Hernquist profile galaxy bulge mass calculator

The Hernquist profile calculator determines the spheroidal mass distribution of elliptical galaxy bulges and central bulges of spirals using the analytical model introduced by Lars Hernquist 1990 in "An Analytical Model for Spherical Galaxies and Bulges" published in The Astrophysical Journal vol. 356, p. 359 (Princeton/Harvard CfA).

Spheroidal density: ρ(r) = (M · a) / (2π) · 1 / [r · (r + a)³], where M is total bulge mass in M_sun, a is the Hernquist scale radius in kpc, r the spherical radius in kpc.

Cumulative mass within r: M(r) = M · r² / (r + a)² (closed-form analytical, essential for fast N-body simulations).

De Vaucouleurs effective radius r_e (projected semi-light radius): r_e ≈ 1.8153 · a, related to the R^(1/4) surface brightness profile of de Vaucouleurs 1948 Annales d'Astrophysique vol. 11, 247.

Physics: the Hernquist model is the 3D analytical counterpart of de Vaucouleurs profile obtained integrating density along line-of-sight; finite total mass (unlike Jaffe 1983 model log-divergent), central cusp ρ ∝ r⁻¹ (compatible with N-body cold dark matter Navarro-Frenk-White 1996), and asymptotic profile ρ ∝ r⁻⁴ for r >> a.

Hernquist gravitational potential: Φ(r) = -G·M / (r + a), very simple (essential for accelerations N-body Tree-code Barnes-Hut 1986).

4 Italian bulge categories:

• Classical bulges (formed by major mergers, M ~ 10¹⁰-10¹² M_sun, M31 Andromeda)
• Pseudo-bulges (formed by secular disk evolution, disk-like, M ~ 10⁹ M_sun, NGC 1097 SB)
• Composite bulges (Milky Way boxy/peanut + classical ~30%)
• Nuclear star clusters (NSC ~10⁶-10⁸ M_sun, separate from bulge, present in 80% late-type INAF SDSS 2014). Italian applications:

• Milky Way bulge M ~ 2×10¹⁰ M_sun a ~ 0.7 kpc studied VVV ESO + APOGEE-2 SDSS-IV INAF Bologna 1936 + INAF Trieste 1898 (boxy/peanut bar bulge)
• M31 Andromeda classical bulge M ~ 4×10¹⁰ M_sun a ~ 1 kpc Hubble HST INAF Padua 1761 + Sagittarius A* central nuclear cluster Sgr A* SMBH 4.3×10⁶ M_sun
• M87 giant elliptical cD bulge M ~ 3×10¹² M_sun a ~ 5 kpc EHT Event Horizon Telescope 2019 INAF + ESO 6.5×10⁹ M_sun BH Virgo cluster (M87* first BH image)
• NGC 4486B nuclear compact elliptical M ~ 10⁹ M_sun a ~ 0.1 kpc Virgo cluster INAF Capodimonte 1812 VST KIDS
• NGC 5128 Centaurus A merger elliptical M ~ 10¹¹ M_sun a ~ 2 kpc radio AGN INAF Bologna
• IC 1101 Abell 2029 cD bulge M ~ 10¹³ M_sun a ~ 50 kpc largest known galaxy. Dynamical model simulations use Hernquist as building block: Gala (Adrian Price-Whelan, ex-Princeton), Galpy (Jo Bovy Toronto), AGAMA (Eugene Vasiliev IoA Cambridge); N-body simulations Eagle 2015 + IllustrisTNG-300 2018 Heidelberg + Magneticum Pillepich use halo Hernquist initial conditions GalIC + DICE. Italian computing: CINECA Casalecchio Bologna 1969 Marconi100 30 PFlops + Leonardo HPC pre-exascale 250 PFlops 2023 (5th TOP500) cosmological N-body bulge high-resolution simulations ART-I-O INAF Trieste 1898 + SISSA Trieste 1978 + Scuola Normale Pisa 1810 (Andrea Ferrara cosmic dawn). Italian X-ray missions bulges: BeppoSAX 1996-2002 ASI Rome 1988 AGN bulges Seyfert NGC 1068 NGC 4151; Athena ESA 2035 X-IFU calorimeter INAF Palermo 1790 25% bulges X-ray hot gas + AGN feedback star formation quenching massive bulges. Italian Physics Nobel Prizes: Marconi 1909 + Fermi 1938 + Segrè 1959 + Carlo Rubbia 1984 (W/Z bosons CERN) + Riccardo Giacconi 2002 (X-ray astronomy, founder Chandra) + Ghez/Genzel 2020 Sgr A*. ICRA La Sapienza Rome 1985 Remo Ruffini central supermassive black holes bulges (M-σ relation Ferrarese-Merritt 2000)

Italian Hernquist bulges table

(Source: Hernquist 1990 ApJ 356, 359 + de Vaucouleurs 1948 r^(1/4) + Kormendy-Ho 2013 ARAA M-σ + EHT 2019 M87* + INAF Rome 2002 data + VVV ESO + APOGEE-2 SDSS-IV INAF Bologna).

Italian Hernquist profile examples

• Milky Way bulge VVV ESO + APOGEE-2 SDSS-IV INAF Bologna 1936 + INAF Trieste 1898: M = 2×10¹⁰ M_sun, a = 0.7 kpc. At r = 2 kpc → M(2) = 2×10¹⁰ · 4 / (2.7)² = 1.1×10¹⁰ M_sun within 2 kpc (54% total mass). Density at r = 2 kpc: ρ = 5.66×10⁷ M_sun/kpc³ = 0.0566 M_sun/pc³. r_e = 1.8153 · 0.7 = 1.27 kpc (matches de Vaucouleurs r^(1/4) brightness profile fit COBE-DIRBE NIR + 2MASS Bologna). Sgr A* SMBH 4.3×10⁶ M_sun (Ghez/Genzel 2020 Nobel) sentinel central nuclear cluster M_NSC ~ 2×10⁷ M_sun separate from bulge
• M31 Andromeda classical bulge Hubble HST INAF Padua 1761: M = 4×10¹⁰ M_sun, a = 1 kpc. At r = 5 kpc → M(5) = 4×10¹⁰ · 25 / 36 = 2.78×10¹⁰ M_sun within 5 kpc (69% total mass). ρ = 5.9×10⁶ M_sun/kpc³ = 5.9×10⁻³ M_sun/pc³. r_e = 1.82 kpc. P3 BH 1.4×10⁸ M_sun (Bender 2005 HST)
• M87 Virgo cD giant elliptical EHT 2019 INAF + ESO Capodimonte VST KIDS 2011-2019: M = 3×10¹² M_sun, a = 5 kpc. At r = 10 kpc → M(10) = 3×10¹² · 100 / 225 = 1.33×10¹² M_sun (44% mass). r_e = 1.8153 · 5 = 9.08 kpc. M87* SMBH 6.5×10⁹ M_sun (EHT first BH image 2019)
• NGC 4486B compact elliptical Virgo VST KIDS INAF Capodimonte 1812: M = 10⁹ M_sun, a = 0.1 kpc nuclear. At r = 0.3 kpc → M(0.3) = 10⁹ · 0.09 / 0.16 = 5.6×10⁸ M_sun (56%). r_e = 0.18 kpc compact. Double nucleus BH 5×10⁸ M_sun (Lauer 1996 HST)
• NGC 5128 Centaurus A radio AGN INAF Bologna 1936: M = 10¹¹ M_sun, a = 2 kpc post-merger. At r = 5 kpc → M(5) = 10¹¹ · 25 / 49 = 5.1×10¹⁰ M_sun. r_e = 3.63 kpc. Centaurus A radio jets AGN 5.5×10⁷ M_sun BH
• IC 1101 Abell 2029 cD largest known galaxy VST KIDS INAF + Athena ESA 2035 INAF Palermo 1790: M = 10¹³ M_sun, a = 50 kpc giant. At r = 100 kpc → M(100) = 10¹³ · 10000 / 22500 = 4.44×10¹² M_sun. r_e = 90.8 kpc. SMBH 4×10¹⁰ M_sun ultramassive
• Comparison with NGC 1399 Fornax cD INAF Bologna 1936 sentinel: M = 5×10¹¹ M_sun, a = 3.5 kpc. At r = 10 kpc → M(10) = 5×10¹¹ · 100 / 182.25 = 2.74×10¹¹ M_sun (55%). r_e = 6.35 kpc. ICRA Sapienza 1985 Remo Ruffini M-σ relation Ferrarese-Merritt 2000 calibration BH 5×10⁸ M_sun

Italian bulges astrophysics

Italian bulges astrophysics reconstructs galaxy dynamical models via the INAF Italian National Institute of Astrophysics Rome 2002 Monte Porzio Catone network (16 observatories, 1,700 researchers): Padua 1761 Galileo + Asiago 1979 Copernico 182 cm bulge spectroscopy; Brera Milan 1762 Schiaparelli; Capodimonte Naples 1812 VST KIDS bulge photometry weak lensing 1,500 deg²; Arcetri Florence 1872 Galileo; Rome 1786; Trieste 1898 Planck CMB hardware + INFN Pisa LFI + bulges spectroscopy SDSS-IV APOGEE-2; Catania-Etna 1880; Bologna Loiano 1936 AGN bulges Seyfert Centaurus A radio + Cassini-Loiano 152 cm; Palermo 1790; Cagliari 1899 SRT 64 m radio.

Italian bulge telescopes: TNG La Palma 1998 3.58 m INAF; LBT Mount Graham Arizona 2005 2×8.4 m INAF 25% (LUCI MOS bulges spectra); VST Cerro Paranal Chile 2011 2.6 m INAF Capodimonte + ESO (KIDS bulge photometry); SRT Sardinia 2013 64 m radio.

Italian bulges missions ASI Rome 1988: Planck Surveyor ESA 2009-2018 CMB σ_8=0.81 H_0=67.4 (INAF Trieste hardware) + bulges integrated luminosity function; BeppoSAX 1996-2002 ASI X-ray AGN bulges Seyfert NGC 1068 NGC 4151; XMM-Newton ESA 1999 Italian RGS Padua-Bologna ICM hot gas elliptical bulges; Chandra NASA 1999 X-ray imaging bulges resolved (Carlo Rubbia 1984 Nobel CERN W/Z + Riccardo Giacconi 2002 Nobel X-ray Chandra founder); JWST 2021 NIRSpec INAF proto-bulges z>5 cosmic dawn (Roberto Maiolino Cambridge ex-INAF Arcetri); Euclid ESA 2023 ASI 25% 15,000 deg² 25 mln galaxy bulge photometry 6 years; Athena X-ray Observatory ESA 2035 X-IFU calorimeter INAF Palermo 1790 25% bulges hot gas + AGN feedback quenching massive ellipticals.

Italian Physics Nobel Prizes: Marconi 1909, Fermi 1938, Segrè 1959, Carlo Rubbia 1984 (W/Z CERN), Riccardo Giacconi 2002 (X-ray Chandra). Andrea Ghez UCLA + Reinhard Genzel MPE 2020 Sgr A*.

ICRA Rome Sapienza 1985 + ICRANet Pescara 2003 Remo Ruffini supermassive bulge BHs (M-σ Ferrarese-Merritt 2000); SISSA Trieste 1978 bulge secular evolution astrophysics; Scuola Normale Pisa 1810 Andrea Ferrara cosmic dawn proto-bulges z>5; CINECA Casalecchio Bologna 1969 Marconi100 30 PFlops + Leonardo HPC pre-exascale 250 PFlops 2023 (5th TOP500) N-body cosmological bulge high-resolution Eagle 2015 + IllustrisTNG-300 2018 Heidelberg + Magneticum Pillepich + ART-I-O INAF + SISSA.

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 + ESFRI ELT-LBT-SKAO ESA Athena 2035 X-IFU INAF Palermo 25%.

Updated July 2026