Interstellar medium (ISM)

The interstellar medium is the gas, dust, magnetic field, and energetic-particle population between stars in a galaxy.

Interstellar gas occupies distinct thermal and ionization regimes because heating, cooling and compression depend on density, radiation and composition. Cold and warm neutral gas can coexist over a finite pressure range, while ionized gas, molecular clouds and hot shock-heated regions add further components. turbulence, magnetic stress and finite cooling times mean that these are dynamically exchanging phases, not necessarily a common global equilibrium.

A useful optically thin model subtracts volumetric heating from radiative cooling. At equilibrium its net cooling vanishes. Pressure-confined temperature perturbations are governed by the derivative of net cooling per unit mass at fixed pressure, so a negative derivative is unstable. The heating law, species abundances and radiation field must be specified before predicting an equilibrium curve. thermal conduction and finite sound-crossing times limit the isobaric approximation.

Dilute shock-heated interstellar plasma with characteristic temperatures around $10^6\mathrm{K}$. supernova energy injection and subsequent cooling are important to its evolution; temperature alone does not determine its volume filling factor.

Diffuse ionized gas near $10^4\mathrm{K}$, maintained by ionization and heating balanced by radiative losses. It is physically distinct from warm gas that remains mostly neutral despite a similar kinetic temperature.

Diffuse mostly neutral atomic gas with temperatures below about $10^4\mathrm{K}$. It can share a pressure range with the cold neutral medium where the thermal equilibrium curve has two stable branches.

The cold predominantly atomic neutral phase of interstellar gas, with characteristic temperatures of order $10^2\mathrm{K}$. It is denser than the warm neutral phase at comparable thermal pressure. Its existence depends on radiation, grain heating and atomic-line cooling; one universal density or temperature does not define every environment.

A superbubble is an expanding cavity of hot gas driven by the combined energy input of stellar winds and multiple supernova explosions. Constant mechanical power and negligible radiative losses lead to a radius proportional to $t^{3/5}$ in a uniform ambient medium.

For constant injected power $\dot{E}$ and ambient mass density ${\rho }_1$, dimensional analysis gives $R=\alpha (\dot{E}/{\rho }_1{)}^{1/5}{t}^{3/5}$. In a similarity solution, set $\eta =r/R$, $\rho ={\rho }_{1}f(\eta )$, $p={\rho }_1{\dot{R}}^{2}g(\eta )$ and $u=\dot{R}h(\eta )$. The strong-shock values on the inner side of $\eta =1$ are $f=(\gamma +1)/(\gamma -1)$ and $g=h=2/(\gamma +1)$.

For point injection of constant power, the outward energy flux through a sphere tending to the source must tend to $\dot{E}$. In the adiabatic superbubble similarity solution, this requiresImposing finite regular central pressure and density with $h=O(\eta )$ would instead give zero flux and fail to supply the prescribed power. A resolved source region or an inner contact-boundary model can replace the point-source condition.

If the similarity profiles account for all kinetic and internal energy inside an adiabatic strong shock, letThe injected energy $\dot{E}t$ equals $4\pi {\rho }_1{\dot{R}}^2{R}^{3}I$. Consequently conservation of energy fixesAn interior region not included in the profiles must have its energy added separately; a radiatively cooled shell changes this energy budget.

Substitution of the adiabatic superbubble similarity solution into spherical mass conservation gives

A supernova remnant is the expanding ejecta and shocked ambient medium left by a supernova explosion.

In the ideal cooled-shell limit with negligible continuing interior pressure, swept-up mass in a uniform medium scales as $R^3$. Constant radial shell momentum then gives $R^3\dot{R}=\mathrm{const}$ and $R\propto t^{1/4}$. This differs from an energy-conserving blast wave or a pressure-driven snowplough.

A self-similar energy-conserving strong spherical explosion into a uniform ambient density has $R_s=\xi (E{t}^2/{\rho }_0{)}^{1/5}$ and $v_s=2R_s/(5t)$. The dimensionless coefficient depends on the equation of state. Radiative losses, nonuniform media or continuous energy injection alter the similarity solution.

A monatomic strong shock has $T_d=3\bar{m}{V}_s^2/(16k_B)$. Together with the Sedov-Taylor blast wave radius this gives $T_d\propto t^{-6/5}$. Extrapolating to ambient temperature or sonic speed estimates when the cold-upstream assumption fails, but the strong-shock law is already becoming inaccurate there. Radiative cooling can end the energy-conserving stage earlier; its transition age additionally needs a cooling function and composition.

During its adiabatic phase, a supernova remnant loses little energy radiatively and expands as a hot, broad shocked flow rather than a thin cooled shell.

A molecular cloud is a cold dense region of the interstellar medium in which hydrogen is predominantly molecular and stars can form.

Rotational carbon-monoxide line emission traces cold molecular clouds because abundant molecular hydrogen is difficult to observe directly. A conversion factor $X_{\mathrm{CO}}$ relates integrated CO intensity to molecular-hydrogen column density.

For one optically thin LTE carbon-monoxide component, $L_{ul}=N_u{A}_{ul}h{\nu }_{ul}$ and $N_J\propto (2J+1)e^{-E_J/(k_{B}T)}$. Thus $C=7A_{32}{\nu }_{32}/(3A_{10}{\nu }_{10})$ for the 3–2/1–0 luminosity ratio. The temperature then determines the $J=1$ population fraction and hence the CO mass. Multiple excitation components, optical depth or subthermal excitation invalidate a single-temperature inference.

Interstellar dust absorbs and scatters short-wavelength starlight and reradiates the absorbed energy thermally in the infrared.

Interstellar extinction is the attenuation of light by absorption and scattering from dust between a source and observer. Its wavelength dependence commonly reddens a source as well as dimming it.

Two optically thin lines from the same upper level have an intrinsic energy-emissivity ratio fixed by their radiative branching rates and photon energies. The upper population cancels, so a common electron-temperature or abundance change cannot alter that ratio. Comparing an optical branch with a weakly attenuated infrared branch can determine differential dust extinction, subject to matched emitting volumes and calibrated emissivities.

The observed colour minus the inferred intrinsic colour in specified photometric passbands; for B and V it is $E(B-V)=(B-V{)}_{\mathrm{obs}}-(B-V{)}_{\mathrm{intrinsic}}$.

An H II region is a volume of ionized hydrogen surrounding one or more sources of ultraviolet photons. Photoionization heating, recombination, gas dynamics, and dust determine its structure and emission spectrum.

An ionization front is the transition between predominantly ionized and predominantly neutral gas. Its motion and thickness depend on photon supply, recombination, gas density, and hydrodynamic response.

In a homogeneous initially neutral medium, an isotropic source with hydrogen-ionizing photon production rate $Q(t)=Q_0{e}^{t/\tau }$ ionizes the displayed comoving inventory when every escaping photon produces one primary ionization and recombinations are absent. The physical front radius is $a(t)X$ in an expanding background. Retarded photon propagation prevents any instantaneous front from travelling faster than light. This photon-budget front is not the equilibrium Strömgren radius, which balances emission against recombinations.

Photoionization equilibrium holds when the local or integrated photoionization rate equals the radiative recombination rate, so the ionized fraction is stationary.

For pure hydrogen illuminated by an isotropic central source with hydrogen-ionizing photon production rate $Q$, an empty cavity of physical radius $r_{\mathrm{in}}$ gives photoionization rate $q=Q⟨\sigma ⟩/(4\pi r_{\mathrm{in}}^2)$ per neutral atom. Local photoionization equilibrium is $(1-\zeta )q={\alpha }_B{n}_H{\zeta }^2$, so $\Theta ={\alpha }_B{n}_H/q$ and the physical quadratic equation root is the displayed rationalized expression. If an unextinguished total recombination-line observed radiative flux $F_{\ell }$ has photon yield $y_{\ell }$, then $Q=4\pi D^2{F}_{\ell }{\lambda }_{\ell }/(y_{\ell }hc)$ and $\Theta =y_{\ell }{\alpha }_B{n}_{H}hc(r_{\mathrm{in}}/D{)}^2/(⟨\sigma ⟩F_{\ell }{\lambda }_{\ell })$. The ratio $r_{\mathrm{in}}/D$, rather than a product of lengths, makes $\Theta$ dimensionless.

When the number of photoionizations equals the number of recombinations, and only those processes exchange thermal energy, equilibrium requires $⟨h\nu -h{\nu }_0{⟩}_{\mathrm{abs}}=⟨E_{\mathrm{rec}}(T)⟩$. Photon absorption and electron capture have their own spectral weights. The density cancels from this event balance, but the actual equilibrium temperature requires the two energy averages.^[ref]

An ionization-bounded nebula contains enough absorbing material to trap the relevant ionizing photons, unlike a density-bounded region that ends before its ionization front. In dust-free hydrogen at steady photoionization equilibrium, the absorbed stellar photon rate equals the integrated Case B recombination rate. Escape, absorption by dust, and competition from other species must be accounted for before identifying the absorbed hydrogen rate with the entire stellar ionizing luminosity.

For fully ionized pure hydrogen of density $n(r)=n_0(r/r_0{)}^{-\beta }$ and constant Case B recombination coefficient $\alpha$, the total recombination rate isThe central integral converges only for $\beta <3/2$; an envelope extending to infinity has finite recombination capacity only for $\beta >3/2$. At $\beta =3/2$, each logarithmic radius interval contributes equally. A cored power-law density profile removes the central singularity.

For the cored power-law density profile that is constant inside $r_c$ and falls as $r^{-\beta }$ outside, with $\beta >3/2$, the maximum integrated Case B recombination rate isA hydrogen-ionizing photon production rate greater than this has no finite radius at photoionization equilibrium. Equality requires an infinite ionized radius. Other core shapes change the coefficient but retain the convergence condition $\beta >3/2$.

The recombination time is the characteristic time for an ionized plasma to recombine after its ionizing source is removed. Uniform hydrogen with electron density $n_e$ and recombination coefficient $\alpha$ has $t_{\mathrm{rec}}=1/(\alpha n_e)$.

Case B recombination excludes direct recombinations to the ground state because their ionizing photons are assumed to be absorbed locally. The coefficient ${\alpha }_B$ therefore counts recombinations to excited states.

The hydrogen-ionizing photon production rate is the number of photons emitted per unit time above the hydrogen ionization threshold:

A spherical blackbody star emits $Q(E_0)=8{\pi }^2{R}_s^2(k{T}_s{)}^{3}F(E_0/k{T}_s)/(h^3{c}^2)$ photons per time above threshold $E_0$. Integrating specific intensity over the outward hemisphere gives flux $\pi B_{\nu }$; dividing the resulting spectral luminosity by photon energy and substituting $x=h\nu /k{T}_s$ proves the formula. The convergent series $F(y)={\sum }_{m\ge 1}{e}^{-my}(y^2/m+2y/m^2+2/m^3)$ follows by integrating the geometric-series expansion termwise.

The ionizing photon escape fraction is the fraction of photons capable of photoionization that leave their source environment. Only these escaping photons contribute to ionizing distant intergalactic medium. Its value can depend on frequency, source geometry, and time.

For $L_{\nu }=A{\nu }^{-s}$ above an ionization threshold ${\nu }_0$, the energy luminosity is $A{\nu }_0^{1-s}/(s-1)$ and the photon rate is $A{\nu }_0^{-s}/(sh)$. Their ratio is the displayed emitted mean energy. For slope two it is twice the threshold energy. Absorption weighting, a high-energy cutoff and secondary ionizations can change the energy deposited per ionization; emitted and locally absorbed spectra must not be conflated.

A Strömgren sphere is the idealized ionized region around a steady source embedded in uniform hydrogen gas, bounded where integrated recombinations consume the source's ionizing-photon supply.

In a uniform sharp-front approximation, an inner ${\mathrm{He}}^{++}$ zone of radius $R_2$ lies inside a ${\mathrm{He}}^{+}$ shell ending at $R_1$, within the hydrogen ionization front $R_H$. For nuclear densities $n_H,n_{He}$, the electron densities are $n_H+2n_{He}$, $n_H+n_{He}$ and $n_H$ in the three zones. Thus the hydrogen recombination integral is $(4\pi /3){\beta }_H{n}_H[n_H{R}_H^3+n_{He}(R_1^3+R_2^3)]$. The helium integrals are $(4\pi /3){\beta }_{HeI}{n}_{He}(n_H+n_{He})(R_1^3-R_2^3)$ and $(4\pi /3){\beta }_{HeII}{n}_{He}(n_H+2n_{He})R_2^3$. These are balanced against species-specific effective photon budgets, not three independently reused cumulative stellar photon counts.

For ionizing-photon rate $Q_{\mathrm{H}}$, uniform pure-hydrogen density $n$, and Case B recombination coefficient ${\alpha }_B$, the Strömgren radius is

Before hydrodynamic expansion matters, photon counting in a uniform initially neutral medium givesThe ionization front approaches the Strömgren radius on the recombination time.

A dusty Strömgren sphere includes ionizing photons absorbed by dust before they can ionize hydrogen. Dust therefore decreases the ionized radius below the dust-free Strömgren radius at fixed source luminosity and gas density.

The neutral-hydrogen 21-centimeter line is the ground-state hyperfine transition of atomic hydrogen. Its brightness and Doppler shift map atomic-gas column density and line-of-sight velocity.

A cosmic ray is a high-energy charged particle propagating through space. Galactic cosmic-ray electrons emit synchrotron and inverse-Compton radiation, while cosmic-ray nuclei produce gamma rays through collisions with gas.

In modern terminology, a primary cosmic ray arriving at Earth's atmosphere is usually a proton or another atomic nucleus, with an electron component as well. Atmospheric collisions generate secondary cosmic rays and air showers, including muons, electrons and photons. A gamma ray can also initiate an air shower, although it is electrically neutral. Consequently, detecting atmospheric ionization does not establish that an incoming primary cosmic ray passed through the detector. Magnetic fields bend charged particles, while gamma rays provide a more direct directional connection to their sources.^[ref]

The Hess cosmic-ray balloon experiment asks how atmospheric shielding and separation from the ground affect ionization. Its power comes from changing the environment around an enclosed detector instead of identifying each incoming particle.

Fermi acceleration transfers bulk-flow energy to fast charged particles through encounters with moving magnetic scattering structures. Scattering is approximately elastic in each structure's rest frame; transformations between these frames produce a systematic energy change.

In first-order Fermi acceleration, converging scattering media produce a mean fractional energy gain linear in their relative speed divided by the particle speed. At a nonrelativistic shock wave, an almost isotropic ultrarelativistic population gains $2(u_1-u_2)/(3c)$ per crossing and $4(u_1-u_2)/(3c)$ per complete return cycle. The distinction between a crossing and a cycle must also be respected when assigning an escape probability.

Diffusive shock acceleration repeatedly returns particles to a shock wave by scattering them in the upstream and downstream plasma. For fluid speeds $u_1>u_2$, particle speed $v\gg u_1$, and shock compression ratio $r=u_1/u_2$, the cycle gain is $⟨\Delta p⟩/p=4(u_1-u_2)/(3v)$ and the cycle escape probability is $4u_2/v$. Their ratio gives $f(p)\propto p^{-3r/(r-1)}$ for the isotropic phase-space one-particle distribution function. For $r=4$ this is $p^{-4}$, corresponding to $dN/dE\propto E^{-2}$ in the ultrarelativistic regime.

For an ultrarelativistic proton of energy $E$ colliding head-on with a photon of energy $ϵ$, the invariant squared collision energy is approximately $(m_p{c}^2{)}^2+4Eϵ$. Producing a final proton and neutral pion requires at least $(m_p{c}^2+m_{\pi }{c}^2{)}^2$. Thus $E_{\mathrm{th}}=m_{\pi }{c}^2(2m_p{c}^2+m_{\pi }{c}^2)/(4ϵ)$. A broad isotropic photon background turns this threshold into an energy-loss feature depending on the photon distribution and source distance, not an absolute ceiling for all arriving protons.

Here a primary cosmic ray means a particle arriving from space before its first atmospheric collision.

A secondary cosmic ray here is a particle produced by an atmospheric interaction of a primary cosmic ray or its descendants. A muon at ground level need not have traveled from the original astronomical source.

An air shower is a cascade of secondary particles produced by an energetic particle interacting in the atmosphere. Both a charged particle and a gamma ray can start one.^[ref]

Victor Franz Hess's 1912 series used two sealed, 3-mm-walled ionization chambers with Theodor Wulf's fiber electrometers. Karl Bergwitz suggested internal plating; Günther & Tegetmeyer improved focusing. Basket-mounted ground comparisons, insulation-loss measurements, simultaneous detectors and barograph-derived heights supplied checks. The highest flight reached 5350 m; its comparisons came from Vienna and after landing, without a launch-site baseline.^[ref]

A third, thin-walled chamber admitted ambient air. Hess multiplied its readings by $750/b$, with pressure $b$ in millimeters of mercury, warning that vessel-wall contributions might violate this scaling. He recorded outside temperature but applied no absolute-temperature correction. Night flights, an eclipse and weather observations constrained alternatives to radiation from above; the eclipse argument addressed directly propagating solar radiation.^[ref]

The electrometers read gradual electric potential loss caused by continuing air ionization. The low-altitude dip and strong high-altitude rise favored an extraterrestrial contribution, without measuring proton composition.^[ref] An ionization chamber integrates electric charge: one energetic particle may create many ion pairs, several particles may contribute to one reading, and a photon can create charged particles in the gas or walls. A voltage-loss reading is therefore neither a count of primary cosmic rays nor a chemical analysis of them.

Domenico Pacini supplied independent evidence: in June 1911 he compared a copper-box-enclosed electroscope at the sea surface and 3 m underwater, more than 300 m offshore above an 8-m-deep seabed. Reduced underwater ionization, stable electrical loss between runs, and a lake comparison supported an external contribution independent of soil radioactive decay. Water provided shielding; this established a different constraint from Hess's altitude comparison.^[ref]

An elementary ionization-chamber charge-loss model makes the inference explicit. Let a positively charged electrode and its electrometer have effective capacitance $C$, electric potential difference $U$, and enclosed gas volume $\mathcal{V}$. Assume constant $C$, a short enough interval that collection conditions remain stable, and no multiplication of charge in the gas. If ionization produces $q$ ion pairs per unit volume per unit time and a fraction $\eta$ is collected, thenHere $e$ is the elementary charge and $I_{\mathrm{leak}}$ is the parasitic electric current through imperfect insulation. Each fully collected ion pair transfers one elementary charge through the measuring circuit; counting its two opposite carriers as two transfers would introduce a factor-of-two error. The relation follows from charge conservation, rather than from a model of the incoming particle's identity.

For an invented example, take $C=2.00\mathrm{pF}$, $\mathcal{V}=2.00\mathrm{L}$, $\eta =1$, a measured fall of $6.00\mathrm{V}{\mathrm{h}}^{-1}$, and a separately determined leakage-equivalent fall of $0.60\mathrm{V}{\mathrm{h}}^{-1}$. The leakage-subtracted fall is $5.40\mathrm{V}{\mathrm{h}}^{-1}$, soThis is about $1.87\times 10^4$ collected ion pairs per second throughout the chamber, not that many primary cosmic rays. Unknown collection losses make $\eta <1$ another uncertainty; changes in capacitance would require differentiating $Q=CU$ with both factors variable.

Gas density is a separate issue. Under a fixed radiation field, fixed gas composition and approximately constant energy deposited per unit gas mass, the volume-based ionization rate scales as $q\propto \rho$. For a vented ionization chamber, the ideal gas law gives the illustrative correctionwhere $T$ is absolute temperature. A pressure-only correction presumes equal temperature. A sealed rigid chamber instead fixes gas number density as it rises: its internal pressure may change with temperature without losing gas. This does not make its insulation, electrometer calibration or collection efficiency insensitive to temperature. Nor can dividing by gas density correct arbitrary radiation from the vessel walls or changing particle energies. Those are distinct effects, not a universal normalization recipe.

An original cosmic-ray altitude attenuation model illustrates how competing origins can produce a dip followed by a rise. It is a hypothetical absorption calculation, not a reconstruction or fit of Hess's measurements, and it excludes air shower production. Put $z$ for height above a flat ground surface and assume gas density $\rho (z)={\rho }_0{e}^{-z/H}$ with constant scale $H$. The overlying atmospheric mass per unit area isLet $G$ denote a ground-origin contribution attenuated by material below the chamber, and $E$ a downward extraterrestrial contribution attenuated by material above it. Assume constant attenuation mass columns ${\Lambda }_g,{\Lambda }_e>0$:These are rates at a fixed detector gas density, in arbitrary common units. The constant $B$ represents residual background. The equation $dI/dX=-I/\Lambda$ gives exponential function attenuation along each propagation direction. Thus $G$ decreases and $E$ increases as $z$ increases; their sum need not be monotone. The different arguments of the two exponential functions encode the entire geometric idea.

The figure chooses $H=8\mathrm{km}$, $G_0=0.8$, $E_{\mathrm{top}}=6$, $B=0.15$, ${\Lambda }_g=0.08X_0$ and ${\Lambda }_e=0.45X_0$, all for illustration. Differentiation givesThe minimum occurs where $E/{\Lambda }_e=G/{\Lambda }_g$, at approximately $1.13\mathrm{km}$ for these invented parameters. At $6\mathrm{km}$ the fixed-gas response is about $1.41$ times its ground value, but multiplying by $\rho (z)/{\rho }_0$ gives an uncorrected vented-chamber response of only $0.66$ times its ground value. A real increase in radiation response per gas mass can therefore be hidden by a decrease in target gas density. This is the reason to distinguish radiation from the amount of gas being irradiated.

The model treats propagation as attenuation without generating new particles. Real primary cosmic rays interact, redistribute energy and produce air showers; an atmospheric profile need not follow either simple exponential function. The model also holds backgrounds and detector response fixed. Its value is to expose the sign of competing effects and the density confound, while Hess's experimental argument depends on his actual comparisons and instrument limitations.