Convective envelopes in rotating OB stars

We study the effects of rotation on the outer convective zones of massive stars. We examine the effects of rotation on the thermal gradient and on the Solberg--Hoiland term by analytical developments and by numerical models. Writing the criterion for convection in rotating envelopes, we show that the effects of rotation on the thermal gradient are much larger and of opposite sign to the effect of the Solberg-Hoiland criterion. On the whole, rotation favors convection in stellar envelopes at the equator and to a smaller extent at the poles. In a rotating 20 Msun star at 94% of the critical angular velocity, there are two convective envelopes, with the bigger one having a thickness of 13.2% of the equatorial radius. In the non-rotating model, the corresponding convective zone has a thickness of only 4.6% of the radius. The occurrence of outer convection in massive stars has many consequences.Comment: 4 pages, 3 figures, accepted by Astronomy and Astrophysic

Visualization and spectral synthesis of rotationally distorted stars

Simple spherical, non-rotating stellar models are inadequate when describing real stars in the limit of very fast rotation: Both the observable spectrum and the geometrical shape of the star deviate strongly from simple models. We attempt to approach the problem of modeling geometrically distorted, rapidly rotating stars from a new angle: By constructing distorted geometrical models and integrating standard stellar models with varying temperature, gravity, and abundances, over the entire surface, we attempt a semi-empirical approach to modeling. Here we present our methodology, and present simple examples of applications.Comment: 7 pages, 5 figures. Contribution to GREAT-ESF Workshop on "Stellar Atmospheres in the Gaia Era", Brussels 23-24 June 201

An Extremely Lithium-Rich Bright Red Giant in the Globular Cluster M3

We have serendipitously discovered an extremely lithium-rich star on the red giant branch of the globular cluster M3 (NGC 5272). An echelle spectrum obtained with the Keck I HIRES reveals a Li I 6707 Angstrom resonance doublet of 520 milli-Angstrom equivalent width, and our analysis places the star among the most Li-rich giants known: log[epsilon(Li)] ~= +3.0. We determine the elemental abundances of this star, IV-101, and three other cluster members of similar luminosity and color, and conclude that IV-101 has abundance ratios typical of giants in M3 and M13 that have undergone significant mixing. We discuss mechanisms by which a low-mass star may be so enriched in Li, focusing on the mixing of material processed by the hydrogen-burning shell just below the convective envelope. While such enrichment could conceivably only happen rarely, it may in fact regularly occur during giant-branch evolution but be rarely detected because of rapid subsequent Li depletion.Comment: 7-page LaTeX file, including 2 encapsulated ps figures + 1 table; accepted for publication in the Astrophysical Journal Letter

PHOEBE 2.0 – Where no model has gone before

phoebe 2.0 is an open source framework bridging the gap between stellar observations and models. It allows to create and fit models simultaneously and consistently to a wide range of observational data such as photometry, spectroscopy, spectrapolarimetry, interferometry and astrometry. To reach the level of precision required by the newest generation of instruments such as Kepler, GAIA and the arrays of large telescopes, the code is set up to handle a wide range of phenomena such as multiplicity, rotation, pulsations and magnetic fields, and to model the involved physics to a new level

Artificial Intelligence Approach to the Determination of Physical Properties of Eclipsing Binaries. I. The EBAI Project

Achieving maximum scientific results from the overwhelming volume of astronomical data to be acquired over the next few decades will demand novel, fully automatic methods of data analysis. Artificial intelligence approaches hold great promise in contributing to this goal. Here we apply neural network learning technology to the specific domain of eclipsing binary (EB) stars, of which only some hundreds have been rigorously analyzed, but whose numbers will reach millions in a decade. Well-analyzed EBs are a prime source of astrophysical information whose growth rate is at present limited by the need for human interaction with each EB data-set, principally in determining a starting solution for subsequent rigorous analysis. We describe the artificial neural network (ANN) approach which is able to surmount this human bottleneck and permit EB-based astrophysical information to keep pace with future data rates. The ANN, following training on a sample of 33,235 model light curves, outputs a set of approximate model parameters (T2/T1, (R1+R2)/a, e sin(omega), e cos(omega), and sin i) for each input light curve data-set. The whole sample is processed in just a few seconds on a single 2GHz CPU. The obtained parameters can then be readily passed to sophisticated modeling engines. We also describe a novel method polyfit for pre-processing observational light curves before inputting their data to the ANN and present the results and analysis of testing the approach on synthetic data and on real data including fifty binaries from the Catalog and Atlas of Eclipsing Binaries (CALEB) database and 2580 light curves from OGLE survey data. [abridged]Comment: 52 pages, accepted to Ap

The Structure of the Homunculus. III. Forming a Disk and Bipolar Lobes in a Rotating Surface Explosion

We present a semi-analytic model for shaping the nebula around eta Carinae that accounts for the simultaneous production of bipolar lobes and an equatorial disk through a rotating surface explosion. Material is launched normal to the surface of an oblate rotating star with an initial kick velocity that scales approximately with the local escape speed. Thereafter, ejecta follow ballistic orbital trajectories, feeling only a central force corresponding to a radiatively reduced gravity. Our model is conceptually similar to the wind-compressed disk model of Bjorkman & Cassinelli, but we modify it to an explosion instead of a steady line-driven wind, we include a rotationally-distorted star, and we treat the dynamics somewhat differently. Continuum-driving avoids the disk inhibition that normally operates in line-driven winds. Our model provides a simple method by which rotating hot stars can simultaneously produce intrinsically bipolar and equatorial mass ejections, without an aspherical environment or magnetic fields. Although motivated by eta Carinae, the model may have generic application to other LBVs, B[e] stars, or SN1987A's nebula. When near-Eddington radiative driving is less influential, our model generalizes to produce bipolar morphologies without disks, as seen in many PNe.Comment: ApJ accepted, 9 page

The nearby eclipsing stellar system delta Velorum - I. Origin of the infrared excess from VISIR and NACO imaging

- Context: The triple stellar system delta Vel system presents a significant infrared excess, whose origin is still being debated. A large infrared bow shock has been discovered using Spitzer/MIPS observations. Although it appears as a significant contributor to the measured IR excess, the possibility exists that a circumstellar IR excess is present around the stars of the system. - Aims: The objective of the present VISIR and NACO observations is to identify whether one of the stars of the delta Vel system presents a circumstellar photometric excess in the thermal IR domain and to quantify it. - Methods: We observed delta Vel using the imaging modes of the ESO/VLT instruments VISIR (in BURST mode) and NACO to resolve the A-B system (0.6" separation) and obtain the photometry of each star. We also obtained one NACO photometry epoch precisely at the primary (annular) eclipse of delta Vel Aa by Ab. - Results: Our photometric measurements with NACO (2.17 mic), complemented by the existing visible photometry allowed us to reconstruct the spectral energy distribution of the three stars. We then compared the VISIR photometry (8.6-12.8 mic) to the expected photospheric emission from the three stars at the corresponding wavelengths. - Conclusions: We can exclude the presence of a circumstellar thermal infrared excess around delta Vel A or B down to a few percent level. This supports the conclusions of Gaspar et al. (2008) that the IR excess of delta Vel has an interstellar origin, although a cold circumstellar disk could still be present. In addition, we derive the spectral types of the three stars Aa, Ab, and B (respectively A2IV, A4V and F8V), and we estimate the age of the system around 400-500 Myr.Comment: 8 pages, 9 figures, A&A, in pres

Surface temperature and synthetic spectral energy distributions for rotationally deformed stars

The spectral energy distribution (SED) of a non-spherical star could differ significantly from the SED of a spherical star with the same average temperature and luminosity. Calculation of the SED of a deformed star is often approximated as a composite of several spectra, each produced by a plane parallel model of given effective temperature and gravity. The weighting of these spectra over the stellar surface, and hence the inferred effective temperature and luminosity, will be dependent on the inclination of the rotation axis of the star with respect to the observer, as well as the temperature and gravity distribution on the stellar surface. Here we calculate the surface conditions of rapidly rotating stars with a 2D stellar structure and evolution code and compare the effective temperature distribution to that predicted by von Zeipel's law. We calculate the composite spectrum for a deformed star by interpolating within a grid of intensity spectra of plane parallel model atmospheres and integrating over the surface of the star. Using this method, we find that the deduced variation of effective temperature with inclination can be as much as 3000 K for an early B star, depending on the details of the underlying model.Comment: 38 pages, 9 figures (AAStex preprint format). Accepted for publication in the Ap

Genetic-Algorithm-based Light Curve Optimization Applied to Observations of the W UMa star BH Cas

I have developed a procedure utilizing a Genetic-Algorithm-based optimization scheme to fit the observed light curves of an eclipsing binary star with a model produced by the Wilson-Devinney code. The principal advantages of this approach are the global search capability and the objectivity of the final result. Although this method can be more efficient than some other comparably global search techniques, the computational requirements of the code are still considerable. I have applied this fitting procedure to my observations of the W UMa type eclipsing binary BH Cassiopeiae. An analysis of V-band CCD data obtained in 1994/95 from Steward Observatory and U- and B-band photoelectric data obtained in 1996 from McDonald Observatory provided three complete light curves to constrain the fit. In addition, radial velocity curves obtained in 1997 from McDonald Observatory provided a direct measurement of the system mass ratio to restrict the search. The results of the GA-based fit are in excellent agreement with the final orbital solution obtained with the standard differential corrections procedure in the Wilson-Devinney code.Comment: 9 pages, 2 figures, 2 tables, uses emulateapj.st

CN and CH Band Strengths of Bright Giants in M3

CN and CH band strengths for ten bright red giants in M3 have been measured from archival spectra obtained with the MMT. A CN-CH band strength anticorrelation is confirmed for the program stars together with other stars for which there is published data. This suggests an anticorrelation between carbon and nitrogen abundances with a constant total abundance of carbon plus nitrogen. However, stars which do not follow the CN-CH anticorrelation are also found. The star III-77, which is found to have both strong CN and CH bands, is the most peculiar among them. While three other stars, VZ194, VZ352, and VZ1420, which show both weak CN and CH band strengths could be AGB stars.Comment: 10 pages, 5 figures, 1 tabl