-
Empirical mass-loss rates and clumping properties of O-type stars in the LMC
Authors:
C. Hawcroft,
L. Mahy,
H. Sana,
J. O. Sundqvist,
M. Abdul-Masih,
S. A. Brands,
L. Decin,
A. deKoter,
J. Puls
Abstract:
We constrain wind parameters of a sample of 18 O-type stars in the LMC, through analysis with stellar atmosphere and wind models including the effects of optically thick clumping. This allows us to determine the most accurate spectroscopic mass-loss and wind structure properties of massive stars at sub-solar metallicity to date and gain insight into the impact of metallicity on massive stellar win…
▽ More
We constrain wind parameters of a sample of 18 O-type stars in the LMC, through analysis with stellar atmosphere and wind models including the effects of optically thick clumping. This allows us to determine the most accurate spectroscopic mass-loss and wind structure properties of massive stars at sub-solar metallicity to date and gain insight into the impact of metallicity on massive stellar winds. Combining high signal to noise (S/N) ratio spectroscopy in the UV and optical gives us access to diagnostics of multiple different physical processes in the stellar wind. We produce synthetic spectra using the stellar atmosphere modelling code FASTWIND, and reproduce the observed spectra using a genetic algorithm based fitting technique. We empirically constrain 15 physical parameters associated with the stellar and wind properties, including temperature, surface gravity, surface abundances, rotation, macroturbulence and wind parameters. We find, on average, mass-loss rates a factor of 4-5 lower than those predicted by Vink et al. 2001, in good agreement with predictions from Bjorklund et al. 2021, and the best agreement with those from Krticka et al. 2018. In the 'weak-wind' regime we find mass-loss rates orders of magnitude below any theoretical predictions. We find a positive correlation of clumping factors (fcl) with effective temperature with an average fcl = 14 +- 8 for the full sample. Above 38 kK an average 46 +- 24% of the wind velocity span is covered by clumps and the interclump density is 10-30% of the mean wind. Below an effective temperature of roughly 38 kK there must be additional light leakage for supergiants. For dwarf stars at low temperatures there is a statistical preference for very low clump velocity spans, however it is unclear if this can be physically motivated as there are no clearly observable wind signatures in UV diagnostics.
△ Less
Submitted 9 July, 2024;
originally announced July 2024.
-
Detection of extragalactic magnetic massive stars
Authors:
S. Hubrig,
M. Schöller,
S. P. Jarvinen,
A. Cikota,
M. Abdul-Masih,
A. Escorza,
R. Jayaraman
Abstract:
Studies of the magnetic characteristics of massive stars have recently received significant attention because they are progenitors of highly magnetised compact objects. Stars initially more massive than about 8M_sun leave behind neutron stars and black holes by the end of their evolution. The merging of binary compact remnant systems produces astrophysical transients detectable by gravitational wa…
▽ More
Studies of the magnetic characteristics of massive stars have recently received significant attention because they are progenitors of highly magnetised compact objects. Stars initially more massive than about 8M_sun leave behind neutron stars and black holes by the end of their evolution. The merging of binary compact remnant systems produces astrophysical transients detectable by gravitational wave observatories. Studies of magnetic fields in massive stars with low metallicities are of particular interest because they provide important information on the role of magnetic fields in the star formation of the early Universe. While several detections of massive Galactic magnetic stars have been reported in the last few decades, the impact of a low-metallicity environment on the occurrence and strength of stellar magnetic fields has not yet been explored. Because of the similarity between Of?p stars in the Magellanic Clouds (MCs) and Galactic magnetic Of?p stars, which possess globally organised magnetic fields, we searched for magnetic fields in Of?p stars in the MCs. Additionally, we observed the massive contact binary Cl* NGC 346 SSN7 in the Small Magellanic Cloud to test the theoretical scenario that the origin of magnetic fields involves a merger event or a common envelope evolution. We obtained and analysed measurements of the magnetic field in four massive Of?p stars in the MCs and the binary Cl* NGC 346 SSN7 using the ESO/VLT FORS2 spectrograph in spectropolarimetric mode. We detected kilogauss-scale magnetic fields in two Of?p-type stars and in the contact binary Cl* NGC 346 SSN7. These results suggest that the impact of low metallicity on the occurrence and strength of magnetic fields in massive stars is low. However, because the explored stellar sample is very small, additional observations of massive stars in the MCs are necessary.
△ Less
Submitted 30 May, 2024;
originally announced May 2024.
-
The colliding-wind binary HD 168112
Authors:
R. Blomme,
G. Rauw,
D. Volpi,
Y. Nazé,
M. Abdul-Masih
Abstract:
Radio surveys of early-type stars have revealed a number of non-thermal emitters. Most of these have been shown to be binaries, where the collision between the two stellar winds is responsible for the non-thermal emission. HD 168112 is a non-thermal radio emitter, whose binary nature has only recently been confirmed spectroscopically. We obtained independent spectroscopic observations to determine…
▽ More
Radio surveys of early-type stars have revealed a number of non-thermal emitters. Most of these have been shown to be binaries, where the collision between the two stellar winds is responsible for the non-thermal emission. HD 168112 is a non-thermal radio emitter, whose binary nature has only recently been confirmed spectroscopically. We obtained independent spectroscopic observations to determine its orbit, in addition to radio observations to see if the thermal or non-thermal nature of the emission changes during the periastron passage. We monitored HD 168112 spectroscopically for a 13 year time span. From these data, we determined the orbital parameters, which we compared to the previous results in the literature. From the spectral index of the radio observations, we found how the nature of the emission changes as the system goes through periastron. Combining our results with other literature data allowed us to further constrain the orbital and stellar parameters. We find HD 168112 to have an orbital period of P = 512.17+0.41-0.11 d, an eccentricity of e = 0.7533+0.0053-0.0124, and a mass ratio close to one. From our spectroscopic modelling, we derived the stellar parameters, but we had difficulty arriving at a spectroscopic mass ratio of one. The radio observations around periastron show only thermal emission, suggesting that most of the synchrotron photons are absorbed in the two stellar winds at that phase. Combining our data with the optical interferometry detection, we could constrain the inclination angle to i ~ 63 deg, and the mass of each component to ~ 26 Msun. We have provided an independent spectroscopic confirmation of the binary nature of HD 168112. Although detected as a non-thermal radio emitter, near periastron the radio emission of this highly eccentric system is thermal and is mainly formed in the colliding-wind region. [abridged]
△ Less
Submitted 3 June, 2024; v1 submitted 6 May, 2024;
originally announced May 2024.
-
X-Shooting ULLYSES: Massive Stars at low metallicity II. DR1: Advanced optical data products for the Magellanic Clouds
Authors:
H. Sana,
F. Tramper,
M. Abdul-Masih,
R. Blomme,
K. Dsilva,
G. Maravelias,
L. Martins,
A. Mehner,
A. Wofford,
G. Banyard,
C. L. Barbosa,
J. Bestenlehner,
C. Hawcroft,
D. John Hillier,
H. Todt,
C. J. K. Larkin,
L. Mahy,
F. Najarro,
V. Ramachandran,
M. C. Ramirez-Tannus,
M. M. Rubio-Diez,
A. A. C. Sander,
T. Shenar,
J. S. Vink,
F. Backs
, et al. (12 additional authors not shown)
Abstract:
Using the medium resolution spectrograph X-shooter, spectra of 235 OB and Wolf-Rayet (WR) stars in sub-solar metallicity environments have been secured. [...]This second paper focuses on the optical observations of 232 Magellanic Clouds targets. It describes the uniform reduction of the UVB (300 - 560 nm) and VIS (550 - 1020 nm) XShootU data as well as the preparation of advanced data products [..…
▽ More
Using the medium resolution spectrograph X-shooter, spectra of 235 OB and Wolf-Rayet (WR) stars in sub-solar metallicity environments have been secured. [...]This second paper focuses on the optical observations of 232 Magellanic Clouds targets. It describes the uniform reduction of the UVB (300 - 560 nm) and VIS (550 - 1020 nm) XShootU data as well as the preparation of advanced data products [...] . The data reduction of the raw data is based on the ESO CPL X-shooter pipeline. We paid particular attention to the determination of the response curves [...] We implemented slit-loss correction, absolute flux calibration, (semi-)automatic rectification to the continuum, and a correction for telluric lines. The spectra of individual epochs were corrected for the barycentric motion, re-sampled and co-added, and the spectra from the two arms were merged into a single flux calibrated spectrum covering the entire optical range with maximum signal-to-noise ratio. [...] We provide three types of data products: (i) two-dimensional spectra for each UVB and VIS exposure; (ii) one-dimensional UVB and VIS spectra before and after response-correction, as well as after applying various processing, including absolute flux calibration, telluric removal, normalisation and barycentric correction; and (iii) co-added flux-calibrated and rectified spectra over the full optical range, for which all available XShootU exposures were combined. For many of the targets, the final signal-to-noise ratio per resolution element is above 200 in both the UVB and the VIS co-added spectra. The reduced data and advanced scientific data products will be made available to the community upon publication of this paper. [...]
△ Less
Submitted 26 February, 2024;
originally announced February 2024.
-
MELCHIORS: The Mercator Library of High Resolution Stellar Spectroscopy
Authors:
P. Royer,
T. Merle,
K. Dsilva,
S. Sekaran,
H. Van Winckel,
Y. Frémat,
M. Van der Swaelmen,
S. Gebruers,
A. Tkachenko,
M. Laverick,
M. Dirickx,
G. Raskin,
H. Hensberge,
M. Abdul-Masih,
B. Acke,
M. L. Alonso,
S. Bandhu Mahato,
P. G. Beck,
N. Behara,
S. Bloemen,
B. Buysschaert,
N. Cox,
J. Debosscher,
P. De Cat,
P. Degroote
, et al. (49 additional authors not shown)
Abstract:
Over the past decades, libraries of stellar spectra have been used in a large variety of science cases, including as sources of reference spectra for a given object or a given spectral type. Despite the existence of large libraries and the increasing number of projects of large-scale spectral surveys, there is to date only one very high-resolution spectral library offering spectra from a few hundr…
▽ More
Over the past decades, libraries of stellar spectra have been used in a large variety of science cases, including as sources of reference spectra for a given object or a given spectral type. Despite the existence of large libraries and the increasing number of projects of large-scale spectral surveys, there is to date only one very high-resolution spectral library offering spectra from a few hundred objects from the southern hemisphere (UVES-POP) . We aim to extend the sample, offering a finer coverage of effective temperatures and surface gravity with a uniform collection of spectra obtained in the northern hemisphere.
Between 2010 and 2020, we acquired several thousand echelle spectra of bright stars with the Mercator-HERMES spectrograph located in the Roque de Los Muchachos Observatory in La Palma, whose pipeline offers high-quality data reduction products. We have also developed methods to correct for the instrumental response in order to approach the true shape of the spectral continuum. Additionally, we have devised a normalisation process to provide a homogeneous normalisation of the full spectral range for most of the objects.
We present a new spectral library consisting of 3256 spectra covering 2043 stars. It combines high signal-to-noise and high spectral resolution over the entire range of effective temperatures and luminosity classes. The spectra are presented in four versions: raw, corrected from the instrumental response, with and without correction from the atmospheric molecular absorption, and normalised (including the telluric correction).
△ Less
Submitted 5 November, 2023;
originally announced November 2023.
-
X-Shooting ULLYSES: massive stars at low metallicity. I. Project Description
Authors:
Jorick S. Vink,
A. Mehner,
P. A. Crowther,
A. Fullerton,
M. Garcia,
F. Martins,
N. Morrell,
L. M. Oskinova,
N. St-Louis,
A. ud-Doula,
A. A. C. Sander,
H. Sana,
J. -C. Bouret,
B. Kubatova,
P. Marchant,
L. P. Martins,
A. Wofford,
J. Th. van Loon,
O. Grace Telford,
Y. Gotberg,
D. M. Bowman,
C. Erba,
V. M. Kalari,
M. Abdul-Masih,
T. Alkousa
, et al. (56 additional authors not shown)
Abstract:
Observations of individual massive stars, super-luminous supernovae, gamma-ray bursts, and gravitational-wave events involving spectacular black-hole mergers, indicate that the low-metallicity Universe is fundamentally different from our own Galaxy. Many transient phenomena will remain enigmatic until we achieve a firm understanding of the physics and evolution of massive stars at low metallicity…
▽ More
Observations of individual massive stars, super-luminous supernovae, gamma-ray bursts, and gravitational-wave events involving spectacular black-hole mergers, indicate that the low-metallicity Universe is fundamentally different from our own Galaxy. Many transient phenomena will remain enigmatic until we achieve a firm understanding of the physics and evolution of massive stars at low metallicity (Z). The Hubble Space Telescope has devoted 500 orbits to observe 250 massive stars at low Z in the ultraviolet (UV) with the COS and STIS spectrographs under the ULLYSES program. The complementary ``X-Shooting ULLYSES'' (XShootU) project provides enhanced legacy value with high-quality optical and near-infrared spectra obtained with the wide-wavelength coverage X-shooter spectrograph at ESO's Very Large Telescope.
We present an overview of the XShootU project, showing that combining ULLYSES UV and XShootU optical spectra is critical for the uniform determination of stellar parameters such as effective temperature, surface gravity, luminosity, and abundances, as well as wind properties such as mass-loss rates in function of Z. As uncertainties in stellar and wind parameters percolate into many adjacent areas of Astrophysics, the data and modelling of the XShootU project is expected to be a game-changer for our physical understanding of massive stars at low Z.
To be able to confidently interpret James Webb Space Telescope (JWST) spectra of the first stellar generations, the individual spectra of low Z stars need to be understood, which is exactly where XShootU can deliver.
△ Less
Submitted 1 June, 2023; v1 submitted 10 May, 2023;
originally announced May 2023.
-
The effects of rotation on massive star spectroscopic observables
Authors:
Michael Abdul-Masih
Abstract:
Rotation is ubiquitous among massive stars. With rotation comes a deformation to the surface geometry, which in turn leads to alterations in the distribution of parameters across the surface including surface gravity, temperature and ionization balance of surface elements. Often, these 3D effects are neglected when analyzing spectra of rapidly rotating massive stars. We aim to determine whether ne…
▽ More
Rotation is ubiquitous among massive stars. With rotation comes a deformation to the surface geometry, which in turn leads to alterations in the distribution of parameters across the surface including surface gravity, temperature and ionization balance of surface elements. Often, these 3D effects are neglected when analyzing spectra of rapidly rotating massive stars. We aim to determine whether neglecting the 3D deformations resulting from rapid rotation has an impact on the final spectroscopic observables, and if so to what degree. Using the SPAMMS code, we generate a grid of synthetic spectra that account for the 3D geometry of rapidly rotating stars and compare them to synthetic spectra generated assuming spherical geometry. Using equivalent width and full width half maximum measurements as proxies, we determine how the measured temperature, helium abundance and projected rotation rates of individual lines in different ionization states vary with rotation rate and inclination. We find that the 3D geometry can have a significant impact on the measured parameters. We show that the temperature is highly dependent on both the rotation rate and the inclination, and that the same system viewed at different inclinations can have measured temperatures that differ by as much as 10\%. We also find that the helium abundance can be underestimated by as much as 60\%, and that lines in different ionization states can have measurable differences in rotation rates. We demonstrate that these differences in rotation rates can be seen in observed data and show that this could allow for an inclination independent measurement of the rotational velocity. Our results indicate that neglecting the 3D effects of rotation can cause significant biases in the measured spectroscopic parameters, and that in many cases, the measured values are more than 3$σ$ away from the true values.
△ Less
Submitted 20 December, 2022;
originally announced December 2022.
-
Distance and age of the massive stellar cluster Westerlund 1. II. The eclipsing binary W36
Authors:
Danilo F. Rocha,
Leonardo A. Almeida,
Augusto Damineli,
Felipe Navarete,
Michael Abdul-Masih,
Gregory N. Mace
Abstract:
Westerlund 1 (Wd 1) is one of the most relevant star clusters in the Milky Way to study massive star formation, although it is still poorly known. Here, we used photometric and spectroscopic data to model the eclipsing binary W36, showing that its spectral type is O6.5 III + O9.5 IV, hotter and more luminous than thought before. Its distance $d_{\rm W36}$ $=$ 4.03$\pm$0.25 kpc agrees, within the e…
▽ More
Westerlund 1 (Wd 1) is one of the most relevant star clusters in the Milky Way to study massive star formation, although it is still poorly known. Here, we used photometric and spectroscopic data to model the eclipsing binary W36, showing that its spectral type is O6.5 III + O9.5 IV, hotter and more luminous than thought before. Its distance $d_{\rm W36}$ $=$ 4.03$\pm$0.25 kpc agrees, within the errors, with three recent Gaia-EDR3-based distances reported in Paper I, Beasor & Davies, and by Negueruela's group. However, they follow different approaches to fix the zero-points for red sources such as those in Wd 1 and to select the best approach, we used an accurate modelling of W36. The weighted mean distance of our parallax (Paper I) and binary distances results in $d_{\rm wd1}$ = 4.05 $\pm$0.20 kpc, with an unprecedented accuracy of 5%. We adopted isochrones based on the Geneva code with supersolar abundances to infer the age of W36B as 6.4 $\pm$ 0.7 Myr. This object seems to be part of the prolific star formation burst represented by OB giants and supergiants that occurred at 7.1 $\pm$ 0.5 Myr ago, which coincides with the recently published PMS isochrone with age 7.2 Myr. Other BA-type luminous evolved stars and Yellow Hypergiants spread in the age range of 8--11 Myr. The four Red Supergiants discussed in paper I represent the oldest population of the cluster with an age of 10.7 $\pm$ 1 Myr. The multiple episodes of star formation in Wd 1 are reminiscent of that reported for the R136/30 Dor LMC cluster.
△ Less
Submitted 10 October, 2022;
originally announced October 2022.
-
The Nature of Unseen Companions in Massive Single-Line Spectroscopic Binaries
Authors:
Hugues Sana,
Michael Abdul-Masih,
Gareth Banyard,
Julia Bodensteiner,
Dominic M. Bowman,
Karan Dsilva,
C. Eldridge,
Matthias Fabry,
Abigail J. Frost,
Calum Hawcroft,
Soetkin Janssens,
Laurent Mahy,
Pablo Marchant,
Norbert Langer,
Timothy Van Reeth,
Koushik Sen,
Tomer Shenar
Abstract:
Massive stars are predominantly found in binaries and higher order multiples. While the period and eccentricity distributions of OB stars are now well established across different metallicity regimes, the determination of mass-ratios has been mostly limited to double-lined spectroscopic binaries. As a consequence, the mass-ratio distribution remains subject to significant uncertainties. Open quest…
▽ More
Massive stars are predominantly found in binaries and higher order multiples. While the period and eccentricity distributions of OB stars are now well established across different metallicity regimes, the determination of mass-ratios has been mostly limited to double-lined spectroscopic binaries. As a consequence, the mass-ratio distribution remains subject to significant uncertainties. Open questions include the shape and extent of the companion mass-function towards its low-mass end and the nature of undetected companions in single-lined spectroscopic binaries. In this contribution, we present the results of a large and systematic analysis of a sample of over 80 single-lined O-type spectroscopic binaries (SB1s) in the Milky Way and in the Large Magellanic Cloud (LMC). We report on the developed methodology, the constraints obtained on the nature of SB1 companions, the distribution of O star mass-ratios at LMC metallicity and the occurrence of quiescent OB+black hole binaries.
△ Less
Submitted 13 September, 2022;
originally announced September 2022.
-
Constraining the overcontact phase in massive binary evolution -- II. Period stability of known O+O overcontact systems
Authors:
Michael Abdul-Masih,
Ana Escorza,
Athira Menon,
Laurent Mahy,
Pablo Marchant
Abstract:
Given that mergers are often invoked to explain many exotic phenomena in massive star evolution, understanding the evolutionary phase directly preceding a merger, the overcontact phase, is of crucial importance. Despite its importance, large uncertainties exist in our understanding of the evolution of massive overcontact binaries. We aim to provide robust observational constraints on the future dy…
▽ More
Given that mergers are often invoked to explain many exotic phenomena in massive star evolution, understanding the evolutionary phase directly preceding a merger, the overcontact phase, is of crucial importance. Despite its importance, large uncertainties exist in our understanding of the evolution of massive overcontact binaries. We aim to provide robust observational constraints on the future dynamical evolution of massive overcontact systems by measuring the rate at which the periods change for a sample of six such objects. Furthermore, we aim to investigate whether the periods of unequal mass systems show higher rates of change than their equal mass counterparts as theoretical models predict. Using archival photometric data from various ground- and space-based missions covering up to ~40 years, we measure the periods of each system over several smaller time spans. We then fit a linear regression through the measured periods to determine the rate at which the period is changing over the entire data set. We find that all of the stars in our sample have very small period changes and that there does not seem to be a correlation with the mass ratio. This implies that the orbital periods for these systems are stable on the nuclear timescale, and that the unequal mass systems may not equalize as expected. When comparing our results with population synthesis distributions, we find large discrepancies between the expected mass ratios and period stabilities. We find that these discrepancies can be mitigated to a degree by removing systems with shorter initial periods, suggesting that the observed sample of overcontact systems may originate from binary systems with longer initial orbital periods.
△ Less
Submitted 2 August, 2022;
originally announced August 2022.
-
SPAMMS: applications and use cases for the 3D spectroscopic analysis technique to study deformed massive stars
Authors:
Michael Abdul-Masih
Abstract:
Whether it be due to rapid rotation or binary interactions, deviations from spherical symmetry are common in massive stars. These deviations from spherical symmetry are known to cause non-uniform distributions of various parameters across the surface including temperature, which can drive internal mixing processes within the envelopes of these massive stars. Despite how common these 3D distortions…
▽ More
Whether it be due to rapid rotation or binary interactions, deviations from spherical symmetry are common in massive stars. These deviations from spherical symmetry are known to cause non-uniform distributions of various parameters across the surface including temperature, which can drive internal mixing processes within the envelopes of these massive stars. Despite how common these 3D distortions are, they are often neglected in spectroscopic analyses. We present a new spectral analysis code called SPAMMS (Spectroscopic PAtch Model for Massive Stars) specifically designed to analyze non-spherical systems. We discuss how the code works and discuss its assumptions. Furthermore, we demonstrate how SPAMMS can be applied to a variety of different types of systems and we show how it can model 3D effects in a way that current analysis techniques are not able to.
△ Less
Submitted 2 August, 2022;
originally announced August 2022.
-
Identifying quiescent compact objects in massive Galactic single-lined spectroscopic binaries
Authors:
L. Mahy,
H. Sana,
T. Shenar,
K. Sen,
N. Langer,
P. Marchant,
M. Abdul-Masih,
G. Banyard,
J. Bodensteiner,
D. M. Bowman,
K. Dsilva,
M. Fabry,
C. Hawcroft,
S. Janssens,
T. Van Reeth,
C. Eldridge
Abstract:
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, s…
▽ More
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.
△ Less
Submitted 15 July, 2022;
originally announced July 2022.
-
An X-ray quiet black hole born with a negligible kick in a massive binary within the Large Magellanic Cloud
Authors:
Tomer Shenar,
Hugues Sana,
Laurent Mahy,
Kareem El-Badry,
Pablo Marchant,
Norbert Langer,
Calum Hawcroft,
Matthias Fabry,
Koushik Sen,
Leonardo A. Almeida,
Michael Abdul-Masih,
Julia Bodensteiner,
Paul A. Crowther,
Mark Gieles,
Mariusz Gromadzki,
Vincent Henault-Brunet,
Artemio Herrero,
Alex de Koter,
Patryk Iwanek,
Szymon Kozłowski,
Daniel J. Lennon,
Jesus Maız Apellaniz,
Przemysław Mroz,
Anthony F. J. Moffat,
Annachiara Picco
, et al. (13 additional authors not shown)
Abstract:
Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray faint binary in the Large Magellanic Cloud. With an orbital period of 10.4-d, it comprises an O-ty…
▽ More
Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray faint binary in the Large Magellanic Cloud. With an orbital period of 10.4-d, it comprises an O-type star of 25 solar masses and an unseen companion of at least nine solar masses. Our spectral analysis excludes a non-degenerate companion at a 5-sigma confidence level. The minimum companion mass implies that it is a black hole. No other X-ray quiet black hole is unambiguously known outside our Galaxy. The (near-)circular orbit and kinematics of VFTS 243 imply that the collapse of the progenitor into a black hole was associated with little or no ejected material or black-hole kick. Identifying such unique binaries substantially impacts the predicted rates of gravitational-wave detections and properties of core-collapse supernovae across the Cosmos.
△ Less
Submitted 15 July, 2022;
originally announced July 2022.
-
Detecting Stripped Stars While Searching for Quiescent Black Holes
Authors:
J. Bodensteiner,
M. Heida,
M. Abdul-Masih,
D. Baade,
G. Banyard,
D. M. Bowman,
M. Fabry,
A. Frost,
L. Mahy,
P. Marchant,
A. Mérand,
M. Reggiani,
Th. Rivinius,
H. Sana,
F. Selman,
T. Shenar
Abstract:
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 bl…
▽ More
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.
△ Less
Submitted 1 July, 2022;
originally announced July 2022.
-
HR 6819 is a binary system with no black hole -- revisiting the source with infrared interferometry and optical integral field spectroscopy
Authors:
A. J. Frost,
J. Bodensteiner,
Th. Rivinius,
D. Baade,
A. Merand,
F. Selman,
M. Abdul-Masih,
G. Banyard,
E. Bordier,
K. Dsilva,
C. Hawcroft,
L. Mahy,
M. Reggiani,
T. Shenar,
M. Cabezas,
P. Hadrava,
M. Heida,
R. Klement,
H. Sana
Abstract:
Two scenarios have been proposed to match the existing observational constraints of the object HR 6819. The system could consist of a close inner B-type giant plus a black hole (BH) binary with an additional Be companion in a wide orbit. Alternatively, it could be a binary composed of a stripped B star and a Be star in a close orbit. Either scenario makes HR 6819 a cornerstone object as the stella…
▽ More
Two scenarios have been proposed to match the existing observational constraints of the object HR 6819. The system could consist of a close inner B-type giant plus a black hole (BH) binary with an additional Be companion in a wide orbit. Alternatively, it could be a binary composed of a stripped B star and a Be star in a close orbit. Either scenario makes HR 6819 a cornerstone object as the stellar BH closest to Earth, or as an example of an important transitional, non-equilibrium phase for Be stars with solid evidence for its nature. We aimed to distinguish between the two scenarios for HR 6819. Both models predict two luminous stars but with very different angular separations and orbital motions. Therefore, the presence of bright sources in the 1-100 milliarcsec (mas) regime is a key diagnostic for determining the nature of the HR 6819 system. We obtained new high-angular resolution data with VLT/MUSE and VLTI/GRAVITY of HR 6819. The MUSE data are sensitive to bright companions at large scales, whilst the interferometric GRAVITY data are sensitive down to separations on mas scales and large magnitude differences. The MUSE observations reveal no bright companion at large separations and the GRAVITY observations indicate the presence of a stellar companion at an angular separation of ~1.2 mas that moves on the plane of the sky over a timescale compatible with the known spectroscopic 40-day period. We conclude that HR 6819 is a binary system and that no BH is present in the system. The unique nature of HR 6819, and its proximity to Earth make it an ideal system for quantitatively characterising the immediate outcome of binary interaction and probing how Be stars form.
△ Less
Submitted 2 March, 2022;
originally announced March 2022.
-
The R136 star cluster dissected with Hubble Space Telescope/STIS. III. The most massive stars and their clumped winds
Authors:
Sarah A. Brands,
Alex de Koter,
Joachim M. Bestenlehner,
Paul A. Crowther,
Jon O. Sundqvist,
Joachim Puls,
Saida M. Caballero-Nieves,
Michael Abdul-Masih,
Florian A. Driessen,
Miriam García,
Sam Geen,
Götz Gräfener,
Calum Hawcroft,
Lex Kaper,
Zsolt Keszthelyi,
Norbert Langer,
Hugues Sana,
Fabian R. N. Schneider,
Tomer Shenar,
Jorick S. Vink
Abstract:
Context: The star cluster R136 inside the LMC hosts a rich population of massive stars, including the most massive stars known. The strong stellar winds of these very luminous stars impact their evolution and the surrounding environment. We currently lack detailed knowledge of the wind structure that is needed to quantify this impact. Aims: To observationally constrain the stellar and wind propert…
▽ More
Context: The star cluster R136 inside the LMC hosts a rich population of massive stars, including the most massive stars known. The strong stellar winds of these very luminous stars impact their evolution and the surrounding environment. We currently lack detailed knowledge of the wind structure that is needed to quantify this impact. Aims: To observationally constrain the stellar and wind properties of the massive stars in R136, in particular the parameters related to wind clumping. Methods: We simultaneously analyse optical and UV spectroscopy of 53 O-type and 3 WNh-stars using the FASTWIND model atmosphere code and a genetic algorithm. The models account for optically thick clumps and effects related to porosity and velocity-porosity, as well as a non-void interclump medium. Results: We obtain stellar parameters, surface abundances, mass-loss rates, terminal velocities and clumping characteristics and compare these to theoretical predictions and evolutionary models. The clumping properties include the density of the interclump medium and the velocity-porosity of the wind. For the first time, these characteristics are systematically measured for a wide range of effective temperatures and luminosities. Conclusions: We confirm a cluster age of 1.0-2.5 Myr and derive an initial stellar mass of $\geq 250 {\rm M}_\odot$ for the most massive star in our sample, R136a1. The winds of our sample stars are highly clumped, with an average clumping factor of $f_{\rm cl}=29\pm15$. We find tentative trends in the wind-structure parameters as a function of mass-loss rate, suggesting that the winds of stars with higher mass-loss rates are less clumped. We compare several theoretical predictions to the observed mass-loss rates and terminal velocities and find that none satisfactorily reproduces both quantities. The prescription of Krtička & Kubát (2018) matches best the observed mass-loss rates.
△ Less
Submitted 22 February, 2022;
originally announced February 2022.
-
Planet Hunters TESS IV: A massive, compact hierarchical triple star system TIC 470710327
Authors:
Nora L. Eisner,
Cole Johnston,
Silvia Toonen,
Abigail J. Frost,
Soetkin Janssens,
Chris J. Lintott,
Suzanne Aigrain,
Hugues Sana,
Michael Abdul-Masih,
Karla Z. Arellano-Córdova,
Paul G. Beck,
Emma Bordier,
Emily Canon,
Ana Escorza,
Mattias Fabry,
Lars Hermansson,
Steve Howell,
Grant Miller,
Shreeya Sheyte,
Safaa Alhassan,
Elisabeth M. L. Baeten,
Frank Barnet,
Stewart. J. Bean,
Mikael Bernau,
David M. Bundy
, et al. (15 additional authors not shown)
Abstract:
We report the discovery and analysis of a massive, compact, hierarchical triple system (TIC 470710327) initially identified by citizen scientists in data obtained by NASA's Transiting Exoplanet Survey Satellite (TESS). Spectroscopic follow-up observations obtained with the HERMES spectrograph, combined with eclipse timing variations (ETVs), confirm that the system is comprised of three OB stars, w…
▽ More
We report the discovery and analysis of a massive, compact, hierarchical triple system (TIC 470710327) initially identified by citizen scientists in data obtained by NASA's Transiting Exoplanet Survey Satellite (TESS). Spectroscopic follow-up observations obtained with the HERMES spectrograph, combined with eclipse timing variations (ETVs), confirm that the system is comprised of three OB stars, with a compact 1.10 d eclipsing binary and a non-eclipsing tertiary on a 52.04 d orbit. Dynamical modelling of the system (from radial velocity and ETVs) reveal a rare configuration wherein the tertiary star (O9.5-B0.5V; 14-17 M$_{\odot}$) is more massive than the combined mass of the inner binary (10.9-13.2 M$_{\odot}$). Given the high mass of the tertiary, we predict that this system will undergo multiple phases of mass transfer in the future, and likely end up as a double neutron star gravitational wave progenitor or an exotic Thorne-Zytkow object. Further observational characterisation of this system promises constraints on both formation scenarios of massive stars as well as their exotic evolutionary end-products.
△ Less
Submitted 14 February, 2022;
originally announced February 2022.
-
Binary-object spectral-synthesis in 3D (BOSS-3D) -- Modelling H-alpha emission in the enigmatic multiple system LB-1
Authors:
L. Hennicker,
N. D. Kee,
T. Shenar,
J. Bodensteiner,
M. Abdul-Masih,
I. El Mellah,
H. Sana,
J. O. Sundqvist
Abstract:
Context: To quantitatively decode the information stored within an observed spectrum, detailed modelling of the physical state and accurate radiative transfer solution schemes are required. In the analysis of stellar spectra, the numerical model often needs to account for binary companions and 3D structures in the stellar envelopes. The enigmatic binary (or multiple) system LB-1 constitutes a perf…
▽ More
Context: To quantitatively decode the information stored within an observed spectrum, detailed modelling of the physical state and accurate radiative transfer solution schemes are required. In the analysis of stellar spectra, the numerical model often needs to account for binary companions and 3D structures in the stellar envelopes. The enigmatic binary (or multiple) system LB-1 constitutes a perfect example of such a complex multi-D problem. Aims: To improve our understanding of the LB-1 system, we directly modelled the phase-dependent H-alpha line profiles of this system. To this end, we developed a multi-purpose binary-object spectral-synthesis code in 3D (BOSS-3D). Methods: BOSS-3D calculates synthetic line profiles for a given state of the circumstellar material. The standard pz-geometry commonly used for single stars is extended by defining individual coordinate systems for each involved object and by accounting for the appropriate coordinate transformations. The code is then applied to the LB-1 system, considering two main hypotheses, a binary containing a stripped star and Be star, or a B star and a black hole with a disc. Results: Comparing these two scenarios, neither model can reproduce the detailed phase-dependent shape of the H-alpha line profiles. A satisfactory match with the observations, however, is obtained by invoking a disc around the primary object in addition to the Be-star disc or the black-hole accretion disc. Conclusions: The developed code can be used to model synthetic line profiles for a wide variety of binary systems, ranging from transit spectra of planetary atmospheres, to post-asymptotic giant branch binaries including circumstellar and circumbinary discs and massive-star binaries with stellar winds and disc systems. For the LB-1 system, our modelling provides strong evidence that each object in the system contains a disc-like structure.
△ Less
Submitted 30 November, 2021;
originally announced November 2021.
-
TESS Eclipsing Binary Stars. I. Short cadence observations of 4584 eclipsing binaries in Sectors 1-26
Authors:
Andrej Prsa,
Angela Kochoska,
Kyle E. Conroy,
Nora Eisner,
Daniel R. Hey,
Luc IJspeert,
Ethan Kruse,
Scott W. Fleming,
Cole Johnston,
Martti H. Kristiansen,
Daryll LaCourse,
Danielle Mortensen,
Joshua Pepper,
Keivan G. Stassun,
Guillermo Torres,
Michael Abdul-Masih,
Joheen Chakraborty,
Robert Gagliano,
Zhao Guo,
Kelly Hambleton,
Kyeongsoo Hong,
Thomas Jacobs,
David Jones,
Veselin Kostov,
Jae Woo Lee
, et al. (22 additional authors not shown)
Abstract:
In this paper we present a catalog of 4584 eclipsing binaries observed during the first two years (26 sectors) of the TESS survey. We discuss selection criteria for eclipsing binary candidates, detection of hither-to unknown eclipsing systems, determination of the ephemerides, the validation and triage process, and the derivation of heuristic estimates for the ephemerides. Instead of keeping to th…
▽ More
In this paper we present a catalog of 4584 eclipsing binaries observed during the first two years (26 sectors) of the TESS survey. We discuss selection criteria for eclipsing binary candidates, detection of hither-to unknown eclipsing systems, determination of the ephemerides, the validation and triage process, and the derivation of heuristic estimates for the ephemerides. Instead of keeping to the widely used discrete classes, we propose a binary star morphology classification based on a dimensionality reduction algorithm. Finally, we present statistical properties of the sample, we qualitatively estimate completeness, and discuss the results. The work presented here is organized and performed within the TESS Eclipsing Binary Working Group, an open group of professional and citizen scientists; we conclude by describing ongoing work and future goals for the group. The catalog is available from http://tessEBs.villanova.edu and from MAST.
△ Less
Submitted 25 October, 2021;
originally announced October 2021.
-
Empirical mass-loss rates and clumping properties of Galactic early-type O supergiants
Authors:
C. Hawcroft,
H. Sana,
L. Mahy,
J. O. Sundqvist,
M. Abdul-Masih,
J. C. Bouret,
S. A. Brands,
A. de Koter,
F. A. Driessen,
J. Puls
Abstract:
We investigate the impact of optically thick clumping on stellar wind diagnostics in O supergiants and constrain wind parameters associated with porosity in velocity space. This is the first time the effects of optically thick clumping have been investigated for a sample of massive hot stars, using models including a full optically thick clumping description. We re-analyse spectroscopic observatio…
▽ More
We investigate the impact of optically thick clumping on stellar wind diagnostics in O supergiants and constrain wind parameters associated with porosity in velocity space. This is the first time the effects of optically thick clumping have been investigated for a sample of massive hot stars, using models including a full optically thick clumping description. We re-analyse spectroscopic observations of a sample of eight O supergiants. Using a genetic algorithm wrapper around the NLTE atmosphere code FASTWIND we obtain simultaneous fits to optical and UV spectra and determine photospheric and wind properties and surface abundances. We provide empirical constraints on a number of wind parameters including the clumping factors, mass-loss rates and terminal wind velocities. Additionally, we establish the first systematic empirical constraints on velocity filling factors and interclump densities. These parameters describe clump distribution in velocity-space and density of the interclump medium in physical-space, respectively. We observe a mass-loss rate reduction of a factor of 3.6 compared to theoretical predictions from Vink et al. (2000), and mass-loss rates within a factor 1.4 of predictions from Björklund et al. (2021). We confirm that including optically thick clumping allows simultaneous fitting of recombination lines and resonance lines, including the unsaturated UV phosphorus lines (Pv 1118-1128), without reducing the phosphorus abundance. We find that, on average, half of the wind velocity field is covered by dense clumps. We also find that these clumps are 25 times denser than the average wind, and that the interclump medium is 3-10 times less dense than the mean wind. The former result agrees well with theoretical predictions, the latter suggests that lateral filling-in of radially compressed gas might be critical for setting the scale of the rarefied interclump matter.
△ Less
Submitted 18 August, 2021;
originally announced August 2021.
-
Constraining the Overcontact Phase in Massive Binary Evolution I. Mixing in V382 Cyg, VFTS 352, and OGLE SMC-SC10 108086
Authors:
Michael Abdul-Masih,
Hugues Sana,
Calum Hawcroft,
Leonardo A. Almeida,
Sarah A. Brands,
Selma E. de Mink,
Stephen Justham,
Norbert Langer,
Laurent Mahy,
Pablo Marchant,
Athira Menon,
Joachim Puls,
Jon Sundqvist
Abstract:
As potential progenitors of several exotic phenomena including gravitational wave sources, magnetic stars, and Be stars, close massive binary systems probe a crucial area of the parameter space in massive star evolution. Despite the importance of these systems, large uncertainties regarding the nature and efficiency of the internal mixing mechanisms still exist. In this work, we aim to provide rob…
▽ More
As potential progenitors of several exotic phenomena including gravitational wave sources, magnetic stars, and Be stars, close massive binary systems probe a crucial area of the parameter space in massive star evolution. Despite the importance of these systems, large uncertainties regarding the nature and efficiency of the internal mixing mechanisms still exist. In this work, we aim to provide robust observational constraints on the internal mixing processes by spectroscopically analyzing a sample of three massive overcontact binaries at different metallicities. Using optical phase-resolved spectroscopic data, we perform an atmosphere analysis using more traditional 1D techniques and using state-of-the-art 3D techniques. We compare and contrast the assumptions and results of each technique and investigate how the assumptions affect the final derived atmospheric parameters. We find that in all three cases, both components of system are highly overluminous indicating either efficient internal mixing of helium or previous non-conservative mass transfer. However, we do not find strong evidence of helium or CNO surface abundance changes usually associated with mixing. Additionally, we find that in unequal mass systems, the measured effective temperature and luminosity of the less massive component places it very close to the more massive component on the Hertzsprung-Russell diagram. These results were obtained independently using both of the techniques mentioned above, which suggests that these measurements are robust. The observed discrepancies between the temperature and the surface abundance measurements when compared to theoretical expectations indicate that unaccounted for additional physical mechanisms may be at play.
△ Less
Submitted 16 April, 2021; v1 submitted 15 April, 2021;
originally announced April 2021.
-
Characterization of the variability in the O+B eclipsing binary HD 165246
Authors:
Cole Johnston,
Nicolas Aimar,
Michael Abdul-Masih,
Dominic. M. Bowman,
Tim R. White,
Calum Hawcroft,
Hugues Sana,
Sanjay Sekaran,
Karan Dsilva,
Andrew Tkachenko,
Conny Aerts
Abstract:
O-stars are known to experience a wide range of variability mechanisms originating at both their surface and their near-core regions. Characterization and understanding of this variability and its potential causes are integral for evolutionary calculations. We use a new extensive high-resolution spectroscopic data set to characterize the variability observed in both the spectroscopic and space-bas…
▽ More
O-stars are known to experience a wide range of variability mechanisms originating at both their surface and their near-core regions. Characterization and understanding of this variability and its potential causes are integral for evolutionary calculations. We use a new extensive high-resolution spectroscopic data set to characterize the variability observed in both the spectroscopic and space-based photometric observations of the O+B eclipsing binary HD~165246. We present an updated atmospheric and binary solution for the primary component, involving a high level of microturbulence ($13_{-1.3}^{+1.0}\,$km\,s$^{-1}$) and a mass of $M_1=23.7_{-1.4}^{+1.1}$~M$_{\odot}$, placing it in a sparsely explored region of the Hertzsprung-Russell diagram. Furthermore, we deduce a rotational frequency of $0.690\pm 0.003\,$d$^{-1}$ from the combined photometric and line-profile variability, implying that the primary rotates at 40\% of its critical Keplerian rotation rate. We discuss the potential explanations for the overall variability observed in this massive binary, and discuss its evolutionary context.
△ Less
Submitted 16 February, 2021;
originally announced February 2021.
-
Detailed evolutionary models of massive contact binaries: I. Model grids and synthetic populations for the Magellanic Clouds
Authors:
Athira Menon,
Norbert Langer,
Selma E. de Mink,
Stephen Justham,
Koushik Sen,
Dorottya Szécsi,
Alex de Koter,
Michael Abdul-Masih,
Hugues Sana,
Laurent Mahy,
Pablo Marchant
Abstract:
The majority of close massive binary stars with initial periods of a few days experience a contact phase, in which both stars overflow their Roche lobes simultaneously. We perform the first dedicated study of the evolution of massive contact binaries and provide a comprehensive prediction of their observed properties. We compute 2790 detailed binary models for the Large and Small Magellanic Clouds…
▽ More
The majority of close massive binary stars with initial periods of a few days experience a contact phase, in which both stars overflow their Roche lobes simultaneously. We perform the first dedicated study of the evolution of massive contact binaries and provide a comprehensive prediction of their observed properties. We compute 2790 detailed binary models for the Large and Small Magellanic Clouds each, assuming a conservative mass transfer. The initial parameter space for both grids span total masses from 20 to 80$\,\textrm{M}_{\odot}$, orbital periods of 0.6 to 2 days and mass ratios of 0.6 to 1.0. We find that models that remain in contact over nuclear timescales evolve towards equal masses, echoing the mass ratios of their observed counterparts. Ultimately, the fate of our nuclear-timescale models is to merge on the main sequence. Our predicted period-mass ratio distributions of O-type contact binaries are similar for both galaxies, and we expect 10 such systems together in both Magellanic Clouds. While we can largely reproduce the observed distribution, we over-estimate the population of equal-mass contact binaries. This situation is somewhat remedied if we also account for binaries that are approaching contact. Our theoretical distributions work particularly well for contact binaries with periods $<$2 days and total masses $\lessapprox45\,\textrm{M}_{\odot}$. We expect stellar winds, non-conservative mass transfer and envelope inflation to have played a role in the formation of the more massive and longer-period contact binaries.
△ Less
Submitted 4 August, 2021; v1 submitted 26 November, 2020;
originally announced November 2020.
-
Tango of celestial dancers: A sample of detached eclipsing binary systems containing g-mode pulsating components. A case study of KIC9850387
Authors:
S. Sekaran,
A. Tkachenko,
M. Abdul-Masih,
A. Prsa,
C. Johnston,
D. Huber,
S. J. Murphy,
G. Banyard,
A. W. Howard,
H. Isaacson,
D. M. Bowman,
C. Aerts
Abstract:
Eclipsing binary systems with components that pulsate in gravity modes ($g$ modes) allow for simultaneous and independent constraints of the chemical mixing profiles of stars. We aim to assemble a sample of $g$-mode pulsators in detached eclipsing binaries with the purpose of finding good candidates for future evolutionary and asteroseismic modelling. In addition, we present a case study of the ec…
▽ More
Eclipsing binary systems with components that pulsate in gravity modes ($g$ modes) allow for simultaneous and independent constraints of the chemical mixing profiles of stars. We aim to assemble a sample of $g$-mode pulsators in detached eclipsing binaries with the purpose of finding good candidates for future evolutionary and asteroseismic modelling. In addition, we present a case study of the eclipsing binary KIC9850387, identified as our most promising candidate. We selected all of the detached eclipsing binaries in the Kepler eclipsing binary catalogue and performed a visual inspection to determine the presence and density of $g$ modes, and the presence of $g$-mode period-spacing patterns in their frequency spectra. We then characterised our sample based on their $g$-mode pulsational parameters and binary and atmospheric parameters. A spectroscopic follow-up of KIC9850387 was then performed. We then performed spectral disentangling followed by atmospheric modelling and abundance analysis for the primary star. We utilised an iterative approach to simultaneously optimise the pulsational and eclipse models, and subsequently performed an analysis of the pressure- ($p$-) and $g$-mode pulsational frequencies. We compiled a sample of 93 Kepler eclipsing binary stars with $g$-mode pulsating components and identified clear $g$-mode period-spacing patterns in the frequency spectra of seven of these systems. We found that the $g$-mode pulsational parameters and the binary and atmospheric parameters of our sample are weakly correlated at best. We find that the eclipsing binary KIC9850387 is a double-lined spectroscopic binary in a near-circular orbit with an intermediate-mass primary, and a solar-like secondary. We find $\ell=1$ and $\ell=2$ period-spacing patterns in the frequency spectrum of KIC9850387 spanning more than ten radial orders each.
△ Less
Submitted 5 October, 2020;
originally announced October 2020.
-
Physics of Eclipsing Binaries. V. General Framework for Solving the Inverse Problem
Authors:
Kyle E Conroy,
Angela Kochoska,
Daniel Hey,
Herbert Pablo,
Kelly M Hambleton,
David Jones,
Joseph Giammarco,
Michael Abdul-Masih,
Andrej Prsa
Abstract:
PHOEBE 2 is a Python package for modeling the observables of eclipsing star systems, but until now has focused entirely on the forward-model -- that is, generating a synthetic model given fixed values of a large number of parameters describing the system and the observations. The inverse problem, obtaining orbital and stellar parameters given observational data, is more complicated and computation…
▽ More
PHOEBE 2 is a Python package for modeling the observables of eclipsing star systems, but until now has focused entirely on the forward-model -- that is, generating a synthetic model given fixed values of a large number of parameters describing the system and the observations. The inverse problem, obtaining orbital and stellar parameters given observational data, is more complicated and computationally expensive as it requires generating a large set of forward-models to determine which set of parameters and uncertainties best represent the available observational data. The process of determining the best solution and also of obtaining reliable and robust uncertainties on those parameters often requires the use of multiple algorithms, including both optimizers and samplers. Furthermore, the forward-model of PHOEBE has been designed to be as physically robust as possible, but is computationally expensive compared to other codes. It is useful, therefore, to use whichever code is most efficient given the reasonable assumptions for a specific system, but learning the intricacies of multiple codes presents a barrier to doing this in practice. Here we present the 2.3 release of PHOEBE (publicly available from http://phoebe-project.org) which introduces a general framework for defining and handling distributions on parameters, and utilizing multiple different estimation, optimization, and sampling algorithms. The presented framework supports multiple forward-models, including the robust model built into PHOEBE itself.
△ Less
Submitted 5 October, 2020; v1 submitted 30 June, 2020;
originally announced June 2020.
-
Is HR 6819 a triple system containing a black hole? -- An alternative explanation
Authors:
J. Bodensteiner,
T. Shenar,
L. Mahy,
M. Fabry,
P. Marchant,
M. Abdul-Masih,
G. Banyard,
D. M. Bowman,
K. Dsilva,
A. J. Frost,
C. Hawcroft,
M. Reggiani,
H. Sana
Abstract:
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 t…
▽ More
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.
△ Less
Submitted 18 June, 2020;
originally announced June 2020.
-
The 'hidden' companion in LB-1 unveiled by spectral disentangling
Authors:
T. Shenar,
J. Bodensteiner,
M. Abdul-Masih,
M. Fabry,
L. Mahy,
P. Marchant,
G. Banyard,
D. M. Bowman,
K. Dsilva,
C. Hawcroft,
M. Reggiani,
H. Sana
Abstract:
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 sequenc…
▽ More
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.
△ Less
Submitted 11 July, 2020; v1 submitted 27 April, 2020;
originally announced April 2020.
-
Spectroscopic patch model for massive stars using PHOEBE II and FASTWIND
Authors:
Michael Abdul-Masih,
Hugues Sana,
Kyle E. Conroy,
Jon Sundqvist,
Andrej Prša,
Angela Kochoska,
Joachim Puls
Abstract:
Massive stars play an important role in the mechanical and chemical evolution of galaxies. Understanding the internal processes of these stars is vital to our understanding of their evolution and eventual end products. Deformations from spherical geometry are common for massive stars; however, the tools that are currently available for the study of these systems are almost exclusively one-dimensio…
▽ More
Massive stars play an important role in the mechanical and chemical evolution of galaxies. Understanding the internal processes of these stars is vital to our understanding of their evolution and eventual end products. Deformations from spherical geometry are common for massive stars; however, the tools that are currently available for the study of these systems are almost exclusively one-dimensional. We present a new spectroscopic analysis tool tailored for massive stars that deviate from spherical symmetry. This code (entitled SPAMMS) is a spectroscopic patch model that takes the three-dimensional surface geometry of the system into account to produce spectral profiles at given phases and orientations. In using the Wilson-Devinney-like code PHOEBE in combination with the non local thermodynamic equilibrium (NLTE) radiative transfer code FASTWIND, we created a three-dimensional mesh that represents the surface geometry of our system and we assigned FASTWIND emergent intensity line profiles to each mesh triangle, which take the local parameters such as temperature, surface gravity, and radius into account. These line profiles were then integrated across the visible surface, where their flux contribution and radial velocity are taken into account, thus returning a final line profile for the visible surface of the system at a given phase. We demonstrate that SPAMMS can accurately reproduce the morphology of observed spectral line profiles for overcontact systems. Additionally, we show how line profiles of rapidly-rotating single stars differ when taking rotational distortion into account, and the effects that these can have on the determined parameters. Finally, we demonstrate the code's ability to reproduce the Rossiter-Mclaughlin and Struve-Sahade effects.
△ Less
Submitted 19 March, 2020;
originally announced March 2020.
-
The Tarantula Massive Binary Monitoring: III. Atmosphere analysis of double-lined spectroscopic systems
Authors:
L. Mahy,
H. Sana,
M. Abdul-Masih,
L. A. Almeida,
N. Langer,
T. Shenar,
A. de Koter,
S. E. de Mink,
S. de Wit,
N. J. Grin,
C. J. Evans,
A. F. J. Moffat,
F. R. N. Schneider,
R. Barbá,
J. S. Clark,
P. Crowther,
G. Gräfener,
D. J. Lennon,
F. Tramper,
J. S. Vink
Abstract:
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 a…
▽ More
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.
△ Less
Submitted 14 December, 2019;
originally announced December 2019.
-
No signature of the orbital motion of a putative 70 solar mass black hole in LB-1
Authors:
Michael Abdul-Masih,
Gareth Banyard,
Julia Bodensteiner,
Dominic M. Bowman,
Karan Dsilva,
Matthias Fabry,
Calum Hawcroft,
Laurent Mahy,
Pablo Marchant,
Gert Raskin,
Maddalena Reggiani,
Hugues Sana,
Tomer Shenar,
Andrew Tkachenko,
Hans Van Winckel,
Lore Vermeylen
Abstract:
Liu et al. (2019) recently reported the detection of a 68 [+11/-13] solar mass (Msun) black hole (BH) paired with an 8.2 [+0.9/-1.2] Msun B-type sub-giant star in the 78.9-day spectroscopic binary system LB-1. Such a black hole is over twice as massive as any other known stellar-mass black hole with non-compact companions2 and its mass approaches those that result from BH-BH coalescences that are…
▽ More
Liu et al. (2019) recently reported the detection of a 68 [+11/-13] solar mass (Msun) black hole (BH) paired with an 8.2 [+0.9/-1.2] Msun B-type sub-giant star in the 78.9-day spectroscopic binary system LB-1. Such a black hole is over twice as massive as any other known stellar-mass black hole with non-compact companions2 and its mass approaches those that result from BH-BH coalescences that are detected by gravitational wave interferometers. Its presence in a solar-like metallicity environment challenges conventional theories of massive binary evolution, stellar winds and core-collapse supernovae, so that more exotic scenarios seem to be needed to explain the existence and properties of LB-1. Here, we show that the observational diagnostics used to derive the BH mass results from the orbital motion of the B-type star, not that of the BH. As a consequence, no evidence for a massive BH remains in the data, therefore solving the existing tension with formation models of such a massive BH at solar metallicity and with theories of massive star evolution in general.
△ Less
Submitted 9 December, 2019;
originally announced December 2019.
-
Clues on the Origin and Evolution of Massive Contact Binaries: Atmosphere Analysis of VFTS 352
Authors:
Michael Abdul-Masih,
Hugues Sana,
Jon Sundqvist,
Laurent Mahy,
Athira Menon,
Leonardo A. Almeida,
Alex De Koter,
Selma E. de Mink,
Stephen Justham,
Norbert Langer,
Joachim Puls,
Tomer Shenar,
Frank Tramper
Abstract:
The massive O4.5 V + O5.5 V binary VFTS 352 in the Tarantula nebula is one of the shortest-period and most massive overcontact binaries known. Recent theoretical studies indicate that some of these systems could ultimately lead to the formation of gravitational waves via black hole binary mergers through the chemically homogeneous evolution pathway. By analyzing ultraviolet-optical phase-resolved…
▽ More
The massive O4.5 V + O5.5 V binary VFTS 352 in the Tarantula nebula is one of the shortest-period and most massive overcontact binaries known. Recent theoretical studies indicate that some of these systems could ultimately lead to the formation of gravitational waves via black hole binary mergers through the chemically homogeneous evolution pathway. By analyzing ultraviolet-optical phase-resolved spectroscopic data, we aim to constrain atmospheric and wind properties that could be later used to confront theoretical predictions from binary evolution. In particular, surface abundances are powerful diagnostics of the evolutionary status, mass transfer and the internal mixing processes. From a set of 32 VLT/FLAMES visual and 8 HST/COS ultraviolet spectra, we used spectral disentangling to separate the primary and secondary components. Using a genetic algorithm wrapped around the NLTE model atmosphere and spectral synthesis code FASTWIND, we perform an 11-parameter optimization to derive the atmospheric and wind parameters of both components, including the surface abundances of He, C, N, O and Si. We find that both components are hotter than expected compared to single-star evolutionary models indicating that additional mixing processes may be at play. However the derived chemical abundances do not show significant indications of mixing when adopting baseline values typical for the system environment.
△ Less
Submitted 3 June, 2019;
originally announced June 2019.
-
Barium & related stars and their white-dwarf companions II. Main-sequence and subgiant stars
Authors:
A. Escorza,
D. Karinkuzhi,
A. Jorissen,
L. Siess,
H. Van Winckel,
D. Pourbaix,
C. Johnston,
B. Miszalski,
G-M. Oomen,
M. Abdul-Masih,
H. M. J. Boffin,
P. North,
R. Manick,
S. Shetye,
J. Mikołajewska
Abstract:
Barium (Ba) dwarfs and CH subgiants are the less-evolved analogues of Ba and CH giants. They are F- to G-type main-sequence stars polluted with heavy elements by a binary companion when the latter was on the Asymptotic Giant Branch (AGB). This companion is now a white dwarf that in most cases cannot be directly detected. We present a large systematic study of 60 objects classified as Ba dwarfs or…
▽ More
Barium (Ba) dwarfs and CH subgiants are the less-evolved analogues of Ba and CH giants. They are F- to G-type main-sequence stars polluted with heavy elements by a binary companion when the latter was on the Asymptotic Giant Branch (AGB). This companion is now a white dwarf that in most cases cannot be directly detected. We present a large systematic study of 60 objects classified as Ba dwarfs or CH subgiants. Combining radial-velocity measurements from HERMES and SALT high-resolution spectra with radial-velocity data from CORAVEL and CORALIE, we determine the orbital parameters of 27 systems. We also derive their masses by comparing their location in the Hertzsprung-Russell diagram with evolutionary models. We confirm that Ba dwarfs and CH subgiants are not at different evolutionary stages and have similar metallicities, despite their different names. Additionally, Ba giants appear significantly more massive than their main-sequence analogues. This is likely due to observational biases against the detection of hotter main-sequence post-mass-transfer objects. Combining our spectroscopic orbits with the Hipparcos astrometric data, we derive the orbital inclinations and the mass of the WD companion for four systems. Since this cannot be done for all systems in our sample yet (but should be with upcoming Gaia data releases), we also analyse the mass-function distribution of our binaries. We can model this distribution with very narrow mass distributions for the two components and random orbital orientation on the sky. Finally, based on BINSTAR evolutionary models, we suggest that the orbital evolution of low-mass Ba systems can be affected by a second phase of interaction along the Red Giant Branch of the Ba star, impacting on the eccentricities and periods of the giants.
△ Less
Submitted 8 April, 2019;
originally announced April 2019.
-
Kepler Eclipsing Binary Stars. VIII. Identification of False Positive Eclipsing Binaries and Re-extraction of New Light Curves
Authors:
Michael Abdul-Masih,
Andrej Prsa,
Kyle Conroy,
Steven Bloemen,
Tabetha Boyajian,
Laurance R. Doyle,
Cole Johnston,
Veselin Kostov,
David W. Latham,
Gal Matijevic,
Avi Shporer,
John Southworth
Abstract:
The Kepler Mission has provided unprecedented, nearly continuous photometric data of $\sim$200,000 objects in the $\sim$105 deg$^{2}$ field of view from the beginning of science operations in May of 2009 until the loss of the second reaction wheel in May of 2013. The Kepler Eclipsing Binary Catalog contains information including but not limited to ephemerides, stellar parameters and analytical app…
▽ More
The Kepler Mission has provided unprecedented, nearly continuous photometric data of $\sim$200,000 objects in the $\sim$105 deg$^{2}$ field of view from the beginning of science operations in May of 2009 until the loss of the second reaction wheel in May of 2013. The Kepler Eclipsing Binary Catalog contains information including but not limited to ephemerides, stellar parameters and analytical approximation fits for every known eclipsing binary system in the Kepler Field of View. Using Target Pixel level data collected from Kepler in conjunction with the Kepler Eclipsing Binary Catalog, we identify false positives among eclipsing binaries, i.e. targets that are not eclipsing binaries themselves, but are instead contaminated by eclipsing binary sources nearby on the sky and show eclipsing binary signatures in their light curves. We present methods for identifying these false positives and for extracting new light curves for the true source of the observed binary signal. For each source, we extract three separate light curves for each quarter of available data by optimizing the signal-to-noise ratio, the relative percent eclipse depth and the flux eclipse depth. We present 289 new eclipsing binaries in the Kepler Field of View that were not targets for observation, and these have been added to the Catalog. An online version of this Catalog with downloadable content and visualization tools is maintained at http://keplerEBs.villanova.edu.
△ Less
Submitted 18 February, 2016;
originally announced February 2016.
-
Kepler Eclipsing Binary Stars. VII. The Catalog of Eclipsing Binaries Found in the Entire Kepler Data-Set
Authors:
Brian Kirk,
Kyle Conroy,
Andrej Prša,
Michael Abdul-Masih,
Angela Kochoska,
Gal Matijevič,
Kelly Hambleton,
Thomas Barclay,
Steven Bloemen,
Tabetha Boyajian,
Laurance R. Doyle,
B. J. Fulton,
Abe Johannes Hoekstra,
Kian Jek,
Stephen R. Kane,
Veselin Kostov,
David Latham,
Tsevi Mazeh,
Jerome A. Orosz,
Joshua Pepper,
Billy Quarles,
Darin Ragozzine,
Avi Shporer,
John Southworth,
Keivan Stassun
, et al. (25 additional authors not shown)
Abstract:
The primary Kepler Mission provided nearly continuous monitoring of ~200,000 objects with unprecedented photometric precision. We present the final catalog of eclipsing binary systems within the 105 square degree Kepler field of view. This release incorporates the full extent of the data from the primary mission (Q0-Q17 Data Release). As a result, new systems have been added, additional false posi…
▽ More
The primary Kepler Mission provided nearly continuous monitoring of ~200,000 objects with unprecedented photometric precision. We present the final catalog of eclipsing binary systems within the 105 square degree Kepler field of view. This release incorporates the full extent of the data from the primary mission (Q0-Q17 Data Release). As a result, new systems have been added, additional false positives have been removed, ephemerides and principal parameters have been recomputed, classifications have been revised to rely on analytical models, and eclipse timing variations have been computed for each system. We identify several classes of systems including those that exhibit tertiary eclipse events, systems that show clear evidence of additional bodies, heartbeat systems, systems with changing eclipse depths, and systems exhibiting only one eclipse event over the duration of the mission. We have updated the period and galactic latitude distribution diagrams and included a catalog completeness evaluation. The total number of identified eclipsing and ellipsoidal binary systems in the Kepler field of view has increased to 2878, 1.3% of all observed Kepler targets. An online version of this catalog with downloadable content and visualization tools is maintained at http://keplerEBs.villanova.edu.
△ Less
Submitted 29 December, 2015;
originally announced December 2015.
-
Kepler Eclipsing Binary Stars. V. Identification of 31 Eclipsing Binaries in the K2 Engineering Data-set
Authors:
Kyle E. Conroy,
Andrej Prša,
Keivan G. Stassun,
Steven Bloemen,
Mahmoud Parvizi,
Billy Quarles,
Tabetha Boyajian,
Thomas Barclay,
Avi Shporer,
David W. Latham,
Michael Abdul-Masih
Abstract:
Over 2500 eclipsing binaries were identified and characterized from the ultra-precise photometric data provided by the Kepler space telescope. Kepler is now beginning its second mission, K2, which is proving to again provide ultra-precise photometry for a large sample of eclipsing binary stars. In the 1951 light curves covering 12 days in the K2 engineering data-set, we have identified and determi…
▽ More
Over 2500 eclipsing binaries were identified and characterized from the ultra-precise photometric data provided by the Kepler space telescope. Kepler is now beginning its second mission, K2, which is proving to again provide ultra-precise photometry for a large sample of eclipsing binary stars. In the 1951 light curves covering 12 days in the K2 engineering data-set, we have identified and determined the ephemerides for 31 eclipsing binaries that demonstrate the capabilities for eclipsing binary science in the upcoming campaigns in K2. Of those, 20 are new discoveries. We describe both manual and automated approaches to harvesting the complete set of eclipsing binaries in the K2 data, provide identifications and details for the full set of eclipsing binaries present in the engineering data-set, and discuss the prospects for application of eclipsing binary searches in the K2 mission.
△ Less
Submitted 14 July, 2014;
originally announced July 2014.