Parker Solution: Transonic Winds & Magnetic Fields

• The Parker solution is the regular transonic wind model derived from isothermal Euler equations in spherical symmetry, balancing pressure, gravity, and expansion.
• It extends beyond stellar wind theory to include applications like the Parker spiral in heliophysics and stochastic formulations in cosmic-ray transport.
• Methodologies such as time-dependent characteristic analysis and stochastic SDE approaches ensure the selection of unique, dynamically stable solutions under realistic conditions.

Searching arXiv for recent and foundational uses of “Parker solution” and closely related Parker constructions. “Parker solution” most commonly denotes the transonic solution of the steady, spherically symmetric wind problem introduced for a thermally driven outflow from a gravitating point mass, but in technical usage the term also sits within a wider Parker-derived family that includes the Parker spiral of the heliospheric magnetic field, stochastic formulations of the Parker transport equation for galactic cosmic rays, and, in adjacent domains, Sweet–Parker reconnection and Parker weighting in short-scan tomography (Keto, 2020, Waters et al., 2012, Biondo et al., 2022, Wawrzynczak et al., 2015, Baty, 2022, Wiegmann et al., 25 Apr 2026). The unifying feature is not a single equation but a set of canonical balances—between pressure and gravity, radial expansion and rotation, advection and diffusion, or redundant and non-redundant measurements—each of which yields a distinguished regular solution.

1. Classical transonic Parker wind

In its standard form, the Parker solution is the steady-state solution of the isothermal Euler and continuity equations in spherical symmetry. With radial velocity $u(r)$, density $\rho(r)$, sound speed $a$, and point-mass potential $GM/r$, continuity gives a constant mass flux,

$\dot M = 4\pi r^2 \rho u,$

while the steady Euler equation reduces to a Bernoulli-type integral. In the nondimensional variables $x=r(a^2/GM)$ and $u\to u/a$, one convenient form is

$$\mathscr L(r,u)=\ln|u|-\tfrac12u^2+2\ln x+\frac{1}{x}=\text{const.}$$

The regular transonic solution is selected by the sonic-point condition

$u_c^2=a^2,\qquad r_c=\frac{GM}{2a^2},$

or equivalently $u_c=1$, $\rho(r)$0 in the normalization of Keto, and $\rho(r)$1 in the normalization of Waters and Proga (Keto, 2020, Waters et al., 2012).

A closely related differential form is

$\rho(r)$2

which is singular at $\rho(r)$3. Requiring a finite derivative there selects the physically regular branch. Waters and Proga further give an explicit isothermal solution in terms of the Lambert $\rho(r)$4 function,

$\rho(r)$5

with the relevant branch fixed by the sonic-point condition $\rho(r)$6 in their normalization (Waters et al., 2012).

The same mathematical structure covers both outward Parker wind and inward Bondi accretion. This suggests that “Parker solution” is best understood, in its narrowest sense, as the regular transonic branch of a Bernoulli problem with a critical point imposed by smoothness rather than by an external boundary prescription alone.

2. Time dependence, characteristics, and stability

The time-dependent formulation clarifies why the transonic Parker solution is not merely an exact steady branch but a dynamically selected one. In the isothermal case, eliminating the density yields the single PDE

$\rho(r)$7

which Keto casts as an initial-value problem and solves by the method of characteristics: $\rho(r)$8 The resulting characteristic flow reconstructs $\rho(r)$9 from Lagrangian labels and shows that a transonic flow develops under a wide range of realistic initial conditions (Keto, 2020).

Keto also reformulates the steady solution in Hamiltonian language by identifying the Bernoulli function with a Lagrangian per unit mass. The second variation is definite in sign for isothermal flows and also for adiabatic flows with $a$0, implying at least linear stability; for polytropic flows with $a$1, the second variation becomes indefinite (Keto, 2020). In the examples given, a uniform subsonic inflow evolves toward the transonic Bondi solution, and a uniform outward speed evolves toward the Parker wind solution, both passing smoothly through the sonic point.

A central selection principle emerges. Any initial profile whose characteristics cross into the “forbidden zone” between the subsonic and supersonic steady-state families is forced onto the unique transonic trajectory, whereas initial data lying entirely on a non-transonic steady branch merely oscillate about that branch (Keto, 2020). A plausible implication is that the Parker solution is not only a regular steady branch but also the attractor compatible with simultaneous inner and outer boundary constraints.

3. Rotation and disc-wind generalizations

Rotation qualitatively changes the topology of Parker-type solutions. In the rotating isothermal model of Lyutikov, working in the corotating frame with base radius $a$2, sound speed $a$3, and angular frequency $a$4, the azimuthal velocity is

$a$5

After eliminating the density through continuity, the radial Mach number $a$6 obeys

$a$7

with

$a$8

The regularity conditions produce two critical radii,

$a$9

where $GM/r$0 is the X-point through which the global transonic solution passes and $GM/r$1 is an O-point associated with a closed subsonic loop. These merge at

$GM/r$2

and for $GM/r$3 there is no real critical point and hence no smooth Parker-type transonic solution (Lyutikov, 2018).

For an equatorial disk wind launched from a Keplerian base, Lyutikov finds that the corona must be relatively cool in order to connect to the critical curve: $GM/r$4, reported in the abstract as “smaller than $GM/r$5” (Lyutikov, 2018). Waters and Proga obtain a related rotating critical-point equation,

$GM/r$6

which can admit two candidate transonic roots for the same physical conditions (Waters et al., 2012). They extend the Parker framework from spherical winds to disc winds by prescribing streamline geometries. In their notation, the effective area scales with a geometry parameter $GM/r$7, with $GM/r$8 for converging/Parker-like geometries and $GM/r$9 for the constant-inclination-angle (CIA) geometry. The lack of adjacent streamline divergence in the CIA model produces sonic points farther out, smaller mass-flux density, and more extended acceleration zones than in the converging model (Waters et al., 2012).

These rotating and geometric extensions preserve the Parker logic—regularity at a critical point still selects the admissible solution—but they replace the single classical sonic point by a richer phase portrait with multiple candidate branches, finite inner cutoffs, and geometry-dependent mass loading.

4. Parker spiral as the magnetic counterpart

The Parker solution also refers, in heliophysics, to the steady solar-wind magnetic geometry known as the Parker spiral. In the simplest analytic form, for a steady radial supersonic wind with constant $\dot M = 4\pi r^2 \rho u,$0,

$\dot M = 4\pi r^2 \rho u,$1

On the equator this becomes

$\dot M = 4\pi r^2 \rho u,$2

yielding the familiar Archimedean winding of field lines (Biondo et al., 2022).

Biondo et al. reconstruct this spiral with the Reverse In-situ and MHD APproach (RIMAP), described as a “hybrid analytical–numerical” procedure. PSP measurements taken between $\dot M = 4\pi r^2 \rho u,$3 and $\dot M = 4\pi r^2 \rho u,$4 AU are analytically back-mapped to an inner boundary at $\dot M = 4\pi r^2 \rho u,$5 using continuity, magnetic-flux conservation, Parker-spiral winding, and a temperature scaling $\dot M = 4\pi r^2 \rho u,$6. These back-mapped profiles are then imposed as time-dependent boundary conditions in a PLUTO MHD calculation on $\dot M = 4\pi r^2 \rho u,$7, $\dot M = 4\pi r^2 \rho u,$8, $\dot M = 4\pi r^2 \rho u,$9 with resolution $x=r(a^2/GM)$0 (Biondo et al., 2022).

The numerical solution recovers the analytic Parker scalings once steady state is reached: $x=r(a^2/GM)$1 and the spiral pitch angle matches the analytic formula. The reconstructed ecliptic-plane maps show fine longitudinal structure down to $x=r(a^2/GM)$2, heliospheric current-sheet crossings near $x=r(a^2/GM)$3, $x=r(a^2/GM)$4, and $x=r(a^2/GM)$5, and an Alfvén radius ranging from $x=r(a^2/GM)$6 to $x=r(a^2/GM)$7. Along the PSP track, agreement between the RIMAP solution and measurements is reported as typically better than $x=r(a^2/GM)$8–$x=r(a^2/GM)$9 for density and speed, while the large-scale HCS signatures are perfectly captured. Along the single field line intersecting the Metis plane of sky at $u\to u/a$0, the density agrees with Metis inversion within measurement errors and the velocity shows Parker-type acceleration, converging to $u\to u/a$1 km/s by $u\to u/a$2 (Biondo et al., 2022).

In this usage, the Parker solution is not a scalar wind profile but the magnetized large-scale geometry generated by that wind. The spiral is therefore the magnetic extension of the transonic outflow solution rather than an unrelated construction.

5. Parker transport equation and stochastic pseudoparticles

A further Parker-derived solution concept appears in cosmic-ray transport. The non-steady heliospheric transport of galactic cosmic ray particles is described by the Parker transport equation for the omni-directional distribution function $u\to u/a$3,

$u\to u/a$4

where $u\to u/a$5 is the solar-wind velocity, $u\to u/a$6 is the drift velocity, and $u\to u/a$7 is the symmetric part of the diffusion tensor in 3D (Wawrzynczak et al., 2015).

Bobik et al. rewrite the equation in backward-time Fokker–Planck form and derive an equivalent set of stochastic differential equations in heliocentric spherical coordinates. In Euler–Maruyama form,

$u\to u/a$8

$u\to u/a$9

$$\mathscr L(r,u)=\ln|u|-\tfrac12u^2+2\ln x+\frac{1}{x}=\text{const.}$$0

$$\mathscr L(r,u)=\ln|u|-\tfrac12u^2+2\ln x+\frac{1}{x}=\text{const.}$$1

with $$\mathscr L(r,u)=\ln|u|-\tfrac12u^2+2\ln x+\frac{1}{x}=\text{const.}$$2 independent Wiener increments and $$\mathscr L(r,u)=\ln|u|-\tfrac12u^2+2\ln x+\frac{1}{x}=\text{const.}$$3 reproducing the second-derivative structure of the transport operator (Wawrzynczak et al., 2015).

The paper stresses the preeminence of the backward over the forward approach. In the forward method, pseudoparticles are launched at the heliospheric boundary at $$\mathscr L(r,u)=\ln|u|-\tfrac12u^2+2\ln x+\frac{1}{x}=\text{const.}$$4 AU and only a small fraction reach the observer at $$\mathscr L(r,u)=\ln|u|-\tfrac12u^2+2\ln x+\frac{1}{x}=\text{const.}$$5 AU; in the backward method, all trajectories start at the observation point and are traced backward until they exit at $$\mathscr L(r,u)=\ln|u|-\tfrac12u^2+2\ln x+\frac{1}{x}=\text{const.}$$6 AU, so $$\mathscr L(r,u)=\ln|u|-\tfrac12u^2+2\ln x+\frac{1}{x}=\text{const.}$$7 of trajectories contribute to the local solution. The implementation uses spherical-coordinate boundary handling, a reflecting inner radial boundary $$\mathscr L(r,u)=\ln|u|-\tfrac12u^2+2\ln x+\frac{1}{x}=\text{const.}$$8 at $$\mathscr L(r,u)=\ln|u|-\tfrac12u^2+2\ln x+\frac{1}{x}=\text{const.}$$9 AU, an outer boundary at $u_c^2=a^2,\qquad r_c=\frac{GM}{2a^2},$0 AU, and the estimator

$u_c^2=a^2,\qquad r_c=\frac{GM}{2a^2},$1

A typical ensemble size is $u_c^2=a^2,\qquad r_c=\frac{GM}{2a^2},$2 trajectories for rigidity $u_c^2=a^2,\qquad r_c=\frac{GM}{2a^2},$3 GV (Wawrzynczak et al., 2015).

This stochastic Parker solution is used to model both Forbush decreases and the 27-day variation of GCR intensity. For the former, the parallel diffusion coefficient is taken as

$u_c^2=a^2,\qquad r_c=\frac{GM}{2a^2},$4

with $u_c^2=a^2,\qquad r_c=\frac{GM}{2a^2},$5, $u_c^2=a^2,\qquad r_c=\frac{GM}{2a^2},$6, $u_c^2=a^2,\qquad r_c=\frac{GM}{2a^2},$7, and a heliolongitudinal disturbance $u_c^2=a^2,\qquad r_c=\frac{GM}{2a^2},$8 for $u_c^2=a^2,\qquad r_c=\frac{GM}{2a^2},$9. For the 27-day variation, the solar-wind speed is modeled as

$u_c=1$0

with $u_c=1$1, $u_c=1$2, and $u_c=1$3 (Wawrzynczak et al., 2015). The stochastic-SDE and finite-difference solutions reproduce the magnitude and time profile of the observed Forbush decrease at rigidity thresholds $u_c=1$4 GV and $u_c=1$5 GV, and for the 27-day variation at $u_c=1$6 GV the two models agree within a few percent in amplitude and phase (Wawrzynczak et al., 2015).

6. Other Parker-derived constructions

The Parker name also labels mathematically distinct constructions outside the original wind problem.

Construction	Core relation	Reported quantitative behavior
Sweet–Parker reconnection	$u_c=1$7, $u_c=1$8	Classical SP recovered for small $u_c=1$9 and $\rho(r)$00; higher $\rho(r)$01 and/or $\rho(r)$02 produce setup-dependent deviations through $\rho(r)$03 (Baty, 2022)
Parker weighting	$\rho(r)$04 for short scans $\rho(r)$05	On the eXplore CT 120, Parker weighting corrects shading of $\rho(r)$06 HU near the phantom boundary and $\rho(r)$07 HU in the interior, with MTF$\rho(r)$08 changing from $\rho(r)$09 to $\rho(r)$10 lp/mm and $\rho(r)$11 curves overlapping within $\rho(r)$12 (Wiegmann et al., 25 Apr 2026)

In the visco-resistive generalization of Sweet–Parker reconnection, the steady-state outflow speed is

$\rho(r)$13

with $\rho(r)$14 and $\rho(r)$15. The usual Sweet–Parker relations are recovered in the limit of small $\rho(r)$16 and $\rho(r)$17, in particular the normalized reconnection rate $\rho(r)$18. In the opposite limit, numerical experiments in coalescence and tilt setups show that the upstream field $\rho(r)$19 is not fixed but varies with $\rho(r)$20 and $\rho(r)$21, producing modified scalings of $\rho(r)$22, $\rho(r)$23, and $\rho(r)$24 (Baty, 2022).

In cone-beam micro-CT, Parker weighting is a redundancy-correction formula for short-scan FDK reconstruction. For a scan of angular extent $\rho(r)$25, Parker’s piecewise $\rho(r)$26-ramp weighting rises smoothly from $\rho(r)$27 to $\rho(r)$28, stays at unity through the central non-redundant region, and ramps back to $\rho(r)$29, with the conjugate-ray constraint

$\rho(r)$30

On the eXplore CT 120 scanner, the reported implementation uses a $\rho(r)$31 short scan with $\rho(r)$32 projections and a precomputed $\rho(r)$33 weight map. The computational cost is negligible, reported as $\rho(r)$34 of the total. Quantitatively, Parker weighting removes deterministic shading without degrading image quality: MTF$\rho(r)$35 is $\rho(r)$36 lp/mm without weighting and $\rho(r)$37 lp/mm with weighting, the detectability index curves overlap within $\rho(r)$38, and the unweighted-minus-weighted residuals show correction of $\rho(r)$39 HU near the phantom boundary and $\rho(r)$40 HU in the homogeneous interior (Wiegmann et al., 25 Apr 2026).

These usages broaden the semantic range of “Parker solution” beyond stellar-wind theory. In strict astrophysical usage, however, the canonical Parker solution remains the regular transonic wind passing through a critical point. The Parker spiral, Parker transport equation, rotating Parker wind, and Parker-type disc winds are best regarded as mathematically and physically linked descendants of that original transonic construction, whereas Sweet–Parker reconnection and Parker weighting are separate Parker-associated formalisms that retain the same emphasis on selecting a unique regular balance under competing constraints.