Abundances in the Galactic bulge: results from planetary nebulae and giant stars

Context. Our understanding of the chemical evolution (CE) of the Galactic bulge requires the determination of abundances in large samples of giant stars and planetary nebulae (PNe). Studies based on high resolution spectroscopy of giant stars in several fields of the Galactic bulge obtained with very large telescopes have allowed important progress. Aims. We discuss PNe abundances in the Galactic bulge and compare these results with those presented in the literature for giant stars. Methods. We present the largest, high-quality data-set available for PNe in the direction of the Galactic bulge (inner-disk/bulge). For comparison purposes, we also consider a sample of PNe in the Large Magellanic Cloud (LMC). We derive the element abundances in a consistent way for all the PNe studied. By comparing the abundances for the bulge, inner-disk, and LMC, we identify elements that have not been modified during the evolution of the PN progenitor and can be used to trace the bulge chemical enrichment history. We then compare the PN abundances with abundances of bulge field giant. Results. At the metallicity of the bulge, we find that the abundances of O and Ne are close to the values for the interstellar medium at the time of the PN progenitor formation, and hence these elements can be used as tracers of the bulge CE, in the same way as S and Ar, which are not expected to be affected by nucleosynthetic processes during the evolution of the PN progenitors. The PN oxygen abundance distribution is shifted to lower values by 0.3 dex with respect to the distribution given by giants. A similar shift appears to occur for Ne and S. We discuss possible reasons for this PNe-giant discrepancy and conclude that this is probably due to systematic errors in the abundance derivations in either giants or PNe (or both). We issue an important warning concerning the use of absolute abundances in CE studies.

H(2) infrared line emission from the ionized region of planetary nebulae

Context. The analysis and interpretation of the H(2) line emission from planetary nebulae have been done in the literature by 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 H(2) emission is produced inside the ionized region of these objects. Aims. The aim of the present work is to calculate and analyze the infrared line emission of H(2) produced inside the ionized region of planetary nebulae using a one-dimensional photoionization code. Methods. The photoionization code Aangaba was improved in order to calculate the statistical population of the H(2) energy levels, as well as the intensity of the H(2) infrared emission lines in the physical conditions typical of planetary nebulae. A grid of models was obtained and the results then analyzed and compared with the observational data. Results. We show that the contribution of the ionized region to the H(2) line emission can be important, particularly in the case of nebulae with high-temperature central stars. This result explains why H(2) 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 in populating the rovibrational levels of the electronic ground state of H(2) molecules. 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 H(2) molecule on the thermal equilibrium of the gas, concluding that, in the ionized region, H(2) 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.

Age distribution of the central stars of galactic disk planetary nebulae

Context. Determination of the ages of central stars of planetary nebulae (CSPN) is a complex problem, and there is presently no single method that can be generally applied. We have developed several methods of estimating the ages of CSPN, based on both the observed nebular properties and some properties of the stars themselves. Aims. Our aim is to estimate the ages and the age distribution of CSPN and to compare the derived results with mass and age determinations of CSPN and white dwarfs based on empirical determinations of these quantities. Methods. We considered a sample of planetary nebulae in the galactic disk, most of which (similar to 69%) are located in the solar neighbourhood, within 3 kpc from the Sun. We discuss several methods of deriving the age distribution of CSPN, namely; (i) the use of an age-metallicity relation that also depends on the galactocentric distance; (ii) the use of an age-metallicity relation obtained for the galactic disk; and (iii) the determination of ages from the central star masses obtained from the observed nitrogen abundances. Results. We estimated the age distribution of CSPN with average uncertainties of 1-2 Gyr, and compared our results with the expected distribution based both on the observed mass distribution of white dwarfs and on the age distribution derived from available mass distributions of CSPN. Based on our derived age distributions, we conclude that most CSPN in the galactic disk have ages under 6 Gyr, and that the age distribution is peaked around 2-4 Gyr.

PLANETARY NEBULAE IN THE INNER MILKY WAY II: THE BULGE-DISK TRANSITION

In this work, a sample of planetary nebulae located in the inner-disk and bulge of the Galaxy is used in order to find the galactocentric distance which better separates these two populations, from the point of view of abundances. Statistical distance scales are used to study the distribution of abundances across the disk-bulge interface. A Kolmogorov-Smirnov test is used to find the distance at which the chemical properties of these regions better separate. The results of the statistical analysis indicate that, on the average, the inner population has lower abundances than the outer. Additionally, for the a-element abundances, the inner population does not follow the disk radial gradient towards the galactic center. Based on our results, we suggest a bulge-disk interface at 1.5 kpc, marking the transition between the bulge and inner-disk of the Galaxy as defined by the intermediate mass population.

PLANETARY NEBULAE IN THE INNER MILKY WAY: NEW ABUNDANCES

The study of planetary nebulae in the inner-disk and bulge gives important information on the chemical abundances of elements such as He, N, O, Ar, Ne, and on the evolution of these abundances, which is associated with the evolution of intermediate-mass stars and the chemical evolution of time Galaxy. We present accurate abundances of the elements He, N, 5, 0, Ar, and Ne for a sample of 54 planetary nebulae located towards the bulge of the Galaxy, for 33 of which the abundances are derived here for the first time. The abundances are obtained based on observations in the optical domain made at the National Laboratory for Astrophysics (LNA, Brazil). The data show a good agreement; with other results in the literature, in the sense that the distribution of the abundances is similar to that of those works.

PLANETARY NEBULAE AND THE CHEMICAL EVOLUTION OF THE MAGELLANIC CLOUDS

The determination of accurate chemical abundances of planetary nebulae (PN) in different galaxies allows us to obtain important constraints on chemical evolution models for these systems. We have a long-term program to derive abundances in the galaxies of the Local Group, particularly the Large and Small Magellanic Clouds. In this work, we present our new results on these objects and discuss their implications in view of recent abundance determinations in the literature. In particular, we obtain distance-independent correlations involving He, N, O, Ne, S, and Ar, and compare the results with data from our own Galaxy and other galaxies in the Local Group. As a result of our observational program, we have a large database of PN in the Galaxy and the Magellanic Clouds, so that we can obtain reliable constraints on the nucleosynthesis processes in the progenitor stars in galaxies of different metallicities.

Molecular chemistry and the missing mass problem in planetary nebulae

Context. Detections of molecular lines, mainly from H-2 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. Aims. 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. Methods. 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. Results. 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(H-2) 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.

Planetary nebulae and determination of the bulge-disk boundary

In this paper, a sample of planetary nebulae in the Galaxy's inner-disk and bulge is used to find the galactocentric distance that optimally separates these two populations in terms of their abundances. Statistical distance scales were used to investigate the distribution of abundances across the disk–bulge interface, while a Kolmogorov–Smirnov test was used to find the distance at which the chemical properties of these regions separate optimally. The statistical analysis indicates that, on average, the inner population is characterized by lower abundances than the outer component. Additionally, for the α-element abundances, the inner population does not follow the disk's radial gradient toward the Galactic Center. Based on our results, we suggest a bulge–disk interface at 1.5 kpc, marking the transition between the bulge and the inner disk of the Galaxy as defined by the intermediate-mass population.

THREE-DIMENSIONAL PHOTOIONIZATION STRUCTURE AND DISTANCES OF PLANETARY NEBULAE. IV. NGC 40

Continuing our series of papers on the three-dimensional (3D) structure and accurate distances of planetary nebulae (PNe), we present here the results obtained for PN NGC 40. Using data from different sources and wavelengths, we construct 3D photoionization models and derive the physical quantities of the ionizing source and nebular gas. The procedure, discussed in detail in the previous papers, consists of the use of 3D photoionization codes constrained by observational data to derive the 3D nebular structure, physical and chemical characteristics, and ionizing star parameters of the objects by simultaneously fitting the integrated line intensities, the density map, the temperature map, and the observed morphologies in different emission lines. For this particular case we combined hydrodynamical simulations with the photoionization scheme in order to obtain self-consistent distributions of density and velocity of the nebular material. Combining the velocity field with the emission-line cubes we also obtained the synthetic position-velocity plots that are compared to the observations. Finally, using theoretical evolutionary tracks of intermediate-and low-mass stars, we derive the mass and age of the central star of NGC 40 as (0.567 +/- 0.06) M(circle dot) and (5810 +/- 600) yr, respectively. The distance obtained from the fitting procedure was (1150 +/- 120) pc.

Metallicity effects on the modified wind momentum of CSPN

Recent investigations on the central stars of planetary nebulae (CSPN) indicate that the masses based on model atmospheres can be much larger than the masses derived from theoretical mass-luminosity relations. Also, the dispersion in the relation between the modified wind momentum and the luminosity depends on the mass spread of the CSPN, and is larger than observed in massive hot stars. Since the wind characteristics probably depend on the metallicity, we analyze the effects on the modified wind momentum by considering the dispersion in this quantity caused by the stellar metallicity. Our CSPN masses are based on a relation between the core mass and the nebular abundances. We conclude that these masses agree with the known mass distribution both for CSPN and white dwarfs, and that the spread in the modified wind momentum can be explained by the observed metallicity variations.

Probing the 2D kinematic structure of early-type galaxies out to three effective radii

We detail an innovative new technique for measuring the two-dimensional (2D) velocity moments (rotation velocity, velocity dispersion and Gauss-Hermite coefficients h(3) and h(4)) of the stellar populations of galaxy haloes using spectra from Keck DEIMOS (Deep Imaging Multi-Object Spectrograph) multi-object spectroscopic observations. The data are used to reconstruct 2D rotation velocity maps. Here we present data for five nearby early-type galaxies to similar to three effective radii. We provide significant insights into the global kinematic structure of these galaxies, and challenge the accepted morphological classification in several cases. We show that between one and three effective radii the velocity dispersion declines very slowly, if at all, in all five galaxies. For the two galaxies with velocity dispersion profiles available from planetary nebulae data we find very good agreement with our stellar profiles. We find a variety of rotation profiles beyond one effective radius, i.e. rotation speed remaining constant, decreasing and increasing with radius. These results are of particular importance to studies which attempt to classify galaxies by their kinematic structure within one effective radius, such as the recent definition of fast- and slow-rotator classes by the Spectrographic Areal Unit for Research on Optical Nebulae project. Our data suggest that the rotator class may change when larger galactocentric radii are probed. This has important implications for dynamical modelling of early-type galaxies. The data from this study are available on-line.

Mapping low- and high-density clouds in astrophysical nebulae by imaging forbidden line emission

Emission line ratios have been essential for determining physical parameters such as gas temperature and density in astrophysical gaseous nebulae. With the advent of panoramic spectroscopic devices, images of regions with emission lines related to these physical parameters can, in principle, also be produced. We show that, with observations from modern instruments, it is possible to transform images taken from density-sensitive forbidden lines into images of emission from high- and low-density clouds by applying a transformation matrix. In order to achieve this, images of the pairs of density-sensitive lines as well as the adjacent continuum have to be observed and combined. We have computed the critical densities for a series of pairs of lines in the infrared, optical, ultraviolet and X-rays bands, and calculated the pair line intensity ratios in the high- and low-density limit using a four- and five-level atom approximation. In order to illustrate the method, we applied it to Gemini Multi-Object Spectrograph (GMOS) Integral Field Unit (GMOS-IFU) data of two galactic nuclei. We conclude that this method provides new information of astrophysical interest, especially for mapping low- and high-density clouds; for this reason, we call it `the ld/hd imaging method`.

The globular cluster kinematics and galaxy dark matter content of NGC 3923

This paper presents further results from our spectroscopic study of the globular cluster (GC) system of the group elliptical NGC 3923. From observations made with the GMOS instrument on the Gemini South Telescope, an additional 50 GC and ultra-compact dwarf (UCD) candidates have been spectroscopically confirmed as members of the NGC 3923 system. When the recessional velocities of these GCs are combined with the 29 GC velocities reported previously, a total sample of 79 GC/UCD velocities is produced. This sample extends to over 6 arcmin (>6 R-e similar to 30 kpc) from the centre of NGC 3923 and is used to study the dynamics of the GC system and the dark matter content of NGC 3923. It is found that the GC system of NGC 3923 displays no appreciable rotation, and that the projected velocity dispersion is constant with radius within the uncertainties. The velocity dispersion profiles of the integrated light and GC system of NGC 3923 are indistinguishable over the region in which they overlap. We find some evidence that the diffuse light and GCs of NGC 3923 have radially biased orbits within similar to 130 arcsec. The application of axisymmetric orbit-based models to the GC and integrated light velocity dispersion profiles demonstrates that a significant increase in the mass-to-light ratio (from M/L-V = 8 to 26) at large galactocentric radii is required to explain this observation. We therefore confirm the presence of a dark matter halo in NGC 3923. We find that dark matter comprises 17.5(-4.5)(+7.3) per cent of the mass within 1 R-e, 41.2(-10.6)(+18.2) per cent within 2 R-e and 75.6(-16.8)(+15.4) per cent within the radius of our last kinematic tracer at 6.9 R-e. The total dynamical mass within this radius is found to be 1.5(-0.25)(+0.4) x 10(12) M-circle dot. In common with other studies of large ellipticals, we find that our derived dynamical mass profile is consistently higher than that derived by X-ray observations, by a factor of around 2.

The qWR star HD 45166 - II. Fundamental stellar parameters and evidence of a latitude-dependent wind

Context. The enigmatic object HD 45166 is a qWR star in a binary system with an orbital period of 1.596 day, and presents a rich emission-line spectrum in addition to absorption lines from the companion star (B7 V). As the system inclination is very small (i = 0.77 degrees +/- 0.09 degrees), HD 45166 is an ideal laboratory for wind-structure studies. Aims. The goal of the present paper is to determine the fundamental stellar and wind parameters of the qWR star. Methods. A radiative transfer model for the wind and photosphere of the qWR star was calculated using the non-LTE code CMFGEN. The wind asymmetry was also analyzed using a recently-developed version of CMFGEN to compute the emerging spectrum in two-dimensional geometry. The temporal-variance spectrum (TVS) was calculated to study the line-profile variations. Results. Abundances and stellar and wind parameters of the qWR star were obtained. The qWR star has an effective temperature of T(eff) = 50 000 +/- 2000 K, a luminosity of log(L/L(circle dot)) = 3.75 +/- 0.08, and a corresponding photospheric radius of R(phot) = 1.00 R(circle dot). The star is helium-rich (N(H)/N(He) = 2.0), while the CNO abundances are anomalous when compared either to solar values, to planetary nebulae, or to WR stars. The mass-loss rate is. M = 2.2 x 10(-7) M(circle dot) yr(-1), and the wind terminal velocity is v(infinity) = 425 km s(-1). The comparison between the observed line profiles and models computed under different latitude-dependent wind densities strongly suggests the presence of an oblate wind density enhancement, with a density contrast of at least 8: 1 from equator to pole. If a high velocity polar wind is present (similar to 1200 km s(-1)), the minimum density contrast is reduced to 4:1. Conclusions. The wind parameters determined are unusual when compared to O-type stars or to typical WR stars. While for WR stars v(infinity)/v(esc) > 1.5, in the case of HD 45166 it is much smaller (v(infinity)/v(esc) = 0.32). In addition, the efficiency of momentum transfer is eta = 0.74, which is at least 4 times smaller than in a typical WR. We find evidence for the presence of a wind compression zone, since the equatorial wind density is significantly higher than the polar wind. The TVS supports the presence of such a latitude-dependent wind and a variable absorption/scattering gas near the equator.

The stellar content of obscured Galactic giant H II regions. VI. W51A

We present K-band spectra of newly born OB stars in the obscured Galactic giant H II region W51A and approximate to 0.8 '' angular resolution images in the J, H, and K(S)-bands. Four objects have been spectroscopically classified as O-type stars. The mean spectroscopic parallax of the four stars gives a distance of 2.0 +/- 0.3 kpc (error in the mean), significantly smaller than the radio recombination line kinematic value of 5.5 kpc or the values derived from maser proper motion observations (6-8 kpc). The number of Lyman continuum photons from the contribution of all massive stars (NLyc approximate to 1.5 x 10(50) s(-1)) is in good agreement with that inferred from radio recombination lines (NLyc = 1.3 x 10(50) s(-1)) after accounting for the smaller distance derived here. We present analysis of archival high angular resolution images (NAOS CONICA at VLT and T-ReCS at Gemini) of the compact region W51 IRS 2. The K(S)-band images resolve the infrared source IRS 2 indicating that it is a very young compact H II region. Sources IRS 2E was resolved into compact cluster (within 660 AU of projected distance) of three objects, but one of them is just bright extended emission. W51d1 and W51d2 were identified with compact clusters of three objects (maybe four in the case of W51d1) each one. Although IRS 2E is the brightest source in the K-band and at 12.6 mu m, it is not clearly associated with a radio continuum source. Our spectrum of IRS 2E shows, similar to previous work, strong emission in Br gamma and He I, as well as three forbidden emission lines of Fe III and emission lines of molecular hydrogen (H(2)) marking it as a massive young stellar object.