-
PDF
- Split View
-
Views
-
Cite
Cite
Janusz Ziółkowski, Andrzej A Zdziarski, A comment on the properties of the matter flow through the first Lagrangian point, Monthly Notices of the Royal Astronomical Society, Volume 508, Issue 3, December 2021, Pages 3495–3498, https://doi.org/10.1093/mnras/stab2826
- Share Icon Share
ABSTRACT
We analyse properties of the mass outflow from the Roche lobe filling component of a semidetached binary system. We follow the approaches published by Paczyński & Sienkiewicz and by Lubow & Shu, which we compare with other simplified approaches. We find that the density of the flow at L1 is orders of magnitude lower than the density on the same equipotential but away from L1. Furthermore, the effective cross-section of the flow, after averaging over its profile of the momentum density, is much lower than some published estimates done without accounting for the averaging. Thus, the use of some simplified formulae for the density and the flow cross-section can lead to overestimates of the accretion rate and of the mass contained in the L1 regions by very large factors unless they are supported by simultaneous integrations of the equations of stellar structure for the outer layers of the donor.
1 INTRODUCTION
The knowledge of the mass transfer rate from the donor to accretor, |$\dot{M}$|, in accreting stellar systems is of major importance for their understanding. By comparing the transfer rate with the accretion rate, we can assess whether the mass transfer is conservative or associated with mass-loss, as well as examine an effect of irradiation on the mass flow. The transfer rate then determines the past and future evolution of accreting systems. Observationally, |$\beta \dot{M}$| can be estimated from a long-term average of the luminosity in an X-ray binary, e.g. Coriat, Fender & Dubus (2012) and Zdziarski, Ziółkowski & Mikołajewska (2019), where β ≤ 1 is the fraction of the transferred mass accreted. However, for understanding of many aspects of these systems, it is important to estimate |$\dot{M}$| theoretically, and compare the two estimates. In semidetached binaries, such estimates can be done based on considering the evolutionary changes of the stellar radius compared to the changes of the Roche lobe, e.g. Webbink, Rappaport & Savonije (1983). If the mass transfer is driven by nuclear expansion of the donor, |$\dot{M}$| is connected to the time-scale of this expansion, see equation (4.16) of Frank, King & Raine (2002). Using both approaches, we can relatively reliably determine the rate of the mass transfer, |$\dot{M}$|.
On the other hand, sometimes it is useful to consider the dynamics of the flow through the vicinity of the inner Lagrangian point of the Roche lobe, L1. This is the case if we need to consider the conditions near the L1 point, e.g. in order to investigate the effects of irradiation. Such analysis is helped by approximate formulae relating the radius excess, ΔR, i.e. the difference between the stellar radius and the radius of the Roche lobe, to the rate of the mass transfer |$\dot{M}$| (Paczyński & Sienkiewicz 1972, hereafter PS72, who used in part unpublished results of Jȩdrzejec 1969; Savonije 1978, 1979, 1983). On the other hand, Lubow & Shu (1975), hereafter LS75, related the physical conditions at L1 to the sound speed and the binary parameters only. Their analysis has been widely applied in several contexts, see, e.g. Lai, Helling & van den Heuvel (2010).
The accuracy of these formulae remains an open question. Does the rate of the mass transfer can be obtained from either the radius excess or the sound speed? We would like to stress that the mechanism responsible for the determination of |$\dot{M}$| is the internal evolutionary engine of the donor. The structure of the outer stellar layers (including the vicinity of L1) just adjusts to carry out whatever |$\dot{M}$| is set by the evolutionary engine. This includes both the radius excess and the gas density at L1. Here, we test the accuracy of the formulae of PS72 and LS75, and find their applicability is not universal.
2 ESTIMATES OF THE FLOW PARAMETERS
2.1 The approach of PS72
Here, we recall and extend some of the results of PS72 and compare them to other expressions in literature. We use the coordinate system in the corotating frame as shown by PS72 in their fig. A1. The x-axis connects the stellar centres, and it is given in units of the semimajor axis (equal to the separation between the stars, A, for a circular orbit). The orbital motion is neglected, and thus the system has an axial symmetry around x, and the y-axis is perpendicular to x in any direction. We will also neglect a small difference between the cross-sections at x = 0 of the star overflowing its Roche lobe, yf in the notation of PS72, and of the equipotential surface which the star fills, ys. This results in a factor of slightly less than unity by which the true |$\dot{M}$| should be multiplied (see table A1 in PS72).
We note that Savonije (1978) repeated the analysis of PS72 in a more precise way. In particular, he took into account the orbital motion of the matter in the binary system. His formulae are more complicated, but the general results are similar to those of PS72. In particular, he obtained the same functional form of the relations given here by equations (4) and (14); specifically, with the same power exponents.
Equation (14) can give us a reasonably accurate estimate of ΔR if we know |$\dot{M}$|. This formula takes into account the physical state of the outflowing matter through the polytrope constant, K, and the polytropic exponent, n. Thus, in order to determine the values of these parameters we need a model of the outer layers of the donor. Importantly, we definitely need these values at the L1 point rather than at the depth ΔR below the surface of the star far from L1 (as practiced by some simplified formulae users). Also, if we use directly equation (10) with estimates of Σ, we need the values of density and pressure at L1 rather than away from it at the depth ΔR.
2.2 The approach of LS75
3 COMPARISON WITH A MODEL OF GX 339–4
As an example, we consider a model of the binary system GX 339–4. This system accretes from an evolved low-mass star on to a black hole, and its binary period is P = 1.7587 d. Specifically, we consider a model slightly modified with respect to the evolutionary model D obtained by Zdziarski et al. (2019), and assume |$\dot{M} \approx 6.1 \times 10^{16}$| g s−1. This |$\dot{M}$| is equal to the value obtained from the evolutionary model by matching the rate of Roche lobe expansion to the stellar expansion. In that model, we have M = 1 M⊙, R = 2.83 R⊙, the accretor mass of MX = 8 M⊙, giving μ = 1/9, and the radius excess of ΔR ≈ 1.22 × 109 cm (a slightly corrected value, following from equation 14). The stellar density at the depth ΔR at the equipotential describing the Roche lobe away from L1 was found to be ρ0 = 3.1 × 10−6 g cm−3 and the temperature was T0 = 1.56 × 104 K. The local value of the polytropic exponent, n, estimated at that depth, is n = 6.43. The matter there is almost fully ionized, and X = 0.74 and Z = 0.014 was assumed, which implies the mean molecular weight of 0.598. Using equation (3) and assuming ideal gas, we obtain the entropy parameter of K ≈ 1.55 × 1013 (cgs), and the sound speed of cs0 ≈ 1.6 × 106 cm s−1. These values of K and n are then assumed to apply to the flow through L1 and be constant within it.
For the above parameters, we can find the velocity and density at L1 from equations (16) and (17), respectively. They are |$v_{\rm \mathit{ L}_1}\approx 7.8\times 10^5$| cm s−1 and |$\rho _{\rm \mathit{ L}_1}\approx 3.4\times 10^{-10}$| g cm−3. Thus, the density at L1 is as much as ≈104 times smaller than the stellar density at the depth ΔR (while the sound speed is by a factor of two smaller).
The geometrical cross-section of the flow can be calculated using equation (1). We get ys ≈ 0.027, which translates into the flow radius of 2.67 × 1010 cm and the total cross-section area of 2.24 × 1021 cm2. We note that ρ v decreases fast with increasing y, especially for a high value of n, which decrease is taken into account in our formulae for either 〈ρv〉 or Σeff, equations (8), (11), respectively. For the example considered here, ρv drops to 10−3 of its value at the point L1 at y ≈ 0.0215.
Now, let us compare the values of different parameters of the flow obtained with the help of different formulae used in literature.
We shall start with the geometrical cross-section of the flow. As was shown above, the analysis of PS72 leads to the radius of the flow equal to 2.67 × 1010 cm. Using estimate of LS75 given by equation (18), we get rf ≈ 1.89 × 1010 cm, i.e. a value smaller by a factor of only ≈1.4. Taking into account the approximate character of both estimates, the agreement is excellent. We have to remember, however, that using equation (18) requires the knowledge of the sound speed at the L1 point. In our case, this knowledge was based on the analysis following the approach of PS72 and on the simultaneous integration of the equations of stellar structure for the outer layers of the donor.
Now, let us compare the values of the flow effective (as opposed to the geometrical one) cross-section area obtained with the different formulae listed in Section 2. The value following from the analysis of PS72 is given by equation (11). For the considered example, Σeff ≈ 2.3 × 1020 cm2. Then, the approximation neglecting the averaging of the flow over the cross-section proposed by Zdziarski et al. (2007), of equation (13), yields Σ0 ≈ 1.5 × 1022 cm2, i.e. almost two orders of magnitude too much. The estimate by LS75, as given by our equation (20), leads to the value Σ0 ≈ 3.56 × 1020 cm2. This value is only by a factor ≈1.5 greater than the value obtained by us with equation (11), which again means that the agreement between the two approaches is excellent. Still, it requires the knowledge of the sound speed at L1.
We next compare the values of the density of the gas at L1 obtained with the different formulae listed in Section 2. The value following from the analysis of PS72 is given by equation (17), which in our case leads to |$\rho _{\rm \mathit{ L}_1}\approx 3.4\times 10^{-10}$| g cm−3. The estimate by LS75, as given by our equation (19), leads to |$\rho _{\rm \mathit{ L}_1}\approx 2.20\times 10^{-10}$| g cm−3. This is by a factor of only ≈1.5 smaller than the value obtained by us with equation (17). This is the same factor 1.5 as found by us while comparing the different estimates of the effective cross-section, which identity is a consequence of the equivalence of equations (19) and (20).
The value of the density of the gas at L1 is certainly one of the important parameters of the flow. The fact that the values of this parameter obtained with two different approaches agree so well is certainly encouraging. It also supports the claim that this estimate is not far from the true value.
Discussing the approximate formulae present in the literature, we should devote attention to popular formulae equivalent to our equation (10) but used in improper way, As an example we may recall equation (9) of Zdziarski et al. (2007). Those authors use equation (10) but they replace the effective cross-section area with the geometrical one and the density of the gas at L1 with the density on the Roche lobe equipotential away from L1. Since, as we have seen, the effective cross-section area might be overestimated this way by two orders of magnitude and the density of the gas even by four orders of magnitude, such use of equation (10) might lead to an error reaching even six orders of magnitude. It is therefore crucial, while using formulae of this type, to use the proper value of the gas density.
Finally, we should comment on the discrepancy between the density of the gas at the Roche equipotential inside the star far from L1 obtained in our calculations and the density obtained with the prescription of LS75, equation (21). Using it, we find the ratio of ≈52. As we saw in our analysis above, that ratio was about four orders of magnitude (3.1 × 10−6 versus 3.4 × 10−10 g cm−3). This is a serious discrepancy, by more than two orders of magnitude. However, we believe that our estimate is closer to the real situation than that of LS75. The reason is the following one. We started with the determination of the radius excess, ΔR. We used for that purpose equation (14). Fortunately, the sensitivity of ΔR to |$\dot{M}$| is very weak. As demonstrated by Zdziarski et al. (2019) (see their fig. 3), increasing |$\dot{M}$| by four orders of magnitude leads to increase of ΔR by a factor of only two. Since we have a relatively reliable estimate of |$\dot{M}$| based on evolutionary considerations, we may trust that ΔR is determined with rather good precision. Knowing the depth of the Roche equipotential below the stellar surface, we may integrate the stellar structure equations (including the heat transfer equation) from the stellar surface inward. At the depth ΔR, we determine, among others, the density of the gas and the gas sound speed. This procedure is rather straightforward and should not raise serious doubts. On the other hand, the LS75 estimate is based on rather qualitative analysis of a very complicated flow in the outer stellar layers.
Finally, we should stress that all the considered descriptions are still quite approximate. Hydrodynamic simulations are needed to get more accurate descriptions.
4 CONCLUSIONS
We have overviewed and extended the treatment of the flow through L1 of PS72, and compared it with other, simplified, methods. We have stressed that the flow density at L1 can be several orders of magnitude lower than that at the same equipotential but away from L1.
Furthermore, one of the approximate estimates of the flow effective cross-section used in literature overestimates the actual one by two orders of magnitude. The other one, given by LS75, leads to result that remains in excellent agreement with our estimate. The same is true about the estimate of the gas density at L1 given by LS75 (but to calculate it correctly we have to know the structure of the outer layers of the donor).
However, we have found that the estimate of the ratio of the density of the gas at the Roche equipotential inside the star far from L1 to that at L1 according to the prescription of LS75 is not accurate.
The overestimates of the density and the cross-section in some formulae then can lead to a gross overestimate of the outflow rate (calculated as a function of the parameters of the flow), as well as of the amount of mass present in the L1 flow. We should, however, remember that outflow rate is determined by the evolutionary engine of the donor.
Generally, we conclude that such simplified formulae can give reasonably accurate estimates if they are supported by simultaneous integrations of the equations of stellar structure for the outer layers of the donor.
ACKNOWLEDGEMENTS
We thank the referees for valuable suggestions. This research has been supported in part by the Polish National Science Centre grants 2015/18/A/ST9/00746 and 2019/35/B/ST9/03944.
DATA AVAILABILITY
There are no new data associated with this article.