Optimization of Starburst99 for Intermediate-Age and Old Stellar Populations

We have incorporated the latest release of the Padova models into the evolutionary synthesis code Starburst99. The Padova tracks were extended to include the full asymptotic giant branch (AGB) evolution until the final thermal pulse over the mass range 0.9 to 5 solar mass. With this addition, Starburst99 accounts for all stellar phases that contribute to the integrated light of a stellar population with arbitrary age from the extreme ultraviolet to the near-infrared. AGB stars are important for ages between 0.1 and 2 Gyr, with their contribution increasing at longer wavelengths. We investigate similarities and differences between the model predictions by the Geneva and the Padova tracks. The differences are particularly pronounced at ages > 1 Gyr, when incompleteness sets in for the Geneva models. We also perform detailed comparisons with the predictions of other major synthesis codes and found excellent agreement. Our synthesized optical colors are compared to observations of old, intermediate-age, and young populations. Excellent agreement is found for the old globular cluster system of NGC 5128 and for old and intermediate-age clusters in NGC 4038/39. In contrast, the models fail for red supergiant dominated populations with sub-solar abundances. This failure can be traced back to incorrect red supergiant parameters in the stellar evolutionary tracks. Our models and the synthesis code are publicly available as version 5.0 of Starburst99 at http://www.stsci.edu/science/starburst99/.Comment: The revised Starburst99 code discussed in this paper will replace the current version 4.0 on our Starburst99 website by December 31, 2004. Accepted for publication in ApJ; 39 pages, 23 figures, 5 table

Integral-Field Spectroscopy of the Post Red Supergiant IRC +10420: evidence for an axi-symmetric wind

We present NAOMI/OASIS adaptive-optics assisted integral-field spectroscopy of the transitional massive hypergiant IRC +10420, an extreme mass-losing star apparently in the process of evolving from a Red Supergiant toward the Wolf-Rayet phase. To investigate the present-day mass-loss geometry of the star, we study the appearance of the line-emission from the inner wind as viewed when reflected off the surrounding nebula. We find that, contrary to previous work, there is strong evidence for wind axi-symmetry, based on the equivalent-width and velocity variations of H$\alpha$ and Fe {\sc ii} $\lambda$6516. We attribute this behaviour to the appearance of the complex line-profiles when viewed from different angles. We also speculate that the Ti {\sc ii} emission originates in the outer nebula in a region analogous to the Strontium Filament of $\eta$ Carinae, based on the morphology of the line-emission. Finally, we suggest that the present-day axisymmetric wind of IRC +10420, combined with its continued blueward evolution, is evidence that the star is evolving toward the B[e] supergiant phase.Comment: 22 pages, 9 figures, accepted for publication in ApJ. B&W-optimized version can be downloaded from http://www.cis.rit.edu/~bxdpci/pubs.htm

Stellar Hydrodynamics in Radiative Regions

We present an analysis of the response of a radiative region to waves generated by a convective region of the star; this wave treatment of the classical problem of ``overshooting'' gives extra mixing relative to the treatment traditionally used in stellar evolutionary codes. The interface between convectively stable and unstable regions is dynamic and nonspherical, so that the nonturbulent material is driven into motion, even in the absence of ``penetrative overshoot.'' These motions may be described by the theory of nonspherical stellar pulsations, and are related to motion measured by helioseismology. Multi-dimensional numerical simulations of convective flow show puzzling features which we explain by this simplified physical model. Gravity waves generated at the interface are dissipated, resulting in slow circulation and mixing seen outside the formal convection zone. The approach may be extended to deal with rotation and composition gradients. Tests of this description in the stellar evolution code TYCHO produce carbon stars on the asymptotic giant branch (AGB), an isochrone age for the Hyades and three young clusters with lithium depletion ages from brown dwarfs, and lithium and beryllium depletion consistent with observations of the Hyades and Pleiades, all without tuning parameters. The insight into the different contributions of rotational and hydrodynamic mixing processes could have important implications for realistic simulation of supernovae and other questions in stellar evolution.Comment: 27 pages, 5 figures, accepted to the Astrophysical Journa

Observational Tests and Predictive Stellar Evolution II: Non-standard Models

We examine contributions of second order physical processes to results of stellar evolution calculations amenable to direct observational testing. In the first paper in the series (Young et al. 2001) we established baseline results using only physics which are common to modern stellar evolution codes. In the current paper we establish how much of the discrepancy between observations and baseline models is due to particular elements of new physics. We then consider the impact of the observational uncertainties on the maximum predictive accuracy achievable by a stellar evolution code. The sun is an optimal case because of the precise and abundant observations and the relative simplicity of the underlying stellar physics. The Standard Model is capable of matching the structure of the sun as determined by helioseismology and gross surface observables to better than a percent. Given an initial mass and surface composition within the observational errors, and no additional constraints for which the models can be optimized, it is not possible to predict the sun's current state to better than ~7%. Convectively induced mixing in radiative regions, seen in multidimensional hydrodynamic simulations, dramatically improves the predictions for radii, luminosity, and apsidal motions of eclipsing binaries while simultaneously maintaining consistency with observed light element depletion and turnoff ages in young clusters (Young et al. 2003). Systematic errors in core size for models of massive binaries disappear with more complete mixing physics, and acceptable fits are achieved for all of the binaries without calibration of free parameters. The lack of accurate abundance determinations for binaries is now the main obstacle to improving stellar models using this type of test.Comment: 33 pages, 8 figures, accepted for publication in the Astrophysical Journa

Chemical Abundances of Planetary Nebulae in the Sagittarius Dwarf Elliptical Galaxy

Spectrophotometry and imaging of the two planetary nebulae He2-436 and Wray16-423, recently discovered to be in the Sagittarius dwarf elliptical galaxy, are presented. Wray16-423 is a high excitation planetary nebula (PN) with a hot central star. In contrast He2-436 is a high density PN with a cooler central star and evidence of local dust, the extinction exceeding that for Wray16-423 by E(B-V)=0.28. The extinction to Wray16-423, (E(B-V)=0.14), is consistent with the extinction to the Sagittarius (Sgr) Dwarf. Both PN show Wolf-Rayet features in their spectra, although the lines are weak in Wray16-423. Images in [O III] and H-alpha+[N II], although affected by poor seeing, yield a diameter of 1.2'' for Wray16-423 after deconvolution; He~2-436 was unresolved. He2-436 has a luminosity about twice that of Wray16-423 and its size and high density suggest a younger PN. In order to reconcile the differing luminosity and nebular properties of the two PN with similar age progenitor stars, it is suggested that they are on He burning tracks The abundance pattern is very similar in both nebulae and shows an oxygen depletion of -0.4 dex with respect to the mean O abundance of Galactic PN and [O/H] = -0.6. The Sgr PN progenitor stars are representative of the higher metallicity tail of the Sgr population. The pattern of abundance depletion is similar to that in the only other PN in a dwarf galaxy companion of the Milky Way, that in Fornax, for which new spectra are presented. However the abundances are larger than for Galactic halo PN suggesting a later formation age. The O abundance of the Sgr galaxy deduced from its PN, shows similarities with that of dwarf ellipticals around M31, suggesting that this galaxy was a dwarf elliptical before its interaction with the Milky Way.Comment: 24 pages, Latex (aas2pp4.sty) including 5 postscript figures. To appear in Ap

Spitzer/IRS spectroscopy of high mass precursors to planetary nebulae

We present Spitzer/IRS observations of a small sample of heavily obscured IRAS sources displaying both the infrared and OH maser emission characteristic of OH/IR stars on the asymptotic giant branch (AGB), but also radio continuum emission typical of ionized planetary nebulae (PNe), the so-called OHPNe. Our observations show that their mid-infrared spectra are dominated by the simultaneous presence of strong and broad amorphous silicate absorption features together with crystalline silicate features, originated in their O-rich circumstellar shells. Out of the five sources observed, three of them are clearly non-variable at infrared wavelengths, confirming their post-AGB status, while the remaining two still show strong photometric fluctuations, and may still have not yet departed from the AGB. One of the non-variable sources in the sample, IRAS 17393-2727, displays a strong [Ne II] nebular emission at 12.8 microns, indicating that the ionization of its central region has already started. This suggests a rapid evolution from the AGB to the PN stage. We propose that these heavily obscured OHPNe represent the population of high mass precursors to PNe in our Galaxy.Comment: To appear in The Astrophysical Journal Letters (scheduled in the 2007 September 1 issue

New optical and near-infrared Surface Brightness Fluctuations models. A primary distance indicator ranging from Globular Clusters to distant galaxies?

We present new theoretical models for Surface Brightness Fluctuations (SBF) both for optical and near-infrared bands in standard ground-based and Hubble Space Telescope filter systems. Simple Stellar Population simulations are adopted. Models cover the age and metallicity ranges from $t=5$ to $15~Gyr$ and from $Z=0.0001$ to 0.04 respectively. Effects due to the variation of the Initial Mass Function and the stellar color-temperature relations are explored. Particular attention is devoted to very bright stars in the color-magnitude diagram and to investigate the effects of mass loss along the Red Giant Branch (RGB) and the Asymptotic Giant Branch (AGB). It is found that $U$ and $B$ bands SBF amplitudes are powerful diagnostics for the morphology of the Horizontal Branch and the Post-AGB stars population. We point out that a careful treatment of mass loss process along the RGB and AGB is fundamental in determining reliable SBF evaluations. The SBF measurements are used to give robust constraints on the evolution of AGB stars, suggesting that mass loss activity on AGB stars should be twice more efficient than on the RGB stars. Our models are able to reproduce the absolute SBF magnitudes of the Galactic Globular Clusters and of galaxies, and their integrated colors. New calibrations of absolute SBF magnitude in $V$, $R$, $I$, and $K$ photometric filters are provided, which appear reliable enough to directly gauge distances bypassing other distance indicators. The SBF technique is also used as stellar population tracer to derive age and metallicity of a selected sample of galaxies of known distances. Finally, {\it SBF color} versus {\it integrated color} diagrams are proposed as particularly useful in removing the well known {\it age-metallicity degeneracy} affecting our knowledge of remote stellar systems.Comment: AJ accepted, 46 pages, 21 figures, 10 tables, uses aastex.cl

The dust envelope of the pre-planetary nebula IRAS19475+3119

We present the spectral energy distribution (SED) of the pre-planetary nebula, IRAS 19475+3119 (I19475), from the optical to the far-infrared. We identify emission features due to crystalline silicates in the ISO SWS spectra of the star. We have fitted the SED of I19475 using a 1-D radiative transfer code, and find that a shell with inner and outer radii of 8.8X10^{16} and 4.4X10^{17}cm, and dust temperatures ranging from about 94K to 46K provide the best fit. The mass of this shell is greater than/equal to 1[34cm^{2}g^{-1}/kappa(100micron)][delta/200]M_Sun, where kappa(100micron) is the 100micron dust mass absorption coefficient (per unit dust mass), and delta is the gas-to-dust ratio. In agreement with results from optical imaging and millimeter-wave observations of CO emission of I19475, our model fits support an r^{-3} density law for its dust shell, with important implications for the interaction process between the fast collimated post-AGB winds and the dense AGB envelopes which results in the observed shapes of PPNs and PNs. We find that the observed JCMT flux at sub-millimeter wavelengths (850micron) is a factor ~ 2 larger than our model flux, suggesting the presence of large dust grains in the dust shell of I19475 which are not accounted for by our adopted standard MRN grain size distribution.Comment: 38 pages, 8 figures. Accepted for publication in Ap

The Ever Changing Circumstellar Nebula Around UW Centauri

We present new images of the reflection nebula surrounding the R Coronae Borealis Star, UW Cen. This nebula, first detected in 1990, has changed its appearance significantly. At the estimated distance of UW Cen, this nebula is approximately 0.6 ly in radius so the nebula cannot have physically altered in only 8 years. Instead, the morphology of the nebula appears to change as different parts are illuminated by light from the central star modulated by shifting thick dust clouds near its surface. These dust clouds form and dissipate at irregular intervals causing the well-known declines in the R Coronae Borealis (RCB) stars. In this way, the central star acts like a lighthouse shining through holes in the dust clouds and lighting up different portions of the nebula. The existence of this nebula provides clues to the evolutionary history of RCB stars possibly linking them to the Planetary Nebulae and the final helium shell flash stars.Comment: To be published in ApJ Letters. 5 pages, 3 figures (2 in color

Shaping bipolar Planetary Nebulae : How mass loss leads to waistline development

Asymptotic Giant Branch (AGB) stars generally have spherically symmetric envelopes, whereas most post-AGB stars and Planetary Nebulae (PNe) show axisymmetric circumstellar envelopes. While various mechanisms for axisymmetric circumstellar structures may explain the shapes of PNe, they do not address how the shape of the circumstellar shell evolves. Here we address the temporal changes in the axisymmetry of AGB star envelopes, and in particular the development of the torus required in the Generalized Interacting Stellar Winds (GISW) model. Assuming (1) an AGB star rotates with sufficient angular speed at the start of the AGB phase; and (2) that the rotational angular momentum of the AGB star is conserved, we demonstrate that some very important observational features of AGB star axisymmetry evolution can be reproduced. We find that, compared to the star's increasing luminosity and decreasing effective temperature, the decreasing mass of the star primarily affects the axisymmetry of the envelope. When a representative mass loss history is adopted, where most of the mass is lost near the end of the AGB phase, the envelope's axisymmetry increases over time, with the strongest increase occurring near the end of the AGB phase. This may naturally explain why most AGB stars have spherically symmetric envelopes, while axisymmetry seems common-place in the post-AGB/PNe phase. The degree of axisymmetry at the end of the AGB phase is found to increase with increasing main sequence mass, and the onset of axisymmetry occurs only after the onset of the superwind (SW) phase, in good agreement with the observations.Comment: 15 pages, 2 figures, accepted by Ap