Radiative Sweet-Parker Reconnection Model

• The paper demonstrates that incorporating optically thin radiative cooling into the classical Sweet-Parker model enables a variable current-sheet length, leading to accelerated reconnection.
• It employs a methodology where the outflow time across the sheet matches the local cooling time, resulting in a compressed, denser layer that is much shorter than the system scale.
• Experimental and numerical studies validate that strong radiative losses boost plasma compression and modify resistivity, significantly increasing reconnection rates beyond classical predictions.

Searching arXiv for the specified paper and closely related radiatively cooled reconnection work. arxiv_search(query="(Chowdhry et al., 18 Aug 2025) radiatively cooled Sweet-Parker current sheet formation under radiative cooling Uzdensky McKinney", max_results=10, sort_by="submittedDate") The radiatively-cooled Sweet–Parker model is a resistive-MHD extension of classical Sweet–Parker reconnection for plasmas in which optically thin radiative losses materially alter the current-sheet energy balance, compression, temperature, and reconnection rate. In its modern form, the framework combines the compression-based theory of Uzdensky and McKinney with a variable-length current-sheet closure motivated by X-point collapse under cooling: strong cooling can accelerate collapse in the inflow direction while arresting or reversing elongation in the outflow direction, so the sheet length need not remain fixed at the system size. The resulting steady state predicts that, when radiative cooling dominates compressional heating, the current sheet can become shorter than the system size and reconnect faster than in the classical Sweet–Parker limit (Chowdhry et al., 18 Aug 2025).

1. Classical basis and radiative generalization

The classical Sweet–Parker model describes steady-state reconnection in collisional plasmas as a long, thin current sheet with a normalized rate scaling as

$\frac{V_{\rm in}}{V_A} \sim S_L^{-1/2},$

where $S_L$ is the Lundquist number based on the sheet length. In that picture, the layer is geometrically constrained by the system scale and energy leaves primarily by advection rather than by photon emission.

The radiatively-cooled extension developed by Uzdensky and McKinney reformulates this picture for strong optically thin cooling in non-relativistic resistive MHD. In the zero-guide-field case, intense cooling lowers the layer temperature, increases the density through pressure balance, raises the Spitzer resistivity, and thereby enhances the reconnection rate relative to incompressible Sweet–Parker scaling. The compression ratio is written as

$A \equiv \frac{n}{n_0},$

and the modified scalings include

$$v_{\rm rec} \sim V_{A0} S_0^{-1/2} A^{1/2}, \qquad \delta \sim L S_0^{-1/2}A^{-1/2}.$$

Pressure balance gives

$$k_B T = \frac{B_0^2}{16\pi n} = \frac{k_B T_{\rm eq}}{A},$$

while the perpendicular Spitzer resistivity is

$$\eta_\perp = C_\eta\, c r_e\, \ln\Lambda\, \theta_e^{-3/2}.$$

This theory also identifies a condition for evolution to a strong-cooling, thermally stable layer when the cooling function satisfies $Q_{\rm rad}(n,T)\sim n^\alpha T^\beta$: the criterion is $\alpha < 1+\beta$. Bremsstrahlung, with $\alpha=2$ and $\beta=1/2$, fails this criterion; external inverse Compton and cyclotron cooling, with $S_L$0 and $S_L$1, satisfy it. In the strong-guide-field limit, compression is suppressed ($S_L$2), so cooling affects the reconnection rate primarily through the temperature dependence of resistivity rather than through density enhancement (Uzdensky et al., 2010).

2. Variable-length reformulation

The principal modification introduced in "Current sheet formation under radiative cooling" is to allow the current-sheet length $S_L$3 to be a dynamical variable rather than an imposed system-scale constant. The motivation is that strong radiative cooling can make the steady-state sheet much shorter than the system size: X-point collapse analysis shows that cooling accelerates collapse along the inflow direction but can arrest or even reverse sheet elongation in the outflow direction.

To incorporate that effect into a steady-state Sweet–Parker closure, the model enforces the condition that the advection or outflow time across the sheet be comparable to the local cooling time,

$S_L$4

This is the paper’s key new assumption. It expresses the requirement that a plasma parcel escape before catastrophic radiative collapse and links $S_L$5 directly to cooling physics rather than to external geometry. Relative to the original Uzdensky–McKinney construction, this change is both quantitative and qualitative: sheet length is no longer prescribed by $S_L$6, and the reconnection rate is affected not only by compression and resistivity but also by cooling-controlled geometric shortening (Chowdhry et al., 18 Aug 2025).

3. X-point collapse with optically thin cooling

The dynamical basis for the variable-length model is a self-similar, Lagrangian MHD treatment of X-point collapse with generic optically thin cooling,

$S_L$7

The similarity ansatz is

$S_L$8

with Lagrangian volume element

$S_L$9

The inflow-direction scale $A \equiv \frac{n}{n_0},$0 tracks sheet thinning, while the outflow-direction scale $A \equiv \frac{n}{n_0},$1 tracks elongation or contraction. The governing ODEs derived from MHD contain a cooling-dependent term involving

$A \equiv \frac{n}{n_0},$2

so the pressure evolution retains the time-integrated effect of radiative losses.

The asymptotic behavior separates weak and strong cooling regimes. For weak cooling, $A \equiv \frac{n}{n_0},$3, collapse proceeds as in the classic case but is accelerated; at late times,

$A \equiv \frac{n}{n_0},$4

with a correction term of order $A \equiv \frac{n}{n_0},$5. For strong cooling, $A \equiv \frac{n}{n_0},$6, the behavior becomes linear,

$A \equiv \frac{n}{n_0},$7

indicating faster collapse. In the same regime, $A \equiv \frac{n}{n_0},$8 can turn over, so that current-sheet elongation is arrested or reversed. For bremsstrahlung, the paper identifies a critical $A \equiv \frac{n}{n_0},$9 of order 10 separating the regime of continued elongation from the regime in which strong cooling stops or reverses it. This dynamical result is the direct precursor to the steady-state assumption that $$v_{\rm rec} \sim V_{A0} S_0^{-1/2} A^{1/2}, \qquad \delta \sim L S_0^{-1/2}A^{-1/2}.$$0 should be set by the competition between escape and cooling rather than by the system scale alone (Chowdhry et al., 18 Aug 2025).

4. Steady-state scalings, geometry, and reconnection rate

In the steady radiatively-cooled state, the central compressibility is characterized by

$$v_{\rm rec} \sim V_{A0} S_0^{-1/2} A^{1/2}, \qquad \delta \sim L S_0^{-1/2}A^{-1/2}.$$1

When radiative cooling dominates compressional heating, the model predicts a current sheet with $$v_{\rm rec} \sim V_{A0} S_0^{-1/2} A^{1/2}, \qquad \delta \sim L S_0^{-1/2}A^{-1/2}.$$2, $$v_{\rm rec} \sim V_{A0} S_0^{-1/2} A^{1/2}, \qquad \delta \sim L S_0^{-1/2}A^{-1/2}.$$3, and a reconnection rate larger than the classical Sweet–Parker value. For bremsstrahlung, $$v_{\rm rec} \sim V_{A0} S_0^{-1/2} A^{1/2}, \qquad \delta \sim L S_0^{-1/2}A^{-1/2}.$$4 and $$v_{\rm rec} \sim V_{A0} S_0^{-1/2} A^{1/2}, \qquad \delta \sim L S_0^{-1/2}A^{-1/2}.$$5, the paper gives the explicit scalings

$$v_{\rm rec} \sim V_{A0} S_0^{-1/2} A^{1/2}, \qquad \delta \sim L S_0^{-1/2}A^{-1/2}.$$6

$$v_{\rm rec} \sim V_{A0} S_0^{-1/2} A^{1/2}, \qquad \delta \sim L S_0^{-1/2}A^{-1/2}.$$7

and

$$v_{\rm rec} \sim V_{A0} S_0^{-1/2} A^{1/2}, \qquad \delta \sim L S_0^{-1/2}A^{-1/2}.$$8

where $$v_{\rm rec} \sim V_{A0} S_0^{-1/2} A^{1/2}, \qquad \delta \sim L S_0^{-1/2}A^{-1/2}.$$9 encodes the relative strength of cooling and $$k_B T = \frac{B_0^2}{16\pi n} = \frac{k_B T_{\rm eq}}{A},$$0 is the classical Sweet–Parker thickness. The physical interpretation given in the paper is that strong radiative losses remove internal thermal pressure faster than compressional heating can restore it, producing a thinner, denser, and shorter layer. Because the Lundquist number based on the actual sheet length is reduced, the normalized rate $$k_B T = \frac{B_0^2}{16\pi n} = \frac{k_B T_{\rm eq}}{A},$$1 increases correspondingly (Chowdhry et al., 18 Aug 2025).

A compact comparison of the three constructions emphasized in the literature is as follows.

Model	Sheet length	Reconnection-rate scaling
Classical Sweet–Parker	$$k_B T = \frac{B_0^2}{16\pi n} = \frac{k_B T_{\rm eq}}{A},$$2	$$k_B T = \frac{B_0^2}{16\pi n} = \frac{k_B T_{\rm eq}}{A},$$3
Uzdensky–McKinney radiative SP	$$k_B T = \frac{B_0^2}{16\pi n} = \frac{k_B T_{\rm eq}}{A},$$4	$$k_B T = \frac{B_0^2}{16\pi n} = \frac{k_B T_{\rm eq}}{A},$$5
Variable-length strong-cooling model	$$k_B T = \frac{B_0^2}{16\pi n} = \frac{k_B T_{\rm eq}}{A},$$6	Increased above classical SP

A common oversimplification is to treat radiative cooling as a correction to the sheet thermodynamics alone. In the variable-length formulation, cooling also controls the large-scale geometry. Another oversimplification is to assume that any strong cooling law supports a stationary fast layer. The earlier stability result for $$k_B T = \frac{B_0^2}{16\pi n} = \frac{k_B T_{\rm eq}}{A},$$7 indicates that not every cooling mechanism admits a steady strong-cooling solution; this is especially relevant for bremsstrahlung-dominated cases, for which the older theory emphasized the possibility of cooling catastrophe rather than smooth steady evolution (Uzdensky et al., 2010).

5. Experimental and numerical realizations

Laboratory experiments and resistive-MHD simulations have provided concrete realizations of the radiatively-cooled regime. In the first experimental study of strongly radiatively-cooled magnetic reconnection, two exploding aluminum wire arrays driven simultaneously on the Z machine with $$k_B T = \frac{B_0^2}{16\pi n} = \frac{k_B T_{\rm eq}}{A},$$8 and $$k_B T = \frac{B_0^2}{16\pi n} = \frac{k_B T_{\rm eq}}{A},$$9 generated a radiatively-cooled reconnection layer with $$\eta_\perp = C_\eta\, c r_e\, \ln\Lambda\, \theta_e^{-3/2}.$$0 and $$\eta_\perp = C_\eta\, c r_e\, \ln\Lambda\, \theta_e^{-3/2}.$$1. X-ray measurements showed a narrow $$\eta_\perp = C_\eta\, c r_e\, \ln\Lambda\, \theta_e^{-3/2}.$$2 FWHM burst at $$\eta_\perp = C_\eta\, c r_e\, \ln\Lambda\, \theta_e^{-3/2}.$$3 after current start, consistent with rapid formation and cooling of the layer. Time-gated images revealed hotspots moving at up to $$\eta_\perp = C_\eta\, c r_e\, \ln\Lambda\, \theta_e^{-3/2}.$$4, interpreted as plasmoids; spectroscopy showed hotspot temperatures of $$\eta_\perp = C_\eta\, c r_e\, \ln\Lambda\, \theta_e^{-3/2}.$$5, much larger than the inflow and bulk-layer temperatures, and these structures generated the majority of the Al K-shell emission at around $$\eta_\perp = C_\eta\, c r_e\, \ln\Lambda\, \theta_e^{-3/2}.$$6 prior to the onset of cooling. The inferred layer thickness, $$\eta_\perp = C_\eta\, c r_e\, \ln\Lambda\, \theta_e^{-3/2}.$$7, matched the Sweet–Parker prediction, while the energy budget was decisively cooling-dominated: $$\eta_\perp = C_\eta\, c r_e\, \ln\Lambda\, \theta_e^{-3/2}.$$8 in the layer and up to $$\eta_\perp = C_\eta\, c r_e\, \ln\Lambda\, \theta_e^{-3/2}.$$9 in hotspots, compared with $Q_{\rm rad}(n,T)\sim n^\alpha T^\beta$0 and $Q_{\rm rad}(n,T)\sim n^\alpha T^\beta$1 (Datta et al., 2024).

Complementary 2D and 3D simulations of the MARZ platform were carried out in GORGON, with radiative losses modeled using non-local thermodynamic equilibrium tables from Spk and, separately, $Q_{\rm rad}(n,T)\sim n^\alpha T^\beta$2 multi-group radiation transport. These simulations reproduced highly collisional, super-Alfvénic ($Q_{\rm rad}(n,T)\sim n^\alpha T^\beta$3) and supersonic ($Q_{\rm rad}(n,T)\sim n^\alpha T^\beta$4–$Q_{\rm rad}(n,T)\sim n^\alpha T^\beta$5) inflows that formed a reconnection layer with $Q_{\rm rad}(n,T)\sim n^\alpha T^\beta$6 and $Q_{\rm rad}(n,T)\sim n^\alpha T^\beta$7. When $Q_{\rm rad}(n,T)\sim n^\alpha T^\beta$8, the layer underwent radiative collapse into a cold, strongly compressed current sheet with accelerated reconnection. At $Q_{\rm rad}(n,T)\sim n^\alpha T^\beta$9, the non-radiative case had width $\alpha < 1+\beta$0, aspect ratio $\alpha < 1+\beta$1, temperature $\alpha < 1+\beta$2, $\alpha < 1+\beta$3, $\alpha < 1+\beta$4, and normalized reconnection rate $\alpha < 1+\beta$5; the locally cooled case had width $\alpha < 1+\beta$6, aspect ratio $\alpha < 1+\beta$7, temperature $\alpha < 1+\beta$8, $\alpha < 1+\beta$9, $\alpha=2$0, and rate $\alpha=2$1; the radiation-transport case had width $\alpha=2$2, aspect ratio $\alpha=2$3, temperature $\alpha=2$4, $\alpha=2$5, $\alpha=2$6, and rate $\alpha=2$7. Plasmoids formed in both regimes, but under strong cooling they were quenched or extinguished before ejection (Datta et al., 2024).

Solar-chromospheric reconnection provides a useful contrast rather than a direct realization of the same mechanism. In a partially ionized low-solar-chromosphere model with a more realistic OPACITY/CHIANTI-based cooling function, low-$\alpha=2$8 cases remained close to single-fluid Sweet–Parker dynamics, whereas a higher-$\alpha=2$9 case with $\beta=1/2$0 developed ion–neutral inflow decoupling and a reconnection rate about three times faster than the Sweet–Parker prediction. The same study found that plasma temperature still increased inside the current sheet and that stronger radiative losses primarily regulated temperature and ionization rather than generically producing the compression-driven fast-reconnection regime emphasized in optically thin HED plasmas (Ni et al., 2018).

6. Assumptions, limitations, and unresolved issues

The radiatively-cooled Sweet–Parker model is built on a restricted but analytically tractable set of assumptions. In the variable-length formulation, the plasma is optically thin and locally cooled with $\beta=1/2$1; the dynamics are treated in non-relativistic MHD; resistivity becomes important only near current-sheet formation or maximum compression; radiation pressure and Compton drag are neglected; and the transition to steady state is assumed to occur when ohmic heating balances radiative cooling. The construction also neglects plasmoid and tearing instabilities and assumes a single-sheet geometry (Chowdhry et al., 18 Aug 2025).

These assumptions delimit the range of applicability. Earlier theory already noted that the existence of a stable strong-cooling layer depends on the structure of the cooling function through $\beta=1/2$2; in that sense, rapid cooling is not by itself sufficient to guarantee a steady radiatively-compressed Sweet–Parker state (Uzdensky et al., 2010). The experimental and numerical work shows that plasmoid formation can coexist with cooling-dominated reconnection, and in simulations strong cooling can even extinguish islands before ejection, which indicates that single-sheet steady-state closures omit dynamically important structure when the layer becomes unstable (Datta et al., 2024). In partially ionized chromospheric plasmas, realistic cooling does not necessarily accelerate reconnection above Sweet–Parker in low-$\beta=1/2$3 regimes, because ion–neutral coupling and thermochemical effects can dominate the response (Ni et al., 2018).

Several unresolved questions follow directly from the existing literature. One concerns the relation between the earlier bremsstrahlung stability caveat and the newer variable-length construction: the latter shows that strong cooling can arrest outflow elongation and motivate a shorter steady-state sheet, while the former emphasized that bremsstrahlung may lead to cooling catastrophe rather than an evolutionary approach to a stable strong-cooling state. This suggests that geometrical shortening and thermal stability are related but distinct issues. Another concerns missing high-energy-density physics: radiative transfer, pair production, strong radiation fields, and kinetic or out-of-equilibrium effects may invalidate the optically thin MHD closure in some astrophysical environments. The broader significance of the model is therefore not that it closes the problem of radiative reconnection, but that it provides a minimal analytic framework in which cooling alters not only dissipation and compression, but also the global geometry of the reconnecting current sheet (Chowdhry et al., 18 Aug 2025).