Source author record

L. Mahy

L. Mahy appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

34works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

34 published item(s)

preprint2022arXiv

Detecting Stripped Stars While Searching for Quiescent Black Holes

While the number of stellar-mass black holes detected in X-rays or as gravitational wave sources is steadily increasing, the known population remains orders of magnitude smaller than predicted by stellar evolution theory. A significant fraction of stellar-mass black holes is expected to hide in X-ray-quiet binaries where they are paired with a "normal" star. Although a handful of such quiescent black hole candidates have been proposed, the majority have been challenged by follow-up investigations. A confusion that emerged recently concerns binary systems that appear to contain a normal B-type star with an unseen companion, believed to be a black hole. On closer inspection, some of these seemingly normal B-type stars instead turn out to be stars stripped of most of their mass through an interaction with their binary companion, which in at least two cases is a rapidly rotating star rather than a compact object. These contaminants in the search for quiescent black holes are themselves extremely interesting objects as they represent a rare phase of binary evolution, and should be given special attention when searching for binaries hosting black holes in large spectroscopic studies.

preprint2022arXiv

Identifying quiescent compact objects in massive Galactic single-lined spectroscopic binaries

Aims. To search for these rare objects, we study 32 Galactic O-type stars that were reported as SB1s in the literature. In our sample we include Cyg X-1, which is known to host an accreting stellar-mass BH, and HD 74194, a supergiant fast X-ray transient, in order to validate our methodology. The final goal is to characterise the nature of the unseen companions to determine if they are MS stars, stripped helium stars, triples, or compact objects such as neutron stars or stellar-mass BHs. Methods. After measuring radial velocities and deriving orbital solutions for all the systems in our sample, we performed spectral disentangling to extract putative signatures of faint secondary companions from the composite spectra. We derived stellar parameters for the visible stars and estimated the mass ranges of the secondary stars using the binary mass function. Variability observed in the photometric TESS light curves was also searched for indications of the presence of putative companions, degenerate or not. Results. In 17 of the 32 systems reported as SB1s, we extract secondary signatures, down to mass ratios of ~0.15. For the 17 newly detected double-lined spectroscopic binaries (SB2s), we derive physical properties of the individual components and discuss why they have not been detected as such before. Among the remaining systems, we identify nine systems with possible NS or low-mass MS companions. For Cyg X-1 and HD 130298, we are not able to extract any signatures for the companions, and the minimum masses of their companions are estimated to be about 7Msun. Our simulations show that secondaries with such a mass should be detectable from our dataset, no matter their nature: MS stars, stripped helium stars or even triples. While this is expected for Cyg X-1, confirming our methodology, our simulations also strongly suggest that HD 130298 could be another candidate to host a stellar-mass BH.

preprint2022arXiv

The Gaia-ESO Survey: The analysis of the hot-star spectra

The Gaia-ESO Survey (GES) is a large public spectroscopic survey that has collected, over a period of 6 years, spectra of ~ 10^5 stars. This survey provides not only the reduced spectra, but also the stellar parameters and abundances resulting from the analysis of the spectra. The GES dataflow is organised in 19 working groups. Working group 13 (WG13) is responsible for the spectral analysis of the hottest stars (O, B and A type, with a formal cut-off of Teff > 7000 K) that were observed as part of GES. We present the procedures and techniques that have been applied to the reduced spectra, in order to determine the stellar parameters and abundances of these stars. The procedure used is similar to that of other working groups in GES. A number of groups (called `Nodes') each independently analyse the spectra, using their state-of-the-art techniques and codes. Specific for the analysis in WG13 is the large temperature range that is covered (Teff = 7000 - 50,000 K), requiring the use of different analysis codes. Most Nodes can therefore only handle part of the data. Quality checks are applied to the results of these Nodes by comparing them to benchmark stars, and by comparing them one to another. For each star the Node values are then homogenised into a single result: the recommended parameters and abundances. Eight Nodes each analysed (part of) the data. In total 17,693 spectra of 6462 stars were analysed, most of them in 37 open star clusters. The homogenisation led to stellar parameters for 5584 stars. Abundances were determined for a more limited number of stars. Elements studied are He, C, N, O, Ne, Mg, Al, Si and Sc. Abundances for at least one of those elements were determined for 292 stars. The hot-star data analysed here, as well as the Gaia-ESO Survey data in general, will be of considerable use in future studies of stellar evolution and open clusters.

preprint2022arXiv

WR 63: A multiple system (O+O)+WR ?

The spectrum of the Wolf-Rayet (WR) star WR 63 contains spectral lines of two different O stars that show regular radial velocity (RV) variations with amplitudes of ~160 and ~225 km/s on a ~4.03 d period. The light-curve shows two narrow eclipses that are 0.2 mag deep on the same period as the RV changes. On the other hand, our data show no significant RV variations for the WR spectral lines. Those findings are compatible with WR 63 being a triple system composed of two non-interacting late O stars orbiting a WR star on a period larger than 1000 days. The amplitude of the WR spectral line-profile variability reaches 7-8% of the line intensity and seems related to a 0.04 mag periodic photometric variation. Large wind density structures are a possible origin of this variability, but our data are not sufficient to verify this. Our analysis shows that, should the three stars be bound, they would be coeval with an age of about 5.9+/-1.4 Myrs. The distance to the O stars is estimated to be 3.4+/-0.5 kpc. Their dynamical masses are 14.3+/-0.1 and 10.3+/-0.1 M_sol. Using rotating, single star evolutionary tracks, we estimate their initial masses to be 18+/-2 and 16+/-2 M_sol for the primary and the secondary, respectively. Regular spectral monitoring is required in the future to detect RV variations of the WR star that would prove that it is gravitationally bound to the close O+OB system and to determine its mass.

preprint2021arXiv

Detailed models of interacting short-period massive binary stars

About a quarter of massive binary stars undergo mass transfer while both stars burn hydrogen at their cores, first on the thermal and then on the nuclear timescale. The nuclear timescale mass transfer leads to observable counterparts: the semi-detached so-called massive Algol binaries. However, comprehensive model predictions for these systems are sparse. We study them using a large grid of ~10,000 detailed binary evolution models calculated with the stellar evolution code MESA, covering initial donor masses between 10-40 M$_{\odot}$ and initial orbital periods above 1.4 d, at a metallicity suitable for the Large Magellanic Cloud (LMC). Our models imply ~30, or ~3% of the ~1,000 core hydrogen burning O-star binaries in the LMC to be currently in the semi-detached phase. Our donor models are up to 25-times more luminous than single stars of identical mass and effective temperature, which agrees with the observed Algols. A comparison of our models with the observed orbital periods and mass ratios implies rather conservative mass transfer in some systems, while very inefficient one in others. This is generally well reproduced by our spin-dependent mass transfer algorithm, except for the lowest considered masses. The observations reflect the slow increase of the surface nitrogen enrichment of the donors during the semi-detached phase all the way to CNO equilibrium. We also investigate the properties of our models after core hydrogen depletion of the donor star, when these models correspond to Wolf-Rayet/helium+OB star binaries. A dedicated spectroscopic survey of massive Algol systems may allow to derive the dependence of the efficiency of thermal timescale mass transfer on the binary parameters, as well as the efficiency of semiconvective mixing in the stellar interior. This would be a crucial step towards reliable binary models up to the formation of supernovae and compact objects.

preprint2021arXiv

The observed multiplicity properties of B-type stars in the Galactic young open cluster NGC 6231

It is well known that massive O-stars are frequently (if not always) found in binary or higher-order multiple systems, but this fact has been less robustly investigated for the lower mass range of the massive stars, represented by the B-type stars. Obtaining the binary fraction and orbital parameter distributions of B-type stars is crucial to understand the impact of multiplicity on the archetypal progenitor of core-collapse supernovae as well as to properly investigate formation channels for gravitational wave progenitors. This work aims to characterise the multiplicity of the B-star population of the young open cluster NGC 6231 through multi-epoch optical spectroscopy of 80 B-type stars. We analyse 31 FLAMES/GIRAFFE observations of 80 B-type stars, monitoring their radial velocities (RVs) and performing a least-squares spectral analysis (Lomb-Scargle) to search for periodicity in those stars with statistically significant variability in their RVs. We constrain an observed spectroscopic binary fraction of $33\pm5$% for the B-type stars of NGC 6231, with a first order bias-correction giving a true spectroscopic binary fraction of $52\pm8$%. Out of 27 B-type binary candidates, we obtained orbital solutions for 20 systems: 15 single-lined (SB1) and 5 double-lined spectroscopic binaries (SB2s). We present these orbital solutions and the orbital parameter distributions associated with them. Our results indicate that Galactic B-type stars are less frequently found in binary systems than their more massive O-type counterparts, but their orbital properties generally resemble those of B- and O-type stars in both the Galaxy and Large Magellanic Cloud.

preprint2020arXiv

Is HR 6819 a triple system containing a black hole? -- An alternative explanation

HR 6819 was recently proposed to be a triple system consisting of an inner B-type giant + black hole binary with an orbital period of 40d and an outer Be tertiary. This interpretation is mainly based on two inferences: that the emission attributed to the outer Be star is stationary, and that the inner star, which is used as mass calibrator for the black hole, is a B-type giant. We re-investigate the properties of HR 6819 by spectral disentangling and an atmosphere analysis of the disentangled spectra to search for a possibly simpler alternative explanation for HR 6819. Disentangling implies that the Be component is not a static tertiary, but rather a component of the binary in the 40-d orbit. The inferred radial velocity amplitudes imply an extreme mass ratio of M_2/M_1 = 15 +/- 3. We infer spectroscopic masses of 0.4$^{+0.3}_{-0.1}$ Msun and 6$^{+5}_{-3}$ Msun for the primary and secondary, which agree well with the dynamical masses for an inclination of i = 32 deg. This indicates that the primary might be a stripped star rather than a B-type giant. Evolutionary modelling suggests that a possible progenitor system would be a tight (P_i ~ 2d) B+B binary system that experienced conservative mass transfer. While the observed nitrogen enrichment of the primary conforms with the predictions of the evolutionary models, we find no indications for the predicted He enrichment. We suggest that HR 6819 is a binary system consisting of a stripped B-type primary and a rapidly-rotating Be star that formed from a previous mass-transfer event. In the framework of this interpretation, HR 6819 does not contain a black hole. Interferometry can distinguish between these two scenarios by providing an independent measurement of the separation between the visible components.

preprint2020arXiv

Properties of OB star-black hole systems derived from detailed binary evolution models

The recent gravitational wave measurements have demonstrated the existence of stellar mass black hole binaries. It is essential for our understanding of massive star evolution to identify the contribution of binary evolution to the formation of double black holes. A promising way to progress is investigating the progenitors of double black hole systems and comparing predictions with local massive star samples such as the population in 30 Doradus in the Large Magellanic Cloud (LMC). Methods. To this purpose, we analyse a large grid of detailed binary evolution models at LMC metallicity with initial primary masses between 10 and 40 Msun, and identify which model systems potentially evolve into a binary consisting of a black hole and a massive main sequence star. We then derive the observable properties of such systems, as well as peculiarities of the OB star component. We find that about 3% of the LMC late O and early B stars in binaries are expected to possess a black hole companion, when assuming stars with a final helium core mass above 6.6 M to form black holes. While the vast majority of them may be X-ray quiet, our models suggest that these may be identified in spectroscopic binaries, either by large amplitude radial velocity variations ( > 50 km s ) and simultaneous nitrogen surface enrichment, or through a moderate radial velocity ( > 10 km/s ) and simultaneously rapid rotation of the OB star. The predicted mass ratios are such that main sequence companions could be excluded in most cases. A comparison to the observed OB+WR binaries in the LMC, Be/X-ray binaries, and known massive BH binaries supports our conclusion. We expect spectroscopic observations to be able to test key assumptions in our models, with important implications for massive star evolution in general, and for the formation of double-black hole mergers in particular.

preprint2020arXiv

The 'hidden' companion in LB-1 unveiled by spectral disentangling

The intriguing binary LS V +22 25 (LB-1) has drawn much attention following claims of it being a single-lined spectroscopic binary with a 79-day orbit comprising a B-type star and a ~70 Msun black hole. Recent analyses have implied that the visible primary star is a stripped He-rich star. However, the nature of the secondary, which was proposed to be a black hole, a neutron star, or a main sequence star, remains unknown. Based on 26 newly acquired spectroscopic observations, we perform an orbital analysis and spectral disentangling of LB-1 to elucidate the nature of the system. Our analysis reveals that LB-1 contains two components of comparable brightness in the optical. The narrow-lined primary, which we estimate to contribute ~55% in the optical, has spectral properties that suggest that it is a stripped star: it has a small spectroscopic mass (~1 Msun) for a B-type star and it is He- and N-rich. The "hidden" secondary, which contributes about 45% of the optical flux, is a rapidly rotating (vsini ~ 300 km/s) B3 V star with a decretion disk -- a Be star. Hence, LB-1 does not contain a compact object. Instead, it is a rare Be binary system consisting of a stripped star (the former mass donor) and a Be star rotating at near its critical velocity (the former mass accretor). This system is a clear example that binary interactions play a decisive role in the production of rapid stellar rotators and Be stars.

preprint2019arXiv

The Tarantula Massive Binary Monitoring: III. Atmosphere analysis of double-lined spectroscopic systems

Accurate stellar parameters of individual objects in binary systems are essential to constrain the effects of binarity on stellar evolution. These parameters serve as a prerequisite to probing existing and future theoretical evolutionary models. We aim to derive the atmospheric parameters of the 31 SB2s in the TMBM sample. This sample, composed of detached, semi-detached and contact systems with at least one of the components classified as an O star, is an excellent test-bed to study how binarity can impact our knowledge of the evolution of massive stars. 32 epochs of FLAMES/GIRAFFE spectra are analysed using spectral disentangling to construct the individual spectra of 62 components. We apply the CMFGEN atmosphere code to determine their stellar parameters and their He, C and N surface abundances. From these properties, we show that the effects of tides on chemical mixing are limited. Components on longer-period orbits show higher nitrogen enrichment at their surface than those on shorter-period orbits, in contrast to expectations of rotational or tidal mixing, implying that other mechanisms play a role in this process. Components filling their Roche lobe are mass donors. They exhibit higher nitrogen content at their surface and rotate more slowly than their companions. By accreting new material, their companions spin faster and are rejuvenated. Their locations in the N-vsini diagram tend to show that binary products are good candidates to populate the two groups of stars (slowly rotating, nitrogen-enriched and rapidly rotating non-enriched) that cannot be reproduced through single-star population synthesis. This sample is the largest sample of binaries to be studied in such a homogeneous way. The study of these objects gives us strong observational constraints to test theoretical binary evolutionary tracks.

preprint2019arXiv

The Tarantula Massive Binary Monitoring. IV. Double-lined photometric binaries

A high fraction of massive stars are found to be binaries but only a few of them are reported as photometrically variable. By studying the populations of SB2 in the 30 Doradus region, we found a subset of them that have photometry from the OGLE project and that display variations in their light curves related to orbital motions. The goal of this study is to determine the dynamical masses and radii of the 26 binary components to investigate the mass-discrepancy problem and to provide an empirical mass-luminosity relation for the LMC. We use the PHOEBE programme to perform a systematic analysis of the OGLE V and I light curves obtained for 13 binary systems in 30 Dor. We adopt Teff, and orbital parameters derived previously to obtain the inclinations of the systems and the parameters of the individual components. Three systems display eclipses in their light curves, while the others only display ellipsoidal variations. We classify two systems as over-contact, five as semi-detached, and four as detached. The two remaining systems have uncertain configurations due to large uncertainties on their inclinations. The fact that systems display ellipsoidal variations has a significant impact on the inclination errors. From the dynamical masses, luminosities, and radii, we provide LMC-based empirical mass-luminosity and mass-radius relations, and we compare them to other relations given for the Galaxy, the LMC, and the SMC. These relations differ for different mass ranges, but do not seem to depend on the metallicity regimes. We also compare the dynamical, spectroscopic, and evolutionary masses of the stars in our sample. While the dynamical and spectroscopic masses agree with each other, the evolutionary masses are systematically higher, at least for stars in semi-detached systems. This suggests that the mass discrepancy can be partly explained by past or ongoing interactions between the stars.

preprint2019arXiv

The young massive SMC cluster NGC 330 seen by MUSE. I. Observations and stellar content

A majority of massive stars are part of binary systems, a large fraction of which will inevitably interact during their lives. Binary-interaction products (BiPs), i.e. stars affected by such interaction, are expected to be commonly present in stellar populations. BiPs are thus a crucial ingredient in the understanding of stellar evolution. We aim to identify and characterize a statistically significant sample of BiPs by studying clusters of 10-40 Myr, an age at which binary population models predict the abundance of BiPs to be highest. One example of such a cluster is NGC 330 in the Small Magellanic Cloud. Using MUSE WFM-AO observations of NGC 330, we resolve the dense cluster core for the first time and are able to extract spectra of its entire massive star population. We develop an automated spectral classification scheme based on the equivalent widths of spectral lines in the red part of the spectrum. We characterize the massive star content of the core of NGC 330 which contains more than 200 B stars, 2 O stars, 6 A-type supergiants and 11 red supergiants. We find a lower limit on the Be star fraction of $32 \pm 3\%$ in the whole sample. It increases to at least $46 \pm 10\%$ when only considering stars brighter than V=17 mag. We estimate an age of the cluster core between 35 and 40 Myr and a total cluster mass of $88^{+17}_{-18} \times 10^3 M_{\odot}$. We find that the population in the cluster core is different than the population in the outskirts: while the stellar content in the core appears to be older than the stars in the outskirts, the Be star fraction and the observed binary fraction are significantly higher. Furthermore, we detect several BiP candidates that will be subject of future studies.

preprint2016arXiv

Apsidal motion in the massive binary HD152218

Massive binary systems are important laboratories in which to probe the properties of massive stars and stellar physics in general. In this context, we analysed optical spectroscopy and photometry of the eccentric short-period early-type binary HD 152218 in the young open cluster NGC 6231. We reconstructed the spectra of the individual stars using a separating code. The individual spectra were then compared with synthetic spectra obtained with the CMFGEN model atmosphere code. We furthermore analysed the light curve of the binary and used it to constrain the orbital inclination and to derive absolute masses of 19.8 +/- 1.5 and 15.0 +/- 1.1 solar masses. Combining radial velocity measurements from over 60 years, we show that the system displays apsidal motion at a rate of (2.04^{+.23}_{-.24}) degree/year. Solving the Clairaut-Radau equation, we used stellar evolution models, obtained with the CLES code, to compute the internal structure constants and to evaluate the theoretically predicted rate of apsidal motion as a function of stellar age and primary mass. In this way, we determine an age of 5.8 +/- 0.6 Myr for HD 152218, which is towards the higher end of, but compatible with, the range of ages of the massive star population of NGC 6231 as determined from isochrone fitting.

preprint2016arXiv

Early-type massive stars in Carina Nebula within the Gaia-ESO Survey

The Gaia-ESO Survey (GES) is obtaining high quality spectra of $\sim$10$^{5}$ stars in our Galaxy, providing an homogeneous and unique overview of all the main components of the Milky Way, its formation history and the evolution of young, mature and ancient Galactic populations. Our group is in charge of the early-type massive stars that define the youngest population in the survey. In this contribution, we present the results of the quantitative spectroscopic analysis of O-type stars in the Carina Nebula within the Gaia-ESO Survey. For this aim, we have used FASTWIND and CMFGEN stellar atmosphere codes, providing stellar parameters for the current sample (GES data release iDR4).

preprint2016arXiv

Observational signatures of past mass-exchange episodes in massive binaries: The case of HD 149 404

Mass and momentum exchanges in close massive binaries play an important role in their evolution, and produce several observational signatures such as asynchronous rotation and altered chemical compositions, that remain after the stars detach again. We investigated these effects for the detached massive O-star binary HD 149404 (O7.5If + ON9.7I, P = 9.81 days), which is thought to have experienced a past episode of case A Roche-lobe overflow (RLOF). Using phase-resolved spectroscopy, we performed the disentangling of the optical spectra of the two stars. The reconstructed primary and secondary spectra were then analysed with the CMFGEN model atmosphere code to determine stellar parameters, such as the effective temperatures and surface gravities, and to constrain the chemical composition of the components. We complemented the optical study with the study of IUE spectra, which we compare to the synthetic binary spectra. The properties of the stars were compared to evolutionary models. We confirmed a strong overabundance in nitrogen ($\left[\rm N/C\right]\sim150\left[\rm N/C\right]_{\odot}$) for the secondary and a slight nitrogen overabundance ($\left[\rm N/C\right]\sim5\left[\rm N/C\right]_{\odot}$) for the primary star. Comparing the two stars, we found evidence for asynchronous rotation, with a rotational period ratio of $0.50 \pm 0.11$. The hypothesis of a past case A RLOF interaction in HD 149404 is most plausible to explain its chemical abundances and rotational asynchronicity. Some of the observed properties, such as the abundance pattern, are clearly a challenge for current case A binary evolution models, however.

preprint2016arXiv

Tracing back the evolution of the candidate LBV HD168625

Aims. We obtained Herschel/PACS imaging and spectroscopic observations of the nebula surrounding the candidate LBV HD168625. By combining these data with optical spectra of the star, we want to constrain the abundances in the nebula and in the star and compare them to trace back its evolution. Methods. We use CMFGEN to determine the fundamental parameters and the CNO abundances of the central star whilst the abundances of the nebula are derived from the emission lines present in the Herschel/PACS spectrum. Result. The FIR images show a nebula composed of an elliptical ring/torus of ejecta with a ESE-WNW axis and of a second perpendicular bipolar structure composed of empty caps/rings. We detect equatorial shells composed of dust and ionized material with different sizes when observed at different wavelengths, and bipolar caps more of less separated from the central star in H$α$ and mid-IR images. This complex global structure seems to show 2 different inclinations: 40° for the equatorial torus and 60° for the bipolar ejections. From the Herschel/PACS spectrum, we determine nebular abundances of $\mathrm{N/H} = 4.1 \pm 0.8 \times~10^{-4}$ and $\mathrm{C/H} \sim 1.6_{-0.35}^{+1.16} \times~10^{-4}$, as well as a mass of ionized gas of $0.17\pm 0.04$M$_{\odot}$ and a neutral hydrogen mass of about $1.0\pm 0.3$M$_{\odot}$ which dominates. Analysis of the central star reveals Teff = 14000K, log g = 1.74 and log (L/L$_{\odot}$) = 5.58. We derive stellar abundances of about $\mathrm{N/H} = 5.0 \pm 1.5 \times~10^{-4}$, $\mathrm{C/H} = 1.4 \pm 0.5 \times~10^{-4}$ and $\mathrm{O/H} = 3.5 \pm 1.0\times~10^{-4}$. All these measurements taken together are compatible with the evolutionary tracks of a star with an initial mass between $28$ and $33$M$_{\odot}$ and with a crit. rotational rate between 0.3 and 0.4 that has lost its material during or just after the BSG phase.

preprint2015arXiv

A spectroscopic investigation of the O-type star population in four Cygnus OB associations. II. Determination of the fundamental parameters

Aims. Having established the binary status of nineteen O-type stars located in four Cygnus OB associations, we now determine their fundamental parameters to constrain their properties and their evolutionary status. We also investigate their surface nitrogen abundances, which we compare with other results from the literature obtained for galactic O-type stars. Methods. Using optical spectra collected for each object in our sample and some UV data from the archives, we apply the CMFGEN atmosphere code to determine their main properties. For the binary systems, we have disentangled the components to obtain their individual spectra and investigate them as if they were single stars. Results. We find that the distances of several presumably single O-type stars seem poorly constrained because their luminosities are not in agreement with the "standard" luminosities of stars with similar spectral types. The ages of these O-type stars are all less than 7 Myrs. Therefore, the ages of these stars agree with those, quoted in the literature, of the four associations, except for CygOB8 for which the stars seem older than the association itself. However, we point out that the distance of certain stars is debatable relative to values found in the literature. The N content of these stars put in perspective with N contents of several other galactic O-type stars seems to draw the same five groups as found in the "Hunter" diagram for the O and B-type stars in the LMC even though their locations are obviously different. We determine mass-loss rates for several objects from the Halpha line and UV spectra. Finally, we confirm the "mass discrepancy" especially for O stars with masses smaller than 30 Msun. .

preprint2015arXiv

The mass of the very massive binary WR21a

We present multi-epoch spectroscopic observations of the massive binary system WR21a, which include the January 2011 periastron passage. Our spectra reveal multiple SB2 lines and facilitate an accurate determination of the orbit and the spectral types of the components. We obtain minimum masses of $64.4\pm4.8 \ M_{\odot}$ and $36.3\pm1.7 \ M_{\odot}$ for the two components of WR21a. Using disentangled spectra of the individual components, we derive spectral types of O3/WN5ha and O3Vz~((f*)) for the primary and secondary, respectively. Using the spectral type of the secondary as an indication for its mass, we estimate an orbital inclination of $i=58.8\pm2.5^{\mathrm{o}}$ and absolute masses of $103.6\pm10.2 \ M_{\odot}$ and $58.3\pm3.7 \ M_{\odot}$, in agreement with the luminosity of the system. The spectral types of the WR21a components indicate that the stars are very young (1$-$2 Myr), similar to the age of the nearby Westerlund 2 cluster. We use evolutionary tracks to determine the mass-luminosity relation for the total system mass. We find that for a distance of 8 kpc and an age of 1.5 Myr, the derived absolute masses are in good agreement with those from evolutionary predictions.

preprint2014arXiv

Epoch-dependent absorption line profile variability in lambda Cep

We present the analysis of a multi-epoch spectroscopic monitoring campaign of the O6Ief star lambda Cep. Previous observations reported the existence of two modes of non-radial pulsations in this star. Our data reveal a much more complex situation. The frequency content of the power spectrum considerably changes from one epoch to the other. We find no stable frequency that can unambiguously be attributed to pulsations. The epoch-dependence of the frequencies and variability patterns are similar to what is seen in the wind emission lines of this and other Oef stars, suggesting that both phenomena likely have the same, currently still unknown, origin.

preprint2014arXiv

Phase resolved X-ray spectroscopy of HDE228766: Probing the wind of an extreme Of+/WNLha star

HDE228766 is a very massive binary system hosting a secondary component, which is probably in an intermediate evolutionary stage between an Of supergiant and an WN star. The wind of this star collides with the wind of its O8 II companion, leading to relatively strong X-ray emission. Measuring the orbital variations of the line-of-sight absorption toward the X-ray emission from the wind-wind interaction zone yields information on the wind densities of both stars. X-ray spectra have been collected at three key orbital phases to probe the winds of both stars. Optical photometry has been gathered to set constraints on the orbital inclination of the system. The X-ray spectra reveal prominent variations of the intervening column density toward the X-ray emission zone, which are in line with the expectations for a wind-wind collision. We use a toy model to set constraints on the stellar wind parameters by attempting to reproduce the observed variations of the relative fluxes and wind optical depths at 1 keV. The lack of strong optical eclipses sets an upper limit of about 68 degrees on the orbital inclination. The analysis of the variations of the X-ray spectra suggests an inclination in the range 54 - 61 degrees and indicates that the secondary wind momentum ratio exceeds that of the primary by at least a factor 5. Our models further suggest that the bulk of the X-ray emission arises from the innermost region of the wind interaction zone, which is from a region whose outer radius, as measured from the secondary star, lies between 0.5 and 1.5 times the orbital separation.

preprint2013arXiv

A spectroscopic investigation of the O-type star population in four Cygnus OB associations. I. Determination of the binary fraction

Context. Establishing the multiplicity of O-type stars is the first step towards accurately determining their stellar parameters. Moreover, the distribution of the orbital parameters provides observational clues to the way that O-type stars form and to the interactions during their evolution. Aims. Our objective is to constrain the multiplicity of a sample of O-type stars belonging to poorly investigated OB associations in the Cygnus complex and for the first time to provide orbital parameters for binaries identified in our sample. Such information is relevant to addressing the issue of the binarity in the context of O-type star formation scenarios. Methods. We performed a long-term spectroscopic survey of nineteen O-type stars. We searched for radial velocity variations to unveil binaries on timescales from a few days up to a few years, on the basis of a large set of optical spectra. Results. We confirm the binarity for four objects: HD193443, HD228989, HD229234 and HD194649. We derive for the first time the orbital solutions of three systems, and we confirm the values of the fourth, showing that these four systems all have orbital periods shorter than 10 days. Besides these results, we also detect several objects that show non-periodic line profile variations in some of their spectral lines. These variations mainly occur in the spectral lines, that are generally affected by the stellar wind and are not likely to be related to binarity. Conclusions. The minimal binary fraction in our sample is estimated to be 21%, but it varies from one OB association to the next. Indeed, 33% of O stars belonging to CygOB1 are binary systems, 0% in CygOB3, 0% in CygOB8, and 33% in CygOB9. Our spectroscopic investigation also stresses the absence of long-period systems among the stars in our sample...

preprint2013arXiv

Three-dimensional orbits of the triple-O stellar system HD 150136

Context. HD 150136 is a triple hierarchical system and a non-thermal radio emitter. It is formed by an O3-3.5 V + O5.5-6 V close binary and a more distant O6.5-7 V tertiary. So far, only the inner orbital properties have been reliably constrained. Aims. To quantitatively understand the non-thermal emission process, accurate knowledge of the physical and orbital properties of the object is crucial. Here, we aim to investigate the orbital properties of the wide system and to constrain the inclinations of the inner and outer binaries, and with these the absolute masses of the system components. Methods. We used the PIONIER combiner at the Very Large Telescope Interferometer to obtain the very first interferometric measurements of HD 150136. We combined the interferometric observations with new and existing high resolution spectroscopic data to derive the orbital solution of the outer companion in the three-dimensional space. Results. The wide system is clearly resolved by PIONIER, with a projected separation on the plane of the sky of about 9 milli-arcsec. The best-fit orbital period, eccentricity, and inclination are 8.2 yr, 0.73 and 108 degr. We constrain the masses of the three stars of the system to 63 +/- 10, 40 +/- 6, and 33 +/- 12 Msun for the O3-3.5 V, O5.5-6 V and O6.5-7 V components. Conclusions. The dynamical masses agree within errors with the evolutionary masses of the components. Future interferometric and spectroscopic monitoring of HD 150136 should allow one to reduce the uncertainties to a few per cent only and to accurately constrain the distance to the system. This makes HD 150136 an ideal system to quantitatively test evolutionary models of high-mass stars as well as the physics of non-thermal processes occurring in O-type systems.

preprint2012arXiv

9 Sgr: uncovering an O-type spectroscopic binary with an 8.6 year period

The O-type object 9 Sgr is a well-known synchrotron radio emitter. This feature is usually attributed to colliding-wind binary systems, but 9 Sgr was long considered a single star. We have conducted a long-term spectroscopic monitoring of this star to investigate its multiplicity and search for evidence for wind-wind interactions. Radial velocities are determined and analysed using various period search methods. Spectral disentangling is applied to separate the spectra of the components of the binary system. We derive the first ever orbital solution of 9 Sgr. The system is found to consist of an O3.5 V((f+)) primary and an O5- 5.5 V((f)) secondary moving around each other on a highly eccentric (e = 0.7), 8.6 year orbit. The spectra reveal no variable emission lines that could be formed in the wind interaction zone in agreement with the expected properties of the interaction in such a wide system. Our results provide further support to the paradigm of synchrotron radio emission from early-type stars being a manifestation of interacting winds in a binary system.

preprint2012arXiv

Evidence for a physically bound third component in HD 150136

Context. HD150136 is one of the nearest systems harbouring an O3 star. Although this system was for a long time considered as binary, more recent investigations have suggested the possible existence of a third component. Aims. We present a detailed analysis of HD 150136 to confirm the triple nature of this system. In addition, we investigate the physical properties of the individual components of this system. Methods. We analysed high-resolution, high signal-to-noise data collected through multi-epoch runs spread over ten years. We applied a disentangling program to refine the radial velocities and to obtain the individual spectra of each star. With the radial velocities, we computed the orbital solution of the inner system, and we describe the main properties of the orbit of the outer star such as the preliminary mass ratio, the eccentricity, and the orbital-period range. With the individual spectra, we determined the stellar parameters of each star by means of the CMFGEN atmosphere code. Results. We offer clear evidence that HD 150136 is a triple system composed of an O3V((f\ast))-3.5V((f+)), an O5.5-6V((f)), and an O6.5-7V((f)) star. The three stars are between 0-3 Myr old. We derive dynamical masses of about 64, 40, and 35 Msun for the primary, the secondary and the third components by assuming an inclination of 49°. It currently corresponds to one of the most massive systems in our galaxy. The third star moves with a period in the range of 2950 to 5500 d on an outer orbit with an eccentricity of at least 0.3. This discovery makes HD 150136 the first confirmed triple system with an O3 primary star. However, because of the long orbital period, our dataset is not sufficient to constrain the orbital solution of the tertiary component with high accuracy.

preprint2012arXiv

Massive open star clusters using the VVV survey II. Discovery of six clusters with Wolf-Rayet stars

Context: The ESO Public Survey "VISTA Variables in the Vía Láctea" (VVV) provides deep multi-epoch infrared observations for an unprecedented 562 sq. degrees of the Galactic bulge, and adjacent regions of the disk. In this survey nearly 150 new open clusters and cluster candidates have been discovered. Aims: This is the second in a series of papers about young, massive open clusters observed using the VVV survey. We present the first study of six recently discovered clusters. These clusters contain at least one newly discovered Wolf-Rayet (WR) star. Methods: Following the methodology presented in the first paper of the series, wide-field, deep JHKs VVV observations, combined with new infrared spectroscopy, are employed to constrain fundamental parameters for a subset of clusters. Results: We affirm that the six studied stellar groups are real young (2-7 Myr) and massive (between 0.8 and 2.2 10^3 Msol) clusters. They are highly obscured (Av ~ 5-24 mag) and compact (1-2 pc). In addition to WR stars, two of the six clusters also contain at least one red supergiant star. We claim the discovery of 8 new WR stars, and 3 stars showing WR-like emission lines which could be classified WR or OIf. Preliminary analysis provides initial masses of ~30-50 Msol for the WR stars. Finally,we discuss the spiral structure of the Galaxy using as tracers the six new clusters together with the previously studied VVV clusters.

preprint2012arXiv

The 2.35 year itch of Cyg OB2 #9 - I. Optical and X-ray monitoring

Context: Nonthermal radio emission in massive stars is expected to arise in wind-wind collisions occurring inside a binary system. One such case, the O-type star Cyg OB2 #9, was proven to be a binary only four years ago, but the orbital parameters remained uncertain. The periastron passage of 2011 was the first one to be observable under good conditions since the discovery of binarity. Aims: In this context, we have organized a large monitoring campaign to refine the orbital solution and to study the wind-wind collision. Methods: This paper presents the analysis of optical spectroscopic data, as well as of a dedicated X-ray monitoring performed with Swift and XMM. Results: In light of our refined orbital solution, Cyg OB2 #9 appears as a massive O+O binary with a long period and high eccentricity; its components (O5-5.5I for the primary and O3-4III for the secondary) have similar masses and similar luminosities. The new data also provide the first evidence that a wind-wind collision is present in the system. In the optical domain, the broad Ha line varies, displaying enhanced absorption and emission components at periastron. X-ray observations yield the unambiguous signature of an adiabatic collision because, as the stars approach periastron, the X-ray luminosity closely follows the 1/D variation expected in that case. The X-ray spectrum appears, however, slightly softer at periastron, which is probably related to winds colliding at slightly lower speeds at that time. Conclusions: It is the first time that such a variation has been detected in O+O systems, and the first case where the wind-wind collision is found to remain adiabatic even at periastron passage.

preprint2011arXiv

A quantitative study of O stars in NGC2244 and the Mon OB2 association

Our goal is to determine the stellar and wind properties of seven O stars in the cluster NGC2244 and three O stars in the OB association MonOB2. These properties give us insight into the mass loss rates of O stars, allow us to check the validity of rotational mixing in massive stars, and to better understand the effects of the ionizing flux and wind mechanical energy release on the surrounding interstellar medium and its influence on triggered star formation. We collect optical and UV spectra of the target stars which are analyzed by means of atmosphere models computed with the code CMFGEN. The spectra of binary stars are disentangled and the components are studied separately. All stars have an evolutionary age less than 5 million years, with the most massive stars being among the youngest. Nitrogen surface abundances show no clear relation with projected rotational velocities. Binaries and single stars show the same range of enrichment. This is attributed to the youth and/or wide separation of the binary systems in which the components have not (yet) experienced strong interaction. A clear trend of larger enrichment in higher luminosity objects is observed, consistent with what evolutionary models with rotation predict for a population of O stars at a given age. We confirm the weakness of winds in late O dwarfs. In general, mass loss rates derived from UV lines are lower than mass loss rates obtained from Ha. The UV mass loss rates are even lower than the single line driving limit in the latest type dwarfs. These issues are discussed in the context of the structure of massive stars winds. The evolutionary and spectroscopic masses are in agreement above 25 Msun but the uncertainties are large. Below this threshold, the few late-type O stars studied here indicate that the mass discrepancy still seems to hold.

preprint2011arXiv

CoRoT observations of O stars: diverse origins of variability

Six O-type stars were observed continuously by the CoRoT satellite during a 34.3-day run. The unprecedented quality of the data allows us to detect even low-amplitude stellar pulsations in some of these stars (HD 46202 and the binaries HD 46149 and Plaskett's star). These cover both opacity-driven modes and solar-like stochastic oscillations, both of importance to the asteroseismological modelling of O stars. Additional effects can be seen in the CoRoT light curves, such as binarity and rotational modulation. Some of the hottest O-type stars (HD 46223, HD 46150 and HD 46966) are dominated by the presence of red-noise: we speculate that this is related to a sub-surface convection zone.

preprint2011arXiv

The two components of the evolved massive binary LZ Cep. Testing the effects of binarity on stellar evolution

Aims. We present an in-depth study of the two components of the binary system LZCep in order to constrain the effects of binarity on the evolution of massive stars. Methods. We use a set of high-resolution, high signal-to-noise ratio optical spectra obtained over the orbital period of the system to perform a spectroscopic disentangling and derive an orbital solution. We subsequently determine the stellar properties of each component through an analysis with the CMFGEN atmosphere code. Finally, with the derived stellar parameters, we model the Hipparcos photometric light curve using the program NIGHTFALL to obtain the inclination and the real stellar masses. Results. LZCep is a O9III+ON9.7V binary. It is as a semi-detached system in which either the primary or the secondary star almost fills up its Roche lobe. The dynamical masses are about 16.0 $M_{\odot}$ (primary) and 6.5 M_{\odot}$ (secondary). The latter is lower than the typical mass of late-type O stars. The secondary component is chemically more evolved than the primary (which barely shows any sign of CNO processing), with strong helium and nitrogen enhancements as well as carbon and oxygen depletions. These properties (surface abundances and mass) are typical of Wolf-Rayet stars, although the spectral type is ON9.7V. The luminosity of the secondary is consistent with that of core He-burning objects. The preferred, tentative evolutionary scenario to explain the observed properties involves mass transfer from the secondary - which was initially more massive - towards the primary. The secondary is now almost a core He-burning object, probably with only a thin envelope of H-rich and CNO processed material. A very inefficient mass transfer is necessary to explain the chemical appearance of the primary. Alternative scenarios are discussed but they suffer from more uncertainties.

preprint2011arXiv

Variability in the CoRoT photometry of three hot O-type stars. HD 46223, HD 46150 and HD 46966

The detection of pulsational frequencies in stellar photometry is required as input for asteroseismological modelling. The second short run (SRa02) of the CoRoT mission has provided photometric data of unprecedented quality and time-coverage for a number of O-type stars. We analyse the CoRoT data corresponding to three hot O-type stars, describing the properties of their light curves and we search for pulsational frequencies, which we then compare to theoretical model predictions. We determine the amplitude spectrum of the data, using the Lomb-Scargle and a multifrequency HMM-like technique. Frequencies are extracted by prewhitening, and their significance is evaluated under the assumption that the light curve is dominated by red noise. We search for harmonics, linear combinations and regular spacings among these frequencies. We use simulations with the same time sampling as the data as a powerful tool to judge the significance of our results. From the theoretical point of view, we use the MAD non-adiabatic pulsation code to determine the expected frequencies of excited modes. A substantial number of frequencies is listed, but none can be convincingly identified as being connected to pulsations. The amplitude spectrum is dominated by red noise. Theoretical modelling shows that all three O-type stars can have excited modes but the relation between the theoretical frequencies and the observed spectrum is not obvious. The dominant red noise component in the hot O-type stars studied here clearly points to a different origin than the pulsations seen in cooler O stars. The physical cause of this red noise is unclear, but we speculate on the possibility of sub-surface convection, granulation, or stellar wind inhomogeneities being responsible.

preprint2010arXiv

Detection of frequency spacings in the young O-type binary HD 46149 from CoRoT photometry

Using the CoRoT space based photometry of the O-type binary HD46149, stellar atmospheric effects related to rotation can be separated from pulsations, because they leave distinct signatures in the light curve. This offers the possibility of characterising and exploiting any pulsations seismologically. Combining high-quality space based photometry, multi-wavelength photometry, spectroscopy and constraints imposed by binarity and cluster membership, the detected pulsations in HD46149 are analyzed and compared with those for a grid of stellar evolutionary models in a proof-of-concept approach. We present evidence of solar-like oscillations in a massive O-type star, and show that the observed frequency range and spacings are compatible with theoretical predictions. Thus, we unlock and confirm the strong potential of this seismically unexplored region in the HR diagram.

preprint2010arXiv

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.

preprint2010arXiv

Plaskett's Star: Analysis of the CoRoT photometric data

The SRa02 of the CoRoT space mission for Asteroseismology was partly devoted to stars belonging to the Mon OB2 association. An intense monitoring was performed on Plaskett's Star (HD47129) and the unprecedented quality of the light curve allows us to shed new light on this very massive, non-eclipsing binary system. We particularly aimed at detecting periodic variability which might be associated with pulsations or interactions between both components. We also searched for variations related to the orbital cycle which could help to constrain the inclination and the morphology of the binary system. A Fourier-based prewhitening and a multiperiodic fitting procedure were applied to analyse the time series and extract the frequencies of variations. We describe the noise properties to tentatively define an appropriate significance criterion, to only point out the peaks at a certain significance level. We also detect the variations related to the orbital motion and study them by using the NIGHTFALL program. The periodogram exhibits a majority of peaks at low frequencies. Among these peaks, we highlight a list of about 43 values, including notably two different sets of harmonic frequencies whose fundamental peaks are located at about 0.07 and 0.82d-1. The former represents the orbital frequency of the binary system whilst the latter could probably be associated with non-radial pulsations. The study of the 0.07d-1 variations reveals the presence of a hot spot most probably situated on the primary star and facing the secondary. The investigation of this unique dataset constitutes a further step in the understanding of Plaskett's Star. These results provide a first basis for future seismic modelling. The existence of a hot region between both components renders the determination of the inclination ambiguous.

preprint2009arXiv

A new investigation of the binary HD 48099

With an orbital period of about 3.078 days, the double-lined spectroscopic binary HD 48099 is, still now, the only short-period O+O system known in the Mon OB2 association. Even though an orbital solution has already been derived for this system, few information are available about the individual stars. We present, in this paper, the results of a long-term spectroscopic campaign. We derive a new orbital solution and apply a disentangling method to recover the mean spectrum of each star. To improve our knowledge concerning both components, we determine their spectral classifications and their projected rotational velocities. We also constrain the main stellar parameters of both stars by using the CMFGEN atmosphere code and provide the wind properties for the primary star through the study of IUE spectra. This investigation reveals that HD 48099 is an O5.5 V((f))+O9 V binary with M_1 sin^3 i = 0.70 M_{\sun} and M_2 sin^3 i = 0.39 M_{\sun}, implying a rather low orbital inclination. This result, combined with both a large effective temperature and log g, suggests that the primary star (v sini ~ 91 km s^-1) is actually a fast rotator with a strongly clumped wind and a nitrogen abundance of about 8 times the solar value.