Lecture 13: Interstellar Shocks: The Rankine-Hugoniot Jump Conditions

ASTR 360 The Interstellar Medium

Context and Why This Matters

Supernova blast waves (the ultimate fate of the massive stars whose formation this course traced in Unit 5, and whose youth as ionizing sources was the subject of Unit 6) are the dominant source of interstellar shocks, and are believed to heat and maintain the hot ionized medium (Lecture 01) and to accelerate the cosmic rays this unit’s final lecture examines.

Learning Goals

  1. Derive the Rankine-Hugoniot jump conditions (conservation of mass, momentum, and energy across a shock front) for an ideal gas.
  2. Derive the strong-shock (high Mach number) limits for density compression and post-shock temperature.
  3. Apply the strong-shock post-shock temperature formula to a real supernova-remnant shock and compare to observed X-ray temperatures.

Opening Phenomenon

The Cassiopeia A supernova remnant, the visible debris of a stellar explosion whose light reached Earth around 1680, has an outer shock front still expanding at roughly 5,000 km/s -- more than ten times the speed of any sound wave that could propagate ahead of it through the surrounding interstellar medium, making this a genuinely supersonic, strong shock whose passage instantaneously and irreversibly heats, compresses, and accelerates the gas it overtakes.

Terms Developed in Context

shock front

A thin transition where supersonic flow is abruptly compressed, decelerated, and heated. In this lecture, the term is not an isolated label: This concept is needed to reason about interstellar shocks: the rankine-hugoniot jump conditions. Conservation of mass, momentum, and energy flux across a steady, planar shock front (in the shock’s own rest frame) gives three coupled algebraic equations (the Rankine-Hugoniot jump conditions) relating the pre-shock (subscript 1) and post-shock (subscript 2) density, pressure, and velocity, parametrized by the upstream Mach number M = v_1/c_{s,1} (the ratio of the pre-shock flow speed to the pre-shock sound speed).

It is tied to the lecture’s observables or evidence in this way: A shock front is a thin transition layer across which gas properties (density, pressure, temperature, velocity) change discontinuously (on a microscopic, collisional-mean-free-path length scale) because the flow ahead of the shock is moving faster than the local sound speed can communicate the disturbance -- exactly the situation for supernova ejecta plowing into the surrounding ISM at thousands of km/s. A second connection is important for using the term in physical reasoning: In the strong-shock limit (M>>1, an excellent approximation for supernova-remnant shocks moving at thousands of km/s into a ~10 km/s-sound-speed warm ISM), these equations simplify dramatically: for a monatomic ideal gas (adiabatic index gamma=5/3), the density compression ratio approaches a fixed value rho2/rho1 -> (gamma+1)/(gamma-1) = 4, independent of the shock’s actual speed, while the post-shock temperature continues to grow with the square of the shock velocity, T2 = 3 mu m_H v_shock^2/(16k). The observational interpretation is constrained by Chandra X-ray Observatory imaging and spectroscopy of Cassiopeia A and Tycho’s supernova remnant reveal exactly the kind of multi-million-kelvin, shock-heated plasma this lecture’s post-shock temperature formula predicts, with observed X-ray temperatures in the few-keV range broadly consistent with strong collisionless-shock heating (Vink 2004; Williams et al. 2016).

Boundary and common confusion. Assuming the predicted post-shock temperature exactly equals the observed X-ray plasma temperature. Real collisionless astrophysical shocks (unlike idealized collisional gas-dynamic shocks) do not necessarily equilibrate electron and ion temperatures instantaneously, so the electron temperature actually measured spectroscopically can be lower than this formula’s ion-temperature-like prediction if electron-ion equilibration is incomplete -- a genuine, actively studied physical subtlety, honestly disclosed here rather than glossed over as an exact match. The term should therefore be applied only after identifying the measured quantity, the model assumption that links it to the desired physical quantity, and the regime in which that link remains valid.

Rankine-Hugoniot conditions

Conservation relations for mass, momentum, and energy across a steady shock discontinuity. In this lecture, the term is not an isolated label: This concept is needed to reason about interstellar shocks: the rankine-hugoniot jump conditions. In the strong-shock limit (M>>1, an excellent approximation for supernova-remnant shocks moving at thousands of km/s into a ~10 km/s-sound-speed warm ISM), these equations simplify dramatically: for a monatomic ideal gas (adiabatic index gamma=5/3), the density compression ratio approaches a fixed value rho2/rho1 -> (gamma+1)/(gamma-1) = 4, independent of the shock’s actual speed, while the post-shock temperature continues to grow with the square of the shock velocity, T2 = 3 mu m_H v_shock^2/(16k).

It is tied to the lecture’s observables or evidence in this way: Chandra X-ray Observatory imaging and spectroscopy of Cassiopeia A and Tycho’s supernova remnant reveal exactly the kind of multi-million-kelvin, shock-heated plasma this lecture’s post-shock temperature formula predicts, with observed X-ray temperatures in the few-keV range broadly consistent with strong collisionless-shock heating (Vink 2004; Williams et al. 2016). A second connection is important for using the term in physical reasoning: This strong-shock post-shock temperature formula directly converts an observationally accessible quantity (shock velocity, measured from proper motions or Doppler-shifted line emission) into a predicted X-ray-emitting plasma temperature -- a direct, testable link between a remnant’s dynamics and its observed X-ray spectrum. The observational interpretation is constrained by The strong compression (density increase) predicted at a shock front is directly connected to why supernova remnants and other shocked regions appear as thin, bright, compressed shells rather than smoothly fading transitions -- exactly the morphology seen in Cas A and Tycho’s remnant’s Chandra/optical imagery.

Boundary and common confusion. Assuming the predicted post-shock temperature exactly equals the observed X-ray plasma temperature. Real collisionless astrophysical shocks (unlike idealized collisional gas-dynamic shocks) do not necessarily equilibrate electron and ion temperatures instantaneously, so the electron temperature actually measured spectroscopically can be lower than this formula’s ion-temperature-like prediction if electron-ion equilibration is incomplete -- a genuine, actively studied physical subtlety, honestly disclosed here rather than glossed over as an exact match. The term should therefore be applied only after identifying the measured quantity, the model assumption that links it to the desired physical quantity, and the regime in which that link remains valid.

Mach number

The flow speed divided by the relevant sound speed, quantifying how strongly supersonic a shock is. In this lecture, the term is not an isolated label: This concept is needed to reason about interstellar shocks: the rankine-hugoniot jump conditions. This strong-shock post-shock temperature formula directly converts an observationally accessible quantity (shock velocity, measured from proper motions or Doppler-shifted line emission) into a predicted X-ray-emitting plasma temperature -- a direct, testable link between a remnant’s dynamics and its observed X-ray spectrum.

It is tied to the lecture’s observables or evidence in this way: The strong compression (density increase) predicted at a shock front is directly connected to why supernova remnants and other shocked regions appear as thin, bright, compressed shells rather than smoothly fading transitions -- exactly the morphology seen in Cas A and Tycho’s remnant’s Chandra/optical imagery. A second connection is important for using the term in physical reasoning: Conservation of mass, momentum, and energy flux across a steady, planar shock front (in the shock’s own rest frame) gives three coupled algebraic equations (the Rankine-Hugoniot jump conditions) relating the pre-shock (subscript 1) and post-shock (subscript 2) density, pressure, and velocity, parametrized by the upstream Mach number M = v_1/c_{s,1} (the ratio of the pre-shock flow speed to the pre-shock sound speed). The observational interpretation is constrained by A shock front is a thin transition layer across which gas properties (density, pressure, temperature, velocity) change discontinuously (on a microscopic, collisional-mean-free-path length scale) because the flow ahead of the shock is moving faster than the local sound speed can communicate the disturbance -- exactly the situation for supernova ejecta plowing into the surrounding ISM at thousands of km/s.

Boundary and common confusion. Assuming the predicted post-shock temperature exactly equals the observed X-ray plasma temperature. Real collisionless astrophysical shocks (unlike idealized collisional gas-dynamic shocks) do not necessarily equilibrate electron and ion temperatures instantaneously, so the electron temperature actually measured spectroscopically can be lower than this formula’s ion-temperature-like prediction if electron-ion equilibration is incomplete -- a genuine, actively studied physical subtlety, honestly disclosed here rather than glossed over as an exact match. The term should therefore be applied only after identifying the measured quantity, the model assumption that links it to the desired physical quantity, and the regime in which that link remains valid.

strong shock limit

The high-Mach-number asymptotic regime in which upstream thermal pressure is negligible in the jump conditions. In this lecture, the term is not an isolated label: This concept is needed to reason about interstellar shocks: the rankine-hugoniot jump conditions. Conservation of mass, momentum, and energy flux across a steady, planar shock front (in the shock’s own rest frame) gives three coupled algebraic equations (the Rankine-Hugoniot jump conditions) relating the pre-shock (subscript 1) and post-shock (subscript 2) density, pressure, and velocity, parametrized by the upstream Mach number M = v_1/c_{s,1} (the ratio of the pre-shock flow speed to the pre-shock sound speed).

It is tied to the lecture’s observables or evidence in this way: A shock front is a thin transition layer across which gas properties (density, pressure, temperature, velocity) change discontinuously (on a microscopic, collisional-mean-free-path length scale) because the flow ahead of the shock is moving faster than the local sound speed can communicate the disturbance -- exactly the situation for supernova ejecta plowing into the surrounding ISM at thousands of km/s. A second connection is important for using the term in physical reasoning: In the strong-shock limit (M>>1, an excellent approximation for supernova-remnant shocks moving at thousands of km/s into a ~10 km/s-sound-speed warm ISM), these equations simplify dramatically: for a monatomic ideal gas (adiabatic index gamma=5/3), the density compression ratio approaches a fixed value rho2/rho1 -> (gamma+1)/(gamma-1) = 4, independent of the shock’s actual speed, while the post-shock temperature continues to grow with the square of the shock velocity, T2 = 3 mu m_H v_shock^2/(16k). The observational interpretation is constrained by Chandra X-ray Observatory imaging and spectroscopy of Cassiopeia A and Tycho’s supernova remnant reveal exactly the kind of multi-million-kelvin, shock-heated plasma this lecture’s post-shock temperature formula predicts, with observed X-ray temperatures in the few-keV range broadly consistent with strong collisionless-shock heating (Vink 2004; Williams et al. 2016).

Boundary and common confusion. Assuming the predicted post-shock temperature exactly equals the observed X-ray plasma temperature. Real collisionless astrophysical shocks (unlike idealized collisional gas-dynamic shocks) do not necessarily equilibrate electron and ion temperatures instantaneously, so the electron temperature actually measured spectroscopically can be lower than this formula’s ion-temperature-like prediction if electron-ion equilibration is incomplete -- a genuine, actively studied physical subtlety, honestly disclosed here rather than glossed over as an exact match. The term should therefore be applied only after identifying the measured quantity, the model assumption that links it to the desired physical quantity, and the regime in which that link remains valid.

compression ratio

The downstream-to-upstream mass-density ratio across a shock. In this lecture, the term is not an isolated label: This concept is needed to reason about interstellar shocks: the rankine-hugoniot jump conditions. In the strong-shock limit (M>>1, an excellent approximation for supernova-remnant shocks moving at thousands of km/s into a ~10 km/s-sound-speed warm ISM), these equations simplify dramatically: for a monatomic ideal gas (adiabatic index gamma=5/3), the density compression ratio approaches a fixed value rho2/rho1 -> (gamma+1)/(gamma-1) = 4, independent of the shock’s actual speed, while the post-shock temperature continues to grow with the square of the shock velocity, T2 = 3 mu m_H v_shock^2/(16k).

It is tied to the lecture’s observables or evidence in this way: Chandra X-ray Observatory imaging and spectroscopy of Cassiopeia A and Tycho’s supernova remnant reveal exactly the kind of multi-million-kelvin, shock-heated plasma this lecture’s post-shock temperature formula predicts, with observed X-ray temperatures in the few-keV range broadly consistent with strong collisionless-shock heating (Vink 2004; Williams et al. 2016). A second connection is important for using the term in physical reasoning: This strong-shock post-shock temperature formula directly converts an observationally accessible quantity (shock velocity, measured from proper motions or Doppler-shifted line emission) into a predicted X-ray-emitting plasma temperature -- a direct, testable link between a remnant’s dynamics and its observed X-ray spectrum. The observational interpretation is constrained by The strong compression (density increase) predicted at a shock front is directly connected to why supernova remnants and other shocked regions appear as thin, bright, compressed shells rather than smoothly fading transitions -- exactly the morphology seen in Cas A and Tycho’s remnant’s Chandra/optical imagery.

Boundary and common confusion. Assuming the predicted post-shock temperature exactly equals the observed X-ray plasma temperature. Real collisionless astrophysical shocks (unlike idealized collisional gas-dynamic shocks) do not necessarily equilibrate electron and ion temperatures instantaneously, so the electron temperature actually measured spectroscopically can be lower than this formula’s ion-temperature-like prediction if electron-ion equilibration is incomplete -- a genuine, actively studied physical subtlety, honestly disclosed here rather than glossed over as an exact match. The term should therefore be applied only after identifying the measured quantity, the model assumption that links it to the desired physical quantity, and the regime in which that link remains valid.

post-shock temperature

The thermal temperature immediately downstream of a shock after bulk kinetic energy has been randomized, subject to species equilibration. In this lecture, the term is not an isolated label: This concept is needed to reason about interstellar shocks: the rankine-hugoniot jump conditions. This strong-shock post-shock temperature formula directly converts an observationally accessible quantity (shock velocity, measured from proper motions or Doppler-shifted line emission) into a predicted X-ray-emitting plasma temperature -- a direct, testable link between a remnant’s dynamics and its observed X-ray spectrum.

It is tied to the lecture’s observables or evidence in this way: The strong compression (density increase) predicted at a shock front is directly connected to why supernova remnants and other shocked regions appear as thin, bright, compressed shells rather than smoothly fading transitions -- exactly the morphology seen in Cas A and Tycho’s remnant’s Chandra/optical imagery. A second connection is important for using the term in physical reasoning: Conservation of mass, momentum, and energy flux across a steady, planar shock front (in the shock’s own rest frame) gives three coupled algebraic equations (the Rankine-Hugoniot jump conditions) relating the pre-shock (subscript 1) and post-shock (subscript 2) density, pressure, and velocity, parametrized by the upstream Mach number M = v_1/c_{s,1} (the ratio of the pre-shock flow speed to the pre-shock sound speed). The observational interpretation is constrained by A shock front is a thin transition layer across which gas properties (density, pressure, temperature, velocity) change discontinuously (on a microscopic, collisional-mean-free-path length scale) because the flow ahead of the shock is moving faster than the local sound speed can communicate the disturbance -- exactly the situation for supernova ejecta plowing into the surrounding ISM at thousands of km/s.

Boundary and common confusion. Assuming the predicted post-shock temperature exactly equals the observed X-ray plasma temperature. Real collisionless astrophysical shocks (unlike idealized collisional gas-dynamic shocks) do not necessarily equilibrate electron and ion temperatures instantaneously, so the electron temperature actually measured spectroscopically can be lower than this formula’s ion-temperature-like prediction if electron-ion equilibration is incomplete -- a genuine, actively studied physical subtlety, honestly disclosed here rather than glossed over as an exact match. The term should therefore be applied only after identifying the measured quantity, the model assumption that links it to the desired physical quantity, and the regime in which that link remains valid.

Evidence

Derivation and Model

Working Equation

\[ \frac{\rho_2}{\rho_1} \xrightarrow{M\gg1} \frac{\gamma+1}{\gamma-1} = 4\ (\gamma{=}5/3), \qquad T_2 = \frac{3\,\mu\, m_H\, v_{\rm shock}^2}{16\,k} \]

Worked Example

  1. Cassiopeia A’s forward shock (v_shock=5,000 km/s, standard published value, Vink 2004 and subsequent Chandra proper-motion studies) gives a predicted strong-shock post-shock (ion) temperature of T2=3.46e+08 K, equivalent to kT2=29.85 keV -- roughly an order of magnitude above the few-keV electron temperature Chandra actually measures from Cas A’s shocked shell, an honest, expected discrepancy (not a computational error) that is itself real evidence for incomplete electron-ion equilibration at this fast, young shock (Ghavamian, Laming & Rakowski 2007).
  2. Tycho’s supernova remnant (v_shock=4,600 km/s, Williams et al. 2016) gives a predicted post-shock ion temperature of T2=2.93e+08 K (kT2=25.27 keV), somewhat lower than Cas A’s, consistent with (and directly following from) its somewhat lower shock velocity through the T2 proportional to v_shock^2 scaling -- and, like Cas A, well above Tycho’s observed few-keV X-ray electron temperature for the same electron-ion equilibration reason.
  3. The strong-shock compression ratio, rho2/rho1=4.0, applies to both remnants’ forward shocks (both are firmly in the strong-shock, high-Mach-number regime given their thousands-of-km/s speeds against a ~10 km/s warm-ISM sound speed), illustrating that the compression ratio saturates at a fixed value even though the two remnants’ predicted temperatures differ substantially.

Common Misconception

Assuming the predicted post-shock temperature exactly equals the observed X-ray plasma temperature. Real collisionless astrophysical shocks (unlike idealized collisional gas-dynamic shocks) do not necessarily equilibrate electron and ion temperatures instantaneously, so the electron temperature actually measured spectroscopically can be lower than this formula’s ion-temperature-like prediction if electron-ion equilibration is incomplete -- a genuine, actively studied physical subtlety, honestly disclosed here rather than glossed over as an exact match.

Active-Learning Guidance

Using T2 proportional to v_shock^2, explain why a supernova remnant’s X-ray temperature is expected to decline as the remnant ages and its forward shock decelerates (as it sweeps up and is slowed by an increasing mass of surrounding interstellar gas), and connect this to why very young remnants like Cas A and Tycho are hotter, harder X-ray sources than older, more evolved remnants.

Lab Connection

Lab 07’s capstone applies this lecture’s shock-jump formulas to Cassiopeia A alongside its cosmic-ray and magnetic-field content (Lecture 14), synthesizing this course’s final unit into one real, multiwavelength object.

Synthesis Questions

References