Fan Guo

The mechanism for accelerating energetic particles in the outer heliosphere is still under debate. The termination shock, turbulence, and magnetic reconnection have all been proposed as viable mechanisms. We investigate compression-driven particle acceleration during magnetic reconnection in the heliospheric current sheet using the MHD–SDE framework that incorporates energetic particle feedback into the MHD evolution. For the heliosheath conditions considered here, we find that reconnection-driven acceleration is inefficient downstream of the termination shock, whereas reconnection near the heliopause can generate stronger compression and produce a power-law differential intensity spectrum J ( ε ) ∝ ε ^−5.4 . However, even the hardest spectrum achieved during magnetic reconnection remains much steeper than the typical J ( ε ) ∝ ε ^−1.5 spectrum observed in the heliosheath. This indicates that magnetic reconnection alone is insufficient to explain the acceleration of energetic particles in the heliosheath.

Multi-messenger, multi-viewpoint, and time-resolved observations of solar flares are now providing unprecedented constraints on particle acceleration sites, energy conversion, and energy transport. The interpretation of current observations, including microwave imaging spectroscopy from EOVSA, hard x-ray (HXR) imaging from Solar Orbiter/STIX, gamma-ray diagnostics from Fermi, and in situ measurements from Parker Solar Probe and Solar Orbiter, collectively demands modeling frameworks that go beyond traditional spatially unresolved, one-zone models or single-mechanism descriptions. This review surveys multiscale and multidimensional modeling approaches, including kinetic, magnetohydrodynamic (MHD), and macroscopic particle models, that are being developed to meet the need. Kinetic simulations reveal that three-dimensional (3D) effects, including field-line chaos and self-generated turbulence, are essential for sustained power-law particle acceleration. MHD simulations now capture flux-rope eruptions, plasmoid-unstable current sheets, and turbulent flare regions in realistic magnetic topologies. Macroscopic models coupling MHD with energetic-particle models produce spatially resolved electron distributions and synthetic HXR and microwave emissions for direct comparison with observations. Despite these advances, outstanding challenges remain in bridging kinetic and global scales, improving MHD simulations and macroscopic particle models, and achieving quantitative model-observation closure.

Certain solar flares exhibit a distinctive candle-flame or cusp-shaped feature above the bright flare arcade visible in extreme ultraviolet (EUV) and X-ray channels sensitive to high-temperature plasma. The presence of a cusp-like structure is generally regarded as a key piece of morphological evidence for magnetic reconnection to power explosive energy release in solar flares. In addition, downward-propagating plasma flows above the flare arcade have often been interpreted as outflows driven by magnetic reconnection. However, the relationship between the observed candle-flame-shaped morphology and the underlying magnetic field geometry for reconnection remains unclear. Likewise, the observed speed of the plasma downflows has been found to be too low compared to the upstream Alfvén speed predicted by reconnection theories. With the help of a recently developed three-dimensional magnetohydrodynamics model, we examine the locations where magnetic topology changes from antiparallel to closed (Y-points) in a candle-flame-shaped flare, compare the observational emission features with synthetic EUV images generated from the model, and analyze their time evolutions. We also investigate the role of projection effects and line-of-sight integration in the measurements of plasma downflow speeds. Our analysis reveals that the Y-points do not necessarily coincide with the apparent cusp tip. Also, the apparent speeds of the supra-arcade downflows, as derived from tracks in the time–distance plots, underestimate the true Alfvén speeds in the reconnection inflow region by at least a factor of 2 up to an order of magnitude.

During near-Sun crossings of the heliospheric current sheet (HCS), Parker Solar Probe (PSP) observed populations of high-energy protons and heavier ions, indicating possible energization by magnetic reconnection up to 10 s–100 s keV nucleon ⁻¹ . Here we study ion acceleration by magnetic reconnection at the HCS. To estimate ion energization, we solve the Parker transport equation coupled to a large-scale 2D MHD reconnection simulation. We find that multiple ion species develop power-law distributions with both spectral index and high-energy cutoff E _max consistent with in situ data. By accounting for the injection physics determined by kinetic simulations, we confirm that the charge-to-mass ratio scales as E _max ∝ ( Q / M ) ^α with α ∼ 0.8–1.1, approximately consistent with PSP measurements in the broader range α ∼ 0.6–1.7. In the limit where ions are injected at the same energy per nucleon, α can be as low as ∼0.3. These findings further support the role of magnetic reconnection in producing high-energy heavy ions at the HCS.

Identifying the physical mechanism driving blazar flares remains a central challenge in high-energy astrophysics. We show that the energy dependence of the standard deviation of the polarization angle variability ( σ _PA ) provides a powerful and robust discriminator of blazar flaring mechanisms. Using particle-in-cell-integrated polarized radiative transfer simulations, we perform the most rigorous statistical analyses of polarization variability to date. We demonstrate that magnetic reconnection and magnetized turbulence imprint qualitatively distinct energy dependence of σ _PA that directly reflects their different magnetic field evolution and particle transport. Reconnection predicts higher σ _PA with higher photon energy until the synchrotron spectral peak, whereas turbulence produces nearly flat σ _PA across the synchrotron spectral component. These trends are resilient to realistic observational limitations. Applying our results to optical and IXPE data of Mrk 421 and 1ES 1959+650, we find strong evidence for reconnection-driven flares embedded in a turbulent blazar zone. Energy-dependent σ _PA emerges as a decisive new probe of particle acceleration in relativistic jets.

In situ observations of energetic particles associated with interplanetary (IP) shocks, known as energetic storm particle (ESP) events, often show time-intensity profiles quite different from predictions of classical diffusive shock acceleration theory. We use numerical simulations, including test-particle simulations and hybrid simulations (with fluid electrons and kinetic protons) for shocks propagating into a turbulent magnetic field to study the ESP intensity-time profiles across a strong IP shock. We find that several types of energetic particle intensity profiles similar to in situ observations can be produced in our simulations. The peak of energetic particle count is often near, but can be shifted from, the locations of the shock front, with the peak usually being observed downstream of the shock. These findings may help understand particle acceleration at both traveling IP shocks and shocks in other heliospheric and astrophysical environments.

Solar flares are among the most dramatic events in the solar system, releasing substantial magnetic energy and accelerating a large number of electrons to high energies. Notably, in certain events, the above-the-looptop region may contain a significant population of nonthermal electrons, both in number and energy. For the first time, we adopt a novel numerical method that combines magnetohydrodynamics with energetic particles incorporating feedback from nonthermal electrons to investigate electron acceleration and transport in solar flares. We find that a large fraction of energetic electrons are accelerated via the current sheet and termination shock regions. Most energetic electrons are concentrated in the above-the-looptop region, carrying a sizable amount of the released energy. We observe that greater feedback of nonthermal electrons leads to steeper energy spectra. The energy density of the nonthermal electrons oscillates due to the periodic impact of magnetic islands into the above-the-looptop region, which may help explain the observed quasiperiodic pulsations. Our simulations provide new insights into the origin of nonthermal electrons and associated emissions in the above-the-looptop region.

On 2022 September 5, a large solar energetic particle (SEP) event was detected by Parker Solar Probe (PSP) and Solar Orbiter (SolO) at heliocentric distances of 0.07 and 0.71 au, respectively. PSP observed an unusual velocity dispersion signature: particles below ∼1 MeV exhibited a normal velocity dispersion, while higher-energy particles displayed an inverse velocity arrival (IVA) feature, with the most energetic particles arriving later than those at lower energies. The maximum energy increased from about 20–30 MeV upstream to over 60 MeV downstream of the shock. The arrival of SEPs at PSP was significantly delayed relative to the expected onset of the eruption. In contrast, SolO detected a typical large SEP event characterized by a regular velocity dispersion at all energies up to 100 MeV. To understand these features, we simulate particle acceleration and transport from the shock to the observers with our newly developed SEP model—Particle ARizona and MIchigan Solver on Advected Nodes. Our results reveal that the IVA and delayed particle onset detected by PSP originate from the time-dependent diffusive shock acceleration processes. After shock passage, PSP’s magnetic connectivity gradually shifted due to its high velocity near perihelion, detecting high-energy SEPs streaming sunward. Conversely, SolO maintained a stable magnetic connection to the strong shock region where efficient acceleration was achieved. These results underscore the importance of spatial and temporal dependence in SEP acceleration at interplanetary shocks and provide new insights to understand SEP variations in the inner heliosphere.

In this study, we analyze an important property of the very local interstellar medium, i.e., turbulence that affects the intensity and shape of the Interstellar Boundary Explorer (IBEX) ribbon. Specifically, we simulate the propagation of the ribbon’s parent ions before they become secondary energetic neutral atoms that can be observed at 1 au by IBEX. We then test how different slab/2D turbulence fractions affect the intensity and shape of the modeled ribbon and compare to the latest IBEX ribbon-separated data. We compare 1D cuts across the modeled and observed ribbon as a function of angle away from the ribbon center. We find that the intensities of the modeled ribbon are larger than most of the data, except near the ecliptic plane. However, the model intensities are sensitive to how we model the neutral solar wind, which forms the source ion population for the ribbon. We then compare the model and data by normalizing the fluxes to their respective peak intensities and find the peaks’ angular distances from the ribbon center. We find that most of the model peaks lie within 3° of the data peaks, which we consider our 1σ uncertainty (half the size of an IBEX pixel). By averaging the instances where the model peaks match the data peaks (within 3°), we find that the mean is slab-50%, with a standard deviation of ±28% and standard error of ±5%.