m physical modeling / IICSM

CONNECTED PHYSICS / STATES OF MATTER

Astrophysics: from radiation to stars and galaxies

Connect observed light to physical properties, derive stellar and orbital models, and solve 20 worked problems spanning stars, planets, gas clouds, and galaxies.

Subject library · 51 guides · derivations & worked examples

Matter pathway: atom → solid → liquid → gas → plasma. Quantum mechanics and quantum field theory provide foundations across the pathway; they are not additional phases. This is a connected modeling route, not a universal heating curve. Actual phases depend on pressure, composition, and kinetics.

1. Infer distance and luminosity from observed light

Definitions & inputs. L is total luminosity (W), F received bolometric flux (W m⁻²), d distance, and p annual parallax angle. m and M are apparent and absolute magnitudes in the same photometric band.

  1. Conservation spreads emitted power over a sphere; flux is power per collecting area, not total power.

    L=∫F dA=4πd2F⇒F=L4πd2L=\int F\,dA=4\pi d^2F\quad\Rightarrow\quad F=\frac{L}{4\pi d^2}
  2. For a small angle, the apparent displacement relative to the solar-system barycentre gives distance. Annual parallax is the semi-angle, not the full six-month displacement.

    p≃1 AUd (radians),d(pc)=1p(arcsec)p\simeq\frac{1\ {\rm AU}}d\ ({\rm radians}),\quad d({\rm pc})=\frac1{p({\rm arcsec})}
  3. The magnitude scale is logarithmic: brighter objects have smaller magnitudes.

    m1−m2=−2.5log⁡10(F1/F2)m_1-m_2=-2.5\log_{10}(F_1/F_2)
  4. Insert inverse-square flux into the magnitude difference. A is extinction in magnitudes in the chosen band.

    M=m(d=10 pc),m−M=5log⁡10[d/(10 pc)]+AM=m(d=10\ {\rm pc}),\quad m-M=5\log_{10}[d/(10\ {\rm pc})]+A

Interpretation. Distance, extinction, and bolometric correction are separate inference steps. A band-limited flux is not automatically the bolometric flux used with total luminosity.

↑ Return to definitions and contents

2. Thermal spectra, stellar radius, and Doppler shifts

Definitions & inputs. Bλ is spectral radiance per wavelength per solid angle; T temperature; σSB the Stefan–Boltzmann constant; R radius; λ0 rest wavelength; vr radial velocity positive for recession.

  1. Planck’s distribution combines photon mode density and Bose occupation. This is radiance per unit wavelength, not per unit frequency.

    Bλ(T)=2hc2λ51ehc/(λkBT)−1B_\lambda(T)=\frac{2hc^2}{\lambda^5}\frac1{e^{hc/(\lambda k_BT)}-1}
  2. Integrate over the outward hemisphere and wavelength, then multiply by the stellar surface area. Effective temperature is defined by this luminosity relation.

    Fsurface=∫0∞πBλdλ=σSBT4,L=4πR2σSBTeff4F_{\rm surface}=\int_0^\infty\pi B_\lambda d\lambda=\sigma_{\rm SB}T^4,\quad L=4\pi R^2\sigma_{\rm SB}T_{\rm eff}^4
  3. Set ∂Bλ/∂λ=0. The nonzero root x≈4.96511 yields Wien’s displacement law; the Bν peak is at a different location.

    x=hc/(λkBT),5(1−e−x)=x,λmax⁡T=2.897771955×10−3 m Kx=hc/(\lambda k_BT),\quad5(1-e^{-x})=x,\quad\lambda_{\max}T=2.897771955\times10^{-3}\ {\rm m\,K}
  4. A small nonrelativistic Doppler shift gives line-of-sight velocity. Gravitational and cosmological redshifts require their own interpretation.

    z=λobs−λ0λ0≃vrc(∣vr∣≪c)z=\frac{\lambda_{\rm obs}-\lambda_0}{\lambda_0}\simeq\frac{v_r}{c}\quad(|v_r|\ll c)

Interpretation. Combining distance, luminosity, color, and spectroscopy constrains stars, but each inference retains model assumptions.

↑ Return to definitions and contents

3. Infer masses from orbital dynamics

Definitions & inputs. M1,M2 are two masses, a is the semimajor axis of their relative separation, P orbital period, and i inclination. G is Newton’s constant.

  1. Subtract the two body equations of motion to obtain the relative orbit.

    r¨=−G(M1+M2)r3r\ddot{\mathbf r}=-\frac{G(M_1+M_2)}{r^3}\mathbf r
  2. For a circular orbit, equate centripetal and gravitational acceleration. The same Kepler law holds for an elliptical orbit with semimajor axis a.

    circular: 4π2aP2=G(M1+M2)a2⇒P2=4π2a3G(M1+M2)\text{circular: }\frac{4\pi^2a}{P^2}=\frac{G(M_1+M_2)}{a^2}\Rightarrow P^2=\frac{4\pi^2a^3}{G(M_1+M_2)}
  3. The centre-of-mass condition relates the two orbital sizes.

    M1a1=M2a2,a=a1+a2M_1a_1=M_2a_2,\quad a=a_1+a_2
  4. The primary’s radial-velocity semiamplitude includes eccentricity e and inclination. It does not alone determine the companion mass without other information.

    K1=2πa1sin⁡iP1−e2K_1=\frac{2\pi a_1\sin i}{P\sqrt{1-e^2}}

Interpretation. Use the relative separation for a, not one component’s distance from the centre of mass.

↑ Return to definitions and contents

4. Derive the equations of stellar structure

Definitions & inputs. Mr is mass enclosed within radius r, ρ density, P total pressure, Lr outward luminosity, ε net local energy generation per unit mass (W kg⁻¹), and κR Rosseland mean opacity (m² kg⁻¹).

  1. A thin spherical shell has volume 4πr²dr.

    dMr=4πr2ρ dr⇒dMrdr=4πr2ρdM_r=4\pi r^2\rho\,dr\Rightarrow\frac{dM_r}{dr}=4\pi r^2\rho
  2. Balance the shell’s pressure force against gravity; pressure decreases outward.

    4πr2dP=−GMrr2dMr⇒dPdr=−GMrρr24\pi r^2dP=-\frac{GM_r}{r^2}dM_r\Rightarrow\frac{dP}{dr}=-\frac{GM_r\rho}{r^2}
  3. Local energy production changes the luminosity passing through each shell. Time-dependent thermal evolution adds entropy terms.

    dLr=ϵ dMr⇒dLrdr=4πr2ρϵdL_r=\epsilon\,dM_r\Rightarrow\frac{dL_r}{dr}=4\pi r^2\rho\epsilon
  4. In optically thick radiative diffusion, insert F=Lr/(4πr²). A convective zone needs a different heat-transport prescription.

    Frad=−16σSBT33κRρdTdr,dTdr=−3κRρLr64πσSBr2T3F_{\rm rad}=-\frac{16\sigma_{\rm SB}T^3}{3\kappa_R\rho}\frac{dT}{dr},\quad\frac{dT}{dr}=-\frac{3\kappa_R\rho L_r}{64\pi\sigma_{\rm SB}r^2T^3}
  5. Illustrative ideal-gas and radiation terms depend on mean particle mass μmp. Degeneracy or strong interactions require other equations of state.

    Pgas=ρkBTμmp,Prad=4σSB3cT4P_{\rm gas}=\frac{\rho k_BT}{\mu m_p},\quad P_{\rm rad}=\frac{4\sigma_{\rm SB}}{3c}T^4

Interpretation. Hydrostatic balance is not thermal balance. These coupled equations must normally be integrated numerically rather than replaced by one uniform-density estimate.

↑ Return to definitions and contents

5. Virial balance and stellar energy timescales

Definitions & inputs. Ω is negative gravitational potential energy, U internal energy, N particle count, μ mean molecular weight in proton-mass units, and L luminosity.

  1. Multiply hydrostatic balance by 4πr³dr and recognize the gravitational binding-energy integral.

    ∫0R4πr3dPdrdr=−∫0MGMrrdMr=Ω\int_0^R4\pi r^3\frac{dP}{dr}dr=-\int_0^M\frac{GM_r}{r}dM_r=\Omega
  2. Integrate the left side by parts; discard surface pressure only under the stated assumption.

    4πR3P(R)−3∫P dV=Ω,P(R)≃0⇒3∫P dV+Ω=04\pi R^3P(R)-3\int P\,dV=\Omega,\quad P(R)\simeq0\Rightarrow3\int P\,dV+\Omega=0
  3. For a monatomic ideal gas, gravitational contraction makes the total energy more negative while raising internal energy.

    U=32∫P dV⇒2U+Ω=0,E=U+Ω=Ω/2U=\tfrac32\int P\,dV\Rightarrow2U+\Omega=0,\quad E=U+\Omega=\Omega/2
  4. Insert the uniform-sphere binding energy and U=3NkB⟨T⟩/2 with N=M/(μmp). This is a particle-weighted mean temperature, not the central temperature.

    ρ=constant:Ω=−3GM25R,⟨T⟩=GMμmp5kBR\rho=\text{constant}:\quad\Omega=-\frac{3GM^2}{5R},\quad\langle T\rangle=\frac{GM\mu m_p}{5k_BR}
  5. The first is a gravitational energy-scale estimate; the second uses burned hydrogen mass fraction fH and conversion efficiency η. Neither assumes L remains constant in a realistic evolution calculation.

    tKH∼GM2RL,tfuel∼ηfHMc2Lt_{\rm KH}\sim\frac{GM^2}{RL},\quad t_{\rm fuel}\sim\frac{\eta f_HMc^2}{L}

Interpretation. Nuclear fuel can sustain luminosity much longer than gravitational contraction alone. Dense stellar remnants require quantum equations of state and sometimes GR.

↑ Return to definitions and contents

6. Radiative transfer, optical depth, and the Eddington scale

Definitions & inputs. Iν is specific intensity; αν absorption coefficient (m⁻¹); jν emission coefficient; κ is a flux-mean mass opacity; τ is optical depth; Sν=jν/αν is the source function.

  1. Along the propagation direction absorption removes intensity and emission adds it.

    dIνds=−ανIν+jν,dτν=ανds\frac{dI_\nu}{ds}=-\alpha_\nu I_\nu+j_\nu,\quad d\tau_\nu=\alpha_\nu ds
  2. For constant source function, solve the first-order equation with an integrating factor. Pure attenuation sets Sν=0.

    dIνdτν=−Iν+Sν⇒Iν(τ)=Iν(0)e−τ+Sν(1−e−τ)\frac{dI_\nu}{d\tau_\nu}=-I_\nu+S_\nu\Rightarrow I_\nu(\tau)=I_\nu(0)e^{-\tau}+S_\nu(1-e^{-\tau})
  3. Radiation carries momentum flux F/c; multiply by opacity per unit mass to get acceleration.

    arad=κF/c=κL4πr2c,ag=GM/r2a_{\rm rad}=\kappa F/c=\frac{\kappa L}{4\pi r^2c},\quad a_g=GM/r^2
  4. Balance the two accelerations. Radius cancels for spherical inverse-square flux.

    arad=ag⇒LEdd=4πGMcκa_{\rm rad}=a_g\Rightarrow L_{\rm Edd}=\frac{4\pi GMc}{\kappa}

Interpretation. The Eddington luminosity is a model-dependent force-balance scale, not an absolute brightness limit for every source geometry.

↑ Return to definitions and contents

7. Gas-cloud collapse: free-fall time and Jeans instability

Definitions & inputs. ρ0 is background mass density, cs is isothermal sound speed, k perturbation wave number, and ω its frequency. μmp is the mean particle mass.

  1. Integrate energy conservation for a shell released from rest; use M=4πρ0R³/3. This is a pressureless dynamical time, not a universal star-formation time.

    r˙2=2GM(1/r−1/R),tff=∫0Rdr2GM(1/r−1/R)=3π32Gρ0\dot r^2=2GM(1/r-1/R),\quad t_{\rm ff}=\int_0^R\frac{dr}{\sqrt{2GM(1/r-1/R)}}=\sqrt{\frac{3\pi}{32G\rho_0}}
  2. Linearize continuity, momentum balance, and Poisson’s equation about the chosen background.

    ∂tδρ+ρ0∇⋅δv=0,∂tδv=−cs2ρ0∇δρ−∇δΦ,∇2δΦ=4πGδρ\partial_t\delta\rho+\rho_0\nabla\cdot\delta\mathbf v=0,\quad\partial_t\delta\mathbf v=-\frac{c_s^2}{\rho_0}\nabla\delta\rho-\nabla\delta\Phi,\quad\nabla^2\delta\Phi=4\pi G\delta\rho
  3. Eliminate velocity and potential, then substitute a plane wave. Negative ω² indicates exponential growth.

    ∂t2δρ=cs2∇2δρ+4πGρ0δρ⇒ω2=cs2k2−4πGρ0\partial_t^2\delta\rho=c_s^2\nabla^2\delta\rho+4\pi G\rho_0\delta\rho\Rightarrow\omega^2=c_s^2k^2-4\pi G\rho_0
  4. The zero-frequency boundary defines the Jeans wavelength. Here Jeans mass means the mass within a sphere of radius λJ/2; other conventions have different factors.

    cs=kBTμmp,λJ=csπGρ0,MJ=4π3ρ0(λJ/2)3c_s=\sqrt{\frac{k_BT}{\mu m_p}},\quad\lambda_J=c_s\sqrt{\frac\pi{G\rho_0}},\quad M_J=\frac{4\pi}3\rho_0(\lambda_J/2)^3

Interpretation. Magnetic support, rotation, turbulent pressure, cooling, and geometry change collapse thresholds. Define the mass convention before comparing values.

↑ Return to definitions and contents

8. Planet transits and orbital inference

Definitions & inputs. Rp and R★ are planet and stellar radii; δ is the fractional flux loss; a orbital semimajor axis; P period; b transit impact parameter in stellar-radius units.

  1. Divide blocked projected area by the stellar disk area.

    Fblocked∝πRp2,F⋆∝πR⋆2⇒δ=(Rp/R⋆)2F_{\rm blocked}\propto\pi R_p^2,\quad F_\star\propto\pi R_\star^2\Rightarrow\delta=(R_p/R_\star)^2
  2. A transit measures a radius ratio; absolute planet radius inherits uncertainty in the star’s radius.

    Rp=R⋆δR_p=R_\star\sqrt\delta
  3. Repeated transits determine period; a stellar mass model then converts it to an orbital scale.

    Mp≪M⋆:a=(GM⋆P24π2)1/3M_p\ll M_\star:\quad a=\left(\frac{GM_\star P^2}{4\pi^2}\right)^{1/3}
  4. Viewing geometry controls whether a transit occurs and its chord length. A shallow dip alone does not uniquely establish a planet.

    b=acos⁡iR⋆(circular orbit)b=\frac{a\cos i}{R_\star}\quad(\text{circular orbit})

Interpretation. Spectroscopy, stellar characterization, multiple transits, and false-positive checks complement the simplified area model.

↑ Return to definitions and contents

9. Galaxies, dynamical mass, and the cosmic distance connection

Definitions & inputs. v is circular speed at radius r, M(<r) enclosed mass, H0 a chosen present expansion rate, z cosmological redshift, and dL luminosity distance.

  1. Balance centripetal and gravitational acceleration. Disks require their actual gravitational geometry for accurate mass inference.

    v2/r=GM(<r)/r2⇒M(<r)=v2r/Gv^2/r=GM(<r)/r^2\Rightarrow M(<r)=v^2r/G
  2. A flat rotation curve does not match a finite central point mass whose Keplerian speed decreases as r⁻¹/².

    v≃constant⇒M(<r)∝rv\simeq\text{constant}\Rightarrow M(<r)\propto r
  3. Low-redshift recession offers an approximate distance estimate. The luminosity-distance definition retains inverse-square notation while incorporating expansion effects.

    cz≃H0d(z≪1),F=L/(4πdL2)cz\simeq H_0d\quad(z\ll1),\qquad F=L/(4\pi d_L^2)
  4. At larger redshift, integrate the expansion history. This expression assumes spatial flatness and standard photon propagation.

    flat FLRW:dL(z)=(1+z)c∫0zdz′H(z′)\text{flat FLRW:}\quad d_L(z)=(1+z)c\int_0^z\frac{dz^\prime}{H(z^\prime)}

Interpretation. Astrophysical observations constrain models rather than directly measuring all intrinsic properties. Geometry, dust, selection effects, and calibration propagate into the result.

↑ Return to definitions and contents

Twenty worked astrophysics examples

Inputs are illustrative reference values, not current catalog measurements. Solar reference values used here are M⊙=1.98847×10³⁰ kg, R⊙=6.957×10⁸ m, and L⊙=3.828×10²⁶ W. One Julian year is 365.25 days. Each example links to its governing derivation.

Example 01. Distance from a 25-milliarcsecond parallax

Definitions & inputs. Annual parallax p=25 mas.

Review the governing derivation ↑

  1. Convert milliarcseconds to arcseconds.

    p=25 mas=0.025 arcsecp=25\ {\rm mas}=0.025\ {\rm arcsec}
  2. Invert the angle in the parsec definition.

    d=1/p=1/0.025=40.0 pcd=1/p=1/0.025=40.0\ {\rm pc}
  3. Convert to SI distance.

    d=40(3.08567758×1016)=1.23427e+18 md=40(3.08567758\times10^{16})=1.23427e+18\ {\rm m}

Interpretation. If parallax uncertainty is large, direct inversion can be a biased distance estimator.

↑ Return to definitions and contents

Example 02. Flux of a solar-luminosity star at 10 pc

Definitions & inputs. L=3.828×10²⁶ W; d=10 pc.

Review the governing derivation ↑

  1. Convert distance to metres.

    d=10 pc=3.08567758×1017 md=10\ {\rm pc}=3.08567758\times10^{17}\ {\rm m}
  2. Divide emitted power by the sphere’s area.

    F=L/(4πd2)F=L/(4\pi d^2)
  3. This is the total energy flux integrated over wavelength.

    F=3.19934e−10 W m−2F=3.19934e-10\ {\rm W\,m^{-2}}

Interpretation. At twice the distance the flux falls to one quarter.

↑ Return to definitions and contents
Inverse-square bolometric flux of a solar-luminosity source with no extinction. Both axes are logarithmic.
Inverse-square bolometric flux of a solar-luminosity source with no extinction. Both axes are logarithmic. Download SVG · Plot data (JSON)

Example 03. Distance from apparent and absolute magnitudes

Definitions & inputs. m=10.0, M=5.0, extinction A=0.0 in the same band.

Review the governing derivation ↑

  1. Insert the observed and intrinsic magnitudes.

    m−M−A=5=5log⁡10[d/(10 pc)]m-M-A=5=5\log_{10}[d/(10\ {\rm pc})]
  2. Divide the modulus by five.

    log⁡10[d/(10 pc)]=1\log_{10}[d/(10\ {\rm pc})]=1
  3. Exponentiate and restore the ten-parsec normalization.

    d=10×101=100 pcd=10\times10^1=100\ {\rm pc}

Interpretation. Positive uncorrected extinction would make this inferred distance too large.

↑ Return to definitions and contents

Example 04. Temperature from a wavelength-spectrum peak

Definitions & inputs. A blackbody Bλ spectrum peaks at λmax=500 nm.

Review the governing derivation ↑

  1. Convert the wavelength.

    λmax⁡=500 nm=5.00×10−7 m\lambda_{\max}=500\ {\rm nm}=5.00\times10^{-7}\ {\rm m}
  2. Use the nonzero Planck-peak solution.

    T=b/λmax⁡,b=2.897771955×10−3 m KT=b/\lambda_{\max},\quad b=2.897771955\times10^{-3}\ {\rm m\,K}
  3. Divide the displacement constant by the wavelength.

    T=5795.54 KT=5795.54\ {\rm K}

Interpretation. A Bν peak cannot be substituted into this wavelength formula.

↑ Return to definitions and contents
Planck spectral radiance per micrometre at three prescribed temperatures. These are blackbody curves, not observed stellar spectra; the vertical scale is absolute, not separately normalized.
Planck spectral radiance per micrometre at three prescribed temperatures. These are blackbody curves, not observed stellar spectra; the vertical scale is absolute, not separately normalized. Download SVG · Plot data (JSON)

Example 05. Radius from luminosity and temperature

Definitions & inputs. L=4L⊙, Teff=5772 K; σSB=5.670374419×10⁻⁸ SI.

Review the governing derivation ↑

  1. Isolate radius before substituting.

    L=4πR2σSBTeff4⇒R=L/(4πσSBTeff4)L=4\pi R^2\sigma_{\rm SB}T_{\rm eff}^4\Rightarrow R=\sqrt{L/(4\pi\sigma_{\rm SB}T_{\rm eff}^4)}
  2. Evaluate with the stated bolometric luminosity.

    R=4(3.828×1026)/(4πσSB57724)=1.3914e+09 mR=\sqrt{4(3.828\times10^{26})/(4\pi\sigma_{\rm SB}5772^4)}=1.3914e+09\ {\rm m}
  3. Compare with R⊙=6.957×10⁸ m.

    R/R⊙=2R/R_\odot=2

Interpretation. At fixed temperature, radius scales as the square root of luminosity.

↑ Return to definitions and contents

Example 06. Radial velocity from an absorption-line shift

Definitions & inputs. λ0=656.28 nm and λobs=656.50 nm.

Review the governing derivation ↑

  1. Compute the redward shift.

    Δλ=656.50−656.28=0.22 nm\Delta\lambda=656.50-656.28=0.22\ {\rm nm}
  2. Form a dimensionless fractional shift.

    z=Δλ/λ0=0.22/656.28z=\Delta\lambda/\lambda_0=0.22/656.28
  3. A positive value means recession under this convention.

    vr≃cz=100.497 km s−1v_r\simeq cz=100.497\ {\rm km\,s^{-1}}

Interpretation. Instrumental, gravitational, and atmospheric line shifts must be separated in precision work.

↑ Return to definitions and contents

Example 07. Total mass of a resolved binary

Definitions & inputs. Relative-orbit semimajor axis a=2 AU; period P=2 Julian years.

Review the governing derivation ↑

  1. Rearrange Kepler’s third law.

    M1+M2=4π2a3/(GP2)M_1+M_2=4\pi^2a^3/(GP^2)
  2. Use consistent SI units.

    a=2(1.495978707×1011) m,P=2(365.25)(86400) sa=2(1.495978707\times10^{11})\ {\rm m},\quad P=2(365.25)(86400)\ {\rm s}
  3. Evaluate and divide by the stated solar mass.

    Mtot=2.00002 M⊙M_{\rm tot}=2.00002\ M_\odot

Interpretation. The total mass does not by itself specify the individual masses.

↑ Return to definitions and contents

Example 08. Orbital scale for a ten-day planet

Definitions & inputs. P=10 days, M★=M⊙, Mp≪M★.

Review the governing derivation ↑

  1. Convert the measured repeat period.

    P=10(86400)=864000 sP=10(86400)=864000\ {\rm s}
  2. Solve Kepler’s law for semimajor axis.

    a=[GM⋆P2/(4π2)]1/3a=[GM_\star P^2/(4\pi^2)]^{1/3}
  3. Express the result in astronomical units.

    a=0.0908389 AUa=0.0908389\ {\rm AU}

Interpretation. An orbital period gives an orbital scale only after specifying the host mass.

↑ Return to definitions and contents

Example 09. Mean density of a solar-size star

Definitions & inputs. M=M⊙=1.98847×10³⁰ kg; R=R⊙=6.957×10⁸ m.

Review the governing derivation ↑

  1. Compute the sphere’s volume.

    V=4πR3/3V=4\pi R^3/3
  2. Divide total mass by volume.

    ρˉ=M/V=3M/(4πR3)\bar\rho=M/V=3M/(4\pi R^3)
  3. The interior density varies greatly even though the mean is modest.

    ρˉ=1409.82 kg m−3\bar\rho=1409.82\ {\rm kg\,m^{-3}}

Interpretation. Do not use the mean density as a central density.

↑ Return to definitions and contents

Example 10. Surface gravity of a solar-size star

Definitions & inputs. M=M⊙, R=R⊙.

Review the governing derivation ↑

  1. Cancel the test mass.

    Fg=GMm/R2⇒g=Fg/m=GM/R2F_g=GMm/R^2\Rightarrow g=F_g/m=GM/R^2
  2. Substitute the solar reference values.

    g=274.208 m s−2g=274.208\ {\rm m\,s^{-2}}
  3. Convert to cgs before computing the spectroscopic log g convention.

    log⁡10[g/(cm s−2)]=4.43808\log_{10}[g/({\rm cm\,s^{-2}})]=4.43808

Interpretation. log g is conventionally quoted using centimetres per second squared.

↑ Return to definitions and contents

Example 11. Central pressure of a uniform-density sphere

Definitions & inputs. Toy sphere: M=M⊙ and R=R⊙; set P(R)=0 and ρ constant.

Review the governing derivation ↑

  1. Insert the uniform enclosed mass into hydrostatic balance.

    Mr=4πρr3/3,dP/dr=−(4πGρ2/3)rM_r=4\pi\rho r^3/3,\quad dP/dr=-(4\pi G\rho^2/3)r
  2. Integrate inward from the zero-pressure surface.

    P(r)=2πGρ23(R2−r2)P(r)=\frac{2\pi G\rho^2}{3}(R^2-r^2)
  3. Replace density by 3M/(4πR³).

    Pc=2πGρ2R23=3GM28πR4P_c=\frac{2\pi G\rho^2R^2}{3}=\frac{3GM^2}{8\pi R^4}
  4. Evaluate the benchmark central pressure.

    Pc=1.34474e+14 PaP_c=1.34474e+14\ {\rm Pa}

Interpretation. This is not the actual solar central pressure; it is a useful hydrostatic solver test.

↑ Return to definitions and contents
Exact pressure profile for the constant-density Newtonian sphere with zero surface pressure. This is a benchmark, not a realistic solar model.
Exact pressure profile for the constant-density Newtonian sphere with zero surface pressure. This is a benchmark, not a realistic solar model. Download SVG · Plot data (JSON)

Example 12. Virial mean temperature of a toy star

Definitions & inputs. M=M⊙, R=R⊙, mean molecular weight μ=0.61.

Review the governing derivation ↑

  1. Express internal energy using the number of freely moving particles.

    2U+Ω=0,U=32(M/μmp)kB⟨T⟩2U+\Omega=0,\quad U=\tfrac32(M/\mu m_p)k_B\langle T\rangle
  2. Substitute the uniform-sphere binding energy and solve.

    Ω=−3GM2/(5R)⇒⟨T⟩=GMμmp/(5kBR)\Omega=-3GM^2/(5R)\Rightarrow\langle T\rangle=GM\mu m_p/(5k_BR)
  3. Evaluate the particle-weighted mean estimate.

    ⟨T⟩=2.81953e+06 K\langle T\rangle=2.81953e+06\ {\rm K}

Interpretation. This temperature is not a prediction of the Sun’s central temperature.

↑ Return to definitions and contents

Example 13. Gravitational contraction timescale

Definitions & inputs. M=M⊙, R=R⊙, L=L⊙.

Review the governing derivation ↑

  1. Divide the gravitational energy scale by the radiated power.

    Eg∼GM2/R,tKH∼Eg/LE_g\sim GM^2/R,\quad t_{\rm KH}\sim E_g/L
  2. Evaluate in seconds.

    tKH=9.90946e+14 st_{\rm KH}=9.90946e+14\ {\rm s}
  3. Convert to million Julian years.

    tKH=31.4012 Myrt_{\rm KH}=31.4012\ {\rm Myr}

Interpretation. For a uniform ideal-gas sphere, the available total binding energy introduces a 3/10 factor; this conventional scale deliberately omits it.

↑ Return to definitions and contents

Example 14. Hydrogen-burning energy-budget timescale

Definitions & inputs. M=M⊙, L=L⊙, burned hydrogen mass fraction fH=0.10 of total initial mass, η=0.007.

Review the governing derivation ↑

  1. Define the fuel fraction as hydrogen actually burned, not a separate core fraction to be multiplied again.

    MH,burned=fHM=0.10MM_{H,\rm burned}=f_HM=0.10M
  2. Use the approximate hydrogen-to-helium conversion efficiency.

    Efuel=ηfHMc2=0.0007Mc2E_{\rm fuel}=\eta f_HMc^2=0.0007Mc^2
  3. Convert the resulting energy-budget time to billions of years.

    tfuel=Efuel/L=10.3558 Gyrt_{\rm fuel}=E_{\rm fuel}/L=10.3558\ {\rm Gyr}

Interpretation. This is a fuel-budget estimate, not a stellar-evolution track.

↑ Return to definitions and contents

Example 15. Eddington scale for ten solar masses

Definitions & inputs. M=10M⊙ and prescribed opacity κ=0.034 m² kg⁻¹.

Review the governing derivation ↑

  1. Write the outward and inward accelerations.

    arad=κL/(4πr2c),ag=GM/r2a_{\rm rad}=\kappa L/(4\pi r^2c),\quad a_g=GM/r^2
  2. Cancel radius and solve for luminosity.

    arad=ag⇒LEdd=4πGMc/κa_{\rm rad}=a_g\Rightarrow L_{\rm Edd}=4\pi GMc/\kappa
  3. Evaluate the force-balance luminosity.

    LEdd=1.47054e+32 W=384153L⊙L_{\rm Edd}=1.47054e+32\ {\rm W}=384153L_\odot

Interpretation. Opacity and emission geometry change the applicable luminosity scale.

↑ Return to definitions and contents

Example 16. Transmission through optical depth two

Definitions & inputs. Incident intensity I0; optical depth τ=2.0; source function S=0.

Review the governing derivation ↑

  1. Separate and integrate the transfer equation.

    dI/dτ=−I⇒∫I0IdI′/I′=−τdI/d\tau=-I\Rightarrow\int_{I_0}^{I}dI^\prime/I^\prime=-\tau
  2. Only about 13.5% is transmitted.

    I/I0=e−2=0.135335I/I_0=e^{-2}=0.135335
  3. Convert the same attenuation into magnitudes.

    A=−2.5log⁡10(I/I0)=2.5ln⁡10τ=2.17147 magA=-2.5\log_{10}(I/I_0)=\frac{2.5}{\ln10}\tau=2.17147\ {\rm mag}

Interpretation. A nonzero source function or scattered light would change the observed intensity.

↑ Return to definitions and contents

Example 17. Free-fall time of a molecular-cloud parcel

Definitions & inputs. Total particle number density n=10⁹ m⁻³ (10³ cm⁻³), mean particle mass μmp with μ=2.33.

Review the governing derivation ↑

  1. Convert number density to mass density.

    ρ=μmpn=2.33mp(109)\rho=\mu m_pn=2.33m_p(10^9)
  2. Use the integrated pressureless collapse time.

    tff=3π/(32Gρ)t_{\rm ff}=\sqrt{3\pi/(32G\rho)}
  3. Convert seconds to million years.

    tff=1.06629 Myrt_{\rm ff}=1.06629\ {\rm Myr}

Interpretation. Pressure, magnetic fields, or turbulence can delay or prevent collapse.

↑ Return to definitions and contents

Example 18. Jeans mass of a ten-kelvin gas cloud

Definitions & inputs. T=10 K, μ=2.33, n=10⁹ m⁻³; ρ=μmpn.

Review the governing derivation ↑

  1. Calculate thermal sound speed, not the adiabatic value with an extra γ factor.

    cs=kBT/(μmp)=188.22 m s−1c_s=\sqrt{k_BT/(\mu m_p)}=188.22\ {\rm m\,s^{-1}}
  2. Set ω²=0 in the dispersion relation.

    λJ=csπ/(Gρ)=0.670362 pc\lambda_J=c_s\sqrt{\pi/(G\rho)}=0.670362\ {\rm pc}
  3. Compute the mass using the stated radius convention.

    MJ=4π3ρ(λJ/2)3=9.08267 M⊙M_J=\frac{4\pi}3\rho(\lambda_J/2)^3=9.08267\ M_\odot

Interpretation. Jeans mass scales as T^(3/2)ρ^(−1/2) at fixed composition; it is not a universal observed core mass.

↑ Return to definitions and contents
Isothermal Jeans mass at total particle density 10³ cm⁻³ and μ=2.33. Mass is measured inside radius λJ/2; magnetic and turbulent support are omitted.
Isothermal Jeans mass at total particle density 10³ cm⁻³ and μ=2.33. Mass is measured inside radius λJ/2; magnetic and turbulent support are omitted. Download SVG · Plot data (JSON)

Example 19. Planet radius from a one-percent transit

Definitions & inputs. Transit depth δ=0.0100, stellar radius R★=R⊙.

Review the governing derivation ↑

  1. Take the square root of the area ratio.

    δ=(Rp/R⋆)2⇒Rp/R⋆=0.0100=0.100\delta=(R_p/R_\star)^2\Rightarrow R_p/R_\star=\sqrt{0.0100}=0.100
  2. Convert the inferred radius to kilometres.

    Rp=0.100R⊙=69570 kmR_p=0.100R_\odot=69570\ {\rm km}
  3. For independent small uncertainties, propagate logarithmic errors; the stellar radius can dominate.

    σRpRp≃(σR⋆/R⋆)2+14(σδ/δ)2\frac{\sigma_{R_p}}{R_p}\simeq\sqrt{(\sigma_{R_\star}/R_\star)^2+\tfrac14(\sigma_\delta/\delta)^2}

Interpretation. A one-percent drop does not imply a one-percent radius ratio.

↑ Return to definitions and contents

Example 20. Enclosed mass from a circular speed

Definitions & inputs. v=200 km s⁻¹ at r=10 kpc.

Review the governing derivation ↑

  1. Convert the measured scales.

    v=2.00×105 m s−1,r=104 pcv=2.00\times10^5\ {\rm m\,s^{-1}},\quad r=10^4\ {\rm pc}
  2. Use circular force balance.

    M(<r)=v2r/GM(<r)=v^2r/G
  3. Evaluate the dynamical mass within the specified radius.

    M(<r)=9.30006e+10 M⊙M(<r)=9.30006e+10\ M_\odot

Interpretation. This estimates total gravitating mass under the model; separating stars, gas, and dark matter requires independent information.

↑ Return to definitions and contents

Astrophysics brings the physics chapters together

Numerical Modeling supplies methods for stellar ODE boundary-value problems, radiative transfer, gravitational N-body motion, and fluid or magnetohydrodynamic PDEs. The analytical examples here provide useful checks on those calculations.

Return to the complete physics pathway →

Notation used throughout

n: number density (m⁻³); N: particle count; ρ: mass density (kg m⁻³); ρc: charge density; p: pressure; T: temperature (K); kB: Boltzmann constant; h, ℏ: Planck constants; β=1/(kBT); μ: chemical potential; f: phase-space distribution; g(r): pair distribution. In solid displacements u is a displacement; in the liquid closure u(r) is pair energy; in fluid equations u is bulk velocity. Subscripts identify phase or species. Every approximation must use consistent SI units or explicitly stated reduced units.