Part III: Solar Magnetic Activity | Chapter 12

Solar Flares

CSHKP model, magnetic reconnection, Sweet-Parker and Petschek rates, and particle acceleration

12.1 The Standard Flare Model (CSHKP)

Derivation 1: Flare Energy from Magnetic Field Annihilation

The CSHKP (Carmichael-Sturrock-Hirayama-Kopp-Pneuman) model describes the standard two-ribbon flare geometry: a vertical current sheet forms above a magnetic arcade, reconnection occurs, and energy is released into the corona and chromosphere.

Step 1. The available magnetic free energy is the excess over the potential field:

$$E_{\text{free}} = \frac{1}{2\mu_0}\int (B^2 - B_p^2) \, dV$$

Step 2. For a volume \(V = L^3\) with field \(B \sim 0.03\) T (300 G):

$$\boxed{E_{\text{free}} \sim \frac{B^2 L^3}{2\mu_0} \sim \frac{(0.03)^2 \times (3 \times 10^7)^3}{2 \times 4\pi \times 10^{-7}} \sim 10^{25} \text{ J (X-class flare)}}$$

This is \(\sim 10^{32}\) erg, consistent with observations of large flares. The flare classification is based on peak soft X-ray flux: C (\(10^{-6}\) W/m\(^2\)), M (\(10^{-5}\)), X (\(10^{-4}\)).

12.2 Sweet-Parker Reconnection

Derivation 2: Sweet-Parker Reconnection Rate

Sweet (1958) and Parker (1957) derived the reconnection rate for a steady-state resistive current sheet of length \(L\) and thickness \(\delta\).

Step 1. Mass conservation: plasma flows in at speed \(v_{\text{in}}\)through area \(L\) and out at \(v_{\text{out}}\) through area \(\delta\):

$$\rho v_{\text{in}} L = \rho v_{\text{out}} \delta$$

Step 2. The outflow is driven by the magnetic tension, accelerating plasma to the Alfven speed: \(v_{\text{out}} = v_A = B/\sqrt{\mu_0\rho}\).

Step 3. The current sheet thickness is set by resistive diffusion balancing the inflow:\(\delta = \eta / v_{\text{in}}\).

Step 4. Combining these:

$$\boxed{\frac{v_{\text{in}}}{v_A} = \frac{1}{\sqrt{R_m}} = \frac{1}{\sqrt{Lv_A/\eta}}}$$$$\boxed{\delta = \frac{L}{\sqrt{R_m}}}$$

For the solar corona: \(R_m \sim 10^{12}\), so \(v_{\text{in}}/v_A \sim 10^{-6}\). This gives a reconnection time \(\tau \sim L/(v_{\text{in}}) \sim 10^6 L/v_A \sim 10^8\) s — far too slow to explain flares (which release energy in minutes). This is the "reconnection rate problem."

12.3 Petschek Reconnection

Derivation 3: Fast Reconnection via Standing Shocks

Petschek (1964) proposed that the current sheet need not extend over the full system length. Instead, slow-mode standing shocks form, processing the inflowing plasma over a much larger area.

Step 1. The Petschek geometry has a compact diffusion region of length \(l \ll L\)with two pairs of standing slow-mode shocks extending from the ends.

Step 2. The maximum reconnection rate:

$$\boxed{\frac{v_{\text{in}}}{v_A} \sim \frac{\pi}{8\ln R_m} \sim 0.01\text{--}0.1}$$

This is fast enough to explain flare timescales. Modern numerical simulations with anomalous resistivity, plasmoid instability, and Hall MHD confirm reconnection rates of 0.01-0.1 of the Alfven speed, consistent with Petschek-like fast reconnection.

12.4 Particle Acceleration

Derivation 4: DC Electric Field Acceleration

Step 1. In the reconnection region, a DC electric field \(E \sim v_{\text{in}}B\)accelerates charged particles along the current sheet:

$$E \sim v_{\text{in}} B \sim 0.01 v_A B = 0.01 \times 10^6 \times 0.01 = 100 \text{ V/m}$$

Step 2. Over a length \(\ell\), a particle gains energy:

$$\boxed{\varepsilon = qE\ell \sim 100 \times 10^7 \text{ eV} = 1 \text{ GeV (for } \ell = 10^7 \text{ m)}}$$

Flares accelerate electrons to 10-100 keV (producing hard X-rays via bremsstrahlung) and protons to MeV-GeV energies. The nonthermal emission carries 10-50% of the total flare energy.

12.5 Thermal and Nonthermal Emission

Derivation 5: The Neupert Effect

The Neupert effect (1968) is the empirical observation that the soft X-ray (SXR) light curve is approximately the time integral of the hard X-ray (HXR) light curve.

Step 1. Accelerated electrons stream down magnetic field lines and produce HXR bremsstrahlung when they hit the dense chromosphere (thick-target model).

Step 2. The deposited energy heats chromospheric plasma, which expands upward (chromospheric evaporation), filling coronal loops with hot (\(\sim 10^7\) K) plasma that emits in soft X-rays.

Step 3. The SXR emission measure grows as the integral of the energy deposition rate:

$$\boxed{F_{\text{SXR}}(t) \propto \int_0^t F_{\text{HXR}}(t') \, dt'}$$

This is a key diagnostic linking the acceleration (impulsive phase, HXR) to the thermal response (gradual phase, SXR), supporting the thick-target model of flare energy deposition.

Numerical Simulation

Solar Flares: Reconnection Rates, Neupert Effect, Current Sheet, Energy Spectrum

Python
script.py150 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

12.6 Full Sweet-Parker Reconnection Derivation

From Mass and Energy Conservation

Step 1: Geometry. Consider a steady-state current sheet of length \(L\) and half-thickness \(\delta\). Plasma flows in at speed \(v_{\text{in}}\) from top and bottom, and exits at speed \(v_{\text{out}}\) from the two ends.

Step 2: Mass conservation. Incompressible flow in the 2D sheet:

$$\rho\,v_{\text{in}}\,(2L) = \rho\,v_{\text{out}}\,(2\delta) \quad\Rightarrow\quad v_{\text{in}} L = v_{\text{out}}\delta \quad\ldots(1)$$

Step 3: Energy conservation. The outflow is accelerated by the magnetic tension (the \(j\times B\) force). The upstream magnetic energy density is converted to kinetic energy:

$$\frac{B^2}{2\mu_0} = \frac{1}{2}\rho\,v_{\text{out}}^2 \quad\Rightarrow\quad v_{\text{out}} = v_A = \frac{B}{\sqrt{\mu_0\rho}} \quad\ldots(2)$$

Step 4: Ohmic diffusion balance. Inside the sheet, the inflow speed must balance the rate of magnetic diffusion across the sheet thickness:

$$v_{\text{in}} = \frac{\eta}{\delta} \quad\Rightarrow\quad \delta = \frac{\eta}{v_{\text{in}}} \quad\ldots(3)$$

Step 5: Combining. From (1): \(\delta = v_{\text{in}}L/v_A\). Equating with (3):

$$\frac{\eta}{v_{\text{in}}} = \frac{v_{\text{in}}L}{v_A} \quad\Rightarrow\quad v_{\text{in}}^2 = \frac{\eta\,v_A}{L}$$

Step 6. Defining \(R_m = Lv_A/\eta\):

$$\boxed{\frac{v_{\text{in}}}{v_A} = \frac{1}{\sqrt{R_m}} \approx 10^{-6} \text{ for coronal } R_m\sim 10^{12}}$$$$\boxed{\delta = \frac{L}{\sqrt{R_m}} \approx 1 \text{ m for } L=10^7 \text{ m}}$$

The reconnection timescale \(\tau_{SP} = L/v_{\text{in}} = \sqrt{\tau_A\tau_d}\) is the geometric mean of the Alfven time and diffusion time. For the corona: \(\tau_{SP} \sim 10^8\) s (~3 years), far too slow for flares observed to release energy in minutes.

12.7 Why Petschek Reconnection Is Fast

Standing Slow-Mode Shocks and Logarithmic Scaling

Step 1. In the Petschek model, the resistive diffusion region has length \(l^* \ll L\). Four standing slow-mode shocks extend from the corners of this compact region, subtending an angle\(\theta\) to the inflow direction.

Step 2. Most of the energy conversion occurs at the shocks, not in the diffusion region. The opening angle \(\theta\) determines the inflow speed: \(v_{\text{in}} \approx v_A\sin\theta\).

Step 3. The maximum rate occurs when the diffusion region length \(l^*\) shrinks to the minimum allowed by the Sweet-Parker scaling within that sub-region. Self-consistently:

$$l^* \sim L\exp\left(-\frac{\pi v_A}{4 v_{\text{in}}}\right)$$

Step 4. Maximizing the inflow rate subject to this constraint:

$$\boxed{\frac{v_{\text{in}}}{v_A}\bigg|_{\max} \approx \frac{\pi}{8\ln R_m} \approx 0.01\text{--}0.03 \text{ for coronal } R_m}$$

The logarithmic dependence on \(R_m\) makes the rate only weakly sensitive to the resistivity, explaining why flare timescales are similar despite large uncertainties in the effective resistivity. Modern simulations with plasmoid instability achieve rates of ~0.01, consistent with observations.

12.8 Power-Law Spectrum from First-Order Fermi Acceleration

Deriving \(f(E) \propto E^{-\delta}\) from Statistical Arguments

Step 1. At each interaction with converging magnetic mirrors (or equivalently, each complete Fermi cycle), a particle gains fractional energy \(\epsilon = \Delta E/E\):

$$E_n = E_0(1+\epsilon)^n \quad\text{after } n \text{ interactions}$$

Step 2. The probability of remaining in the acceleration region after \(n\)interactions is \(P_n = (1-P_{\text{esc}})^n\) where \(P_{\text{esc}}\) is the escape probability per cycle.

Step 3. Eliminating \(n\): from \(n = \ln(E/E_0)/\ln(1+\epsilon)\):

$$N(>E) \propto P_n \propto \left(\frac{E}{E_0}\right)^{\ln(1-P_{\text{esc}})/\ln(1+\epsilon)}$$

Step 4. For small \(\epsilon\) and \(P_{\text{esc}}\):

$$\boxed{f(E) = \frac{dN}{dE} \propto E^{-\delta}, \quad \delta = 1 + \frac{P_{\text{esc}}}{\epsilon}}$$

For the first-order Fermi process at a shock with compression ratio \(r\):\(\epsilon \propto (r-1)/r\) and \(P_{\text{esc}} \propto 1/r\), giving\(\delta = (r+2)/(2(r-1))\) for the energy spectrum. For a strong shock (\(r=4\)):\(\delta = 2\), producing the classic \(E^{-2}\) spectrum. In flares, the observed hard X-ray spectral index of 3-5 corresponds to an electron spectrum of\(\delta \approx 2\text{--}4\), broadly consistent.

12.9 CSHKP Flare Model Diagram

The standard CSHKP model shows the geometry of a two-ribbon flare with reconnection in a vertical current sheet above the magnetic arcade.

ChromospherePost-flare loopsCurrent SheetX-pointInflowPlasmoidRibbonRibbone- beamsHXRHXR

CSHKP standard flare model: magnetic reconnection at the X-point in the current sheet releases energy. Accelerated electrons stream down field lines producing HXR at chromospheric footpoints (ribbons). Heated plasma fills post-flare loops visible in SXR and EUV.

Extended Simulation: Reconnection Rates & GOES Light Curves

Extended: Reconnection Models, GOES Flare Simulation, Current Sheets, Fermi Spectra

Python
script.py141 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Practice Problems

Problem 1:Calculate the Sweet-Parker reconnection rate for a current sheet in the solar corona. Given: Lundquist number $S = R_m = 10^{12}$, Alfven speed $v_A = 1000\;\text{km/s}$, and sheet length $L = 10^9\;\text{cm}$. Find the inflow speed and reconnection timescale.

Solution:

1. Sweet-Parker reconnection rate: $M_A = v_{\text{in}}/v_A = S^{-1/2}$.

2. $M_A = (10^{12})^{-1/2} = 10^{-6}$.

3. Inflow speed: $v_{\text{in}} = M_A \times v_A = 10^{-6} \times 10^8\;\text{cm/s} = 100\;\text{cm/s} = 1\;\text{m/s}$.

4. Sweet-Parker sheet thickness: $\delta = L/\sqrt{S} = 10^9/10^6 = 10^3\;\text{cm} = 10\;\text{m}$.

5. Reconnection timescale: $\tau_{SP} = L/v_{\text{in}} = 10^9/100 = 10^7\;\text{s} \approx 116\;\text{days}$.

6. This is far too slow to explain solar flares ($\tau_{\text{flare}} \sim 10^2\text{-}10^3\;\text{s}$). The Sweet-Parker rate is $\sim 10^{-6}$, but observations require $M_A \sim 0.01\text{-}0.1$. This discrepancy motivates fast reconnection models (Petschek, plasmoid instability, turbulent reconnection).

Problem 2:Compare the Petschek reconnection rate with the Sweet-Parker rate from Problem 1. The Petschek maximum rate is $M_A^{\text{Petschek}} = \pi/(8\ln S)$. Calculate the speedup factor and the resulting flare timescale.

Solution:

1. Petschek rate: $M_A^P = \frac{\pi}{8\ln S} = \frac{\pi}{8\ln(10^{12})} = \frac{\pi}{8 \times 27.63}$.

2. $M_A^P = \frac{3.1416}{221.0} = 0.0142$.

3. Speedup factor: $M_A^P / M_A^{SP} = 0.0142 / 10^{-6} = 1.42 \times 10^4$.

4. Petschek inflow speed: $v_{\text{in}} = 0.0142 \times 10^8 = 1.42 \times 10^6\;\text{cm/s} = 14.2\;\text{km/s}$.

5. Reconnection timescale: $\tau_P = L/v_{\text{in}} = 10^9 / 1.42 \times 10^6 = 704\;\text{s} \approx 12\;\text{min}$.

6. This is consistent with observed impulsive phase durations ($\sim 10$ min). Petschek reconnection achieves fast rates via slow-mode shocks that extend from a compact diffusion region, converting magnetic energy over a large volume without requiring all plasma to pass through the thin resistive layer. Modern simulations show plasmoid instability achieves a universal rate $M_A \sim 0.01$ independent of $S$.

Problem 3:Estimate the total magnetic energy available in a solar active region with field strength $B = 300\;\text{G}$ and volume $V = (3 \times 10^9\;\text{cm})^3$. What GOES flare class could this energy produce?

Solution:

1. Magnetic energy density: $u_B = \frac{B^2}{8\pi}$ (CGS).

2. $u_B = \frac{(300)^2}{8\pi} = \frac{9 \times 10^4}{25.13} = 3581\;\text{erg/cm}^3$.

3. Volume: $V = (3 \times 10^9)^3 = 2.7 \times 10^{28}\;\text{cm}^3$.

4. Total magnetic energy: $E_B = u_B \times V = 3581 \times 2.7 \times 10^{28} = 9.67 \times 10^{31}\;\text{erg}$.

5. Only the free (non-potential) magnetic energy can be released; typically 10-30% of $E_B$. Available energy: $E_{\text{free}} \sim 0.1\text{-}0.3 \times 10^{32} \approx 10^{31}\;\text{erg}$.

6. GOES classification: $10^{31}\;\text{erg}$ corresponds to an X-class flare (X-class flares release $10^{31}\text{-}10^{32}\;\text{erg}$). For comparison: C-class $\sim 10^{29}$, M-class $\sim 10^{30}$, X-class $\sim 10^{31}\text{-}10^{32}\;\text{erg}$. The largest recorded flares ($\sim 10^{33}\;\text{erg}$) require even larger active regions.

Problem 4:Calculate the Alfven speed in the solar corona for a magnetic field $B = 50\;\text{G}$ and electron density $n_e = 10^9\;\text{cm}^{-3}$. How does this compare with the sound speed at $T = 2 \times 10^6\;\text{K}$?

Solution:

1. Alfven speed: $v_A = \frac{B}{\sqrt{4\pi\rho}}$ (CGS), where $\rho = n_e m_p$ for a hydrogen plasma.

2. Mass density: $\rho = 10^9 \times 1.673 \times 10^{-24} = 1.673 \times 10^{-15}\;\text{g/cm}^3$.

3. $v_A = \frac{50}{\sqrt{4\pi \times 1.673 \times 10^{-15}}} = \frac{50}{\sqrt{2.103 \times 10^{-14}}} = \frac{50}{1.450 \times 10^{-7}}$.

4. $v_A = 3.45 \times 10^8\;\text{cm/s} = 3450\;\text{km/s}$.

5. Sound speed: $c_s = \sqrt{2k_BT/m_p} = \sqrt{2 \times 1.381 \times 10^{-16} \times 2 \times 10^6 / 1.673 \times 10^{-24}} = 1.82 \times 10^7\;\text{cm/s} = 182\;\text{km/s}$.

6. Plasma beta: $\beta = c_s^2/v_A^2 = (182/3450)^2 = 0.0028$. The corona is strongly magnetically dominated ($\beta \ll 1$), meaning magnetic forces control the plasma dynamics. Alfven waves propagate ~19 times faster than sound waves, and the magnetic field determines coronal loop structure, flare dynamics, and CME acceleration.

Problem 5:A solar flare releases $E = 10^{32}\;\text{erg}$ in $\tau = 600\;\text{s}$. If the reconnection region has dimensions $L = 5 \times 10^9\;\text{cm}$ and $B = 200\;\text{G}$, estimate the required reconnection rate $M_A$ and compare with Sweet-Parker and Petschek predictions.

Solution:

1. Power: $P = E/\tau = 10^{32}/600 = 1.67 \times 10^{29}\;\text{erg/s}$.

2. Poynting flux into the reconnection region: $P \approx \frac{B^2}{4\pi}v_{\text{in}} L^2$ (energy flux $\times$ area).

3. Energy density: $B^2/(4\pi) = (200)^2/(4\pi) = 3183\;\text{erg/cm}^3$.

4. Required inflow: $v_{\text{in}} = \frac{P}{(B^2/4\pi)L^2} = \frac{1.67 \times 10^{29}}{3183 \times (5 \times 10^9)^2} = \frac{1.67 \times 10^{29}}{7.96 \times 10^{22}} = 2.10 \times 10^6\;\text{cm/s}$.

5. Alfven speed: $v_A = B/\sqrt{4\pi\rho} \approx 4400\;\text{km/s}$ (for $n_e = 10^9$). Reconnection rate: $M_A = v_{\text{in}}/v_A = 21/4400 = 0.0048 \approx 0.005$.

6. Comparison: Sweet-Parker ($M_A \sim 10^{-6}$) is 5000$\times$ too slow. Petschek ($M_A \sim 0.014$) and plasmoid-mediated reconnection ($M_A \sim 0.01$) are consistent with the required rate. This confirms that fast reconnection mechanisms are essential to explain solar flare energetics.

← SunspotsNext: Coronal Mass Ejections →
Share:XRedditLinkedIn
Rate this chapter: