# RADIATIVELY–DRIVEN OUTFLOWS AND AVOIDANCE OF COMMON–ENVELOPE EVOLUTION IN CLOSE BINARIES

## 1. INTRODUCTION

The question of what happens to a compact object that is fed mass at rates far higher than its Eddington limit has a long history (Shakura & Sunyaev, 1973; Kafka & Mészáros, 1976; Begelman, 1979). In the context of accreting binary systems, this problem is particularly acute because of the possibility of common–envelope (CE) evolution at such rates. That is, the accreting component may be unable either to accept or to expel the mass at a sufficiently high rate to avoid the formation of an envelope engulfing the entire binary system. The frictional drag of this envelope can shrink the binary orbit drastically. If the resulting release of orbital energy is enough to unbind the envelope, the binary will emerge from the common envelope with a smaller separation; if not, the binary components may coalesce. CE evolution is probably required for the formation of binaries such as cataclysmic variables, in which the binary separation is far smaller than the radius of the accreting white dwarf’s red giant progenitor. However, it is in general an open question whether CE evolution occurs in any given binary.

This question is thrown into sharp relief by recent work on the evolution of the low–mass X–ray binary Cygnus X–2 (King & Ritter, 1999), which has a period of 9.84 d. The rather precise spectroscopic information found by Casares, Charles, & Kuulkers (1998), together with the observed effective temperature of the secondary, shows that this star has a mass definitely below $`0.7\mathrm{M}_{}`$ and yet a luminosity of order $`150\mathrm{L}_{}`$. King & Ritter (1999) consider several possible explanations and show that the only viable one is that Cygnus X–2 is a product of early massive Case B evolution. Here ‘Case B’ means that the mass–losing star has finished core hydrogen–burning, and is expanding across the Hertzsprung gap: ‘early’ means that the stellar envelope is radiative rather than convective, and ‘massive’ that the helium core is non–degenerate; see Kippenhahn & Weigert, 1967. In Cygnus X–2 an initially more massive ($`M_{2i}3.5\mathrm{M}_{}`$) secondary transferred mass on a thermal timescale ($`10^6`$ yr) to the neutron star. This idea gives a satisfying fit to the present observed properties of Cyg X–2, as well as a natural explanation for the large white dwarf companion masses found in several millisecond pulsar binaries with short orbital periods. CE evolution cannot have occurred, as Cyg X–2’s long orbital period means that there was far too little orbital energy available for the CE mechanism to have ejected so much mass. Thus an inescapable feature of this picture is that the neutron star is evidently able to eject essentially all of the matter ($`23\mathrm{M}_{}`$) transferred to it at highly super–Eddington rates $`10^6\mathrm{M}_{}`$ yr<sup>-1</sup>. Indeed, the neutron star mass in Cyg X–2 is rather close to the canonical value of $`1.4\mathrm{M}_{}`$. The aim of this paper is to determine under what conditions such expulsion can occur without the system going into a common envelope.

## 2. EXPULSION BY RADIATION PRESSURE

There are essentially two views as to the fate of matter dumped onto a compact object at a highly super–Eddington rate. In spherically–symmetric, dissipative accretion of an electron–scattering medium, the luminosity generated by infall down to radius $`R`$ will reach the Eddington limit at a radius

$$R_{\mathrm{ex}}\left(\frac{\dot{M}_{\mathrm{tr}}}{\dot{M}_{\mathrm{Edd}}}\right)R_S,$$
(1)   
where $`\dot{M}_{\mathrm{tr}}`$ is the mass infall rate at large radius (i.e. the mass transfer rate from the companion star in our case), $`\dot{M}_{\mathrm{Edd}}=L_{\mathrm{Edd}}/c^2`$ is the Eddington accretion rate, and $`R_S`$ is the Schwarzschild radius (Begelman, 1979). This is also the “trapping radius”, below which photon diffusion outward cannot overcome the advection of photons inward. If the compact object is a black hole, the radiation generated in excess of the Eddington limit can thus be swept into the black hole, and lost. If the compact object is a neutron star, however, radiation pressure building up near the star’s surface must resist inflow in excess of $`\dot{M}_{\mathrm{Edd}}`$, causing the stalled envelope to grow outward. This situation would lead to the formation of a common envelope.

The outcome may be very different if the accretion flow has even a small amount of angular momentum. Shakura & Sunyaev (1973) suggested that super–Eddington flow in an accretion disk would lead to the formation of a strong wind perpendicular to the disk surface, which could carry away most of the mass. Such a model (an “Adiabatic Inflow-Outflow Solution,”, or ADIOS) was elaborated by Blandford & Begelman (1999: hereafter BB99), who considered radiatively inefficient accretion flows in general. BB99 recalled that viscous transfer of angular momentum also entails the transfer of energy outward. If the disk were unable to radiate efficiently (as would be the case at $`R<R_{\mathrm{tr}}`$), the energy deposited in the material well away from the inner boundary would unbind it, leading to the creation of powerful wind. BB99 described a family of self-similar models in which the mass inflow rate decreases inward as $`\dot{M}r^n`$ with $`0<n<1`$. The exact value of $`n`$ depends on the physical processes depositing energy and angular momentum in the wind. If these are very efficient (e.g., mediated by highly organized magnetic torques) $`n`$ could be close to zero, in which case little mass would be lost. However, if the wind is produced inefficiently, $`n`$ would have to be close to 1 and the mass flux reaching the central parts of the accretion disk would be much smaller than the mass transferred from the secondary. For example, two-dimensional hydrodynamical simulations of the evolution of a non-radiative viscous torus (Stone, Pringle & Begelman 1999) show the development of a convectively driven circulation with little mass reaching the central object, and $`n1`$. The development of this strong mass loss is generic and is not related to the assumption of self-similarity. In effect, what is happening is that the energy liberated by a small fraction of the mass reaching the deep gravitional potential serves to unbind the majority of the matter which is weakly bound at large distances.

While the specific details of mass loss from super–Eddington flow have not been worked out (in particular, radiation-dominated convection is poorly understood), it is reasonable to assume that the wind will be produced inefficiently, with $`n1`$ as in the case of hydrodynamic convection. We also assume that most of the matter will be blown away from $`R_{\mathrm{ex}}`$. Applying equation (1), we find

$$R_{\mathrm{ex}}1.3\times 10^{14}\dot{m}_{\mathrm{tr}}\mathrm{cm},$$
(2)   
where $`\dot{m}_{\mathrm{tr}}`$ is the mass transfer rate expressed in $`\mathrm{M}_{}`$ yr<sup>-1</sup>. Note that $`R_{\mathrm{ex}}`$ is independent of the mass of the compact accretor. Since we restrict attention to electron scattering opacity only, we require that hydrogen should be strongly ionized at $`R_{\mathrm{ex}}`$. This is ensured by requiring the radiation temperature near $`R_{\mathrm{ex}}`$ to exceed $`T_H10^4`$ K. Since the luminosity emerging from $`R_{\mathrm{ex}}`$ is close to the Eddington limit, we require $`L_{\mathrm{Edd}}4\pi R_{\mathrm{ex}}^2\sigma T_H^4`$, which is satisfied if $`R_{\mathrm{ex}}/R_S10^7m_1^{1/2}`$, or equivalently (using equation)

$$\dot{M}_{\mathrm{tr}}10^7\dot{M}_{\mathrm{Edd}}m_1^{1/2}2\times 10^2m_1^{1/2}\mathrm{M}_{}\mathrm{yr}^1$$
(3)   
where $`m_1=M_1/\mathrm{M}_{}`$ is the mass of the compact accretor (black hole or neutron star).

CE evolution will be avoided if $`R_{\mathrm{ex}}`$ is smaller than the accretor’s Roche lobe radius $`R_1`$. If the accretor is the less massive star (as will generally hold in cases of interest) we can use standard formulae to write

$$r_1=1.9m_1^{1/3}P_\mathrm{d}^{2/3},$$
(4)   
where $`r_1=R_1/\mathrm{R}_{}`$ and $`P_\mathrm{d}`$ is the orbital period measured in days. Combining with equation (2) gives

$$\dot{M}_{\mathrm{tr}}10^3m_1^{1/3}P_\mathrm{d}^{2/3}\mathrm{M}_{}\mathrm{yr}^1.$$
(5)   
This form of the limit can be compared directly with observation if we have estimates of the transfer rate, orbital period and the accretor mass. For more systematic study it is useful to replace the dependence on the accretor’s Roche lobe by that on its companion’s. Thus, since the mass transfer rate is specified by properties of the companion star, which is assumed to fill its Roche lobe radius $`R_2`$, we eliminate $`R_1`$ from the condition $`R_{\mathrm{ex}}R_1`$ by using the relation

$$\frac{R_1}{R_2}=\left(\frac{m_1}{m_2}\right)^{0.45},$$
(6)   
(cf King et al., 1997) where $`M_2=m_2\mathrm{M}_{}`$ is the companion mass. Writing $`R_2=r_2\mathrm{R}_{}`$ we finally get the limit

$$\dot{M}_{\mathrm{tr}}5\times 10^4m_1^{0.45}m_2^{0.45}r_2\mathrm{M}_{}\mathrm{yr}^1$$
(7)   
on the mass transfer rate if CE evolution is not to occur.

## 3. AVOIDANCE OF COMMON ENVELOPE EVOLUTION

By specifying the nature of the companion star we fix $`m_2,r_2`$ and $`\dot{M}_{\mathrm{tr}}`$ in (7), and so can examine whether CE evolution is likely in any given case. Rapid mass transfer occurs if the companion star is rather more massive than the accretor, since then the act of transferring mass shrinks the donor’s Roche lobe. The mass transfer proceeds on a dynamical or thermal timescale depending on whether the donor star’s envelope is largely convective or radiative (e.g. Savonije, 1983). In the first case, CE evolution is quite likely to ensue, as the mass transfer rate rises to very high values. However, even in this case it is worth checking the inequality (7) in numerical calculations, as the e–folding time for the mass transfer is $`t_e(H/R_2)t_M`$, where $`H`$ is the stellar scaleheight and $`t_M`$ is the mass transfer timescale set by whatever process (e.g., nuclear evolution) brought the donor into contact with its Roche lobe initially. For main–sequence and evolved stars we have $`H/R_210^4,10^2`$ respectively. Thus $`t_e`$ may be long enough that the companion mass is exhausted before (7) is violated.

Thermal–timescale mass transfer is rather gentler, and offers the possibility of avoidance of CE evolution. In addition to the case mentioned above, thermal–timescale mass transfer will also occur if the donor star is crossing the Hertzsprung gap and has not yet developed a convective envelope (i.e., is not close to the Hayashi line), even if it is the less massive star. Detailed calculations (Kolb, 1998) show that in both cases the mass transfer rate is given roughly by

$$\dot{M}_{\mathrm{tr}}\frac{M_2}{t_{\mathrm{KH}}},$$
(8)   
where

$$t_{\mathrm{KH}}=3\times 10^7\frac{m_2^2}{r_2l_2}\mathrm{yr}$$
(9)   
was the Kelvin–Helmholtz time of the star when it left the main sequence, and $`L_2=l_2\mathrm{L}_{}`$ was its luminosity. (Note that by definition the donor is not in thermal equilibrium, so an originally main–sequence donor will develop a non–equilibrium structure as mass transfer proceeds.) The condition of a radiative envelope requires a main–sequence mass $`m_21`$, so we may take

$$r_2m_2^{0.8},l_2m_2^3.$$
(10)   
Inserting in (9) and (8) we find

$$\dot{M}_{\mathrm{tr}}3\times 10^8m_2^{2.8},$$
(11)   
so comparing with (7) we require

$$m_253m_1^{0.18}$$
(12)   
and thus (from 11)

$$\dot{M}_{\mathrm{tr},\mathrm{max}}2\times 10^3m_1^{0.51}\mathrm{M}_{}\mathrm{yr}^1.$$
(13)   
Hence we expect CE evolution to be avoided in thermal–timescale mass transfer from a main–sequence star, or from a Hertzsprung gap star, provided that it has a radiative envelope. This is in agreement with the assumption of no CE evolution in Cyg X–2 made by King & Ritter (1999), where the initial donor mass was about $`3.5\mathrm{M}_{}`$.

## 4. CONCLUSIONS

We have derived a general criterion for the avoidance of common–envelope evolution in a binary in which the accretor is a neutron star or a black hole. This shows that thermal–timescale mass transfer from a main–sequence star is unlikely to lead to CE evolution, as is mass transfer from a Hertzsprung gap star, provided that the envelope is radiative. The first possibility allows the early massive Case B evolution inferred by King & Ritter (1999) for the progenitor of Cyg X–2. SS433 may be an example of the second possibility, with a fairly massive donor star. We will discuss this possibility in detail in a future paper.

The considerations of this paper suggest that common–envelope evolution with a neutron–star or black–hole accretor generally requires an evolved donor with a deep convective envelope. This represents a slight restriction on some of the routes invoked in the possible formation of Thorne–$`\dot{\mathrm{Z}}`$ytkow objects.

This research was carried out at the Institute for Theoretical Physics and supported in part by the National Science Foundation under Grant No. PHY94–07194. ARK gratefully acknowledges support by the UK Particle Physics and Astronomy Research Council through a Senior Fellowship. MCB acknowledges support from NSF grant AST95–29170 and a Guggenheim Fellowship.