Jeongbhin Seo

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.

We investigate binary mergers of galaxy clusters, the formation of shocks, and the resulting radio relics using three-dimensional simulations. The initial setup consists of two idealized spherical subclusters with a mass ratio below 3, each permeated by turbulent magnetic fields, and we follow their merger with a high-order accurate magnetohydrodynamic (MHD) code. In parallel, we track the acceleration of cosmic-ray electrons (CRe) via diffusive shock acceleration (DSA) at merger-driven shocks, together with radiative cooling and Fermi-II (turbulent) acceleration in the postshock region, employing a newly developed Eulerian Fokker–Planck solver. Synchrotron emission is computed from the simulated CRe distribution and magnetic fields. In this paper, we detail these numerical approaches and present the first results obtained with them. Two prominent axial shocks emerge along the merger axis; the shock ahead of the heavier subcluster systematically attains a higher Mach number, although it is more compact, than that ahead of the lighter subcluster. Turbulent magnetic fields, which are both inherited from the initial conditions and amplified during the merger, produce patchy, fine-scale structures in the radio surface brightness. Because of the combined effects of turbulent acceleration, spatially nonuniform magnetic fields, and the curved geometry of merger shocks, the volume-integrated radio spectra show deviations from the canonical power-law steepening expected for a planar shock with a uniform field. Reacceleration of preexisting fossil CRe enhances the surface brightness. Our results highlight the coupled roles of merger dynamics, MHD turbulence, and CRe physics in shaping the observed properties of radio relics in cluster outskirts.

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.

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.

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.

Recent observations and simulations indicate that solar flares undergo extremely complex 3D evolution, making 3D particle transport models essential for understanding electron acceleration and interpreting flare emissions. In this study, we investigate this problem by solving Parker’s transport equation with 3D MHD simulations of solar flares. By examining energy conversion in the 3D system, we evaluate the roles of different acceleration mechanisms, including reconnection current sheet (CS), termination shock (TS), and supra-arcade downflows (SADs). We find that large-amplitude turbulent fluctuations are generated and sustained in the 3D system. The model results demonstrate that a significant number of electrons are accelerated to hundreds of keV and even a few MeV, forming power-law energy spectra. These energetic particles are widely distributed, with concentrations at the TS and in the flare looptop region, consistent with results derived from recent hard X-ray (HXR) and microwave (MW) observations. By selectively turning particle acceleration on or off in specific regions, we find that the CS and SADs effectively accelerate electrons to several hundred keV, while the TS enables further acceleration to MeV. However, no single mechanism can independently account for the significant number of energetic electrons observed. Instead, the mechanisms work synergistically to produce a large population of accelerated electrons. Our model provides spatially and temporally resolved electron distributions in the whole flare region and at the flare footpoints, enabling synthetic HXR and MW emission modeling for comparison with observations. These results offer important insights into electron acceleration and transport in 3D solar flare regions.

Relativistic jets of radio galaxies (RGs) are possible sources of ultra-high-energy cosmic rays (UHECRs). Recent studies combining relativistic hydrodynamic simulations with Monte Carlo particle transport have demonstrated that UHECRs can be accelerated to energies beyond 10²⁰ eV through shocks, turbulence, and relativistic shear in jet-induced flows of Fanaroff–Riley type RGs. The resulting time-asymptotic UHECR spectrum is well modeled by a double power law with an “extended” exponential cutoff, primarily shaped by relativistic shear acceleration. In this study, we adopt this novel source spectrum and simulate the propagation of UHECRs from nearby RGs using the CRPropa code. We focus on Virgo A (Vir A), Centaurus A (Cen A), Fornax A (For A), and Cygnus A (Cyg A), expected to be the most prominent UHECR sources among RGs. We then analyze the energy spectrum and mass composition of UHECRs arriving at Earth. We find that, due to the extended high-energy tail in the source spectrum, UHECRs from Vir A, which has a higher Lorentz factor, exhibit a higher flux at the highest energies and a lighter mass composition at Earth compared to those from Cen A and For A with lower Lorentz factors. Despite Cyg A having an even higher Lorentz factor, the large distance limits its contribution. With a small number of nearby prominent RGs, our findings suggest that if RGs are the major sources of UHECRs, the energy spectrum and mass composition of observed UHECRs would exhibit hemispheric differences between the Northern and Southern skies at the highest energies.

Magnetic reconnection and turbulence in magnetically-dominated environments have been proposed as important nonthermal particle acceleration mechanisms that generate high energy particles and associated emissions. While the acceleration to high energy that produces the power-law energy distribution has drawn strong interest, recent studies actively discuss pre-acceleration, or injection, to a sufficient energy for a sustained and prolonged Fermi-like acceleration. The injection process is important for determining the fraction of nonthermal particles and energy partition between thermal and nonthermal particles. We review recent advances in understanding the injection mechanisms responsible for populating these nonthermal power-law spectra, and conclude with an outlook for studies and applications of injection models.