Relativistic effects on radiative ejection of coronae in variable Xray sources
Nicolaus Copernicus Astronomical Center, Polish Academy of Sciences, ul. Bartycka 18, 00716 Warszawa, Poland
email: mbhupe@camk.edu.pl; wlodek@camk.edu.pl
Received: 1 July 2013
Accepted: 13 March 2014
Context. Optically thin coronae around neutron stars suffering an Xray burst can be ejected as a result of rapid increase in stellar luminosity. In general relativity, radiation pressure from the central luminous star counteracts gravitational attraction more strongly than in Newtonian physics. However, motion near the neutron star is very effectively impeded by the radiation field.
Aims. To explore the mechanisms leading to ejection of accretion disk coronae
Methods. We perform a general relativistic calculation of the motion of a test particle in a spherically symmetric radiation field.
Results. At every radial distance from the star larger than that of the innermost stable circular orbit (ISCO), and any initial luminosity of the star, there exists a luminosity change which leads to coronal ejection. The luminosity required to eject from the system the inner parts of the optically thin neutronstar corona is very high in the presence of radiation drag and always close to the Eddington luminosity. Outer parts of the corona, at a distance of 20 R_{G} or more, will be ejected by a subEddington outburst. Mildly fluctuating luminosity will lead to dissipation in the plasma and may explain the observed Xray temperatures of coronae in low mass Xray binaries. At large radial distances from the star (3 × 10^{3} R_{G} or more) the results do not depend on whether or not PoyntingRobertson drag is included in the calculation.
Key words: stars: neutron / stars: winds, outflows / Xrays: binaries / Xrays: bursts / scattering / accretion, accretion disks
© ESO, 2014
1. Introduction
Most low mass Xray binaries (LMXBs) vary in luminosity on many timescales. Particularly rapid and luminous variations are exhibited by the numerous Xray bursters, which are thought to be neutron stars undergoing a thermonuclear explosion on their surface, yielding an Eddington luminosity at maximum. Sometimes the maximum flux corresponds to superEddington luminosities, as in the pulsating neutron star “LMC transient” A0535668 (at a firm distance of 50 kpc), which is thought to attain an isotropic flux L_{∞} = 1.2 × 10^{39} erg/s = 6.9 L_{Edd} for a 1.4 M_{⊙} star (Bradt & McClintock 1983). One may ask whether the variability of Xray luminosity has any influence on the state of circumstellar matter. The observed LMXBs are thought to be powered by accretion occurring through an optically thick disk. In this paper we are considering the response of accretion flows in the optically thin regions of the system to Xray variability of the source.
We are considering a corona in a variable Xray source. In this context, corona means optically thin plasma “above the surface” of the disk. The effects of luminosity change on the motion of optically thin accreting matter in Newtonian dynamics are already known, an impulsive increase in the stellar luminosity of the star may lead to coronal ejection if the luminosity is increased by one half or more of the difference between Eddington luminosity and the initial one Kluźniak (2013). However, that result was derived analytically taking no account of radiation drag. Walker & Meszaros (1989) pointed out the effect of PoyntingRobertson drag Robertson (1937) for accretion flows around high luminosity neutron stars. Proper inclusion of radiation drag requires a numerical solution in general relativity (GR), and we are presenting such solutions in this paper. Recently, Stahl et al. (2013) discussed the problem for superEddington outbursts in initially nonradiating stars. Here, we consider a wide range of initial and final luminosities.
We are attempting to understand luminosity effects on the corona by modeling testparticle motion around the star. The dynamics are described by equations of motion of a particle moving in a spherically symmetric radiation field in the Schwarzschild metric, while interaction with the radiation by momentum absorption with a crosssection whose numerical value will correspond to the Thomson crosssection times the mass of the particle expressed in units of proton mass. This assures that the conventional Eddington luminosity, L_{Edd}, will balance gravity exactly for hydrogen plasma at infinity (i.e., in the Newtonian limit). The calculation can be carried over to other compositions, and other crosssections for photon absorption, by suitably redefining the Eddington luminosity.
The radiation, assumed to be originating on the surface of a neutron star, will have two important effects on the motion: first, it will counteract gravity by transfering the radial component of (a part of) its momentum; second, it will exert a drag on the moving particle. The calculation is performed in GR, and includes all nonvanishing components of the radiation stress tensor.
The paper is organized in the following manner. In Sect. 2 we shall briefly explain the equations of motion that describe the trajectories of test particles. In Sect. 3.1, in order to isolate the gravitational effects of GR, we shall present results of a simplified calculation without radiation drag, and in Sects. 3.2, 3.3 we will describe the test particle behavior determined from the complete equations of motion (e.o.m.), which include radiation drag. In Sect. 4 we discuss the results and present the conclusions in Sect. 5.
2. Equations of motion
We performed all the calculations in the Schwarzschild metric, using spherical polar coordinates (r,θ,φ). For a spherically symmetric radiation field, motion of a testparticle is restricted to one plane and we choose it to be the equatorial plane (θ = π/ 2). Abramowicz et al. (1990) performed a rigorous analysis of purely radial motion of a test particle in the combined gravity and isotropic radiation fields of a spherical, nonrotating, compact star, and we shall use their stress energy tensor of radiation, T^{(μ)(ν)}, calculated in a stationary observer’s tetrad assuming isotropic emission from the surface of star, see also Stahl et al. (2013). Numerical solutions of the equations of motion of test particle trajectories in a steady radiation field have been obtained e.g., by Bini et al. (2009) and Stahl et al. (2012) in the Schwarzschild metric, and by Oh et al. (2010); Bini et al. (2011) in the Kerr metric.
To describe in detail the GR effects of variable luminosity on the motion of optically thin plasma, and specifically the effects of impulsive changes in luminosity on testparticle trajectories, we use two sets of equations in the present paper. One includes all components of the radiation stress tensor, including the terms responsible for radiation drag, and the other, a simplified set, describes the trajectories of test particle without considering radiation drag.
2.1. Complete equations of motion (“with drag”)
We performed the computations using dimensionless coordinates scaled by gravitational radius R_{G} = GM/c^{2}, i.e., the dimensionless radial position, x, stellar radius, X, and proper time, τ, The equations of motion (e.o.m.) are (e.g., Stahl et al. 2012, 2013, see also Appendix A) where we found it convenient to abbreviate the metric coefficient as B ≡ g_{tt} = (1 − 2 /x), and where and the time component of the particle four velocity is given by (5)In all the computations the luminosity at infinity, L_{∞}, is represented in units of Eddington luminosity as k (Eqs. (1), (2)). One can also display the results in terms of the luminosity at the surface of star, L(R). These are related through (6)To obtain x(τ) and φ(τ) we integrated Eqs. (1) and (2) using the code of Stahl et al. (2012) and Wielgus et al. (2012), relying on the DormandPrince method, which is a fourthorder accuracy, adaptive stepsize RungeKutta type integration method. In all our simulations we assumed a fixed radius of neutron star of X = 5. At this stellar radius, the marginally stable radius (ISCO) is outside the surface of the neutron star (Kluzniak & Wagoner 1985).
2.2. Equations of motion without drag
The terms responsible for radiation drag are those which explicitly involve spatial components of velocity and are proportional to a nonvanishing component of the radiation stress tensor. For example, the last term in Eq. (2) does not contribute to radiation drag because it is not proportional to the radiation stress tensor, but all the other terms on the right hand side of this equation do. On the other hand, the first term on the right hand side of Eq. (1) is proportional to the stress tensor, but it does not involve spatial components of velocity, and so it is not a radiation drag term. Neglecting the terms responsible for radiation drag, one obtains simplified equations of motion involving only one component of the radiation stress tensor,
3. Coronal ejection
If a particle is moving in a circular orbit, it can become unbound under a sudden increase of luminosity. For a test particle in Newtonian dynamics, moving initially in a Keplerian circular orbit under initial luminosity k_{1}, the minimum luminosity change required to unbind the particle is (k_{Edd} − k_{1})/2. Thus, in the Newtonian limit the minimum luminosity required to eject a test particle is L_{NL} = (1 + k_{1})L_{Edd}/ 2. An impulsive change of luminosity to this value will also lead to coronal ejection if the radiation drag is neglected Kluźniak (2013).
We will be contrasting conditions under which particles escape to infinity with those in which the particles remain in the system. The results presented below are based on numerical solutions for the trajectories. Therefore, we need a numerical criterion for deciding which outgoing trajectories are bounded, and which ones are unbounded. We adopted one^{1} in which we computed the Newtonian effective specific energy of the test particle at x = 3000. By the effective specific energy we mean the specific energy in the potential −GM_{eff}/r, with M_{eff} = M(1 − L_{∞}/L_{Edd}); this prescription takes account of the radial radiation pressure term. If the Newtonian effective specific energy is positive at x = 3000, the trajectory is deemed to be unbounded (i.e., extending to infinity), and if it is negative, the trajectory is deemed bounded. In practice, at such a large distance the trajectory of those particles which have reached it is always very nearly radial, and the influence of drag is negligible there.
To isolate the effects of drag from gravitational effects in GR we shall discuss two cases, in the first case we shall omit effects of radiation drag (Eqs. (7), (8)), and in the second case we shall use the complete e.o.m. (Eqs. (1), (2)) that include also the radiation drag terms. Our results show convergence of the two cases at large radii (Fig. 3). Interestingly, radiation drag dominates GR corrections, and already at a radius of about 10^{2} R_{G} the nondrag GR corrections appear to be negligible. At r = 10^{3} R_{G} the fully relativistic solution (including drag) differs from the Newtonian limit by just a few percent. In the following we give the value of luminosities at infinity scaled by the Eddington luminosity. Thus, a “luminosity of k_{1}” signifies L_{∞} = k_{1}L_{Edd}.
Fig. 1 Minimum luminosity that will eject to infinity a particle initially in circular orbit at x_{esc}. Plot compares solutions of the full e.o.m. (solid green curve, Eqs. (1), (2)) with solutions of the simplified equations that neglect the effects of radiation drag ([ND], thin red line, Eqs. (7), (8)). Initial luminosity is k_{1} = L_{1}/L_{Edd} = 0, stellar radius X = 5. The Newtonian limit is also shown ([NL], dotdashed line), as well as the radius of the Eddington capture sphere (ECS). See text for details. 

Open with DEXTER 
Fig. 2 Radius of the escape sphere when drag is neglected. The minimum luminosity, k_{esc}, that unbinds a particle in circular orbit at x_{esc} is shown as a function of the initial orbital radius, for two values of initial luminosity, k_{1} = 0.0, and 0.36. Stellar radius is X = 5.0. Constant k_{esc} lines represent Newtonian limit (NL) for k_{1} = 0.0,0.36. If there were no drag, particles from all circular orbits of radii in the range 6 ≤ x(0) <x_{esc}(k_{esc}) would escape the system, and the initial radius in the outgoing trajectory, x(0), would be its periastron as well. 

Open with DEXTER 
Fig. 3 Convergence of results using Eqs. (1), (2) and (6), (7), with the Newtonian limit. The ratio is plotted of the minimum luminosity sufficient to unbind a particle, initially in a circular GR orbit at k_{1} = 0, to the corresponding value in the Newtonian limit. Uppermost (thick green) curve corresponds to Eqs. (1) and (2), and red continuous curve to Eqs. (7) and (8). 

Open with DEXTER 
3.1. GR effects without drag
We consider an initial luminosity, k_{1} which at some point in the evolution of the system will be changed impulsively to a different value k_{2}. To contrast the effects of the Schwarzschild spacetime with the Newtonian case, in this subsection we assume the absence of radiation drag. Thus a test particle can move in a Keplerian orbit at a fixed radial distance from the star, as long as the luminosity is held fixed. It will continue to move on nearby orbits under very small changes in the luminosity. But if the luminosity change is large enough, the particle can escape the system.
We investigate for what luminosity change k_{2} − k_{1} the particle will become unbound (neglecting radiation drag). The initial conditions correspond to a circular orbit, at a particular radius around the star of luminosity k_{1}. To find the trajectory we integrate Eqs. (7), (8). Since the effective radiative force diminishes more rapidly with increasing distance from the center of a variable Xray source than the gravitational attraction, in the absence of drag a smaller luminosity change is required to eject the particle from a circular orbit close to the ISCO (x = 6 for a nonrotating neutron star) than from a more distant one. This is clear in Fig. 2. The plotted value, k_{esc}(x_{esc}), corresponds to the minimum luminosity change k_{esc} − k_{1} sufficient to unbind a particle that was in a circular orbit at x = x_{esc} at luminosity k_{1}, assuming that drag is neglected. Also shown (constant k_{esc} lines in the figure) is the corresponding Newtonian limit of Kluźniak (2013).
It turns out that for any given pair of initial and final luminosities there is a sphere, inside of which the luminosity change is sufficiently large for all orbits to become unbounded–all particles within the sphere abandon their previously stable circular orbits for trajectories extending to infinity. On the other hand, particles previously orbiting outside the sphere remain bound. However, these results are drastically altered by the radiation drag (Sect. 3.2). Another way of looking at Fig. 2 is that it presents the escape sphere radii (in motion without drag) as a function of the final luminosity L_{esc}, for two different initial luminosities.
Fig. 4 Testparticle trajectory beginning at r = 15R_{G} for a constant k = L_{2}/L_{Edd} = 0.048. Instead of maintaining a circular orbit the test particle follows a spiral trajectory and accretes on the surface of star. 

Open with DEXTER 
3.2. GR effects with drag
Now we shall integrate Eqs. (1) and (2) to compute the trajectories of test particles suffering the influence of the full radiation tensor. Radiation drag is automatically included in the formalism, and the angular momentum and energy of the test particle are no longer constants of motion–in this case we cannot get an analytical expression like in Newtonian dynamics by integrating the equations. In a steady radiation field, a circular orbit cannot be maintained–drag will cause the orbiting test particle to loose angular momentum and energy, and will cause it to follow a spiral trajectory (Walker & Meszaros 1989), cf. Fig. 4.
If initially a test particle is moving under the influence of luminosity k_{1}, a sudden change to a different luminosity k_{2} will change its trajectory. For an initially nonluminous star (k_{1} = 0) the discussion is straightforward. The initial trajectory may be taken to correspond to a stable circular orbit, and the parameter k in Eqs. (1) and (2) may be taken to be nonzero starting at any moment. Again, for any radius of the initial circular orbit, there is a minimum luminosity k_{2} = L_{esc}/L_{Edd} which makes the testparticle unbound (Fig. 1). Note that now k_{esc} is a monotonically decreasing function of x. For comparison, also plotted in Fig. 1 are k_{esc} in the Newtonian limit (dotdashed curve, Kluźniak 2013), and in the GR calculation with radiation drag neglected (Eqs. (7), (8) and Sect. 3.1). Also shown (magenta dashed curve) is the effective Eddington luminosity, k_{eff}(x), i.e., that luminosity^{2} at which a test particle located at x may be in a static equilibrium (Phinney 1987). The inverse function x(k_{eff}) gives the radius of the Eddington capture sphere (Stahl et al. 2012, 2013), discussed below.
In the general case (k_{1} ≠ 0), the choice of initial conditions is less obvious. Since we are interested in the change of the motion of accretion disk corona, resulting from the change in luminosity, let us first consider circular orbits. The initial condition for these calculations corresponds to the motion of a test particle instantaneously in a circular orbit appropriate for stellar luminosity k_{1}. Effectively, we assume that in a steady luminosity field the system will adjust to the presence of radiation drag, resulting in nearly circular motion of the accreting fluid.
Imagine a particle executing a circular orbit around the star when it shines with luminosity k_{1}. This would only be possible if some external agency (e.g., the accretion disk) continually replenished the angular momentum and energy lost to radiation drag. The star is assumed to impulsively change its luminosity to k_{2}. We imagine the external agency to cease operation precisely at that instant when the spherical wavefront corresponding to the change in luminosity arrives at the orbit (for instance, the coupling between the disc and the corona is disrupted by the additional radiation pressure). Future hydrodynamic calculations will no doubt elucidate the true state of affair–ours is an exploratory calculation aiming to qualitatively identify the main features in the response of a corona to variable illumination by the central star.
Fig. 5 Solid curves: variation with the final luminosity of minimum escape radius, x = x_{esc}(k_{2}), after a sudden change to k_{2} from k_{1} = 0, or 0.36. The same curves also show the minimum final luminosity, k_{2} = k_{esc}(x) required to eject to infinity a particle previously orbiting at x. Stellar radius is X = 5.0. Dotted curve: variation of ECS radius with k_{2} (see text). 

Open with DEXTER 
Technically, we integrated Eqs. (1) and (2) with a value of k corresponding to the final luminosity, with the initial condition corresponding to a radius and velocity of a circularorbit solution of Eqs. (7) and (8) at luminosity k_{1}. The particle then follows a trajectory defined by its initial velocity and position, as well as the three stellar parameters M, R, k_{2}.
We find that at any particular radius, x, there is a minimum luminosity k_{2} = k_{esc} that will cause a particle satisfying the initial conditions to escape to infinity. A final luminosity value larger than k_{esc} will also lead to the expulsion of the particle. The k_{esc}(x) curve is monotonically decreasing. Figure 5 illustrates this behavior. The inverse function, x_{esc}(k_{2}) = x(k_{esc}) is the minimum radius of circular orbits from which the particles escape (minimum escape radius), for a given value of k_{2}. For a given value of k_{2}, the sphere at x_{esc} now divides space into an outer region, x ≥ x_{esc}, from which the particles escape, and an inner region, x<x_{esc}, the particles from which remain bound to the system (although they may leave the region x<x_{esc}).
Clearly, radiation drag is responsible for capturing all particles that are initially within the sphere of escape radius x_{esc}. At lower luminosities, the particles will actually be accreted onto the central star. However, at about Eddington luminosities they will lose all their momentum before settling on the star, and will levitate above the surface of the star in a state of equilibrium on a spherical surface concentric with the star, the socalled Eddington capture sphere (ECS) (Stahl et al. 2012, 2013). The ECS exists for luminosity k at a certain x satisfying X ≤ x< ∞, when the equation k_{eff}(x) = k has a solution (cf. footnotes 2 and 3).
This behaviour is illustrated in Fig. 6, where two trajectories are shown for motion starting with the same initial condition (instantaneously circular orbit at x = 8.0 under luminosity k_{1} = 0.36), but occurring for two different values of luminosities k_{2}, both satisfying k_{2}<k_{esc}(x), so that both motions are bounded. The two trajectories differ qualitatively in their termination point: at final luminosity k_{2} = 0.6 the particle is accreted with nonzero velocity onto the surface of the star, while at k_{2} = 0.888 the ECS is present, on which the nonescaping particles come to rest. The latter case is an example of particles being trapped outside their initial position, radiation drag causing here a net displacement outwards, somewhat paradoxically. The ECS luminosityradius relation is shown in Fig. 5. The quasiparadoxical net displacement outwards occurs only for radii to the left of the intersection of the ECS curve with the k_{esc}(x) curve, and only for values of k_{2} between the two curves. In Fig. 5 the minimum escape radius has been presented for simulations with two values of initial luminosity, k_{1} = 0.0 and 0.36. We also plotted variation of ECS radius with the final luminosity, k_{2}. For a fixed initial luminosity k_{1}, the escape sphere and the ECS coincide at a particular radius (Stahl et al. 2013). In Fig. 5 for initial luminosity k_{1} = 0.0 and 0.36, the minimum escape radius and the ECS radius are the same at x_{esc} = 8.412 and 14.365, and corresponding luminosities k_{2} = 0.873 and 0.928, respectively.
Figure 7 is helpful in illustrating the concepts of the minimum escape radius x_{esc}, and of the ECS. Here, we started with an initial luminosity of k_{1} = 0.36, impulsively changed to k_{2} = 0.930. Two trajectories are illustrated tangent to circles of radii x_{1} = 12.0 and x_{2} = 14.0. The initial velocities correspond to instantaneously circular motion under the influence of k_{1} = 0.36, as described above. In this example the minimum escape radius is x_{esc} = 13.0, and clearly x_{1}<x_{esc}, while x_{2}>x_{esc}. In the latter case [x(0) = x_{2}] the particle enters an escaping trajectory under the impulsive change of luminosity, while in the former [x(0) = x_{1}] it is captured by the ECS at x = 14.927.
Fig. 6 Two qualitatively different trajectories with the same initial condition for two different luminosities: k_{2} = 0.6 (thin blue trajectory) and k_{2} = 0.888 (violet trajectory). For k_{2} = 0.888 an ECS is present at x_{esc} = 9.44. The radius of star X = 5.0. The initial conditions are taken for x = 8.0 and k_{1} = 0.36 (cf. Fig. 5). 

Open with DEXTER 
Fig. 7 Trajectories corresponding to the k_{1} = L_{1}/L_{Edd} = 0.36 curve in Fig. 5. Final luminosity k_{2} = 0.930 is fixed in both trajectories. The stellar surface is shown at X = 5.0, and the ECS at x_{ECS} = 14.927 (magenta), the minimum escape radius is x_{esc} = 13.0 (middle circle in black). Bounded (solid violet line) trajectory starts at x = 12.0, unbounded (black dashed) trajectory starts at x = 14.0. 

Open with DEXTER 
Fig. 8 Trajectories corresponding to the k_{1} = 0.36 curve in Fig. 5 for L_{2}/L_{Edd} = 0.918. Black dashed trajectory is unbounded. The bounded trajectory (blue) ends on an ECS, which in this case is within the sphere of minimum escape radius (black circle). 

Open with DEXTER 
Two nonoverlapping trajectories for the case of nonzero initial luminosity (k_{1} = 0.36) are illustrated in Fig. 8. The final luminosity is taken to be k_{2} = 0.918, and now the ECS is inside the sphere of minimum escape radius. The black dashed trajectory illustrates an unbounded trajectory starting from an instantaneously circular orbit of radius exceeding x_{esc}. The test particle on the blue continuous trajectory is captured on the ECS. If the final luminosity is smaller than the one corresponding to an ECS at the stellar radius (at X = 5), the nonescaping particles (i.e., ones with x(0) <x_{esc}) will accrete on the star^{3}.
3.3. Luminosity effect on accreting particle
In previous section we discussed effects of sudden luminosity increase on a test particle which is initially moving in a circular orbit. In this section we describe a more realistic situation in which we suddenly increased the luminosity when the test particle is following a spiral trajectory. In this case we used as the initial condition a velocity corresponding to a trajectory initially tangent to a circular orbit. We chose initial radial distance, x = 20.0 at which the particle would be moving in a circular orbit if there were no radiation drag under luminosity k_{1} = 0.36. The particle will follow a spiral trajectory (Fig. 4), loosing its angular momentum, and finally either it will accrete on the surface of star or it will be captured on the ECS, depending on the luminosity.
We are interested in a test particle which is following a spiral trajectory under low luminosities, and while it does we suddenly increased the luminosity so that the particle can be ejected. Interestingly a higher luminosity is required to eject a test particle which is moving in a spiral trajectory at a particular radial distance than for a test particle moving in a circular orbit at the same distance. We can also say that for a particular final luminosity if test particles will be ejected from circular orbits for all x>x_{esc}, a test particle in a spiral trajectory will be ejected only if it is at a larger radial distance than x_{esc}, Fig. 10 (with the exception of x_{esc} = 6M).
The higher luminosity required to eject an accreting test particle is a consequence of a radial inward velocity component in case of spiral trajectory. From Fig. 10 it is clear that it requires a very much higher luminosity at outer regions for an infalling (“accreting”) test particle as compared to test particle initially following a circular orbit at the same radial distance. It is also clear from Fig. 9 that azimuthal velocity is lower for the particle following spiral trajectory as compared to particle following circular orbit. This will make particle even more tightly bound in the system when it is moving in a spiral trajectory. Thus it requires higher luminosity to eject the test particle. The outward radial four velocity component increases after sudden increase in luminosity and this causes the ejection of the test particle.
We know from GR that if relativistic specific energy of a geodesic particle e> 1, the particle is unbound and if e< 1, the particle is bound. GR specific energy is given by the covariant time component of four velocity, (9)If the GR specific energy for test particle in spiral trajectory is e_{s}(x) and for circular orbit e_{c}(x), then from Fig. 11 we see that 1 − e_{c}(x) < 1 − e_{s}(x). This shows how tightly the test particle is bound. In general, while higher luminosity is required to eject the test particle that is initially in a spiral orbit than would be the case for a circular orbit, this is not the case at the location of the ISCO (r = 6M), where Eddington luminosity is required (and suffices) to eject the particle, whether orbiting in the ISCO, or in a spiral trajectory intersecting the ISCO.
4. Discussion
Fig. 9 Variation of azimuthal component of four velocity, dφ/ dτ, with radial distance. From top to bottom: solid black curve corresponds to dφ/ dτ for the test particle in a sequence of circular orbits at different radii in the range, thick dashed (blue) curve, dφ/ dτ for a test particle following spiral trajectory for initial luminosity k_{1} = 0.36; red dashed dotted curve for k_{1} = 0.48; and magenta curve for k_{1} = 0.60. 

Open with DEXTER 
Fig. 10 Variation of the minimum luminosity required to eject the test particle decreases with radial distance. Dashed (blue) curve is showing the variation in case of an “accreting” test particle in a spiral trajectory and solid black curve is showing the variation for a circular orbit initial condition (Fig. 5). In both cases the initial luminosity k_{1} = 0.36 is fixed, (magenta) shortdashed curve is representing the ECS for final luminosity. 

Open with DEXTER 
Fig. 11 Relativistic specific energy as a function of radial distance for test particles in circular orbits and one in a spiral trajectory starting at an initial radial distance of x = 20.0. In both curves the luminosity is fixed at k_{1} = 0.36. In case of the circular orbits the drag terms are neglected, while the spiral trajectory is a solution of the full e.o.m. with the drag terms included. 

Open with DEXTER 
Radiation increases so strongly towards the source in general relativity that at first sight it seems easy to unbind matter orbiting a compact luminous star. For instance, a modest increase in luminosity would be sufficient to expel matter from the ISCO if there were no radiation drag, e.g., a change from 0 to 0.283 L_{Edd}, or from 0.360 L_{Edd} to 0.589 L_{Edd} would suffice, as can be seen in Fig. 2. However, the same radiation very strongly impedes the motion of matter moving in the optically thin regions illuminated by the star by exerting a drag that reduces the angular momentum and energy inherent in the motion. When this effect is properly taken into account, the conclusion is reversed: it is, in fact, very difficult to expel matter orbiting at the ISCO, typically this is only possible when superEddington luminosities are attained. However, the effects of drag fall off very quickly with distance to the radiation source, so that an Eddington outburst of a lower luminosity source (say, initial luminosity <0.3 L_{Edd}), such as an atoll source is sufficient to clear out all test particles orbiting in the optically thin region at r ≳ 10R_{G} (Figs. 1, 5)
The numerical results obtained in the paper apply to the motion of test particles. However, they should be also applicable, at least qualitatively, to optically thin plasma, e.g., the coronae of accretion disks. The radiation front following an outburst of the source moves at the speed of light appropriate for the medium, which is always larger than both the sound speed and the speed of the orbiting particles along any trajectory they may follow. All parts of the optically thin plasma will feel the influence of radiation before any significant hydrodynamic interaction occurs between the different regions at various radial distance from the source.
First let us consider plasma outside the spherical surface of minimum escape radius. Assuming that the plasma is distributed axisymmetrically, we can rule out intersection of outgoing trajectories from a ring of plasma originally located at a given radial distance from the source. Further, as the effects of drag fall off with distance, the motion of the outer rings of plasma will not be impeded by that of the inner rings. However, it is true that the inner rings will at first travel faster than the outer rings, because the outward velocity shortly after the passage of the radiation front scales with initial orbital velocity. This will lead to compressive heating of the plasma, thus robbing it of some kinetic energy and making the escape of the plasma less likely. The expected effect of replacing test particles with plasma is to increase the minimum luminosity change required to make the trajectories unbounded or, equivalently, to increase the minimum escape radius x_{esc} for a given final luminosity.
By the same arguments, plasma will remain bound to the system in those regions from which test particles cannot escape. Thus, we arrive at a robust conclusion that even in Eddington luminosity outbursts radiation drag makes it difficult for matter orbiting near the ISCO to escape the system. However, that region will be temporarily cleared of any optically thin plasma while the plasma follows the initial outwardsgoing part of its bounded trajectory.
As we have seen in Figs. 6–8, the bounded trajectories can cover a wide range of distance. It seems obvious that under the conditions corresponding to the noncircular bounded trajectories discussed in this paper, colliding shells of plasma will undergo compressive and dissipative heating. The temperatures involved are expected to be quite high, a fraction of the virial energy being involved in the transfer of kinetic to thermal energy (Kluźniak 2013). Thus, we would like to suggest that strong illumination alone of optically thin plasma may provide a sufficient heating mechanism to explain the origin of the observed high temperatures in the Xray emitting coronae of LMXBs.
5. Conclusions
We have numerically solved in full GR the equations of motion of a test particle moving in a strong radiation field. The motion of test particles, and most likely also of plasma, around a gravitating luminous body is strongly affected by rapid changes of its luminosity. For any initial luminosity and circular orbit, a sufficiently large increase in luminosity will unbind the particle. The effects of GR, and especially of radiation drag, make this “escape” luminosity a strong function of distance for the first several tens of gravitational radii. This GR result is in contrast with the Newtonian case, where the luminosity leading to ejection is independent of radius.
In particular, close to the ISCO radiation drag is so strong, that it takes an about Eddington luminosty outburst to eject particles to infinity. If the initial luminosity is a sizable fraction of Eddington the final luminosity required to eject the particle may be superEddington. This may explain why Xray bursters typically seem to be rebuilding their inner accretion disks shortly after the outburst. However, the effects of radiation drag drop very rapidly with distance, and already at radii ~ 20R_{G} a subEddington burst is typically sufficient to eject test particles to infinity.
We expect that similar conclusions will also hold for the motion of plasma, or hot gas, and we intend to verify these expectations in future work with a hydro code.
This criterion differs in some respects from the one adopted by Stahl et al. (2013).
The ECS radius ranges from R to infinity for luminosity in the range (Stahl et al. 2013).
Acknowledgments
We thank Wenfei Yu for illuminating conversations in the early stages of this project. Research supported in part by Polish NCN grants UMO34 2011/01/B/ST9/05439, 2013/08/A/ST9/00795.
References
 Abramowicz, M. A., Ellis, G. F. R., & Lanza, A. 1990, ApJ, 361, 470 [NASA ADS] [CrossRef] (In the text)
 Bini, D., Jantzen, R. T., & Stella, L. 2009, Class. Quant. Grav., 26, 055009 [NASA ADS] [CrossRef] (In the text)
 Bini, D., Geralico, A., Jantzen, R. T., Semerák, O., & Stella, L. 2011, Class. Quant. Grav., 28, 035008 [NASA ADS] [CrossRef] (In the text)
 Bradt, H. V. D., & McClintock, J. E. 1983, ARA&A, 21, 13 [NASA ADS] [CrossRef] (In the text)
 Kluźniak, W. 2013, A&A, 551, A70 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Kluzniak, W., & Wagoner, R. V. 1985, ApJ, 297, 548 [NASA ADS] [CrossRef] (In the text)
 Oh, J. S., Kim, H., & Lee, H. M. 2010, Phys. Rev. D, 81, 084005 [NASA ADS] [CrossRef] (In the text)
 Phinney, E. S. 1987, in Superluminal Radio Sources, eds. J. A. Zensus, & T. J. Pearson, 301 (In the text)
 Robertson, H. P. 1937, MNRAS, 97, 423 [NASA ADS] (In the text)
 Stahl, A., Wielgus, M., Abramowicz, M., Kluźniak, W., & Yu, W. 2012, A&A, 546, A54 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Stahl, A., Kluźniak, W., Wielgus, M., & Abramowicz, M. 2013, A&A, 555, A114 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Walker, M. A., & Meszaros, P. 1989, ApJ, 346, 844 [NASA ADS] [CrossRef] (In the text)
 Wielgus, M., Stahl, A., Abramowicz, M., & Kluźniak, W. 2012, A&A, 545, A123 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
Appendix A
Here, following Stahl et al. (2013), we present an abbreviated derivation of Eqs. (1), (2). The acceleration of the particle is in general given by (A.1)where ds^{2} = M^{2}dτ^{2}. Taking the force of radiation to be (A.2)we obtain the equation of motion (A.3)We determine the radiation flux F^{i} from the expression (A.4)where is the projection tensor, and T^{jk} is the stressenergy tensor for the radiation. Thus, the flux is (A.5)with (A.6)Expressing ϵ in tetrad components, we obtain Eq. (4). We take the tetrad components of the stressenergy tensor from Abramowicz et al. (1990): 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
All Figures
Fig. 1 Minimum luminosity that will eject to infinity a particle initially in circular orbit at x_{esc}. Plot compares solutions of the full e.o.m. (solid green curve, Eqs. (1), (2)) with solutions of the simplified equations that neglect the effects of radiation drag ([ND], thin red line, Eqs. (7), (8)). Initial luminosity is k_{1} = L_{1}/L_{Edd} = 0, stellar radius X = 5. The Newtonian limit is also shown ([NL], dotdashed line), as well as the radius of the Eddington capture sphere (ECS). See text for details. 

Open with DEXTER  
In the text 
Fig. 2 Radius of the escape sphere when drag is neglected. The minimum luminosity, k_{esc}, that unbinds a particle in circular orbit at x_{esc} is shown as a function of the initial orbital radius, for two values of initial luminosity, k_{1} = 0.0, and 0.36. Stellar radius is X = 5.0. Constant k_{esc} lines represent Newtonian limit (NL) for k_{1} = 0.0,0.36. If there were no drag, particles from all circular orbits of radii in the range 6 ≤ x(0) <x_{esc}(k_{esc}) would escape the system, and the initial radius in the outgoing trajectory, x(0), would be its periastron as well. 

Open with DEXTER  
In the text 
Fig. 3 Convergence of results using Eqs. (1), (2) and (6), (7), with the Newtonian limit. The ratio is plotted of the minimum luminosity sufficient to unbind a particle, initially in a circular GR orbit at k_{1} = 0, to the corresponding value in the Newtonian limit. Uppermost (thick green) curve corresponds to Eqs. (1) and (2), and red continuous curve to Eqs. (7) and (8). 

Open with DEXTER  
In the text 
Fig. 4 Testparticle trajectory beginning at r = 15R_{G} for a constant k = L_{2}/L_{Edd} = 0.048. Instead of maintaining a circular orbit the test particle follows a spiral trajectory and accretes on the surface of star. 

Open with DEXTER  
In the text 
Fig. 5 Solid curves: variation with the final luminosity of minimum escape radius, x = x_{esc}(k_{2}), after a sudden change to k_{2} from k_{1} = 0, or 0.36. The same curves also show the minimum final luminosity, k_{2} = k_{esc}(x) required to eject to infinity a particle previously orbiting at x. Stellar radius is X = 5.0. Dotted curve: variation of ECS radius with k_{2} (see text). 

Open with DEXTER  
In the text 
Fig. 6 Two qualitatively different trajectories with the same initial condition for two different luminosities: k_{2} = 0.6 (thin blue trajectory) and k_{2} = 0.888 (violet trajectory). For k_{2} = 0.888 an ECS is present at x_{esc} = 9.44. The radius of star X = 5.0. The initial conditions are taken for x = 8.0 and k_{1} = 0.36 (cf. Fig. 5). 

Open with DEXTER  
In the text 
Fig. 7 Trajectories corresponding to the k_{1} = L_{1}/L_{Edd} = 0.36 curve in Fig. 5. Final luminosity k_{2} = 0.930 is fixed in both trajectories. The stellar surface is shown at X = 5.0, and the ECS at x_{ECS} = 14.927 (magenta), the minimum escape radius is x_{esc} = 13.0 (middle circle in black). Bounded (solid violet line) trajectory starts at x = 12.0, unbounded (black dashed) trajectory starts at x = 14.0. 

Open with DEXTER  
In the text 
Fig. 8 Trajectories corresponding to the k_{1} = 0.36 curve in Fig. 5 for L_{2}/L_{Edd} = 0.918. Black dashed trajectory is unbounded. The bounded trajectory (blue) ends on an ECS, which in this case is within the sphere of minimum escape radius (black circle). 

Open with DEXTER  
In the text 
Fig. 9 Variation of azimuthal component of four velocity, dφ/ dτ, with radial distance. From top to bottom: solid black curve corresponds to dφ/ dτ for the test particle in a sequence of circular orbits at different radii in the range, thick dashed (blue) curve, dφ/ dτ for a test particle following spiral trajectory for initial luminosity k_{1} = 0.36; red dashed dotted curve for k_{1} = 0.48; and magenta curve for k_{1} = 0.60. 

Open with DEXTER  
In the text 
Fig. 10 Variation of the minimum luminosity required to eject the test particle decreases with radial distance. Dashed (blue) curve is showing the variation in case of an “accreting” test particle in a spiral trajectory and solid black curve is showing the variation for a circular orbit initial condition (Fig. 5). In both cases the initial luminosity k_{1} = 0.36 is fixed, (magenta) shortdashed curve is representing the ECS for final luminosity. 

Open with DEXTER  
In the text 
Fig. 11 Relativistic specific energy as a function of radial distance for test particles in circular orbits and one in a spiral trajectory starting at an initial radial distance of x = 20.0. In both curves the luminosity is fixed at k_{1} = 0.36. In case of the circular orbits the drag terms are neglected, while the spiral trajectory is a solution of the full e.o.m. with the drag terms included. 

Open with DEXTER  
In the text 