Cooling flow Created 2026-10-20 Updated 2026-10-20
An inward flow driven by loss of thermal pressure as hot gas undergoes radiative cooling. A classical unheated cluster model predicts substantial central cooling and deposition of cold material. active-galactic-nucleus feedback, transport and other heating can reduce the actual inflow; an X-ray cooling luminosity is not automatically a measured mass-deposition rate.
In an optically thin plasma, one convention defines the radiative cooling power per volume by , where counts hydrogen nuclei. A convention instead has . Thus the tabulated functions differ even for the same gas. With total particle density , a monatomic-gas radiative gas cooling time is when heating is negligible.[ref]
Gas-rich galaxy merger Created 2026-10-06 Updated 2026-10-25
A gas-rich galaxy merger involves a substantial reservoir of gas capable of radiative cooling. Dissipation and gravitational torques permit gas inflow and enhanced star formation, unlike the nearly collisionless stellar evolution in a dry galaxy merger. A surviving or replenished gas reservoir can subsequently form a galactic disk.
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.
When collisions thermalize an optically thin two-level atom, the upper fraction is and the emitted power per atom is . Its logarithmic temperature slope is . Radiative cooling still occurs, but collision frequency no longer sets the emitted power as in the low-density limit.
Lava lake Created 2026-10-06 Updated 2026-10-14
A lava lake is a persistent exposed body of molten or partly molten magma. Surface radiative cooling, thermal convection, latent heat, and internal radiogenic heating can jointly control formation and melting of its crust.
Work in the rest frame of the shock. The normal momentum flux is conserved, so . In the strong-shock limit , while mass conservation gives . The given Strong-shock Rankine-Hugoniot conditions therefore yield
Using the downstream ideal gas law gives
For a monatomic perfect gas, and the coefficient is . This is the immediate post-shock temperature, before substantial radiative cooling; ionization, dissociation or unequal electron and ion temperatures require additional physical modeling.
In a spherical self-gravitating equilibrium with positive mass density, away from the center. Thus can locally be replaced by , so both and can be expressed as functions of . Hydrostatic equilibrium further gives
Since decreases outward wherever , eliminating gives a barotropic equation of state along this equilibrium. This is a relation describing its structure; it does not require all thermodynamic perturbations of an arbitrary material to be barotropic.
For , integrate the polytropic enthalpy balance :
For a finite-radius zero-pressure boundary with and , and therefore . With the gravitational potential zero at infinity, an isolated sphere has . It is an integration constant tied to that boundary, not necessarily a zero of the gravitational potential.
Put and , with . The preceding relation gives . Substitution into the spherical gravitational Poisson equation yields the Lane-Emden equation
Let be the first zero of the regular Lane-Emden function and . Then the Lane-Emden mass formula and radius are
For a common fixed polytropic index with a finite first zero, their scalings are
The usual regular zero-pressure models have a finite first zero for ; the model has infinite radius, so a finite-radius formula cannot simply be applied to it.
For an ideal gas of fixed composition, . Consequently homologous polytropes at fixed central temperature have and both , recovering . The common central temperature, composition and dimensionless structure are the relevant assumptions. Keeping fixed instead gives the different polytropic mass-radius relation, and no finite index makes its generic fixed- exponent exactly one.
For adiabatic monatomic hydrogen, and . At a common fixed specific entropy, hence common , the same relations give
Here , and the omitted constant in the last relation depends on . A population with different cloud entropies need not follow that mass sequence.
Efficient radiative cooling in molecular clouds prevents compression from forcing the temperature upward as steeply as in an adiabatic monatomic gas. Rotational transitions of molecules and radiation by dust provide low-temperature cooling channels; which dominates depends on composition, opacity and temperature. It is therefore possible to assemble dense massive cold clouds while disposing of compressional heat. This explains the direction of the idealized comparison, rather than proving a universal cloud-mass cutoff.
For a common isothermal sound speed , a spherical isothermal equilibrium has equal to times a dimensionless structural factor. Holding that factor fixed gives the homologous isothermal mass-radius relation
For example, the singular isothermal sphere has . A regular finite cloud requires a pressure-confined isothermal sphere or another outer boundary: an unconfined isothermal sphere does not end at a finite vacuum radius. Equal temperature alone, while allowing arbitrary truncation or density contrast, does not make the coefficient universal.
Once swept-up matter dominates the ejecta and radiative cooling is negligible, the Sedov-Taylor blast wave has its shock radius fixed by dimensional analysis:
Here is a dimensionless constant of order unity fixed by the adiabatic index. Treating the ambient matter as hydrogen gives , while . Thus
For a monatomic gas, the usual similarity normalization instead gives about . An estimate of a few parsecs is therefore appropriate. Negligible ejecta mass does not by itself guarantee this energy-conserving stage: substantial radiation losses or a strongly nonuniform ambient mass density would invalidate the assumed similarity law.
Immediately upstream, free fall from rest at infinity gives . Treat the surface shock as locally planar and stationary, with negligible upstream pressure. The Strong-shock Rankine-Hugoniot conditions for a monatomic gas, , give and . momentum conservation then gives
Use the source's gas-constant convention . The immediate postshock temperature is
Here adiabatic means that the jump does not lose appreciable heat to the surroundings; the shock is nevertheless irreversible and raises entropy. The gas still has velocity immediately downstream, so setting all upstream kinetic energy equal to downstream internal energy would give the wrong coefficient. Later settling, radiative cooling or changing ionization can alter the temperature; they are outside this immediate, fixed-composition perfect-gas jump.
A massive main-sequence donor loses material in a stellar wind. Gravitational capture by the neutron star can produce a nearly radial, supersonic infall when the captured material has little specific angular momentum. A magnetic field may channel the final infall toward polar accretion columns, where a standing shock decelerates the flow. Large angular momentum instead favours an accretion disk; strictly radial accretion is an idealization.
For mean particle mass , the ideal gas law is . Combining it with b gives
For Newtonian free fall from a distance large compared with the stellar radius , . Therefore
TODO: a unique numerical temperature and characteristic photon band require source data not supplied here. The held question gives but gives neither nor a composition fixing . The mass alone does not fix the free-fall speed or the temperature. For example, doubling at the same mass and composition halves both quantities in the last box.
A normalization, rather than an additional claimed datum, makes the scale explicit:
Here is the characteristic particle energy; the mean translational energy of a nonrelativistic thermal particle is . If one additionally chose and , this ideal one-temperature estimate would correspond to gamma-ray-scale thermal photon energies. Those are conditional choices, not information recovered from the source. Real accreting neutron stars commonly emit X-rays: electron-ion disequilibrium, radiative cooling, radiation pressure, reprocessing and relativistic infall invalidate a direct identification of this ideal ion-heating estimate with the escaping photon spectrum.
Methane reaching the stratosphere is oxidized and is an important source of stratospheric water vapor. That water can feed HOx radical family chemistry and alter radiative cooling and polar-cloud conditions. Methane also converts chlorine atoms into hydrogen chloride, opposing some chlorine-driven ozone loss. These effects compete: increasing methane does not have one universal ozone sign at every altitude and latitude.
For the albedo-to-forcing conversion, Earth intercepts sunlight over a disk but distributes the global mean over its full sphere, giving incident mean flux . The absorbed shortwave flux is , and hence
The cooling magnitude is for the assumed one-year effect of one mass unit of injection. Albedo does not change the sunlight incident at the top of the atmosphere; it changes the absorbed/net incoming flux. The minus sign therefore accompanies the radiative cooling.
galactic feedback couples star formation and black-hole accretion back to the gas supply. Stellar radiation, winds, supernovae and cosmic rays inject energy and momentum; active-galactic-nucleus feedback adds radiation, jets and winds powered by accretion. Heating, turbulence or gas ejection can suppress later star formation, while compression can sometimes promote it. Negative feedback permits self-regulation rather than unrestricted gas consumption.
Define the mass-loading factor . If the launch speed scales with the escape speed and the injected energy per stellar mass is fixed, conservation of energy suggests ; fixed injected momentum instead suggests . Coupling efficiencies, radiative losses and recycling change these ideal scalings. Thermal energy injected into dense gas can radiate away before driving an outflow, so an energy budget alone does not demonstrate effective feedback.
Radio-filled X-ray cavities provide a mechanical diagnostic. For a cavity containing a relativistic fluid with , its enthalpy is ; dividing by an estimated age gives a mechanical-power estimate. Chandra X-ray Observatory observations of the Perseus cluster revealed fronts and ripples associated with radio bubbles, interpreted as weak shocks and pressure waves that can transport jet energy into the cooling atmosphere.[ref]
The test requires ambient density and temperature estimates, radio/X-ray spatial comparison, projected cavity volumes and an age model. Brightness edges must be distinguished from absorption and cool rims; wave energy is not automatically thermalized locally. These controls determine whether inferred mechanical heating can actually offset radiative cooling.
Galactic star formation is limited by how much gas becomes cold, shielded and self-gravitating, and by how long it remains in that state. A useful local estimate is
The star-formation efficiency per free-fall time absorbs the fact that turbulence, magnetic support and stellar feedback prevent all eligible gas from collapsing in one free-fall time of a uniform sphere. Rotation and disk stability influence the supply of dense clouds, while gravitational collapse and gravitational fragmentation operate inside them.
For galaxy-averaged measurements, the Kennicutt–Schmidt law has , with the classic mixed spiral/starburst sample giving . This empirical relation was obtained using H-alpha and far-infrared star-formation tracers together with atomic-hydrogen and CO gas measurements.[ref] Dust attenuation, contamination by older stellar heating, the stellar initial mass function and the CO-to-molecular-gas conversion must be controlled. A disk-averaged total-gas exponent is not automatically a small-scale cloud law.
Young stars then photoionize gas, drive winds and eventually inject supernova energy and momentum. These processes can disrupt clouds, sustain random motions or expel gas; radiative cooling and fresh inflow rebuild the cold supply. Large-scale star formation is a regulated conversion of gas, rather than unimpeded free fall of the entire disk. A long gas-depletion time is compatible with short galactic dynamical times in the dense fraction.
First specify the perturbation: take fixed gas density, fixed ionic abundance and the same temperature-independent heating, so thermal stability requires at an equilibrium. Let , be the collider density and the atomic density. A two-level optically thin balance is
The critical density of an atomic transition is approximately . At low density, each rare collisional excitation is followed by an escaping photon, so
The isochoric thermal stability of a collisionally excited line therefore gives
At the linear restoring coefficient vanishes; the cooling maximum is marginal to first order and is not robust two-sided stability.
At high density, collisional de-excitation competes successfully with radiative decay, and collisions establish local thermodynamic equilibrium level populations. Detailed balance gives , where . The cooling is then
Thus the high-density fixed-density branch is stable for finite , becoming weakly restoring as the upper fraction saturates at high temperature. It is incorrect to retain the low-density factor in the radiative cooling law once populations are thermalized.[ref]
The isobaric thermal stability of a two-level coolant is a distinct question. If the intended heating is a fixed rate per particle, then and stability tests the cooling per particle. The low-density boundary becomes ; the high-density cooling per particle remains proportional to , so remains stable. If instead heating is constant per volume, the isobaric low-density boundary is , and the high-density stability condition is . The source gives no pressure constraint or explicit heating normalization: the boxed result is the fixed-density answer, with these conditional alternatives stated to avoid asserting a universal boundary.
Figure 1.
Low-density collisional line cooling has a temperature maximum, whereas high-density LTE emission rises toward saturation
. Source. Details.
A protostar forms when a self-gravitating cloud fragment can contract faster than pressure support and energy loss can arrest it. A useful local scale is the Jeans length, , with characteristic Jeans mass and free-fall time
for an initially uniform sphere. Geometry, turbulence and magnetic support change the detailed threshold; none of these estimates asserts that every cloud contracts as a uniform sphere.
While cooling is efficient, compression is nearly isothermal. Then decreases as the density rises, allowing fragmentation. Rotation becomes increasingly important during contraction because a parcel conserving specific angular momentum encounters a centrifugal scale . The resulting accretion disc must transport angular momentum before matter can reach the central object. Infall, a dense accreting interior and an outer envelope can therefore coexist; a protostar is not initially an isolated fixed-mass star.
Increasing optical depth eventually slows radiative cooling. The central gas then becomes less compressible and can form a pressure-supported first core. Molecular hydrogen dissociation absorbs compressional energy and can trigger a second collapse before a compact stellar core is established. Subsequent accretion supplies a luminosity of order
The retained fraction depends on shock and disc radiation, so this expression is an energy scale rather than a guaranteed luminosity all emerging from the stellar surface.
A central accreting core can be nearly hydrostatic on its short sound-crossing time while its mass, entropy and radius evolve on longer times. A low-mass object is often largely convective because of its opacity and steep radiative temperature gradient; deuterium burning can temporarily contribute to support and moderate contraction. Neither convection nor the timing of deuterium ignition is universal across all masses and accretion histories.
When the main accretion phase ends, gravitational contraction remains an energy source. The virial theorem explains how contraction both heats the interior and supplies radiation: for an ideal-gas-supported object with negligible surface-pressure correction, , so only part of the released binding energy appears as increased thermal energy. The corresponding Kelvin-Helmholtz time is
Fully convective pre-main-sequence evolution approximately follows a Hayashi track, with decreasing radius and luminosity at comparatively weakly changing effective temperature. Objects developing radiative interiors can follow Henyey tracks. Eventually sustained core hydrogen burning balances the luminosity and the star reaches the main sequence. Massive objects can ignite hydrogen while substantial accretion continues, so “accretion ends, then hydrogen ignites” is not a universal ordering.
A useful structural distinction is therefore between an accretion-powered protostar and a largely fixed-mass pre-main-sequence contraction model. An infrared-bright envelope or disc traces obscured formation and reprocessing of radiation; a schematic evolutionary track should not replace the actual mass- and entropy-dependent stellar-structure problem.
A spatially uniform pressure is a useful leading approximation when the sound speed makes the acoustic crossing time much shorter than the radiative cooling and thermal conduction times. The motions must remain at low Mach number, and gravity and other stresses must be negligible. Sound then redistributes pressure before appreciable thermal evolution occurs. Small pressure corrections can still drive the slow motion; imposing exactly zero pressure gradient in the full momentum equation at every order would be a stronger assumption.
The equation of state gives the specific internal energy of an ideal gas as . With , the continuity equation is . Differentiating therefore yields
Insert this into the internal energy equation, whose compressional term is , and use spatially uniform pressure, so . The result is
For temporally constant pressure put , so the mass density is . The continuity equation now gives . Multiplication of the last equation by , with , gives
This is the isobaric conductive cooling in mass coordinates transformation before changing coordinates.
There is a boundary assumption implicit in the proposed Lagrangian coordinate. For a mass integral from a fixed origin, integration of the continuity equation gives
Thus is material when the base mass flux vanishes, as appropriate for the stationary cold reservoir below. If the origin permits mass flux, replace by ; this corrected label has zero material derivative. In a material mass coordinate, and . In particular
The conduction coefficient therefore becomes , with the same as the cooling coefficient. With , the equation is
In the hot gas there is no volumetric radiative cooling. The initial temperature is uniform and the cold boundary introduces no finite mass length. The remaining nonlinear diffusion equation is invariant under , with the same temperature scale. Its natural similarity solution is therefore
The far-field condition also recovers the initial temperature at every fixed as . Substitution gives the ordinary differential equation
For , this reduces to . Integrating first for gives . The two boundary conditions fix , hence the error-function isobaric cooling front is
Here the error function has normalization , as printed in the original PDF; the converted TeX's is incorrect. The following original sketch shows the temperature profiles. The cooling layer broadens as , and at fixed positive mass label the temperature falls as time increases, while the remote hot gas remains at .
Figure 1.
Isobaric conductive cooling profiles at increasing scaled times
. Details.
For , the magnitude of the heat flux into the cold reservoir is
Although the thermal conductivity tends to zero there, the physical temperature gradient diverges: their product has this finite limit. The stationary reservoir has no incoming advective energy flux and radiates the conducted energy immediately. Since , the radiated power per unit area is
Thus the requested exponent is . The initial divergence is the usual idealization of an instantaneous temperature discontinuity in a diffusion equation.
If the gas can cool appreciably below the virial temperature, radiative cooling removes thermal energy and pressure support. In a dark matter halo the gas then contracts, dissipating more energy as it falls. Efficient condensation requires the radiative gas cooling time to be short enough compared with the relevant dynamical or assembly time; merely having an available low-temperature transition does not guarantee that the gas reaches it quickly.
The collisionless dark matter cannot lose comparable energy through radiation and remains extended. Gas with appreciable conserved angular momentum stops radial collapse when rotation supports it, often forming a disk; lower-angular-momentum gas reaches a more compact central region. Cold dense gas can fragment into self-gravitating clouds and form stars if its gravitational instability overcomes remaining support. Stellar feedback subsequently reheats or expels gas and regulates the conversion. Efficient cooling enables central baryonic condensation and star formation; angular momentum and feedback determine the resulting galaxy.
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.