Accreting raindrop with quadratic drag Created 2026-10-14 Updated 2026-10-18
For a spherical drop of constant mass density , radius and downward velocity , suppose its volume grows at rate by collecting stationary mist and its upward quadratic drag is , with . The drop mass is . Its incoming mass accretion rate and the variable-mass system momentum law give the displayed equations. The accreted matter must be accelerated from rest, producing the additional term alongside drag.
Accretion column Created 2026-10-07 Updated 2026-11-05
An accretion column carries an inward accretion flow onto a localized region of a star. In the strong-field slender approximation, the magnetic field controls the flux-tube area while longitudinal pressure and gravity determine acceleration. The signed outward-coordinate velocity is negative; the inward speed and mass accretion rate are usually defined positive.
Angular-momentum flux Created 2026-10-05 Updated 2026-10-14
The rate at which angular momentum crosses a surface, including advection and stresses. With positive inward mass accretion rate, the net inward flux in a steady accretion disk is .
Black-hole accretion Created 2026-09-24 Updated 2026-11-05
Black-hole accretion is the inflow of matter onto a black hole, converting some of its gravitational binding energy into heat, electromagnetic radiation and outflows. Matter with substantial angular momentum generally forms an accretion disk; internal stresses transfer angular momentum outward so that matter can move inward. The released energy can be radiated before the material crosses the event horizon. In a thin-disk model the innermost stable circular orbit marks the transition toward rapid plunging motion, although the detailed inner flow depends on the accretion state.[ref]
The relation defines the radiative efficiency of black-hole accretion, where is emitted luminosity, the mass accretion rate and the speed of light. A radiatively inefficient accretion flow carries much of its available energy inward or into outflows instead of emitting it, so luminosity alone need not measure the mass supply.
Two integrations of the steady Keplerian viscous diffusion equation give this general transported quantity. fixes the mass accretion rate, while fixes the additive angular-momentum flux or inner torque. A zero-torque inner boundary condition gives ; a nonaccreting constant-torque disk has .
Keplerian accretion disk Created 2026-09-24 Updated 2026-11-05
A Keplerian accretion disk orbits a dominant central mass with angular velocity and specific angular momentum . In an axisymmetric steady thin disk without wind losses, take the inward mass accretion rate to be constant. Conservation of angular momentum gives
The zero-torque inner boundary condition at fixes this constant to . Since Keplerian rotation gives , the viscous torque in an accretion disk is . Equating the two expressions yields
Far outside the inner edge this approaches . Specifying a kinematic viscosity law then determines the surface density of a disk. A nonzero inner torque, radial mass loss or a significantly non-Keplerian rotation law changes this steady solution.[ref]
For a steady, spherically symmetric radial inflow, conservation of mass makes the mass accretion rate independent of radius:
Atomic hydrogen gives . The Newtonian free-fall speed from rest at infinity is , obtained by setting the specific mechanical energy to zero. Hence the spherical free-fall accretion rate is
Using the given values gives , , and
The area factor assumes a full spherical inflow with negligible angular momentum. With only a local number density and speed, an arbitrary nonspherical inflow's total rate is not fixed; one instead integrates over the inflow surface. This is a flux estimate, rather than the pressure-selected Bondi accretion solution.
Following a parcel. Let be a fluid pathline, with . The multivariable chain rule gives the material derivative,
The first term measures change at fixed position; the advection term accounts for the parcel moving through spatial gradients. The derivative is also called the Lagrangian derivative.
The forces on a parcel. A material volume has mass . The pressure forces on its boundary sum to , while the gravitational acceleration is . Newton's second law therefore gives the Euler equations for an inviscid fluid,
This neglects viscosity and other body forces. The left side is the parcel's acceleration, not just the time derivative at a fixed spatial location.
Bernoulli conservation. The claimed streamline result needs steady flow and a barotropic fluid, , or an equivalent enthalpy relation along each isentropic streamline. These hypotheses are necessary; the unsteady Euler equation alone does not imply the stated conservation law. Define the barotropic enthalpy function
Dot the steady Euler equation with , using . This gives
Thus the Bernoulli function is
Every term has units of energy per mass, . Different streamlines can have different constants.
The spherical accretion equation. Write the inward speed as , so the velocity vector is and the positive mass accretion rate is
Assume a prescribed central potential and negligible gas self-gravity. The steady radial Euler equation is
For a barotropic equation of state, , with . Eliminate the density derivative to obtain
A smooth sonic point requires the right side to vanish where , hence
For a nonisothermal gas this is an implicit relation involving the local sound speed. In the isothermal case is constant and this is the Bondi sonic point.
The isothermal integral. For , choose the enthalpy reference so that . The reservoir is at rest at infinity and , making its Bernoulli constant zero. Thus
or
The logarithm has a dimensionless argument. Constant temperature during compression generally requires heat removal, so this is Isothermal Bondi accretion, rather than adiabatic accretion with .
The critical accretion rate. Select the smooth reservoir-connected transonic branch. Evaluate the last relation at and :
Mass conservation then gives the isothermal Bondi accretion rate,
The units are . Equivalently , with .
The inner-boundary selection of isothermal Bondi flow is important: a standard absorbing accretor inside the sonic radius selects this critical solution. The steady equations and alone also admit subcritical subsonic branches supported by a different inner boundary; they do not by themselves fix a unique rate. The usual answer presumes transonic capture.
Figure 1.
Original numerical transonic isothermal accretion profiles, normalized by the sonic radius and outer density; the speed becomes supersonic inside the sonic point
. Details.
Take for inward radial velocity, and assume a constant specific-heat ratio . The isentropic flow has and adiabatic sound speed . Steady spherical mass conservation and the radial Euler equation for fluid motion give
Eliminating the mass density derivative gives , or
Differentiating gives the corresponding sound speed equation:
At a smooth sonic point, . A finite derivative requires the numerator of the flow equation to vanish too. Thus
This regularity requirement concerns a smooth transonic branch, rather than an arbitrary singular solution reaching unit Mach number.
Integrating the momentum equation using the specific enthalpy yields the Bernoulli equation
For a finite positive mass accretion rate, its large-radius solution has
These expressions satisfy both the flow equations asymptotically and show and . In particular agrees with the first boxed equation. Combining the Bernoulli equation with sonic regularity additionally gives , so a finite positive sonic point occurs in the range identified below.
For the small-radius inward free fall branch, neglecting enthalpy against in the Bernoulli equation predicts . The continuity equation then fixes , and consequently . The neglected-to-leading ratio is
It tends to zero precisely when
This makes the assumed supersonic flow and leading free fall self-consistent. At equality, the ratio is of constant order rather than vanishing; above it, the proposed asymptotic balance fails. The isothermal limit also permits free fall, with logarithmic rather than power-law enthalpy.
For the stellar funnel, write the magnetic scalar potential as , where and . Its magnetic dipole field has and , so a magnetic-field-line equation is
The edge of the polar flux tube meets the stellar surface at . Thus, near the pole, , and its transverse radius is . This proves the dipolar flux-tube area law
Equivalently magnetic flux conservation gives and . The stipulated small cross section keeps the tube close to the polar axis, where its direction is approximately radial; the approximation is not valid all the way to arbitrarily large .
Outside the star the prescribed field is a gradient, hence and . Its Lorentz force density is therefore . The magnetic-pressure and magnetic-tension parts cancel in this exact current-free dipole. Field-aligned velocity also has , consistent with a steady prescribed field. The following calculation uses the narrow polar funnel approximation.
Now steady mass conservation reads . With , replace the spherical factor 2 by 3:
The Bernoulli equation is unchanged because the magnetic field supplies no work or force. On a gravity-dominated free fall branch, , but and . Hence
The same bound follows from regular transonic accretion in a power-law tube: and the reservoir-matched Bernoulli equation gives . The stellar surface cuts off the mathematical small- limit; the criterion refers to a funnel that becomes gravity dominated and crosses its sonic point before reaching that surface.
For an equatorial circular orbit in the Paczyński-Wiita potential, radial force balance gives
Its specific angular momentum is . Perturbing the circular radius at fixed gives the radial epicyclic frequency squared . Therefore
Circular motion is radially stable for , marginal at , and unstable between and . Thus the inner edge of the nearly circular thin accretion disk is , joining a plunging region of a black-hole accretion disk inside the innermost stable circular orbit. This conclusion concerns the circular thin-disk model, not an assertion that all gas disappears inside this radius.
Take the inward mass accretion rate as positive, and write the positive outward viscous torque in an accretion disk as
Steady conservation of angular momentum makes the net inward flux independent of radius. Impose a zero-torque inner boundary condition at , appropriate to an approximately stress-free transition to plunging motion. Then
This is the steady accretion disk with arbitrary rotation law angular-momentum balance.
The local viscous dissipation integrated through the full disk thickness is . It is the power per unit planar area summed over both radiating faces. Since
substitution gives the two-face dissipation of a Paczyński-Wiita disk:
In particular vanishes at the stress-free inner edge and is nonnegative outside it because increases there. A one-face radiative flux is ; the total accretion luminosity is therefore , with no additional factor of two.
For a direct evaluation put . Then
The first integral is , using . For the second,
so it equals . Hence
This is the radiative efficiency of black-hole accretion within this pseudo-Newtonian model.
Independently, the specific orbital energy of a circular particle is
Thus and . The energy released per unit accreted mass is exactly , confirming the integrated accretion luminosity. The Paczyński-Wiita circular orbit reproduces the Schwarzschild inner radius but gives an approximate binding energy; the efficiency derived here is for the stated potential.
Write the inward radial velocity as , with . Steady spherical mass conservation and the radial Euler equations for an inviscid fluid give
For a barotropic fluid, define the sound speed and the barotropic enthalpy function by . Choosing , integration gives the Bernoulli equation
The continuity equation gives . Eliminating mass density from the momentum equation produces the barotropic spherical accretion equation
The sonic point is where the inward Mach number equals one. A smooth solution cannot cross the vanishing coefficient on the left with a nonzero right-hand side, so a regular sonic radius must obey
Conversely a continuous flow that is subsonic at infinity and supersonic inward must cross unit Mach number; smooth passage selects a simultaneous zero of the two sides, hence a regular sonic point. The associated critical accretion branch, rather than an arbitrary subsonic solution with smaller flux, is selected by these outer and inner requirements.
For the polytropic equation of state, set to distinguish its exponent from any separately specified perturbation exponent. Assume and first take the usual positive index . Then
For the attractive mixed potential, assume . The sonic condition and Bernoulli equation give
Hence polytropic accretion in a mixed inverse-power potential reduces to
For all nonzero coefficients in the left side are nonnegative for , with at least one strictly positive: there is no finite positive root. For the constant term is negative and the leading coefficient positive, giving exactly one positive root. Define and . Then
This is the finite-radius transonic range for the two positive attractive terms.
The critical point has real distinct crossing slopes. To check this rather than just solve the sonic equalities, put , and . Differentiating the barotropic spherical accretion equation at its simultaneous zero gives
Its discriminant is . The sonic Bernoulli equation says , so this discriminant is strictly positive for . The accretion crossing is the slope for which , namely the smaller root. The two slopes have opposite signs of this Mach-number derivative. This is the regular transonic saddle, not a singular vertical tangent.
The corresponding mass accretion rate follows without an additional arbitrary mass density normalization:
One may equivalently replace by . With a nonzero term, inner free fall has and continuity gives . Thus increases without bound precisely in the positive-index range , consistent with the required inner supersonic branch.
The assumptions on the parameters matter. If and , the finite sonic radius instead requires , the usual Bondi accretion range; is the zero-radius limiting critical case. If and , the condition remains . With arbitrary signs of , one must use the quadratic and retain only positive roots with and real admissible crossing slopes; attraction was not explicitly given as a sign condition in the source.
If positive polytropic index was not intended as an implicit physical convention, there is a further mathematical range. For , positive compressibility also permits , giving . Both potential terms attractive then give a positive sonic root as well:
The same mass accretion rate formula applies; its crossing-slope discriminant is positive. The range has and no physical acoustic sonic point, has constant pressure and zero sound speed, and does not define the printed equation of state. Thus the usual answer assumes and positive attractive coefficients; these conventions should not be confused with additional printed restrictions.
Put . In steady state, is constant, so
The inward mass accretion rate and outward viscous torque in an accretion disk are
Use the positive mass accretion rate , and retain , . Define , which is positive under the condition in (b). Matching the Bernoulli function to the reservoir and using the polytropic equation of state gives
The transonic spherical accretion rate in a power-law potential is therefore
The combination has dimensions of length to the power , so this expression has dimensions of mass per time. The sonic point selects the flux that connects the subsonic reservoir to the inward supersonic transonic branch.
For the endpoint limits of power-law spherical accretion, hold fixed. As , , so . Hence
This is also obtained directly from an isothermal equation of state: the Bernoulli function becomes , so . For , , it recovers the isothermal Bondi accretion rate.
For , and means . Since ,
Thus
The finite limiting flux accompanies and ; it does not assert a finite-radius sonic point at the endpoint. For this is the familiar limit .
For completeness, the printed range has the endpoint . At , and the flux is independent of , namely . For , and , so as . These are the corresponding extended endpoint limits.
For the finite-mass similarity solution above, the total disk mass is
Substitute , so and . This gives
The requested linear increase of mass is therefore correct. It is supplied by the infinite-radius boundary, rather than by a source-free isolated finite disk. The inward mass accretion rate is
which tends to zero at the origin and to at infinity.
For Keplerian rotation, the specific angular momentum is . The total disk angular momentum would be
This integral diverges at , corresponding to infinite radius. There is no finite constant total angular momentum for the nonzero printed similarity solution. The seemingly time-independent last expression is an infinite quantity; it cannot prove a finite conservation law.
The problem remains visible under a cutoff. With , write , where and is the exponential integral. Then
Even the finite change after subtracting a common logarithmic divergence is not zero. Direct integration of the disk diffusion equation gives the boundary identity
for this profile the outer expression has a nonzero limit, giving the same finite change rate. An inner zero torque alone does not impose zero outer angular-momentum flux.
More generally, the printed family can have finite angular momentum at both ends only if and , hence it is trivial. Thus the request for a nonzero finite constant is incompatible with the stated viscosity exponent, similarity ansatz and infinite radial domain. A physical finite-domain or altered similarity problem would need different boundary data or assumptions; silently assigning a finite value to the divergent integral does not repair it.
The specific orbital energy released in a black-hole accretion flow reaching an inner radius is of order . For a circular Keplerian orbit, the binding energy is more precisely , with the other half of the potential-energy decrease retained as orbital kinetic energy. Multiplying by the mass accretion rate gives the accretion luminosity scale
For an inner radius of order the Schwarzschild radius, , this becomes
This is an energy-release scaling, not unit radiative efficiency of black-hole accretion. The actual inner orbit, stresses and advected energy determine the coefficient and the emitted luminosity; the innermost stable circular orbit of a Schwarzschild black hole lies outside its horizon.
Take to mean the magnitude of the inward radial velocity. In the Gaussian units used for the magnetic pressure, approximate the surface density of a disk by . The mass accretion rate and free-fall speed then give
At the magnetospheric truncation radius, balance the magnetic pressure of the magnetic dipole field with the specified ram pressure:
Thus
Discarding the convention-dependent numerical factor gives
The magnetic pressure grows inward as , faster than the ram pressure , so the field disrupts the flow inside this radius. The numerical coefficient depends on the vertical averaging and on treating the inflow as radial free fall, as stipulated here; a detailed thin disk boundary model need not have the same coefficient.
The magnetic torque transfers angular momentum from the faster rotator to the slower one. When , the neutron star supplies angular momentum to the inner accretion disk and possible propeller flow, so it tends to spin down. When the inner accretion disk rotates faster than the star, the magnetic connection and accreted angular momentum tend to spin it up.
For slowly varying mass accretion rate and magnetic dipole moment, these tendencies favour a magnetic spin equilibrium with the corotation radius near the magnetospheric truncation radius:
The precise balance depends on accreted and ejected angular momentum and on the extent of magnetic coupling across the accretion disk; corotation is a scaling estimate, not an exact zero-torque theorem.
Use for inward mass accretion rate, and for magnetic torque adding angular momentum to the disk. Combining the continuity equation with conservation of angular momentum identifies the inward mass flux as . In a steady Keplerian accretion disk, this is the supplied constant , so
The boundary condition is . Since , the viscous torque in an accretion disk is . Hence the steady Keplerian accretion disk with an inner torque has
In particular, . For a physical steady solution with this boundary condition, .
For , the surface density of a disk rises from a small inner value toward : this is the usual Keplerian accretion disk with an approximately zero-torque inner boundary condition. For , the inner surface density of a disk is large and decreases approximately as , giving a torque-dominated accretion disk. This approximation applies where ; sufficiently far out, any fixed finite again approaches the same constant asymptote.
Figure 1.
Steady surface-density profiles with weak and strong inner torques
. Details.
The constant net inward angular-momentum flux is
For , net angular momentum moves outward. The strong-torque profile is therefore decretion disk-like in its stress-dominated structure, but the specified mass flux remains inward: a genuine net decretion disk requires an outward mass flux, not merely large .
At , the Maxwell stress tensor estimate gives
Using the magnetospheric truncation radius scaling in part (b), all dimensional parameters cancel:
If the displayed is retained and the scaling radius is assigned unit coefficient, the estimate is ; retaining and exact free-fall speed throughout gives . Those prefactors are not controlled by the scaling argument. In particular, it does not generate a parametrically large .
The transonic accretion in a dipolar flux tube values are
The mass density follows from for the same polytropic equation of state. If is the total loaded area at the surface, mass conservation at the sonic point gives
For two equal polar caps with surface angular radius , , so the total mass accretion rate is
For only one loaded polar cap, divide this expression by two. Expressing the result first in also covers a bundle whose loading angle is defined at another radius. The smooth isothermal limit has , , and .
Take to be the inward radial speed and select the regular critical solution for a gas reservoir of density and adiabatic sound speed . Other mathematically subsonic steady solutions need a different inner boundary condition. For the polytropic equation of state, let , so
Mass conservation and the Bernoulli equation give
The gas is at rest at infinity, so its kinetic and gravitational terms vanish there.
Differentiate these two conserved quantities, using . Eliminating the density derivative yields the polytropic Bondi accretion equation
At a regular sonic point, both sides vanish:
The Bernoulli equation at that point reduces to . For , equivalently ,
Since , the critical density is . Evaluate the conserved mass accretion rate at the sonic point:
where
The isothermal limit gives , agreeing with Isothermal Bondi accretion.
For , and the Newtonian critical radius tends to zero while diverges. The limiting critical solution approaches Mach number one only as ; it has no finite-radius Newtonian sonic transition. This is a failure of the Newtonian inner approximation, not an obstruction to accretion onto a black hole. A physical Schwarzschild black hole has a finite event horizon, and the regular relativistic inflow becomes supersonic outside it.
The singular-looking factor in the rate has a finite limit. Writing , its logarithm is , so the critical Bondi accretion rate for gamma equals five thirds is
This is also the leading cold-reservoir limit of relativistic spherical accretion; an exact rate for relativistically hot gas requires the relativistic equations, not this Newtonian expression.
For completeness, the finite relativistic sonic point can be seen without assigning it an arbitrary inner radius. Let denote the relativistic squared sound speed and . For an isentropic relativistic gas with , its dimensionless specific enthalpy is . Conservation of rest-mass flux and relativistic Bernoulli energy gives the critical relations
For and , expansion gives , hence
Here the asymptotic relativistic and Newtonian sound speeds agree to leading order. This finite-radius transition explains why the gas can reach the horizon despite the Newtonian result.