Accreting Black Holes Skewing and Bending the Optical Emission from Massive Wolf-Rayet Companions -- A Case Study of IC10 X-1

Sayantan Bhattacharya,Dimitris M. Christodoulou,Andre-Nicolas Chene,Silas G. T. Laycock,Breanna A. Binder,Demosthenes Kazanas

DOI: https://doi.org/10.1093/mnras/stad2094

2023-07-07

Abstract:We present a statistical analysis of the He ii 4686 emission line in the spectra of the black hole and Wolf-Rayet (WR) star of the high-mass X-ray binary IC10 X-1. This line is visibly skewed, and the third moment (skewness) varies with the binary's orbital phase. We describe a new method of extracting such weak/faint features lying barely above a noisy continuum. Using the moments of these features, we have been able to decompose these skewed lines into two symmetric Gaussian profiles as a function of the orbital phase. The astrophysical implications of this decomposition are significant due to the complex nature of wind-accretion stream interactions in such binary systems. Previous studies have already shown a 0.25 phase lag in the radial velocity curve of the star and the X-ray eclipse, which indicates that the He ii emitters might be in the stellar wind, hence not tracing the star's orbital motion. Results from this work further suggest the existence of two separate emitting regions, one in the stellar wind in the shadow of the WR star, and another in the accretion stream that impacts the black hole's outer accretion disk; and the observed skewed He ii lines can be reproduced by superposition of the two corresponding time-dependent Gaussian emission profiles.

High Energy Astrophysical Phenomena,Solar and Stellar Astrophysics

What problem does this paper attempt to address?

### Problems the paper attempts to solve

This paper aims to study the wind - accretion flow interaction between black holes (BH) and Wolf - Rayet stars (WR stars) by analyzing the He II 4686 emission line of high - energy X - ray binaries (HMXB) in the IC10 X - 1 system. Specifically, the paper attempts to solve the following problems:

1. **Skewness of the He II 4686 emission line**:
- This emission line shows obvious skewness in the spectrum, and its skewness (i.e., the third - order moment) changes with the binary orbital phase. The paper proposes a new method for extracting these weak features from the noisy background and decomposing them into two symmetric Gaussian distributions.

2. **Physical origin of the emission line**:
- By decomposing the skewed He II 4686 emission line, the paper attempts to reveal the physical sources of these emission lines. Research shows that there may be two independent emission regions: one is located in the shadow region of the WR star wind, and the other is located in the accretion flow that impacts the accretion disk outside the black hole.

3. **Dynamics of the binary system**:
- Previous studies have found that there is a 0.25 - phase lag between the radial velocity curve and the X - ray eclipse, indicating that the He II emission line may come from the stellar wind rather than directly tracing the orbital motion of the star. The paper further explores this phenomenon and proposes a physical model of the two emission regions.

4. **Statistical analysis method**:
- The paper introduces a new data - fitting method for processing weak signals in the noisy background. This method can more accurately extract and analyze the statistical moments of the He II 4686 emission line, thereby better understanding the dynamic processes of the binary system.

### Main conclusions

- **Emission line decomposition**: By decomposing the skewed He II 4686 emission line into two symmetric Gaussian distributions, researchers can more clearly identify two different emission regions.
- **Physical significance**: These two emission regions correspond to the shadow region of the WR star wind and the accretion flow that impacts the accretion disk outside the black hole respectively, which provides important clues for understanding the complex interaction of the wind - accretion flow.
- **Dynamic model**: The model proposed in the paper helps to explain the observed skewness of the He II 4686 emission line and its change with the orbital phase, providing a theoretical basis for further research on the dynamics of the binary system.

### Formula display

- **Decomposition of Gaussian distribution**:
\[
0=\delta_1 + \delta_2
\]
\[
2V = V_1+\delta_2^2+V_2+\delta_2^2
\]
\[
2S=\delta_1(3V_1 + \delta_1^2)+\delta_2(3V_2+\delta_2^2)
\]
\[
2K = 3V_1^2+6V_1\delta_1^2+\delta_1^4+3V_2^2+6V_2\delta_2^2+\delta_2^4
\]
where \(\delta_i=\mu_i-\mu\).

- **Solution of cubic equation**:
\[
6x^3+3(K - 3)x - S^2=0
\]
\[
\mu_{\pm}=\mu\pm c
\]
\[
c = \sqrt{Vx}=D\sqrt{x}
\]
\[
V_{\pm}=V\left(1 - x\pm\frac{S}{3\sqrt{x}}\right)
\]

- **Separation distance and bimodal index**:
\[
D_{\mu}=\frac{\Delta\mu}{2D}=\frac{c}{D}=\sqrt{x}
\]
\[
D_V=\frac{V_1 - V_2}{2V}=\frac{S}{3\sqrt{x}}=\frac{S}{3D_{\mu}}
\]