Bondi accretion rate calculator

The Bondi accretion rate calculator determines the amount of gas captured in free fall by a compact object (black hole, neutron star, white dwarf, protostar) immersed in a uniform gaseous medium per the formulation by Austrian cosmologist Hermann Bondi (1919-2005) published 1952 in Monthly Notices of the Royal Astronomical Society "On spherically symmetrical accretion" at King's College London. The classic Bondi-Hoyle-Lyttleton formula: Ṁ_Bondi = 4π · λ_c · G² · M² · ρ_∞ / c_s³, with G gravitational constant (6.674×10⁻¹¹), M object mass kg, ρ_∞ ambient gas density kg/m³, c_s sound speed m/s, λ_c polytropic factor (γ = 5/3 → λ = 0.25; γ = 4/3 → λ = 0.71; γ = 1 isothermal → λ ≈ 1.12). The Bondi radius: r_B = G · M / c_s² (gravitational sphere of influence).

Free fall velocity at Bondi radius: v_ff = √(2GM/r_B) = √2 · c_s. Numerically for 1 M_sun in 10 K molecular cloud (c_s = 0.2 km/s) ρ = 10⁻²⁰ kg/m³: r_B = 7×10¹⁵ m = 50,000 AU; Ṁ ≈ 10⁻⁶ M_sun/year (protostellar formation).

For supermassive AGN BHs M = 10⁹ M_sun in 10⁷ K intracluster medium (c_s = 500 km/s) ρ = 10⁻²⁵ kg/m³: r_B = 50 kpc; Ṁ ≈ 10⁻³ M_sun/year.

Relative velocity of object vs medium (Bondi-Hoyle-Lyttleton 1944) modifies formula: Ṁ_BHL = 4π G² M² ρ / (v_rel² + c_s²)^{3/2}.

Italian astrophysical Bondi applications:

• Sagittarius A* Milky Way INAF Bologna Loiano 1936 (BH 4×10⁶ M_sun, ρ = 10⁻²⁰ kg/m³ Bondi r_B = 0.1 pc, Ṁ_Bondi = 10⁻⁶ M_sun/year → observed 10⁻⁹ M_sun/year → "starved" radiative feedback 1,000×)
• microquasar Cygnus X-1 INAF Brera Milan 1762 (BH 21 M_sun accreting from O9.7 Iab supergiant via 2,000 km/s stellar wind)
• X-ray binaries Centaurus X-3, Hercules X-1, Cygnus X-3 (Riccardo Giacconi Uhuru 1971, Roche-lobe overflow accretion)
• AGN M87 EHT Event Horizon Telescope 2019 first image (BH 6.5×10⁹ M_sun, relativistic jet 5,000 ly, Ṁ_Bondi ≈ 10⁻³ M_sun/year)
• protostars T-Tauri ρ Ophiuchi INAF stellar formation (Ṁ ≈ 10⁻⁶ M_sun/year)
• Bondi-Hoyle-Lyttleton "wind accretion" symbiotic stars INAF Padua 1761. Italian accretion theory:

• Remo Ruffini ICRA International Center Relativistic Astrophysics Rome 1985 Sapienza University ("blackholic energy" Bondi-Eddington combined optimization, GRB BdHN Binary-driven Hypernova model)
• Riccardo Giacconi Genoa 1931-2018 Nobel Physics 2002
• Roberto Maiolino Cambridge ex-INAF Arcetri Florence 1872 first galaxies JWST + AGN feedback. Andrea Ghez UCLA + Reinhard Genzel MPE 2020 Nobel Prize "discovery of supermassive compact object center Milky Way" Sgr A* tracking S2 stars 16 years measured at Keck and VLT (Italian ESO INAF 25% participation)

Italian Bondi astrophysics tables

(Source: Bondi 1952 MNRAS 112, 195 + Bondi-Hoyle-Lyttleton 1944 + Frank, King, Raine "Accretion Power" Cambridge 2002 + INAF Rome 2002).

Italian Bondi calculation examples

• Sagittarius A* Milky Way 8.3 kpc INAF Bologna Loiano 1936: M_BH = 4.0×10⁶ M_sun = 8×10³⁶ kg. Interstellar gas ρ_∞ ≈ 10⁻²⁰ kg/m³ (n_H = 100/cm³), c_s = 500 km/s. r_B = G·M/c_s² = 2.1×10¹⁵ m = 0.07 pc. Ṁ_Bondi = 10⁻⁶ M_sun/year theoretical. Observed L_X ≈ 10²⁶ W → Ṁ_real ≈ 10⁻⁹ M_sun/year (1,000× sub-Bondi → "starved", ADAF feedback Narayan-Yi 1994)
• Microquasar Cygnus X-1 1,860 pc INAF Brera Milan 1762: BH 21 M_sun + O9.7 Iab supergiant. Stellar wind v_w = 2,000 km/s, ρ_w ≈ 10⁻⁹ kg/m³. Bondi-Hoyle-Lyttleton: Ṁ_BHL ≈ 10⁻⁷ M_sun/year, L_X = 0.1 c² Ṁ ≈ 6×10²⁹ W
• M87 SMBH Virgo 16.4 Mpc EHT 2019: M = 6.5×10⁹ M_sun, ρ_∞ ≈ 10⁻²⁵ kg/m³ (intracluster). Ṁ_Bondi ≈ 10⁻³ M_sun/year. Relativistic jet 5,000 ly = 1.5 kpc, M87* power 6×10³⁵ W detected radio LBA Italian INAF
• Quasar 3C 273 Virgo 749 Mpc INAF TNG La Palma 1998: M_BH = 8.9×10⁸ M_sun, ρ ≈ 10⁻²³ kg/m³ Shakura-Sunyaev 1973 disk. Ṁ_Bondi ≈ 0.1 M_sun/year
• T-Tauri protostar ρ Ophiuchi giant molecular cloud 130 pc INAF Arcetri Florence 1872: M = 0.5 M_sun, T = 10 K (c_s = 0.2 km/s), ρ_∞ ≈ 10⁻¹⁹ kg/m³. r_B = 2.5×10¹⁵ m = 17,000 AU. Ṁ_Bondi ≈ 10⁻⁶ M_sun/year (1 million year formation)
• Hercules X-1 NS X-ray INAF Trieste 1898: NS 1.4 M_sun Roche-lobe overflow from 2 M_sun giant. Ṁ ≈ 10⁻⁹ M_sun/year, L_X ≈ 10³⁰ W Eddington bound 1.76×10³¹ W
• IMBH Holmberg IX X-1 super-Eddington INAF + NICER 2017 ASI-NASA: M_BH ≈ 10 M_sun, ρ_local ≈ 10⁻¹⁰ kg/m³ Roche-lobe overflow 5 M_sun giant. Ṁ_BHL ≈ 10⁻⁵ M_sun/year, L_X ≈ 10⁴⁰ erg/s super-Eddington outflow

Italian accretion astrophysics

Italian accretion astrophysics:

• ICRA International Center Relativistic Astrophysics Rome Sapienza 1985 + ICRANet Pescara 2003 Remo Ruffini director, GRB BdHN Binary-driven Hypernova model
• INAF Italian National Institute of Astrophysics Rome 2002 Monte Porzio Catone (16 observatories, 1,700 researchers). Observatories: Padua 1761 Galileo + Asiago Cima Ekar 1979 Copernico 182 cm; Brera Milan 1762 Schiaparelli; Capodimonte Naples 1812; Arcetri Florence 1872; Rome 1786; Trieste 1898; Cagliari 1899 SRT 64 m radio Sardinia 2013; Catania-Etna 1880; Bologna Loiano 1936; Palermo 1790. Telescopes: TNG National Galileo Telescope La Palma 1998 3.58 m INAF; LBT Mount Graham Arizona 2005 2×8.4 m INAF 25%; VST Cerro Paranal Chile 2011 2.6 m INAF + ESO
• Italian X-ray missions ASI 1988: BeppoSAX 1996-2002 GRB970228; XMM-Newton ESA 1999 RGS; NICER NASA 2017 ASI; IXPE 2021 NASA-ASI INAF Brera + INFN Pisa; Athena ESA 2035 X-IFU INAF Palermo
• Italian Physics Nobel Prizes: Marconi 1909, Fermi 1938, Segrè 1959, Carlo Rubbia 1984, Riccardo Giacconi 2002
• Andrea Ghez UCLA + Reinhard Genzel MPE 2020 Nobel Prize Sgr A* S2 star tracking 16 years Keck VLT
• SISSA Trieste 1978 astrophysics MHD; CINECA Bologna 1969 Marconi100 30 PFlops accretion disk simulations. Italian R&D tax bonuses 2025: R&D Tax Credit L. 160/2019 up to 25%; Patent Box L. 178/2020 50%; PNRR Mission 4 Component 2 €6.9B

Updated July 2026