Showing 1–50 of 196 results for author: Rauw, G

A family of binaries with an extreme mass ratio

Authors: Yael Naze, Gregor Rauw, Piotr A. Kolaczek-Szymanski, Nikolay Britavskiy, Jonathan Labadie-Bartz

Abstract: Multiplicity is ubiquitous among massive stars. While the stellar components usually display similar masses, some binaries with extremely low mass ratios were also observed. Some of them are primordial, while others arise from binary interactions. The identification of systems with extreme mass ratios brings valuable information, notably on the origin of fast rotation in massive stars. We identify… ▽ More Multiplicity is ubiquitous among massive stars. While the stellar components usually display similar masses, some binaries with extremely low mass ratios were also observed. Some of them are primordial, while others arise from binary interactions. The identification of systems with extreme mass ratios brings valuable information, notably on the origin of fast rotation in massive stars. We identify new short-period systems with extreme mass ratios through the detection of eclipses and reflection effects. The physical properties of a dozen newly identified cases were precisely evaluated through high-quality photometry and spectroscopy. In addition to characterizing these binaries, we found a clear signature of apsidal motion in one system, and three other systems display long-term shifts in eclipse times. All systems we reported here are composed of a massive star and a cool low-mass companion. They are therefore primordial cases. This doubles the known number of these systems in the Galaxy. In this context, it is important to note that most massive stars in these systems, as well as in previous systems reported in the literature, rotate fast (supersynchronous compared to the orbital motion). The high incidence of fast rotation in these nascent binaries provides strong constraints for star formation models. △ Less

Submitted 27 October, 2025; v1 submitted 17 October, 2025; originally announced October 2025.

Comments: accepted by A&A

Fluorescent Fe K line emission of gamma Cas stars II. Predictions for magnetic interactions and white-dwarf accretion scenarios

Authors: Gregor Rauw

Abstract: About 12 percent of the early-type Be stars, so-called gamma Cas stars, exhibit an unusually hard and bright thermal X-ray emission that could result from accretion onto a white dwarf companion or from magnetic interactions between the Be star and its decretion disk. Exploring the full power of high-resolution X-ray spectroscopy of gamma Cas stars requires comparison of observations of the fluores… ▽ More About 12 percent of the early-type Be stars, so-called gamma Cas stars, exhibit an unusually hard and bright thermal X-ray emission that could result from accretion onto a white dwarf companion or from magnetic interactions between the Be star and its decretion disk. Exploring the full power of high-resolution X-ray spectroscopy of gamma Cas stars requires comparison of observations of the fluorescent Fe Kalpha emission lines near 6.4 keV with synthetic profiles of this line complex computed in the framework of the magnetic interaction and the accreting WD scenarios. For the latter, we further distinguish between accretion onto a non-magnetic and a magnetic WD. Our models account for different reservoirs of reprocessing material: the Be circumstellar decretion disk, the Be photosphere, an accretion disk around the WD companion, a magnetically channelled accretion flow and the WD photosphere. We find considerably different line properties for the different scenarios. For a non-magnetic accreting WD, the global Fe Kalpha complex is extremely broad, reaching a full width of 140 eV, whilst it is ~ 40 eV for the magnetic star-disk interaction and the magnetic accreting WD cases. In the magnetic star-disk interaction, the line centroid follows the orbital motion of the Be star, whereas it moves along with the WD in the case of an accreting WD. For gamma Cas, given the 15 times larger amplitude of the WD orbital motion, the shift in position for an accreting WD should be easily detectable with high-resolution spectrographs such as Resolve on XRISM, but remains essentially undetectable for the magnetic star-disk interaction. Upcoming high-resolution spectroscopy of the fluorescent Fe Kalpha emission lines in the X-ray spectra of gamma Cas stars will thus allow to distinguish between the competing scenarios. △ Less

Submitted 18 October, 2025; v1 submitted 3 October, 2025; originally announced October 2025.

Comments: Accepted for publication in A&A

Journal ref: A&A 703, A188 (2025)

Apsidal motion and proximity effects in the massive binary BD+60 497

Authors: Gregor Rauw, Piotr A. Kolaczek-Szymanski, Yael Naze, Lucas Nys

Abstract: The eccentric short-period O-star binary BD+60 497 is an interesting laboratory in which to study tidal interactions in massive binary systems, notably via the detection and characterisation of apsidal motion. The rate of apsidal motion in such systems can help constrain their age and provide insight into the degree of mass concentration in the interior of massive stars. We used spectroscopic data… ▽ More The eccentric short-period O-star binary BD+60 497 is an interesting laboratory in which to study tidal interactions in massive binary systems, notably via the detection and characterisation of apsidal motion. The rate of apsidal motion in such systems can help constrain their age and provide insight into the degree of mass concentration in the interior of massive stars. We used spectroscopic data collected over two decades to reconstruct the individual spectra of the stars and to establish their epoch-dependent radial velocities. An orbital solution, explicitly accounting for apsidal motion was adjusted to the data. Space-borne photometric time series were analysed with Fourier methods and with binary models. We derived a rate of apsidal motion of $6.15^{+1.05}_{-1.65}$ degree/yr, which suggests an age of $4.13^{+0.42}_{-1.37}$ Myr. The disentangled spectra unveiled a curious change in the spectral properties of the secondary star between the epochs 2002-2003 and 2018-2022 with the secondary spectrum appearing to be of an earlier spectral type over recent years. Photometric data show variability at the 6 mmag level on the period of the binary system, which is hard to explain in terms of proximity effects. Whilst the rate of apsidal motion agrees well with theoretical expectations, the changes in the reconstructed secondary spectrum hint at a highly non-uniform surface temperature distribution for this star. Different effects are discussed that could contribute to the photometric variations. The current most-likely explanation is a mix of proximity effects and tidally excited oscillations △ Less

Submitted 2 October, 2025; v1 submitted 16 September, 2025; originally announced September 2025.

Comments: Accepted for publication in A&A

Journal ref: A&A 702, A274 (2025)

The multi-faceted variability of HD 192639: Stochastic behaviour, regularities, and an exceptional mass-ejection event

Authors: Gregor Rauw, Yaël Nazé, Charles-Antoine Gilon

Abstract: Spectroscopic and photometric variability is widespread among O-type supergiants. It is linked to various phenomena affecting the star and its circumstellar environment, thereby providing direct information concerning them. To investigate such connections, we decided to revisit the prototypical O7.5 Iabf supergiant HD 192639. High-cadence spectroscopic monitoring was performed simultaneously with… ▽ More Spectroscopic and photometric variability is widespread among O-type supergiants. It is linked to various phenomena affecting the star and its circumstellar environment, thereby providing direct information concerning them. To investigate such connections, we decided to revisit the prototypical O7.5 Iabf supergiant HD 192639. High-cadence spectroscopic monitoring was performed simultaneously with intensive space-borne photometric observations. The data were analysed with several methods to characterise the variability. Besides the usual stochastic, low-frequency photometric variability, our observations reveal the presence of recurrent - but transient - modulations on a timescale of about five days. The same signal is present in the spectroscopic data and was already seen two decades ago. This stability suggests that this timescale corresponds to the stellar rotation. Furthermore, our observations unveil, for the first time, an unusually strong dimming event in the light curve associated with absorption and emission changes in H I and He I lines. This unprecedented trough corresponds to an episodic ejection of a rather large amount of mass (its column density being comparable to that of the steady wind). While rare, such an event could hint at an overlooked aspect of mass loss in massive stars. △ Less

Submitted 3 August, 2025; originally announced August 2025.

Comments: Accepted for publucation in A&A

Another one (BH+OB pair) bites the dust

Authors: Yael Naze, Gregor Rauw

Abstract: Most (or possibly all) massive stars reside in multiple systems. From stellar evolution models, numerous systems with an OB star coupled to a black hole would be expected to exist. There have been several claimed detections of such pairs in recent years and this is notably the case of HD96670. Using high-quality photometry and spectroscopy in the optical range, we revisited the HD96670 system. We… ▽ More Most (or possibly all) massive stars reside in multiple systems. From stellar evolution models, numerous systems with an OB star coupled to a black hole would be expected to exist. There have been several claimed detections of such pairs in recent years and this is notably the case of HD96670. Using high-quality photometry and spectroscopy in the optical range, we revisited the HD96670 system. We also examined complementary X-ray observations to provide a broader view of the system properties. The TESS light curves of HD96670 clearly show eclipses, ruling out the black hole companion scenario. This does not mean that the system is not of interest. Indeed, the combined analysis of photometric and spectroscopic data indicates that the system most likely consists of a O8.5 giant star paired with a stripped-star companion with a mass of ~4.5Msol, a radius of ~1Rsol, and a surface temperature of ~50kK. While several B+sdOB systems have been reported in the literature, this would be the first case of a Galactic system composed of an O star and a faint stripped star. In addition, the system appears brighter and harder than normal OB stars in the X-ray range, albeit less so than for X-ray binaries. The high-energy observations provide hints of phase-locked variations, as typically seen in colliding wind systems. As a post-interaction system, HD96670 actually represents a key case for probing binary evolution, even if it is not ultimately found to host a black hole. △ Less

Submitted 11 March, 2025; originally announced March 2025.

Comments: accepted for publication in A&A

The X-ray Integral Field Unit at the end of the Athena reformulation phase

Authors: Philippe Peille, Didier Barret, Edoardo Cucchetti, Vincent Albouys, Luigi Piro, Aurora Simionescu, Massimo Cappi, Elise Bellouard, Céline Cénac-Morthé, Christophe Daniel, Alice Pradines, Alexis Finoguenov, Richard Kelley, J. Miguel Mas-Hesse, Stéphane Paltani, Gregor Rauw, Agata Rozanska, Jiri Svoboda, Joern Wilms, Marc Audard, Enrico Bozzo, Elisa Costantini, Mauro Dadina, Thomas Dauser, Anne Decourchelle , et al. (257 additional authors not shown)

Abstract: The Athena mission entered a redefinition phase in July 2022, driven by the imperative to reduce the mission cost at completion for the European Space Agency below an acceptable target, while maintaining the flagship nature of its science return. This notably called for a complete redesign of the X-ray Integral Field Unit (X-IFU) cryogenic architecture towards a simpler active cooling chain. Passi… ▽ More The Athena mission entered a redefinition phase in July 2022, driven by the imperative to reduce the mission cost at completion for the European Space Agency below an acceptable target, while maintaining the flagship nature of its science return. This notably called for a complete redesign of the X-ray Integral Field Unit (X-IFU) cryogenic architecture towards a simpler active cooling chain. Passive cooling via successive radiative panels at spacecraft level is now used to provide a 50 K thermal environment to an X-IFU owned cryostat. 4.5 K cooling is achieved via a single remote active cryocooler unit, while a multi-stage Adiabatic Demagnetization Refrigerator ensures heat lift down to the 50 mK required by the detectors. Amidst these changes, the core concept of the readout chain remains robust, employing Transition Edge Sensor microcalorimeters and a SQUID-based Time-Division Multiplexing scheme. Noteworthy is the introduction of a slower pixel. This enables an increase in the multiplexing factor (from 34 to 48) without compromising the instrument energy resolution, hence keeping significant system margins to the new 4 eV resolution requirement. This allows reducing the number of channels by more than a factor two, and thus the resource demands on the system, while keeping a 4' field of view (compared to 5' before). In this article, we will give an overview of this new architecture, before detailing its anticipated performances. Finally, we will present the new X-IFU schedule, with its short term focus on demonstration activities towards a mission adoption in early 2027. △ Less

Submitted 15 February, 2025; originally announced February 2025.

Comments: 44 pages, 14 figures, accepted for publication in Experimental Astronomy

Caught in flagrante delicto: evidence for past mass transfer in massive binaries?

Authors: Gregor Rauw

Abstract: Many massive binary systems undergo mass and angular momentum transfer over the course of their evolution. This kind of interaction is expected to deeply affect the properties of the mass donor and mass gainer and to leave various observational signatures. The most common smoking guns of a past mass transfer episode are notably rapid rotation of the mass gainer and altered surface chemical abundan… ▽ More Many massive binary systems undergo mass and angular momentum transfer over the course of their evolution. This kind of interaction is expected to deeply affect the properties of the mass donor and mass gainer and to leave various observational signatures. The most common smoking guns of a past mass transfer episode are notably rapid rotation of the mass gainer and altered surface chemical abundances of the stripped mass donor star. Quantitative observational studies of evolved massive binaries are crucial to gain insight into poorly constrained parameters of binary evolution models such as the fraction of mass lost by the mass donor that is actually accreted by the mass gainer. Yet, drawing conclusions about a past mass transfer episode requires a detailed analysis of all aspects of a binary system which sometimes leads to unexpected results. In this contribution, we review the existing observational evidence for past mass exchange events in massive main-sequence and post main-sequence binaries. △ Less

Submitted 18 October, 2024; originally announced November 2024.

Comments: Accepted for publication in the Bulletin de la Société Royale des Sciences de Liège, Proceedings of the 41st Liège Astrophysical Colloquium, held in Liège 15-19 July 2024, edited by Y. Nazé, & E. Bozzo

X-raying the wind-wind collisions in HD 168112 and HD 167971

Authors: Gregor Rauw, Ronny Blomme, Yaël Nazé, Delia Volpi

Abstract: The O-type long-period binary HD 168112 and triple HD 167971 star systems have been known for several decades for their non-thermal synchrotron radio emission. This emission arises from relativistic electrons accelerated in the hydrodynamic shocks of the wind collisions in these systems. Such wind collisions are expected to produce a strong X-ray emission that varies as a function of orbital phase… ▽ More The O-type long-period binary HD 168112 and triple HD 167971 star systems have been known for several decades for their non-thermal synchrotron radio emission. This emission arises from relativistic electrons accelerated in the hydrodynamic shocks of the wind collisions in these systems. Such wind collisions are expected to produce a strong X-ray emission that varies as a function of orbital phase. In wide eccentric binaries, such as our targets, the X-ray emission arises from an adiabatic plasma and its intensity should scale as the inverse of the orbital separation. We present a set of XMM-Newton observations of these systems which help us gain insight into the properties of their wind interactions. △ Less

Submitted 18 October, 2024; originally announced October 2024.

Comments: in French language. Accepted for publication in the Bulletin de la Société Royale des Sciences de Liège, Proceedings of the 41st Liège Astrophysical Colloquium, held in Liège 15-19 July 2024, edited by Y. Nazé, & E. Bozzo

How does the solar chromospheric activity look like under different inclination angles?

Authors: G. Vanden Broeck, S. Bechet, G. Rauw, F. Clette

Abstract: Chromospheric plages are distributed between mid-latitude and the Equator and never close to the Poles. Therefore, we suspect that the inclination angle of the solar rotation axis has an impact on the observable chromospheric emission. We reproduce the solar images from any inclination in order to study the effect of the inclination axis on the solar variability by using direct observations of the… ▽ More Chromospheric plages are distributed between mid-latitude and the Equator and never close to the Poles. Therefore, we suspect that the inclination angle of the solar rotation axis has an impact on the observable chromospheric emission. We reproduce the solar images from any inclination in order to study the effect of the inclination axis on the solar variability by using direct observations of the Sun in the Ca II K line. More than 2700 days of observations since the beginning of the Ca II K observations with USET, in July 2012, were used in our analysis. For each observation day, we produce synoptic maps to map the entire solar surface during a full solar rotation. Then by choosing a given inclination, we generate solar-disk views, representing the segmented brightest structures of the chromosphere (plages and enhanced network), as seen under this inclination. The area fraction are extracted from the masks for each inclination and we compare the evolution of those time series to quantify the impact of the inclination angle. We find a variation of the area fraction between an Equator-on view and a Pole-on view. Our results show an important impact of the viewing angle on the detection of modulation due to the solar rotation. With the dense temporal sampling of USET data, the solar rotation is detectable up to an inclination of about $|i| = 70^{\circ}$ and the solar-cycle modulation is clearly detected for all inclinations, though with a reduced amplitude in polar views. When applying a sparse temporal sampling typical for time series of solar-like stars, the rotational modulation is no longer detected, whatever the inclination. On the other hand, we find that the activity-cycle modulation remains detectable as long as the sampling contains at least 20 observations per year and the cycle amplitude reaches at least 30\% of the solar-cycle amplitude. △ Less

Submitted 15 October, 2024; originally announced October 2024.

Comments: Submitted to A&A

Relationship between TIGRE solar S-index and USET Ca II K full disk images

Authors: G. Vanden Broeck, S. Bechet, F. Clette, G. Rauw, K. -P. Schröder, M. Mittag

Abstract: Full disk observations of the solar chromosphere in the Ca II K line represent a valuable dataset for studies of solar magnetic activity. The S-index is widely used to investigate the magnetic activity of stars, however, its connection to the coverage of stellar magnetic structure is still poorly understood. We use the archives of full disk Ca II K images taken by the Royal Observatory of Belgium… ▽ More Full disk observations of the solar chromosphere in the Ca II K line represent a valuable dataset for studies of solar magnetic activity. The S-index is widely used to investigate the magnetic activity of stars, however, its connection to the coverage of stellar magnetic structure is still poorly understood. We use the archives of full disk Ca II K images taken by the Royal Observatory of Belgium with the USET to derive the area fraction of the brightest chromospheric structures over the last decade. These data allowed us to study the end of the solar cycle 24 and the beginning of solar cycle 25. The brightest regions of the solar surface were then segmented using an algorithm based on an intensity threshold. We computed the area fraction over the solar disk and compared it with the S-index from TIGRE. For the detection of periodic modulations, we applied a discrete Fourier power spectrum method to both datasets. A tight linear relationship was found between the USET area fraction and the TIGRE S-index, with an improved correlation obtained in the low-activity regime by considering the enhanced network. In both time series, we detected the modulation caused by the rotation of bright structures on the solar disk. However, this detection is constrained in the case of TIGRE due to its observation strategy. We studied the correlation between the disk coverage with chromospheric structures and the variability of the S-index on an overlapping period of ten years. We concluded that the disk coverage index is a good proxy for the S-index and will be useful in future studies of the magnetic activity of solar-type stars. The USET area fraction dataset is most appropriate for evaluating the solar rotation period and will be used in future works to analyze the impact of the inclination of the stellar rotation axis on the detectability of such periodic modulations in solar-type stars. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: Accepted for publication in Astronomy & Astrophysics

X-raying the zeta Tau binary system

Authors: Yael Naze, Christian Motch, G. Rauw, Myron A. Smith, Jan Robrade

Abstract: The Be star zeta Tau was recently reported to be a gamma Cas analog; that is, it displays an atypical (bright and hard) X-ray emission. The origin of these X-rays remains debated.The first X-ray observations indicated a very large absorption of the hot plasma component (N_H~ 10^{23}/cm^2). This is most probably related to the edge-on configuration of the zeta Tau disk. If the X-ray emission arises… ▽ More The Be star zeta Tau was recently reported to be a gamma Cas analog; that is, it displays an atypical (bright and hard) X-ray emission. The origin of these X-rays remains debated.The first X-ray observations indicated a very large absorption of the hot plasma component (N_H~ 10^{23}/cm^2). This is most probably related to the edge-on configuration of the zeta Tau disk. If the X-ray emission arises close to the companion, an orbital modulation of the absorption could be detected as the disk comes in and out of the line of sight. New XMM-Newton data were obtained to characterize the high-energy properties of zeta Tau in more detail. They are complemented by previous Chandra and SRG/eROSITA observations as well as by optical spectroscopy and TESS photometry. The high-quality XMM-Newton data reveal the presence of a faint soft X-ray emission, which appears in line with that recorded for non-gamma Cas Be stars. In addition, zeta Tau exhibits significant short-term variability at all energies, with larger amplitudes at lower frequencies (``red noise''), as is found in X-ray data of other gamma Cas stars. Transient variability (softness dip, low-frequency signal) may also be detected at some epochs. In addition, between X-ray exposures, large variations in the spectra are detected in the 1.5-4.keV energy band. They are due to large changes in absorption toward the hottest (9keV) plasma. These changes are not correlated with either the orbital phase or the depth of the shell absorption of the Halpha line. These observed properties are examined in the light of proposed gamma Cas models. △ Less

Submitted 11 July, 2024; v1 submitted 21 June, 2024; originally announced June 2024.

Comments: accepted for publication by A&A

Ups and downs in the X-ray emission of the colliding wind binaries HD 168112 and HD 167971

Authors: Gregor Rauw, Ronny Blomme, Yael Naze, Delia Volpi, Sorenza Fernandez-Vera

Abstract: The long-period O-star binary system HD 168112 and the triple O-star system HD 167971 are well-known sources of non-thermal radio emission that arises from a colliding wind interaction. The wind-wind collisions in these systems should result in phase-dependent X-ray emissions. The presence of a population of relativistic electrons in the wind interaction zone could affect the properties of the X-r… ▽ More The long-period O-star binary system HD 168112 and the triple O-star system HD 167971 are well-known sources of non-thermal radio emission that arises from a colliding wind interaction. The wind-wind collisions in these systems should result in phase-dependent X-ray emissions. The presence of a population of relativistic electrons in the wind interaction zone could affect the properties of the X-ray emission and make it deviate from the behaviour expected for adiabatic shocks. We investigate the X-ray emission of these systems with the goals of quantifying the fraction of the X-ray flux arising from wind interactions and determining whether these emissions follow the predictions for adiabatic wind-wind collisions. Six X-ray observations were collected with XMM-Newton. Three observations were scheduled around the most recent periastron passage of HD 168112. Spectra and light curves were analysed and compared with simple predictions of model calculations for X-ray emission from colliding wind systems. The X-ray emission of HD 168112 varies as the inverse of the orbital separation, as expected for an adiabatic wind interaction zone. The relative contribution of intrinsic X-ray emission from wind-embedded shocks varies between 38% at periastron to 81% at apastron. The wind-wind collision zone remains adiabatic even around periastron passage. The X-ray emission of HD 167971 displays variations on the orbital timescale of the inner eclipsing binary. The existing data of this system do not allow us to probe variations on the timescale of the outer orbit. Shock modification due to the action of relativistic electrons does not seem to be efficiently operating in the HD 168112 system. In the existing observations, a significant part of the emission of HD 167971 must arise in the inner eclipsing binary. The origin of this emission is as yet unclear. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: Accepted for publication in Astronomy & Astrophysics

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.

Comments: 13 pages, 7 figures, accepted for publication in Astronomy & Astrophysics. Corrected values of Omega on p. 9

Phase-resolved XMM-Newton observations of the massive post-RLOF system HD 149404

Authors: G. Rauw, S. Lizin, S. Rosu, E. Mossoux

Abstract: We investigated the X-ray emission of HD 149404, a 9.81-day period O-star binary in a post-Roche lobe overflow evolutionary stage. X-ray emission of O-star binaries consists of the intrinsic emission of the individual O stars and a putative additional component arising from the wind-wind interaction. Phase-locked variations in the X-ray spectra can be used to probe the properties of the stellar wi… ▽ More We investigated the X-ray emission of HD 149404, a 9.81-day period O-star binary in a post-Roche lobe overflow evolutionary stage. X-ray emission of O-star binaries consists of the intrinsic emission of the individual O stars and a putative additional component arising from the wind-wind interaction. Phase-locked variations in the X-ray spectra can be used to probe the properties of the stellar winds of such systems. XMM-Newton observations of HD 149404 collected at two conjunction phases and a quadrature phase were analysed. X-ray spectra were extracted and flux variations as a function of orbital phase were inferred. The flux ratios were analysed with models considering various origins for the X-ray emission. The highest and lowest X-ray fluxes are observed at conjunction phases respectively with the primary and secondary star in front. The flux variations are nearly grey with only marginal energy dependence. None of the models accounting for photoelectric absorption by homogeneous stellar winds perfectly reproduces the observed variations. Whilst the overall X-ray luminosity is consistent with a pure intrinsic emission, the best formal agreement with the observed variations is obtained with a model assuming pure wind-wind collision X-ray emission. The lack of significant energy-dependence of the opacity most likely hints at the presence of optically thick clumps in the winds of HD 149404. △ Less

Submitted 22 March, 2024; originally announced March 2024.

Comments: Accepted for publication in Astronomy & Astrophysics

Tracing the evolution of short-period binaries with super-synchronous fast rotators

Authors: N. Britavskiy, M. Renzo, Y. Nazé, G. Rauw, P. Vynatheya

Abstract: The initial distribution of rotational velocities of stars is still poorly known, and how the stellar spin evolves from birth to the various end points of stellar evolution is an actively debated topic. Binary interactions are often invoked to explain the existence of extremely fast-rotating stars ($v\sin\,i$ $\gtrsim$ 200 $km\,s^{-1}$). The primary mechanisms through which binaries can spin up st… ▽ More The initial distribution of rotational velocities of stars is still poorly known, and how the stellar spin evolves from birth to the various end points of stellar evolution is an actively debated topic. Binary interactions are often invoked to explain the existence of extremely fast-rotating stars ($v\sin\,i$ $\gtrsim$ 200 $km\,s^{-1}$). The primary mechanisms through which binaries can spin up stars are tidal interactions, mass transfer, and possibly mergers. To evaluate these scenarios, we investigated in detail the evolution of three known fast-rotating stars in short-period spectroscopic and eclipsing binaries, namely HD 25631, HD 191495, and HD 46485, with primaries of masses of 7, 15, and 24 $M_{\odot}$, respectively, with companions of $\sim1\,M_\odot$ and orbital periods of less than 7 days. These systems belong to a recently identified class of binaries with extreme mass ratios, whose evolutionary origin is still poorly understood. We evaluated in detail three scenarios that could explain the fast rotation observed in these binaries: it could be primordial, a product of mass transfer, or the result of a merger within an originally triple system. We computed grids of single and binary MESA models varying tidal forces and initial binary architectures to investigate the evolution and reproduce observational properties of these systems. We find that, because of the extreme mass-ratio between binary components, tides have a limited impact, regardless of the prescription used, and that the observed short orbital periods are at odds with post-mass-transfer scenarios. The most likely scenario to form such young, rapidly rotating, and short-period binaries is primordial rotation, implying that the observed binaries are pre-interaction ones. These binaries show that the initial spin distribution of massive stars can have a wide range of rotational velocities. △ Less

Submitted 20 January, 2024; originally announced January 2024.

Comments: Accepted for publication in the Astronomy & Astrophysics journal. 20 pages, 17 figures. Reproducible at https://github.com/NikolayBritavskiyAstro/fast_rotating_binaries; the input and output MESA data are available at https://zenodo.org/records/10479754

Fluorescent Fe K line emission of gamma Cas stars I. Do gamma Cas stars host propelling neutron stars?

Authors: Gregor Rauw

Abstract: Gamma~Cas stars are early-type Be stars that exhibit an unusually hard and bright thermal X-ray emission. One of the proposed scenarios to explain these properties postulates the existence of a neutron star companion in the propeller stage, during which the magnetosphere of a rapidly rotating neutron star repels infalling material. To test this model, we examined the fluorescent Fe K$α$ emission l… ▽ More Gamma~Cas stars are early-type Be stars that exhibit an unusually hard and bright thermal X-ray emission. One of the proposed scenarios to explain these properties postulates the existence of a neutron star companion in the propeller stage, during which the magnetosphere of a rapidly rotating neutron star repels infalling material. To test this model, we examined the fluorescent Fe K$α$ emission line at $\sim 6.4$\,keV in the X-ray spectra of $γ$~Cas stars, which offers a powerful diagnostic of both the primary source of hard X-rays and the reprocessing material. We computed synthetic line profiles of the fluorescent Fe K$α$ emission line in the framework of the propelling neutron star scenario. Two reservoirs of material contribute to the fluorescence in this case: the Be circumstellar decretion disk and a shell of cool material that surrounds the shell of X-ray-emitting plasma around the putative propelling neutron star. We analysed the synthetic line profiles and expected equivalent widths of the lines for three well-studied $γ$~Cas stars. The predicted line strengths fall short of the observed values by at least an order of magnitude. Pushing the model parameters to reproduce the observed line strengths led to column densities towards the primary X-ray source that exceed the observationally determined values by typically a factor of 20, and would further imply a higher X-ray luminosity than observed. The strengths of the observed Fe K$α$ fluorescent emission lines in $γ$~Cas stars are inconsistent with the expected properties of a propeller scenario as proposed in the literature. △ Less

Submitted 19 December, 2023; originally announced December 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Colliding winds in WR21 and WR31 -- I. The X-ray view

Authors: Yael Naze, Gregor Rauw, Rachel Johnson, Eric Gosset, Jennifer L. Hoffman

Abstract: WR21 and WR31 are two WR+O binaries with short periods, quite similar to the case of V444 Cyg. The XMM-Newton observatory has monitored these two objects and clearly revealed phase-locked variations as expected from colliding winds. The changes are maximum in the soft band (0.5--2.keV, variations by a factor 3--4) where they are intrinsically linked to absorption effects. The increase in absorptio… ▽ More WR21 and WR31 are two WR+O binaries with short periods, quite similar to the case of V444 Cyg. The XMM-Newton observatory has monitored these two objects and clearly revealed phase-locked variations as expected from colliding winds. The changes are maximum in the soft band (0.5--2.keV, variations by a factor 3--4) where they are intrinsically linked to absorption effects. The increase in absorption due to the dense WR wind is confirmed by the spectral analysis. The flux maximum is however not detected exactly at conjunction with the O star in front but slightly afterwards, suggesting Coriolis deflection of the collision zone as in V444 Cyg. In the hard band (2.--10. keV), the variations (by a factor of 1.5--2.0) are much more limited. Because of the lower orbital inclinations, eclipses as observed for V444 Cyg are not detected in these systems. △ Less

Submitted 12 September, 2023; v1 submitted 1 September, 2023; originally announced September 2023.

Comments: accepted for publication by MNRAS

Extreme mass ratios and fast rotation in three massive binaries

Authors: Yael Naze, Nikolay Britavskiy, Gregor Rauw, Jonathan Labadie-Bartz, S. Simon-Diaz

Abstract: The origin of rapid rotation in massive stars remains debated, although binary interactions are now often advocated as a cause. However, the broad and shallow lines in the spectra of fast rotators make direct detection of binarity difficult. In this paper, we report on the discovery and analysis of multiplicity for three fast-rotating massive stars: HD25631 (B3V), HD191495 (B0V), and HD46485 (O7V)… ▽ More The origin of rapid rotation in massive stars remains debated, although binary interactions are now often advocated as a cause. However, the broad and shallow lines in the spectra of fast rotators make direct detection of binarity difficult. In this paper, we report on the discovery and analysis of multiplicity for three fast-rotating massive stars: HD25631 (B3V), HD191495 (B0V), and HD46485 (O7V). They display strikingly similar TESS light curves, with two narrow eclipses superimposed on a sinusoidal variation due to reflection effects. We complement these photometric data by spectroscopy from various instruments (X-Shooter, Espadons, FUSE...), to further constrain the nature of these systems. The detailed analyses of these data demonstrates that the companions of the massive OB stars have low masses (~1Msol) with rather large radii (2-4 Rsol) and low temperatures (<15 kK). These companions display no UV signature, which would exclude a hot subdwarf nature, but disentangling of the large set of X-Shooter spectra of HD25631 revealed the typical signature of chromospheric activity in the companion's spectrum. In addition, despite the short orbital periods (P=3-7d), the fast-rotating OB-stars still display non-synchronized rotation and all systems appear young (<20Myr). This suggests that, as in a few other cases, these massive stars are paired in those systems with non-degenerate, low-mass PMS companions, implying that fast rotation would not be a consequence of a past binary interactions in their case. △ Less

Submitted 7 August, 2023; v1 submitted 4 August, 2023; originally announced August 2023.

Comments: accepted for publication by MNRAS

Revisiting the orbital motion of WR 138

Authors: Gregor Rauw, Yaël Nazé, Eric Gosset

Abstract: The optical spectrum of WR 138 exhibits emission lines typical of a WN6o star and absorption lines from a rapidly-rotating OB star. Using a large set of spectroscopic data, we establish a new orbital solution of the WN6o star based on the radial velocities of highly-ionized nitrogen lines. We show that the WN6o star moves on a 4.3 yr orbit with a comparatively low eccentricity of 0.16. The radial… ▽ More The optical spectrum of WR 138 exhibits emission lines typical of a WN6o star and absorption lines from a rapidly-rotating OB star. Using a large set of spectroscopic data, we establish a new orbital solution of the WN6o star based on the radial velocities of highly-ionized nitrogen lines. We show that the WN6o star moves on a 4.3 yr orbit with a comparatively low eccentricity of 0.16. The radial velocities of the OB star display considerable scatter. Our best estimates of the velocities of He I absorption lines result in a mass-ratio of $m_{\rm WN6o}/m_{\rm OB} = 0.53 \pm 0.09$. We disentangle the spectra of the two stars and derive a projected rotational velocity of $v\,\sin{i} = 350 \pm 10$ km s$^{-1}$ for the OB star. Our best orbital parameters, combined with the Gaia parallax of WR 138, are at odds with a previous interferometric detection of the companion, suggesting that there is either a bias in this detection or that WR 138 is actually a triple system. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comments: accepted for publication by New Astronomy

New insight into the massive eccentric binary HD 165052: self-consistent orbital solution, apsidal motion, and fundamental parameters

Authors: Sophie Rosu, Edwin A. Quintero, Gregor Rauw, Philippe Eenens

Abstract: HD165052 is a short-period massive eccentric binary system that undergoes apsidal motion. As the rate of apsidal motion is directly related to the internal structure constants of the binary components, its study allows to get insight into the internal structure of the stars. We use medium- and high-resolution spectroscopic observations of HD165052 to provide constraints on the fundamental properti… ▽ More HD165052 is a short-period massive eccentric binary system that undergoes apsidal motion. As the rate of apsidal motion is directly related to the internal structure constants of the binary components, its study allows to get insight into the internal structure of the stars. We use medium- and high-resolution spectroscopic observations of HD165052 to provide constraints on the fundamental properties of the binary system and the evolutionary state of its components. We apply a spectral disentangling code to reconstruct artefact-free spectra of the individual stars and derive the radial velocities (RVs) at the times of the observations. We perform the first analysis of the disentangled spectra with the non-LTE model atmosphere code CMFGEN to determine the stellar properties. We derive the first self-consistent orbital solution of all existing RV data, including those reported in the literature, accounting for apsidal motion. We build, for the very first time, dedicated stellar evolution tracks with the Clés code requesting the theoretical effective temperatures and luminosities to match those obtained from our spectroscopic analysis. The binary system HD165052, consisting of an O6.5V((f)) primary and an O7V((f)) secondary, displays apsidal motion at a rate of (11.30+0.64-0.49)$°$yr$^{-1}$. Evolutionary masses are compared to minimum dynamical masses to constrain the orbital inclination. Evolutionary masses Mev,P=24.8$\pm$1.0M$_\odot$ and Mev,S=20.9$\pm$1.0M$_\odot$ and radii Rev,P=7.0+0.5-0.4R$_\odot$ and Rev,S=6.2+0.4-0.3R$_\odot$ are derived, and the inclination is constrained to 22.1$°\le i\le 23.3°$. Theoretical apsidal motion rates, derived assuming an age of 2.0+/-0.5 Myr for the binary, are in agreement with the observational determination. The agreement with theoretical apsidal motion rates enforces the inferred values of the evolutionary stellar masses and radii. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: 17 pages. arXiv admin note: text overlap with arXiv:2205.11207

Diving into the magnetosphere of the Of?p star HD108

Authors: Gregor Rauw, Yaël Nazé, Asif ud-Doula, Coralie Neiner

Abstract: We analyse optical and X-ray spectroscopy of the Of?p star HD108, known for its strong dipolar magnetic field and its optical line profile variability with a timescale of $54 \pm 3$ yrs, interpreted as the stellar rotation period. Optical emission lines have now recovered from their minimum emission state reached in 2007 - 2008. The variations of the equivalent width of the H$α$ emission provide c… ▽ More We analyse optical and X-ray spectroscopy of the Of?p star HD108, known for its strong dipolar magnetic field and its optical line profile variability with a timescale of $54 \pm 3$ yrs, interpreted as the stellar rotation period. Optical emission lines have now recovered from their minimum emission state reached in 2007 - 2008. The variations of the equivalent width of the H$α$ emission provide constraints on the inclination of the rotation axis ($i$) and the obliquity of the magnetic axis ($β$). The best agreement between model and observations is found for ($i$, $β$) pairs with $i + β\simeq 85^{\circ}$ and $i \in [30^{\circ},55^{\circ}]$. The Balmer emission lines display stochastic variability at the $\sim 5$ % level on timescales of a few days. TESS photometry unveils transient modulations on similar timescales in addition to prominent red noise variations. A Chandra X-ray observation of December 2021, when the star was at a higher emission level, indicates a slight increase of the flux and a spectral hardening compared to the August 2002 XMM-Newton observation, taken near minimum emission state. Magnetohydrodynamic simulations are used to compute synthetic X-ray spectra. With our current best estimate of the $\dot{M}_{B=0}$ mass-loss rate, the simulated X-ray luminosity and spectral energy distribution agree very well with the observations. Finally, the radial velocities vary on a period of 8.5 years with a peak-to-peak amplitude of 10 - 11 km s$^{-1}$, suggesting orbital motion with an unseen companion of at least 4 M$_{\odot}$. △ Less

Submitted 2 March, 2023; originally announced March 2023.

Comments: Accepted for publication in MNRAS

The Athena X-ray Integral Field Unit: a consolidated design for the system requirement review of the preliminary definition phase

Authors: Didier Barret, Vincent Albouys, Jan-Willem den Herder, Luigi Piro, Massimo Cappi, Juhani Huovelin, Richard Kelley, J. Miguel Mas-Hesse, Stéphane Paltani, Gregor Rauw, Agata Rozanska, Jiri Svoboda, Joern Wilms, Noriko Yamasaki, Marc Audard, Simon Bandler, Marco Barbera, Xavier Barcons, Enrico Bozzo, Maria Teresa Ceballos, Ivan Charles, Elisa Costantini, Thomas Dauser, Anne Decourchelle, Lionel Duband , et al. (274 additional authors not shown)

Abstract: The Athena X-ray Integral Unit (X-IFU) is the high resolution X-ray spectrometer, studied since 2015 for flying in the mid-30s on the Athena space X-ray Observatory, a versatile observatory designed to address the Hot and Energetic Universe science theme, selected in November 2013 by the Survey Science Committee. Based on a large format array of Transition Edge Sensors (TES), it aims to provide sp… ▽ More The Athena X-ray Integral Unit (X-IFU) is the high resolution X-ray spectrometer, studied since 2015 for flying in the mid-30s on the Athena space X-ray Observatory, a versatile observatory designed to address the Hot and Energetic Universe science theme, selected in November 2013 by the Survey Science Committee. Based on a large format array of Transition Edge Sensors (TES), it aims to provide spatially resolved X-ray spectroscopy, with a spectral resolution of 2.5 eV (up to 7 keV) over an hexagonal field of view of 5 arc minutes (equivalent diameter). The X-IFU entered its System Requirement Review (SRR) in June 2022, at about the same time when ESA called for an overall X-IFU redesign (including the X-IFU cryostat and the cooling chain), due to an unanticipated cost overrun of Athena. In this paper, after illustrating the breakthrough capabilities of the X-IFU, we describe the instrument as presented at its SRR, browsing through all the subsystems and associated requirements. We then show the instrument budgets, with a particular emphasis on the anticipated budgets of some of its key performance parameters. Finally we briefly discuss on the ongoing key technology demonstration activities, the calibration and the activities foreseen in the X-IFU Instrument Science Center, and touch on communication and outreach activities, the consortium organisation, and finally on the life cycle assessment of X-IFU aiming at minimising the environmental footprint, associated with the development of the instrument. Thanks to the studies conducted so far on X-IFU, it is expected that along the design-to-cost exercise requested by ESA, the X-IFU will maintain flagship capabilities in spatially resolved high resolution X-ray spectroscopy, enabling most of the original X-IFU related scientific objectives of the Athena mission to be retained. (abridged). △ Less

Submitted 28 November, 2022; v1 submitted 30 August, 2022; originally announced August 2022.

Comments: 48 pages, 29 figures, Accepted for publication in Experimental Astronomy with minor editing

The X-ray emission of Be+stripped star binaries

Authors: Yael Naze, Gregor Rauw, Myron A. Smith, Christian Motch

Abstract: Using observations from Chandra, Swift and XMM-Newton, we investigate the high-energy properties of all known (18) Be+sdO systems as well as 7 additional Be binaries suspected to harbour stripped stars. The observed X-ray properties are found to be similar to those observed for other Be samples. The vast majority of these systems (15 out of 25) display very faint (and soft) X-ray emission, and six… ▽ More Using observations from Chandra, Swift and XMM-Newton, we investigate the high-energy properties of all known (18) Be+sdO systems as well as 7 additional Be binaries suspected to harbour stripped stars. The observed X-ray properties are found to be similar to those observed for other Be samples. The vast majority of these systems (15 out of 25) display very faint (and soft) X-ray emission, and six others are certainly not bright X-ray sources. Only two systems display gamma-Cas characteristics (i.e. bright and hard X-rays), and one of them is a new detection: HD37202 (zeta Tau). It presents an extremely hard spectrum, due to a combination of high temperature and high absorption (possibly due to its high inclination). In parallel, it may be noted that the previously reported cyclic behaviour of this Be star has disappeared in recent years. Instead, shorter cycles and symmetric line profiles are observed for the Halpha line. It had been recently suggested that the peculiar X-ray emissions observed in gamma-Cas stars could arise from a collision between the disk of a Be star and the wind of its hot, stripped-star companion. The small fraction of gamma-Cas analogs in this sample, as well as the properties of the known companions of the gamma-Cas cases (low mass or not extremely hot, contrary to predictions), combined to the actual stripped-star and colliding-wind empirical knowledge, make the disk-wind collision an unlikely scenario to explain the gamma-Cas phenomenon. △ Less

Submitted 25 August, 2022; v1 submitted 8 August, 2022; originally announced August 2022.

Comments: accepted for publication by MNRAS

The X-ray emission of gamma Cassiopeiae during the 2020-2021 disc eruption

Authors: Gregor Rauw, Yaël Nazé, Christian Motch, Myron A. Smith, Joan Guarro Fló, Raimundo Lopes de Oliveira

Abstract: gamma Cas is known for its hard and intense X-ray emission that could trace accretion by a compact companion, wind interaction with a hot sub-dwarf companion, or magnetic interaction between the star and its Be decretion disc. These scenarios should lead to diverse dependences of the hard X-ray emission on disc density. We collected X-ray observations of gamma Cas during an episode of enhanced dis… ▽ More gamma Cas is known for its hard and intense X-ray emission that could trace accretion by a compact companion, wind interaction with a hot sub-dwarf companion, or magnetic interaction between the star and its Be decretion disc. These scenarios should lead to diverse dependences of the hard X-ray emission on disc density. We collected X-ray observations of gamma Cas during an episode of enhanced disc activity around January 2021. We investigate the variations in the disc properties using time series of dedicated optical spectroscopy and existing broadband photometry. Epoch-dependent Doppler maps of the H-alpha, H-beta, and He I 5876 emission lines are built to characterise the emission regions in velocity space. We analyse 4 XMM-Newton observations taken at key phases of the enhanced disc activity episode. Archival data are used to study the long-term correlation between optical and X-ray emission. Optical spectroscopy unveils an increase in the radial extent of the emission regions during the episode of enhanced disc activity, whilst no increase in the V-band flux is recorded. Doppler maps do not reveal any stable feature in the disc resulting from the putative action of the companion on the outer parts of the Be disc. No increase in the hard emission is observed in relation to the enhanced disc activity. However, at two occasions, the soft X-ray emission of gamma Cas is strongly attenuated, suggesting more efficient obscuration by a large flaring Be disc. There is a strong correlation between the long-term variations in the X-ray flux and in the V-band photometry. The observed behaviour of gamma Cas suggests no direct link between the properties of the outer regions of the Be disc and the hard X-ray emission, but favours a link between the level of X-ray emission and the properties of the inner part of the Be disc. These results thus disfavour an accretion or colliding wind scenario. △ Less

Submitted 17 June, 2022; originally announced June 2022.

Comments: Accepted for publication in A&A

Journal ref: A&A 664, A184 (2022)

Eight years of TIGRE robotic spectroscopy: Operational experience and selected scientific results

Authors: J. N. González-Perez, M. Mittag, J. H. M. M. Schmitt, K. -P. Schröder, D. Jack, G. Rauw, Y. Nazé

Abstract: TIGRE (Telescopio Internacional de Guanajuato Robótico Espectroscópico) has been operating in fully robotic mode in the Observatory La Luz (Guanajuato, Mexico) since the end of 2013. With its sole instrument, HEROS, an échelle spectrograph with a spectral resolution R~20000, TIGRE has collected more than 48000 spectra of 1151 different sources with a total exposure time of more than 11000 hours in… ▽ More TIGRE (Telescopio Internacional de Guanajuato Robótico Espectroscópico) has been operating in fully robotic mode in the Observatory La Luz (Guanajuato, Mexico) since the end of 2013. With its sole instrument, HEROS, an échelle spectrograph with a spectral resolution R~20000, TIGRE has collected more than 48000 spectra of 1151 different sources with a total exposure time of more than 11000 hours in these eight years. Here we briefly describe the system and the upgrades performed during the last years. We present the statistics of the weather conditions at the La Luz Observatory, emphasizing the characteristics that affect the astronomical observations. We evaluate the performance and efficiency of TIGRE, both optical and operational, and describe the improvements of the system implemented to optimize the telescope's performance and meet the requirements of the astronomer in terms of timing constraints for the observations and the quality of the spectra. We describe the actions taken to slow down the optical efficiency loss due to the aging of the optical surfaces as well as the upgrades of the scheduler and the observing procedures to minimize the time lost due to interrupted observations or observations that do not reach the required quality. Finally, we highlight a few of the main scientific results obtained with TIGRE data. △ Less

Submitted 6 June, 2022; originally announced June 2022.

Comments: 26 pages, 14 figures, accepted in Frontiers in Astronomy and Space Sciences

Apsidal motion in massive eccentric binaries: The case of CPD-41$^\circ$7742, and HD 152218 revisited

Authors: Sophie Rosu, Gregor Rauw, Yaël Nazé, Eric Gosset, Christiaan Sterken

Abstract: This paper is part of a study of the apsidal motion in close eccentric massive binary systems, which aims to constrain the internal structure of the stars. We focus on the binary CPD-41°7742 and revisit HD152218. Independent studies of CPD-41°7742 in the past showed large discrepancies in the longitude of periastron of the orbit, hinting at the presence of apsidal motion. We perform a consistent a… ▽ More This paper is part of a study of the apsidal motion in close eccentric massive binary systems, which aims to constrain the internal structure of the stars. We focus on the binary CPD-41°7742 and revisit HD152218. Independent studies of CPD-41°7742 in the past showed large discrepancies in the longitude of periastron of the orbit, hinting at the presence of apsidal motion. We perform a consistent analysis of all observational data, explicitly accounting for the apsidal motion. We use spectroscopic and photometric observations of CPD-41°7742 to infer values for the fundamental stellar and binary parameters. We apply a disentangling method to the spectra to simultaneously derive the RVs at the times of observation and reconstruct the individual stellar spectra. These are analysed by means of CMFGEN to determine the stellar properties. We determine the apsidal motion rate in two ways: We complement our RVs with those reported in the literature and we use the phase shifts between the primary and secondary eclipses. The light curves are analysed by means of Nightfall to constrain the orbital inclination and stellar masses. Stellar structure and evolution models are constructed with Cles. Different prescriptions for the mixing are adopted. Photometric data of HD152218 are analysed and stellar evolution models are built. The binary CPD-41°7742 displays apsidal motion at a rate of 15.38+0.42-0.51° yr-1. Our analysis of the observational data of CPD-41°7742 explicitly accounting for the apsidal motion allows us to explain the discrepancy in periastron longitudes. The age estimates are in good agreement with estimates obtained for other massive binaries in NGC 6231. This study confirms the need for enhanced mixing in the stellar evolution models to reproduce the observational stellar properties; this points towards larger convective cores than usually considered. △ Less

Submitted 25 May, 2022; v1 submitted 23 May, 2022; originally announced May 2022.

Comments: arXiv admin note: text overlap with arXiv:2202.02012

Journal ref: A&A 664, A98 (2022)

X-ray emission of massive stars and their winds

Authors: Gregor Rauw

Abstract: Most types of massive stars display X-ray emission that is affected by the properties of their stellar winds. Single non-magnetic OB stars have an X-ray luminosity that scales with their bolometric luminosity and their emission is thought to arise from a distribution of wind-embedded shocks. The lack of significant short-term stochastic variability indicates that the winds consist of a large numbe… ▽ More Most types of massive stars display X-ray emission that is affected by the properties of their stellar winds. Single non-magnetic OB stars have an X-ray luminosity that scales with their bolometric luminosity and their emission is thought to arise from a distribution of wind-embedded shocks. The lack of significant short-term stochastic variability indicates that the winds consist of a large number of independent fragments. Detailed variability studies unveiled a connection between the photosphere and the wind: well-studied O-type stars exhibit a ~ 10% modulation of their emission on timescales consistent with the rotation period, and a few early B-type pulsators display ~ 10% modulations of their X-ray flux with the pulsation period. Unlike OB stars, their evolved descendants (WR and LBV stars) lack a well-defined relation between their X-ray and bolometric luminosities, and several subcategories of objects remain undetected. These properties most likely stem from the combined effects of wind optical depth and wind velocity. Magnetic OB stars display an enhanced X-ray emission frequently modulated by the rotation of the star. These properties are well explained by the magnetically confined wind shock model and an oblique magnetic rotator configuration. Some massive binaries display phase-dependent excess emission arising from the collision between the winds of the binary components. Yet, the majority of the massive binaries do not show such an emission, probably as a consequence of radiative cooling of the shock-heated plasma. Finally, a growing subset of the Be stars, the so-called gamma Cas stars, feature an unusually hard and strong thermal X-ray emission that varies over a wide range of timescales. Several scenarios have been proposed to explain these properties, but the origin of the phenomenon remains currently one of the major unsolved puzzles in stellar X-ray astrophysics. △ Less

Submitted 31 March, 2022; originally announced March 2022.

Comments: Invited chapter of the Handbook of X-ray and Gamma-ray Astrophysics

Apsidal motion in massive eccentric binaries in NGC 6231: The case of HD 152219

Authors: S. Rosu, G. Rauw, M. Farnir, M. -A. Dupret, A. Noels

Abstract: The measurement of the apsidal motion in close eccentric massive binary systems provides essential information to probe the internal structure of the stars that compose the system. Following the determination of the fundamental stellar and binary parameters, we make use of the tidally induced apsidal motion to infer constraints on the internal structure of the stars composing the binary system HD1… ▽ More The measurement of the apsidal motion in close eccentric massive binary systems provides essential information to probe the internal structure of the stars that compose the system. Following the determination of the fundamental stellar and binary parameters, we make use of the tidally induced apsidal motion to infer constraints on the internal structure of the stars composing the binary system HD152219. The extensive set of spectroscopic, photometric, and radial velocity observations allows us to constrain the fundamental parameters of the stars together with the rate of apsidal motion of the system. Stellar structure and evolution models are further built with the Clés code testing different prescriptions for the internal mixing occurring inside the stars. The effect of stellar rotation axis misalignment with respect to the normal to the orbital plane on our interpretation of the apsidal motion in terms of internal structure constants is investigated. Made of an O9.5 III primary star (M1 = 18.64+/-0.47M${_\odot}$, R1 = 9.40+0.14-0.15R${_\odot}$, Teff,1 = 30900+/-1000 K) and a B1-2 V-III secondary star (M2 = 7.70+/-0.12M${_\odot}$, R2 = 3.69+/-0.06R${_\odot}$, Teff,2 = 21697+/-1000 K), the binary system HD152219 displays apsidal motion at a rate (1.198+/-0.300)°yr-1. The weighted-average mean of the internal structure constant of the binary system is inferred: k2 = 0.00173+/-0.00052. For the Clés models to reproduce the k2-value of the primary star, a significant enhanced mixing is required, notably through the turbulent mixing, but at the cost that other stellar parameters cannot be reproduced simultaneously. The difficulty to reproduce the k2-value simultaneously with the stellar parameters as well as the incompatibility between the age estimates of the primary and secondary stars are indications that some physics of the stellar interior are still not completely understood. △ Less

Submitted 4 February, 2022; originally announced February 2022.

Comments: 30 pages

Journal ref: A&A 660, A120 (2022)

X-ray response to disk evolution in two gamma-Cas stars

Authors: Yael Naze, Gregor Rauw, Terrence Bohlsen, Bernard Heathcote, Padric Mc Gee, Paulo Cacella, Christian Motch

Abstract: The Halpha emission of a set of southern gamma-Cas stars was monitored since 2019, with the aim of detecting transition events and examining how their peculiar X-ray emission would react in such cases. Two stars, HD119682 and V767Cen, were found to display slowly decreasing disk emissions. These decreases were not perfectly monotonic and several temporary and limited re-building events were observ… ▽ More The Halpha emission of a set of southern gamma-Cas stars was monitored since 2019, with the aim of detecting transition events and examining how their peculiar X-ray emission would react in such cases. Two stars, HD119682 and V767Cen, were found to display slowly decreasing disk emissions. These decreases were not perfectly monotonic and several temporary and limited re-building events were observed. For HD119682, the emission component in Halpha disappeared in mid-July 2020. In X-rays, the X-ray flux was twice smaller than recorded two decades ago but of a similar level as observed a decade ago. The X-ray flux decreased over the campaign by 30%, but the hardness remained similar in datasets of all epochs. In particular, the gamma-Cas character remained as clear as before even when there was no trace of disk emission in the Halpha line. For V767Cen, the full disappearance of disk emission in Halpha never occurred. We followed closely a disk rebuilding event, but no significant change in flux or hardness was detected. These behaviours are compared to those of other gamma-Cas stars and their consequences on the X-ray generation are discussed. △ Less

Submitted 1 February, 2022; originally announced February 2022.

Comments: accepted by MNRAS

Science with the Ultraviolet Explorer (UVEX)

Authors: S. R. Kulkarni, Fiona A. Harrison, Brian W. Grefenstette, Hannah P. Earnshaw, Igor Andreoni, Danielle A. Berg, Joshua S. Bloom, S. Bradley Cenko, Ryan Chornock, Jessie L. Christiansen, Michael W. Coughlin, Alexander Wuollet Criswell, Behnam Darvish, Kaustav K. Das, Kishalay De, Luc Dessart, Don Dixon, Bas Dorsman, Kareem El-Badry, Christopher Evans, K. E. Saavik Ford, Christoffer Fremling, Boris T. Gansicke, Suvi Gezari, Y. Goetberg , et al. (31 additional authors not shown)

Abstract: UVEX is a proposed medium class Explorer mission designed to provide crucial missing capabilities that will address objectives central to a broad range of modern astrophysics. The UVEX design has two co-aligned wide-field imagers operating in the FUV and NUV and a powerful broadband medium resolution spectrometer. In its two-year baseline mission, UVEX will perform a multi-cadence synoptic all-sky… ▽ More UVEX is a proposed medium class Explorer mission designed to provide crucial missing capabilities that will address objectives central to a broad range of modern astrophysics. The UVEX design has two co-aligned wide-field imagers operating in the FUV and NUV and a powerful broadband medium resolution spectrometer. In its two-year baseline mission, UVEX will perform a multi-cadence synoptic all-sky survey 50/100 times deeper than GALEX in the NUV/FUV, cadenced surveys of the Large and Small Magellanic Clouds, rapid target of opportunity followup, as well as spectroscopic followup of samples of stars and galaxies. The science program is built around three pillars. First, UVEX will explore the low-mass, low-metallicity galaxy frontier through imaging and spectroscopic surveys that will probe key aspects of the evolution of galaxies by understanding how star formation and stellar evolution at low metallicities affect the growth and evolution of low-metallicity, low-mass galaxies in the local universe. Such galaxies contain half the mass in the local universe, and are analogs for the first galaxies, but observed at distances that make them accessible to detailed study. Second, UVEX will explore the dynamic universe through time-domain surveys and prompt spectroscopic followup capability will probe the environments, energetics, and emission processes in the early aftermaths of gravitational wave-discovered compact object mergers, discover hot, fast UV transients, and diagnose the early stages of stellar explosions. Finally, UVEX will become a key community resource by leaving a large all-sky legacy data set, enabling a wide range of scientific studies and filling a gap in the new generation of wide-field, sensitive optical and infrared surveys provided by the Rubin, Euclid, and Roman observatories. This paper discusses the scientific potential of UVEX, and the broad scientific program. △ Less

Submitted 17 January, 2023; v1 submitted 30 November, 2021; originally announced November 2021.

Comments: 69 pages, 43 figures

Velocity monitoring of gamma-Cas stars reveals their binarity status

Authors: Yael Naze, Gregor Rauw, Stefan Czesla, Myron A. Smith, Jan Robrade

Abstract: The binary status of gamma-Cas stars has been discussed while theoretically examining the origin of their peculiar X-ray emission. However, except in two cases, no systematic radial velocity monitoring of these stars had been undertaken yet to clarify their status. We now fill this gap using TIGRE, CARMENES, and UVES high-resolution spectroscopy. Velocities were determined for 16 stars, revealing… ▽ More The binary status of gamma-Cas stars has been discussed while theoretically examining the origin of their peculiar X-ray emission. However, except in two cases, no systematic radial velocity monitoring of these stars had been undertaken yet to clarify their status. We now fill this gap using TIGRE, CARMENES, and UVES high-resolution spectroscopy. Velocities were determined for 16 stars, revealing shifts and/or changes in line profiles. The orbit of six new binaries could be determined: the long periods (80-120d) and small velocity amplitudes (5-7km/s) suggest low mass companions (0.6-1 M$_{\odot}$). The properties of the known gamma-Cas binaries appear similar to those of other Be systems, with no clear-cut separation between them. One of the new systems is a candidate for a rare case of quadruple system involving a Be star. Five additional gamma-Cas stars display velocity variations compatible with the presence of companions, but no orbital solution could yet be formally established for them hence they only receive the status of "binary candidate". △ Less

Submitted 18 November, 2021; originally announced November 2021.

Comments: accepted by MNRAS

Red noise and pulsations in evolved massive stars

Authors: Yael Naze, Gregor Rauw, Eric Gosset

Abstract: We examine high-cadence space photometry taken by the Transiting Exoplanet Survey Satellite (TESS) of a sample of evolved massive stars (26 Wolf-Rayet stars and 8 Luminous Blue Variables or candidate LBVs). To avoid confusion problems, only stars without bright Gaia neighbours and without evidence of bound companions are considered. This leads to a clean sample, whose variability properties should… ▽ More We examine high-cadence space photometry taken by the Transiting Exoplanet Survey Satellite (TESS) of a sample of evolved massive stars (26 Wolf-Rayet stars and 8 Luminous Blue Variables or candidate LBVs). To avoid confusion problems, only stars without bright Gaia neighbours and without evidence of bound companions are considered. This leads to a clean sample, whose variability properties should truly reflect the properties of the WR and LBV classes. Red noise is detected in all cases and its fitting reveals characteristics very similar to those found for OB-stars. Coherent variability is also detected for 20% of the WR sample. Most detections occur at moderately high frequency (3--14/d), hence are most probably linked to pulsational activity. This work doubles the number of WRs known to exhibit high-frequency signals. △ Less

Submitted 14 January, 2021; originally announced January 2021.

Comments: accepted by MNRAS - the arxiv version has some figures in low-resolution

LIFELINE: The program for the simulation of the X-ray line profiles in massive colliding wind binaries

Authors: E. Mossoux, G. Rauw

Abstract: The study of the X-ray line profiles produced by massive colliding wind binaries is a powerful tool for the characterisation of the stellar winds. We built a self-consistent program for the computation of line profiles named LIFELINE. The resulting theoretical profiles can be compared to the line profile that will be observed with future high-resolution X-ray spectrographs to retrieve the characte… ▽ More The study of the X-ray line profiles produced by massive colliding wind binaries is a powerful tool for the characterisation of the stellar winds. We built a self-consistent program for the computation of line profiles named LIFELINE. The resulting theoretical profiles can be compared to the line profile that will be observed with future high-resolution X-ray spectrographs to retrieve the characteristics of the stellar winds generating them. We considered a grid of 780 O-type binaries and computed, for each of them, the wind velocity distribution of each star, taking the impact of the radiation pressure and gravity force of the companion star into account. We then computed the characteristics of the wind shock region and followed the emitted photons towards the observer to compute their absorption. Finally, the Fe K line profiles near 6.7keV were constructed from the distribution of the photons as a function of the radial velocities of their emitting region. LIFELINE can be used to compare the theoretical line profiles to the observed ones or to compute theoretical profiles for a new binary system. We highlight the results for three systems. While the line profiles created in adiabatic wind collision regions are quite simple, the line profiles arising from regions in the radiative regime, as found in short-period binaries, are more sophisticated notably because of the Coriolis effect on the shape of the shock. The predicted differences in line morphology between systems with different wind properties are quite significant, allowing a detailed comparison between the theoretical profiles and those that will be observed with future high-resolution X-ray spectrometers. △ Less

Submitted 30 November, 2020; originally announced November 2020.

Comments: 19 pages, 10 figures, accepted by A&A

Journal ref: A&A 646, A89 (2021)

Similar but different: the varied landscape of Onfp/Oef stars variability

Authors: Gregor Rauw, Yael Naze

Abstract: The Oef category gathers rapidly rotating and evolved O-stars displaying a centrally reversed He II 4686 emission line. The origin of the variability of their photospheric and wind spectral lines is debated, with rotational modulation or pulsations as the main contenders. To shed new light on this question, we analysed high-quality and high-cadence TESS photometric time series for five Oef stars.… ▽ More The Oef category gathers rapidly rotating and evolved O-stars displaying a centrally reversed He II 4686 emission line. The origin of the variability of their photospheric and wind spectral lines is debated, with rotational modulation or pulsations as the main contenders. To shed new light on this question, we analysed high-quality and high-cadence TESS photometric time series for five Oef stars. We also collected a new time series of spectra for one target (lambda Cep) which had been the subject of specific debates in the last years. These observations reveal the variety of Oef behaviours. While space-based photometric data reveal substantial red noise components in all targets, only zeta Pup seems to display a long-lived periodicity. In our sample, stars exhibit a dominant signal at low frequencies but it appears relatively short-lived. This is reminiscent of rotational modulations by transient photospheric spots, though this scenario is challenged by the case of HD 14442, whose 1.230 d$^{-1}$ signal significantly exceeds the critical rotational frequency. In parallel, no evidence of persistent p-mode non-radial pulsations is found in either photometry or spectroscopy of the stars, only temporary excitation of g-mode pulsations could offer an alternative explanation for the dominant signals. Finally, the revised luminosities of the stars using GAIA-DR2 show that they are not all supergiants as zeta Pup. The question then arises whether the Oef peculiarity denotes a homogeneous class of objects after all. △ Less

Submitted 23 October, 2020; originally announced October 2020.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society

Journal ref: MNRAS, 500, 2096-2111 (2021)

The Massive Binary System 9 Sgr Revisited: New Insights into Disentangling Methods

Authors: Edwin A. Quintero, Philippe Eenens, Gregor Rauw

Abstract: Disentangling techniques are often needed to obtain the spectra of the individual components of binary or multiple systems. A thorough analysis of the shift-and-add algorithm of Marchenko, Moffat, & Eenens (1998) reveals that in many cases the line fluxes are poorly reproduced and spurious wings appear. The causes of these discrepancies are discussed and a new disentangling package, QER20, is pres… ▽ More Disentangling techniques are often needed to obtain the spectra of the individual components of binary or multiple systems. A thorough analysis of the shift-and-add algorithm of Marchenko, Moffat, & Eenens (1998) reveals that in many cases the line fluxes are poorly reproduced and spurious wings appear. The causes of these discrepancies are discussed and a new disentangling package, QER20, is presented which significantly reduces these errors and vastly increases the performance. When applied to the massive binary 9 Sgr, our new code yields line fluxes which are notably different from those previously published and lead us to revise the spectral classification to slightly earlier subtypes: O3V((f +)) for the primary and O5V((f)) for the secondary. We show that with the MME98 algorithm the classification of massive stars in binaries can be off by several subtypes whilst there are no such errors when the QER20 package is used. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: Accepted for publication in Astronomical Notes, 12 pages, 10 figures

Journal ref: Astronomische Nachrichten, 341(6-7), 628-637 (2020)

Apsidal motion in the massive binary HD 152248 -- Constraining the internal structure of the stars

Authors: S. Rosu, A. Noels, M. -A. Dupret, G. Rauw, M. Farnir, S. Ekström

Abstract: Apsidal motion in massive eccentric binaries offers precious information about the internal structure of the stars. This is especially true for twin binaries consisting of two nearly identical stars. We make use of the tidally induced apsidal motion in the twin binary HD152248 to infer constraints on the internal structure of the O7.5 III-II stars composing this system. We build stellar evolution… ▽ More Apsidal motion in massive eccentric binaries offers precious information about the internal structure of the stars. This is especially true for twin binaries consisting of two nearly identical stars. We make use of the tidally induced apsidal motion in the twin binary HD152248 to infer constraints on the internal structure of the O7.5 III-II stars composing this system. We build stellar evolution models with the code Clés assuming different prescriptions for the internal mixing occurring inside the stars. We identify the models that best reproduce the observationally determined present-day properties of the components of HD152248, as well as their $k_2$, and the apsidal motion rate of the system. We analyse the impact of some poorly constrained input parameters, including overshooting, turbulent diffusion, and metallicity. We further build 'single' and 'binary' GENEC models that account for stellar rotation to investigate the impacts of binarity and rotation. We discuss some effects that could bias our interpretation of the apsidal motion in terms of the internal structure constant. Reproducing the observed $k_2$ value and rate of apsidal motion simultaneously with the other stellar parameters requires a significant amount of internal mixing or enhanced mass-loss. The results suggest that a single-star evolution model is sufficient to describe the physics inside this binary system. Qualitatively, the high turbulent diffusion required to reproduce the observations could be partly attributed to stellar rotation. Higher-order terms in the apsidal motion are negligible. Only a very severe misalignment of the rotation axes could significantly impact the rate of apsidal motion, but such a high misalignment is highly unlikely in such a binary system. We infer an age estimate of $5.15\pm0.13$ Myr for the binary and initial masses of $32.8\pm0.6$ M$_\odot$ for both stars. △ Less

Submitted 18 September, 2020; originally announced September 2020.

Comments: 16 pages, A&A accepted

Journal ref: A&A 642, A221 (2020)

TESS lightcurves of gamma-Cas stars

Authors: Yael Naze, Gregor Rauw, Andrzej Pigulski

Abstract: gamma-Cas stars constitute a subgroup of Be stars showing unusually hard and bright X-ray emission. In search for additional peculiarities, we analyzed the TESS lightcurves of 15 gamma-Cas analogs. Their periodograms display broad frequency groups and/or narrow isolated peaks, often superimposed over red noise. The detected signals appear at low frequencies, with few cases of significant signals b… ▽ More gamma-Cas stars constitute a subgroup of Be stars showing unusually hard and bright X-ray emission. In search for additional peculiarities, we analyzed the TESS lightcurves of 15 gamma-Cas analogs. Their periodograms display broad frequency groups and/or narrow isolated peaks, often superimposed over red noise. The detected signals appear at low frequencies, with few cases of significant signals beyond 5/d (and all of them are faint). The signal amplitudes, and sometimes the frequency content, change with time, even in the absence of outburst events. On the basis of their optical photometric variability, gamma-Cas stars reveal no distinctive behaviour and thus appear similar to Be stars in general. △ Less

Submitted 19 August, 2020; originally announced August 2020.

Comments: accepted by MNRAS - the arxiv version has figures in low-resolution

Journal ref: MNRAS, 498, 3171-3183 (2020)

Irregular emission cycles in the Oe star HD 60848

Authors: Gregor Rauw, Yael Naze, Fran Campos, Joan Guarro Flo, Umberto Sollecchia

Abstract: We analyse the spectroscopic and photometric variability of the Oe star HD 60848 over the last twenty five years. The spectra reveal recurrent, but irregular cycles of increased circumstellar emission lines. These cycles are highly asymmetric displaying a slow increase over about 6 years, followed by a fast decay within about 6 months. Our analysis focuses on the most recent cycle (2013 - 2020). T… ▽ More We analyse the spectroscopic and photometric variability of the Oe star HD 60848 over the last twenty five years. The spectra reveal recurrent, but irregular cycles of increased circumstellar emission lines. These cycles are highly asymmetric displaying a slow increase over about 6 years, followed by a fast decay within about 6 months. Our analysis focuses on the most recent cycle (2013 - 2020). The equivalent width and velocity separation of the emission peaks indicate variations of the outer disk radius by a factor ~ 2.2, although the variability appears more complex than expected from first principle relations for optically thin Keplerian disks. We observe a time delay between the variations of the strengths of He I 5876 on the one hand and H-alpha and H-beta on the other hand. We interpret this behaviour in a two-step disk growth scenario, where the disk first expands radially before its density increases. A difference in behaviour is also seen between H-alpha and the H I Paschen lines, with the latter displaying a more symmetric cycle, similar to the photometric variability. The rather fast decays of the H-alpha emission observed in 2001, 2009 and 2018 - 2019 suggest that the strong radiation field and early spectral type of the star lead to a faster dissipation of the disk than in later-type Be stars, as theoretically expected. We discuss X-ray observations of the star both during a high and a low-emission state. The X-ray spectrum is soft at both epochs, and the X-ray fluxes are only marginally different and remain consistent with the canonical Lx/Lbol relation of O-type stars. These results indicate that the circumstellar decretion disk of HD 60848 has essentially no impact on the star's X-ray emission, and that the latter most likely arises inside the stellar wind. △ Less

Submitted 17 August, 2020; originally announced August 2020.

Comments: New Astronomy, accepted

Journal ref: New Astronomy, 83, id 101462 (14p, 2021)

Exocomets: A spectroscopic survey

Authors: I. Rebollido, C. Eiroa, B. Montesinos, J. Maldonado, E. Villaver, O. Absil, A. Bayo, H. Canovas, A. Carmona, Ch. Chen, S. Ertel, Th. Henning, D. P. Iglesias, R. Launhardt, R. Liseau, G. Meeus, A. Moór, A. Mora, J. Olofsson, G. Rauw, P. Riviere-Marichalar

Abstract: While exoplanets are now routinely detected, the detection of small bodies in extrasolar systems remains challenging. Since the discovery of sporadic events interpreted as exocomets (Falling Evaporating Bodies) around $β$ Pic in the early 80s, only $\sim$20 stars have been reported to host exocomet-like events. We aim to expand the sample of known exocomet-host stars, as well as to monitor the hot… ▽ More While exoplanets are now routinely detected, the detection of small bodies in extrasolar systems remains challenging. Since the discovery of sporadic events interpreted as exocomets (Falling Evaporating Bodies) around $β$ Pic in the early 80s, only $\sim$20 stars have been reported to host exocomet-like events. We aim to expand the sample of known exocomet-host stars, as well as to monitor the hot-gas environment around stars with previously known exocometary activity. We have obtained high-resolution optical spectra of a heterogeneous sample of 117 main-sequence stars in the spectral type range from B8 to G8. The data have been collected in 14 observing campaigns expanding over 2 years from both hemispheres. We have analysed the Ca ii K&H and Na i D lines in order to search for non-photospheric absorptions originated in the circumstellar environment, and for variable events that could be caused by outgassing of exocomet-like bodies. We have detected non-photospheric absorptions towards 50% of the sample, attributing a circumstellar origin to half of the detections (i.e. 26% of the sample). Hot circumstellar gas is detected in the metallic lines inspected via narrow stable absorptions, and/or variable blue-/red-shifted absorption events. Such variable events were found in 18 stars in the Ca ii and/or Na i lines; 6 of them are reported in the context of this work for the first time. In some cases the variations we report in the Ca ii K line are similar to those observed in $β$ Pic. While we do not find a significant trend with the age or location of the stars, we do find that the probability of finding CS gas in stars with larger vsin i is higher. We also find a weak trend with the presence of near-infrared excess, and with anomalous ($λ$ Boo-like) abundances, but this would require confirmation by expanding the sample. △ Less

Submitted 24 March, 2020; originally announced March 2020.

Comments: Accepted in A&A on March 23rd 2020

Journal ref: A&A 639, A11 (2020)

Search for non-thermal X-ray emission in the colliding wind binary Cyg OB2 #8A

Authors: E. Mossoux, J. M. Pittard, G. Rauw, Y. Nazé

Abstract: Cyg OB2 #8a is a massive O-type binary displaying strong non-thermal radio emission. Owing to the compactness of this binary, emission of non-thermal X-ray photons via inverse Compton scattering is expected. We first revised the orbital solution for Cyg OB2 #8a using new optical spectra. We then reduced and analysed X-ray spectra obtained with XMM-Newton, Swift, INTEGRAL, and NuSTAR. The analysis… ▽ More Cyg OB2 #8a is a massive O-type binary displaying strong non-thermal radio emission. Owing to the compactness of this binary, emission of non-thermal X-ray photons via inverse Compton scattering is expected. We first revised the orbital solution for Cyg OB2 #8a using new optical spectra. We then reduced and analysed X-ray spectra obtained with XMM-Newton, Swift, INTEGRAL, and NuSTAR. The analysis of the XMM-Newton and Swift data allows us to better characterise the X-ray emission from the stellar winds and colliding winds region at energies below 10 keV. We confirm the variation of the broad-band light curve of Cyg OB2 #8a along the orbit with, for the first time, the observation of the maximum emission around phase 0.8. The minimum ratio of the X-ray to bolometric flux of Cyg OB2 #8a remains well above the level expected for single O-type stars, indicating that the colliding wind region is not disrupted during the periastron passage. The analysis of the full set of publicly available INTEGRAL observations allows us to refine the upper limit on the non-thermal X-ray flux of the Cyg OB2 region between 20 and 200 keV. Two NuSTAR observations (phases 0.028 and 0.085) allow us to study the Cyg OB2 #8a spectrum up to 30 keV. These data do not provide evidence of the presence of non-thermal X-rays, but bring more stringent constraints on the flux of a putative non-thermal component. Finally, we computed, thanks to a new dedicated model, the anisotropic inverse Compton emission generated in the wind shock region. The theoretical non-thermal emission appears to be compatible with observational limits and the kinetic luminosity computed from these models is in good agreement with the unabsorbed flux observed below 10 keV. △ Less

Submitted 20 March, 2020; originally announced March 2020.

Comments: 12 pages, 8 figures, A&A accepted

Journal ref: A&A 636, A109 (2020)

Let there be more variability in two gamma Cas stars

Authors: Yael Naze, Andrzej Pigulski, Gregor Rauw, Myron Smith

Abstract: We investigate the short-term optical variability of two gamma Cas analogs, pi Aqr and BZ Cru, thanks to intensive ground-based spectroscopic and space-borne photometric monitorings. For both stars, low-amplitude (mmag) coherent photometric variability is detected. The associated signals display long-term amplitude variations, as in other Be stars. However, these signals appear at high frequencies… ▽ More We investigate the short-term optical variability of two gamma Cas analogs, pi Aqr and BZ Cru, thanks to intensive ground-based spectroscopic and space-borne photometric monitorings. For both stars, low-amplitude (mmag) coherent photometric variability is detected. The associated signals display long-term amplitude variations, as in other Be stars. However, these signals appear at high frequencies, especially in pi Aqr, indicating p-modes with a high degree l, a quite unusual feature amongst Be stars. While BZ Cru presents only low-level spectral variability, without clear periodicity, this is not the case of pi Aqr. In this star, the dominant photometric frequencies, near ~12/d, are confirmed spectroscopically in separate monitorings taken during very different disk activity levels ; the spectroscopic analysis suggests a probable tesseral nature for the mode. △ Less

Submitted 28 February, 2020; originally announced February 2020.

Comments: accepted by MNRAS - the arxiv version has some figures in low-resolution

Journal ref: MNRAS, 494, 958-974 (2020)

Three discoveries of gamma Cas analogs from dedicated XMM-Newton observations of Be stars

Authors: Yael Naze, Christian Motch, Gregor Rauw, Shami Kumar, Jan Robrade, Raimundo Lopes de Oliveira, Myron A. Smith, Jose M. Torrejon

Abstract: In the last years, a peculiarity of some Be stars - their association with unusually hard and intense X-ray emission - was shown to extend beyond a mere few cases. In this paper, we continue our search for new cases by performing a limited survey of 18 Be stars using XMM-Newton. The targets were selected either on the basis of a previous X-ray detection (Exosat, ROSAT, XMM-slew survey) without spe… ▽ More In the last years, a peculiarity of some Be stars - their association with unusually hard and intense X-ray emission - was shown to extend beyond a mere few cases. In this paper, we continue our search for new cases by performing a limited survey of 18 Be stars using XMM-Newton. The targets were selected either on the basis of a previous X-ray detection (Exosat, ROSAT, XMM-slew survey) without spectral information available, or because of the presence of a peculiar spectral variability. Only two targets remain undetected in the new observations and three other stars only display faint and soft X-rays. Short-term and/or long-term variations were found in one third of the sample. The spectral characterization of the X-ray brightest 13 stars of the sample led to the discovery of three new gamma Cas (HD44458, HD45995, V558Lyr), bringing the total to 25 known cases, and another gamma Cas candidate (HD120678), bringing the total to 2. △ Less

Submitted 14 February, 2020; v1 submitted 13 February, 2020; originally announced February 2020.

Comments: accepted for publication in MNRAS

Journal ref: MNRAS, 493, 2511-2517 (2020)

Apsidal motion in the massive binary HD152248

Authors: S. Rosu, G. Rauw, K. E. Conroy, E. Gosset, J. Manfroid, P. Royer

Abstract: The eccentric massive binary HD152248 (also known as V1007 Sco), which hosts two O7.5 III-II(f) stars, is the most emblematic eclipsing O-star binary in the very young and rich open cluster NGC6231. Its properties render the system an interesting target for studying tidally induced apsidal motion. Measuring the rate of apsidal motion in such a binary system gives insight into the internal structur… ▽ More The eccentric massive binary HD152248 (also known as V1007 Sco), which hosts two O7.5 III-II(f) stars, is the most emblematic eclipsing O-star binary in the very young and rich open cluster NGC6231. Its properties render the system an interesting target for studying tidally induced apsidal motion. Measuring the rate of apsidal motion in such a binary system gives insight into the internal structure and evolutionary state of the stars composing it. A large set of optical spectra was used to reconstruct the spectra of the individual binary components and establish their radial velocities using a disentangling code. Radial velocities measured over seven decades were used to establish the rate of apsidal motion. We furthermore analysed the reconstructed spectra with the CMFGEN model atmosphere code to determine stellar and wind properties of the system. Optical photometry was analysed with the Nightfall binary star code. A complete photometric and radial velocity model was constructed in PHOEBE 2 to determine robust uncertainties. We find a rate of apsidal motion of $(1.843^{+0.064}_{-0.083})°$ yr$^{-1}$. The photometric data indicate an orbital inclination of $(67.6^{+0.2}_{-0.1})°$ and Roche-lobe filling factors of both stars of about 0.86. Absolute masses of $29.5^{+0.5}_{-0.4}$M$_\odot$ and mean stellar radii of $15.07^{+0.08}_{-0.12}$R$_\odot$ are derived for both stars. We infer an observational value for the internal structure constant of both stars of $0.0010\pm0.0001$. Our in-depth analysis of the massive binary HD152248 and the redetermination of its fundamental parameters can serve as a basis for the construction of stellar evolution models to determine theoretical rates of apsidal motion to be compared with the observational one. In addition, the system hosts two twin stars, which offers a unique opportunity to obtain direct insight into the internal structure of the stars. △ Less

Submitted 4 February, 2020; originally announced February 2020.

Comments: 22 pages, A&A accepted

Journal ref: A&A 635, A145 (2020)

Closing gaps to our origins. The UV window into the Universe

Authors: Ana I. Gomez de Castro, Martin A. Barstow, Fréderic Baudin, Stefano Benetti, Jean Claude Bouret, Noah Brosch, Domitilla de Martino, Giulio del Zanna, Chris Evans, Miriam García, Boris Gaensicke, Carolina Kehrig, Jon Lapington, Alain Lecavelier des Etangs, Giampiero Naletto, Yael Nazé, Coralie Neiner, Jonathan Nichols, Marina Orio, Isabella Pagano, Gregor Rauw, Steven Shore, Gagik Tovmasian, Asif ud-Doula, Kevin France , et al. (1 additional authors not shown)

Abstract: The investigation of the emergence of life is a major endeavour of science. Astronomy is contributing to it in three fundamental manners: (1) by measuring the chemical enrichment of the Universe, (2) by investigating planet formation and searching for exoplanets with signatures of life and, (3) by determining the abundance of aminoacids and the chemical routes to aminoacid and protein growth in as… ▽ More The investigation of the emergence of life is a major endeavour of science. Astronomy is contributing to it in three fundamental manners: (1) by measuring the chemical enrichment of the Universe, (2) by investigating planet formation and searching for exoplanets with signatures of life and, (3) by determining the abundance of aminoacids and the chemical routes to aminoacid and protein growth in astronomical bodies. This proposal deals with the first two. In the Voyage to 2050, the world-wide scientific community is getting equipped with large facilities for the investigation of the emergence of life in the Universe (i.e. VLT, JWST, ELT, GMT, TMT, ALMA, FAST, VLA, ATHENA, SKA) including the ESA's CHEOPS, PLATO and ARIEL missions. This white paper is a community effort to call for the development of a large ultraviolet optical observatory to gather fundamental data for this investigation that will not be accessible through other ranges of the electromagnetic spectrum. A versatile space observatory with UV sensitivity a factor of 50-100 greater than existing facilities will revolutionize our understanding of the pathway to life in the Universe. △ Less

Submitted 2 November, 2019; originally announced November 2019.

Comments: White Paper submitted to the European Space Agency call for Voyage 2050. This WP proposes the development of a European Ultraviolet Visible Observatory (EUVO) to study the origins of life and it is largely based in arXiv 1306.3358 (ApSS,354,229)

Surprises in the simultaneous X-ray and optical monitoring of pi Aqr

Authors: Yael Naze, G. Rauw, M. Smith

Abstract: To help constrain the origin of the peculiar X-ray emission of gamma Cas stars, we conducted a simultaneous optical and X-ray monitoring of pi Aqr in 2018. At that time, the star appeared optically bright and active, with a very strong Halpha emission. Our monitoring covers three 84d orbital cycles, allowing us to probe phase-locked variations as well as longer-term changes. In the new optical dat… ▽ More To help constrain the origin of the peculiar X-ray emission of gamma Cas stars, we conducted a simultaneous optical and X-ray monitoring of pi Aqr in 2018. At that time, the star appeared optically bright and active, with a very strong Halpha emission. Our monitoring covers three 84d orbital cycles, allowing us to probe phase-locked variations as well as longer-term changes. In the new optical data, the radial velocity variations seem to span a smaller range than previously reported, which might indicate possible biases. The X-ray emission is variable, but without any obvious correlation with orbital phase or Halpha line strength. Furthermore, the average X-ray flux and the relative range of flux variations are similar to those recorded in previous data, although the latter data were taken when the star was less bright and its disk had nearly entirely disappeared. Only the local absorption component in the X-ray spectrum appears to have strengthened in the new data. This absence of large changes in X-ray properties despite dramatic disk changes appears at odds with previous observations of other gamma Cas stars. It also constrains scenarios proposed to explain the gamma Cas phenomenon. △ Less

Submitted 20 November, 2019; v1 submitted 24 October, 2019; originally announced October 2019.

Comments: accepted for publication in A&A

Journal ref: A&A 632, A23 (2019)

A Gaia view of the two OB associations Cygnus OB2 and Carina OB1: The signature of their formation process

Authors: Beomdu Lim, Yael Naze, Eric Gosset, Gregor Rauw

Abstract: OB associations are the prime star forming sites in galaxies. However the detailed formation process of such stellar systems still remains a mystery. In this context, identifying the presence of substructures may help tracing the footprints of their formation process. Here, we present a kinematic study of the two massive OB associations Cygnus OB2 and Carina OB1 using the precise astrometry from t… ▽ More OB associations are the prime star forming sites in galaxies. However the detailed formation process of such stellar systems still remains a mystery. In this context, identifying the presence of substructures may help tracing the footprints of their formation process. Here, we present a kinematic study of the two massive OB associations Cygnus OB2 and Carina OB1 using the precise astrometry from the Gaia Data Release 2 and radial velocities. From the parallaxes of stars, these OB associations are confirmed to be genuine stellar systems. Both Cygnus OB2 and Carina OB1 are composed of a few dense clusters and a halo which have different kinematic properties: the clusters occupy regions of 5-8 parsecs in diameter and display small dispersions in proper motion, while the halos spread over tens of parsecs with a 2-3 times larger dispersions in proper motion. This is reminiscent of the so-called "line width-size" relation of molecular clouds related to turbulence. Considering that the kinematics and structural features were inherited from those of their natal clouds would then imply that the formation of OB associations may result from structure formation driven by supersonic turbulence, rather than from the dynamical evolution of individual embedded clusters. △ Less

Submitted 9 September, 2019; originally announced September 2019.

Comments: 20 pages, 12 figures, 4 tables, accepted for publication in MNRAS

Journal ref: MNRAS, 490, 440-454 (2019)

Evolution of the disk of pi Aqr: from near-disappearance to a strong maximum

Authors: Y. Naze, G. Rauw, J . Guarro Flo, A. De Bruin, O. Garde, O. Thizy, F. Houpert, E. Pollmann, C. J. Sawicki, M. Leonardi, M. Moll, C. T. Quandt, P. Berardi, T. Lester, P. Fosanelli, A. Favaro, J. N. Terry, K. Graham, B. Mauclaire, T. Bohlsen, M. Pujol, E. Bertrand, E. Bryssinck, V. Desnoux, P. Lailly , et al. (8 additional authors not shown)

Abstract: Some Be stars display important variability of the strength of the emission lines formed in their disk. This is notably the case of pi Aqr. We present here the recent evolution of the Be disk in this system thanks to spectra collected by amateur spectroscopists since the end of 2013. A large transition occurred: the emission linked to the Be disk nearly disappeared in January 2014, but the disk ha… ▽ More Some Be stars display important variability of the strength of the emission lines formed in their disk. This is notably the case of pi Aqr. We present here the recent evolution of the Be disk in this system thanks to spectra collected by amateur spectroscopists since the end of 2013. A large transition occurred: the emission linked to the Be disk nearly disappeared in January 2014, but the disk has recovered, with a line strength now reaching levels only seen during the active phase of 1950--1990. In parallel to this change in strength occurs a change of disk structure, notably involving the disappearance of the strong asymmetry responsible for the V/R modulation. △ Less

Submitted 21 June, 2019; originally announced June 2019.

Comments: accepted by New Astronomy

Journal ref: New Astronomy (vol 73, id 101279, 2019)

Variations on a theme -- the puzzling behaviour of Schulte 12

Authors: Yael Naze, Gregor Rauw, Stefan Czesla, Laurent Mahy, Fran Campos

Abstract: One of the first massive stars detected in X-rays, Schulte 12 has remained a puzzle in several aspects. In particular, its extreme brightness both in the visible and X-ray ranges is intriguing. Thanks to Swift and XMM-Newton observations covering ~5000d, we report the discovery of a regular 108d modulation in X-ray flux of unknown origin. The minimum in the high-energy flux appears due to a combin… ▽ More One of the first massive stars detected in X-rays, Schulte 12 has remained a puzzle in several aspects. In particular, its extreme brightness both in the visible and X-ray ranges is intriguing. Thanks to Swift and XMM-Newton observations covering ~5000d, we report the discovery of a regular 108d modulation in X-ray flux of unknown origin. The minimum in the high-energy flux appears due to a combination of increased absorption and decreased intrinsic emission. We examined in parallel the data from a dedicated spectroscopic and photometric monitoring in the visible and near-IR domains, complemented by archives. While a similar variation timescale is found in those data, they do not exhibit the strict regular clock found at high energies. Changes in line profiles cannot be related to binarity but rather correspond to non-radial pulsations. Considering the substantial revision of the distance of Schulte 12 from the second GAIA data release, the presence of such oscillations agrees well with the evolutionary status of Schulte 12, as it lies in an instability region of the HR diagram. △ Less

Submitted 3 June, 2019; originally announced June 2019.

Comments: accepted by A&A

Journal ref: A&A 627, A99 (2019)

Quest for the tertiary component in Cyg OB2 #5

Authors: Gregor Rauw, Yaël Nazé, Fran Campos

Abstract: The Cyg OB2 #5 system is thought to consist of a short-period (6.6 d) eclipsing massive binary orbited by an OB-star orbiting with a period of ~6.7 yr; these stars in turn are orbited by a distant early B-star with a period of thousands of years. However, while the inner binary has been studied many times, information is missing on the other stars, in particular the third star whose presence was i… ▽ More The Cyg OB2 #5 system is thought to consist of a short-period (6.6 d) eclipsing massive binary orbited by an OB-star orbiting with a period of ~6.7 yr; these stars in turn are orbited by a distant early B-star with a period of thousands of years. However, while the inner binary has been studied many times, information is missing on the other stars, in particular the third star whose presence was indirectly postulated from recurrent modulations in the radio domain. Besides, to this date, the X-ray light curve could not be fully interpreted, for example in the framework of colliding-wind emission linked to one of the systems. We obtained new optical and X-ray observations of Cyg OB2 #5, which we combined to archival data. We performed a thorough and homogeneous investigation of all available data, notably revisiting the times of primary minimum in photometry. In the X-ray domain, XMM-Newton provides scattered exposures over ~5000 d whilst Swift provides a nearly continuous monitoring for the last couple of years. Although the X-ray light curve reveals clear variability, no significant period can be found hence the high-energy emission cannot be explained solely in terms of colliding winds varying along either the short or intermediate orbits. The optical data reveal for the first time clear signs of reflex motion. The photometry indicates the presence of a 2366 d (i.e. 6.5 yr) period while the associated radial velocity changes are detected at the 3 sigma level in the systemic velocity of the He II 4686 emission line. With the revised period, the radio light curve is interpreted consistently in terms of a wind interaction between the inner binary and the tertiary star. From these optical and radio data, we derive constraints on the physical properties of the tertiary star and its orbit. △ Less

Submitted 16 May, 2019; originally announced May 2019.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A, 627, A2 (21p, 2019)

BRITE photometry of the massive post-RLOF system HD149404

Authors: G. Rauw, A. Pigulski, Y. Nazé, A. David-Uraz, G. Handler, F. Raucq, E. Gosset, A. F. J. Moffat, C. Neiner, H. Pablo, A. Popowicz, S. M. Rucinski, G. A. Wade, W. Weiss, K. Zwintz

Abstract: HD149404 is an evolved non-eclipsing O-star binary that has previously undergone a Roche lobe overflow interaction. Understanding some key properties of the system requires a determination of the orbital inclination and of the dimensions of the components. The BRITE-Heweliusz satellite was used to collect photometric data of HD149404. Additional photometry was retrieved from the SMEI archive. Thes… ▽ More HD149404 is an evolved non-eclipsing O-star binary that has previously undergone a Roche lobe overflow interaction. Understanding some key properties of the system requires a determination of the orbital inclination and of the dimensions of the components. The BRITE-Heweliusz satellite was used to collect photometric data of HD149404. Additional photometry was retrieved from the SMEI archive. These data were analysed using a suite of period search tools. The orbital part of the lightcurve was modelled with the nightfall binary star code. The Gaia-DR2 parallax of HD149404 was used to provide additional constraints. The periodograms reveal a clear orbital modulation of the lightcurve with a peak-to-peak amplitude near 0.04 mag. The remaining non-orbital part of the variability is consistent with red noise. The lightcurve folded with the orbital period reveals ellipsoidal variations, but no eclipses. The minimum when the secondary star is in inferior conjunction is deeper than the other minimum due to mutual reflection effects between the stars. Combined with the Gaia-DR2 parallaxes, the photometric data indicate an orbital inclination in the range of 23° to 31° and a Roche lobe filling factor of the secondary larger than or equal to 0.96. The luminosity of the primary star is consistent with its present-day mass, whereas the more evolved secondary appears overluminous for its mass. We confirm that the primary's rotation period is about half the orbital period. Both features most probably stem from the past Roche lobe overflow episode. △ Less

Submitted 8 November, 2018; originally announced November 2018.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 621, A15 (2019)