Chemical composition and origin of nebulae around Luminous Blue Variables

We use the analysis of the heavy element abundances (C, N, O, S) in circumstellar nebulae around Luminous Blue Variables to infer the evolutionary phase in which the material has been ejected. (1) We discuss the different effects that may have changed the gas composition of the nebula since it was ejected (2) We calculate the expected abundance changes at the stellar surface due to envelope convection in the red supergiant phase. If the observed LBV nebulae are ejected during the RSG phase, the abundances of the LBV nebulae require a significantly smaller amount of mass to be lost than assumed in evolutionary models. (3) We calculate the changes in the surface composition during the main sequence phase by rotation induced mixing. If the nebulae are ejected at the end of the MS-phase, the abundances in LBV nebulae are compatible with mixing times between 5 x 10^6 and 1 x 10^7 years. The existence of ON stars supports this scenario. (4) The predicted He/H ratio in the nebulae are significantly smaller than the current observed photospheric values of their central stars. Combining various arguments we show that the LBV nebulae are ejected during the blue SG phase and that the stars have not gone through a RSG phase. The chemical enhancements are due to rotation induced mixing, and the ejection is possibly triggered by near-critical rotation. During the ejection, the outflow was optically thick, which resulted in a large effective radius and a low effective temperature. This also explains the observed properties of LBV dust.Comment: 18 pages, 4 figures, to be published in The Astrophysical Journal, April 20, 200

The dusty Nebula surrounding HR Car: a Spitzer view

We present mid-IR observations of the Galactic Luminous Blue Variable (LBV) HR Car and its associated nebula carried out with the Spitzer Space Telescope using both IRAC and IRS, as part of a GTO program aimed to study stellar ejecta from evolved stars. Our observations reveal a rich mid-IR spectrum of the inner nebula showing both solid state and atomic gas signatures. Strong low-excitation atomic fine structure lines such as $26.0 \mu$m [\ion{Fe}{2}] and $34.8 \mu$m [\ion{Si}{2}], indicate, for the first time, the presence of a PDR in this object class. While the physics and chemistry of the low-excitation gas appears to be dominated by photodissociation, a possible contribution due to shocks can be inferred from the evidence of gas phase Fe abundance enhancement. The presence of amorphous silicates, inferred from the observed characteristic broad feature at $10 \mu$m located in the inner nebula, suggests that dust has formed during the LBV outburst. This is in contrast with the detection of crystalline dust in other probably more evolved Galactic LBVs, which is similar to the crystalline dust observed in red supergiants. This has been considered to be evidence of dust production during evolutionary phases prior to the outburst.Comment: 27 pages, 6 figures. accepted by Ap

Closed-form sums for some perturbation series involving associated Laguerre polynomials

Infinite series sum_{n=1}^infty {(alpha/2)_n / (n n!)}_1F_1(-n, gamma, x^2), where_1F_1(-n, gamma, x^2)={n!_(gamma)_n}L_n^(gamma-1)(x^2), appear in the first-order perturbation correction for the wavefunction of the generalized spiked harmonic oscillator Hamiltonian H = -d^2/dx^2 + B x^2 + A/x^2 + lambda/x^alpha 0 0, A >= 0. It is proved that the series is convergent for all x > 0 and 2 gamma > alpha, where gamma = 1 + (1/2)sqrt(1+4A). Closed-form sums are presented for these series for the cases alpha = 2, 4, and 6. A general formula for finding the sum for alpha/2 = 2 + m, m = 0,1,2, ..., in terms of associated Laguerre polynomials, is also provided.Comment: 16 page

Star formation history in the SMC: the case of NGC602

Deep HST/ACS photometry of the young cluster NGC 602, located in the remote low density "wing" of the Small Magellanic Cloud, reveals numerous pre-main sequence stars as well as young stars on the main sequence. The resolved stellar content thus provides a basis for studying the star formation history into recent times and constraining several stellar population properties, such as the present day mass function, the initial mass function and the binary fraction. To better characterize the pre-main sequence population, we present a new set of model stellar evolutionary tracks for this evolutionary phase with metallicity appropriate for the Small Magellanic Cloud (Z = 0.004). We use a stellar population synthesis code, which takes into account a full range of stellar evolution phases to derive our best estimate for the star formation history in the region by comparing observed and synthetic color-magnitude diagrams. The derived present day mass function for NGC 602 is consistent with that resulting from the synthetic diagrams. The star formation rate in the region has increased with time on a scale of tens of Myr, reaching $0.3-0.7 \times 10^{-3} M_\odot yr^{-1}$ in the last 2.5 Myr, comparable to what is found in Galactic OB associations. Star formation is most complete in the main cluster but continues at moderate levels in the gas-rich periphery of the nebula.Comment: 24 pages. Accepted for publication in A

Ionization structure in the winds of B[e] supergiants: I. Ionization equilibrium calculations in a H plus He wind

The non-spherically symmetric winds of B[e] supergiants are investigated. An empirical density distribution is chosen that accounts for the density concentrations and ratios derived from observations, and our model winds are assumed to contain only hydrogen and helium. We first calculate the approximate ionization radii for H and He and compare the results with the ionization fractions calculated from the more accurate ionization balance equations. We find that winds with a r^-2 density distribution turn out to reach a constant ionization fraction as long as the wind density is low, i.e. in polar direction. For the high density equatorial regions, however, we find that the winds become neutral just above the stellar surface of the hot and massive B[e] supergiants forming a disk-like neutral region. In such a disk molecules and dust can form even very near the hot central star.Comment: 11 pages, 8 figures, accepted for publication in A&

The red stellar population in NGC 1569

We present HST NICMOS photometry of the resolved stellar population in NGC 1569. The CMD in the F110W and F160W photometric bands contains ~2400 stars with a formal photometric error < 0.1 mag down to J~23.5 and H~22.5. We describe the data processing which is required to calibrate the instrumental peculiarities of NICMOS. Two different packages for PSF-fitting photometry are used to strengthen the photometric results in the crowded stellar field of NGC 1569. The resulting CMD is discussed in terms of the major evolutionary properties of the resolved stellar populations. For a distance modulus of 26.71 and a reddening E(B-V)=0.56, our CMD samples stars down to ~0.8 Mo, corresponding to look-back times > 15 Gyr. This is clear indication of SF activity spanning an entire Hubble time. The metallicity of the reddest RGB stars is in better agreement with Z=0.004 as measured in HII regions, than with Z=0.0004 as expected from the stellar ages. The presence of - yet undetected - very metal-poor stars embedded in the stellar distribution around J=22.75 and J-H=1.15 is, however, not ruled out. The youngest stars (< 50 Myr) are preferentially found around the two central super star clusters, whereas the oldest population has a more uniform spatial distribution. A SFR per unit area of 1 Mo yr*(-1) kpc*(-2) and a mass formed in stars of ~ 1.4x10*6 Mo in the last 50 Myr are derived from the CMD. The NIR CMD places strong constraints on the lower limit of the onset of SF in NGC 1569. The exceptionally high crowding in the NICMOS images of NGC 1569 is a challenge for the photometric analysis. As a result, optical and NIR images of NGC 1569 sample different populations and cannot be cross-correlated. Nevertheless, we demonstrate the consistency of the SF histories derived from the optical and NIR CMDs.Comment: 41 pages including 1 table (Latex) and 14 figures (postscript). Accepted for publication in the Astronomical Journal, March 2001 issu

Optical-uv spectrum and proper motion of the middle-aged pulsar b1055-52

PSRB1055-52 is a middle-aged (~535 kyr) radio, X-ray, and gamma-ray pulsar showing X-ray thermal emission from the neutron star (NS) surface. A candidate optical counterpart to PSRB1055-52 was proposed by Mignani and coworkers based on Hubble Space Telescope (HST) observations performed in 1996, in one spectral band only. We report on HST observations of this field carried out in 2008, in four spectral bands. The astrometric and photometric analyses of these data confirm the identification of the proposed candidate as the pulsar's optical counterpart. Similarly to other middle-aged pulsars, its optical-UV spectrum can be described by the sum of a power-law (PLO) component, presumably emitted from the pulsar magnetosphere, and a Rayleigh-Jeans (RJ) component emitted from the NS surface. The spectral index of the PLO component, alpha_O=1.05+/-0.34, is larger than for other pulsars with optical counterparts. The RJ component, with the brightness temperature TO=(0.66+/-0.10) d_350**2 R_O,13**-2 MK (where d_350 and R_O,13 are the distance to the pulsar in units of 350 pc and the radius of the emitting area in units of 13 km), shows a factor of 4 excess with respect to the extrapolation of the X-ray thermal component into the UV-optical. This hints that the RJ component is emitted from a larger, colder area, and suggests that the distance to the pulsar is smaller than previously thought. From the absolute astrometry of the HST images we measured the pulsar coordinates with a position accuracy of 0.15". From the comparison with previous observations we measured the pulsar proper motion, mu = 42+/-5 mas/yr, which corresponds to a transverse velocity V_t = (70+/-8) d_350 km/s.Comment: 12 pages, 5 figures, accepted for publication on Astrophysical Journal, (Fig1a available at http://www.mssl.ucl.ac.uk/~rm2/

The initial stellar mass function from random sampling in hierarchical clouds II: statistical fluctuations and a mass dependence for starbirth positions and times

Observed variations in the slope of the initial stellar mass function are shown to be consistent with a model in which the protostellar gas is randomly sampled from hierarchical clouds at a rate proportional to the square root of the local density. RMS variations in the IMF slope around the Salpeter value are +/- 0.4 when only 100 stars are observed, and +/- 0.1 when 1000 stars are observed. The hierarchical-sampling model also reproduces the tendency for massive stars to form closer to the center of a cloud, at a time somewhat later than the formation time of the lower mass stars. The assumed density dependence for the star formation rate is shown to be appropriate for turbulence compression, magnetic diffusion, gravitational collapse, and clump or wavepacket coalescence. The low mass flattening in the IMF comes from the inability of gas to form stars below the thermal Jeans mass at typical temperatures and pressures. Consideration of heating and cooling processes indicate why the thermal Jeans mass should be nearly constant in normal environments, and why it might increase in some starburst regions. The steep IMF in the extreme field is not explained by the model, but other origins are suggested.Comment: 21 pages, 8 figures, scheduled for ApJ vol. 515, April 10, 199