Astrophysics and Stellar Physics

Introduction to Astrophysics

Astrophysics applies physics laws—gravity, thermodynamics, quantum mechanics, electromagnetism, nuclear physics, relativity—to understand the formation, structure, evolution, and ultimate fate of celestial objects and the universe as a whole. From the Sun's nuclear furnace to black hole event horizons, from the cosmic microwave background to the large-scale web of galaxies, astrophysics has revealed a universe far more vast, violent, and strange than pre-telescopic imagination could conceive. The observable universe spans ~93 billion light-years across, contains ~2 trillion galaxies, and is 13.8 billion years old—all established through physical reasoning from observations of electromagnetic radiation, gravitational waves, cosmic rays, and neutrinos reaching Earth.

Key milestones: Cecilia Payne-Gaposchkin (1925)—Sun's composition mostly hydrogen, not rock-like; Hans Bethe (1939)—solar energy from pp chain nuclear fusion; discovery of pulsars (1967, Jocelyn Bell)—rotating neutron stars; first quasar identification (1963, Maarten Schmidt 3C 273); discovery of CMB (1965, Penzias and Wilson)—relic Big Bang radiation; Hubble Space Telescope deep field images (1996) revealing thousands of galaxies in a blank patch of sky; first gravitational wave detection (2015, LIGO); first black hole image (2019, Event Horizon Telescope). Modern multi-messenger astronomy combining gravitational waves, neutrinos, and photons across all wavelengths provides unprecedented physical diagnostics of cosmic events.

Stellar Structure and Evolution

Stellar Equilibrium and the Main Sequence

Stars spend most of their lives in hydrostatic equilibrium—pressure gradient force balancing gravitational self-attraction: dP/dr = -GM(r)rho(r)/r^2. For main-sequence stars burning hydrogen in cores, energy transport by radiation (in radiative zones) and convection (when temperature gradient is super-adiabatic) carries luminosity outward. Mass-luminosity relationship L ∝ M^4 (approximately) on the main sequence—massive stars are vastly more luminous but burn their fuel far faster (lifetime t ∝ M/L ∝ M^-3: a 10 M_sun star lives ~10 Myr, a 1 M_sun star ~10 Gyr). The Hertzsprung-Russell diagram—luminosity versus surface temperature (spectral type)—organises stars into main sequence (hydrogen burning), red giants (shell hydrogen burning), horizontal branch (core helium burning), and white dwarf sequence, providing a visual summary of stellar evolution. Stellar populations: Population I (metal-rich, young disk stars like the Sun); Population II (metal-poor, old halo and globular cluster stars formed from pristine Big Bang composition); Population III (hypothetical first metal-free massive stars formed at z~20—sought through JWST as gravitational lensing observations).

Stellar Death and Compact Objects

Stellar fate depends on initial mass. Low mass stars (M < ~8 M_sun): after giant branch evolution, stellar envelope expelled as planetary nebula; core collapses to white dwarf (~0.6 M_sun of C-O, or O-Ne for intermediate mass); degenerate electron pressure provides support against gravity. Massive stars (M > ~8 M_sun): core silicon burning builds iron core to Chandrasekhar mass 1.4 M_sun; electron capture begins, degeneracy pressure fails, core collapses in ~0.1 seconds to ~10 km neutron star (supported by neutron degeneracy and nuclear repulsion)—gravitational energy release (~3×10^46 J, more than 1% of solar rest mass energy) powers core-collapse supernova. If remnant mass >~3 M_sun (Tolman-Oppenheimer-Volkoff limit), nothing can prevent collapse to black hole. Pair instability supernovae in very massive stars (140-260 M_sun): electron-positron pair production softens equation of state triggering runaway nuclear burning, completely disrupting the star with no remnant.

Black Holes

Black Hole Physics

Black holes are regions of spacetime where gravity is too strong for anything—including light—to escape from within the event horizon (Schwarzschild radius r_s = 2GM/c^2 = 3 km × M/M_sun). Stellar mass black holes (~3-100 M_sun) form from massive star collapse; supermassive black holes (10^6-10^10 M_sun) reside at galaxy centres—including Sgr A* (4×10^6 M_sun, imaged by EHT 2022) at the Milky Way centre. Rotating Kerr black holes are characterised by mass M and spin a—frame dragging (ergosphere region outside event horizon where spacetime rotation forces co-rotation) enables the Penrose process extracting rotational energy. Hawking radiation: quantum field theory in curved spacetime predicts black holes radiate thermally at T_H = hbar*c^3/(8pi*G*M*k_B)—stellar mass black holes have T_H = 6×10^-8 K (undetectable), primordial mini black holes have higher T_H and would be evaporating today. The information paradox—whether Hawking evaporation loses information (violating quantum mechanics)—remains a major unsolved problem at the quantum gravity frontier.

Galaxies and Cosmology

Galaxies—gravitationally bound systems of 10^7–10^12 stars, gas, dust, and dark matter—span ellipticals, spirals (Milky Way, Andromeda), and irregulars. Galaxy rotation curves (flat rather than Keplerian fall-off at large radii) indicate ~85% of mass in non-luminous dark matter halos—the most compelling observational evidence for dark matter. The Milky Way contains ~250-400 billion stars, 100,000 light-year diameter disk, ~50 kpc dark matter halo. Large-scale structure of the universe: galaxies cluster in groups, clusters (Virgo cluster, Coma cluster), superclusters (Laniakea with 10^17 M_sun) and cosmic web filaments/sheets/voids—quantified by the two-point correlation function and power spectrum measured by SDSS, 2dFGRS, and DESI galaxy surveys. The DESI Baryon Acoustic Oscillation survey (2024 results from 6 million galaxies) measured dark energy with unprecedented precision, finding potential hints of dynamic dark energy evolution.

Examples and Applications

Example 1: Gravitational Wave Astronomy

LIGO, Virgo, and KAGRA gravitational wave interferometers detect spacetime distortions h = Delta L / L ~ 10^-21 from binary compact object mergers—equivalent to measuring 1/1000 of a proton diameter change in 4km arms. GW150914 (September 2015)—first detection, binary black hole merger ~30+36 M_sun at 1.3 billion light-years producing peak luminosity ~3.6×10^49 W (surpassing all stars in observable universe combined briefly). GW170817 (2017)—neutron star merger detected simultaneously in gravitational waves and gamma-ray burst (GRB), followed by kilonova optical/IR counterpart as r-process elements (gold, platinum, europium) synthesised. The O3/O4 observing runs have detected 90+ compact binary events in the gravitational wave catalogue (GWTC-3). Future: LISA space-based interferometer (3 spacecraft, 2.5 million km arms) targeting supermassive black hole mergers, EMRIs, and stochastic gravitational wave background; Einstein Telescope (underground, 10 km L-shaped); Cosmic Explorer (40 km arms) for Earth-based next generation.

Example 2: The James Webb Space Telescope

JWST (launched December 2021, at L2 Lagrange point 1.5 million km from Earth) is a 6.5m primary mirror infrared space telescope—the successor to HST—operating 0.6-28 micrometers providing diffraction-limited resolution in near-to-mid-IR. Its capabilities: detecting galaxies at z>12 (less than 400 million years after Big Bang) observed as UV photons redshifted into IR; characterising exoplanet atmospheres via transmission spectroscopy of transits—detecting CO2, H2O, CH4, SO2 in hot Jupiter and super-Earth atmospheres with unprecedented detail (TRAPPIST-1 system atmospheric observations); imaging protostellar jets and disk formation (Herbig-Haro objects) in unprecedented clarity; and imaging Solar System objects (Europa ocean plumes confirmed, Neptune ring details). JWST's early universe observations found unexpectedly massive, bright, early galaxies (Maisie's Galaxy at z=11.4, CEERS-1749 at z=17 claimed) challenging galaxy formation models by having higher stellar masses earlier than predicted.

Example 3: Neutron Star Equation of State

Neutron stars (M ~1.2-2.0 M_sun, R ~10-13 km, central density ~5-10 times nuclear saturation density) are the densest objects with well-measured properties—probing the equation of state (pressure vs. density) of cold dense matter inaccessible to terrestrial experiments. Multi-messenger constraints: GW170817 tidal deformability measurement (Lambda < 800) constrains the equation of state at 1-2x nuclear density; NICER X-ray timing mission measures pulsars' radii by modelling thermal X-ray surface emission pulsation shapes—achieved ~1 km precision radius measurements for PSR J0030+0451 (R=13.0 km) and PSR J0740+6620 (heaviest neutron star measured, 2.08 M_sun, R=12.4 km). Maximum mass (twice solar) constraints from Shapiro delay measurements of pulsars in binary systems disfavour purely nucleonic models at high density—possibly requiring quark matter cores or hyperons with strong repulsive interactions. Understanding dense matter at neutron star conditions has implications for heavy-ion physics (RHIC/LHC quark-gluon plasma) and the r-process.

Example 4: Exoplanet Atmospheric Physics

Comparative planetology of exoplanet atmospheres is now a mature JWST-era field. Atmospheric characterisation techniques: primary eclipse transmission spectroscopy (stellar light filtered through planet atmosphere during transit reveals absorption features at specific wavelengths); secondary eclipse thermal emission spectroscopy (planet-only thermal emission subtracted when planet passes behind star); phase curve monitoring (full orbital thermal emission variation mapping east-west temperature contrast). Hot Jupiters have dayside temperatures of 1000-3000 K, efficient day-night heat redistribution (eastward equatorial superrotation observed via phase curve offsets), and primarily clear H/He atmospheres with metal hydrides, alkalis, and Ti/VO at high temperatures. TRAPPIST-1 b and c JWST thermal emission measurements constrained surface temperatures inconsistent with thick CO2 atmospheres. Detectability of biosignature gases (O2, O3, CH4, N2O) in temperate rocky planet atmospheres from the TRAPPIST-1 system requires ~100 JWST transit observations—achievable within JWST's 20-year mission lifetime if prioritised.

Example 5: Type Ia Supernovae and Dark Energy

Type Ia supernovae—thermonuclear explosions of white dwarfs at or near the Chandrasekhar mass—are cosmological standard candles (intrinsic luminosity ~5×10^9 L_sun) with standardisable peak brightness via the Phillips relation (brighter supernovae decline more slowly after peak). The Supernova Cosmology Project (Perlmutter) and High-Z SN Search Team (Riess, Schmidt) measured SN Ia distances at z~0.3-1.0 in the 1990s, finding distant supernovae dimmer than expected for decelerating expansion—discovering cosmic acceleration and dark energy (Nobel 2011). Current constraints: dark energy density parameter OmegaLambda ~0.69, dark matter OmegaDM ~0.26, baryon matter Omega_b ~0.05 (Lambda-CDM concordance cosmology). The Hubble tension—H0 measurements from early universe CMB (~67.4 km/s/Mpc, Planck 2018) persistently disagree with late universe distance ladder (73.4 km/s/Mpc, SH0ES Riess team)—at ~5 sigma tension suggesting either systematic errors or new physics (early dark energy, extra relativistic species, modified gravity).

Example 6: Active Galactic Nuclei and Quasars

Active galactic nuclei (AGN) are powered by accretion of matter onto supermassive black holes (SMBH) in galaxy centres—releasing gravitational energy as accretion disk radiation, relativistic jets, and ionisation-driven winds. Quasars (quasi-stellar objects)—the most luminous AGN—can outshine their entire host galaxy by 100-1000× at peak, reaching L ~10^47 erg/s from accretion rates of ~10 M_sun/year. The unified model of AGN: a central SMBH with accretion disk surrounded by a dense dusty torus—orientation to the line of sight determines the observational class (Seyfert 1/2, quasar, blazer, radio galaxy). M87's SMBH (6.5×10^9 M_sun, 55 billion light-years distant) imaged by EHT 2019 as a bright accretion ring around a dark shadow matching general relativistic prediction for the photon orbit size. AGN feedback—enormous energy input into host galaxy interstellar medium—quenches star formation and establishes the observed M_SMBH-sigma* correlation between black hole mass and host galaxy bulge velocity dispersion, coupling SMBH and galaxy evolution.

Example 7: Cosmic Rays and High-Energy Astrophysics

Cosmic rays—charged particles (protons 90%, helium 9%, heavier nuclei, electrons, positrons) spanning energies from ~GeV to ~3×10^20 eV (ultra-high energy cosmic rays, UHECRs)—impinge on Earth from all directions. Below ~10^15 eV (the knee), cosmic rays are Galactic in origin (accelerated by supernova remnant shocks); above the ankle (~10^18.5 eV), they are extragalactic. UHECRs above ~6×10^19 eV (GZK limit: proton interactions with CMB photons producing pions, limiting range to ~100 Mpc) should show anisotropy toward sources—Auger Observatory observations correlate with Centaurus A and starburst galaxies within 100 Mpc. Gamma-ray astronomy (Fermi-LAT, CTA): gamma-ray bursts (GRBs—most energetic transients in the universe, ~10^44-10^47 J in seconds), gamma-ray emission from SNRs (supernova remnant shock acceleration), and pairs of photons from Markarian 421/501 blazar jets probe Lorentz invariance violation at quantum gravity energy scales. IceCube neutrino telescope (1 km^3 ice detector at South Pole) detected astrophysical neutrino flux and identified NGC 1068 Seyfert galaxy as a high-energy neutrino source.

Example 8: The Milky Way Structure

Mapping the Milky Way's structure from inside requires indirect methods—parallaxes, RR Lyrae standard candles, pulsars, masers. Gaia space astrometry mission (2013-present) has measured parallaxes and proper motions of 1.5 billion stars with sub-microarcsecond precision—revealing the Milky Way's spiral arm structure, the warp and flare of the outer disk, numerous stellar streams (remnants of disrupted globular clusters and dwarf galaxies), and the Gaia-Sausage-Enceladus merger event ~10 Gyr ago when a Milky-Way mass-ratio dwarf galaxy was accreted, depositing retrograde halo stars and sequentially truncating the thick disk formation. The Galactic Centre region—within 200 pc of Sgr A*—shows dense stellar populations (Nuclear Star Cluster, ~3×10^7 M_sun), young massive star clusters (Arches, Quintuplet) formed ~3 Myr ago from compressed gas clouds, and stellar orbits (S2, with ~16-year period just 100 AU from Sgr A*) providing the most precise measurement of Sgr A* mass and General Relativistic precession in strong gravity fields.

Try it live

Everything above runs in your browser — open SPH Fluid and change the parameters while it is running. Nothing is installed, nothing is uploaded, the whole model lives in one tab.