---
title: Kglobal Model for Macroscale Reconnection
url: https://www.emergentmind.com/topics/kglobal-model
type: topic
---

# Kglobal Model for Macroscale Reconnection

Searching arXiv for recent papers on the kglobal model in magnetic reconnection.
The **kglobal model** is a computational model for **macroscale magnetic reconnection** that combines an MHD-like fluid backbone with **guiding-center macro-particles** to study non-thermal particle energization without resolving kinetic scales such as the ion and electron gyro-radii, inertial lengths, or Debye length. It was introduced to address regimes in which fully kinetic PIC simulations are computationally prohibitive, while pure MHD lacks self-consistent particle acceleration. Across its published development, kglobal has been used to study **energetic electron production**, **large-scale parallel electric fields and return currents**, **simultaneous ion and electron energization**, **upstream slow shocks**, and **system-size control of maximum particle energy** during reconnection [1907.07554], [2401.14500], [2109.06351], [2512.13394].

## 1. Definition and physical scope

kglobal is described as a **macroscale reconnection model for energetic electron production**, and later as a framework extended to include **particle ions** as well as electrons. Its basic premise is that in large reconnection systems the dominant energization channels are controlled by **large-scale magnetic geometry**, especially **Fermi reflection in large-scale magnetic fields**, rather than by kinetic-scale boundary layers localized near x-lines and separatrices [1907.07554], [2401.14500].

The model therefore **eliminates kinetic scales** and evolves the plasma self-consistently on fluid scales while advancing energetic particles as **guiding-center macro-particles**. In this formulation, particles move across the magnetic field via $\mathbf{E}\times\mathbf{B}$ drift and along the field at their parallel velocity, while their gyromotion is not resolved. This design makes the model suitable for systems such as solar flares, large magnetotail reconnection regions, and heliospheric current sheets, where the scale separation between the global energy-release region and kinetic plasma scales is extreme [1907.07554], [2508.00119], [2410.16539].

A central consequence of this ordering is that kglobal is not intended as a full kinetic first-principles description of microscopic processes. Instead, it is intended to retain the reconnection dynamics most relevant to **large-scale particle acceleration**, including **flux-rope formation**, **island mergers**, **pressure anisotropy feedback**, and **parallel electric-field effects** that remain important after kinetic layers are removed [1907.07554], [2401.14500].

## 2. Core formulation and dynamical variables

The original kglobal model contains three plasma components: **ion fluid**, **electron fluid**, and **particle electrons**. In the later ion-enabled formulation, a fourth component, **particle ions**, is added. No species conversion is allowed after initialization [2401.14500].

In the electron-only formulation, particle electrons are advanced in the guiding-center limit with conserved magnetic moment
\[
\mu_{ep} = \frac{p_{ep,\perp}^2}{2B} = \text{const.}
\]
and parallel momentum evolution
\[
\frac{d p_{ep,\parallel}}{dt} = p_{ep,\parallel}\mathbf{v}_E\cdot \boldsymbol{\kappa} -\frac{\mu_e}{\gamma_e}\mathbf{b}\cdot\nabla B - eE_{\parallel}.
\]
The terms represent **curvature/Fermi acceleration**, **mirror force**, and **parallel electric-field acceleration**. The same guiding-center structure is later used for particle ions, with
\[
\mu_{ip}=\frac{p_{ip,\perp}^2}{2B}=\text{const.}
\]
and
\[
\frac{d p_{ip,\parallel}}{dt} = p_{ip,\parallel}\mathbf{v}_E\cdot \boldsymbol{\kappa} -\mu_i\,\mathbf{b}\cdot\nabla B +eE_{\parallel}.
\]
This establishes a unified macroscale representation of both ion and electron energization [2401.14500].

The ion fluid obeys a continuity equation and an MHD-like momentum equation coupled to particle stresses and the electromagnetic field. In the original model the momentum equation includes ion-fluid stress, electron-fluid stress, particle-electron stress, the Lorentz force, and the field-aligned electric force. The magnetic field evolves through Faraday’s law using the perpendicular electric field,
\[
\frac{\partial \mathbf{B}}{\partial t} = -c\,\boldsymbol{\nabla\times}\mathbf{E}_{\perp}, \qquad \mathbf{E}_{\perp}=-\frac{1}{c}\mathbf{v}_i\times\mathbf{B}.
\]
The parallel electric field is determined from electron pressure and stress terms rather than from kinetic boundary-layer dynamics [2401.14500].

The model uses **quasi-neutrality** and a **zero-parallel-current closure**. In the original formulation,
\[
n_{ef}=n_i-n_{ep},
\]
and
\[
n_{ef}v_{ef,\parallel}=n_iv_{i,\parallel}-n_{ep}v_{ep,\parallel},
\]
with the ordering
\[
\frac{V_d}{C_A}=\frac{J}{neC_A}\sim \frac{d_i}{L}\ll 1.
\]
In the ion-extended version, these relations become
\[
n_{ef}=n_{if}+n_{ip}-n_{ep},
\]
and
\[
n_{ef}v_{ef,\parallel}=n_{if}v_{if,\parallel}+n_{ip}v_{ip,\parallel}-n_{ep}v_{ep,\parallel}.
\]
These closures encode the macroscale assumption that current-drift scales have been ordered out, while field-aligned flows remain self-consistent [2401.14500].

## 3. Large-scale parallel electric field and energy-conserving closure

A major development of kglobal was the explicit inclusion of a **large-scale parallel electric field** produced by **magnetic-field-aligned gradients in electron pressure**. This extension distinguishes between two different classes of field-aligned electric field: localized **kinetic-scale parallel electric fields** near the x-line and separatrices, and a **macroscale $E_{\parallel}$** that acts over the reconnection exhaust and is dynamically important for transport, early-stage heating, and return-current closure [1907.07554].

In the upgraded formulation, the large-scale parallel electric field is
\[
E_{\parallel}=\frac{-1}{n_ie}\left( \boldsymbol{B} \cdot \boldsymbol{\nabla} \left( \frac{m_en_c v_{\parallel c}^2}{B} \right) + \boldsymbol{b}\cdot \boldsymbol{\nabla} P_{c} + \boldsymbol{b} \cdot \boldsymbol{\nabla} \cdot \boldsymbol{T}_h\right),
\]
where $n_i=n_c+n_h$, $n_c$ and $n_h$ are cold and hot electron densities, $v_{\parallel c}$ is the cold-electron parallel flow, $P_c$ is the cold-electron scalar pressure, $\boldsymbol{b}=\boldsymbol{B}/B$, and $\boldsymbol{T}_h$ is the hot-electron gyrotropic stress tensor. The latter is written as
\[
\boldsymbol{T}_{h}=T_{eh\parallel}\boldsymbol{bb}+P_{eh\perp}(\boldsymbol{I}-\boldsymbol{bb}),
\]
with
\[
P_{eh\perp}=\int d\boldsymbol{p}_{e}\frac{p_{e\perp}^2}{2m_e \gamma_e}f, \qquad T_{eh\parallel}=\int d\boldsymbol{p}_{e}\frac{p_{e\parallel}^2}{m_e \gamma_e}f.
\]
The normalization of this field is
\[
E_\parallel \sim \frac{m_e C_{Ae}^2}{eL_0} = \frac{m_i C_A^2}{eL_0},
\]
and
\[
\frac{E_\parallel}{E_\perp}\sim \frac{d_i}{L_0}\ll 1.
\]
Accordingly, $E_\parallel$ is retained for motion along the field but neglected in Faraday’s law for the evolution of the magnetic field [1907.07554].

This extension also introduces **cold-electron return currents** in open systems. The cold-electron parallel speed is constrained by
\[
v_{||c}=\frac{1}{n_c}(n_iv_{||i}-n_hv_{||h}),
\]
so that escaping hot-electron current is balanced by a compensating cold-electron flow. The model therefore represents large-scale current closure and suppression of hot-electron escape without reintroducing kinetic boundary layers [1907.07554].

A defining feature of this formulation is that the added $E_\parallel$ is derived in an **energy-conserving** manner. The resulting conservation law is
\[
W_{MHD}+W_c+W_h=constant,
\]
with cold-electron energy
\[
W_c=\frac{m_en_cv_{\parallel c}^2}{2}+\frac{1}{\Gamma -1}P_c.
\]
A later ion-inclusive formulation extends this conservation structure to
\[
\int d\mathbf{x}\, \frac{\partial}{\partial t} \left(W_{if}+W_{ip}+W_{ef}+W_{ep}+B^2/8\pi\right)=0,
\]
up to small errors associated with formally neglected terms. This places energy conservation at the center of the model’s reduced-physics design [1907.07554], [2401.14500].

## 4. Extension from electron-only to ion-electron energization

The paper “A Computational Model for Ion and Electron Energization during Macroscale Magnetic Reconnection” extends kglobal from an **electron-only kinetic-fluid hybrid** to a model with **both ions and electrons** as particle species [2401.14500]. The principal structural change is the inclusion of **particle-ion inertia** in the fluid momentum equation.

In the original model, electron inertia could be neglected in the ion-fluid momentum balance. This is no longer possible once particle ions are added. The upgraded fluid-ion momentum equation becomes
\[
\begin{split}
m_i&\frac{\partial\, (n_{i}\mathbf{v}_{if})}{\partial t} = - \boldsymbol{\nabla \cdot}\mathbb{T}_{if} -(\boldsymbol{\nabla \cdot}\mathbb{T}_{ef})_\perp + \mathbf{J}\boldsymbol{\times}\mathbf{B}/c  \\
&- \boldsymbol{\nabla\cdot}\mathbb{T}_{ip} -(\boldsymbol{\nabla \cdot}\mathbb{T}_{ep})_{\perp} +en_iE_{\parallel}\mathbf{b} -\mathbf{I}_{i\parallel},
\end{split}
\]
with the correction term
\[
\mathbf{I}_{i\parallel} = \frac{\partial}{\partial t}m_in_{ip}(v_{ip\parallel}-v_{if\parallel})\mathbf{b}.
\]
This term accounts for the difference in parallel inertia between fluid and particle ions [2401.14500].

The particle-ion stress tensor is gyrotropic,
\[
\mathbb{T}_{ip} = P_{ip\parallel}\mathbf{b}\mathbf{b} + P_{ip\perp}(\mathbb{I}-\mathbf{b}\mathbf{b}) + m_in_{ip}\mathbf{v}_{ip}\mathbf{v}_{ip},
\]
with
\[
P_{ip\perp}=\int d\mathbf{p}_i \frac{p^2_{i\perp}}{2m_i}f, \qquad P_{ip\parallel}=\int d\mathbf{p}_i \frac{p^2_{i\parallel}}{m_i}f.
\]
The electron-fluid and ion-fluid pressures use adiabatic closures,
\[
\frac{d}{dt}\left(\frac{P_{if}}{n_{if}^{\gamma}}\right)=0, \qquad \frac{d}{dt}\left(\frac{P_{ef}}{n_{ef}^{\gamma}}\right)=0,
\]
with $\gamma=5/3$ [2401.14500].

This extension broadens the model’s scope from energetic electron production to the **simultaneous non-thermal energization of ions and electrons**. A plausible implication is that kglobal becomes suitable not only for electron-acceleration problems, but also for **energy partition** studies in large reconnection systems where proton and electron heating evolve differently.

## 5. Dominant acceleration channels and species dependence

Across the kglobal literature, the dominant acceleration mechanism is repeatedly identified as **Fermi reflection** associated with **contracting magnetic structures**, **reconnection exhausts**, and especially **flux-rope or magnetic-island mergers**. In the model’s physical picture, particles reflect from moving magnetic structures whose motion is tied to Alfvénic reconnection outflows, so the energization rate can be proportional to particle energy [1907.07554], [2512.13394], [2508.00119].

In the study of maximum particle energy, one merger of two equal flux ropes of radius $r_0$ is described by the field-line shortening
\[
s_0 = 4\pi r_0, \qquad s = 2\sqrt{2}\,\pi r_0 = \frac{s_0}{\sqrt{2}},
\]
which yields a $\sqrt{2}$ increase in $v_\parallel$ and therefore approximately
\[
W \rightarrow 2W.
\]
After $N$ comparable mergers,
\[
W = W_i 2^N.
\]
The number of mergers is linked to the effective system size through hyper-resistive control of the smallest island width,
\[
w_0 \sim L S_\nu^{-1/4},
\]
with
\[
S_\nu=\frac{C_A L^3}{\nu},
\]
so that when the largest island reaches system scale,
\[
2^N \sim \left(\frac{L}{w_0}\right)^2 \sim S_\nu^{1/2},
\]
and therefore
\[
W_{\max}\sim W_i S_\nu^{1/2}.
\]
The simulations vary $S_\nu$ from $1.2\times10^7$ to $6.1\times10^9$ and find that the inferred merger count rises from about $N\approx 7.95$ to $N\approx 9.64$, while the high-energy cutoff shifts to higher energies as $S_\nu$ increases [2512.13394].

The model has also been used to explain why **proton heating and energization exceed those of electrons** when the upstream temperatures of the two species are equal. In that study, particles entering a reconnection exhaust gain a velocity increment of order
\[
\Delta v \sim 2C_A,
\]
so that for a single reflection
\[
\Delta \mathcal{E} = \frac{1}{2}m\left[(v_{\parallel,0}+2C_A)^2-v_{\parallel,0}^2\right] = 2mC_A v_{\parallel,0}+2mC_A^2.
\]
For protons, the initial gain scales like
\[
\Delta \mathcal{E}_i \sim m_i C_A^2.
\]
For electrons, the initial energy gain is smaller and scales as
\[
\Delta \mathcal{E}_e \sim \left(\beta_{e0}\frac{m_e}{m_i}\right)^{1/2} m_i C_A^2 \ll m_i C_A^2.
\]
The simulations report late-time temperature increments nearly insensitive to mass ratio,
\[
\Delta T_e \approx 0.070,\ 0.079,\ 0.071,\qquad \Delta T_i \approx 0.313,\ 0.306,\ 0.314
\]
for $m_i/m_e=25,\ 100,\ 400$, respectively. This is interpreted as evidence that the relative heating is controlled more by the early injection energy and the large-scale Fermi process than by the artificial mass ratio [2508.00119].

Guide-field strength modifies the efficiency of this process. A strong guide field increases the curvature radius of reconnected field lines and suppresses Fermi reflection, reducing energy gain. In the maximum-energy study the guide field is held fixed so that the principal scaling remains the system-size dependence $W_{\max}\propto S_\nu^{1/2}$ [2512.13394].

## 6. Reconnection structures, shocks, and benchmarks

kglobal has been used not only for non-thermal tails, but also for **shock formation** and **anisotropic MHD-like wave physics**. In “Slow Shock Formation Upstream of Reconnecting Current Sheets,” slow shocks are documented in the upstream region of simulations with the **kglobal kinetic macroscale simulation model**. During multi-island reconnection, the formation and merging of flux ropes drives plasma flows and pressure disturbances in the upstream region. These disturbances steepen into slow shocks that propagate along the reconnecting component of the magnetic field and satisfy the expected **Rankine-Hugoniot jump conditions**. Plasma heating arises from both compression across the shock and the parallel electric field that develops to maintain charge neutrality in a kinetic system. The shocks are weaker at lower plasma $\beta$, where shock steepening is slow, and their contribution to electron heating is described as relatively minor compared with **Fermi reflection** and the parallel electric fields that bound the reconnection outflow [2109.06351].

The model has also been benchmarked against wave phenomena that probe its anisotropic closures. In the large-scale $E_\parallel$ paper, kglobal accurately captures the damping of **electron acoustic modes** in a plasma with hot kinetic electrons and cold fluid electrons. For the perturbed distribution,
\[
\partial_t \tilde{f} + v_{||}\nabla_{||}\tilde{f} - \frac{e}{m_e}\tilde{E}_{||}\partial_{v_{||}}f_0 = 0,
\]
the perturbed parallel electric field is
\[
\tilde{E}_\parallel=-\frac{1}{n_ie} \left( \frac{5}{3} \nabla_\parallel T_{c}\tilde{n}_{h} + \nabla_\parallel \tilde{T}_{h} \right),
\]
and the dispersion relation is
\[
\frac{n_{0c}}{n_{0h}}=Z'(\zeta)\left(\frac{5}{6}\frac{T_{0c}}{T_{0h}}+\zeta^2\right).
\]
The reported damping rates agree very well with the linear theory. The same paper benchmarks the suppression of hot-electron transport by comparing kglobal with the PIC code **p3d**, finding that **kglobal with $E_\parallel$** matches p3d well over most of the domain, whereas **kglobal without $E_\parallel$** spreads much more rapidly [1907.07554].

The ion-inclusive formulation is benchmarked using a **circularly polarized Alfvén wave in anisotropic plasma** and the **firehose instability**. The predicted phase speed is
\[
C_p = C_A \sqrt{1 - \frac{4\pi(P_{\parallel} - P_{\perp})}{B^2}} \equiv \alpha C_A,
\]
and the firehose growth rate is
\[
\gamma = k C_A |\alpha| - \nu k^4,
\]
with $\nu = 5.0\times 10^{-5}$ and a test case using $\alpha^2=-0.2$. The measured phase speeds and growth rates are reported to agree well with theory, and no spurious short-wavelength instabilities appear [2401.14500].

## 7. Applications, observational interpretation, and limitations

The kglobal model has been used as an interpretive framework for direct spacecraft observations. In a Parker Solar Probe study of the near-Sun heliospheric current sheet, kglobal is described as a **global reconnection modeling framework** upgraded to handle **ion energization as well as electron energization**, with all kinetic scales eliminated. Simulations of a Harris-like equilibrium use upstream parameters based directly on observations: $B_1 = 510\ \mathrm{nT}$, $B_2 = 470\ \mathrm{nT}$, $n_1 = 1500\ \mathrm{cm^{-3}}$, $n_2 = 3000\ \mathrm{cm^{-3}}$, proton temperature $T_p = 46\ \mathrm{eV}$, electron temperature $T_e = 91\ \mathrm{eV}$, and guide fields $B_g = 0.2\,B_{\rm rec}$ and $B_g = 0.3\,B_{\rm rec}$ [2410.16539].

The hybrid Alfvén speed is
\[
C_{Ah} = \left[\frac{B_1 B_2 (B_1 + B_2)}{4\pi m_i (n_1 B_2 + n_2 B_1)}\right]^{1/2},
\]
with
\[
C_{Ah} \approx 232\ \mathrm{km\ s^{-1}},
\qquad
m_i C_{Ah}^2 \approx 0.56\ \mathrm{keV}.
\]
Despite this modest magnetic energy per particle, the simulations produce **protons up to $\sim 500\ \mathrm{keV}$** for $B_g = 0.2\,B_{\rm rec}$ and **electrons up to $\sim 100\ \mathrm{keV}$**, while the observed proton spectrum between about $67$ and $527\ \mathrm{keV}$ is fitted by
\[
\frac{dj}{dE} \propto E^{-\gamma}, \qquad \gamma \approx 5.1 \pm 0.12.
\]
The $B_g = 0.3\,B_{\rm rec}$ simulation gives a softer spectrum with power-law index $\sim -5$ and maximum proton energy $\sim 200\ \mathrm{keV}$, and is identified as the best match to the observed slope. The study argues that a guide field near $\sim 0.25\,B_{\rm rec}$ would likely match both slope and cutoff more closely [2410.16539].

These applications reinforce a recurring interpretive theme: kglobal attributes extended power-law tails and large maximum energies to **repeated Fermi acceleration in merging magnetic islands or flux ropes**, while allowing pressure anisotropy and firehose feedback to regulate the reconnection dynamics. This suggests that the model is especially effective when the dominant physics is **macroscale island-driven acceleration** rather than kinetic boundary-layer processes or strong turbulence.

The model’s limitations are stated explicitly in the literature. kglobal is **2D** in several published applications, removes kinetic scales, uses **guiding-center dynamics**, employs **hyper-resistivity** and, in some studies, an artificial **proton-to-electron mass ratio of 25**. The upper energy cutoff can also be controlled by the **size of the simulation domain**, so simulated maxima are not always pure physics cutoffs. The model does not support **plasma waves** that require violation of charge neutrality, and it is not designed to capture fine kinetic boundary layers or short-wavelength plasma oscillations [1907.07554], [2410.16539].

A common misconception is that kglobal is simply a cheaper PIC surrogate. The published descriptions do not support that characterization. Instead, kglobal is a **reduced-kinetic macroscale model** built around a specific claim: that the high-energy particle population in large reconnection systems is controlled mainly by **large-scale Fermi processes**, **parallel electric fields from pressure gradients**, and **pressure-anisotropy feedback**, rather than by directly resolving kinetic boundary layers [1907.07554], [2401.14500].

Source: https://www.emergentmind.com/topics/kglobal-model