A VLT/FLAMES survey for massive binaries in Westerlund 1 III. The WC9d binary W239 and implications for massive stellar evolution

Context. There is growing evidence that a treatment of binarity amongst OB stars is essential for a full theory of stellar evolution. However the binary properties of massive stars - frequency, mass ratio & orbital separation - are still poorly constrained. Aims. In order to address this shortcoming we have undertaken a multiepoch spectroscopic study of the stellar population of the young massive cluster Westerlund 1. In this paper we present an investigation into the nature of the dusty Wolf-Rayet star and candidate binary W239. Methods. To accomplish this we have utilised our spectroscopic data in conjunction with multi-year optical and near-IR photometric observations in order to search for binary signatures. Comparison of these data to synthetic non-LTE model atmosphere spectra were used to derive the fundamental properties of the WC9 primary. Results. We found W239 to have an orbital period of only similar to 5.05 days, making it one of the most compact WC binaries yet identified. Analysis of the long term near-IR lightcurve reveals a significant flare between 2004-6. We interpret this as evidence for a third massive stellar component in the system in a long period (> 6 yr), eccentric orbit, with dust production occuring at periastron leading to the flare. The presence of a near-IR excess characteristic of hot (similar to 1300 K) dust at every epoch is consistent with the expectation that the subset of persistent dust forming WC stars are short (< 1 yr) period binaries, although confirmation will require further observations. Non-LTE model atmosphere analysis of the spectrum reveals the physical properties of the WC9 component to be fully consistent with other Galactic examples. Conclusions. The simultaneous presence of both short period Wolf-Rayet binaries and cool hypergiants within Wd 1 provides compelling evidence for a bifurcation in the post-Main Sequence evolution of massive stars due to binarity. Short period O+OB binaries will evolve directly to the Wolf-Rayet phase, either due to an episode of binary mediated mass loss - likely via case A mass transfer or a contact configuration - or via chemically homogenous evolution. Conversely, long period binaries and single stars will instead undergo a red loop across the HR diagram via a cool hypergiant phase. Future analysis of the full spectroscopic dataset for Wd 1 will constrain the proportion of massive stars experiencing each pathway; hence quantifying the importance of binarity in massive stellar evolution up to and beyond supernova and the resultant production of relativistic remnants.

Search for very high energy gamma-ray emission from pulsar-pulsar wind nebula systems with the MAGIC telescope

The MAGIC collaboration has searched for high-energy gamma-ray emission of some of the most promising pulsar candidates above an energy threshold of 50 GeV, an energy not reachable up to now by other ground-based instruments. Neither pulsed nor steady gamma-ray emission has been observed at energies of 100 GeV from the classical radio pulsars PSR J0205+6449 and PSR J2229+6114 (and their nebulae 3C58 and Boomerang, respectively) and the millisecond pulsar PSR J0218+4232. Here, we present the flux upper limits for these sources and discuss their implications in the context of current model predictions.

The embedded cluster DBSB 48 in the nebula Hoffleit 18: Comparison with Trumpler 14

We derive fundamental parameters of the embedded cluster DBSB 48 in the southern nebula Hoffleit 18 and the very young open cluster Trumpler 14, by means of deep JHK(s) infrared photometry. We build colour-magnitude and colour-colour diagrams to derive reddening and age, based on main sequence and pre-main sequence distributions. Radial stellar density profiles are used to study cluster structure and guide photometric diagram extractions. Field-star decontamination is applied to uncover the intrinsic cluster sequences in the diagrams. Ages are inferred from K-excess fractions. A prominent pre-main sequence population is present in DBSB 48, and the K-excess fraction f(K) = 55 +/- 6% gives an age of 1.1 +/- 0.5 Myr. A mean reddening of A(Ks) = 0.9 +/- 0.03 was found, corresponding to A(v) = 8.2 +/- 0.3. The cluster CMD is consistent with the far kinematic distance of 5 kpc for Hoffleit 18. For Trumpler 14 we derived similar parameters as in previous studies in the optical, in particular an age of 1.7 +/- 0.7 Myr. The fraction of stars with infrared excess in Trumpler 14 is f(K) = 28 +/- 4%. Despite the young ages, both clusters are described by a King profile with core radii R-core = 0.46 +/- 0.05 pc and R-core = 0.35 +/- 0.04 pc, respectively, for DBSB 48 and Trumpler 14. Such cores are smaller than those of typical open clusters. Small cores are probably related to the cluster formation and/or parent molecular cloud fragmentation. In DBSB 48, the magnitude extent of the upper main sequence is Delta K-s approximate to 2 mag, while in Trumpler 14 it is Delta K-s approximate to 5 mag, consistent with the estimated ages. (c) 2008 Elsevier B.V. All rights reserved.

Near-infrared integral field spectroscopy of the Homunculus nebula around eta Carinae using Gemini/CIRPASS

This work presents the first integral field spectroscopy of the Homunculus nebula around eta Carinae in the near-infrared spectral region (J band). We confirmed the presence of a hole on the polar region of each lobe, as indicated by previous near-IR long-slit spectra and mid-IR images. The holes can be described as a cylinder of height (i.e. the thickness of the lobe) and diameter of 6.5 and 6.0 x 10(16) cm, respectively. We also mapped the blue-shifted component of He I lambda 10830 seen towards the NW lobe. Contrary to previous works, we suggested that this blue-shifted component is not related to the Paddle but it is indeed in the equatorial disc. We confirmed the claim of N. Smith and showed that the spatial extent of the Little Homunculus matches remarkably well the radio continuum emission at 3 cm, indicating that the Little Homunculus can be regarded as a small H II region. Therefore, we used the optically thin 1.3 mm radio flux to derive a lower limit for the number of Lyman-continuum photons of the central source in eta Car. In the context of a binary system, and assuming that the ionizing flux comes entirely from the hot companion star, the lower limit for its spectral type and luminosity class ranges from O5.5 III to O7 I. Moreover, we showed that the radio peak at 1.7 arcsec NW from the central star is in the same line-of-sight of the `Sr-filament` but they are obviously spatially separated, while the blue-shifted component of He I lambda 10830 may be related to the radio peak and can be explained by the ultraviolet radiation from the companion star.

Modelling the warm H(2) infrared emission of the Helix nebula cometary knots

Molecular hydrogen emission is commonly observed in planetary nebulae. Images taken in infrared H(2) emission lines show that at least part of the molecular emission is produced inside the ionized region. In the best studied case, the Helix nebula, the H(2) emission is produced inside cometary knots (CKs), comet-shaped structures believed to be clumps of dense neutral gas embedded within the ionized gas. Most of the H(2) emission of the CKs seems to be produced in a thin layer between the ionized diffuse gas and the neutral material of the knot, in a mini-photodissociation region (mini-PDR). However, PDR models published so far cannot fully explain all the characteristics of the H(2) emission of the CKs. In this work, we use the photoionization code AANGABA to study the H(2) emission of the CKs, particularly that produced in the interface H(+)/H(0) of the knot, where a significant fraction of the H(2) 1-0 S(1) emission seems to be produced. Our results show that the production of molecular hydrogen in such a region may explain several characteristics of the observed emission, particularly the high excitation temperature of the H(2) infrared lines. We find that the temperature derived from H(2) observations, even of a single knot, will depend very strongly on the observed transitions, with much higher temperatures derived from excited levels. We also proposed that the separation between the H alpha and [N II] peak emission observed in the images of CKs may be an effect of the distance of the knot from the star, since for knots farther from the central star the [N II] line is produced closer to the border of the CK than H alpha.

The peculiar isolated neutron star in the Carina Nebula Deep XMM-Newton and ESO-VLT observations of 2XMM J104608.7-594306

While fewer in number than the dominant rotation-powered radio pulsar population, peculiar classes of isolated neutron stars (INSs) which include magnetars, the ROSAT-discovered "Magnificent Seven" (M7), rotating radio transients (RRATs), and central compact objects in supernova remnants (CCOs) - represent a key element in understanding the neutron star phenomenology. We report the results of an observational campaign to study the properties of the source 2XMM J104608.7-594306, a newly discovered thermally emitting INS. The evolutionary state of the neutron star is investigated by means of deep dedicated observations obtained with the XMM-Newton Observatory, the ESO Very Large Telescope, as well as publicly available gamma-ray data from the Fermi Space Telescope and the AGILE Mission. The observations confirm previous expectations and reveal a unique type of object. The source, which is likely within the Carina Nebula (N-H = 2.6x10(21) cm(-2)), has a spectrum that is both thermal and soft, with kT(infinity) = 135 eV. Non-thermal (magnetospheric) emission is not detected down to 1% (3 sigma, 0.1-12 keV) of the source luminosity. Significant deviations (absorption features) from a simple blackbody model are identified in the spectrum of the source around energies 0.6 keV and 1.35 keV. While the former deviation is likely related to a local oxygen overabundance in the Carina Nebula, the latter can only be accounted for by an additional spectral component, which is modelled as a Gaussian line in absorption with EW = 91 eV and sigma = 0.14 keV (1 sigma). Furthermore, the optical counterpart is fainter than m(V) = 27 (2 sigma) and no gamma-ray emission is significantly detected by either the Fermi or AGILE missions. Very interestingly, while these characteristics are remarkably similar to those of the M7 or the only RRAT so far detected in X-rays, which all have spin periods of a few seconds, we found intriguing evidence of very rapid rotation, P = 18.6ms, at the 4 sigma confidence level. We interpret these new results in the light of the observed properties of the currently known neutron star population, in particular those of standard rotation-powered pulsars, recycled objects, and CCOs. We find that none of these scenarios can satisfactorily explain the collective properties of 2XMM J104608.7-594306, although it may be related to the still poorly known class of Galactic anti-magnetars. Future XMM-Newton data, granted for the next cycle of observations (AO11), will help us to improve our current observational interpretation of the source, enabling us to significantly constrain the rate of pulsar spin down.

From stellar nebula to planets: The refractory components

Context. To date, calculations of planet formation have mainly focused on dynamics, and only a few have considered the chemical composition of refractory elements and compounds in the planetary bodies. While many studies have been concentrating on the chemical composition of volatile compounds (such as H2O, CO, CO2) incorporated in planets, only a few have considered the refractory materials as well, although they are of great importance for the formation of rocky planets. Aims. We computed the abundance of refractory elements in planetary bodies formed in stellar systems with a solar chemical composition by combining models of chemical composition and planet formation. We also considered the formation of refractory organic compounds, which have been ignored in previous studies on this topic. Methods. We used the commercial software package HSC Chemistry to compute the condensation sequence and chemical composition of refractory minerals incorporated into planets. The problem of refractory organic material is approached with two distinct model calculations: the first considers that the fraction of atoms used in the formation of organic compounds is removed from the system (i.e., organic compounds are formed in the gas phase and are non-reactive); and the second assumes that organic compounds are formed by the reaction between different compounds that had previously condensed from the gas phase. Results. Results show that refractory material represents more than 50 wt % of the mass of solids accreted by the simulated planets with up to 30 wt % of the total mass composed of refractory organic compounds. Carbide and silicate abundances are consistent with C/O and Mg/Si elemental ratios of 0.5 and 1.02 for the Sun. Less than 1 wt % of carbides are present in the planets, and pyroxene and olivine are formed in similar quantities. The model predicts planets that are similar in composition to those of the solar system. Starting from a common initial nebula composition, it also shows that a wide variety of chemically different planets can form, which means that the differences in planetary compositions are due to differences in the planetary formation process. Conclusions. We show that a model in which refractory organic material is absent from the system is more compatible with observations. The use of a planet formation model is essential to form a wide diversity of planets in a consistent way.

From stellar nebula to planetesimals

Context. Solar and extrasolar comets and extrasolar planets are the subject of numerous studies in order to determine their chemical composition and internal structure. In the case of planetesimals, their compositions are important as they govern in part the composition of future planets. Aims. The present works aims at determining the chemical composition of icy planetesimals, believed to be similar to present day comets, formed in stellar systems of solar chemical composition. The main objective of this work is to provide valuable theoretical data on chemical composition for models of planetesimals and comets, and models of planet formation and evolution. Methods. We have developed a model that calculates the composition of ices formed during the cooling of the stellar nebula. Coupled with a model of refractory element formation, it allows us to determine the chemical composition and mass ratio of ices to rocks in icy planetesimals throughout in the protoplanetary disc. Results. We provide relationships for ice line positions (for different volatile species) in the disc, and chemical compositions and mass ratios of ice relative to rock for icy planetesimals in stellar systems of solar chemical composition. From an initial homogeneous composition of the nebula, a wide variety of chemical compositions of planetesimals were produced as a function of the mass of the disc and distance to the star. Ices incorporated in planetesimals are mainly composed of H2O, CO, CO2, CH3OH, and NH3. The ice/rock mass ratio is equal to 1 ± 0.5 in icy planetesimals following assumptions. This last value is in good agreement with observations of solar system comets, but remains lower than usual assumptions made in planet formation models, taking this ratio to be of 2–3.

Rotational velocities of single and binary O-type stars in the Tarantula Nebula

Rotation is a key parameter in the evolution of massive stars, affecting their evolution, chemical yields, ionizing photon budget, and final fate. We determined the projected rotational velocity, υ e sin i, of ~330 O-type objects, i.e. ~210 spectroscopic single stars and ~110 primaries in binary systems, in the Tarantula nebula or 30 Doradus (30 Dor) region. The observations were taken using VLT/FLAMES and constitute the largest homogeneous dataset of multi-epoch spectroscopy of O-type stars currently available. The most distinctive feature of the υ e sin i distributions of the presumed-single stars and primaries in 30 Dor is a low-velocity peak at around 100 km s^-1. Stellar winds are not expected to have spun-down the bulk of the stars significantly since their arrival on the main sequence and therefore the peak in the single star sample is likely to represent the outcome of the formation process. Whereas the spin distribution of presumed-single stars shows a well developed tail of stars rotating more rapidly than 300 km s^-1, the sample of primaries does not feature such a high-velocity tail. The tail of the presumed-single star distribution is attributed for the most part - and could potentially be completely due - to spun-up binary products that appear as single stars or that have merged. This would be consistent with the lack of such post-interaction products in the binary sample, that is expected to be dominated by pre-interaction systems. The peak in this distribution is broader and is shifted toward somewhat higher spin rates compared to the distribution of presumed-single stars. Systems displaying large radial velocity variations, typical for short period systems, appear mostly responsible for these differences.

Search for VHE gamma-ray emission from Geminga pulsar and nebula with the MAGIC telescopes

The Geminga pulsar, one of the brighest gamma-ray sources, is a promising candidate for emission of very-high-energy (VHE > 100 GeV) pulsed gamma rays. Also, detection of a large nebula have been claimed by water Cherenkov instruments. We performed deep observations of Geminga with the MAGIC telescopes, yielding 63 hours of good-quality data, and searched for emission from the pulsar and pulsar wind nebula. We did not find any significant detection, and derived 95% confidence level upper limits. The resulting upper limits of 5.3 × 10^(−13) TeV cm^(−2)s^(−1) for the Geminga pulsar and 3.5 × 10^(−12) TeV cm^(−2)s^(−1) for the surrounding nebula at 50 GeV are the most constraining ones obtained so far at VHE. To complement the VHE observations, we also analyzed 5 years of Fermi-LAT data from Geminga, finding that the sub-exponential cut-off is preferred over the exponential cut-off that has been typically used in the literature. We also find that, above 10 GeV, the gamma-ray spectra from Geminga can be described with a power law with index softer than 5. The extrapolation of the power-law Fermi-LAT pulsed spectra to VHE goes well below the MAGIC upper limits, indicating that the detection of pulsed emission from Geminga with the current generation of Cherenkov telescopes is very difficult.

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

Context. The analysis and interpretation of the H(2) line emission from planetary nebulae have been done in the literature by assuming that the molecule survives only in regions where the hydrogen is neutral, as in photodissociation, neutral clumps, or shocked regions. However, there is strong observational and theoretical evidence that at least part of the H(2) emission is produced inside the ionized region of these objects. Aims. The aim of the present work is to calculate and analyze the infrared line emission of H(2) produced inside the ionized region of planetary nebulae using a one-dimensional photoionization code. Methods. The photoionization code Aangaba was improved in order to calculate the statistical population of the H(2) energy levels, as well as the intensity of the H(2) infrared emission lines in the physical conditions typical of planetary nebulae. A grid of models was obtained and the results then analyzed and compared with the observational data. Results. We show that the contribution of the ionized region to the H(2) line emission can be important, particularly in the case of nebulae with high-temperature central stars. This result explains why H(2) emission is more frequently observed in bipolar planetary nebulae (Gatley's rule), since this kind of object typically has hotter stars. Collisional excitation plays an important role in populating the rovibrational levels of the electronic ground state of H(2) molecules. Radiative mechanisms are also important, particularly for the upper vibrational levels. Formation pumping can have minor effects on the line intensities produced by de-excitation from very high rotational levels, especially in dense and dusty environments. We included the effect of the H(2) molecule on the thermal equilibrium of the gas, concluding that, in the ionized region, H(2) only contributes to the thermal equilibrium in the case of a very high temperature of the central star or a high dust-to-gas ratio, mainly through collisional de-excitation.

Distinguishing post-AGB impostors in a sample of pre-main sequence stars

Context. A sample of 27 sources, cataloged as pre-main sequence stars by the Pico dos Dias Survey (PDS), is analyzed to investigate a possible contamination by post-AGB stars. The far-infrared excess due to dust present in the circumstellar envelope is typical of both categories: young stars and objects that have already left the main sequence and are suffering severe mass loss. Aims. The two known post-AGB stars in our sample inspired us to seek for other very likely or possible post-AGB objects among PDS sources previously suggested to be Herbig Ae/Be stars, by revisiting the observational database of this sample. Methods. In a comparative study with well known post-AGBs, several characteristics were evaluated: (i) parameters related to the circumstellar emission; (ii) spatial distribution to verify the background contribution from dark clouds; (iii) spectral features; and (iv) optical and infrared colors. Results. These characteristics suggest that seven objects of the studied sample are very likely post-AGBs, five are possible post-AGBs, eight are unlikely post-AGBs, and the nature of seven objects remains unclear.

Detection of high-velocity material from the wind-wind collision zone of Eta Carinae across the 2009.0 periastron passage

We report near-infrared spectroscopic observations of the Eta Carinae massive binary system during 2008-2009 using the CRIRES spectrograph mounted on the 8m UT 1 Very Large Telescope (VLT Antu). We detect a strong, broad absorption wing in He I lambda 10833 extending up to -1900 km s(-1) across the 2009.0 spectroscopic event. Analysis of archival Hubble Space Telescope/Space Telescope Imaging Spectrograph ultraviolet and optical data identifies a similar high-velocity absorption (up to -2100 km s(-1)) in the ultraviolet resonance lines of Si IV lambda lambda 1394, 1403 across the 2003.5 event. Ultraviolet resonance lines from low-ionization species, such as Si II lambda lambda 1527, 1533 and CII lambda lambda 1334, 1335, show absorption only up to -1200 km s(-1), indicating that the absorption with velocities -1200 to -2100 km s(-1) originates in a region markedly more rapidly moving and more ionized than the nominal wind of the primary star. Seeing-limited observations obtained at the 1.6m OPD/LNA telescope during the last four spectroscopic cycles of Eta Carinae (1989-2009) also show high-velocity absorption in He I lambda 10833 during periastron. Based on the large OPD/LNA dataset, we determine that material with velocities more negative than -900 km s(-1) is present in the phase range 0.976 <= phi <= 1.023 of the spectroscopic cycle, but absent in spectra taken at phi <= 0.947 and phi >= 1.049. Therefore, we constrain the duration of the high-velocity absorption to be 95 to 206 days (or 0.047 to 0.102 in phase). We propose that the high-velocity absorption component originates in shocked gas in the wind-wind collision zone, at distances of 15 to 45 AU in the line-of-sight to the primary star. With the aid of three-dimensional hydrodynamical simulations of the wind-wind collision zone, we find that the dense high-velocity gas is along the line-of-sight to the primary star only if the binary system is oriented in the sky such that the companion is behind the primary star during periastron, corresponding to a longitude of periastron of omega similar to 240 degrees-270 degrees. We study a possible tilt of the orbital plane relative to the Homunculus equatorial plane and conclude that our data are broadly consistent with orbital inclinations in the range i = 40 degrees-60 degrees.

BONA FIDE, STRONG-VARIABLE GALACTIC LUMINOUS BLUE VARIABLE STARS ARE FAST ROTATORS: DETECTION OF A HIGH ROTATIONAL VELOCITY IN HR CARINAE

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 Si lambda lambda 4088-4116. To match their observed line profiles from 2009 May, a high rotational velocity of nu(rot) similar or equal to 150 +/- 20 km s(-1) is needed (assuming an inclination angle of 30 degrees), implying that HR Car rotates at similar or equal to 0.88 +/- 0.2 of its critical velocity for breakup (nu(crit)). 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 nu(crit) is reached. To the left of this strip, a forbidden zone with nu(rot)/nu(crit) > 1 is present, explaining why no LBVs are detected in this zone. Since dormant/ex LBVs like P Cygni and HD 168625 have low nu(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.