Researcher profile

H. Sana

H. Sana contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
30works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

30 published item(s)

preprint2022arXiv

A spectroscopic multiplicity survey of Galactic Wolf-Rayet stars: II. The northern WNE sequence

Most massive stars reside in multiple systems that will interact over the course of their lifetime. Classical Wolf-Rayet (WR) stars represent the final end stages of stellar evolution at the upper-mass end. As part of a homogeneous, magnitude-limited ($V\leq12$) spectroscopic survey of northern Galactic WR stars, this paper aims to establish the observed and intrinsic multiplicity properties of the early-type nitrogen-rich WR population (WNE). We obtained high-resolution spectroscopic time series of the complete magnitude-limited sample of 16 WNE stars observable with the 1.2 m Mercator telescope at La Palma, typically providing a time base of about two to eight years. We measured relative radial velocities (RVs) using cross-correlation and used RV variations to flag binary candidates. Adopting a peak-to-peak RV variability threshold of 50 km/s as a criterion found an observed multiplicity fraction of 0.44$\pm$0.12. Using an updated Monte Carlo method with a Bayesian framework, we calculated the three-dimensional likelihood for the intrinsic binary fraction, the maximum period, and the power-law index for the period distribution for the WNE population. We also re-derived multiplicity parameters for the Galactic WC population. We found an intrinsic multiplicity fraction of $0.56\substack{+0.20 \\ -0.15}$ for the parent WNE population. For the Galactic WC population, we re-derive an intrinsic multiplicity fraction of $0.96\substack{+0.04 \\ -0.22}$. The derived multiplicity parameters for the WNE population are quite similar to those derived for main-sequence O binaries but differ from those of the WC population. The significant shift in the WC period distribution towards longer periods is too large to be explained via expansion of the orbit due to stellar winds, and we discuss possible implications of our results.

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

Modeling overcontact binaries, I. The effect of tidal deformation

In the realm of massive stars, strong binary interaction is commonplace. One extreme case are overcontact systems, which is expected to be part of the evolution of all stars evolving towards a merger, and is hypothesized to play a role in the formation of binary black holes. However, important simplifications are made to model the evolution of overcontact binaries, namely that tidal or rotational deformation is frequently ignored. Yet, both observation and theory show that overcontact stars are heavily tidally deformed, leaving a potentially important effect on the outer layers unaccounted for in models. In this work we develop the methodology to represented tidally deformed stars. Using numerical methods, we compute the structure correction factors to the 1D spherical stellar structure equations due to the binary Roche potential, and compare them to existing results and the structure corrections of single rotating stars. We implement the new structure correction factors into the stellar evolution code MESA and explore several case studies. We compare the differences between our simulations when no rotation is included, when we treat rotation using single star corrections (i.e. only accounting for centrifugal deformation) or when we use tidal deformation. We find that ignoring rotation in deformed detached eclipsing binaries can produce a radius discrepancy of up to 5%. The difference between tidal and single star centrifugal distortion models is more benign at 1%, showing that single rotating star models are a suitable approximation of tidally deformed stars in a binary system. In overcontact configurations, we find a similar 5% variation in surface properties as a result of tidal distortion with respect to non-rotating models, showing that it is inappropriate to model binary stars that fill their Roche lobe significantly, as non-rotating.

preprint2022arXiv

On the origin of close massive binaries in the M17 star-forming region

Spectroscopic multiplicity surveys of O stars in young clusters and OB associations have revealed that a large portion ($\sim$ 70%) of these massive stars (M$_{i}$ $\gt$ 15 $M_{\odot}$) belong to close and short-period binaries (physical separation d $\lt$few au). Despite the recent and significant progress, the formation mechanisms leading to such close massive multiple systems remain to be elucidated. As a result, young massive close binaries (or higher-order multiple systems) are unique laboratories to figure out the pairing mechanism of high-mass stars. We present the first VLTI/GRAVITY observations of six young O stars in the M17 star-forming region ($\lesssim$ 1 Myr) and two additional foreground stars. From the interferometric model fitting of visibility amplitudes and closure phases, we search for companions and measure their positions and flux ratios. Combining the resulting magnitude difference with atmosphere models and evolutionary tracks, we further constrain the masses of the individual components. All of the six high-mass stars are in multiple systems, leading to a multiplicity fraction (MF) of 100%, yielding a 68% confidence interval of 94-100%. We detect a total number of 9 companions with separations up to 120 au. Including previously identified spectroscopic companions, the companion fraction of the young O stars in our sample reaches 2.3$\pm$0.6. The derived masses span a wide range from 2.5 to 50 $M_{\odot}$, with a great tendency towards high-mass companions. While based on a modest sample, our results clearly indicate that the origin of the high degree of multiplicity is rooted in their star formation mechanism. No clear evidence for one of the competing concepts of massive star formation (core accretion or competitive accretion) could be found. However, our results are compatible with migration as a scenario for the formation of close massive binaries.

preprint2022arXiv

Properties of the Be-type stars in 30 Doradus

The evolutionary status of Be-type stars remains unclear, with both single-star and binary pathways having been proposed. Here, VFTS spectroscopy of 73 Be-type stars, in the spectral-type range, B0--B3, is analysed to estimate projected rotational velocities, radial velocities and stellar parameters. They are found to be rotating faster than the corresponding VFTS B-type sample but simulations imply that their projected rotational velocities are inconsistent with them all rotating at near critical velocities. The de-convolution of the projected rotational velocities estimates leads to a mean rotational velocity estimate of 320-350 km/s, approximately 100 km/s larger than that for the corresponding B-type sample. There is a dearth of targets with rotational velocities less than 0.4 of the critical velocity, with a broad distribution reaching up to critical rotation. Our best estimate for the mean or median of the rotational velocitiy is 0.68 of the critical velocity. Rapidly-rotating B-type stars are more numerous than their Be-type counterparts, whilst the observed frequency of Be-type stars identified as binary systems is significantly lower than that for normal B-type stars, consistent with their respective radial-velocity dispersions. The semi-amplitudes for the Be-type binaries are also smaller. Similar results are found for a SMC Be-type sample centred on NGC346 with no significant differences being found between the two samples. These results are compared with the predictions of single and binary stellar evolutionary models for Be-type stars. Assuming that a single mechanism dominated the production of classical Be-type stars, our comparison would favour a binary evolutionary history.

preprint2021arXiv

A relation between the radial velocity dispersion of young clusters and their age: Evidence for hardening as the formation scenario of massive close binaries

The majority of massive stars ($>8$ $\rm{M_{\odot}}$) in OB associations are found in close binary systems. Nonetheless, the formation mechanism of these close massive binaries is not understood yet. Using literature data, we measured the radial-velocity dispersion ($σ_\mathrm{RV}$) as a proxy for the close binary fraction in ten OB associations in the Galaxy and the Large Magellanic Cloud, spanning an age range from 1 to 6 Myrs. We find a positive trend of this dispersion with the cluster's age, which is consistent with binary hardening. Assuming a universal binary fraction of $f_\mathrm{bin}$ = 0.7, we converted the $σ_\mathrm{RV}$ behavior to an evolution of the minimum orbital period $P_\mathrm{cutoff}$ from $\sim$9.5 years at 1 Myr to $\sim$1.4 days for the oldest clusters in our sample at $\sim$6 Myr. Our results suggest that binaries are formed at larger separations, and they harden in around 1 to 2 Myrs to produce the period distribution observed in few million year-old OB binaries. Such an inward migration may either be driven by an interaction with a remnant accretion disk or with other young stellar objects present in the system. Our findings constitute the first empirical evidence in favor of migration as a scenario for the formation of massive close binaries.

preprint2021arXiv

Analytic, dust-independent mass-loss rates for red supergiant winds initiated by turbulent pressure

Context. Red supergiants are observed to undergo vigorous mass-loss. However, to date, no theoretical model has succeeded in explaining the origins of these objects' winds. This strongly limits our understanding of red supergiant evolution and Type II-P and II-L supernova progenitor properties. Aims. We examine the role that vigorous atmospheric turbulence may play in initiating and determining the mass-loss rates of red supergiant stars. Methods. We analytically and numerically solve the equations of conservation of mass and momentum, which we later couple to an atmospheric temperature structure, to obtain theoretically motivated mass-loss rates. We then compare these to state-of-the-art empirical mass-loss rate scaling formulae as well as observationally inferred mass-loss rates of red supergiants. Results. We find that the pressure due to the characteristic turbulent velocities inferred for red supergiants is sufficient to explain the mass-loss rates of these objects in the absence of the normally employed opacity from circumstellar dust. Motivated by this initial success, we provide a first theoretical and fully analytic mass-loss rate prescription for red supergiants. We conclude by highlighting some intriguing possible implications of these rates for future studies of stellar evolution, especially in light of the lack of a direct dependence on metallicity.

preprint2021arXiv

Binary-object spectral-synthesis in 3D (BOSS-3D) -- Modelling H-alpha emission in the enigmatic multiple system LB-1

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.

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

Probing the low-mass end of the companion mass function for O-type stars

Almost all massive stars are part of a binary system. Given the wide range of separations at which these companions are found, several observational techniques have been adopted to characterize them, but contrasts greater than 4 in the H-band have never been reached between 0&#34;.1 and 1&#34;. We used VLT/SPHERE to observe simultaneously with the IRDIS and IFS sub-systems 18 O-type stars within 6 kpc and ages between 1-5 Myrs to probe the existence of stellar companions in the angular separation range from 0&#34;.15 to 6&#34; down to very low mass ratios. The IFS YJH- band observations have allowed us to probe the presence of sub-solar companions in a 1&#34;.7x1&#34;.7 field-of-view down to magnitude limits of deltaH=10 at 0&#34;.4. In the wider 12&#34;x12&#34; IRDIS field-of-view, we reached contrasts of deltaK=12 at 1&#34;, enabling us to look for even fainter companions. This paper presents five newly discovered intermediate (<1&#34;) separation companions, three of which are smaller than 0.2M_sun. If confirmed by future analyses of proper motions, these new detections represent the lowest-mass companions ever found around O-type stars. Assuming that all sources detected within 1&#34; are physically bound, the observed fraction of companions for O-type stars between 0&#34;.15 and 0&#34;.9 is 0.39+/-0.15, whereas it increases to 1.6+/-0.3 in the separation range from 0&#34;.9 to 6&#34;. These findings clearly support the notion that massive stars form almost exclusively in multiple systems, and that larger AO-assisted coronagraphic surveys are crucial in placing constraints on the multiplicity properties of massive star companions in regions of the parameter space that have previously gone unexplored, and demonstrate that the companion mass function is populated down to the lowest stellar masses.

preprint2021arXiv

The CubeSpec space mission. I. Asteroseismology of massive stars from time series optical spectroscopy: science requirements and target list prioritisation

Context. There is currently a niche for providing high-cadence, high resolution, time-series optical spectroscopy from space, which can be filled by using a low-cost cubesat mission. The Belgian-led ESA CubeSpec mission is specifically designed to provide space-based, low-cost spectroscopy with specific capabilities that can be optimised for a particular science need. Approved as an ESA in-orbit demonstrator, the CubeSpec satellite&#39;s primary science objective will focus on obtaining high-cadence, high resolution optical spectroscopic data to facilitate asteroseismology of pulsating massive stars. Aims. In this first paper, we aim to search for pulsating massive stars suitable for the CubeSpec mission, specifically $β$ Cep stars, which typically require time series spectroscopy to identify the geometry of their pulsation modes. Methods. Based on the science requirements needed to enable asteroseismology of massive stars with the capabilities of CubeSpec&#39;s spectrograph, we combine a literature study for pulsation with the analysis of recent high-cadence time series TESS photometry to classify the variability for stars brighter than V < 4 mag and between O9 and B3 in spectral type. Results. Among the 90 stars that meet our magnitude and spectral type requirements, we identify 23 promising $β$ Cep stars with high-amplitude (non-)radial pulsation modes with frequencies below 7 d$^{-1}$. Using further constraints on projected rotational velocities, pulsation amplitudes and number of pulsation modes, we devise a prioritised target list for the CubeSpec mission according to its science requirements and the potential of the targets for asteroseismology. The full target catalogue further provide a modern TESS-based review of line profile and photometric variability properties among bright O9-B3 stars.

preprint2021arXiv

The inner circumstellar dust of the red supergiant Antares as seen with VLT/SPHERE/ZIMPOL

The processes by which red supergiants lose mass are not fully understood thus-far and their mass-loss rates lack theoretical constraints. The ambient surroundings of the nearby M0.5 Iab star Antares offers an ideal environment to obtain detailed empirical information on the outflow properties at its onset, and hence indirectly, on the mode(s) of mass loss. We present and analyse optical VLT/SPHERE/ZIMPOL polarimetric imaging with angular resolution down to 23 milli-arcsec, sufficient to spatially resolve both the stellar disk and its direct surroundings. We detect a conspicuous feature in polarised intensity that we identify as a clump containing dust, which we characterise through 3D radiative transfer modelling. The clump is positioned behind the plane of the sky, therefore has been released from the backside of the star, and its inner edge is only 0.3 stellar radii above the surface. The current dust mass in the clump is $1.3^{+0.2}_{-1.0} \times 10^{-8}$ M$_{\odot}$, though its proximity to the star implies that dust nucleation is probably still ongoing. The ejection of clumps of gas and dust makes a non-negligible contribution to the total mass lost from the star which could possibly be linked to localised surface activity such as convective motions or non-radial pulsations.

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.

preprint2021arXiv

Uncovering astrometric black hole binaries with massive main-sequence companions with Gaia

The hunt for compact objects is on. Rarely seen massive binaries with a compact object are a crucial phase in the evolution towards compact object mergers. In Gaia data release 3 (DR3), the first Gaia astrometric orbital solutions for binary sources will become available. We investigate how many black holes (BH) with massive main-sequence dwarf companions (OB+BH binaries) are expected to be detected as binaries in DR3 and at the end of the nominal 5-yr mission (DR4). We estimate the fraction of identifiable OB+BH binaries and discuss the distributions of the masses of both components and the orbital periods. We study the impact of different BH-formation scenarios. Using tailored models for the massive star population, which assume a direct collapse and no kick upon BH formation (the fiducial case), we estimate the fraction of OB+BH systems that Gaia will detect as binaries. A distance distribution according to that of the second Alma Luminous Star catalogue (ALSII) is assumed. We investigate how many of the systems detected as binaries are identifiable as OB+BH binaries, using a method based on astrometric data. In the fiducial case we conservatively estimate that 77% of the OB+BH binaries in ALSII will be detected as binaries in DR3, of which 89% are identifiable as OB+BH binaries. This leads to a total of around 190 OB+BH binaries, a 20-fold increase in the known sample of OB+BH binaries, covering an uncharted parameter space of long-period binaries. The size and properties of the identifiable OB+BH population will contain crucial observational constraints to improve our understanding of BH formation. In DR4, the detected fraction will increase to 85%, of which 82% will be identifiable. Hence, an additional ~5 systems could be identified, which are expected to have either very short or long periods. The fractions become smaller for different BH-formation scenarios. (truncated)

preprint2020arXiv

Investigating the lack of main-sequence companions to massive Be stars

About 20% of all B-type stars are classical Be stars. The Be phenomenon is strongly correlated with rapid rotation, the origin of which remains unclear. It may be rooted in single- or binary-star evolution. In the framework of the binary channel, the initially more massive star transfers mass and angular momentum to the original secondary, which becomes a Be star. The system then evolves into a Be binary with a post-main-sequence companion, which may later be disrupted in a supernova event. Hence, if the binary channel dominates the formation of Be stars, one may expect a strong lack of close Be binaries with main sequence (MS) companions. Through an extensive, star-by-star review of the literature of a magnitude-limited sample of Galactic early-type Be stars, we investigate whether Be binaries with MS companions are known to exist. Our sample is constructed from the BeSS database and cross-matched with all available literature on the individual stars. Out of an initial list of 505 Be stars, we compile a final sample of 287 Galactic Be stars earlier than B1.5 with V<=12 mag. Out of those, 13 objects were reported as Be binaries with known post-MS companions and 11 as binaries with unknown, uncertain or debated companions. We find no confirmed reports of Be binaries with MS companions. For the remaining 263 targets, no significant reports of multiplicity exist in the literature, implying that they are either Be binaries with faint companions, or truly single. The clear lack of reported MS companions to Be stars, which stands in contrast to the high number of detected B+B MS binaries, strongly supports the hypothesis that early-type Be stars are binary interaction products that spun up after mass and angular momentum transfer from a companion star. Taken at face value, our results may suggest that a large majority of the early-type Be stars have formed through binary mass-transfer.

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

Multiplicity of the red supergiant population in the young massive cluster NGC 330

The multiplicity properties of massive stars are one of the important outstanding issues in stellar evolution. Quantifying the binary statistics of all evolutionary phases is essential to paint a complete picture of how and when massive stars interact with their companions, and to determine the consequences of these interactions. We investigate the multiplicity of an almost complete census of red supergiant stars (RSGs) in NGC 330, a young massive cluster in the SMC. Using a combination of multi-epoch HARPS and MUSE spectroscopy, we estimate radial velocities and assess the kinematic and multiplicity properties of 15 RSGs in NGC 330. Radial velocities are estimated to better than +/-100 m/s for the HARPS data. The line-of-sight velocity dispersion for the cluster is estimated as 3.20 +0.69-0.52 km/s. When virial equilibrium is assumed, the dynamical mass of the cluster is log (M{dyn} /M{sun}) = 5.20+/-0.17, in good agreement with previous upper limits. We detect significant radial velocity variability in our multi-epoch observations and distinguish between variations caused by atmospheric activity and those caused by binarity. The binary fraction of NGC 330 RSGs is estimated by comparisons with simulated observations of systems with a range of input binary fractions. In this way, we account for observational biases and estimate the intrinsic binary fraction for RSGs in NGC 330 as f{RSG} = 0.3+/-0.1 for orbital periods in the range 2.3< log P [days] <4.3, with q>0.1. Using the distribution of the luminosities of the RSG population, we estimate the age of NGC 330 to be 45+/-5 Myr and estimate a red straggler fraction of 50%. We estimate the binary fraction of RSGs in NGC 330 and conclude that it appears to be lower than that of main-sequence massive stars, which is expected because interactions between an RSG and a companion are assumed to effectively strip the RSG envelope.

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 &#39;hidden&#39; 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 &#34;hidden&#34; 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.

preprint2020arXiv

The high-energy emission from HD~93129A near periastron

We conducted an observational campaign towards one of the most massive and luminous colliding wind binaries in the Galaxy, HD~93129A, close to its periastron passage in 2018. During this time the source was predicted to be in its maximum of high-energy emission. Here we present our data analysis from the X-ray satellites \textit{Chandra} and \textit{NuSTAR} and the $γ$-ray satellite \textit{AGILE}. High-energy emission coincident with HD~93129A was detected in the X-ray band up to $\sim$18~keV, whereas in the $γ$-ray band only upper limits were obtained. We interpret the derived fluxes using a non-thermal radiative model for the wind-collision region. We establish a conservative upper limit for the fraction of the wind kinetic power that is converted into relativistic electron acceleration, $f_\mathrm{NT,e} < 0.02$. In addition, we set a lower limit for the magnetic field in the wind-collision region as $B_\mathrm{WCR} > 0.3$~G. We also argue a putative interpretation of the emission from which we estimate $f_\mathrm{NT,e} \approx 0.006$ and $B_\mathrm{WCR} \approx 0.5$~G. We conclude that multi-wavelength, dedicated observing campaigns during carefully selected epochs are a powerful tool for characterising the relativistic particle content and magnetic field intensity in colliding wind binaries.

preprint2020arXiv

The VLT-FLAMES Tarantula Survey XXXII. Low-luminosity late O-type stars -- classification, main physical parameters, and silicon abundances

Analysis of late O-type stars observed in the Large Magellanic Cloud (LMC) by the VLT-FLAMES Tarantula Survey (VFTS) revealed a discrepancy between the physical properties estimated from model-atmosphere analysis and those expected from their morphological classifications. Here we revisit the analysis of 32 of these puzzling objects using new hydrogen-helium-silicon FASTWIND models and a different fitting approach to re-evaluate their physical properties. Our new analysis confirms that these stars indeed have properties that are typical of late O-type dwarfs. We also present the first estimates of silicon abundances for O-type stars in the 30 Dor clusters NGC 2060 and NGC 2070, with a weighted mean abundance for our sample of 7.05 +/- 0.03. Our values are about 0.20 dex lower than those previously derived for B-type stars in the LMC clusters N 11 and NGC 2004 using TLUSTY models. Various possibilities (e.g. differences in the analysis methods, effects of microturbulence, and real differences between stars in different clusters) were considered to account for these results. We also used our grid of FASTWIND models to reassess the impact of using the Galactic classification criteria for late O-type stars in the LMC by scrutinising their sensitivity to different stellar properties. At the cool edge of the O star regime the HeII 4686/HeI 4713 ratio used to assign luminosity class for Galactic stars can mimic giants or bright giants in the LMC, even for objects with high gravities (log_g > 4.0 dex). We argue that this line ratio is not a reliable luminosity diagnostic for late O-type stars in the LMC, and that the SiIV 4989/HeI4026 ratio is more robust for these types.

preprint2020arXiv

The Wolf-Rayet binaries of the nitrogen sequence in the Large Magellanic Cloud: spectroscopy, orbital analysis, formation, and evolution

Massive Wolf-Rayet (WR) stars dominate the radiative and mechanical energy budget of galaxies and probe a critical phase in the evolution of massive stars prior to core-collapse. It is not known whether core He-burning WR stars (classical WR, cWR) form predominantly through wind-stripping (w-WR) or binary stripping (b-WR). With spectroscopy of WR binaries so-far largely avoided due to its complexity, our study focuses on the 44 WR binaries / binary candidates of the Large Magellanic Cloud (LMC, metallicity Z~0.5 Zsun), identified on the basis of radial velocity variations, composite spectra, or high X-ray luminosities. Relying on a diverse spectroscopic database, we aim to derive the physical and orbital parameters of our targets, confronting evolution models of evolved massive stars at sub-solar metallicity, and constraining the impact of binary interaction in forming them. Spectroscopy is performed using the Potsdam Wolf-Rayet (PoWR) code and cross-correlation techniques. Disentanglement is performed using the code Spectangular or the shift-and-add algorithm. Evolutionary status is interpreted using the Binary Population and Spectral Synthesis (BPASS) code, exploring binary interaction and chemically-homogeneous evolution. No obvious dichotomy in the locations of apparently-single and binary WN stars on the Hertzsprung-Russell diagram is apparent. According to commonly used stellar evolution models (BPASS, Geneva), most apparently-single WN stars could not have formed as single stars, implying that they were stripped by an undetected companion. Otherwise, it must follow that pre-WR mass-loss/mixing (e.g., during the red supergiant phase) are strongly underestimated in standard stellar evolution models.

preprint2020arXiv

VLT/X-shooter spectroscopy of massive young stellar objects in the 30 Doradus region of the Large Magellanic Cloud

The process of massive star ($M\geq8~M_\odot$) formation is still poorly understood. Observations of massive young stellar objects (MYSOs) are challenging due to their rarity, short formation timescale, large distances, and high circumstellar extinction. Here, we present the results of a spectroscopic analysis of a population of MYSOs in the Large Magellanic Cloud (LMC). We took advantage of the spectral resolution and wavelength coverage of X-shooter (300-2500 nm), mounted on the European Southern Observatory Very Large Telescope, to detect characteristic spectral features in a dozen MYSO candidates near 30 Doradus, the largest starburst region in the Local Group hosting the most massive stars known. The X-shooter spectra are strongly contaminated by nebular emission. We used a scaling method to subtract the nebular contamination from our objects. We detect H$α,β$, [O I] 630.0 nm, Ca II infrared triplet, [Fe II] 1643.5 nm, fluorescent Fe II 1687.8 nm, H$_2$ 2121.8 nm, Br$γ$, and CO bandhead emission in the spectra of multiple candidates. This leads to the spectroscopic confirmation of 10 candidates as bona fide MYSOs. We compare our observations with photometric observations from the literature and find all MYSOs to have a strong near-infrared excess. We compute lower limits to the brightness and luminosity of the MYSO candidates, confirming the near-infrared excess and the massive nature of the objects. No clear correlation is seen between the Br$γ$ luminosity and metallicity. Combining our sample with other LMC samples results in a combined detection rate of disk features such as fluorescent Fe II and CO bandheads which is consistent with the Galactic rate (40\%). Most of our MYSOs show outflow features.

preprint2020arXiv

Why binary interaction does not necessarily dominate the formation of Wolf-Rayet stars at low metallicity

Classical Wolf-Rayet (WR) stars are massive, hydrogen depleted, post main-sequence stars that exhibit emission-line dominated spectra. For a given metallicity Z, stars exceeding a certain initial mass M_single(Z) can reach the WR phase through intrinsic mass-loss (single-star channel). Stars of lower masses can reach the WR phase via binary mass transfer (binary channel). It is commonly assumed that the binary channel dominates the formation of WR stars in environments with low Z such as the SMC and LMC. However, their reported WR binary fractions of 30-40% are comparable to that of the Galaxy. Here, we explain this apparent contradiction by considering the minimum initial mass M_spec(Z) needed for the stripped product to appear as a WR star. We calibrate M_spec(Z) using the lowest-luminosity WR stars in the Clouds and the Galaxy. A range of M_single(Z) values are explored using various evolution codes. We estimate the additional contribution of the binary channel by considering the interval [M_spec(Z), M_single(Z)], which characterises the initial-mass range in which binaries can form additional WR stars. Results: The WR-phenomenon ceases below luminosities of logL = 4.9, 5.25, and 5.6 [Lsun] in the Galaxy, LMC, and SMC, which translates to He-star masses of 7.5, 11, 17 Msun and initial masses of M_spec = 18, 23, 37 Msun. Stripped stars with lower initial masses in the respective galaxies would tend to not appear as WR stars. M_single lies in the range 20-30, 30-60, and > 40 Msun for the Galaxy, LMC, and SMC. We find that that the additional contribution of the binary channel is a non-trivial function of Z that cannot be conclusively claimed to be monotonically increasing with decreasing Z. Hence, one should not a-priori expect that binary interactions become increasingly important in forming WR stars at low Z, or that the WR binary fraction grows with decreasing Z.

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.

preprint2019arXiv

The young stellar content of the giant HII regions M8, G333.6-0.2, and NGC6357 with VLT/KMOS

Context: The identification and characterisation of populations of young massive stars in (giant) HII regions provides important constraints on i) the formation process of massive stars and their early feedback on the environment, and ii) the initial conditions for population synthesis models predicting the evolution of ensembles of stars. Aims: We identify and characterise the stellar populations of the following young giant HII regions: M8, G333.6-0.2, and NGC6357. Methods: We acquired H- and K-band spectra of around 200 stars using The K-band KMOS on the ESO Very Large Telescope. The targets for M8 and NGC6357 were selected from the MYStIX project, which combines X-ray observations with near-infrared and mid-infrared data. For G333.6-0.2, the sample selection is based on the near-infrared colours combined with X-ray data. We introduce an automatic spectral classification method in order to obtain temperatures and luminosities for the observed stars. We analyse the stellar populations using their photometric, astrometric, and spectroscopic properties and compared the position of the stars in the Hertzprung-Russell diagram with stellar evolution models to constrain their ages and mass ranges. Results: We confirm the presence of candidate ionising sources in the three regions and report new ones, including the first spectroscopically identified O stars in G333.6-0.2. In M8 and NGC6357, two populations are identified: i) OB main-sequence stars ($M > 5~\rm{M_{\odot}}$) and ii) pre-main sequence stars ($M\approx0.5-5~\rm{M_{\odot}}$). The ages of the clusters are $\sim$1-3~Myr, $< 3$~Myr, and $\sim$0.5-3~Myr for M8, G333.6-0.2, and NGC6357, respectively. We show that MYStIX selected targets have $>$ 90\% probability of being members of the HII region, whereas a selection based on near infrared (NIR) colours leads to a membership probability of only $\sim$70\%.

preprint2018arXiv

Response to comment on &#34;An excess of massive stars in the local 30 Doradus starburst&#34;

Farr and Mandel reanalyse our data, finding initial-mass-function slopes for high mass stars in 30 Doradus that agree with our results. However, their reanalysis appears to underpredict the observed number of massive stars. Their technique results in more precise slopes than in our work, strengthening our conclusion that there is an excess of massive stars above $30\,\mathrm{M}_\odot$ in 30 Doradus.