First orbital solution for the non-thermal emitter Cyg OB2 #9

After the first detection of its binary nature, the spectroscopic monitoring of the non-thermal radio emitter Cyg OB2 #9 (P=2.4yrs) has continued, doubling the number of available spectra of the star. Since the discovery paper of 2008, a second periastron passage has occurred in February 2009. Using a variety of techniques, the radial velocities could be estimated and a first, preliminary orbital solution was derived from the HeI5876 line. The mass ratio appears close to unity and the eccentricity is large, 0.7--0.75. X-ray data from 2004 and 2007 are also analyzed in quest of peculiarities linked to binarity. The observations reveal no large overluminosity nor strong hardness, but it must be noted that the high-energy data were taken after the periastron passage, at a time where colliding wind emission may be low. Some unusual X-ray variability is however detected, with a 10% flux decrease between 2004 and 2007. To clarify their origin and find a more obvious signature of the wind-wind collision, additional data, taken at periastron and close to it, are needed.Comment: 15 pages, 4 figures, accepted by Ap

A Radial Velocity Survey of the Cygnus OB2 Association

We conducted a radial velocity survey of the Cygnus OB2 Association over a 6 year (1999 - 2005) time interval to search for massive close binaries. During this time we obtained 1139 spectra on 146 OB stars to measure mean systemic radial velocities and radial velocity variations. We spectroscopically identify 73 new OB stars for the first time, the majority of which are likely to be Association members. Spectroscopic evidence is also presented for a B3Iae classification and temperature class variation (B3 - B8) on the order of 1 year for Cygnus OB2 No. 12. Calculations of the intial mass function with the current spectroscopic sample yield Gamma = -2.2 +/- 0.1. Of the 120 stars with the most reliable data, 36 are probable and 9 are possible single-lined spectroscopic binaries. We also identify 3 new and 8 candidate double-lined spectroscopic binaries. These data imply a lower limit on the massive binary fraction of 30% - 42%. The calculated velocity dispersion for Cygnus OB2 is 2.44 +/- km/s, which is typical of open clusters. No runaway OB stars were found.Comment: 56 pages, 23 figures, 5 tables, accepted for publication in the Astrophysical Journa

Observed Consequences of Presupernova Instability in Very Massive Stars

This chapter concentrates on the deaths of very massive stars, the events leading up to their deaths, and how mass loss affects the resulting death. The previous three chapters emphasized the theory of wind mass loss, eruptions, and core collapse physics, but here we emphasize mainly the observational properties of the resulting death throes. Mass loss through winds, eruptions, and interacting binaries largely determines the wide variety of different types of supernovae that are observed, as well as the circumstellar environments into which the supernova blast waves expand. Connecting these observed properties of the explosions to the initial masses of their progenitor stars is, however, an enduring challenge and is especially difficult for very massive stars. Superluminous supernovae, pair instability supernovae, gamma ray bursts, and "failed" supernovae are all end fates that have been proposed for very massive stars, but the range of initial masses or other conditions leading to each of these (if they actually occur) are still very certain. Extrapolating to infer the role of very massive stars in the early universe is essentially unencumbered by observational constraints and still quite dicey.Comment: 39 pages, 5 figures, to appear as chapter in the book "Very Massive Stars in the Local Universe", ed. J. Vin

A high angular resolution survey of massive stars in Cygnus OB2 : JHK adaptive optics results from the Gemini Near-Infrared Imager

We present results of a high angular resolution survey of massive OB stars in the Cygnus OB2 association that we conducted with the NIRI camera and ALTAIR adaptive optics system of the Gemini North telescope. We observed 74 O- and early B-type stars in Cyg OB2 in the JHK infrared bands in order to detect binary and multiple companions. The observations are sensitive to equal-brightness pairs at separations as small as 0.08 \arcsec, and progressively fainter companions are detectable out to Δ K = 9 mag at a separation of 2 arcsec. This faint contrast limit due to readnoise continues out to 10 arcsec near the edge of the detector. We assigned a simple probability of chance alignment to each companion based upon its separation and magnitude difference from the central target star and upon areal star counts for the general star field of Cyg OB2. Companion stars with a field membership probability of less than 1% are assumed to be physical companions. This assessment indicates that 47% of the targets have at least one resolved companion that is probably gravitationally bound. Including known spectroscopic binaries, our sample includes 27 binary, 12 triple, and 9 systems with four or more components. These results confirm studies of high mass stars in other environments that find that massive stars are born with a high multiplicity fraction. The results are important for the placement of the stars in the H-R diagram, the interpretation of their spectroscopic analyses, and for future mass determinations through measurement of orbital motion

Variable millimetre radiation from the colliding-wind binary Cygnus OB2 #8A

Context. Massive binaries have stellar winds that collide. In the colliding-wind region, various physically interesting processes occur, leading to enhanced X-ray emission, non-thermal radio emission, as well as non-thermal X-rays and gamma-rays. Non-thermal radio emission (due to synchrotron radiation) has so far been observed at centimetre wavelengths. At millimetre wavelengths, the stellar winds and the colliding-wind region emit more thermal free-free radiation, and it is expected that any non-thermal contribution will be difficult or impossible to detect. Aims. We aim to determine if the material in the colliding-wind region contributes substantially to the observed millimetre fluxes of a colliding-wind binary. We also try to distinguish the synchrotron emission from the free-free emission. Methods. We monitored the massive binary Cyg OB2 #8A at 3 mm with the NOrthern Extended Millimeter Array (NOEMA) interferometer of the Institut de Radioastronomie Millimétrique (IRAM). The data were collected in 14 separate observing runs (in 2014 and 2016), and provide good coverage of the orbital period. Results. The observed millimetre fluxes range between 1.1 and 2.3 mJy, and show phase-locked variability, clearly indicating that a large part of the emission is due to the colliding-wind region. A simple synchrotron model gives fluxes with the correct order of magnitude, but with a maximum that is phase-shifted with respect to the observations. Qualitatively this phase shift can be explained by our neglect of orbital motion on the shape of the colliding-wind region. A model using only free-free emission results in only a slightly worse explanation of the observations. Additionally, on the map of our observations we also detect the O6.5 III star Cyg OB2 #8B, for which we determine a 3 mm flux of 0.21 ± 0.033 mJy. Conclusions. The question of whether synchrotron radiation or free-free emission dominates the millimetre fluxes of Cyg OB2 #8A remains open. More detailed modelling of this system, based on solving the hydrodynamical equations, is required to give a definite answer

The Tarantula Massive Binary Monitoring: I. Observational campaign and OB-type spectroscopic binaries

© ESO, 2017.Context. Massive binaries play a crucial role in the Universe. Knowing the distributions of their orbital parameters is important for a wide range of topics from stellar feedback to binary evolution channels and from the distribution of supernova types to gravitational wave progenitors, yet no direct measurements exist outside the Milky Way. Aims. The Tarantula Massive Binary Monitoring project was designed to help fill this gap by obtaining multi-epoch radial velocity (RV) monitoring of 102 massive binaries in the 30 Doradus region. Methods. In this paper we analyze 32 FLAMES/GIRAFFE observations of 93 O- and 7 B-type binaries. We performed a Fourier analysis and obtained orbital solutions for 82 systems: 51 single-lined (SB1) and 31 double-lined (SB2) spectroscopic binaries. Results. Overall, the binary fraction and orbital properties across the 30 Doradus region are found to be similar to existing Galactic samples. This indicates that within these domains environmental effects are of second order in shaping the properties of massive binary systems. A small difference is found in the distribution of orbital periods, which is slightly flatter (in log space) in 30 Doradus than in the Galaxy, although this may be compatible within error estimates and differences in the fitting methodology. Also, orbital periods in 30 Doradus can be as short as 1.1 d, somewhat shorter than seen in Galactic samples. Equal mass binaries (q> 0.95) in 30 Doradus are all found outside NGC 2070, the central association that surrounds R136a, the very young and massive cluster at 30 Doradus's core. Most of the differences, albeit small, are compatible with expectations from binary evolution. One outstanding exception, however, is the fact that earlier spectral types (O2-O7) tend to have shorter orbital periods than later spectral types (O9.2-O9.7). Conclusions. Our results point to a relative universality of the incidence rate of massive binaries and their orbital properties in the metallicity range from solar (Z) to about half solar. This provides the first direct constraints on massive binary properties in massive star-forming galaxies at the Universe's peak of star formation at redshifts z ~ 1 to 2 which are estimated to have Z 0.5 Z

Dissecting a supernova impostor's circumstellar medium: MUSEing about the SHAPE of eta Carinae's outer ejecta

Aims. The role of episodic mass loss is one of the outstanding questions in massive star evolution. The structural inhomogeneities and kinematics of their nebulae are tracers of their mass-loss history. We conduct a three-dimensional morpho-kinematic analysis of the ejecta of η Car outside its famous Homunculus nebula. Methods. We carried out the first large-scale integral field unit observations of η Car in the optical, covering a field of view of 1′× 1′ centered on the star. Observations with the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope (VLT) reveal the detailed three-dimensional structure of η Car’s outer ejecta. Morpho-kinematic modeling of these ejecta is conducted with the code SHAPE. Results. The largest coherent structure in η Car’s outer ejecta can be described as a bent cylinder with roughly the same symmetry axis as the Homunculus nebula. This large outer shell is interacting with the surrounding medium, creating soft X-ray emission. Doppler velocities of up to 3000 km s-1 are observed. We establish the shape and extent of the ghost shell in front of the southern Homunculus lobe and confirm that the NN condensation can best be modeled as a bowshock in the orbital/equatorial plane. Conclusions. The SHAPE modeling of the MUSE observations provides a significant gain in the study of the three-dimensional structure of η Car’s outer ejecta. Our SHAPE modeling indicates that the kinematics of the outer ejecta measured with MUSE can be described by a spatially coherent structure, and that this structure also correlates with the extended soft X-ray emission associated with the outer debris field. The ghost shell immediately outside the southern Homunculus lobe hints at a sequence of eruptions within the time frame of the Great Eruption from 1837–1858 or possibly a later shock/reverse shock velocity separation. Our 3D morpho-kinematic modeling and the MUSE observations constitute an invaluable dataset to be confronted with future radiation-hydrodynamics simulations. Such a comparison may shed light on the yet elusive physical mechanism responsible for η Car-like eruptions

The relation between the most-massive star and its parental star cluster mass

We present a thorough literature study of the most-massive star, m_max, in several young star clusters in order to assess whether or not star clusters are populated from the stellar initial mass function (IMF) by random sampling over the mass range 0.01 < m < 150 M_sol without being constrained by the cluster mass, M_ecl. The data reveal a partition of the sample into lowest mass objects (M_ecl < 10^2 M_sol), moderate mass clusters (10^2 M_sol < M_ecl < 10^3 M_sol) and rich clusters above 10^3 M_sol. Additionally, there is a plateau of a constant maximal star mass (m_max ~ 25 M_sol) for clusters with masses between 10^3 M_sol and 4 10^3 M_sol. Statistical tests of this data set reveal that the hypothesis of random sampling from the IMF between 0.01 and 150 M_sol is highly unlikely for star clusters more massive than 10^2 M_sol with a probability of p ~ 2 10^-7 for the objects with M_ecl between 10^2 M_sol and 10^3 M_sol and p ~ 3 10^-9 for the more massive star clusters. Also, the spread of m_max values at a given M_ecl is smaller than expected from random sampling. We suggest that the basic physical process able to explain this dependence of stellar inventory of a star cluster on its mass may be the interplay between stellar feedback and the binding energy of the cluster-forming molecular cloud core. Given these results, it would follow that an integrated galactic initial mass function (IGIMF) sampled from such clusters would automatically be steeper in comparison to the IMF within individual star clusters.Comment: 26 pages, 13 figures, 3 tables, accepted for publication in MNRA