Search

Accurate mass-loss rates are essential for meaningful stellar evolutionary models. For massive single stars with initial masses between 8 - 30\msun the implementation of cool supergiant mass loss in stellar models strongly affects the resulting evolution, and the most commonly used prescription for these cool-star phases is that of de Jager. Recently, we published a new \mdot\ prescription calibrated to RSGs with initial masses between 10 - 25\msun, which unlike previous prescriptions does not over estimate \mdot\ for the most massive stars. Here, we carry out a comparative study to the MESA-MIST models, in which we test the effect of altering mass-loss by recomputing the evolution of stars with masses 12-27\msun\ with the new \mdot-prescription implemented. We show that while the evolutionary tracks in the HR diagram of the stars do not change appreciably, the mass of the H-rich envelope at core-collapse is drastically increased compared to models using the de Jager prescription. This increased envelope mass would have a strong impact on the Type II-P SN lightcurve, and would not allow stars under 30\msun\ to evolve back to the blue and explode as H-poor SN. We also predict that the amount of H-envelope around single stars at explosion should be correlated with initial mass, and we discuss the prospects of using this as a method of determining progenitor masses from supernova light curves., 7 pages, 3 figures, accepted for publication in ApJ

Neutron star X-ray binaries (NS XRBs) associated with supernova remnants (SNRs) are the youngest X-ray binaries that can provide insights into the early evolution of X-ray binaries and the formation properties of neutron stars. There are an increasing number of NS XRBs discovered to be harbored in SNRs in our and nearby galaxies. In this work, we perform binary population synthesis calculations to simulate the population of NS XRBs associated with SNRs for different types of companions, including Roche-lobe overfilling main-sequence stars, Be stars and supergiants. We estimate their birthrates and present the distributions of orbital parameters and companion mass for each type of companions. Our calculations show that the majority of them are Be X-ray binaries (BeXRBs) and that a few BeXRBs are expected to be associated with SNRs in a Milky Way-type galaxy., Comment: Accepted for publication in ApJ

A failed core-collapse supernova from a non-rotating progenitor can eject mass due to a weakening of gravity associated to neutrino emission by the protoneutron star. This mechanism yields observable transients and sets an upper limit to the mass of the black hole (BH) remnant. Previous global simulations of this mechanism have included neutrino losses parametrically, however, with direct implications for the ejecta mass and energy. Here we evolve the inner supernova core with a spherically-symmetric, general-relativistic neutrino radiation-hydrodynamic code until BH formation. We then use the result in a Newtonian code that follows the response of the outer layers of the star to the change in gravity and resolves the surface pressure scale height. We find that the dense-matter equation of state (EOS) can introduce a factor $\sim 2$ variation in gravitational mass lost to neutrinos, with a stiff EOS matching previous parametric results, and a soft EOS yielding lower ejecta masses and energies by a factor of several. This difference is caused primarily by the longer time to BH formation in stiffer EOSs. With a soft EOS, our red and yellow supergiant progenitors fail to unbind mass if hydrogen recombination energy is not included. Using a linear ramp in time for mass-energy lost to neutrinos (with suitable parameters) yields a stellar response within $\sim 10\%$ of that obtained using the detailed history of neutrino losses. Our results imply quantitative but not qualitative modifications to previous predictions for shock breakout, plateau emission, and final BH masses from these events., Comment: Accepted by ApJ with minor corrections. New equation (6) quantifies the difference between baryonic and gravitational mass in the pre-supernova progenitor

Analysis of the recent long exposure Chandra X-ray observation of the early-type O star zeta Pup shows clear variability with a period previously reported in optical photometric studies. These 813 ks of HETG observations taken over a roughly one year time span have two signals of periodic variability: a high significance period of 1.7820 +/- 0.0008 day, and a marginal detection of periodic behavior close to either 5 day or 6 day period. A BRITE-Constellation nanosatellite optical photometric monitoring, using near-contemporaneous observations to the Chandra data, confirms a 1.78060 +/- 0.00088 day period for this star. The optical period coincides with the new Chandra period within their error ranges, demonstrating a link between these two wavebands and providing a powerful lever for probing the photosphere/wind connection in this star. The phase lag of the X-ray maximum relative to the optical maximum is approximately phi=0.45, but consideration of secondary maxima in both datasets indicates possibly two hot spots on the star with an X-ray phase lag of phi=0.1 each. The details of this periodic variation of the X-rays are probed by displaying a phased and trailed X-ray spectrum and by constructing phased light curves for wavelength bands within the HETG spectral coverage, ranging down to bands encompassing groups of emission lines. We propose that the 1.78 day period is the stellar rotation period and explore how stellar bright spots and associated co-rotating interacting regions or CIRs could explain the modulation of the optical and X-ray output for this star and their phase difference., 17 pages; 10 figures

The distribution of stars in the Hertzsprung Russell diagram (HRD) for a stellar conglomeration represents a snapshot of its evolving stellar population. Some of the supergiant stars may transit the HRD from blue to red and then again to blue during their late evolutionary stages, as exemplified by the progenitor of SN 1987A. Others may transit a given part of the HRD more than twice in a "blue loop" and end up as red supergiants before they explode. Since stars in blue loops spend a considerable part of their lives there, these stages may change the relative number of modeled supergiants in the HRD. Their lifetimes in turn depend upon the initial mass of the star, how convection in its interior is modeled, and how much mass loss takes place during its evolution. The observed ratio of the number of blue to red supergiants and yellow to red supergiants sensitively test the stellar evolution theory. We compare modeled number ratios of these supergiants with observed data from the Large Magellanic Cloud as it has a metallicity very similar to that of the environment of SN 2013ej. We successfully model these by taking into account moderate (exponential) convective overshooting. We explore its effect on the final radius and mass of the star prior to core collapse. The radius differs dramatically with overshoot. These factors controlling pre-supernova structure may affect the post-explosion optical/IR light curves and spectral development., Comment: 14 pages, 3 figures, 3 tables, accepted at the ApJ for publication

Cassiopeia A (Cas A) is one of the best studied young Galactic supernova remnants. While providing a rare opportunity to study in detail the remnant of a Type IIb supernova, questions remain regarding the nature of its progenitor, its mass-loss history, and its pre-SN evolution. Here we present an optical investigation of the circumstellar environment around Cas A and find clumpy and filamentary Halpha emission nebulosities concentrated 10-15 pc (10-15 arcminutes) to the north and east. First reported by Minkowski as a faint H II region, these nebulosities exhibit distinct morphological and spectroscopic properties relative to the surrounding diffuse emissions. Compared to neighboring H II regions, these nebulae show stronger [N II] 6548, 6583 A and [S II] 6716, 6731 A emissions relative to Halpha. We show that Cas A's highest-velocity ejecta knots are interacting with some of the closest projected emission nebulae, thus providing strong evidence that these nebulae lie at the same distance as the remnant. We interpret these surrounding nebulosities to be the remains of the progenitor's red supergiant wind which accumulated against the southern edge of a large extended H II region located north of Cas A. Our findings are consistent with the view that Cas A's progenitor underwent considerable mass-loss, first from a fast main-sequence wind, then from a slower, clumpy red supergiant wind, and finally from a brief high-velocity wind, like that from a yellow supergiant., Comment: 19 pages, 9 figures, accepted to ApJ

Magnetic configurations, stable on the long term, appear to exist in various evolutionary phases, from Main-Sequence stars to white dwarfs and neutron stars. The large scale ordered nature of these fields, often approximately dipolar, and their scaling according to the flux conservation scenario favor the model of a fossil field (Duez et al. 2010). We make some first estimates of the magnetic coupling between the stellar cores and the outer layers in red giants and supergiants. Analytical expressions of the truncation radius of the field coupling are established for a convective envelope and for a rotating radiative zone with horizontal turbulence. The timescales of the internal exchanges of angular momentum are considered. Numerical estimates are made on the basis of recent model grids. The direct magnetic coupling of the core to the extended convective envelope of red giants and supergiants appears unlikely. However, we find that the intermediate radiative zone is fully coupled to the core during the He-burning and later phases. This coupling is able to produce a strong spin down of the core of red giants and supergiants, also leading to relatively slowly rotating stellar remnants, like white dwarfs and pulsars. Some angular momentum is also transferred to the outer convective envelope of red giants and supergiants during the He-burning phase and later., Comment: 11 pages, accepted for publication in ApJ

We present the evolution of massive star progenitors of supernovae of type IIP. We take the example of the nearby and well-studied SN 2013ej. We explore how convective overshoot affects the stellar structure, surface abundances, and effective temperature of massive stars, using the Modules for Experiments in Stellar Astrophysics (MESA). In particular models with moderate overshoot ($f$ = 0.02 to 0.031) show the presence of blue loops in the Hertzsprung-Russell diagram with a red to blue [$log_{10}(T_{eff}/\rm K)$ from $< 3.6$ to $> 4.0$] excursion and transition back to red, during core helium burning phase. Models with overshoot outside this range of $f$ values kept the star in the red supergiant state throughout the post helium ignition phases. The surface CNO abundance shows enrichment post-main-sequence and again around the time when helium is exhausted in core. These evolutionary changes in surface CNO abundance are indistinguishable in the currently available observations due to large observational uncertainties. However, these observations may distinguish between the ratios of surface nitrogen to oxygen at different evolutionary stages of the star. We also compare the effects of convective overshoot on various parameters related to likelihood of explosion of a star as opposed to collapse to a black hole. These parameters are the compactness parameter, M$_4$, and $��_4$. Combination $��_4 \times $M$_4$, and $��_4$ have similar variation with $f$ and both peak out at $f$ = 0.032. We find that all of our 13 M$_{\odot}$ models are likely to explode., 17 pages, 8 Figures, 2 Tables, Accepted for publication in The Astrophysical Journal

The hot massive luminous blue variables (LBVs) represent an important evolutionary phase of massive stars. Here, we report the discovery of a new LBV -- LAMOST J0037+4016 in the distant outskirt of the Andromeda galaxy. It is located in the south-western corner (a possible faint spiral arm) of M31 with an unexpectedly large projection distance of $\sim$ 22 kpc from the center. The optical light curve shows a 1.2 mag variation in $V$ band and its outburst and quiescence phases both last over several years. The observed spectra indicate an A-type supergiant at epoch close to the outburst phase and a hot B-type supergiant with weak [Fe II] emission lines at epoch of much dimmer brightness. The near-infrared color-color diagram further shows it follows the distribution of Galactic and M31 LBVs rather than B[e] supergiants. All the existing data strongly show that LAMOST J0037+4016 is an LBV. By spectral energy distribution fitting, we find it has a luminosity ($4.42 \pm 1.64$)$\times 10^5$ $L_{\odot}$ and an initial mass $\sim 30$ $M_{\odot}$, indicating its nature of less luminosity class of LBV., Comment: 7 pages, 4 figures, 3 tables, accepted by ApJL

Supergiant X-ray binaries usually comprise a neutron star accreting from the wind of a OB supergiant companion. They are classified as classical systems and the supergiant fast X-ray transients (SFXTs). The different behavior of these sub-classes of sources in X-rays, with SFXTs displaying much more pronounced variability, is usually (at least) partly ascribed to different physical properties of the massive star clumpy stellar wind. In case of SFXTs, a systematic investigation of the effects of clumps on flares/outbursts of these sources has been reported by Bozzo et al. (2017) exploiting the capabilities of the instruments on-board XMM-Newton to perform a hardness-resolved spectral analysis on timescales as short as a few hundreds of seconds. In this paper, we use six XMM-Newton observations of IGR J18027-2016 to extend the above study to a classical supergiant X-ray binary and compare the findings with those derived in the case of SFXTs. As these observations of IGR J18027-2016 span different orbital phases, we also study its X-ray spectral variability on longer timescales and compare our results with previous publications. Although obtaining measurements of the clump physical properties from X-ray observations of accreting supergiant X-ray binaries was already proven to be challenging, our study shows that similar imprints of clumps are found in the X-ray observations of the supergiant fast X-ray transients and at least one classical system, i.e. IGR J18027-2016. This provides interesting perspectives to further extend this study to many XMM-Newton observations already performed in the direction of other classical supergiant X-ray binaries., 13 pages, 6 figures, 2 tables; Accepted for publication in ApJ

High cadence ultraviolet, optical and near-infrared photometric and low-resolution spectroscopic observations of the peculiar Type II supernova (SN) 2018hna are presented. The early phase multiband light curves exhibit the adiabatic cooling envelope emission following the shock breakout up to ~14 days from the explosion. SN~2018hna has a rise time of $\sim$\,88 days in the V-band, similar to SN 1987A. A $\rm^{56}Ni$ mass of ~0.087$\pm$0.004 $\rm M_{\odot}$ is inferred for SN 2018hna from its bolometric light curve. Hydrodynamical modelling of the cooling phase suggests a progenitor with a radius ~50 $\rm R_{\odot}$, a mass of ~14-20 $\rm M_{\odot}$ and explosion energy of ~1.7-2.9$\rm \times$ $\rm 10^{51}\ erg$. The smaller inferred radius of the progenitor than a standard red supergiant is indicative of a blue supergiant progenitor of SN 2018hna. A sub-solar metallicity (~0.3 $\rm Z_{\odot}$) is inferred for the host galaxy UGC 07534, concurrent with the low-metallicity environments of 1987A-like events., 8 pages, 5 figures, Accepted for Publication in ApJ Letters

Recent observations of supernovae (SNe) just after the explosion suggest that a good fraction of SNe have the confined circumstellar material (CSM) in the vicinity, and the pre-SN enhanced mass loss may be a common property. The physical mechanism of this phenomenon is still unclarified, and the energy deposition into the envelope has been proposed as a possible cause of the confined CSM. In this work, we have calculated the response of the envelope to various types of sustained energy deposition starting from a few years before the core collapse. We have further investigated how the resulting progenitor structure would affect appearance of the ensuing supernova. While it has been suspected that a super-Eddington energy deposition may lead to a strong and/or eruptive mass loss to account for the confined CSM, we have found that a highly super-Eddington energy injection into the envelope changes the structure of the progenitor star substantially, and the properties of the resulting SNe become inconsistent with usual SNe. This argument constrains the energy budget involved in the possible stellar activity in the final years to be at most one order of magnitude higher than the Eddington luminosity. Such an energy generation however would not dynamically develop a strong wind in the time scale of a few years. We therefore propose that a secondary effect (e.g., pulsation or binary interaction) triggered by the moderate envelope inflation, which is caused by sub-Eddington energy injection, likely induces the mass loss., Comment: 18 pages, 14 figures, accepted for publication in ApJ

IGR J16493-4348 is an eclipsing supergiant high-mass X-ray binary (sgHMXB), where accretion onto the compact object occurs via the radially outflowing stellar wind of its early B-type companion. We present an analysis of the system's X-ray variability and periodic modulation using pointed observations (2.5-25 keV) and Galactic bulge scans (2-10 keV) from the Rossi X-ray Timing Explorer (RXTE) Proportional Counter Array (PCA), along with Swift Burst Alert Telescope (BAT) 70-month snapshot (14-195 keV) and transient monitor (15-50 keV) observations. The orbital eclipse profiles in the PCA bulge scans and BAT light curves are modeled using asymmetric and symmetric step and ramp functions. We obtain an improved orbital period measurement of 6.7828 $\pm$ 0.0004 days from an observed minus calculated (O-C) analysis of mid-eclipse times derived from the BAT transient monitor and PCA scan data. No evidence is found for the presence of a strong photoionization or accretion wake. We refine the superorbital period to 20.067 $\pm$ 0.009 days from the discrete Fourier transform (DFT) of the BAT transient monitor light curve. A pulse period of 1093.1036 $\pm$ 0.0004 s is measured from a pulsar timing analysis using pointed PCA observations spanning $\sim$1.4 binary orbits. We present pulse times of arrival (ToAs), circular and eccentric timing models, and calculations of the system's Keplerian binary orbital parameters. We derive an X-ray mass function of $f_{x}(M)$ $=$ 13.2$^{+2.4}_{-2.5}$ $M_{\odot}$ and find a spectral type of B0.5 Ia for the supergiant companion through constraints on the mass and radius of the donor. Measurements of the eclipse half-angle and additional parameters describing the system geometry are provided., Comment: 21 pages, 13 figures, 9 tables, Accepted for publication in ApJ on 2018 November 9

Multiple stellar populations are observed in almost all globular-clusters, but the origin of this phenomenon is still debated. We investigate the role cool supergiants may have played. To do this, we combine two investigative methods: state-of-the-art massive stellar evolution and calculations of the hydrodynamic structure of the cluster-gas. This approach allows us to study how star-formation in young massive clusters depends on the energy- and mass-input of the first-generation of stars, while predicting the chemical composition of the second-generation. We find that the presence of massive (9-500 M$_{\odot}$) metal-poor supergiants in the young cluster leads to a star-formation episode within the first 4 Myr of the cluster's lifetime, that is, before the first core-collapse supernovae explode or the gas is expelled. The stellar winds accumulate in the cluster center, forming the second-generation there. Its composition is predicted to show variations in O & Na abundances, consistently with observations. The abundance of helium is, similarly to other scenarios involving massive stars, higher than what is referred from observations. Supposing dynamical removal of stars from the outskirts of the cluster, or applying a top-heavy initial-mass-function, we can predict a number ratio of the second-generation as high as 20-80%. The effect of metallicity is shown to be important, as the most luminous supergiants are only predicted at low-metallicity, thus limiting $-$ but not excluding $-$ the extent of a polluted second-generation at high-metallicity. These massive stars becoming black-holes suggests globular-clusters hosting gravitational-wave progenitors. Our scenario predicts a correlation between the mass of the cluster and the extent of the multiple population phenomenon., Accepted by ApJ, in Press

Understanding the mechanism of core-collapse supernova explosions requires knowledge of the nuclear equation of state (EoS). Recent multi-dimensional numerical simulations indicate that explosions are possible. Nevertheless, it is not yet fully understood which equation of state is realized in the proto-neutron star formed during SN explosions. We examine the EoS sensitivity of the relic supernova neutrino spectrum as a probe of the nuclear EoS. This sensitivity arises largely from the contribution to neutrino emission from failed supernovae. We consider a variety of astrophysical scenarios, which include different progenitor masses for a successful explosion, the cosmological star formation rate, starbursts, quiescent star formation, and the metallicity dependence of the initial mass function. We find that the EoS signature remains robust under a variety of conditions. We demonstrate the viability of future neutrino detectors to distinguish the nuclear EoS via the relic supernova neutrino spectrum.

The photon-scattering winds of M-giants absorb parts of the chromospheric emission lines and produce self-reversed spectral features in high resolution {\it HST}/GHRS spectra. These spectra provide an opportunity to assess fundamental parameters of the wind, including flow and turbulent velocities, the optical depth of the wind above the region of photon creation, and the star's mass-loss rate. This paper is the last paper in the series `GHRS Observations of Cool, Low-Gravity Stars'; the last several have compared empirical measurements of spectral emission lines with models of the winds and mass-loss of K-giant and supergiants. We have used the Sobolev with Exact Integration (SEI) radiative transfer code, along with simple models of the outer atmosphere and wind, to determine and compare the wind characteristics of the two M-giant stars, $\gamma$~Cru (M3.5III) and $\mu$~Gem (M3IIIab), with previously derived values for low-gravity K-stars. The analysis specifies the wind parameters and calculates line profiles for the \ion{Mg}{2} resonance lines, in addition to a range of unblended \ion{Fe}{2} lines. Our line sample covers a large range of wind opacities and, therefore, probes a range of heights in the atmosphere. Our results show that $\mu$~Gem has a slower and more turbulent wind then $\gamma$~Cru. Also, $\mu$~Gem has weaker chromosphere, in terms of surface flux, with respect to $\gamma$~Cru. This suggests that $\mu$~Gem is more evolved than $\gamma$~Cru. Comparing the two M-giants in this work with previously studied K-giant and supergiant stars ($\alpha$~Tau, $\gamma$~Dra, $\lambda$~Vel) reveals that the M-giants have slower winds than the earlier giants, but exhibit higher mass-loss rates. Our results are interpreted in the context of the winds being driven by Alfv\'{e}n waves., Comment: Accepted for publications on ApJ, 17 pages, 7 figures

Very recently, the optical bands of the 2Σ+ ← X2Σ+system of AlO have been identified in the red supergiant star VYCMa. In an effort to explain the origin of this transition, we used state-of-the-art quantum chemical calculations with proven high accuracy to compute the lowest singlet and triplet electronic states of the AlOH and HAlO isomers as well as their equilibrium geometry and electronic properties. Our calculated potential energy surfaces implicate the three singlet electronic states 21A', 31A', and 11A'' in the photodissociation of the [Al,O,H] system. Only AlO, H, Al, and OH products can occur through the photodissociation of [Al,O,H]; AlH and O are not allowed. For the photodissociation of AlOH, the AlO product can occur only in its excited states AlO(2Π) and AlO(2Σ+).

The long-standing puzzle of reversed-C bisectors of spectral lines is addressed with the tools of the third-signature plot and bisector mapping. The F8 supergiant γ Cyg, which shows strongly reversed bisectors, is analyzed and found to have a monotonic decline in velocities with height in its photosphere, consistent with a convective overshoot region similar to what is seen for stars on the cool side of the granulation boundary. Further, a flux deficit is derived that is not unlike those for stars on the cool side of the granulation boundary, although the γ Cyg flux deficit is wider and peaks at a higher velocity. The reversed-C bisector shape is caused by a combination of a steep decline in velocities with height, as shown in the third-signature plot, with a flux deficit spanning only a fraction of the red wing of the line profiles, in contrast to cooler stars where the deficit extends over most of the red wing. Apparently no unusual velocity fields or other bizarre behaviors are needed to explain the granulation boundary; it stems more simply from the continuous changes of normal convective overshoot across the HR diagram.

The yellow supergiant content of nearby galaxies provides a critical test of massive star evolutionary theory. While these stars are the brightest in a galaxy, they are difficult to identify because a large number of foreground Milky Way stars have similar colors and magnitudes. We previously conducted a census of yellow supergiants within M31 and found that the evolutionary tracks predict a yellow supergiant duration an order of magnitude longer than we observed. Here we turn our attention to the SMC, where the metallicity is 10x lower than that of M31, which is important as metallicity strongly affects massive star evolution. The SMC's large radial velocity (~160 km/s) allows us to separate members from foreground stars. Observations of ~500 candidates yielded 176 near-certain SMC supergiants, 16 possible SMC supergiants, along with 306 foreground stars and provide good relative numbers of yellow supergiants down to 12Mo. Of the 176 near-certain SMC supergiants, the kinematics predicted by the Besancon model of the Milky Way suggest a foreground contamination of >4%. After placing the SMC supergiants on the H-R diagram and comparing our results to the Geneva evolutionary tracks, we find results similar to those of the M31 study: while the locations of the stars on the H-R diagram match the locations of evolutionary tracks well, the models over-predict the yellow supergiant lifetime by a factor of ten. Uncertainties about the mass-loss rates on the main-sequence thus cannot be the primary problem with the models., Accepted by the ApJ

Pulsations driven by partial ionization of hydrogen in the envelope are often considered important for driving winds from red supergiants (RSGs). In particular, it has been suggested by some authors that the pulsation growth rate in a RSG can be high enough to trigger an unusually strong wind (or a super-wind), when the luminosity to mass ratio becomes sufficiently large. Using both hydrostatic and hydrodynamic stellar evolution models with initial masses ranging from 15 to 40 \Msun, we investigate 1) how the pulsation growth rate depends on the global parameters of supergiant stars, and 2) what would be the consequences of a pulsation-driven super-wind, if it occurred, for the late stages of massive star evolution. We suggest that such a super-wind history would be marked by a runaway increase, followed by a sudden decrease, of the winds mass loss rate. The impact on the late evolution of massive stars would be substantial, with stars losing a huge fraction of their H-envelope even with a significantly lower initial mass than previously predicted. This might explain the observed lack of Type II-P supernova progenitors having initial mass higher than about 17 \Msun. We also discuss possible implications for a subset of Type IIn supernovae.

Ultra-long gamma ray bursts (ULGRBs) are a distinct class of GRBs characterized by durations of several thousands of seconds, about two orders of magnitude longer than those of standard long GRBs (LGRBs). The driving engine of these events has not been uncovered yet, and ideas range from magnetars, to tidal disruption events, to extended massive stars, such as blue super giants (BSG). BSGs, a possible endpoint of stellar evolution, are attractive for the relatively long free-fall times of their envelopes, allowing accretion to power a long-lasting central engine. At the same time, their large radial extension poses a challenge to the emergence of a jet. Here we perform an end-to-end simulation aimed at assessing the viability of BSGs as ULGRB progenitors. The evolution to core collapse of a BSG star model is calculated with the MESA code. We then compute the accretion rate for the fraction of envelope material with enough angular momentum to circularize and form an accretion disk, and input the corresponding power into a jet which we evolve through the star envelope with the FLASH code. Our simulation shows that the jet can emerge, and the resulting light curves resemble those observed in ULGRBs, with durations $T_{90}$ ranging from $\approx$4000~s to $\approx10^4$ s depending on the viewing angle., Comment: Matches version in press in ApJ

We present the results of modeling dust SEDs across the SMC with the aim of mapping the distribution of PAHs in a low-metallicity environment. Using Spitzer Survey of the SMC (S3MC) photometry from 3.6-160 um over the main star-forming regions of the Wing and Bar along with spectral mapping from 5-38 um from the Spitzer Spectroscopic Survey of the SMC (S4MC) in selected regions, we model the dust SED and emission spectrum to determine the fraction of dust in PAHs across the SMC. We use the regions of overlapping photometry and spectroscopy to test the reliability of the PAH fraction as determined from SED fits alone. The PAH fraction in the SMC is low compared to the Milky Way and variable--with relatively high fractions (q_PAH~1-2%) in molecular clouds and low fractions in the diffuse ISM (=0.6%). We use the map of PAH fraction across the SMC to test a number of ideas regarding the production, destruction and processing of PAHs in the ISM. We find weak or no correlation between the PAH fraction and the distribution of carbon AGB stars, the location of supergiant H I shells and young SN remnants, or the turbulent Mach number. We find that the PAH fraction is correlated with CO intensity, peaks in the dust surface density and the molecular gas surface density as determined from 160 um emission. The PAH fraction is high in regions of active star-formation, as predicted by its correlation with molecular gas, but is suppressed in H II regions. Because the PAH fraction in the diffuse ISM is generally very low--in accordance with previous work on modeling the SED of the SMC--and the PAH fraction is relatively high in molecular regions, we suggest that PAHs are destroyed in the diffuse ISM of the SMC and/or PAHs are forming in molecular clouds. We discuss the implications of these observations for our understanding of the PAH life cycle, particularly in low-metallicity galaxies., Comment: accepted for publication in ApJ, abstract abridged

We report the identification of a source coincident with SN 2009kr in HST pre-explosion images. The object appears to be a single point source with an intrinsic colour V-I = 1.1 and M_V = -7.6. If this is a single star it would be a yellow supergiant of log L/L_{sol} \sim 5.1 and a mass of 15 (+5/-4) M_{sol}. The spatial resolution does not allow us yet to definitively determine if the progenitor object is a single star, a binary system, or a compact cluster. We show that the early lightcurve is flat, similar to IIP SNe, but that the the spectra are somewhat peculiar, displaying unusual P-Cygni profiles. The evolution of the expanding ejecta will play an important role in understanding the progenitor object., Comment: Submitted to ApJ Letters

We present the first three-dimensional (3D) simulations of the large-scale mixing that takes place in the shock-heated stellar layers ejected in the explosion of a 15.5 solar-mass blue supergiant star. The outgoing supernova shock is followed from its launch by neutrino heating until it breaks out from the stellar surface more than two hours after the core collapse. Violent convective overturn in the post-shock layer causes the explosion to start with significant asphericity, which triggers the growth of Rayleigh-Taylor (RT) instabilities at the composition interfaces of the exploding star. Deep inward mixing of hydrogen (H) is found as well as fast-moving, metal-rich clumps penetrating with high velocities far into the H-envelope of the star as observed, e.g., in the case of SN 1987A. Also individual clumps containing a sizeable fraction of the ejected iron-group elements (up to several 0.001 solar masses) are obtained in some models. The metal core of the progenitor is partially turned over with Ni-dominated fingers overtaking oxygen-rich bullets and both Ni and O moving well ahead of the material from the carbon layer. Comparing with corresponding 2D (axially symmetric) calculations, we determine the growth of the RT fingers to be faster, the deceleration of the dense metal-carrying clumps in the He and H layers to be reduced, the asymptotic clump velocities in the H-shell to be higher (up to ~4500 km/s for the considered progenitor and an explosion energy of 10^{51} ergs, instead of, 15 pages, 8 figures, 30 eps files; significantly extended and more figures added after referee comments; accepted by The Astrophysical Journal

We investigate the relationship between R Coronae Borealis (RCB) stars and hydrogen-deficient carbon (HdC) stars by measuring precise 16O/18O ratios for five cool RCB stars. The 16O/18O ratios are derived by spectrum synthesis from high-resolution (R=50,000) K-band spectra. Lower limits to the 16O/17O and 14N/15}N ratios as well as Na and S abundances (when possible) are also given. RCB stars in our sample generally display less 18O than HdC stars - the derived 16O/18O ratios range from 3 to 20. The only exception is the RCB star WX CrA, which seems to be a HdC-like star with 16O/18O=0.3. Our result of a higher 16O/18O ratio for the RCB stars must be accounted for by a theory of the formation and evolution of HdC and RCB stars. We speculate that a late dredge-up of products of He-burning, principally 12C and 16O, may convert a 18O-rich HdC star into a 18O-poor RCB star as the H-deficient star begins its final evolution from a cool supergiant to the top of the white dwarf cooling track., Accepted for publication in the Astrophysical Journal; 21 pages, 3 Tables, 11 Figures

We present the analysis of near infrared (NIR), adaptive optics (AO) Subaru and archived HST imaging data of a region near the northern middle lobe (NML) of the Centaurus A (Cen A) jet, at a distance of $\sim15$ kpc north-east (NE) from the center of NGC5128. Low-pass filtering of the NIR images reveals strong -- $>3\sigma$ above the background mean -- signal at the expected position of the brightest star in the equivalent HST field. Statistical analysis of the NIR background noise suggests that the probability to observe $>3\sigma$ signal at the same position, in three independent measurements due to stochastic background fluctuations alone is negligible ($\leq10^{-7}\%$) and, therefore, that this signal should reflect the detection of the NIR counterparts of the brightest HST star. An extensive photometric analysis of this star yields $V-I$, visual-NIR, and NIR colors expected from a yellow supergiant (YSG) with an estimated age $\sim10^{+4}_{-3}$ Myr. Furthermore, the second and third brighter HST stars are, likely, also supergiants in Cen A, with estimated ages $\sim16^{+6}_{-3}$ Myr and $\sim25^{+15}_{-9}$ Myr, respectively. The ages of these three supergiants are in good agreement with the ages of the young massive stars that were previously found in the vicinity and are thought to have formed during the later phases of the jet-HI cloud interaction that appears to drive the star formation (SF) in the region for the past $\sim100$ Myr., Comment: 11 pages, 6 figures, 2 tables, accepted for publication in ApJ

Fifty archival radio observations of the supergiant binary Cyg OB2 #5 using the Very Large Array over 20 years are re-examined to determine the location and character of the previously detected variable radio emission. The radio emission from the system consists of a primary component that is associated with the binary, and a non-thermal source (NE) that has been ascribed to a wind-collision region (WCR) between the stellar winds of the binary and that of a B-type star (Star D) to the NE. NE shows no evidence of variation in 23 epochs where it is resolved separately from the primary radio component, demonstrating that the variable emission arises in the primary. Since NE is non-variable, the radio flux from the primary can now be well determined for the first time, most especially in observations that do not resolve both the primary and NE components. The variable radio emission from the primary has a period of 6.7+/-0.3 years which is described by a simple model of a non-thermal source orbiting within the stellar wind envelope of the binary. Such a model implies the presence of a third, unresolved stellar companion (Star C) orbiting the 6.6-day binary with a period of 6.7 years. The variable non-thermal emission arises from either a WCR between Star C and the binary system, or possibly from Star C directly. The model gives a mass-loss rate of 3.4 x 10^{-5} solar mass/yr for Cyg OB2 #5, unusually high for an Of supergiant and comparable to that of WR stars, and consistent with an unusually strong He I 1.083-micron emission line, also redolent of WR stars. An examination of radial velocity observations suggests reflex motion of the binary due to Star C. The natures of NE and Star D are also examined. (abridged), Comment: Accepted for publication in ApJ Dec 1, 2009. 13 pages, 11 figures, 4 tables

The circumstellar envelope of the hypergiant star IRC+10420 has been traced as far out in SiO J=2-1 as in CO J = 1-0 and CO J = 2-1, in dramatic contrast with the centrally condensed (thermal) SiO- but extended CO-emitting envelopes of giant and supergiant stars. Here, we present an observation of the circumstellar envelope in SiO J=1-0 that, when combined with the previous observation in {\sioii}, provide more stringent constraints on the density of the SiO-emitting gas than hitherto possible. The emission in SiO peaks at a radius of $\sim$2\arcsec\ whereas that in SiO J=2-1 emission peaks at a smaller radius of $\sim$1\arcsec, giving rise to their ring-like appearances. The ratio in brightness temperature between SiO J=1-0 and SiO J=2-1 decreases from a value well above unity at the innermost measurable radius to about unity at radius of $\sim$2\arcsec, beyond which this ratio remains approximately constant. Dividing the envelope into three zones as in models for the CO J = 1-0 and CO J = 2-1 emission, we show that the density of the SiO-emitting gas is comparable with that of the CO-emitting gas in the inner zone, but at least an order of magnitude higher by comparison in both the middle and outer zones. The SiO-emitting gas therefore originates from dense clumps, likely associated with the dust clumps seen in scattered optical light, surrounded by more diffuse CO-emitting interclump gas. We suggest that SiO molecules are released from dust grains due to shock interactions between the dense SiO-emitting clumps and the diffuse CO-emitting interclump gas., Comment: Accepted for publication in ApJ

The inclusion of rotationally induced mixing in stellar evolution can alter the structure and composition of pre-supernova stars. We survey the effects of progenitor rotation on nucleosynthetic yields in Population III and II supernovae (SNe) using the new adaptive mesh refinement code CASTRO. We examine piston-driven spherical explosions in 15, 25, and 40 M{sub sun} stars at Z = 0 and 10{sup -4} Z{sub sun} with three explosion energies and two rotation rates. Rotation in the Z = 0 models resulted in primary nitrogen production and a stronger hydrogen burning shell which led all models to die as red supergiants (in contrast to the blue supergiant progenitors made without rotation). On the other hand, the Z = 10{sup -4} Z{sub sun} models that included rotation ended their lives as compact blue stars. Because of their extended structure, the hydrodynamics favors more mixing and less fallback in the metal-free stars than the Z = 10{sup -4} models. As expected, higher energy explosions produce more enrichment and less fallback than do lower energy explosions, and at constant explosion energy, less massive stars produce more enrichment and leave behind smaller remnants than do more massive stars. We compare our nucleosynthetic yields to themore » chemical abundances in the three most iron-poor stars yet found and reproduce the abundance pattern of one, HE 0557-4840, with a zero metallicity, 15 M{sub sun}, 2.4 x 10{sup 51} erg SN. A Salpeter IMF-averaged integration of our yields for Z = 0 models with explosion energies of 2.4 x 10{sup 51} erg or less is in good agreement with the abundances observed in larger samples of extremely metal-poor (EMP) stars, provided 15 M{sub sun} stars are included. Since the abundance patterns of EMP stars likely arise from a representative sample of progenitors, our yields suggest that 15-40 M{sub sun} core-collapse SNe with moderate explosion energies contributed the bulk of the metals to the early universe.« less

We present HST/NICMOS photometry, and low-resolution K-band spectra of the GLIMPSE9 stellar cluster. The newly obtained color-magnitude diagram shows a cluster sequence with H-Ks =1 mag, indicating an interstellar extinction Aks=1.6\pm0.2 mag. The spectra of the three brightest stars show deep CO band-heads, which indicate red supergiants with spectral type M1-M2. Two 09-B2 supergiants are also identified, which yield a spectrophotometric distance of 4.2\pm0.4 kpc. Presuming that the population is coeval, we derive an age between 15 and 27 Myr, and a total cluster mass of 1600\pm400 Msun, integrated down to 1 Msun. In the vicinity of GLIMPSE9 are several HII regions and SNRs, all of which (including GLIMPSE 9) are probably associated with a giant molecular cloud (GMC) in the inner galaxy. GLIMPSE9 probably represents one episode of massive star formation in this GMC. We have identified several other candidate stellar clusters of the same complex.

We report optical observations of the luminous blue variable (LBV) HR Carinae which show that the star has reached a visual minimum phase in 2009. More importantly, we detected absorptions due to SiIV λλ4088-4116. To match their observed line profiles from 2009 May, a high rotational velocity of v rot ≃ 150 ∼ 20 km s-1 is needed (assuming an inclination angle of 30°), implying that HR Car rotates at ≃0.88 ∼ 0.2 of its critical velocity for breakup (vcrit). Our results suggest that fast rotation is typical in all strong-variable, bona fide galactic LBVs, which present S-Dor-type variability. Strong-variable LBVs are located in a well-defined region of the HR diagram during visual minimum (the "LBV minimum instability strip"). We suggest this region corresponds to where v crit is reached. To the left of this strip, a forbidden zone with vrot/vcrit>1 is present, explaining why no LBVs are detected in this zone. Since dormant/ex LBVs like P Cygni and HD 168625 have low v rot, we propose that LBVs can be separated into two groups: fast-rotating, strong-variable stars showing S-Dor cycles (such as AG Car and HR Car) and slow-rotating stars with much less variability (such as P Cygni and HD 168625). We speculate that supernova (SN) progenitors which had S-Dor cycles before exploding (such as in SN 2001ig, SN 2003bg, and SN 2005gj) could have been fast rotators. We suggest that the potential difficulty of fast-rotating Galactic LBVs to lose angular momentum is additional evidence that such stars could explode during the LBV phase., Instituto de Astrofísica de La Plata, Facultad de Ciencias Astronómicas y Geofísicas

Ultraluminous X-ray sources are non-nuclear point sources in external galaxies with LX = 2 × 1039-1041 erg s–1, and thus are possibly stellar black holes with special radiation mechanisms or intermediate mass black holes of 100-105 M ☉. To measure their dynamical mass, one needs to identify the secondary in the optical followed by photometric and spectroscopic monitoring. The counterpart light is usually contaminated, sometimes even dominated, by the X-ray irradiated accretion disk, which complicates the derivation of the properties of the secondary itself. Here, we report a study on M101-ULX1 with 26 Hubble Space Telescope (HST) observations and 33 X-ray observations over 16 years. There were three multi-band HST observations at different optical states over eight years, the spectral energy distributions (SEDs) of which cannot be fitted by the previously claimed secondary of a B supergiant or any other single star. All SEDs can be decomposed into the same stellar component plus a second component due to X-ray irradiation with different contributions at different optical states. The secondary in ULX1 is most likely a WR star as revealed by the SED decompositions and confirmed by the presence of He II λ4686 emission in the counterpart spectrum. Five X-ray outbursts and two optical peaks are revealed from X-ray and optical light curves, and optical peaks are identified as the optical response to X-ray outbursts. The correlated optical and X-ray variability analysis is suggestive of ellipsoidal modulation with a half-amplitude of 0.05 mag and a period in the range of 1.5-4.5 days. Future photometric and spectroscopic observations are required to confirm the orbital period and measure the dynamical mass for ULX1.

We report on Suzaku observations of four unidentified sources from the INTEGRAL and Swift BAT Galactic plane surveys. All the sources have a large neutral hydrogen column density and are likely members of an emerging class of heavily absorbed high mass X-ray binary (HMXB) first identified in INTEGRAL observations. Two of the sources in our sample are approximately constant flux sources, one source shows periodic variation and one source exhibited a short, bright X-ray outburst. The periodicity is transient, suggesting it is produced by a neutron star in an elliptical orbit around a stellar wind source. We analyze the flaring source in several segments to look for spectral variation and discuss the implications of the findings for the nature of the source. We conclude that all four sources in our sample can be identified with the emerging class of highly absorbed HMXBs, that one is a newly identified transient X-ray pulsar and that at least one is a newly identified supergiant fast X-ray transient (SFXT)., 22 pages, 5 figures, submitted to ApJ

We explore the sensitivity of massive stars to neutrino magnetic moments. We find that the additional cooling due to the neutrino magnetic moments bring about qualitative changes to the structure and evolution of stars in the mass window 7 Msun < M < 18 Msun, rather than simply changing the time scales for the burning. We describe some of the consequences of this modified evolution: the shifts in the threshold masses for creating core-collapse supernovae and oxygen-neon-magnesium white dwarfs and the appearance of a new type of supernova in which a partial carbon-oxygen core explodes within a massive star. The resulting sensitivity to the magnetic moment is at the level of (2-4) * 10^{-11} ��_B., 13 pages, 9 figures; a relevant reference added, typos corrected, results unchanged; submitted to ApJ

Observations of the carbon-bearing molecules CO, HCN, CS, HNC, CN, and HCO+ have been conducted toward the circumstellar envelope of the oxygen-rich red supergiant star, VY Canis Majoris (VY CMa), using the Arizona Radio Observatory (ARO). CO and HCN were also observed toward the O-rich shells of NML Cyg, TX Cam, IK Tau, and W Hya. Rotational transitions of these species at 1 mm, 0.8 mm, and 0.4 mm were measured with the ARO Submillimeter Telescope, including the J = 6 → 5 line of CO at 691 GHz toward TX Cam and W Hya. The ARO 12 m was used for 2 mm and 3 mm observations. Four transitions were observed for HCO+ in VY CMa, the first definitive identification of this ion in a circumstellar envelope. Molecular line profiles from VY CMa are complex, indicating three separate outflows: a roughly spherical flow and separate red- and blueshifted winds, as suggested by earlier observations. Spectra from the other sources appear to trace a single outflow component. The line data were modeled with a radiative transfer code to establish molecular abundances relative to H2 and source distributions. Abundances for CO derived for these objects vary over an order of magnitude, f ~ 0.4-5 × 10–4, with the lower values corresponding to the supergiants. For HCN, a similar range in abundance is found (f ~ 0.9-9 × 10–6), with no obvious dependence on the mass-loss rate. In VY CMa, HCO+ is present in all three outflows with f ~ 0.4-1.6 × 10–8 and a spatial extent similar to that of CO. HNC is found only in the red- and blueshifted components with [HCN]/[HNC] ~ 150-190, while [CN]/[HCN] ~ 0.01 in the spherical flow. All three velocity components are traced in CS, which has a confined spatial distribution and f ~ 2-6 × 10–7. These observations suggest that carbon-bearing molecules in O-rich shells are produced by a combination of photospheric shocks and photochemistry. Shocks may play a more prominent role in the supergiants because of their macroturbulent velocities.

A luminous optical transient (OT) that appeared in NGC 300 in early 2008 had a maximum brightness, M_V ~ -12 to -13, intermediate between classical novae and supernovae. We present ground-based photometric and spectroscopic monitoring and adaptive-optics imaging of the OT, as well as pre- and post-outburst space-based imaging with HST and Spitzer. The optical spectrum at maximum showed an F-type supergiant photosphere with superposed emission lines of hydrogen, Ca II, and [Ca II], similar to the spectra of low-luminosity Type IIn "supernova impostors" like SN 2008S, as well as cool hypergiants like IRC +10420. The emission lines have a complex, double structure, indicating a bipolar outflow with velocities of ~75 km/s. The luminous energy released in the eruption was ~10^47 ergs, most of it emitted in the first 2 months. By registering new HST images with deep archival frames, we have precisely located the OT site, and find no detectable optical progenitor brighter than broad-band V magnitude 28.5. However, archival Spitzer images reveal a bright, non-variable mid-IR pre-outburst source. We conclude that the NGC 300 OT was a heavily dust-enshrouded luminous star, of ~10-15 Msun, which experienced an eruption that cleared the surrounding dust and initiated a bipolar wind. The progenitor was likely an OH/IR source which had begun to evolve on a blue loop toward higher temperatures, but the precise cause of the outburst remains uncertain., Revised version, accepted by ApJ Letters. 7 pages, 4 figures

We have identified 17 A-type stars in the Galactic Plane that have mid-IR excesses at 8 micron. From the observed colors in the (r'-H_alpha)-(r'-i') plane, we first identified 23050 early A-type main sequence (MS) star candidates in the Isaac Newton Photometric H-Alpha Survey (IPHAS) point source database that are located in Spitzer GLIMPSE Galactic Plane fields. Imposing the requirement that they be detected in all seven 2MASS and IRAC bands led to a sample of 2692 candidate A-type stars with fully sampled 0.6 to 8 micron SEDs. Optical classification spectra of 18 of the IPHAS candidate A-type MS stars showed that all but one could be well fitted using main sequence A-type templates, with the other being an A-type supergiant. Out of the 2692 A-type candidates 17 (0.6%) were found to have 8-micron excesses above the expected photospheric values. Taking into account non-A-Type contamination estimates, the 8-micron excess fraction is adjusted to ~0.7%. The distances to these sources range from 0.7-2.5 kpc. Only 10 out of the 17 excess stars had been covered by Spitzer MIPSGAL survey fields, of which 5 had detectable excesses at 24 micron. For sources with excesses detected in at least two mid-IR wavelength bands, blackbody fits to the excess SEDs yielded temperatures ranging from 270 to 650 K, and bolometric luminosity ratios L_IR/L* from 2.2x10^{-3}-1.9x10^{-2}, with a mean value of 7.9x10^{-3} (these bolometric luminosities are lower limits as cold dust is not detectable by this survey). Both the presence of mid-IR excesses and the derived bolometric luminosity ratios are consistent with many of these systems being in the planet-building transition phase between the early protoplanetary disk phase and the later debris disk phase.

A new circumstellar molecule, the radical AlO (X 2Σ+), has been detected toward the envelope of the oxygen-rich supergiant star VY Canis Majoris (VY CMa) using the Arizona Radio Observatory (ARO). The N = 7 → 6 and 6 → 5 rotational transitions of AlO at 268 and 230 GHz were observed at 1 mm using the ARO Submillimeter Telescope (SMT) and the N = 4 → 3 line was detected at 2 mm using the ARO 12 m telescope. Based on the shape of the line profiles, AlO most likely arises from the dust-forming region in the spherical outflow of VY CMa, as opposed to the blue or redshifted winds, with a source size of θ s ~ 05. Given this source size, the column density of AlO was found to be N tot ~ 2 × 1015 cm–2 for T rot ~ 230 K, with a fractional abundance, relative to H2, of ~10–8. Gas-phase thermodynamic equilibrium chemistry is the likely formation mechanism for AlO in VY CMa, but either shocks disrupt the condensation process into Al2O3, or chemical "freezeout" occurs. The species therefore survives further into the circumstellar envelope to a radius of r ~ 20 R *. The detection of AlO in VY CMa is additional evidence of an active gas-phase refractory chemistry in oxygen-rich envelopes, and suggests such objects may be fruitful sources for other new oxide identifications.

We present 26 point-sources discovered with Chandra within 200" (~20kpc) of the center of the barred supergiant galaxy NGC 1365. The majority of these sources are high-mass X-ray binaries, containing a neutron star or a black hole accreting from a luminous companion at a sub-Eddington rate. Using repeat Chandra and XMM-Newton as well as optical observations, we discuss in detail the natures of two highly-variable ultraluminous X-ray sources (ULXs): NGC 1365 X1, one of the most luminous ULXs known since the ROSAT era, which is X-ray variable by a factor of 30, and NGC 1365 X2, a newly discovered transient ULX, variable by a factor of >90. Their maximum X-ray luminosities (3-5 x 10^40 erg/s, measured with Chandra) and multiwavelength properties suggest the presence of more exotic objects and accretion modes: accretion onto intermediate mass black holes (IMBHs) and beamed/super-Eddington accretion onto solar-mass compact remnants. We argue that these two sources have black-hole masses higher than those of the typical primaries found in X-ray binaries in our Galaxy (which have masses of, 18 pages, accepted by ApJ

We present a measurement of metallicity in the Galactic center Quintuplet Cluster made using quantitative spectral analysis of two Luminous Blue Variables (LBVs). The analysis employs line-blanketed NLTE wind/atmosphere models fit to high-resolution near-infrared spectra containing lines of H, HeI, SiII, MgII, and FeII. We are able to break the H/He ratio vs. mass-loss rate degeneracy found in other LBVs and to obtain robust estimates of the He content of both objects. Our results indicate solar iron abundance and roughly twice solar abundance in the alpha-elements. These results are discussed within the framework of recent measurements of oxygen and carbon composition in the nearby Arches Cluster and iron abundances in red giants and supergiants within the central 30 pc of the Galaxy. The relatively large enrichment of alpha-elements with respect to iron is consistent with a history of more nucleosynthesis in high mass stars than the Galactic disk.

The observed X-ray source temperature distributions in OB stellar winds, as determined from high energy resolution Chandra observations, show that the highest temperatures occur near the star, and then steadily decrease outward through the wind. To explain this unexpected behavior, we propose a shock model concept that utilizes a well-known magnetic propulsion mechanism; the surface ejection of "diamagnetic plasmoids" into a diverging external magnetic field. This produces rapidly accelerating self-contained structures that plow through an ambient wind and form bow shocks that generate a range in X-ray temperatures determined by the plasmoid-wind relative velocities. The model free parameters are the plasmoid initial Alfven speed, the initial plasma-beta of the external medium, and the divergence rate of the external field. These are determined by fitting the predicted bow shock temperatures with the observed OB supergiant X-ray temperature distribution. We find that the initial external plasma-beta has a range between 0 and 2, and the assumed radially-decreasing external magnetic field strength that scales as r^{-S} has a value of S lying between 2 and 3. Most importantly, the initial plasmoid Alfven speed is found to be well-constrained at a value of 0.6 times the terminal velocity, which appears to represent an upper limit for all normal OB stars. This intriguing new limit on OB magnetic properties, as derived from Chandra observations, emphasizes the need for further studies of magnetic propulsion mechanisms in these stars., Comment: Accepted for publication in ApJ Letters, 4 pages, 4 figures (color), 1 table

A survey of the 12C/13C ratio toward circumstellar envelopes has been conducted at millimeter wavelengths using the facilities of the Arizona Radio Observatory (ARO). The ratios were obtained for a sample of local C- and O-rich asymptotic giant branch and supergiant stars from observations of the 12C and 13C isotopologues of CO and CN, respectively. The J = 1 → 0 transitions of both molecules were observed at λ = 3 mm using the ARO 12 m telescope, while the J = 2 → 1 lines of the two species were measured using the ARO Sub-Millimeter Telescope (SMT) at λ = 1 mm. The 12C/13C ratios were determined from the CO data by modeling both transitions simultaneously with a circumstellar radiative transfer code, which can account for the high opacities present in the emission from this species. In the case of CN, the hyperfine structure was used to evaluate opacity effects. Ratios obtained independently from CO and CN are in good agreement. For the C-rich envelopes, the ratios fall in the range 12C/13C ~ 25–90, while the O-rich shells have values of 10-35. Ratios of 12C/13C ~ 3–14 are found for the supergiant stars, with the exception of VY CMa, where the values lie in the range 25–46. All ratios obtained in this study are ≤ 89, the solar value, suggesting that substantial carbon-13 enrichment may be currently occurring in the local interstellar medium. A qualitative model was constructed based on first and third dredge-up convective mixing that can reproduce the observed ratios. Substantial mixing of H-burning products must occur to explain the ratios in the O-rich objects, while a wide range of 12C/13C values can be generated by only a few percent mixing of He-burning ashes in the C-rich case. The 12C/13C ratios obtained in this study should help improve stellar yield models and contribute to the understanding of Galactic chemical evolution.

IGR J17544-2619 and XTE J1739-302 are considered the prototypical sources of the new class of High Mass X-ray Binaries, the Supergiant Fast X-ray Transients (SFXTs).These sources were observed during bright outbursts on 2008 March 31 and 2008 April 8, respectively, thanks to an on-going monitoring campaign we are performing with Swift, started in October 2007. Simultaneous observations with XRT and BAT allowed us to perform for the first time a broad band spectroscopy of their outbursts. The X-ray emission is well reproduced with absorbed cutoff powerlaws, similar to the typical spectral shape from accreting pulsars. IGR J17544-2619 shows a significantly harder spectrum during the bright flare (where a luminosity in excess of 1E36 erg/s is reached) than during the long-term low level flaring activity (1E33-1E34 erg/s), while XTE J1739-302 displayed the same spectral shape, within the uncertainties, and a higher column density during the flare than in the low level activity. The light curves of these two SFXTs during the bright flare look similar to those observed during recent flares from other two SFXTs, IGRJ11215-5952 and IGRJ16479-4514, reinforcing the connection among the members of this class of X-ray sources.

Supergiant Fast X-ray Transients (SFXTs) are a new class of HMXBs discovered thanks to the monitoring of the Galactic plane performed with the INTEGRAL satellite in the last 5 years. These sources display short outbursts (significantly shorter than typical Be/X-ray binaries) with a peak luminosity of a few 1E36 erg/s. The quiescent level, measured only in a few sources, is around 1E32 erg/s. We are performing a monitoring campaign with Swift of four SFXTs (IGRJ16479-4514, XTEJ1739-302, IGRJ17544-2619 and AXJ1841.0-0536/IGRJ18410-0535). We report on the first four months of Swift observations, started on 2007 October 26. We detect a low level X-ray activity in all four SFXTs which demonstrates that these transient sources accrete matter even outside their outbursts. This fainter X-ray activity is composed of many flares with a large flux variability, on timescales of thousands of seconds. The lightcurve variability is also evident on larger timescales of days, weeks and months, with a dynamic range of more than one order of magnitude in all four SFXTs. The X-ray spectra are typically hard, with an average 2-10 keV luminosity during this monitoring of about 1E33-1E34 erg/s. We detected pulsations from the pulsar AXJ1841.0-0536, with a period of 4.7008+/-0.0004 s. This monitoring demonstrates that these transients spend most of the time accreting matter, although at a much lower level (~100-1000 times lower than during the bright outbusts), and that the true quiescence, characterized by a soft spectrum and a luminosity of a few 1E32 erg/s, observed in the past only in a couple of members of this class, is probably a very rare state., Accepted for publication in ApJ

We present an optical/near-infrared search for a counterpart to the perplexing radio transient GCRT J1745-3009, a source located ~1 degree from the Galactic Center. Motivated by some similarities to radio bursts from nearby ultracool dwarfs, and by a distance upper limit of 70 pc for the emission to not violate the 1e12 K brightness temperature limit for incoherent radiation, we searched for a nearby star at the position of GCRT J1745-3009. We found only a single marginal candidate, limiting the presence of any late-type star to >1 kpc (spectral types earlier than M9), >200 pc (spectral types L and T0-T4), and >100 pc (spectral types T4-T7), thus severely restricting the possible local counterparts to GCRT J1745-3009. We also exclude any white dwarf within 1 kpc or a supergiant star out to the distance of the Galactic Center as possible counterparts. This implies that GCRT J1745-3009 likely requires a coherent emission process, although whether or not it reflects a new class of sources is unclear., 10 pages, 5 figures. Accepted for publication in the Astrophysical Journal

Millimeter-wave observations of PN, CP, and HCP have been carried out toward circumstellar envelopes of evolved stars using the Arizona Radio Observatory (ARO). HCP and PN have been identified in the carbon-rich source CRL 2688 via observations at 1 mm using the Submillimeter Telescope (SMT) and 2-3 mm with the Kitt Peak 12 m. An identical set of measurements were carried out toward IRC +10216, as well as observations of CP at 1 mm. PN was also observed toward VY Canis Majoris (VY CMa), an oxygen-rich supergiant star. The PN and HCP line profiles in CRL 2688 and IRC +10216 are roughly flat topped, indicating unresolved, optically thin emission; CP, in contrast, has a distinct "U" shape in IRC +10216. Modeling of the line profiles suggests abundances, relative to H2, of f(PN) ~ (3–5) × 10−9 and f(HCP) ~ 2 × 10−7 in CRL 2688, about an order of magnitude higher than in IRC +10216. In VY CMa, f(PN) is ~4 × 10−8. The data in CRL 2688 and IRC +10216 are consistent with LTE formation of HCP and PN in the inner envelope, as predicted by theoretical calculations, with CP a photodissociation product at larger radii. The observed abundance of PN in VY CMa is a factor of 100 higher than LTE predictions. In IRC +10216, the chemistry of HCP/CP mimics that of HCN/CN and suggests an N2 abundance of f ~ 1 × 10−7. The chemistry of phosphorus appears active in both carbon- and oxygen-rich envelopes of evolved stars.

In this paper we survey the theory of wind accretion in high mass X-ray binaries hosting a magnetic neutron star and a supergiant companion. We concentrate on the different types of interaction between the inflowing wind matter and the neutron star magnetosphere that are relevant when accretion of matter onto the neutron star surface is largely inhibited; these include the inhibition through the centrifugal and magnetic barriers. Expanding on earlier work, we calculate the expected luminosity for each regime and derive the conditions under which transition from one regime to another can take place. We show that very large luminosity swings (~10^4 or more on time scales as short as hours) can result from transitions across different regimes. The activity displayed by supergiant fast X-ray transients, a recently discovered class of high mass X-ray binaries in our galaxy, has often been interpreted in terms of direct accretion onto a neutron star immersed in an extremely clumpy stellar wind. We show here that the transitions across the magnetic and/or centrifugal barriers can explain the variability properties of these sources as a results of relatively modest variations in the stellar wind velocity and/or density. According to this interpretation we expect that supergiant fast X-ray transients which display very large luminosity swings and host a slowly spinning neutron star are characterized by magnetar-like fields, irrespective of whether the magnetic or the centrifugal barrier applies. Supergiant fast X-ray transients might thus provide a new opportunity to detect and study magnetars in binary systems., Accepted for publication in ApJ. 16 pages, 6 figures

We report the first Very Long Baseline Array (VLBA) observations of 43 GHz v=1, J=1-0 SiO masers in the circumstellar envelope of the M-type semi-regular supergiant variable star AH Sco at 2 epochs separated by 12 days in March 2004. These high-resolution VLBA images reveal that the distribution of SiO masers is roughly on a persistent elliptical ring with the lengths of the major and minor axes of about 18.5 and 15.8 mas, respectively, along a position angle of 150^{\circ}. And the red-shifted masers are found to be slightly closer to the central star than the blue-shifted masers. The line-of-sight velocity structure of the SiO masers shows that with respect to the systemic velocity of -6.8 km/s the higher velocity features are closer to the star, which can be well explained by the simple outflow or infall without rotation kinematics of SiO masers around AH Sco. Study of proper motions of 59 matched features between two epochs clearly indicates that the SiO maser shell around AH Sco was undergoing an overall contraction to the star at a velocity of ~13 km/s at a distance of 2.26 kpc to AH Sco. Our 3-dimensional maser kinematics model further suggests that such an inward motion is very likely due to the gravitation of the central star. The distance to AH Sco of 2.26$\pm$0.19 kpc obtained from the 3-dimensional kinematics model fitting is consistent with its kinematic distance of 2.0 kpc., Comment: 23 pages, 4 figures, 4 tables, accepted for publication in ApJ

We present the results of a search for the JK = 10→ 00 transition of PH3 (phosphine) at 267 GHz toward several circumstellar envelopes using the Arizona Radio Observatory 10 m Submillimeter Telescope (SMT). In the carbon-rich shells of IRC +10216 and CRL 2688, we have detected emission lines exactly at the PH3 frequency. Toward the oxygen-rich supergiant VY Canis Majoris, only an upper limit was obtained, while in the evolved carbon-rich proto-planetary nebula CRL 618, the transition is contaminated by vibrationally excited HC3N (ν7 = 4). The line shape in IRC +10216 appears to consist of two distinct components: a flat-topped profile with a width of ~28 km s−1, as is typical for this source, and a narrower feature approximately 4 km s−1 wide. The narrow component likely arises from the inner envelope (r < 8R*) where the gas has not reached the terminal expansion velocity, or it is nonthermal emission. Based on the broader component, the abundance of PH3 with respect to H2 is estimated to be 5 × 10−8 in a region with a radius of r < 150R*. If the narrower component is thermal, it implies a phosphine abundance of ~5 × 10−7 close to the stellar photosphere (r < 8R*). In CRL 2688, the PH3 abundance is less constrained, with plausible values ranging from 3 × 10−8 to 4 × 10−7, assuming a spherical distribution. Phosphine appears to be present in large concentrations in the inner envelope of C-rich AGB stars, and thus may function as a parent molecule for other phosphorus species.