Escape, capture, and levitation of matter in Eddington outbursts
^{1}
Physics DepartmentGothenburg University,
41296
Göteborg,
Sweden
email: gusstaad@student.gu.se, marek.abramowicz@physics.gu.se
^{2}
Nicolaus Copernicus Astronomical Center, ul. Bartycka 18,
00716
Warszawa,
Poland
email:
wlodek@camk.edu.pl
^{3}
Institute of Micromechanics and Photonics,
ul. św A. Boboli 8,
02525
Warszawa,
Poland
email:
maciek.wielgus@gmail.com
Received:
28
March
2013
Accepted:
26
May
2013
Context. An impulsive increase in luminosity by one half or more of the Eddington value will lead to ejection of all optically thin plasma from Keplerian orbits around the radiating star, if gravity is Newtonian and the PoyntingRobertson drag is neglected. Radiation drag may bring some particles down to the stellar surface. On the other hand, general relativistic calculations show that gravity may be balanced by a sufficiently intense radiation field at a certain distance from the star. Aims. We investigate the motion of test particles around highly luminous stars to determine conditions under which plasma may be ejected from the system. Results. In Einstein’s gravity, if the outburst is close to the Eddington luminosity, all test particles orbiting outside an “escape sphere” will be ejected from the system, while all others will be captured from their orbits onto the surface of another sphere, which is well above the stellar surface, and may even be outside the escape sphere, depending on the value of luminosity. Radiation drag will bring all the captured particles to rest on this “Eddington capture sphere”, where they will remain suspended in an equilibrium state as long as the local flux of radiation does not change and remains at the effective Eddington value.
Key words: stars: neutron / stars: winds, outflows / Xrays: binaries / scattering / accretion, accretion disks
© ESO, 2013
1. Introduction
It is widely believed that in luminous compact sources, such as the accreting neutron stars in lowmass Xray binaries, PoyntingRobertson drag will tend to increase the accretion rate as the luminosity of the source increases. Not intending to disagree that in steady sources radiation drag may enhance the accretion rate, here we wish to point out that a sudden increase in luminosity may in fact have the opposite effect in optically thin plasma: atoms (or ions) subjected to increased radiation pressure will tend to move out, away from the stellar surface, PoyntingRobertson drag notwithstanding.
In Newtonian dynamics, the effects of increased radiation pressure are straightforward to compute. A simple application of orbital mechanics reveals that an increase in luminosity by one half of the Eddington value^{1} is sufficient to eject optically thin plasma from its Keplerian orbits at any initial distance from the star.
In Einstein’s general relativity (GR) this remains qualitatively true at a large enough distance from the star – in the optically thin regime a sudden increase in luminosity of sufficient magnitude will unbind the orbiting fluid. However, closer to the star, competing effects of gravity and radiation pressure may lead to very nonintuitive behavior, particularly when combined with the effects of radiation drag. We show in an exact GR calculation that if the luminosity at infinity is close to the Eddington value, optically thin matter in the vicinity of a compact star is neither ejected from the system nor accreted onto the star. Instead, it will collapse onto the surface of an imaginary sphere, whose radius is a function of the stellar radius and the luminosity (Eq. (3)). Rather than orbit the star, fluid on this Eddington capture sphere remains at rest in a state of equilibrium (Stahl et al. 2012), suspended above the stellar radius, in what is to all intents and purposes a state of levitation.
2. The Eddington capture sphere
In Einstein’s general relativity, the redshifted luminosity L(r) of an isotropic star decreases according to (1)The effective Eddington luminosity at r, i.e., the luminosity at which radiation pressure balances gravity, is (2)For a given luminosity L_{∞}, a static balance between the radiation pressure force and gravitational pull is then achieved at the radius at which L(r) = L_{eff}(r) (Phinney 1987). The formula for this unique radius, r_{Edd}(L_{∞}), can be written in the form (3)We use standard Schwarzschild coordinates, with the Schwarzschild radius R_{Schw} = 2GM/c^{2}, and L_{∞}/L_{Edd} the stellar luminosity at infinity in units of the Eddington value, L_{Edd} = 4πGMm_{p}c/σ, where σ is the Thomson crosssection. Strictly speaking the formula is true for an optically thin shell of hydrogen plasma, but similar formulae hold for other compositions of the fluid, as well as other photon momentum absorption processes. Abramowicz et al. (1990) showed that r_{Edd} is a position of stable equilibrium for test particles (with proton mass, and Thomson crosssection for photon momentum absorption) moving radially in the combined gravitational and radiation fields of a spherical, nonrotating, isotropically radiating compact star.
It is now known that the sphere at r_{Edd} is also a position of equilibrium with respect to nonradial perturbations, and that in fact particles from a wide class of trajectories will be captured on this sphere, because of radiation drag (see e.g., Bini et al. 2009; Oh et al. 2011; Stahl et al. 2012,for work in the Schwarzschild metric). Hence, we refer to it as the Eddington capture sphere (ECS). For Kerr metric results, see e.g., Oh et al. (2010); Bini et al. (2011).
3. Equations of motion
We carry out all calculations in the Schwarzschild metric using the standard Schwarzschild spherical coordinates, and extend the treatment of Abramowicz et al. (1990) to nonradial motion; we use their stressenergy tensor of radiation T^{(i)(k)}, which was calculated in the stationary observer tetrad assuming isotropic emission from the stellar surface. An explicit form of the stressenergy tensor is given in the Appendix, along with a summary of the derivation of the equations of motion. Spherical symmetry of the problem assures that for any set of initial conditions the trajectory of a test particle is confined to a single plane, which we take to be the equatorial plane θ = π/2.
From this point on we take G = 1 = c, and scale all radii with M (i.e., one half of the Schwarzschild radius). The proper time is given by the Schwarzschild metric restricted to the equatorial plane (4)We introduce B = 1−2/r, and denote the fourvelocity by u^{i}, so u^{t} = dt/dτ. The angular frequency observed at infinity is Ω = (dφ/dτ)/u^{t}, when discussing circular orbits we use its Keplerian value Ω_{K}(r) = 1/r^{3/2} ≡ v_{K}(r)/r. In the figures, we find it convenient to use v^{φ} ≡ rΩ, scaled with the coordinate speed of light, at those (initial) moments when the radial velocity is zero, and v^{r} ≡ dr/dt when the azimuthal velocity is zero. The stellar radius is denoted by R.
We make the standard assumption of efficient electromagnetic coupling between protons and electrons and, as is usual in the discussion of Eddington luminosity, take the mass of the test particle to be equal to that of the proton. In so doing, we neglect possible effects of charge separation discussed by Walker (1988). The effects of PoyntingRobertson drag in steady disk accretion were considered by Walker & Mészáros (1989). We consider a timedependent problem. On the assumption that the test particle is absorbing momentum (from the radiation field) with the Thomson crosssection in its rest frame, thus suffering a restframe force of (σ/c) × flux, the particle trajectory is described by two coupled, second order differential equations: with (7)The parameter k is the stellar surface luminosity in Eddington units, (8)Because BT^{(r)(t)} = πI(R)(1−2/R)(R^{2}/r^{2}), Eqs. (5), (6) allow a static solution r(τ) = r_{Edd} (Eq. (3)), provided that (9)By integrating Eqs. (5) and (6), we are able to find the particle orbits, such as the ones illustrated in Fig. 1
All numerical results presented below were obtained with the DormandPrince method, which is a fourthorder accuracy, adaptive RungeKutta type integrator.
Fig. 1 Trajectories of particles leaving their circular orbits following an impulsive increase in luminosity from zero to that value (k = 1.45 for R = 5) at which the escape sphere and the Eddington capture sphere coincide. Particles initially in circular orbits at (or outside) the escape sphere leave the system on outward bound trajectories (thin dashed red line). Those from initial orbits inside the escape sphere, turn back and are captured by the ECS, on which they remain at rest in their final position (thick dots). 

Open with DEXTER 
Fig. 2 Escape speeds normalized to the speed of light, as a function of the starting radius, for R = 5,L_{∞} = 0.9 L_{Edd}. Solid line corresponds to radial motion. The Eddington capture sphere is at the position marked with a vertical dashed line. 

Open with DEXTER 
4. Escape from the vicinity of a luminous star
In the absence of radiation, it is easy enough to determine whether or not a test particle will leave the system. A particle is bound or unbound depending on its angular momentum and energy. An unbound test particle is certain to escape to infinity, if only it is outward going, i.e., if u^{r}(0) ≥ 0, while incoming particles may, or may not, be captured by the star depending on their impact parameter.
The presence of radiation drag complicates matters, as energy and angular momentum are no longer constants of the motion. One needs to compute the actual trajectory for given initial conditions to find out whether or not the particle leaves the system. In this paper we only consider trajectories with u^{r}(0) ≥ 0. For examples of inward bound trajectories see, e.g., Stahl et al. (2012).
We have integrated the equations of motion to find trajectories that are not bounded^{2} and found the minimum (“escape”) velocities required for escape from the system. These depend on the initial velocity vector, with the highest value of the escape speed corresponding to initial velocities tangent to a circle concentric with the spherical star (“purely azimuthal” initial motion), and the lowest for radial motion (v^{φ} = 0). These two limiting cases are shown in Fig. 2 as a function of the radius when the stellar radius and luminosity are R = 5, and L_{∞} = 0.9 L_{Edd} or k = 1.5.
Technically, the criterion we used to decide whether or not a test particle escapes to infinity was to integrate the equations of motion until the particle turned around (u^{r} < 0), or conversely, reached r = 1500 R_{Schw}. In the latter case, motion was always radial to a high accuracy, and radiation drag negligible. The particle was deemed to be on a bounded trajectory if its (relativistic) specific energy at r = 1500 R_{Schw} was less than unity. If the specific energy was ≥1, the particle was classified as escaping the system.
Except at the very highest luminosities, it is rather hard for particles to escape from the ECS. Figure 3 shows the values of escape velocities from the ECS, as a function of the luminosity parameter k of Eq. (8).
Fig. 3 Escape speeds (normalized to the speed of light) from the Eddington capture sphere vs. luminosity for R = 5. The whole range of luminosities is shown, for which there is an ECS. 

Open with DEXTER 
In the limit of L_{∞} → L_{Edd}, the ECS becomes infinitely large, r_{Edd} → ∞, and the escape velocity goes to zero. In the opposite (lower) limit of Eq. (9), r_{Edd} = R. In this case, the escape speed from the stellar surface, which now coincides with the ECS, exceeds one half of the speed of light (Fig. 4) for all R ≤ 5. This is not an unusually high value, if one recalls (Kluźniak & Wagoner 1985) that in the Schwarzschild metric the highest Keplerian orbital speed is c/2 (in the marginally stable orbit, a.k.a. the ISCO, at r = 6), however, at first sight it may seem surprising that the escape speed remains so high even in the presence of the effective Eddington luminosity of radiation, whose pressure balances gravity, Eq. (2). As a reminder we note that in this context the escape speed has a somewhat different meaning than usually. It is not the speed necessary for the particle to be unbound, i.e., its energy to be larger than its rest mass. Rather, it is that speed which is sufficient to “overcome” radiation drag, so that the orbit is not bounded, i.e., a speed sufficient for the particle to escape to infinity in spite of continual loss of angular momentum to the radiation field.
5. Impulsive change in luminosity
The question we are addressing in this paper is whether or not an optically thin corona will be ejected from a system undergoing an impulsive increase in luminosity to the effective Eddington value, and what becomes of the particles which are not ejected. Astrophysically, such a situation may be encountered, e.g., in Xray bursters (see Sect. 6). In Newtonian mechanics and gravity, this question was answered by Kluźniak (2013), who found that in the absence of radiation drag all particles will be ejected if the luminosity increases above a critical value, equal to (L + L_{Edd})/2, where L was the luminosity prior to the increase. To solve the problem in full GR with radiation drag included, we integrate Eqs. (5) and (6).
Consider the trajectory of a test particle in a circular orbit at r about a star whose luminosity will undergo an impulsive change from zero to L(R). Thus, the initial conditions correspond to a particle with v^{r} = 0, and v^{φ} = v_{K} appropriate for that circular orbit, v_{K} being the Keplerian speed defined following Eq. (4). We begin the integration at that moment in which the luminosity at r jumps from zero to L(r) of Eq. (1), and find that this initial v^{φ} is higher than the azimuthal escape speed for all r greater than some critical value, which is a function of R and L(R), while it is lower for all radii below that critical value. Accordingly, we introduce the escape sphere dividing space into an outer region, from which all test particles will be hurled out to infinity by the impulsive change in luminosity, and the inner region, between the stellar surface and the escape sphere, from which the particles do not escape the system (cf. Fig. 5).
Fig. 4 Escape speeds from the Eddington capture sphere when it coincides with the stellar surface, as a function of the stellar radius, for purely radial and azimuthal initial velocities. 

Open with DEXTER 
Fig. 5 Radii of the escape sphere and the Eddington capture sphere as a function of the stellar luminosity (R = 5). The two spheres coincide at the critical value of k = 1.45, or L_{∞} = 0.873 L_{Edd}, corresponding to r_{Edd} = 8.40. 

Open with DEXTER 
The fate of particles within the escape sphere is rather interesting. Because of radiation drag, for luminosities below the lower limit of Eq. (9), the particles would have ended up on the stellar surface. However, for the luminosities considered here they all come to rest on the ECS. This is illustrated in Fig. 1, which shows an escaping trajectory of a particle previously in circular orbit on what is now the escape sphere, as well as bounded trajectories from circular orbits of lower radii.
As it turns out, the escape sphere may be inside or outside the ECS. Thus, for a high enough luminosity, all particles initially outside the escape sphere escape to infinity, while those inside it migrate outside to settle on the ECS. Figure 5 shows the radii of the escape sphere and of the capture sphere, as a function of the luminosity parameter. The two coincide for a critical value of the luminosity parameter, which we find to be k_{crit} = 1.45 for R = 5, the orbits shown in Fig. 1 correspond to this case.
For particles that are initially well within the escape sphere, the migration time to the ECS is rather short. This is illustrated in Fig. 6. For a 2 M_{⊙} star the unit of time is 2 GM_{⊙}/c^{3} ≈ 10^{5} s. Thus, all times shown in the figure correspond to a fraction of a second, and this is how long the high luminosity has to be sustained for the particles to collect on the ECS.
Fig. 6 Infall time onto the Eddington capture sphere from various initial radii within the escape sphere. Note that the tick marks on the vertical scale are separated by 10^{4} GM/c^{3}. 

Open with DEXTER 
6. Discussion
While our computations were performed for test particles, the main results should also hold for plasma, in which internal dissipation of energy may occur. Clearly, the particles that eventually settle on the ECS lose all their kinetic energy (and angular momentum), and it makes little difference whether some of it is first redistributed to other particles.
As the radiation front of increased flux travels at the speed of light, all plasma particles are accelerated and the outbound particles do not run into other particles departing from larger radii. The particles that are on unbounded trajectories very quickly approach radial motion, because of radiation drag. Very little plasma compression and dissipation is likely to occur in this situation. In any case, the corona is assumed to be radiatively inefficient, so any dissipation of energy does not lead to binding plasma that is unbound – the total energy of the plasma remains positive. However, there are two serious limitations of the current work, which should be overcome in future calculations. We have assumed spherical symmetry, while the most luminous, and interesting, astronomical sources accrete in a disklike geometry. Further, in the microscopic domain we have neglected the transfer of momentum between the scattered (or reemitted) photons and the test particle – only the momentum transferred by the absorbed photons was included in the calculation. Thus the radius of the escape sphere is likely to have been underestimated. On the other hand, the position of the ECS is insensitive to this assumption, particles at rest on the ECS radiate isotropically, and the only momentum transfer is from the absorption process included in our calculation.
The future will tell whether the theoretical finding presented here have any application to Xray bursts, which are spherically symmetric to a high degree, and how they affect the estimates of neutron star radius, as in, e.g., Damen et al. (1990). As the precise value of the stellar radius and mass is not known, it is hard at present to translate the observed flux from the neutron star to a precise value of the luminosity k parameter. Should a direct observational determination of k turn out to be possible, a value for the stellar radius would immediately follow (cf. Eq. (8)). However, we note that Xray bursters are thought to attain Eddington luminosity, e.g., Lewin et al. (1993), and a recent burst in the neutron star 2S 0918549 has exhibited an interesting achromatic oscillation, which has been interpreted as an accretion disk response to superEddington fluxes in an Xray burst (in ’t Zand et al. 2011). Wielgus et al. (2012) point out that, in principle, a luminosity oscillation may be related to the appearance and disappearance of the ECS. Very interesting luminosity variations were also reported in the Z sources, or rather, in the Z phase of bright LMXBs (Lin et al. 2009), and in some microquasars – of course, these are precisely the sources where spherical symmetry is broken, and more work is needed before any reliable conclusions can be reached on the relevance to these sources of the phenomena discussed here.
Perhaps it is worth noting that superEddington luminosity may also be present among some accreting black holes, notably ULXs (Feng & Soria 2011; Gladstone et al. 2011), but also in supermassive black holes, (e.g., Kollmeier et al. 2011; Nobuta et al. 2012), and microquasars, such as V4641 Sgr in outburst (Revnivtsev et al. 2002).
7. Conclusions
We have shown in a GR calculation including radiation drag that a sudden change in luminosity from zero to about the (effective) Eddington value leads to ejection of all test particles outside a certain sphere (the escape sphere). In Newtonian physics all particles outside the stellar surface are either ejected, or not, depending on the luminosity jump (Kluźniak 2013), when radiation drag is neglected. In Einstein’s gravity, the redshifted luminosity drops with radius, and this effect, combined with radiation drag, introduces a radius dependence to the problem, thus allowing some particle to be left behind in the system. At the high final luminosities considered here, the nonescaping particles are not accreted by the star. Instead, they lose all kinetic energy and angular momentum and migrate to the ECS, typically within a fraction of a second, where they remain at rest as long as the local flux of radiation corresponds to the effective Eddington value.
We expect our work to be applicable to the optically thin coronae of variable Xray sources. Future work will show how the position of the escape sphere depends on the initial luminosity of the system (Mishra et al. 2013).
More precisely, an impulsive increase in luminosity by one half of the difference between the most recent steady value and the Eddington luminosity is sufficient to unbind a nonrelativistic test particle (Kluźniak 2013).
Acknowledgments
We thank Saul Rappaport and Wenfei Yu for inspiring discussions. In addition, Saul Rappaport’s help was invaluable in validating our codes. Research supported in part by Polish NCN grants UMO2011/01/B/ST9/05439, and N N203 511238, and the Czech grant MSM 4781305903, as well as by a Swedish VR grant.
References
 Abramowicz, M. A., Ellis, G. F. R., & Lanza, A. 1990, ApJ, 361, 470 [NASA ADS] [CrossRef] [Google Scholar]
 Bini, D., Jantzen, R. T., & Stella, L. 2009, Class. Q. Grav., 26, 055009 [NASA ADS] [CrossRef] [Google Scholar]
 Bini, D., Geralico, A., Jantzen, R. T., Semerák, O., & Stella, L. 2011, Class. Q. Grav., 28, 035008 [NASA ADS] [CrossRef] [Google Scholar]
 Damen, E., Wijers, R. A. M. J., van Paradijs, J., et al. 1990, A&A, 237, 103 [NASA ADS] [Google Scholar]
 Feng, H., & Soria, R. 2011, New Astron. Rev., 55, 166 [NASA ADS] [CrossRef] [Google Scholar]
 Gladstone, J. C., Roberts, T. P., & Done, C. 2011, Astron. Nachr., 332, 345 [NASA ADS] [CrossRef] [Google Scholar]
 in ’t Zand, J. J. M., Galloway, D. K., & Ballantyne, D. R. 2011, A&A, 525, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kluźniak, W. 2013, A&A, 551, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kluźniak, W., & Wagoner, R. V. 1985, ApJ, 297, 548 [NASA ADS] [CrossRef] [Google Scholar]
 Kollmeier, J. A., Onken, C., Kochanek, C. S., et al. 2011, ApJ, 648, 128 [NASA ADS] [CrossRef] [Google Scholar]
 Lewin, W. H. G., van Paradijs, J., & Taam, R. E. 1993, Space Sci. Rev., 62, 223 [NASA ADS] [CrossRef] [Google Scholar]
 Lin, D., Remillard, R. A., & Homan, J. 2009, ApJ, 696, 1257 [NASA ADS] [CrossRef] [Google Scholar]
 Oh, J. S., Kim, H., & Lee, H. M. 2010, Phys. Rev. D, 81, 084005 [NASA ADS] [CrossRef] [Google Scholar]
 Oh, J. S., Kim, H., & Lee, H. M. 2011, New Astron., 16, 3, 183 [Google Scholar]
 Mishra, B., & Kluzniak, W. 2013, A&A, submitted [Google Scholar]
 Nobuta, K., Akiyama, M., Ueda, Y., et al. 2012, ApJ, 761, 143 [NASA ADS] [CrossRef] [Google Scholar]
 Phinney, E. S. 1987, in Superluminal Radio Sources, eds. J. A. Zensus, & T. J. Pearson (Cambridge: Cambridge University Press), 12 [Google Scholar]
 Revnivtsev, M., Gilfanov, M., Churazov, E., & Sunyaev, R. 2002, A&A, 391, 1013 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Stahl, A., Wielgus, M., Abramowicz, M., Kluźniak, W., & Yu, W. 2012, A&A, 546, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Walker, M. A. 1988, ApJ, 330, L47 [NASA ADS] [CrossRef] [Google Scholar]
 Walker, M. A., & Mészáros, P. 1989, ApJ, 346, 844 [NASA ADS] [CrossRef] [Google Scholar]
 Wielgus, M., Stahl, A., Abramowicz, M., & Kluźniak, W. 2012, A&A, 545, 123 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A
Here, we present an abbreviated derivation of Eqs. (5), (6). The acceleration of the particle is in general given by (A.10)where ds^{2} = M^{2}dτ^{2}. Taking the force of radiation to be (A.11)we obtain the equation of motion (A.12)We determine the radiation flux F^{i} from the expression (A.13)
where is the projection tensor, and T^{jk} is the stressenergy tensor for the radiation. Thus, the flux is (A.14)with (A.15)Expressing ε ≡ ϵ/ [πI(r)] in tetrad components, we obtain Eq. (7). We take the tetrad components of the stressenergy tensor from Abramowicz et al. (1990): (A.16)(A.17)(A.18)(A.19)with all other components except T^{(θ)(θ)} vanishing. The T^{(θ)(θ)}component, when it appears, will however always be multiplied by u^{θ} = 0, so we may ignore this component for the present problem. In the equations above, α = α(r) is the viewing angle and I(r) is the intensity of radiation. In our dimensionless coordinates, they are given by (A.20) (A.21)
All Figures
Fig. 1 Trajectories of particles leaving their circular orbits following an impulsive increase in luminosity from zero to that value (k = 1.45 for R = 5) at which the escape sphere and the Eddington capture sphere coincide. Particles initially in circular orbits at (or outside) the escape sphere leave the system on outward bound trajectories (thin dashed red line). Those from initial orbits inside the escape sphere, turn back and are captured by the ECS, on which they remain at rest in their final position (thick dots). 

Open with DEXTER  
In the text 
Fig. 2 Escape speeds normalized to the speed of light, as a function of the starting radius, for R = 5,L_{∞} = 0.9 L_{Edd}. Solid line corresponds to radial motion. The Eddington capture sphere is at the position marked with a vertical dashed line. 

Open with DEXTER  
In the text 
Fig. 3 Escape speeds (normalized to the speed of light) from the Eddington capture sphere vs. luminosity for R = 5. The whole range of luminosities is shown, for which there is an ECS. 

Open with DEXTER  
In the text 
Fig. 4 Escape speeds from the Eddington capture sphere when it coincides with the stellar surface, as a function of the stellar radius, for purely radial and azimuthal initial velocities. 

Open with DEXTER  
In the text 
Fig. 5 Radii of the escape sphere and the Eddington capture sphere as a function of the stellar luminosity (R = 5). The two spheres coincide at the critical value of k = 1.45, or L_{∞} = 0.873 L_{Edd}, corresponding to r_{Edd} = 8.40. 

Open with DEXTER  
In the text 
Fig. 6 Infall time onto the Eddington capture sphere from various initial radii within the escape sphere. Note that the tick marks on the vertical scale are separated by 10^{4} GM/c^{3}. 

Open with DEXTER  
In the text 