Preprint
Article

A New Perspective to the Relation Between Accretion Luminosity and Mass Accretion Rate

This version is not peer-reviewed.

Submitted:

15 July 2024

Posted:

16 July 2024

You are already at the latest version

Abstract
In this paper I show that, opposite to what is generally assumed, when a source accretes matter at the so-called Eddington mass accretion rate M˙ed, the generated accretion luminosity Lacc is smaller than Led, the Eddington accretion luminosity that can stop accretion. I demonstrate my thesis using pieces of informations that, one by one, are well known in the accretion theory but that, for some reason, they are not considered altogether. Namely, I use the following facts: i) Lacc, in the ideal case, is equal to the rate at which kinetic energy is deposited into the accreting source; ii) radiation pressure decreases the velocity of the infalling matter; iii) the generally accepted formula Lacc=GMM˙/R is valid only if the matter is in free fall. To give a quantitative demonstration of my thesis I use two accretion models: the first is the very simple spherical accretion of zero-temperature, fully ionized hydrogen; then I use a model developed some years ago where the energy transfer between matter and radiation is taken into account. The fact that M˙ed generates a Lacc<Led can have consequences in all the cases in which one would need to model accretion with M˙>M˙ed to reproduce the observations.
Keywords: 
;  ;  

1. Introduction

Accretion is an important phenomenon in astrophysics and can be the only source of luminosity for those objects that do not have an internal energy production, like neutron stars or black holes. It is also important during the formation of a star because in this case too the accreting source can irradiate only the energy, or a fraction of it, due to the infalling matter that releases its kinetic energy when its fall stops.
The amount of kinetic energy that can be radiated away depends strongly on how the motion is stopped and on the efficiency of converting this energy in luminosity. In any case, energy conservations requires that if all the kinetic energy is converted in accretion luminosity, L acc in the following, we have
L acc ( k ) = 1 2 v f 2 M ˙
M ˙ being the mass accretion rate and v f the final velocity of the infalling matter at R, the radial position where the energy conversion happens; ( k ) stands for kinetic because this definition of L acc is based on the velocity of the particles. Note that, throughout this paper, L acc ( k ) or L acc alone, is the accretion luminosity at R.
Equation (1) is seldom, if ever, found in papers on accretion because L acc is, essentially always, written in this way
L acc = G M M ˙ R
with M mass and R radius, in the sense described above, of the central source. When matter is in free fall, then v f 2 = v ff 2 = 2 G M / R , with v ff final free-fall velocity, and the two above equations are equivalent: the potential energy per unit mass, G M / R , is entirely converted in kinetic energy.
The accretion luminosity exerts a force [1] that alters the motion of the matter reducing its velocity. The question arises whether the radiation pressure can ultimately halt accretion, and for which M ˙ this happens. In analogy with what derived by Eddington [2] who found that the radiation pressure due to thermonuclear reactions in the interior of a star can make the stellar atmosphere evaporate1, the L acc that stops accretion is named L ed . For spherical accretion of ionized hydrogen, it is easy to derive [4]
L ed = 4 π c m σ · G M = 3.3 × 10 4 M M L
where c is the velocity of light, m is the proton mass and σ is the Thomson cross section for electrons. If we assume that Equation (2) holds for any M ˙ , we can define the Eddington mass accretion rate, M ˙ ed , as
M ˙ ed = 4 π c m σ · R = 1.1 × 10 3 R R M yrs 1 .
The fact that L ed can be reached with a finite value of M ˙ has important implications in astrophysics, for instance in topics like high-mass star formation [e.g., [5], accretion onto compact objects or growth of super massive black holes [e.g., [6]. In these cases, the observational results sometimes require to model the accretion with high rates, such that M ˙ > M ˙ ed . Indeed, many works have been developed to overcome the limitation on M ˙ , see, e.g., Donnan et al. [7] and references therein.
Now, the aim of this paper is to show that Equation (2) is not valid when radiation pressure is present. Clearly, G M / R measures the depth of the gravitational potential well but the effect of the radiation pressure is to reduce the velocity of the matter. It would make no sense, otherwise, to look for the value of M ˙ that stops accretion. Thus, v f < v ff , or v f < 2 G M / R , while, as said, Equation (2) holds only if v f = v ff . In turn, this means that L acc ( k ) < L acc . The inequality might make the reader think that energy is not conserved; on the contrary, in the next sections I will show that assuming Equation (2) violates the energy conservation because it does not take into account that by decelerating the matter, radiation makes a work. The force acting on each particle is not only the accelerating gravitational force, but also the decelerating radiative force. The latter term is also due to gravitation so, in some sense, gravity is responsible of both forces.
The idea that gravitational force seems generated by a mass smaller than M, because of the radiation pressure, is not new: for instance, it can be found in Taam et al. [8] who introduced an effective mass, smaller than M, to take into account the radiative deceleration; however, they overlooked the consequence of this force reduction. If M ˙ , and then the number of falling particles, increases more work is necessary to slow down the matter. In turn, this means that more gravitational energy, in the form of radiation, is used to decelerate the matter, and the fraction of this energy that is deposited at R is reduced. It is then clear that the radiation pressure acts as a negative feedback, it eats kinetic energy which can be seen as a reduction of gravitational potential energy. As a consequence, a complete accretion stop can happen, if ever, only for M ˙ > M ˙ ed .
This paper is so organized: In Section 2 the two models I used are presented and discussed. In Section 3 the results of the second model, more realistic than the first one, in terms of velocity of the particles and accretion luminosity, are presented. Finally, in Section 4, the results are discussed and conclusions are given.

2. Methods

To demonstrate my thesis, I will exploit two accretion models: the first is very simple and for this reason analytical formulae can be derived, but at the price of being a bit unphysical. It is, however, often used to give a quick demonstration that L acc = L ed when M ˙ = M ˙ ed [4,9]. I will use it to show easily the consequences of adopting Equation (1) instead of Equation (2). The second model is more accurate because it takes into account the energy transfer between matter and radiation.

2.1. A Simple Spherical Accretion Model

This model is unphysical because it considers the interaction between matter and radiation but only in dynamical terms, ignoring that the radiation is absorbed by the infalling matter. As written, its utility comes from the fact that the equation of motion can be solved analytically.
A particle of mass m falling from a radius r experiences two forces: i) the gravitational force G M m / r 2 down to the surface radius R that we suppose much smaller than r ; and 2) the radiation pressure due to Thomson scattering. In the classical treatment of radiation [1], the radiation momentum is completely absorbed by matter at a rate W σ where W is the energy density of the electromagnetic field and σ is the Thomson cross section. Since σ depends on the inverse of the square of the particle mass, this force acts mainly on the electrons. On the contrary, the gravitational force is larger for protons: as a consequence the flow should be treated as a two fluids plasma but assuming that we have only fully ionized hydrogen, we can treat the problem as one single fluid. In fact, the electric field arising from the interactions between protons and electrons keeps the two components coupled [10]. The force acting on each particle is then σ I / c where I = c W is the intensity field and c is the light speed.
Depending on how we define the accretion luminosity, the equation of motion can be written in two different ways:
F tot = σ I c G M m r 2 = { σ L acc 4 π c r 2 G M m r 2 = σ G M M . 4 π c r 2 R G M m r 2 ( 5 ) σ L acc ( k ) 4 π c r 2 G M m r 2 = σ v f 2 M . 8 π c r 2 R G M m r 2 ( 6 )
where m is the proton mass and, when ignoring absorption, we can set I ( r ) = L acc / 4 π r 2 .
Equation (4) follows immediately from Equation (5) if the condition F tot = 0 is imposed. However, it should be noted that F tot = 0 means that the velocity is constant, not necessarily that v = 0 everywhere, so that we have to verify that v = 0 is indeed an allowed solution.
The problem is that if v = 0 , clearly L acc = 0 as well: the condition M ˙ = M ˙ ed , through Equation (2), implies L acc = L ed and, at same time, implies that L acc = 0 . Thus, the combination of Equation (5) and Equation (2) does not admit the solution v = 0 . The inconsistency comes from what has been already said: Equation (2) is valid only for matter in free fall and this is not the case. The only conclusion we can draw from setting F tot = 0 in Equation (5) is that v = const 0 .
If we adopt Equation (1) and solve Equation (6), it is immediate to derive
1 2 m v 2 = 1 r 1 2 σ M ˙ 4 π c G M m
from which
v ( r ) = 2 G M r R v f 2 μ r
with μ = M ˙ / M ˙ ed . Thus
v f 2 = 2 G M R 1 1 + μ = v ff 2 1 1 + μ
where v ff is the final velocity for matter in free-fall: as expected, v f = v ff when μ = 0 . The solution v = 0 is acceptable asymptotically, it is approached as μ .
By combining Equation (1) with Equation (9), after some algebra one derives
L acc ( k ) = L ed μ 1 + μ
showing that L ed , like v = 0 , is asymptotically reached for μ .
In this way we can see that increasing M ˙ increases L acc ( k ) as well. But taking into account the negative feedback of the radiation pressure, that is feeded by the gravitational potential well, the increase in luminosity does not scale linearly with M ˙ .

2.1.1. The Case L 0

In Equation (6), the luminosity produced internally by the accreting source, named in the following L , has been assumed to be zero. In case L 0 , the equations change but the final result is, qualitatively, not different from the previous case, provided that L < L ed .
Equation (6) becomes
d v 2 d r = σ 4 π c m r 2 v f 2 M ˙ + 2 σ L 4 π c m r 2 2 G M r 2 = R μ v f 2 r 2 + 2 G M r 2 Λ 2 G M r 2
with Λ = L / L ed ; Equation (8) is now
v ( r ) = 2 G M r R v f 2 μ r 2 G M r Λ
and
v f 2 = v ff 2 1 Λ 1 + μ ,
so that
v ( r ) = 2 G M r 1 Λ 1 + μ .
Thus, accretion is still possible with any M ˙ provided that Λ < 1 . When Λ = 1 , the radiation pressure due to L that clearly does not depend on M ˙ , stops accretion.
The accretion luminosity is easily derived:
L acc ( k ) = 1 2 v f 2 M ˙ = L ed μ ( 1 Λ ) 1 + μ
with L acc = 0 for Λ = 1 .

2.2. A More Realistic Accretion Model

In the previous section we have seen the consequences of defining the accretion luminosity following Equation (1) or Equation (2), by using a simple model. This model is based on three assumptions: i) zero-temperature flow, so that no interaction between particles is taken into account; ii) no energy exchange between matter and radiation; iii) L ( r ) decreases with radius only by the geometrical factor 4 π r 2 , which implies a constant otpical depth τ = 0 at all radii.
Now, assumption i) makes the model not too much realistic even if not unphysical. But assumptions ii) and iii) are not consistent with the fact that the interaction between radiation and matter must happen in order to justify the equation of motion given either by Equation (5) or by Equation (6).
Now, assumption i) makes the model not too much realistic even if not unphysical. But assumptions ii) and iii) are not consistent with the fact that the interaction between radiation and matter must happen in order to justify the equation of motion given either by Equation (5) or by Equation (6).
In this section, thus, I will employ a more realistic accretion model. In literature one can find accurate accretion models that describe the physics of the accretion process, e.g., Fukue [11] or Murakami et al. [12], just to make a few of examples. Adapting these models to my needs, however, is not simple because in many cases the use of Equation (2) is not written explicitly, so that it is not immediate to find out how to modify a given model to adopt Equation (1). And, as I am going to show, it can happen that Equation (2) and Equation (1) are unintentionally adopted at the same time.
For these reasons, I use a model developed fifty years ago by Maraschi et al. [13], MTR in the following, that while considering the exchange between matter and radiation, so assumption ii) above is no longer present, it is still quite simple.
The set of equations used by MRT is:
E r = 1 ρ r U 3 ( Euler s equation ) ,
4 π r 2 ρ v = M ˙ ( continuity equation ) ,
r ( 4 π r 2 F ) = M ˙ E r ( energy balance between radiation and gas ) ,
F = c m σ 1 ρ r U 3 + 4 3 U v ( transfer equation for photons ) .
where the multiplying 4/3 factor in front of U v was added by Kafka and Mészáros [14]; note also that I have assumed a negative velocity, so the sign in the continuity equation and in the term 4 / 3 U v is opposite to MRT’s set.
Equation (18) gives
F = M ˙ 4 π r 2 ( E + ψ ) = M ˙ 4 π r 2 E + G M R
where the constant ψ is found assuming that particles’ energy goes to zero as r . Now, if we analyze carefully Equation (20), we will see that it implies Equation (1) even if MTR assumed explicitly that Equation (2) holds. In fact, Energy E is given by
E = 1 2 v 2 G M r
while the accretion luminosity is L = 4 π r 2 F ; for r = R we have
L acc = M ˙ E ( R ) + G M R = M ˙ 1 2 v f 2 G M R + G M R = 1 2 M ˙ v f 2 = L acc ( k )
so that, in MTR’s model there is the hidden adoption of Equation (1). Of course, if one assumes that v f = v ff , that is, that the final velocity is equal to the free-fall final velocity, the two formulations of L acc give the same result. But MRT derived for their model the velocity profile v ( r ) , finding that for any μ 0 , it is always v f < v ff , which is expected. Then we have
L acc ( k ) = 1 2 M ˙ v f 2 < 1 2 M ˙ v ff 2 = G M M ˙ R
which shows that MTR’s model was not internally consistent. Of course, my aim is not to show the inconsistency of MTR’s work, I make these considerations only to explain why I adopt a simple model instead of a complex and much more realistic, accretion model. Only in this way I can go deep in the physics of the solution to read behind the lines what are all the assumptions, intentionally or not, of the model.
Now, for r we have
L M ˙ 1 2 v 2 + G M R G M M ˙ R
where the equal sign holds when v 0 for r . By comparing Equation (23) with Equation (24), we have to conclude either that two equations are not consistent, or that L ( r ) increase with r, which looks unphysical. The reason of the inconsistency comes from neglecting the absorption of the radiation, but this problem becomes evident only if Equation (1) is assumed, because with Equation (2) it is clear that L G M M ˙ R .
The set of Equations (17) — (19) cannot be solved under the assumption iii) (i.e., τ ( r ) = 0 r , see the beginning of this section). We have to take into account the absorption of the radiation by writing
F = M ˙ 4 π r 2 E + G M R e τ
with
τ = σ m R r ρ d r = σ m M ˙ 4 π R r 1 v r 2 d r , τ r = σ m M ˙ 4 π 1 v r 2
( τ / r will be used later on in this section).
By combining the different equations, it is possible to arrive to a second order differential equation in v: all the steps are detailed in Appendix A where the general case L 0 is treated. Here I present the final equation for the case L = 0 :
2 V z 2 + 1 z 2 μ 4 1 2 e τ 1 2 V 2 μ e τ 2 V 2 1 z 1 V z + 1 z 3 V 2 μ z 1 + μ z 3 e τ c μ 2 v ff z + V 3 2 z V + μ 2 z 3 V e τ c μ v ff z 2 V c μ v ff z V 2 = 0
with z = r / R , V = v R 2 G M and, for typical neutron stars,
c v ff 1.840 R 10 km 1 2 M M 1 2
so that Equation (26) becomes
τ = c v ff μ 1 z 1 V z 2 d z 1.840 μ R 10 km 1 2 M M 1 2 1 z 1 V z 2 d z
with τ > 0 because V < 0 .
When μ = 0 Equation (27) turns into
2 V z 2 1 2 z 2 V 2 V z 1 z 3 V = 0
and it is easy to verify that the free-fall velocity V = 1 / z is a solution.
Let’s now see if a solution V = cos t , say V nf (no force), is possible. Equation (27) becomes
2 μ z 1 + V nf μ e τ c μ 2 v ff z + V nf 3 2 z V nf μ 2 e τ c μ v ff z 2 V nf c μ v ff z V nf 2 = 0
A constant solution, by definition, must be valid for any z, in particular also for z = 1 where the previous equation simplifies in
2 μ 1 + V nf μ c μ 2 v ff + V nf 3 2 V nf μ c μ v ff V nf 1 = 0
from which
V nf 3 + c μ 2 v ff V nf 2 + 3 μ 2 2 μ V nf + c μ v ff = 0 .
It is convenient to plot this cubic equation to look for the solutions: this is done in Figure 1 for a set of values of μ . The case μ = 0 is not shown because it corresponds to matter in free fall and V cannot be constant; for μ = 0.01 , a curve not shown in the figure, a solution for V nf is not found because the curve, in the range shown, never intersects the line y = 0 . A value of V nf that is solution of the cubic equation starts appearing somewhere in the range 0.01 < μ 0.05 . But for μ 2 , the range of values explored in the next section, Equation (27) admits only V ( z ) < 0 z , so that V nf > 0 is not acceptable: as a consequence, Equation (27) does not admit a solution V = const and, in turn, this means that the total force acting on each particle can never be zero.

3. Results

Equation (27) has been solved using the solver Radau of the routine solve_ivp, part of the Python package scipy [15], using an iterative approach2: first, Equation (8), that in adimensional variables becomes
V = 1 z ( 1 + μ ) ,
is used in Equation (29) to derive a first estimate of the optical depth, namely
τ = 2 c v ff μ 1 + μ 1 z 1 .
With the above τ ( z ) , Equation (27) is solved with boundary conditions
V = 1 z ( 1 + μ ) V z = 1 2 1 1 + μ 1 z 1.5
where z = 100
Once V ( z ) and V / z are derived, τ ( z ) is computed again but now directly from Equation (29). The new curve τ ( z ) is used in Equation (27) and the procedure is iterated. Given two iterations, i and i + 1 , the iterations is stopped once δ i + 1 = | ( τ ( z ) i + 1 τ ( z ) i ) | / τ ( z ) i + 1 < 10 3 .
This approach works for μ < 1 . When μ 1 , δ i decreases until reaches a minimum, then it oscillates between two values, one small and one large. In these cases, I stop the procedure once δ i + 1 > δ i , and the solution corresponding to iteration i is used as solution.
I have also tried another approach in which, given the new solution of V ( z ) , Equation (27) is solved starting from z = 1 and going to larger z. The optical depth is computed at runtime through
τ 0 = 0 τ i = τ i 1 + c v ff μ 1 z i 1 2 V i 1 + 1 z i 2 V i z i z i 1 / 2
This approach failed because sometime the routine solve_ivp calls the integrating function twice with the same3z, but the second time the first derivative of V jumps to high values and V ( z ) diverges.
Five values of μ are chosen for this work: 0.1, i.e., one tenth of M ˙ ed , 0.5, 1.0, 1.5 and 2.0. In Table 1, for each μ < 1 , column Iteration gives how many iterations are necessary to reach convergence ( δ i + 1 < 10 3 ); for the other values of μ , the two last iterations are reported, once δ starts oscillating between two values. In these case, the iterations are stopped when δ i + 1 > δ i and the adopted τ ( z ) is the one corresponding at the smallest δ : for instance, when μ = 1 the last decreasing value is for iteration 9, and the adopted τ is 4.0739.
In Figure 2, the curves τ ( z ) are shown. Because of the lack of convergence of τ , it is possible that τ is a bit overestimated for μ 1 .
The solutions V ( z ) are shown in the left panel of Figure 3: at large distances they are similar, approaching the asymptotic law given by Equation (34). At short radii, z < 2 , when μ 1 a maximum in | V | starts appearing, clearly visible in the right panel of Figure 3 where a zoom-in of the left panel for small z is shown.
This means that infalling particles are accelerated up to very close to the surface of the accreting object, and then, at z 2 , they are strongly decelerated. In a complete thermodynamical treatment, this behaviour would be modelled with a schock, causing a strong increase in temperature. The presence of a shock was already foreseen by Shapiro and Salpeter [16], and even at similar distances; it is, however, surprising that the possibility of a shock is suggested even with a simple, zero-temperature approach as in this case. Clearly, lacking a proper thermodynamical set of equations, we can only note that the strong deceleration looks like a shock, but nothing more can be said.
Finally, in Figure 4 the luminosity Λ is shown, where
Λ ( z ) = μ V 2 1 z + 1 e τ
As expected, Λ ( z ) follows the trend of V ( z ) modulated by the exponential of τ . As said before, however, it is likely that the optical depth is overestimated at large z, so that, especially for μ > 1 , the luminosity might be underestimated.
Figure 4 gives a quantitative representation of what has been said at the end of the Introduction: for high M ˙ , the work done by the radiation pressure to decelerate the infall matter makes the emerging luminosity be more and more faint. This is because even if the number of particles increases with M ˙ , the kinetic energy per particle is more and more reduced.

4. Discussion

We have seen that Equation (20) implies Equation (1), so it could seem that the aim of this paper, to promote the use of L acc ( k ) instead of the more common L acc , is not new: even if not explicitly written, we can find that equation already in a fifty-years-old model. Moreover, Equation (2) is often not explicitly written in modern accretion models, as if it were not used.
However, Equation (2) is at the base of the following equivalence, used in pratically all the papers on accretion:
μ M ˙ M ˙ ed = L acc L ed Λ .
No matter how L acc is written, the above equivalence clearly states that when μ = 1 then Λ = 1 , and once L acc = L ed accretion is stopped.
On the contrary, if one follows Equation (1), then
μ M ˙ M ˙ ed > L acc ( k ) L ed Λ ( k ) ;
as a consequence, not only μ = 1 does not imply that the accretion luminosity is equal to L ed , also, it is to be demonstrated that L ed can ever be reached with a finite M ˙ .
I have derived Equation (40) first using a basic accretion model that has the only advantage to give an analytical expression for the velocity law v ( r ) .
Once v ( r ) is known, L acc ( k ) can be computed and one can see that L ed is an asymptotic value reached only for M ˙ . In this model, when μ = 1 we have Λ ( k ) = 0.5 .
I have then adopted a second accretion model developed by MTR. This model is more realistic than the first one, but with a strong assumption: the falling particles do not interact each other during the motion. This is equivalent to say that the flow is at zero-temperature. Such an assumption is, of course, less than optimal and the results obtained are, with no doubt, of limited value. However, it is important to stress that this model is not unphysical, no physical law is violated, it is only not too much realistic. It can be seen as an ideal experiment of thermodynamics that cannot be executed in real world but that can give indications on how a system evolves. The model used here can give insights on what happens in the accretion flow when Equation (1) is used.
I understand that the reader can feel uncomfortable if someone claims that what has been assumed in the last sixty years, Equation (2) or Equation (39), is not correct. But my conclusion is based on Equation (1): first, it is based on energy conservation; second, it was already used in the past, as shown before. Equation (2) violates the energy conservation because it does not take into account that the radiation pressure makes work on the infalling matter, changing the equation of motion. As a consequnce, also the energy balance must be adapted. It is clear that the depth of the gravitational potential well remains G M / R , independent on M ˙ . On the other hand, the radiation pressure acts as a counter-reaction whose amplitude is a fraction of G M / R , and one has to demostrate that this fraction can indeed be 1.
All the accretion models show that the matter is not in free fall, being decelerated by the radiation pressure, thus Equation (2) can no longer be valid and Equation (40) holds. For any developed accretion model, it has to be demonstrated that radiation pressure can halt accretion with a finite M ˙ , this assumption must be proved and, in general, for μ = 1 L ed is not reached.

Funding

This research received no external funding.

Data Availability Statement

All the python routine used in this work are freely available at https://drive.google.com/drive/folders/0ADGdDSs8TC5WUk9PVA.

Acknowledgments

Preliminary versions of this paper were read by Prof. Francesco Strafella and Dr. Patrick Hennebelle, who I kindly thank for their suggestions. I wish also to thank Ms. Anne Zhang, Managing Editor of Universe, for being patient with all the delays in the submission of the work.

Conflicts of Interest

The author declares no conflicts of interest.

Appendix A

In this section, the system of Equations (16)–(19) will be used to derive the second order differential equation in V, Equation (27), in the more general case of L , the internal luminosity of the accreting source, different from zero. By setting in the following L = 0 , all the equations given in Section 2.2 can be recovered.
First, Equation (25) becomes
F = M ˙ 4 π r 2 E + G M R e τ + L 4 π r 2 e τ
with τ given by Equation (26).
Then, I follow MTR’s procedure: first, combining Equation (A1) with Equation (19) we have
M ˙ 4 π r 2 E + G M R e τ + L 4 π r 2 e τ = c m σ 1 ρ r U 3 + 4 3 U v ,
and using Equation (16)
M ˙ 4 π r 2 E + G M R e τ + L 4 π r 2 e τ = c m σ E r + 4 3 U v
from which
U = 3 4 M ˙ 4 π r 2 v E + G M R e τ + 3 4 L 4 π r 2 v e τ 3 4 c m σ v E r .
To eliminate U, Equation (A4) is derived with respect to r and Equation (16) is used again. But before that, I write E as a function of v, differently from MTR who wrote v as a function of E,
E = 1 2 v 2 G M r , E r = v v r + G M r 2
so that
U r = 3 ρ E r = 3 M ˙ 4 π r 2 v E r = 3 M ˙ 4 π r 2 v r + 3 M ˙ 4 π r 2 v G M r 2 ,
3 4 M ˙ 4 π r E e τ r 2 v = 3 4 M ˙ 4 π r 2 e τ 1 2 v r + 1 2 σ M ˙ 4 π m r 2 v r + G M r v 2 v r + G M r v σ M ˙ 4 π m r 2 v + 3 G M r 2 v
3 4 G M M ˙ 4 π R r e τ r 2 v = 3 4 M ˙ 4 π G M R e τ r 2 v σ M ˙ 4 π m r 2 v 1 v v r 2 r
3 4 L 4 π r e τ r 2 v = 3 4 L 4 π e τ r 2 v σ M ˙ 4 π m r 2 v 1 v v r 2 r
and
3 4 c m σ r 1 v E r = 3 4 c m σ 2 v r 2 G M r 2 v 2 v r 2 G M r 3 v .
Thus
3 M ˙ 4 π r 2 v r + 3 M ˙ 4 π r 2 v G M r 2 = = 3 4 M ˙ 4 π r 2 e τ 1 2 v r + 1 2 σ M ˙ 4 π m r 2 v r + G M r v 2 v r G M r v σ M ˙ 4 π m r 2 v + 3 G M r 2 v + + 3 4 M ˙ 4 π G M R e τ r 2 v σ M ˙ 4 π m r 2 v 1 v v r 2 r + 3 4 L 4 π e τ r 2 v σ M ˙ 4 π m r 2 v 1 v v r 2 r 3 4 c m σ 2 v r 2 G M r 2 v 2 v r 2 G M r 3 v ,
and after expanding all the multiplications
3 4 c m σ 2 v r 2 + 3 M ˙ 4 π r 2 v r + 3 M ˙ 4 π r 2 v G M r 2 3 4 c m σ G M r 2 v 2 v r 3 4 c m σ 2 G M r 3 v 3 4 M ˙ 4 π r 2 e τ 1 2 v r + 3 4 M ˙ 4 π r 2 e τ 1 2 σ M ˙ 4 π m r 2 + 3 4 M ˙ 4 π r 2 e τ v r 3 4 M ˙ 4 π r 2 e τ G M r v 2 v r + 3 4 M ˙ 4 π r 2 e τ G M r v σ M ˙ 4 π m r 2 v + 3 4 M ˙ 4 π r 2 e τ 3 G M r 2 v 3 4 M ˙ 4 π G M R e τ r 2 v σ M ˙ 4 π m r 2 v + 3 4 M ˙ 4 π G M R e τ r 2 v 1 v v r + 3 4 M ˙ 4 π G M R e τ r 2 v 2 r + 3 4 L 4 π e τ r 2 v σ M ˙ 4 π m r 2 v + 3 4 L 4 π e τ r 2 v 1 v v r + 3 4 L 4 π e τ r 2 v 2 r = 0 .
After simplifying the term 4 σ / 3 c m we get
2 v r 2 + 4 σ c m M ˙ 4 π r 2 v r G M r 2 v 2 v r σ c m M ˙ 4 π r 2 e τ 1 2 v r σ c m M ˙ 4 π r 2 e τ G M r v 2 v r + σ c m M ˙ 4 π G M R e τ r 2 v 1 v v r + + σ c m L 4 π e τ r 2 v 1 v v r + 4 σ c m M ˙ 4 π r 2 v G M r 2 2 G M r 3 v σ c m M ˙ 4 π r 2 e τ 1 2 σ M ˙ 4 π m r 2 + σ c m M ˙ 4 π r 2 e τ v r + + σ c m M ˙ 4 π r 2 e τ G M r v σ M ˙ 4 π m r 2 v σ c m M ˙ 4 π r 2 e τ 3 G M r 2 v σ c m M ˙ 4 π G M R e τ r 2 v σ M ˙ 4 π m r 2 v + + σ c m M ˙ 4 π G M R e τ r 2 v 2 r σ c m L 4 π e τ r 2 v σ M ˙ 4 π m r 2 v + σ c m L 4 π e τ r 2 v 1 v v r + σ c m L 4 π e τ r 2 v 2 r = 0
that can be cast in this form
2 v r 2 + 4 σ M ˙ 4 π c m r 2 G M r 2 v 2 σ M ˙ 4 π c m r 2 e τ 1 2 σ M ˙ 4 π c m r 2 e τ G M r v 2 + σ M ˙ 4 π c m G M R e τ r 2 v 1 v + + σ c m L 4 π e τ r 2 v 1 v v r + 4 σ M ˙ 4 π c m r 2 v G M r 2 2 G M v r 3 σ M ˙ 4 π c m r 2 e τ 1 2 σ M ˙ 4 π m r 2 + + σ M ˙ 4 π c m r 2 e τ v r + σ M ˙ 4 π c m r 2 e τ G M r v σ M ˙ 4 π m r 2 v σ M ˙ 4 π c m r 2 e τ 3 G M r 2 v + σ M ˙ 4 π c m G M R e τ r 2 v σ M ˙ 4 π m r 2 v + σ M ˙ 4 π c m G M R e τ r 2 v 2 r σ c m L 4 π e τ r 2 v σ M ˙ 4 π m r 2 v + σ c m L 4 π e τ r 2 v 2 r = 0 .
Now we use the relation
σ M ˙ 4 π c m = R μ
so that
2 v r 2 + 4 R μ r 2 G M r 2 v 2 R μ r 2 e τ 1 2 R μ r 2 e τ G M r v 2 + μ G M e τ r 2 v 1 v + σ c m L 4 π e τ r 2 v 1 v v r + + 4 R μ r 2 v G M r 2 2 G M v r 3 R μ r 2 e τ 1 2 R c μ r 2 + R μ r 2 e τ v r R μ r 2 e τ G M r v R c μ r 2 v + R μ r 2 e τ 3 G M r 2 v μ G M e τ r 2 v R c μ r 2 v + μ G M e τ r 2 v 2 r σ c m L 4 π e τ r 2 v R c μ r 2 v + σ c m L 4 π e τ r 2 v 2 r = 0 ,
and finally
2 v r 2 + 1 r R μ r 4 1 2 e τ G M r v 2 μ G M r e τ v 2 R r 1 + σ c m L 4 π e τ r v 2 v r + + 2 r 2 v G M r 2 R μ r 1 + R μ r 3 e τ 1 2 R c μ r + v 3 G M r v + R μ r 2 v e τ G M r R c μ r 2 v + μ c r v 2 R σ c m L 4 π e τ r 3 v R c μ r v 2 = 0 .
It might be of some interest to note that the equation for E derived by MTR contains the term ( E / r ) 2 , while in the above equation a quadratic term in v / r is not present.
At this point, the following dimensionless variables are introduced: z = r / R , V = v R 2 G M , obtaining
2 v r 2 = 2 G M R 1 R 2 2 V z 2 ; 1 r v r = 2 G M R 1 R 2 1 z V z
R μ r 4 1 2 e τ = μ z 4 1 2 e τ ,
G M r v 2 = 1 2 z V 2 ,
μ G M r e τ v 2 R r 1 = μ e τ 2 z V 2 1 z 1 ,
σ c m L 4 π e τ r v 2 = L L ed e τ 2 z V 2 = Λ e τ 2 z V 2
with Λ = L / L ed
2 r 2 v G M r 2 R μ r 1 = 1 R 2 z 3 V 2 G M R 2 μ z 1 ,
R μ r 3 e τ 1 2 R c μ r + v 3 G M r v = μ R 2 z 3 e τ 2 G M R c μ 2 v ff z + V 3 2 z V
where
v ff = 2 G M R
is the final free-fall velocity;
R μ r 2 v e τ G M r R c μ r 2 v + μ c r v 2 R = μ 2 R 2 z 3 V 2 G M R e τ c μ v ff z 2 V c μ v ff z V 2 ,
and, finally,
σ c m L 4 π e τ r 3 v R c μ r v 2 = 2 G M R Λ e τ 2 R 2 z 3 V c μ z V v ff 2 .
All the terms can be collected and the common factor 2 G M / R × 1 / R 2 can be cancelled out, obtaining
2 V z 2 + 1 z 2 μ 4 1 2 e τ 1 2 V 2 μ e τ 2 V 2 1 z 1 + Λ e τ 2 V 2 V z + 1 z 3 V 2 μ z 1 + μ z 3 e τ c μ 2 v ff z + V 3 2 z V + μ 2 z 3 V e τ c μ v ff z 2 V c μ v ff z V 2 Λ e τ 2 z 3 V c μ z V v ff 2 = 0 .
Equation (27) corresponds to Equation (A28) once Λ is set to zero. The procedure to compute the solution V ( z ) is the same outlined in Section 3: the initial τ ( z ) is derived assuming the velocity law given by Equation (14) that, in the new variables, becomes
V = 1 Λ z ( 1 + μ )
Thus, Equation (35) is now
τ = 2 c v ff μ 1 + μ 1 + Λ 1 z 1 .
Equation (A28) has been resolved for two values of μ , 0.1 and 1.0, and two values of Λ , 0.1 and 0.5. If Λ 1 , V changes sign so that Equation (17) and Equation (19) must be changed as well. However, the case of reverted motion is out of the scope of this paper and will not discussed.
In Figure A1 the solutions of V ( Z ) are reported for the four cases studied. The two panels are quite similar to those in Figure 3. The main difference is that Λ does not depend on μ , nor on V f 2 , so that, at the surface of the accreting object, the radiation pressure is higher than in the case of Λ = 0 causing a stronger deceleration.
Figure A1. Left panel: velocity V, in units of free-fall velocity 2 G M / R vs. normalized radial distance for two values of M ˙ , in units of M ˙ ed , and two values of L , in units of L ed . Right panel: zoom-in of the left figure to emphasize the maximum in | V | at short radii.
Figure A1. Left panel: velocity V, in units of free-fall velocity 2 G M / R vs. normalized radial distance for two values of M ˙ , in units of M ˙ ed , and two values of L , in units of L ed . Right panel: zoom-in of the left figure to emphasize the maximum in | V | at short radii.
Preprints 112270 g0a1
The luminosity profile, given by
Λ ( z ) = μ V 2 1 z + 1 + Λ e τ ,
is shown in the two panels of Figure A2. For μ = 1.0 and Λ = 0.5 , the optical depth increases so fast that the luminosity goes to very small values just above the surface of the accreting object. However, as written in Section 2.2, the procedure used to solve Equation (A28) causes a likely overestimated of τ ( r ) : a more rigorous treatment of accretion, without the zero-temperature assumption, would probably give a more gentle transition in V when z approaches 1, and, consequently, a less severe attenuation of Λ at high μ . In any case, at least qualitatively, we can conclude that μ = 1 does not stop accretion.
Figure A2. Left panel: luminosity Λ , in units of L ed , vs. normalized radial distance for different values of M ˙ in units of M ˙ ed . Right panel: zoom-in of the left figure at short radii.
Figure A2. Left panel: luminosity Λ , in units of L ed , vs. normalized radial distance for different values of M ˙ in units of M ˙ ed . Right panel: zoom-in of the left figure at short radii.
Preprints 112270 g0a2

References

  1. Landau, L.D.; Lifshitz, E.M. The classical theory of fields; 1971.
  2. Eddington, A.S. On the radiative equilibrium of the stars. Mon. Not. R. Astron. Soc. 1916, 77, 16–35. [Google Scholar] [CrossRef]
  3. Abramowicz, M.A. Super-Eddington black hole accretion: Polish doughnuts and slim disks. Growing Black Holes: Accretion in a Cosmological Context; Merloni, A.; Nayakshin, S.; Sunyaev, R.A., Eds., 2005, pp. 257–273, [arXiv:astro-ph/astro-ph/0411185]. [CrossRef]
  4. Frank, J.; King, A.R.; Raine, D.J. Accretion power in astrophysics; 1985.
  5. Zinnecker, H.; Yorke, H.W. Toward Understanding Massive Star Formation. Ann. Review of Astron. & Astroph. 2007, arXiv:astro-ph/0707.1279]45, 481–563. [Google Scholar] [CrossRef]
  6. Inayoshi, K.; Visbal, E.; Haiman, Z. The Assembly of the First Massive Black Holes. Ann. Review of Astron. & Astroph. 2020, arXiv:astro-ph.GA/1911.05791]58, 27–97. [Google Scholar] [CrossRef]
  7. Donnan, F.R.; Hernández Santisteban, J.V.; Horne, K.; Hu, C.; Du, P.; Li, Y.R.; Xiao, M.; Ho, L.C.; Aceituno, J.; Wang, J.M.; Guo, W.J.; Yang, S.; Jiang, B.W.; Yao, Z.H. Testing super-eddington accretion on to a supermassive black hole: reverberation mapping of PG 1119+120. Mon. Not. R. Astron. Soc. 2023, arXiv:astro-ph.GA/2302.09370]523, 545–567. [Google Scholar] [CrossRef]
  8. Taam, R.E.; Fu, A.; Fryxell, B.A. Accretion in Wind-driven X-Ray Sources. The Astroph. Journal 1991, 371, 696. [Google Scholar] [CrossRef]
  9. Eckart, A.; Schödel, R.; Straubmeier, C. The black hole at the center of the Milky Way; 2005.
  10. Maraschi, L.; Reina, C.; Treves, A. On Spherical Accretion near the Eddington Luminosity. Astron. & Astroph. 1974, 35, 389. [Google Scholar]
  11. Fukue, J. Bondi Accretion onto a Luminous Object. Pub. of the Astroph. Soc. of Japan 2001, 53, 687–692. [Google Scholar] [CrossRef]
  12. Murakami, M.; Nishihara, K.; Hanawa, T. Self-Similar Gravitational Collapse of Radiatively Cooling Spheres. The Astroph. Journal 2004, 607, 879–889. [Google Scholar] [CrossRef]
  13. Maraschi, L.; Reina, C.; Treves, A. On Spherical Accretion near the Eddington Luminosity. Astron. & Astroph. 1974, 35, 389. [Google Scholar]
  14. Kafka, P.; Mészáros, P. How fast can a black hole eat? General Relativity and Gravitation 1976, 7, 841–846. [Google Scholar] [CrossRef]
  15. Virtanen, P.; Gommers, R.; Oliphant, T.E.; Haberland, M.; Reddy, T.; Cournapeau, D.; Burovski, E.; Peterson, P.; Weckesser, W.; Bright, J.; van der Walt, S.J.; Brett, M.; Wilson, J.; Millman, K.J.; Mayorov, N.; Nelson, A.R.J.; Jones, E.; Kern, R.; Larson, E.; Carey, C.J.; Polat, İ.; Feng, Y.; Moore, E.W.; VanderPlas, J.; Laxalde, D.; Perktold, J.; Cimrman, R.; Henriksen, I.; Quintero, E.A.; Harris, C.R.; Archibald, A.M.; Ribeiro, A.H.; Pedregosa, F.; van Mulbregt, P.; SciPy 1. 0 Contributors. SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python. Nature Methods 2020, 17, 261–272. [Google Scholar] [CrossRef] [PubMed]
  16. Shapiro, S.L.; Salpeter, E.E. Accretion onto neutron stars under adiabatic shock conditions. The Astroph. Journal 1975, 198, 671–682. [Google Scholar] [CrossRef]
1
On the validity of this attribution, see Abramowicz [3].
2
The routine code is available for download, see paragraph Data Availability Statement toward the end of the paper.
3
It can also happen that z i < z i 1 , but this case is handled easily.
Figure 1. The intersections of the straight y = 0 (black line) with the set of curves, drawn with different μ , give the value V nf that solve Equation (33). As explained in the text, an acceptable solution for our problem requires that V nf intersects the line y = 0 for negative values. For μ 2.0 this does not happen: Equation (27) does not admit a solution V = const .
Figure 1. The intersections of the straight y = 0 (black line) with the set of curves, drawn with different μ , give the value V nf that solve Equation (33). As explained in the text, an acceptable solution for our problem requires that V nf intersects the line y = 0 for negative values. For μ 2.0 this does not happen: Equation (27) does not admit a solution V = const .
Preprints 112270 g001
Figure 2. Optical depth τ vs. normalized radial distance for different values of M ˙ in units of M ˙ ed .
Figure 2. Optical depth τ vs. normalized radial distance for different values of M ˙ in units of M ˙ ed .
Preprints 112270 g002
Figure 3. Left panel: velocity V, in units of free-fall velocity 2 G M / R vs. normalized radial distance for different values of M ˙ in units of M ˙ ed . Right panel: zoom-in of the left figure to emphasize the maximum in | V | at short radii.
Figure 3. Left panel: velocity V, in units of free-fall velocity 2 G M / R vs. normalized radial distance for different values of M ˙ in units of M ˙ ed . Right panel: zoom-in of the left figure to emphasize the maximum in | V | at short radii.
Preprints 112270 g003
Figure 4. Left panel: luminosity Λ , in units of L ed , vs. normalized radial distance for different values of M ˙ in units of M ˙ ed . Right panel: zoom-in of the left figure at short radii.
Figure 4. Left panel: luminosity Λ , in units of L ed , vs. normalized radial distance for different values of M ˙ in units of M ˙ ed . Right panel: zoom-in of the left figure at short radii.
Preprints 112270 g004
Table 1. Here I report, for each value of μ , the number of iterations; τ i + 1 and τ i , the optical depth at the largest distance from the accreting object, where i corresponds to the number of iterations; the relative difference δ i + 1 of the optical deptth at infinity, between two successive iterations.
Table 1. Here I report, for each value of μ , the number of iterations; τ i + 1 and τ i , the optical depth at the largest distance from the accreting object, where i corresponds to the number of iterations; the relative difference δ i + 1 of the optical deptth at infinity, between two successive iterations.
μ Iteration τ i + 1 τ i δ i + 1
0.1 2 0.3518 0.3518 8.741 · 10 5
0.5 4 1.9420 1.9419 9.982 · 10 6
1.0 9 4.0648 4.0739 2.240 · 10 3
10 4.0748 4.0648 2.446 · 10 3
1.5 2 7.1799 6.1722 1.403 · 10 1
3 5.7619 7.1799 2.461 · 10 1
2.0 1 16.175 11.471 2.908 · 10 1
2 7.6956 16.175 1.102
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

Downloads

66

Views

19

Comments

0

Subscription

Notify me about updates to this article or when a peer-reviewed version is published.

Email

Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2025 MDPI (Basel, Switzerland) unless otherwise stated