Parker Solar Wind Model

Parker solar wind model is the foundational hydrodynamic and magnetohydrodynamic theory formulated in 1958 by American astrophysicist Eugene Newman Parker demonstrating that a million-kelvin stellar corona with high thermal conductivity cannot maintain static hydrostatic equilibrium out to large distances, and must instead expand continuously from subsonic velocities near the coronal base through a transonic transition point—the Parker sonic radius (r_s)—to form a supersonic solar wind streaming into interplanetary space. Published in The Astrophysical Journal after initial resistance from reviewers who favored Sydney Chapman's static conductive corona or Ludwig Biermann's corpuscular radiation notes, Parker's 1958 paper simultaneously predicted that solar rotation winds the frozen-in radial coronal magnetic field into an Archimedean spiral—the Parker spiral—across the heliosphere. Within astronomy and space physics, the Parker solar wind model was confirmed in situ between 1959 and 1962 by the Soviet Luna 1 and NASA Mariner 2 spacecraft, governs rotational magnetic braking (Skumanich spin-down) of main-sequence stars, drives space weathering and atmospheric stripping across unmagnetized planets such as Mars and Venus, and is directly sampled inside the Alfvén surface by NASA's Parker Solar Probe (launched in 2018).

In 1957, British geophysicist Sydney Chapman modeled the solar corona—heated above the photosphere (whose baseline temperature T_eff = 5,772 K is set by the Eddington-Barbier relation) to T₀ ≈ 1.5 × 10⁶ K by magnetohydrodynamic wave dissipation and nanoflares—as a static, spherically symmetric atmosphere in hydrostatic equilibrium (v = 0). Because thermal conductivity in a fully ionized hydrogen plasma scales steeply with temperature via the Spitzer law (κ_th = κ₀ T^5/2), heat conduction in a static corona with no heat sinks at infinity requires the radial heat flux −4π r² κ₀ T^5/2 (dT / dr) to remain constant, yielding Chapman's slow radial temperature decay:

T(r) = T₀ (r₀ / r)^2/7

Substituting this conductive temperature profile into the static equation of hydrostatic equilibrium, dP / dr = − (G M_⊙ ρ) / r² with the ideal gas law P = (ρ / [μ m_p]) k_B T (where μ ≈ 0.60 is the mean molecular weight of fully ionized solar plasma), and integrating from the coronal base (r = r₀, pressure P₀) to radial infinity (r → ∞) gives the asymptotic pressure at infinity P_∞:

P(r → ∞) = P₀ exp[ − (7 / 5) (G M_⊙ μ m_p) / (r₀ k_B T₀) ]

Evaluating this expression for the solar corona (T₀ = 1.5 × 10⁶ K, r₀ ≈ 1.05 R_⊙, P₀ ≈ 3 × 10⁻³ Pa) yields an asymptotic static coronal pressure of P_∞ ≈ 10⁻⁵ to 10⁻⁶ Pa (n k_B T ~ 10⁻⁴ dyn cm⁻²). However, the total pressure of the local Milky Way interstellar medium (thermal gas plus magnetic field and cosmic rays) is only P_ISM ≈ 10⁻¹³ Pa (~10⁻¹² dyn cm⁻²)—seven to eight orders of magnitude smaller than P_∞! Eugene Parker recognized that because the interstellar medium cannot provide the confining boundary pressure P_∞ required to hold a static coronal atmosphere in place, the solar corona cannot remain static (v = 0) and must instead expand outward dynamically as a hydrodynamic wind.

Conservation of Mass and Momentum in an Isothermal Corona

To capture the core physics induced by high Spitzer thermal conductivity, Parker (1958) analyzed a steady-state (∂/∂t = 0), spherically symmetric radial outflow with velocity v(r) and mass density ρ(r) in an isothermal corona of constant temperature T₀ and isothermal sound speed c_s:

c_s = ( k_B T₀ / [μ m_p] )^1/2 ≈ 128 km s⁻¹ × ( T₀ / [1.2 × 10⁶ K] )^1/2

with ideal gas pressure P(r) = ρ(r) c_s². The conservation of mass (continuity equation) across spherical shells of area 4π r² fixes the constant solar mass-loss rate Ṁ:

4π r² ρ(r) v(r) = Ṁ = constant ⇒ (1 / ρ)(dρ / dr) = − (1 / v)(dv / dr) − (2 / r)

The radial Euler momentum conservation equation balancing convective acceleration, outward thermal pressure gradient, and inward solar gravity is:

ρ v (dv / dr) = − (dP / dr) − (G M_⊙ ρ / r²) = − c_s² (dρ / dr) − (G M_⊙ ρ / r²)

Dividing through by ρ and substituting (1 / ρ)(dρ / dr) from the continuity equation eliminates density ρ(r) and produces Parker's solar wind differential equation:

( v − [c_s² / v] ) (dv / dr) = (2 c_s² / r) − (G M_⊙ / r²)

Or equivalently, written in terms of the radial Mach number M_a(r) = v(r) / c_s:

( 1 / v ) (dv / dr) ( M_a² − 1 ) = (2 / r) [ 1 − (r_s / r) ]

The Parker Sonic Radius (r_s) and Transonic Topology

The right-hand side of Parker's momentum equation vanishes at a unique radial distance where the thermal expansion term 2c_s² / r balances solar gravity G M_⊙ / r². This radius is the Parker sonic radius (r_s):

r_s = G M_⊙ / (2 c_s²) = (G M_⊙ μ m_p) / (2 k_B T₀) ≈ 5.75 R_⊙ × ( 1.0 × 10⁶ K / T₀ )

For a coronal temperature of T₀ = 1.5 × 10⁶ K (c_s ≈ 144 km s⁻¹), the sonic point lies at r_s ≈ 3.8 R_⊙ above the center of the Sun:

• Inside the Sonic Radius (r < r_s): The right-hand side is negative (gravity exceeds spherical geometry divergence). For the wind to accelerate outward (dv / dr > 0) from near-zero velocity at the coronal base (v(r₀) ≪ c_s), the left-hand bracket (M_a² − 1) must also be negative (v < c_s, subsonic flow).
• At the Sonic Radius (r = r_s): For the velocity gradient dv / dr to remain finite and smooth as the right-hand side passes through zero at r = r_s, the flow speed must simultaneously equal the sound speed: v(r_s) = c_s (M_a = 1).
• Beyond the Sonic Radius (r > r_s): Both (r − r_s) and (M_a² − 1) become positive, so the wind continues to accelerate (dv / dr > 0) into a supersonic solar wind (v > c_s).

Integrating Parker's equation analytically yields the transcendental velocity topology:

( v / c_s )² − ln[ ( v / c_s )² ] = 4 ln( r / r_s ) + 4 ( r_s / r ) − 3

Among the four mathematical branches passing through the (r, v) phase plane—(I) double-valued unphysical loops, (II) purely subsonic "solar breeze" solutions that reproduce Chapman's finite pressure P_∞ > 0 at infinity, (III) supersonic-at-base accretion/infall solutions (the inverse of spherical Bondi accretion), and (IV) the unique transonic Parker solar wind solution—only Branch IV starts subsonically at the photosphere/corona (v ≪ c_s), passes smoothly through v(r_s) = c_s, and achieves vanishing thermal pressure P(r → ∞) → 0 at large heliocentric distances (where v(r) ≈ 2 c_s [ln(r / r_s)]^1/2 ≈ 400–800 km s⁻¹ at 1 astronomical unit).

Because the ionized solar wind plasma has high electrical conductivity, Alfvén's theorem of flux freezing locks the coronal magnetic field lines into the radially expanding wind while their footpoints remain anchored in the rotating solar photosphere (sidereal angular velocity Ω_⊙ ≈ 2.865 × 10⁻⁶ rad s⁻¹, period ≈ 25.38 days). In the equatorial plane (θ = π/2), a plasma parcel launched from r₀ with radial speed v_r traces an Archimedean spiral—the Parker spiral—with magnetic field components:

B_r(r) = B₀ (r₀ / r)², B_φ(r) = − B₀ (r₀ / r)² [ Ω_⊙ (r − r₀) / v_r ]

The local spiral pitch angle ψ between the interplanetary magnetic field vector and the radial direction satisfies:

tan ψ = |B_φ / B_r| = Ω_⊙ (r − r₀) / v_r

At Earth's orbit (r = 1 au = 1.496 × 10⁸ km) for v_r ≈ 430 km s⁻¹, the Parker spiral angle is ψ ≈ 45° ("garden-hose angle"), tightening to ψ ≈ 84° at the orbit of Neptune (30 au). Furthermore, because magnetic tension forces the outflowing wind to co-rotate with the star out to the Alfvén radius r_A ≈ 15–20 R_⊙ (where wind speed equals the Alfvén wave speed v_A = B / √[μ₀ρ]), the solar wind carries away specific angular momentum j = Ω_⊙ r_A² (dJ / dt = −[2/3] Ṁ Ω_⊙ r_A²), explaining why solar-type stars spin down with age as Ω(t) ∝ t^−1/2 (the Skumanich law, 1972) even though the Sun's present mass-loss rate is only Ṁ_⊙ ≈ 2 × 10⁻¹⁴ M_⊙ yr⁻¹ (compared to radiation-driven winds in Luminous Blue Variables and Wolf-Rayet stars). In April 2021, NASA's Parker Solar Probe flew beneath the Alfvén critical surface at r = 18.8 R_⊙ (and subsequently down to 9.86 R_⊙), observing magnetic "switchbacks" and Alfvénic turbulence that heat and accelerate the nascent solar wind.

1. Parker, E. N. (1958). "Dynamics of the Interplanetary Gas and Magnetic Fields." The Astrophysical Journal, 128, 664–676. https://doi.org/10.1086/146579
2. Neugebauer, M., & Snyder, C. W. (1966). "Mariner 2 Observations of the Solar Wind: 1. Average Properties." Journal of Geophysical Research, 71(19), 4469–4484. https://doi.org/10.1029/JZ071i019p04469
3. Weber, E. J., & Davis, L., Jr. (1967). "The Angular Momentum of the Solar Wind." The Astrophysical Journal, 148, 217–227. https://doi.org/10.1086/149138
4. Kasper, J. C., Klein, K. G., Lichko, E., Huang, J., Chen, C. H. K., Badman, S. T., ... & Zank, G. P. (2021). "Parker Solar Probe Enters the Magnetically Dominated Solar Corona." Physical Review Letters, 127(25), 255101. https://doi.org/10.1103/PhysRevLett.127.255101