Open Access
Issue
A&A
Volume 692, December 2024
Article Number A36
Number of page(s) 15
Section Stellar atmospheres
DOI https://doi.org/10.1051/0004-6361/202451681
Published online 29 November 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1 Introduction

The investigation of the soft X-ray emission (below 10 keV) from massive (made of O-type and/or Wolf–Rayet stars) binaries constitutes a highly relevant approach to investigate shock physics, providing some insight into stellar winds. In these systems, the strong stellar winds are likely to collide, leading to strong shocks responsible for plasma heating up to several tens of 106 K. The properties of this thermal X-ray emission intimately depend on the shock parameters, such as pre-shock velocity and wind density. Even though in many cases the wind interaction is basically following a standard behavior well described by seminal theoretical works (Stevens et al. 1992; Pittard & Stevens 1997; Pittard & Parkin 2010), one can expect some specific cases to happen. This is typically what may happen for very asymmetric systems where one wind is much stronger than the other, potentially leading to wind-photosphere collisions. As a relevant tracer of phenomena occurring in such systems, soft X-ray spectroscopy constitutes a valid tool to characterize their behavior. The present study focuses on HD 93205, located in Trumpler 16, a stellar cluster targeted by several observation campaigns in X-rays over the past decades.

HD 93205 is an O-type system displaying a well-established double-lined spectroscopic binary nature (Walborn 1971, 1973). Conti & Walborn (1976) classified HD 93205 as O3V + O8V, and they derived a first orbital solution with a period of 6.08 d. The orbital solution was refined by Morrell et al. (2001), pointing to a quite eccentric orbit (see Table 1). HD 93205 is one of the earliest type stars identified in the Galaxy with a well-constrained orbital solution. The spectral type of the primary star in HD 93205 was later on revised as O3.5 V((f)) by Sota et al. (2014).

Many previous studies reported on the detection of HD 93205 in X-rays with a moderate X-ray luminosity of ~1032 erg s−1 (see e.g., Seward et al. 1979; Seward & Chlebowski 1982; Chlebowski & Garmany 1991; Corcoran et al. 1995). Clear phase-locked X-ray modulations were noticed by Corcoran (1996) using ROSAT observations. Using five X-ray observations from XMM-Newton, Antokhin et al. (2003) showed that the X-ray flux changes as a function of the separation between two binary components and found the source spectrum to be soft. Some X-ray variability, suspected to be related to colliding winds, were also reported by Albacete Colombo et al. (2003), Albacete-Colombo et al. (2008), Antokhin et al. (2008) and Gagné et al. (2011). Despite the numerous, previous X-ray variability reports, a detailed spectral and time analysis of the wind interaction in HD 93205 is lacking. In particular, the strong asymmetry of the system points to a potentially complex situation, where the primary wind may collide with the photosphere, at least at some orbital phases given the eccentricity of this short period system. Making use of a long and consistent time series of X-ray observations, our aim is to characterize its main properties in light of what has been learned over the past decades on the physical circumstances ruling the production of soft X-rays in massive O-type binaries, and especially in short period ones.

This paper is organized as follows. Section 1 describes the target and summarizes the outcome of previous works relevant for our purposes. Section 2 presents the X-ray observations and the processing of data used for the present study. The X-ray spectral analysis and light curve are detailed in Sects. 3 and 4, respectively. The X-ray luminosity of the system is discussed in Sect. 5. The physical discussion on the origin of the X-ray excess is present in Sect. 6, before presenting a physically valid scenario for HD 93205 in Sect. 7. Finally, our conclusions are drawn in Sect. 8.

Table 1

Main orbital parameters of HD 93205 as given by Morrell et al. (2001).

2 X-ray monitoring of HD 93205 with XMM-Newton

X-ray data of HD 93205 obtained with XMM-Newton (Jansen et al. 2001) from 2000 July to 2019 December has been analyzed. A log of total 31 epochs of utilized X-ray observations has been given in Table A.1. These observations were obtained for several other massive stars in Trumpler 16 cluster as main target and HD 93205 was observed as a field star by XMM-Newton using different configurations of the three European Photon Imaging Camera (EPIC) instruments, viz. MOS1, MOS2, and PN. Most of the observations were performed with either thick or medium optical blocking filters. The timing and spectral products have been extracted from data using the latest calibration files with SAS v20.0.0.

The SAS tasks EPCHAIN and EMCHAIN were used to process the raw EPIC Observation Data Files (ODF) for PN and MOS observations, respectively. The list of event files were generated considering only the events designated with pattern 0–4 for PN and 0–12 for MOS data using task EVSELECT. No pileup was noticed in any of the data set upon examining with EPATPLOT. The EPIC background is sometimes affected by high background X-ray photons and/or soft proton flares. In order to eliminate such instances, full-frame light curves were generated considering the single-events (PATTERN = 0) in > 10 keV energy range for MOS and 10-12 keV band for PN. The events with count rate more than 0.2 counts s−1 for MOS and 0.4 counts s−1 for PN background light curves were excluded. Rest of the events formed the good time intervals and were utilized for further analysis.

Figure 1 shows the X-ray image extracted from the observation ID 0311990101 in 0.3–12.0 keV energy range. The source coordinates were precisely determined using the source detection algorithm in SAS named as EDETECT_CHAIN. The location of HD 93205 has been detected at the coordinates RA (J2000) = 10:44:33.84 and Dec (J2000) = -59:44:13.92. HD 93205 has a neighboring source HD 93204 detected as close as ~18.3 arcsec at the position RA (J2000) = 10:44:32.40 and Dec (J2000) = −59:44:28.68 as shown in Fig. 1. HD 93204 is an O5.5 V((f)) star emitting in X-rays (Sota et al. 2014). We compared the location of both of these sources with the coordinates determined from the high resolution Chandra Carina catalog (Broos et al. 2011). Both the positions matched very well confirming our source detection approach. A circular region centered at the estimated location of HD 93205 with radius 25 arcsec has been selected to extract the light curves and spectra of the target source from EPIC CCDs. HD 93204 lies within the chosen circular region around HD 93205, therefore, the counts from another circular region of radius 8 arcsec centered at the location of HD 93204 were subtracted from the actual source region. This procedure helped to remove the X-ray contamination from the neighboring source to HD 93205 X-ray products.

The background light curve and spectrum were extracted from a 15 arcsec radius circular region located at the source-free positions surrounding HD 93205. The proximity of background region to the main target ensures roughly same off-axis angle of these two regions on the CCD plane. Chandra images were also investigated to verify that no other X-ray emitting object is present in our chosen background region due to its better resolution and sensitivity than XMM-Newton-EPIC. In order to make corrections in the extracted light curves to consider the good time intervals, dead time, exposure, point-spread function, and background subtraction, the task EPICLCCORR has been used. The source as well as the background spectra were generated by EVSELECT task. ARFGEN and RMFGEN were employed to calculate the dedicated ancillary response files (ARFs) and response matrix files (RMFs), respectively, necessary for the spectral fitting. The extracted spectra were backscaled using the task BACKSCALE. The energy channels in each of the extracted source spectrum were grouped using GRPPHA to have a minimum of ten counts per spectral bin. Further temporal and spectral analyses were performed using HEASoft version v6.29c.

thumbnail Fig. 1

X-ray image of the Eta Carinae region extracted from observation ID 0311990101 in the 0.3–12 keV energy range. The position of HD 93205 as been highlighted with a white circle of 25 arcsec radius. The X-ray flux of neighboring source HD 93204 has been removed from a circular region of 8 arcsec radius. An estimate of background was done from a circular region (cyan color) having radius 15 arcsec lying in the source-free area of the detector.

thumbnail Fig. 2

MOS1 and MOS2 spectra of HD 93205 jointly fitted with the two-temperature thermal plasma emission model from the observation ID 0145740301. The details of the spectral fitting are provided in Table A.2.

3 X-ray spectral analysis

The X-ray spectra of HD 93205 extracted consistently from XMM-Newton-EPIC displays typical features of an optically thin plasma heated in hot stars’ winds above 106 K. The obtained spectra were relatively soft with maximum emission located between 0.8–0.9 keV as shown in Fig. 2. There were no significant counts with energy more than 5 keV above background level. Therefore, the further spectral energy analysis was restricted to 0.3–5.0 keV energy range.

In order to characterize the X-ray spectral properties of HD 93205, and search for any variability with time, we considered optically thin thermal emission models. Such an approach is standard for the modeling of the X-ray emission from massive binaries (Pittard & Parkin 2010; De Becker 2015; Arora et al. 2019, 2021, 2024). We made use of the APEC emission model available in the X-ray spectral fitting package XSPEC (v.12.12.0). We fitted the spectrum with multiple APEC components but using only two emission components turned out to be sufficient to model the X-ray spectrum of HD 93205. The composite model was of the form PHABS(LSM)*PHABS(LOCAL)*(APEC+APEC). The multiplicative photoelectric absorption component PHABS has been used to model the effects of interstellar absorption and local absorption by the stellar wind material on the emitted X-rays. The interstellar column of absorbing matter was fixed to value of NHism=0.24×1022$N_H^{ism} = 0.24 \times {10^{22}}$ cm−2 obtained from Jenkins (2019). However, local column (NHlocal)$\left( {N_H^{local}} \right)$ was kept as a free parameter during the fitting procedure. The reduced χ2 minimization was used to test the goodness of fit. The model assumed solar abundances from Anders & Grevesse (1989) while performing the spectral fitting.

Independent modeling of the spectra obtained at different epochs is not straightforward, as it may lead to a significant dispersion in the best-fit parameters, depending upon the quality of individual data-sets. The PN spectra has more number of counts than other EPIC detectors, therefore, we fitted only PN spectra initially with the two-temperature composite model to explore the parameter space and achieve a first diagnostic of the full time series, without freezing any parameter (apart from interstellar absorption). Upon modeling the PN spectra obtained at each epoch separately, it was found that the average values of two temperatures corresponding to soft and hard emission were 0.16±0.01 (=Kt1) and 0.60±0.02 (=kT2), respectively. The achieved value of reduced χ2 was below 1.3 for all of the data sets. However, for few of the spectra the temperature corresponding to hard emission fitted at higher value around 1.5 keV which better modeled the high energy tail but overall it resulted in poor reduced χ2 than low kT2 model and it resulted in increased normalization constant of soft component (norm1) at the expense of normalization value from hard component (norm2). As it did not improve the physical relevance of the spectral fitting, we decided to proceed further with kT1 =0.2 keV and kT2=0.6 keV. These values are consistent with the ones estimated by Antokhin et al. (2008) and Ryspaeva & Kholtygin (2020).

We note that similar results were obtained for MOS1 and MOS2 spectra, so we decided to focus on the joint fitting of all EPIC spectra at each epoch. As the next step, we decided to freeze some parameters to constrain the exploration of the parameter space and achieve a consistent set of solutions for the full time series. All emission model temperatures were frozen (kT1 = 0.2 keV and kT2 = 0.6 keV), keeping the local absorption column and the two normalization parameters free. Following this approach, we obtained a consistent modeling series for 31 epochs, including cases where the number of counts were too low to obtain adequate results in our first unconstrained attempts. All our results are presented in Table A.2. Reduced χ2 values range between 0.9 and 1.5.

The MOS1 and MOS2 X-ray spectra of HD 93205 shown in Fig. 2 has been jointly modeled. In order to see the evolution of various spectral parameters over the complete orbit of HD 93205, we estimated its binary orbital phase according to the ephemeris provided by Morrell et al. (2001) as HJD = 2450499.089+6.080E. Zero phase corresponds to the time of periastron passage in this eccentric binary system. The best-fit parameters NHlocal$N_H^{local}$, norm1 and norm2 have been phase phase-folded in Fig. 3. The ratio norm1/norm2 plotted in the bottom panel of this figure provides hints about variation in hardness of spectrum at several epochs of observation. The observed and ISM-corrected X-ray flux was also estimated in 0.3–5.0 keV (broad, FBobs,ism${\rm{F}}_B^{obs,ism}$), 0.3–2.0 keV (soft, FSobs,ism${\rm{F}}_S^{obs,ism}$), and 2.0–5.0 keV (hard, FHobs,ism${\rm{F}}_H^{obs,ism}$) energy bands. The variation of X-ray flux as a function of the orbital phase is shown in Fig. 4 in different energy bands. The observed trend suggests that in the broad and soft energy band there is a variation in flux by about a factor of ~1.6 between its lower and higher value achieved after apastron and close to periastron, respectively. However, in hard energy band the ratio is ~2.4. The phase-locked X-ray flux variations are clearly in agreement with the modulations shown in the Fig. 3.

4 X-ray timing analysis

In order to verify the determination of orbital period of HD 93205, we have performed a Fourier transform (FT) of the flux measurements provided in Table A.2. We used both Lomb–Scargle (Lomb 1976; Scargle 1982; Horne & Baliunas 1986) and CLEAN (Roberts et al. 1987) periodograms. Both algorithms were applied on the flux values obtained at 31 epochs.

The Lomb–Scargle algorithm is effective in searching for periodic signals in unevenly time sampled light curves and as shown in upper panel of Fig. 5. It consists of several peaks along with one dominant peak at frequency 0.16445±0.00004 cycles d−1. In order to remove several other peaks which might be originating due to the aliasing, the CLEAN power spectrum was also calculated with a loop-gain of 0.1 and 100 iterations. The CLEANed periodogram showed a clear, distinguished peak at exactly the same frequency as of the Lomb-Scargle power spectrum in the bottom panel of Fig. 5. The peak corresponds to a period of 6.0807±0.0013 days in full consistency with the orbital period obtained by Morrell et al. (2001) through high-quality spectroscopic solution. The present period determination shows relevance of investigating massive stars multiplicity from X-ray time series analysis carried over a considerably long time baseline to the classical orbit determination techniques (e.g., see Arora et al. 2024).

HD 93205 was properly exposed mainly in MOS1 and MOS2, therefore providing us with adequate time series for light curve analysis. The background-subtracted X-ray light curves as obtained from MOS1 and MOS2 are shown in Fig. 6. These light curves were extracted in broad, soft, and hard energy bands. While plotting these light curves, we have considered the count rates estimated from the observations obtained with thick optical blocking filters as given in Table A.1. It has to be noted that the difference in the detector effective area is significant below 0.6–0.7 keV for MOS1. As noticed in Sect. 3, most of the X-ray emission from HD 93205 comes from the soft energy range therefore it becomes important to account for detector efficiency while analysing variation in raw count rate obtained from X-ray light curves. We could see clear distinction in the count rates taken with different filters for the same instrument. However, the detector effective area is taken into picture by using dedicated ARF files while spectral fitting, facilitating analysis of all data points together in Fig. 4. The light curves in Fig. 6 were orbital phase folded and each point corresponds to average count rate from individual observation ID.

The X-ray variability is seen in all the energy band light curves with an order of magnitude lower count rate in the hard band. The phase-locked modulation in count rate is almost similar in the broad and soft energy bands where it becomes maximum around orbital phase zero and drops gradually to minimum around apastron. This behavior is more systematic in broad and soft energy bands. The count rate enhanced by ~60% in these energy band from apastron to periastron, although variation in hard band is more than 200% for MOS1 light curves which has clearer variation than MOS2 light curves. Finally, it is emphasized that soft band accounts for more than 90% of the total counts detected in broad energy band in good agreement with the soft X-ray spectrum discussed in Sect. 3.

thumbnail Fig. 3

Variation of the local equivalent H-column density (NHlocal)$\left( {{\rm{N}}_H^{local}} \right)$, normalization constants (norm1 and norm2), and the ratio normı/norm2 corresponding to two thermal plasma emission components with the orbital phase of HD 93205 (see Table A.2).

5 X-ray luminosity

The fluxes corrected for interstellar absorption were converted to X-ray luminosities, assuming a distance of 2350 pc (Bailer-Jones et al. 2021). Given the flux values across the orbit quoted in Table A.2, we determined luminosities (LX) ranging between 5.2 × 1032 and 8.3 × 1032 erg s−1.

A comparison with the X-ray luminosity reported by Gagné et al. (2011) shows that their value (~3.6 × 1032erg s−1) is lower than the luminosity range reported on above by about a factor 1.5. Their Chandra measurement spanned an orbital phase range (0.21–0.38) away from the maximum of the light curve. We attribute this discrepancy to two reasons. First, Gagné et al. (2011) worked in the 0.5–8.0 keV energy range, while we used the 0.3–5.0 keV band. Since most of the X-ray emission from HD 93205 is soft, our estimate of the X-ray luminosity accounts for a non-negligible contribution that was not measured by Gagné et al. (2011). Second, our modeling points to a major contribution from a plasma with a temperature of about 0.2 keV with a weak contribution from a plasma at 0.6 keV (versus 0.23 and 0.74 keV for Gagné et al. 2011). Such a spectrum is especially bright in the low energy part of the X-ray spectrum that is also the most affected by interstellar absorption. Because of the slightly softer spectrum considered in our study, the correction for that interstellar absorption is more pronounced in our case, leading to ISM-corrected luminosities somewhat higher than reported by Gagné et al. (2011).

For O-type stars, the intrinsic X-ray emission from individual stellar winds is expected to be about 10”7 times the bolometric luminosity (Lbol). This luminosity ratio is valid for not too evolved O-type stars (Owocki et al. 2013; De Becker 2013). Given the values of parameters quoted in Table 2, the measured X-ray luminosity points to an X-ray over-luminosity factor in the range of 1.8–2.8. This translates into an excess luminosity LX,excess = LXLX,PLX,S = 2.3–5.4 × 1032 erg s−1 (where LX,P and LX,S are the individual X-ray luminosities from the primary and secondary winds, respectively). For the O3 primary, the thick wind may be opaque enough to lead to the luminosity ratio to be slightly lower than 10−7. Our estimate of LX,excess may thus be seen as a conservative lower limit. This excess emission produced on top of individual winds constitutes the core of our discussion in Sect. 6.

thumbnail Fig. 4

Modulation of the (a) observed (Fobs, left) and (b) ISM-corrected (Fism, right) X-ray flux in broad, soft, and hard energy bands from HD 93205 obtained after X-ray spectral fitting as a function of the binary orbital phase (see Table A.2).

thumbnail Fig. 5

Lomb–Scargle and CLEANed power spectra of observed flux (FBobs)$\left( {F_B^{obs}} \right)$ from HD 93205 in the 0.3-5.0 keV energy range. The frequency of the peak with the highest power is also mentioned.

6 Origin of the X-ray emission excess

6.1 Wind interaction

In a massive binary system, a significant excess X-ray emission is usually attributable to colliding winds (Stevens et al. 1992; Pittard & Parkin 2010). However, given the highly asymmetric nature of HD 93205, the primary wind may completely overwhelm the secondary one, therefore preventing such a wind-wind collision (WWC) to occur. This opens up the possibility for a wind-photosphere collision (WPC) to develop in the system.

We searched for stable solutions of the stagnation point on the basis of the ram pressure balance condition. Given the small size of the orbit, stellar winds do not reach their terminal velocity before colliding. We assumed a classical β-velocity law (with β = 0.8) V=V(1RD)β$V = {V_\infty }{\left( {1 - {R \over D}} \right)^\beta }$(1)

for both winds (where R and D are the stellar radius and distance available to the wind to be radiatively accelerated before colliding, respectively), along with wind parameters quoted in Table 2. For the evolution of the stellar separation along the orbit, we assumed parameters taken from Table 1. By a comparison between minimum masses and expected absolute masses given by Martins et al. (2005), we deduce an inclination of about 58° (averaging values from both stars), in fair agreement with the 55° suggested by Morrell et al. (2001). This results in an absolute separation of 58.8 R, that has to be modulated by the eccentricity depending on the orbital phase. Our analysis resulted in a lack of stable stagnation point between the two stars at any orbital phase, pointing to the very likely scenario of a collision between the primary’s (stronger) wind and the secondary’s photosphere. However, we note that different assumptions on the wind parameters (including in particular a downward revision of the primary mass loss rate) may allow a stable WWC to form quite close to the secondary’s surface. This deserves to be qualified by a more detailed discussion including the influence of radiative effects.

Table 2

Adopted parameters for both components of HD 93205, based on Martins et al. (2005) and Muijres et al. (2012).

thumbnail Fig. 6

Folded X-ray light curves of HD 93205 as observed by XMM-Newton using (a) MOS1 (left) and (b) MOS2 (right) in broad, soft, and hard energy bands. The average count rate of each individual light curve observed with thick optical blocking filter has been plotted here.

6.2 Influence of radiative effects

The dynamics of the incoming primary wind is very likely to be affected by radiative effects. First of all, in strongly asymmetric systems such as HD 93205, one may be dealing with sudden radiative braking (Gayley et al. 1997). This process considers that the primary wind is suddenly decelerating while getting very close to the secondary’s surface. This may prevent the WPC to occur, leading alternatively to a WWC close to the secondary’s surface. The study by Gayley et al. (1997) consists in a one-dimensional approach that therefore ignores orbital effects and assumes axial symmetry about the line ofsep-aration between the stars. Such a symmetry would be broken in case of significant orbital distortion that tends to increase as the orbital period decreases. A valid indicator of this would be the orbital-to-wind velocity ratio: the greater the ratio, the stronger is the departure from axial symmetry. Despite of this limitation, this approach has the benefit to provide us with some analytical tools to investigate semi-quantitatively the influence of this radiative effect on our system at every orbital phase. A full three-dimensional analysis is well beyond the scope of the the present study.

In order to check for the relevance of the WPC and WWC scenarios, we made use of the formalism presented by Gayley et al. (1997). Unless specified differently, all parameters in equations below are expressed in cgs units. We shall first define the braking coefficient (Cb), Cb=(VPvth)(LSκe2πVP2cD)1αPP/vα${C_b} = \left( {{{{V_P}} \over {{v_{th}}}}} \right){\left( {{{{L_S}{\kappa _e}} \over {2\pi V_P^2cD}}} \right)^{1 - \alpha }}P_{P/v}^{ - \alpha }$(2)

where VP is the primary wind velocity, υth is the plasma thermal velocity, LS is the secondary bolometric luminosity, κe is the free electron cross section per gram, D is the stellar separation, α is the index of the line list power distribution of the standard CAK wind theory (Castor et al. 1975). PP/V is the ratio of the primary wind momentum flux to the secondary radiative momentum flux (ṀP C VP/LS).

The so-called braking radius (xb) is expressed as follows, xb=[ 1+(dη)(1α)/(1+α) ]1${x_b} = {\left[ {1 + {{\left( {{d \over \eta }} \right)}^{(1 - \alpha )/(1 + \alpha )}}} \right]^{ - 1}}$(3)

with η=1α1+α(11+αkCb)1/(1α)d$\eta = {{1 - \alpha } \over {1 + \alpha }}{\left( {{1 \over {1 + \alpha }}k{C_b}} \right)^{1/(1 - \alpha )}}d$(4)

where k is the force multiplier parameter from the CAK theory, and d is the stellar separation expressed in units of the secondary’s stellar radius (d = D/RS).

According to Gayley et al. (1997), radiative braking cannot prevent photospheric collision provided the following condition is fulfilled, Fbr=xbd<1${F_{br}} = {x_b}d < 1$(5)

where Fbr will be referred to as the braking parameter. Adopting α = 0.6 and k = 0.15 (values valid for O-type stars adopted by Gayley et al. 1997), along with κe = 0.37ɡ−1 (not specified by Gayley et al. 1997, but determined from their results), and making use of parameters quoted in Tables 1 ad 2, we checked for the validity of the condition expressed by Eq. (5) along the full orbit. As shown in Fig. 7, the braking parameter is lower than one over the full orbit (solid line). In order to check for the robustness of this result, we repeated the calculation adopting other values for k and α. Fbr tends to go above 1.0 when α goes above 0.6, especially close to apastron. For k = 0.2, Fbr is above 1.0 between orbital phases 0.3 and 0.7 (dashed line in Fig. 7). We note that Gayley et al. (1997) considered k values above 0.15 in the case of a Wolf–Rayet wind, due to the higher metallicity as compared to their O-type progenitors, while Abbott (1982) proposed values of 0.174 and 0.178 for O-stars with effective temperature of 40000 and 50 000 K, respectively (more typical of the HD 93205’s primary). The latter values lead to a braking parameter greater than one in a significant part of the orbit.

Another radiative effect worth considering is radiative inhibition (Stevens & Pollock 1994). This process considers the mutual influence of the radiation pressure of both stars in a massive binary system on their wind material. Basically, the incoming primary wind is likely to be slower toward the direction of the secondary due to the continuous contribution of the secondary radiation field on the primary wind equation of motion. Such a process is known to be able to reduce substantially the pre-shock velocity of a colliding wind, therefore significantly reducing the post-shock temperature of the X-ray emitting plasma. In case of substantial inhibition, the effective wind velocity could drop by 50% (see examples in Stevens & Pollock 1994).

Finally, a likely radiative effect is that of self-regulated shocks (Parkin & Sim 2013). In this case, the X-ray emission produced by the post-shock plasma ionizes the incoming wind ahead of the shock, leading to significant inhibition of the wind acceleration. The effect is twofold: (i) it can reduce the pre-shock velocity, therefore reducing the post-shock plasma temperature, and (ii) the ionization close to the shock in the pre-shock flow can make the line-driving mechanism much less efficient, thus potentially suppressing sudden radiative breaking. As a result, this process has the potential to extend the range of orbital separations where a WPC occurs. Consequently, even cases with CAK parameters in favor of sudden radiative breaking around apastron (Fig. 7, dotted line) may actually lead to a WPC along the full orbit.

Summarizing, radiative effects play a key role in the nature and properties of the wind interaction. The detailed analysis of their combined action would require a full hydro-radiative modeling, based on robust stellar and wind parameters for both components of the system. However, this can’t prevent us from exploring potential scenarios expected to be constrained by our X-ray measurements. According to our analysis, the WPC scenario is worth considering. However, depending on the adopted initial parameters, one cannot fully reject the idea of a WWC. This is especially true if one reduces the primary’s wind velocity and/or mass loss rate, and alternative assumptions are adopted for the line-driving parameters affecting the occurrence of radiative braking. Finally, as a result of the significant eccentricity of the system, the pre-shock conditions are changing along the orbit. We thus also have to consider a hybrid collision scenario (HCP) where the system is switching between a WPC (around periastron) and a WWC scenario (around apastron).

thumbnail Fig. 7

Plot of the braking parameter (Fbr) as a function of the orbital phase. The solid and dotted line stand for two different assumptions for the k parameter. The horizontal dashed line stands for the critical value above which the WPC scenario is not valid.

6.3 Wind-photosphere collision

6.3.1 X-ray emission from a WPC

In the WPC case, a shock warped about the secondary star is formed. The pre-shock flow is made up by the incoming primary wind, and the heated post-shock gas is between the shock front and the secondary’s surface. As a first step, we achieved some quantitative predictions of the X-ray luminosity from the WPC (LX,WPC) based on the approach proposed by Usov (1992, their Eq. (81)). The predicted X-ray luminosity is expressed as follows: LX,WPC=5.8×1032M˙P,52RS,123D134VP,8(erg s1)${L_{X,WPC}} = 5.8 \times {10^{32}}{{\dot M_{P, - 5}^2R_{S,12}^3} \over {D_{13}^4{V_{P,8}}}}\quad \left( {erg{{\rm{s}}^{ - 1}}} \right)$(6)

where P,−5 is the primary wind mass loss rate in units of 10−5 M yr−1, RS,12 is the secondary’s radius in units of 1012 cm, and VP,8 is the pre-shock velocity of the primary’s wind in units of 108 cm s−1. The separation D13 (in units of 1013 cm) is determined as a function of the orbital phase on the basis of the parameters provided in Table 1.

We stress that Eq. (6) established by Usov (1992) considers the thermal X-ray emission produced by the post-shock plasma is continuum free-free emission only, and completely neglects the contribution from emission lines. Equation (6) thus underestimates the actual X-ray luminosity by a factor that depends on the post-shock plasma temperature, and thus on the pre-shock velocity, as expressed below (Stevens et al. 1992): kT1.17VP,82.$kT \approx 1.17V_{P,8}^2.$(7)

We quantified the relative importance of emission lines making use of thermal emission models in XSPEC. We measured the theoretical X-ray flux between 0.3 and 10.0 keV for a series of thermal emission models with consistent normalization parameter, for plasma temperatures (kT) ranging between 0.25 to 10.0 keV (by steps of 0.25 keV). We repeated the procedure using two models offered in XSPEC: (i) apec, accounting for both free-free continuum and emission lines (total emission), and (ii) nlapec, accounting free-free continuum with no contribution from emission lines. We computed the total-to-continuum flux ratio at each temperature, as shown in Fig. 8. The contribution from spectral lines clearly dominates at low plasma temperature, and drops quickly to reach a level of about 10% of the continuum contribution above 2 keV. This correction has thus to be provided to LX,WPC depending on the post-shock temperature.

The pre-shock velocity is varying significantly as a function of the orbital phase, as a consequence of the eccentricity of the short period orbit where the primary wind collides the photosphere while still in its acceleration region. Figure 9 shows the pre-shock velocity at the position of the secondary (middle panel, solid line), along with the post-shock plasma temperature (lower panel, solid line) as a function of the orbital phase. For illustration purpose, dashed lines for the velocity and temperature illustrate the hypothetical situation of a significant drop in the pre-shock velocity due to high amplitude radiative inhibition (see Sect. 6.2).

Based on the considerations on the pre-shock velocity, Eq. (6) allows us to predict the X-ray luminosity from the WPC region, as shown in Fig. 10. According to Fig. 8, the range of temperatures displayed in the bottom panel of Fig. 9 suggests a quite hard X-ray emission, dominated by the free-free continuum with a contribution from spectral lines of some additional 10%. We note however that these temperature are clearly quite excessive regarding the plasma temperature values we reported on in Sect. 3. This point will be addressed in Sect. 7.

thumbnail Fig. 8

Plot of the total-to-continuum X-ray flux ratio as a function of plasma temperature.

thumbnail Fig. 9

Plot of the stellar separation (upper panel), the pre-shock velocity (middle panel), and the post-shock temperature (lower panel) as a function of the orbital phase. Solid lines stand for a first approximation of the pre-shock velocity while dashed lines illustrate the case of strong radiative inhibition.

6.3.2 Post-shock temperature

According to our discussion in Sect. 6, HD 93205 may give rise to a WPC in a significant part of its orbit, or even along its full orbit. Focusing on this scenario, we explored the question of the likely shape of the shock front to investigate in particular the impact of shock obliquity on the distribution of the pre-shock velocities (and resulting post-shock temperatures) along the shock surface. Our methodology is presented in Appendix A.

On the basis of normal pre-shock velocities determined as a function of θ (Eqs. (A.6) to (A.8)), we determined the expected post-shock temperature along the shock surface using Eq. (7). We repeated the calculation considering both periastron and apastron, in both un-inhibited and strongly radiatively inhibited (wind velocity dropped by a factor 2) cases. We restricted the calculation up to off-axis angle (θ) values corresponding to x = 0 (see Fig. A.1). This explains why the range of θ values displayed in Fig. 11 is different for the two considered orbital phases (the shock is seen under a lower off-axis angle at apastron as compared to periastron). Depending on the orbital phase, on the considered θ value, and on the level of radiative inhibition, the range of normal pre-shock velocities is quite high, that is, from about 500 to 3200 km s−1. This leads to an even wider range of post-shock temperatures, from about 0.5 keV to 14keV.

Because of the strong effect of obliquity along the shock front, the resulting average post-shock temperature is significantly lower than the one usually measured at the apex of the shock (along the line of centers). In an attempt to achieve a view of the averaged post-shock temperature more typical of the full emitting shock, we computed the surface of revolution of the shock curve (Eq. (A.1)) between x = 0 and x = RS + Δ, yielding a value of Ashock =830R2${A_{shock{\rm{ }}}} = 830\,{\rm{R}}_ \odot ^2$ (by steps of δ x = 0.01 R). The surface averaged kT was computed using this relation, kT=iAi(kT)iiAi=1Ashock iAi(kT)i$\langle kT\rangle = {{\sum\limits_i {{A_i}} {{(kT)}_i}} \over {\sum\limits_i {{A_i}} }} = {1 \over {{A_{shock{\rm{ }}}}}}\sum\limits_i {{A_i}} {(kT)_i}$(8)

where Ai is the annular element of surface at position xi. As each surface element increases proportionally to y as one moves from the shock apex to more off-axis contributions, the relative weight of cooler plasma is high enough to drop the average temperature significantly below that post-shock gas close to the line of centers. As a result, 〈kT〉 goes from 3.0 to 8.0 keV in the non-inhibited case at periastron and apastron, respectively. Assuming severe radiative inhibition (see above), these values drop down to 0.8 and 2.0 keV. In addition, our basic approach considers only the emission for the positive side of the x-axis. On the side of the back hemisphere of the secondary, some emission is also expected with a high obliquity, that is, with a lower plasma temperature. However, we stress that for negative values of x, the primary wind flow line gets quickly close to the asymptote of the shock surface, leading to a simple advection of pre-shock material along the shock without feeding any further the shock physics at the origin of the X-ray emission. This asymptotic flow case is especially quickly reached for shorter separations. We thus anticipate that the X-ray contributions that are not explicitly considered above should mainly lead to a slight drop of the average plasma temperature away from periastron. In addition, away from the line of centers, the influence of the secondary wind can be sufficiently strong to contribute to the shaping of the shock front. Dealing with these aspects is however beyond the scope of this paper, as this would require a full hydrodynamical and radiative modeling of the system.

Finally, the likely wind velocity parameter space can encompass lower values if for instance the assumed terminal velocity is too high. In addition, the assumed β parameter of the velocity law (Eq. (1)) has some influence on the results. For instance, for β = 1.0, 〈kT〉 goes down to 2.4 (resp. 7.6) at periastron (resp. apastron) for the un-inhibited case. For severe inhibition (see above), these values drop by a factor 4.

thumbnail Fig. 10

Plot of the predicted X-ray luminosity (LX,WPC) as a function of the orbital phase. The solid and dashed line has the same meaning as in Fig. 9.

thumbnail Fig. 11

Plot of the expected normal pre-shock velocity as a function of the angular separation from the line of centers (θ) at periastron (left) and apastron (right), in the framework of the WPC scenario. Solid and dashed lines stand for the non-inhibited and strongly radiatively inhibited cases, respectively.

thumbnail Fig. 12

Velocity curves (top panels) and post-shock temperatures (bottom panels) assuming two different values for the β parameter for both winds: 0.8 (solid curves) and 1.0 (dashed curves). The primary star is set at D = 0, and the secondary is at the other end of the horizontal axis. The vertical dotted lines delimit the maximum region available for wind acceleration between the stellar bodies along the line of centers.

6.4 Wind-wind collision

The case of WWC has been extensively discussed by seminal papers in this field (see e.g., Stevens et al. 1992; Pittard & Stevens 2002; Pittard & Parkin 2010). The occurrence of WWC is expected to significantly affect the X-ray properties of the system as compared to the WPC scenario. The main reason lies in the fact that such a WWC would occur in the acceleration region of the primary wind. As such, a WWC would occur further away from the secondary as compared to the shock discussed in Sect. 6.3; the primary wind pre-shock velocity would be lower, leading to significantly lower post-shock temperatures. This is illustrated in Fig. 12 where the velocity curves are plotted for both winds. Besides, a WWC would harbor two shocks, including one on the secondary’s side. Given the proximity with the secondary star, this secondary shock would be characterized by low pre-shock velocities, resulting in a quite soft X-ray emission component from the post-shock secondary plasma.

Depending on the position of the WWC between the two stars, different ranges of pre-shock velocities could be covered. It is also noticeable the the small size of the orbit and its quite large eccentricity lead the distance between photospheres to range between ~14.7 R at periastron and ~58.2 R at apastron (i.e., a factor ~4). This results in substantial differences in pre-shock velocities as a function of the orbital phase. If strong radiative inhibition was considered, curves in the upper panels of Fig. 12 would be significantly dropped, resulting in a severe drop of the post-shock temperature curves (lower panels). Additionally, the velocity curves plotted in Fig. 12 are valid along the line of centers. Accounting for obliquity due to the WWC warped about the secondary, one may expect the average X-ray emitting plasma temperature on the primary side to be dropped by about a factor 2. The obliquity is however expected to be less severe on the secondary’s side, resulting in an obliquity reduction factor still lower than 1, but significantly larger than 0.5. This would overall constitute a way to predict plasma temperature closer to those measured from the X-ray emission in Sect. 3.

thumbnail Fig. 13

Illustration of the orbit (dashed line) of HD 93205 in the orbital plane with the primary set at coordinates (x = 0, y = 0) (left) and as a projected orbit assuming two likely inclination angles (right). Solid circles represent the primary and secondary stellar bodies. Stellar bodies and orbit sizes are to scale. Orbital phase increases counterclockwise (resp. clockwise) along the orbit in top right (resp. bottom right) illustration.

6.5 Hybrid wind collision

Depending on adopted values for stellar wind parameters and on the action of efficient radiative effects, a stable stagnation point between the two wind may exist at some orbital phases, especially away from periastron. Figure 7 (along with the discussion in Sect. 6.2) also tells us that radiative braking may prevent the existence of WPC around apastron, over a part of the orbit that depends sensitively on adopted parameters. This scenario would therefore consider that the nature of the wind interaction would change over timescales shorter than the orbit. It is however very difficult to predict how such a wind collision system would relax from such a change of configuration. One can just anticipate a rather complex situation whose accurate description is far beyond the scope of this study. On the one hand, this would require a detailed time-dependent modeling of the hydrodynamic of the wind interaction, and on the other hand our current data are not sufficient to provide any insight into this topic.

7 A potential scenario for HD 93205

7.1 Orbital modulation

According to the results reported in Sect. 4, the X-ray emission presents a clear orbital modulation. The maximum-to-minimum LX,excess for the broad band (see Sect. 5) is about 1.6. This means that the X-ray emission from the wind interaction region drops by a factor slightly greater than 1/3 from maximum to minimum. This phase-locked X-ray variability is in principle attributable to a variation of the emission (due to the significant eccentricity) modulated by orientation effects.

The intrinsic X-ray emission from the wind interaction region is expected to peak close to periastron, when the emission measure of the emitting plasma is at maximum. In the framework of the WPC scenario, the expected minimum occurs indeed at apastron (Fig. 10). Assuming a WWC scenario, one would expect a similar behavior (see e.g., Pittard & Parkin 2010). Even though our light curves display a peak close to periastron, we also notice a low state spanning orbital phase between about 0.5 and 1.0, with a steep rise at periastron followed by a shallower decrease of the X-ray flux. This asymmetric behavior certainly points to a modulation by orientation effects.

The discussion of orientation effects in HD 93205 requires having a close look at the orbit. Using the orbital elements published by Morrell et al. (2001), we plotted the orbit of the system in Fig. 13. As the longitude of the ascending node is not known for HD 93205, we set it to 0°. This however doesn’t affect our discussion as this parameter only rules the rotation of the projected orbit on the sky plane. The left part of the figure shows the system in the orbital plane. The right part illustrates the projected orbit making two different assumptions on the actual inclination, as the available information doesn’t allow to lift the ambiguity between i = 58° and i′ = 180° − i = 122°. For both assumed inclination angles, the secondary is in the background at periastron and in the foreground at apastron.

One has to consider that the secondary stellar body is likely to hide a part of the X-ray emission region depending on the orbital phase, therefore leading to a drop of the X-ray flux. In addition, the X-ray emission is supposed to travel through an orbital phase dependent absorption column. The measured X-ray maximum occurs right after periastron, while the secondary’s surface exposed to the primary wind is the most visible. The inclination prevents the secondary’s surface to be eclipsed by the primary, even though one may expect some significant absorption by the primary wind. At minimum, after apastron, a significant part of the shocked region is occulted by the secondary’s body. Given the results reported on here, it seems that the impact of the occultation is more severe than that of the absorption by the primary wind material. This suggests that, even though one could expect the X-ray emitting region to be somewhat extended, a significant part of the X-ray emitting wind interaction is confined in a narrow region close to the secondary’s surface.

thumbnail Fig. 14

Plot of the cooling parameter (ζcooling) as a function of the off-axis angle (θ) for periastron (green) and apastron (blue), assuming the WPC scenario. The solid and dashed lines have the same meaning as in Fig. 9.

7.2 Shock types

Whatever the nature of the interaction (WPC or WWC), the nature of the shock (radiative vs adiabiatic) is also important to understand the X-ray emission from HD 93205. We evaluated the cooling parameter (ζcooling) as defined by Stevens et al. (1992). Assuming the WPC scenario, we calculated ζcooling =D12V84M˙7${\zeta _{cooling{\rm{ }}}} = {{{D_{12}}V_8^4} \over {{{\dot M}_{ - 7}}}}$(9)

where D12 is the separation from the primary to a specific point of the shock front (P, in Fig. A.1) in units of 1012 cm, V8 is the pre-shock velocity normal to the front at the same point in units of 108cm s−1, and −7 is the mass loss rate of the primary in units of 10−7 M yr−1. A ζcooling significantly above 1.0 points to an adiabatic regime, while a value close to or lower than 1.0 it typical of a radiative regime. Figure 14 shows ζcooling as a function of the off-axis angles (θ). At apastron, ζcooling is greater than 1.0 for θ values lower than 7.5° (resp. 4.0°) in the non-inhibited (resp. severely radiatively inhibited) case. Sufficiently away from the line of centers, shock obliquity reduces sufficiently the pre-shock velocity to significantly affect ζcooling. At periastron, ζcooling is on the order of (or even lower than) 1.0 over the whole shock front, even in the un-inhibited case.

The WWC scenario involves shocks occurring a bit further away from the secondary’s photosphere. The pre-shock velocity on the primary’s side is thus lower than in the case of the WPC scenario, and on the secondary’s side a much lower pre-shock velocity is expected given the short available distance for wind acceleration (see Fig. 12). For both shocks, ζcooling is thus very low, given its strong dependence on the velocity. It is thus very reasonable to consider that both shocks will be radiative at any orbital phase, and for most (if not all) off-axis angles. The only exception could happen close to apastron, near the line of centers.

Finally, in the framework of HWC scenario, radiative braking could at best potentially prevent a photospheric collision to happen near apastron. Apastron would thus behave as a WWC system (where shocks are more likely radiative), and periastron would be more compliant with the WPC scenario (with radiative shocks). The main result is that the HWC scenario would very likely result in a fully radiative regime along the full orbit, whatever the position on the shock front.

The prevalence of the radiative nature of the shocks has inevitably consequences on the X-ray emission. The occurrence of the thin-shell instability is expected to alter the shock front, leading to an irregularly corrugated structure (Vishniac 1994; Walder & Folini 1998). The resulting obliquity of the shocks causes a drastic drop of the pre-shock velocity, leading to significantly lower post-shock temperatures (Kee et al. 2014). This could explain the quite low post-shock temperature, that it, about 0.2 and 0.6 keV, measured across the orbit. On top of the likely radiative inhibition (see Sect. 6.2), this additional obliquity should severely affect the energy injection of the primary wind into the shock front, therefore reducing substantially the heating of the X-ray emitting plasma.

Coming again to the WPC scenario that should be valid along at least a major part of the orbit, Fig. 10 shows a predicted X-ray luminosity ranging between 6 × 1031 (at apastron) and 1033 erg s−1 (at periastron), in the radiatively inhibited case. If we account for the significant drop of the pre-shock velocity due to the thin-shell instability, one has to consider a much lower average post-shock plasma temperature. Typically, at the maximum plasma temperature we measured in Sect. 3, that is, ~0.6 keV, we have to multiply that predicted X-ray luminosity by a factor of about 5 to account for significant line emission on top of the continuum (see Fig. 8 and discussion in Sect. 6.3). This leads to an X-ray luminosity ranging between 3 × 1032 (at apastron) and 5 × 1033 erg s−1 (at periastron). The former value is compliant with our measurements (Sect. 5), while the latter appears to be quite excessive. The lack of significant rise of the X-ray luminosity at peristron could be attributable to a significant photoelectic absorption by the primary wind, but it may also result from the radiative nature of the shocks, more pronounced close to periastron. According to 2D simulations of collision between equal expanding winds by Kee et al. (2014), such radiative shocks should be less efficient at producing X-rays, with a reduction of the X-ray luminosity by a factor up to a few tens in the most severe cases (see e.g., De Becker et al. 2004; Antokhin et al. 2004). The suppression of thermal X-rays in radiative shocks was also addressed by Steinberg & Metzger (2018). More recent hydrodynamic simulations of Kashi (2020) for colliding wind structure of the massive binary system HD 166734 (07.5If + 09I(f), period ≃ 34.538 d, eccentricity ≃ 0.618) along its binary orbit have also suggested formation of instabilities in the nonadiabatic interaction of unequal winds. The resulting X-ray emission is seen to be reduced as a consequence of accretion of primary wind onto secondary close to the periastron.

8 Summary and conclusions

We investigated the X-ray emission from the eccentric massive binary HD 93205 on the basis of a time series of XMM-Newton observations spanning about two decades. The X-ray spectrum is quite soft, with a typical plasma temperature not higher than 0.6 keV. We identified a periodic variability of the X-ray flux with a period of about 6.08 d, in full agreement with the orbital solution published by Morrell et al. (2001).

The stellar wind of the primary component (O3.5 V((ff))) clearly dominates that of the secondary (O8 V), inhibiting the formation of a regular wind-wind interaction in between the two stars. Taking into account radiative effects such as sudden radiative breaking, radiative inhibition, and the potential role played by self-regulated shocks, we conclude that HD 93205 very likely undergoes a wind-photosphere interaction along most of its orbit. We caution, however, that this statement is supported by the one-dimensional approach used to address this topic, and we cannot reject the idea of a potential revision of this conclusion upon a more detailed three-dimensional analysis, which is out of the scope of the present paper. The only exception is close to apastron, where one cannot exclude a wind-wind interaction to occur. One could thus be dealing with a hybrid wind collision, switching between a wind-photosphere collision and a wind-wind collision depending on the orbital phase. The orbital modulation of the X-ray emission is interpreted in terms of the variation of the intrinsic emission modulated by orientation effects (photoelectric absorption and occultation of a part of the X-ray emitting region by a secondary stellar body). The X-ray light curve peaks close to periastron, with a steep rise at peri-astron followed by a shallower decay. The low state of the X-ray emission occurs between phases 0.5 and 1.0, while the secondary is in front. The latter configuration implies that a significant part of the X-ray emitting interaction region is located on the opposite side of the secondary stellar body (close to the line of centers, between the two stars). This is in agreement with our prediction of an interaction region (WPC of WWC) very close to the secondary’s surface.

Despite the short separation in the system allowing at least the primary wind to hit the secondary’s photosphere before it reaches its terminal velocity, predicted pre-shock velocities are too high to be compliant with the very soft emission revealed by the X-ray monitoring. Even considering a quite severe radiative inhibition dropping the pre-shock velocity by a factor 2, the predicted X-ray emission is harder than measured. This discrepancy could be explained by the radiative nature of the shock. In such shocks, instabilities significantly alter the shape of the shock front, leading to some additional obliquity in the incidence of the wind flow onto the shock surface.

Our study clearly demonstrates that HD 93205 is a rather challenging laboratory to investigate shock physics in short period, asymmetric binary systems, where various radiative effects and radiative shocks concur to display an interesting behavior in soft X-rays. We encourage readers to consider this target as a test case for state-of-the-art future modeling including, self-consistently, all the physics that we considered in our discussion.

Acknowledgements

We thank the referee for the careful reading of the manuscript and giving us constructive comments and suggestions. BA acknowledges the grant provided by Wallonia Brussels International to carry out this work. This research has made use of observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and NASA.

Appendix A Shock shape in the framework of the WPC scenario

As a first approximation of the shock shape, we made use of the shock stand-off distance and shape equation for a shock produced by a flow hitting a sphere (Billig 1967; Cardona & Lago 2023). In cartesian coordinates with the x-axis along the line of centers of the binary, the equation expressing y as a function of x is the following, y=Rctanζ[ (xRSΔRccotan2ζ+1)21 ]1/2$y = {{{R_c}} \over {\tan \zeta }}{\left[ {{{\left( {{{x - {R_S} - \Delta } \over {{R_c}{{{\mathop{\rm cotan}\nolimits} }^2}\zeta }} + 1} \right)}^2} - 1} \right]^{1/2}}$(A.1)

where RS is the secondary stellar radius, Rc is the shock curvature radius, Δ is the stand-off distance (separation between the shock apex and the stellar surface along the line of centers), and ζ is defined as ζ=arcsin (1Ma)$\zeta = \arcsin \left( {{1 \over {{M_a}}}} \right)$(A.2)

where Ma is the sonic pre-shock Mach number of the primary wind.

The expressions for Rc and Δ are the following: Rc=1.143RS exp (0.54(Ma1)1.2)${R_c} = 1.143{R_S}\exp \left( {{{0.54} \over {{{\left( {{M_a} - 1} \right)}^{1.2}}}}} \right)$(A.3) Δ=0.143RS exp (3.2Ma2)$\Delta = 0.143{R_S}\exp \left( {{{3.2} \over {M_a^2}}} \right)$(A.4)

Assuming a wind temperature that is 0.5 times the effective temperature, the sound speed (expressed as vs=0.1Tw,P${\v _s} = 0.1\sqrt {{T_{w,P}}} $, assuming a perfect gas with mean molecular weight on the order of 1.0) is on the order of 15 km s−1. Given the expected range of pre-shock velocities as a function of orbital phase, Ma ranges between 170 and 220. In that hypersonic regime, the shape and position of the shock are quite insensitive to changes of pre-shock sonic Mach number. Along the full orbit, we obtain Rc = 9.75 R and Δ = 1.22 R. We note that these values are still valid in case of severe radiative inhibition leading to a drop by a factor 2 of the primary wind velocity (as considered in Fig. 9 and Fig. 10).

thumbnail Fig. A.1

Illustration of the shock shape (red curve) produced by the WPC at an arbitrary orbital phase. Details required to compute the normal component of the pre-shock velocity are shown. The line of centers of the system is along the x-axis.

An illustration of the shock shape is presented in Fig. A.1. The pre-shock velocity vector is represented for a specific flow line coming from the primary and hitting the shock at point P. The distance from the primary star to point P is easily expressed from the geometry of the system, s=(DX)2+y2$s = \sqrt {{{(D - X)}^2} + {y^2}} $(A.5)

Knowing s, the pre-shock velocity VP was computed using a standard β velocity law (Eq. 1).

We made use of the geometry shown in Fig. A.1 (based on Eq. A.1) to determine for every point of the shock curve the amplitude of the normal component of the pre-shock velocity (VP,n) that counts for the shock physics at work in the system. With angles θ and ψ defined as illustrated in Fig. A.1, we note that ϕ=π/2ψ+θ$\phi = \pi /2 - \psi + \theta $(A.6)

and VP,n=VPcosϕ${V_{P,n}} = {V_P}\cos \phi $(A.7)

The angle θ is determined for any point P from the x and y coordinates related through Eq. A.1, θ= arctan (yDx)$\theta = \arctan \left( {{y \over {D - x}}} \right)$(A.8)

Table A.1

Log of XMM-Newton observations of HD 93205.

Table A.2

Best-fit parameters obtained from spectral fitting of HD 93205 as observed from XMM-Newton-EPIC.

References

  1. Abbott, D. C. 1982, ApJ, 259, 282 [Google Scholar]
  2. Albacete-Colombo, J. F., Damiani, F., Micela, G., Sciortino, S., & Harnden, F. R., J. 2008, A&A, 490, 1055 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Albacete Colombo, J. F., Méndez, M., & Morrell, N. I. 2003, MNRAS, 346, 704 [NASA ADS] [CrossRef] [Google Scholar]
  4. Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 [Google Scholar]
  5. Antokhin, I. I., Rauw, G., Vreux, J. M., & van der Hucht, K. A. 2003, in Astronomical Society of the Pacific Conference Series, Vol. 305, Magnetic Fields in O, B and A Stars: Origin and Connection to Pulsation, Rotation and Mass Loss, eds. L. A. Balona, H. F. Henrichs, & R. Medupe, 383 [NASA ADS] [Google Scholar]
  6. Antokhin, I. I., Owocki, S. P., & Brown, J. C. 2004, ApJ, 611, 434 [Google Scholar]
  7. Antokhin, I. I., Rauw, G., Vreux, J. M., van der Hucht, K. A., & Brown, J. C. 2008, A&A, 477, 593 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Arora, B., Pandey, J. C., & De Becker, M. 2019, MNRAS, 487, 2624 [Google Scholar]
  9. Arora, B., Pandey, J. C., De Becker, M., et al. 2021, AJ, 162, 257 [NASA ADS] [CrossRef] [Google Scholar]
  10. Arora, B., De Becker, M., & Pandey, J. C. 2024, A&A, 687, A34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Demleitner, M., & Andrae, R. 2021, AJ, 161, 147 [Google Scholar]
  12. Billig, F. S. 1967, J. Spacecr. Rockets, 4, 822 [NASA ADS] [CrossRef] [Google Scholar]
  13. Broos, P. S., Townsley, L. K., Feigelson, E. D., et al. 2011, ApJS, 194, 2 [NASA ADS] [CrossRef] [Google Scholar]
  14. Cardona, V., & Lago, V. 2023, Phys. Lett. A, 491, 129185 [NASA ADS] [CrossRef] [Google Scholar]
  15. Castor, J. I., Abbott, D. C., & Klein, R. I. 1975, ApJ, 195, 157 [Google Scholar]
  16. Chlebowski, T., & Garmany, C. D. 1991, ApJ, 368, 241 [NASA ADS] [CrossRef] [Google Scholar]
  17. Conti, P. S., & Walborn, N. R. 1976, ApJ, 207, 502 [NASA ADS] [CrossRef] [Google Scholar]
  18. Corcoran, M. F. 1996, in Revista Mexicana de Astronomia y Astrofísica Conference Series, 5, eds. V. Niemela, N. Morrell, P. Pismis, & S. Torres-Peimbert, 54 [Google Scholar]
  19. Corcoran, M. F., Swank, J., Rawley, G., et al. 1995, in Revista Mexicana de Astronomia y Astrofísica Conference Series, 2, eds. V. Niemela, N. Morrell, & A. Feinstein, 97 [Google Scholar]
  20. De Becker, M. 2013, New A, 25, 7 [CrossRef] [Google Scholar]
  21. De Becker, M. 2015, MNRAS, 451, 1070 [Google Scholar]
  22. De Becker, M., Rauw, G., Pittard, J. M., et al. 2004, A&A, 416, 221 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Gagné, M., Fehon, G., Savoy, M. R., et al. 2011, ApJS, 194, 5 [CrossRef] [Google Scholar]
  24. Gayley, K. G., Owocki, S. P., & Cranmer, S. R. 1997, ApJ, 475, 786 [Google Scholar]
  25. Horne, J. H., & Baliunas, S. L. 1986, ApJ, 302, 757 [Google Scholar]
  26. Jansen, F., Lumb, D., Altieri, B., et al. 2001, A&A, 365, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Jenkins, E. B. 2019, ApJ, 872, 55 [Google Scholar]
  28. Kashi, A. 2020, MNRAS, 492, 5261 [NASA ADS] [CrossRef] [Google Scholar]
  29. Kee, N. D., Owocki, S., & ud-Doula, A. 2014, MNRAS, 438, 3557 [Google Scholar]
  30. Lomb, N. R. 1976, Ap&SS, 39, 447 [Google Scholar]
  31. Martins, F., Schaerer, D., & Hillier, D. J. 2005, A&A, 436, 1049 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Morrell, N. I., Barbá, R. H., Niemela, V. S., et al. 2001, MNRAS, 326, 85 [NASA ADS] [CrossRef] [Google Scholar]
  33. Muijres, L. E., Vink, J. S., de Koter, A., Müller, P. E., & Langer, N. 2012, A&A, 537, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Owocki, S. P., Sundqvist, J. O., Cohen, D. H., & Gayley, K. G. 2013, MNRAS, 429, 3379 [NASA ADS] [CrossRef] [Google Scholar]
  35. Parkin, E. R., & Sim, S. A. 2013, ApJ, 767, 114 [NASA ADS] [CrossRef] [Google Scholar]
  36. Pittard, J. M., & Stevens, I. R. 1997, MNRAS, 292, 298 [NASA ADS] [Google Scholar]
  37. Pittard, J. M., & Stevens, I. R. 2002, A&A, 388, L20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Pittard, J. M., & Parkin, E. R. 2010, MNRAS, 403, 1657 [NASA ADS] [CrossRef] [Google Scholar]
  39. Roberts, D. H., Lehar, J., & Dreher, J. W. 1987, AJ, 93, 968 [Google Scholar]
  40. Ryspaeva, E., & Kholtygin, A. 2020, Res. Astron. Astrophys., 20, 108 [CrossRef] [Google Scholar]
  41. Scargle, J. D. 1982, ApJ, 263, 835 [Google Scholar]
  42. Seward, F. D., & Chlebowski, T. 1982, ApJ, 256, 530 [NASA ADS] [CrossRef] [Google Scholar]
  43. Seward, F. D., Forman, W. R., Giacconi, R., et al. 1979, ApJ, 234, L55 [NASA ADS] [CrossRef] [Google Scholar]
  44. Sota, A., Maíz Apellániz, J., Morrell, N. I., et al. 2014, ApJS, 211, 10 [Google Scholar]
  45. Steinberg, E., & Metzger, B. D. 2018, MNRAS, 479, 687 [NASA ADS] [Google Scholar]
  46. Stevens, I. R., & Pollock, A. M. T. 1994, MNRAS, 269, 226 [Google Scholar]
  47. Stevens, I. R., Blondin, J. M., & Pollock, A. M. T. 1992, ApJ, 386, 265 [Google Scholar]
  48. Usov, V. V. 1992, ApJ, 389, 635 [Google Scholar]
  49. Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2001, A&A, 369, 574 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Vishniac, E. T. 1994, ApJ, 428, 186 [Google Scholar]
  51. Walborn, N. R. 1971, ApJ, 167, L31 [CrossRef] [Google Scholar]
  52. Walborn, N. R. 1973, ApJ, 179, 517 [NASA ADS] [CrossRef] [Google Scholar]
  53. Walder, R., & Folini, D. 1998, A&A, 330, L21 [NASA ADS] [Google Scholar]

All Tables

Table 1

Main orbital parameters of HD 93205 as given by Morrell et al. (2001).

Table 2

Adopted parameters for both components of HD 93205, based on Martins et al. (2005) and Muijres et al. (2012).

Table A.1

Log of XMM-Newton observations of HD 93205.

Table A.2

Best-fit parameters obtained from spectral fitting of HD 93205 as observed from XMM-Newton-EPIC.

All Figures

thumbnail Fig. 1

X-ray image of the Eta Carinae region extracted from observation ID 0311990101 in the 0.3–12 keV energy range. The position of HD 93205 as been highlighted with a white circle of 25 arcsec radius. The X-ray flux of neighboring source HD 93204 has been removed from a circular region of 8 arcsec radius. An estimate of background was done from a circular region (cyan color) having radius 15 arcsec lying in the source-free area of the detector.

In the text
thumbnail Fig. 2

MOS1 and MOS2 spectra of HD 93205 jointly fitted with the two-temperature thermal plasma emission model from the observation ID 0145740301. The details of the spectral fitting are provided in Table A.2.

In the text
thumbnail Fig. 3

Variation of the local equivalent H-column density (NHlocal)$\left( {{\rm{N}}_H^{local}} \right)$, normalization constants (norm1 and norm2), and the ratio normı/norm2 corresponding to two thermal plasma emission components with the orbital phase of HD 93205 (see Table A.2).

In the text
thumbnail Fig. 4

Modulation of the (a) observed (Fobs, left) and (b) ISM-corrected (Fism, right) X-ray flux in broad, soft, and hard energy bands from HD 93205 obtained after X-ray spectral fitting as a function of the binary orbital phase (see Table A.2).

In the text
thumbnail Fig. 5

Lomb–Scargle and CLEANed power spectra of observed flux (FBobs)$\left( {F_B^{obs}} \right)$ from HD 93205 in the 0.3-5.0 keV energy range. The frequency of the peak with the highest power is also mentioned.

In the text
thumbnail Fig. 6

Folded X-ray light curves of HD 93205 as observed by XMM-Newton using (a) MOS1 (left) and (b) MOS2 (right) in broad, soft, and hard energy bands. The average count rate of each individual light curve observed with thick optical blocking filter has been plotted here.

In the text
thumbnail Fig. 7

Plot of the braking parameter (Fbr) as a function of the orbital phase. The solid and dotted line stand for two different assumptions for the k parameter. The horizontal dashed line stands for the critical value above which the WPC scenario is not valid.

In the text
thumbnail Fig. 8

Plot of the total-to-continuum X-ray flux ratio as a function of plasma temperature.

In the text
thumbnail Fig. 9

Plot of the stellar separation (upper panel), the pre-shock velocity (middle panel), and the post-shock temperature (lower panel) as a function of the orbital phase. Solid lines stand for a first approximation of the pre-shock velocity while dashed lines illustrate the case of strong radiative inhibition.

In the text
thumbnail Fig. 10

Plot of the predicted X-ray luminosity (LX,WPC) as a function of the orbital phase. The solid and dashed line has the same meaning as in Fig. 9.

In the text
thumbnail Fig. 11

Plot of the expected normal pre-shock velocity as a function of the angular separation from the line of centers (θ) at periastron (left) and apastron (right), in the framework of the WPC scenario. Solid and dashed lines stand for the non-inhibited and strongly radiatively inhibited cases, respectively.

In the text
thumbnail Fig. 12

Velocity curves (top panels) and post-shock temperatures (bottom panels) assuming two different values for the β parameter for both winds: 0.8 (solid curves) and 1.0 (dashed curves). The primary star is set at D = 0, and the secondary is at the other end of the horizontal axis. The vertical dotted lines delimit the maximum region available for wind acceleration between the stellar bodies along the line of centers.

In the text
thumbnail Fig. 13

Illustration of the orbit (dashed line) of HD 93205 in the orbital plane with the primary set at coordinates (x = 0, y = 0) (left) and as a projected orbit assuming two likely inclination angles (right). Solid circles represent the primary and secondary stellar bodies. Stellar bodies and orbit sizes are to scale. Orbital phase increases counterclockwise (resp. clockwise) along the orbit in top right (resp. bottom right) illustration.

In the text
thumbnail Fig. 14

Plot of the cooling parameter (ζcooling) as a function of the off-axis angle (θ) for periastron (green) and apastron (blue), assuming the WPC scenario. The solid and dashed lines have the same meaning as in Fig. 9.

In the text
thumbnail Fig. A.1

Illustration of the shock shape (red curve) produced by the WPC at an arbitrary orbital phase. Details required to compute the normal component of the pre-shock velocity are shown. The line of centers of the system is along the x-axis.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.