Astrophysical cooling function Created 2026-10-06 Updated 2026-10-07
An astrophysical cooling function packages radiative losses of a plasma into a temperature-dependent coefficient. In the hydrogen-density-squared convention, the volume energy loss is , so the coefficient has units . Its dependence on ionization state, chemical abundances and radiation field must be specified; a low-density collisional-equilibrium curve is not universal.
The Mach number is , comparing the flow speed relative to a disturbance with the local sound speed. For a body moving supersonically, spherical acoustic disturbances emitted earlier have a tangent envelope with Mach cone half-angle . Signals cannot propagate arbitrarily far upstream. Compressive characteristics can converge and steepen into a shock wave, within which microscopic dissipation prevents a multivalued fluid state; supersonic motion alone does not require every flow to contain a shock.
Work in the shock frame with steady normal velocities . Integrating conservation of mass, normal momentum and total energy over a thin pillbox, neglecting gravity across the thin layer, gives the Rankine-Hugoniot conditions for a perfect gas:
The energy equation includes pressure work through specific enthalpy. The shock is adiabatic if negligible heat escapes the layer, but its irreversible dissipation means it is not isentropic.
Put and . Momentum conservation gives . Substitution in the energy equation produces the unshocked root and the compressive root
Since , the strong-shock limit has negligible upstream thermal pressure compared with ram pressure. The Strong-shock Rankine-Hugoniot conditions give
For unchanged composition, the ideal gas law implies , so
For , the compression tends to four and , where is mean particle mass.
With constant specific heat , entropy production in a perfect-gas shock is measured by the change in specific entropy
The strong-shock expression is therefore
The exact jump is positive for every compressive shock with , in accordance with the Second law of thermodynamics; the asymptotic expression is only for . There is a normalization issue in the printed description: this logarithmic formula with constant describes entropy per unit mass. The true entropy density is , so its difference is , rather than the same logarithmic jump with a constant heat capacity per unit volume. A material parcel's irreversible entropy increase should be assessed using .
Downstream gas is denser and hotter. For optically thin collisional emission the loss per unit volume is typically , whereas stored thermal energy scales approximately as . Compression can therefore reduce the radiative gas cooling time; excitation and ionization at the higher temperature can open additional cooling channels. The actual efficiency still depends on the astrophysical cooling function, composition and optical depth. In a radiative shock with separated scales, the temperature first rises across a narrow dissipative layer, then decreases through a much broader cooling region of characteristic length . Subsonic cooling usually brings further compression with approximately constant pressure. The immediate microscopic shock thickness and the macroscopic cooling length are distinct.
Let state be the fully cooled gas, with and the same composition. Set , the isothermal sound speed squared. Across the entire shock plus cooling layer, mass and momentum are still conserved, but gas energy alone is not. With ,
The nontrivial isothermal shock therefore has
If the upstream Mach number were defined with instead, this would be . Retaining the distinction avoids losing the factor .
The gas's specific entropy rises irreversibly at the hot shock and then decreases as escaping radiation removes heat. Since the final temperature equals the initial temperature,
This decrease for the gas is consistent with the Second law of thermodynamics: the emitted radiation and its eventual surroundings carry entropy away. It does not make the microscopic compressive shock reversible. The figure illustrates both the hot entropy jump and the subsequent cooling; its spatial cooling profile is schematic because no cooling law has been specified.
Figure 1.
Temperature and specific entropy across an adiabatic shock followed by radiative cooling
. Details.
For a monatomic ideal gas, write the thermal energy density as . At low optical depth, radiative loss uses an astrophysical cooling function, often written , while contains heating. At fixed density the net radiative gas cooling time is approximately when the denominator is positive; a cooling flow closer to constant pressure must include enthalpy and compressional work. Comparing this time with the free-fall time of a uniform sphere determines whether pressure can be lost during gravitational contraction.
Hot ionized gas cools by bremsstrahlung and recombination; collisional excitation followed by photon escape supplies atomic and metal-line cooling at lower temperatures. molecular line cooling through rotational or vibrational lines and dust thermal emission become important in cold dense gas. Abundances, ionization state, density-dependent collisional de-excitation, and photon trapping all change the cooling function. Optically thick gas requires radiative transfer, not a simple local density-squared loss law.
Heating comes from photoionization, grain photoelectric heating driven by ultraviolet light, cosmic rays, shocks, turbulent dissipation and compression. photoionization heating deposits the photon energy above the ionization threshold as electron kinetic energy. A thermal equilibrium satisfies ; its stability depends on how the net loss changes under the relevant pressure or density constraint, not merely on equality. Different stable balances permit a multiphase interstellar medium. Spectral-line ratios constrain temperature and ionization, but interpreting them requires density diagnostics, interstellar extinction corrections and a radiation-field model.
Galaxy formation combines gravitational assembly of dark-matter haloes with dissipative baryonic physics. Primordial density fluctuations supply a statistical seed distribution. In linear theory, each Fourier mode evolves through a growth factor and a scale-dependent transfer function; the power spectrum describes the variance by scale. Radiation pressure, free streaming and the time of horizon entry influence the small-scale transfer, so one should not replace the spectrum by a single universal growing power of time.
When fluctuations become nonlinear, overdense regions collapse, virialize and build larger structures through accretion and mergers. Cold dark matter permits a broad hierarchy of bound haloes. A spherical-collapse threshold and a mass-variance prescription, such as the Press-Schechter formalism, organize an approximate halo abundance; departures from spherical collapse and correlated environments limit its literal interpretation.
Gas falling into a halo is compressed and may be shock-heated toward a virial scale
Unlike collisionless dark matter, it can radiate away thermal energy. A rough cooling time is
where the denominator specifies a chosen optically thin astrophysical cooling function convention. Cooling is effective for condensation when this time is short compared with the relevant dynamical or supply time. Atomic transitions, molecular processes, ionization state, density and metal abundance matter. A stable virial shock is not obligatory for all accretion: cold dense inflows can supply gas in regimes where shock-heated gas would cool rapidly.
The distinction between a gravitationally assembled massive halo and a radiatively condensing luminous component is fundamental. In the classic two-stage picture, non-dissipative clustering builds potential wells while cooling allows gas to concentrate much further than the halo component.[ref]
Angular momentum inherited from tidal torques and later accretion helps gas settle into a rotationally supported disc after cooling; disc size depends on angular-momentum retention and transport. Locally self-gravitating gas fragments and forms stars. Stellar radiation, winds and supernovae then heat gas, drive turbulence, enrich it with heavy elements and can expel it from shallow potential wells. Metal enrichment can also accelerate later cooling, so feedback is more than a single positive or negative energy term.
Galaxy growth is usefully expressed as a gas-budget competition:
Here is the promptly returned stellar mass fraction; unlike the simplification in Question1, a general evolutionary model need not set it to zero. The equation separates newly supplied, converted, expelled and recycled gas rather than assuming that halo mass uniquely fixes stellar mass.
Accreting central black holes can provide radiative and mechanical feedback. Heating a halo, ejecting gas or maintaining a hot atmosphere can reduce fresh fuel; removal of a satellite's gas supply and ram pressure in a surrounding medium provide additional environmental routes. Quiescent galaxy assembly can still add stellar mass through dry mergers, so cessation of star formation is not identical to cessation of all mass assembly.
Mergers and disc instabilities can rearrange morphology and trigger bursts, while smooth accretion can rebuild discs. Cosmological expansion changes the supply and merger environment; accelerated late-time expansion also slows linear growth. A predictive model must connect these mechanisms to the stellar mass function, gas fractions, star-formation histories, sizes, chemical abundances and clustering. Those observables test different pieces of the assembly-and-feedback picture; none alone implies a one-to-one instantaneous relation between galaxy luminosity and total halo mass.
For an optically thin low-density gas in collisional ionization equilibrium, the astrophysical cooling function collects losses per pair of hydrogen nuclei into . Its units are ; the time unit is missing from the printed guide to the vertical-axis labels. For the primordial atomic cooling curve, only hydrogen and helium provide bound-state coolants.
Figure 1.
Schematic primordial cooling function with hydrogen and helium features, metal enrichment and photoionization effects
. Details.
Near , collisional excitation of hydrogen becomes effective and subsequent photon emission removes thermal energy. The excitation rate is exponentially suppressed below the atomic energy threshold. Atomic line cooling and collisional ionization generate a strong hydrogen feature around a few ; helium excitation and ionization produce further structure around . Radiative recombination also contributes. Once the gas is highly ionized, these bound-state losses weaken. At –, thermal bremsstrahlung dominates, with approximately apart from slowly varying factors. The figure is an original qualitative sketch anchored to the two supplied values, not an atomic-rate calculation.
Increasing galactic metallicity adds many ionic transitions, producing metal-line cooling and substantially enhancing cooling, particularly in the – range. Metals also provide low-energy fine-structure transitions below . The change is not a uniform vertical shift: the positions and strengths of features depend on the elemental abundances and ionization state.
Photoionization removes bound electrons even where collisions alone would leave atoms neutral, often suppressing the hydrogen and helium line peaks and changing the metal-ion population. It also adds photoionization heating. The net thermal loss is then , with possible heating-cooling equilibrium near . A single density-independent no longer describes every irradiated cloud: the answer also depends on radiation intensity, spectrum, shielding and density. The illustrative photoionized curve denotes altered cooling alone, not the net loss including heating.
Below , primordial gas can use molecular line cooling, especially rotational and vibrational transitions of molecular hydrogen and, in suitable chemical conditions, hydrogen deuteride. Molecule formation and protection against photodissociation are essential. Enriched gas can additionally use metal-line cooling, other molecular species and dust thermal emission, with gas-dust energy exchange important at high density. At high cosmological redshift, residual free electrons can transfer energy to the cosmic microwave background through Compton cooling by the cosmic microwave background if the gas is hotter than the radiation. Expansion can also cool gas adiabatically, but neither that process nor Compton cooling by the cosmic microwave background is an atomic loss law. The atomic threshold is a cooling-channel limitation, not a universal minimum gas temperature.
In the Lambda-CDM model, nearly Gaussian primordial density contrasts grow under gravity within an expanding universe containing cold dark matter, ordinary baryons and a cosmological constant. The cold dark matter is effectively collisionless and has negligible pressure on galactic scales. Before recombination, baryons are coupled to the photon fluid: radiation pressure and acoustic oscillations prevent their perturbations from behaving like pressureless matter. After recombination they can fall into the gravitational potentials already established by dark matter, subject to gas pressure and the Jeans mass.
For small density contrasts, evolution is linear. On pressure-free scales the growing mode is multiplied by the linear growth factor , with during matter domination. The cosmological density power spectrum can be written
Here is the cosmological transfer function. For nearly scale-invariant initial conditions, the large-scale matter spectrum behaves approximately as , whereas well inside the matter-radiation equality scale it falls approximately as , until the microscopic dark-matter cutoff matters. This fall of the dimensional does not imply less fluctuation power on every smaller mass scale: the power per logarithmic wavenumber is , and the smoothed matter density variance is obtained by integrating it against a mass-dependent window. Over the relevant cold-dark-matter hierarchy, smaller mass windows generally have larger variance.
The hierarchical galaxy formation picture follows: fluctuations on small mass scales typically reach the nonlinear collapse threshold first, while larger objects assemble later through accretion and dark-matter halo mergers. It is a statistical ordering, not a claim that every small object precedes every large rare peak. When becomes order unity, the linear growth factor is no longer a solution for the local density. Collisionless dark matter develops multistream regions and bound dark matter halos; phase mixing and violent relaxation redistribute orbital energies, and virialized structures approximately obey the virial theorem.
Baryons have an additional nonlinear route. Infall and shocks convert bulk kinetic energy into thermal energy, with characteristic virial temperature
Unlike collisionless dark matter, the gas can lose this energy through radiative cooling. The optically thin gas cooling time is the thermal-energy density divided by the radiative loss rate, for example
The density convention in must agree with the denominator; the astrophysical cooling function can also be defined using instead. The cooling criterion for galaxy formation compares this time with the collapse or supply time. Rapidly cooling gas loses pressure support, contracts and can form stars. Slowly cooling gas remains in a hot atmosphere. Stable virial shocks are not obligatory in every low-mass system: gas can also arrive in cold streams and cool while being accreted.
Angular momentum prevents indefinite radial contraction. Tidal torque theory supplies an initial halo spin, and later mergers change it. If gas radiates energy while retaining much of its specific angular momentum, it settles into a rotationally supported galactic disk rather than reaching the centre. The relation explains why modest halo spin can set a disk radius much smaller than its virial radius of a dark-matter halo. Torques, bars and gravitational encounters can transport angular momentum outwards and feed central concentrations; radiative cooling alone does not remove it.
Galaxy mergers alter stellar structure as well as assembling mass. A major galaxy merger can strongly disturb or destroy an existing galactic disk, randomizing stellar orbits and creating a spheroid through violent relaxation. A gas-rich galaxy merger also permits dissipation, inflow and a burst of star formation; gas left over or accreted afterwards can rebuild a galactic disk. Minor galaxy mergers add stars to outer components, thicken disks and grow bulges. Dry galaxy mergers add stellar mass and can increase size without much new star formation. Halo merging therefore does not imply instantaneous merging of its galaxies: satellite orbital decay requires Chandrasekhar dynamical friction and can take a substantial time.
The atomic and molecular cooling thresholds for galaxy formation supply a lower characteristic scale. Primordial atomic gas cools inefficiently below roughly because electronic excitation is suppressed. The corresponding halo mass is of order at a redshift of order ten, with approximate dependence at fixed threshold temperature. Molecular hydrogen can cool gas at hundreds of kelvin and permit smaller early objects, provided it forms and survives dissociating radiation. Metal-line cooling changes these thresholds after enrichment. Thus the atomic threshold is not an absolute minimum mass for all stellar systems.
At the other end, sufficiently massive dark matter halos have high virial temperatures and low-density hot gas. Above the strong atomic-line-cooling interval, thermal bremsstrahlung has . At comparable halo gas density, , while depends mainly on formation density. Cooling therefore becomes less able to condense all the gas within the available time. This upper galaxy mass from gas cooling argument selects galaxy-sized condensations, broadly halo masses around – in simple low-redshift estimates, rather than single luminous galaxies containing every baryon in a group or cluster. Its numerical scale depends on epoch, metallicity and gas profile. Subsequent galaxy mergers can build larger stellar systems; cooling is not an absolute upper bound on their final mass.
Finally, a Press-Schechter halo mass function has many low-mass objects and a steep high-mass cutoff, while the luminosity function of galaxies also has a faint component and a bright cutoff, often described by a Schechter function. The two shapes are related through the halo-to-galaxy luminosity mapping, not by identifying luminosity with total halo mass. In the idealized one-central-galaxy, no-scatter limit,
If and , then . A constant conversion gives similar shapes. Actual stellar feedback, inefficient low-temperature cooling and reionization suppress faint galaxies relative to small haloes, while long cooling times and active-galactic-nucleus feedback suppress luminosity at large halo masses. Satellites, scatter, stellar populations and dust further affect the correspondence. The halo hierarchy supplies the gravitational framework; cooling, angular momentum and feedback determine which parts become luminous galaxies.
For an optically thin, low-density plasma, the astrophysical cooling function packages collisional radiation losses and temperature-dependent ion fractions into the coefficient multiplying . Because the volume loss rate has units , the coefficient has units . The energy-density unit printed in the numerical hint cannot be the unit of this coefficient. Interpret the quoted logarithmic values in the dimensionally consistent cooling-coefficient unit.
Assume collisional ionization equilibrium, a primordial hydrogen-helium mixture and no external photoheating. In the specified range, the primordial atomic cooling curve has the following features. Just above , thermal electrons begin to excite neutral hydrogen efficiently; subsequent line emission, especially Lyman-alpha emission, causes a steep rise. The excitation rate contains a threshold factor of order . Hydrogen line cooling is strong near a few times ; the supplied value at provides a useful low-temperature label.
As hydrogen becomes ionized, neutral-hydrogen line cooling declines. Helium excitation and ionization produce a further shoulder or peak around . Collisional excitation, collisional ionization and radiative recombination all contribute: excitation photons remove electron kinetic energy, ionization consumes it, and recombination produces free-bound radiation. Once hydrogen and helium are almost fully stripped, their bound-state cooling disappears and the curve falls into a relatively inefficient interval. At high temperature, thermal bremsstrahlung dominates, with an approximate tail and a weak Gaunt-factor correction.
The sketch uses the two supplied numerical labels and a qualitative hydrogen-helium interpolation; it is not a tabulated atomic-rate calculation.
Figure 1.
Qualitative primordial atomic cooling curve with hydrogen and helium line features, a bremsstrahlung tail, and an illustrative metal-enriched comparison
. Details.
Metal-line cooling raises the cooling coefficient markedly over much of –, because heavier elements supply many ions and excitation transitions after hydrogen and helium have lost their bound electrons. It also broadens and reshapes the line-cooling peaks. Fine-structure lines can permit cooling below the hydrogen atomic threshold; molecular hydrogen can likewise cool metal-free gas below that threshold, but lies outside the requested temperature range. At sufficiently high temperatures thermal free-free emission again dominates the continuum. The enhanced curve in the figure is a schematic comparison, not a numerical claim about a specified metallicity. Thus metals generally shorten the cooling time and extend the temperature range of efficient cooling; ionization state, abundance and radiation field determine the actual curve.
Photoionization heating Created 2026-10-07 Updated 2026-10-14
Absorption of an ionizing photon leaves energy in a photoelectron above the binding energy. Thermalization supplies heat to the gas. Photoionization heating competes with the astrophysical cooling function, can establish thermal equilibrium, and can raise the pressure enough to inhibit baryon accretion into shallow dark-matter haloes.