Interstellar Medium Modulation of Nonlinear Kinetic Alfvén Morphology in Structured Galactic Environments

TL;DR

This paper develops a spatially dependent framework for nonlinear KA solitons in the magnetized interstellar medium by coupling a diffuse WIM background with embedded H II regions, SWBs, and SNRs. Using a multi-component analytic ISM model and κ-distributed electrons/positrons, it derives location-dependent KdV coefficients that govern KA dispersion and nonlinearity, enabling Galactic maps of soliton existence, amplitude, and width. The study identifies exclusion zones in high-$\beta$ H II regions and ultra-low-$\beta$ SNR interiors, while SWB and SNR shells provide favorable environments for KA solitons, with sharp property variations near structure boundaries. These results connect large-scale ISM structure to ion-kinetic dissipation, offering potential observables for radio scattering and pulsar scintillation, and providing a physically motivated sub-grid framework for galaxy-scale simulations of wave dissipation and transport.

Abstract

We present a spatially dependent framework for the existence and propagation of nonlinear kinetic Alfvén (KA) structures in the interstellar medium (ISM). Using a multi-component analytical model that incorporates the diffuse warm ionized medium together with localized H II regions, supernova remnants (SNR), and stellar-wind bubbles (SWB), we derive location-dependent coefficients governing KA dispersion and nonlinearity. The reductive perturbation method is applied to obtain Korteweg-de Vries (KdV) equations, enabling the characterization of solitons under realistic astrophysical conditions. Numerical analysis demonstrates how superthermality, plasma $β$, temperature, and density gradients modulate soliton amplitude, width, and stability. Our results reveal distinct exclusion zones (EZs) for KA solitons in high-$β$ HII regions and SWB/SNR interiors, as well as ultra low-$β$ regions near central pulsar wind nebulae. While H II regions exhibit simple Gaussian-driven depletions, the complex ``hole-and-shell" morphologies of SWBs and SNRs imprint sharp spatial variations and discontinuities on soliton properties. This study establishes a direct link between macroscopic ISM morphology, ion-kinetic scale dissipation, and the emergence of coherent Alfvénic activity, with implications for radio scattering, pulsar scintillation, and fine-scale signatures in astrophysical observations.

Figures (6)

• Figure 1: Contour plots of key ISM parameters in the local galactic disk plane near embedded astronomical structures. Radius of magenta dashed circles represent the characteristic size of the embedded structures. Panel (a) shows $\log_{10}$ of total magnetic field strength $B_0$($\mu$G), featuring diamagnetic depletion in the H II region and shell like compressions at the SWB and SNR. Panel (b) shows $\log_{10}$ of total ion number density $n_{i0}$ (cm$^{-3}$). Panel (c) shows $\log_{10}$ of electron temperature $T_{e0}$ (K). Panel (d) shows $\log_{10}$ of plasma $\beta$. High-$\beta$ ($\beta > 1$) regions are masked in cyan, while ultra low-$\beta$ ($\beta < m_e/m_i$) regions are masked in green. The H II region is centered at $(R, Z) \approx (4,0)\,$kpc, SWB at $(7,0)\,$kpc and SNR at $(10,0)\,$kpc. Note that the background field in panel (a) is normalized to $\mu$G to match the 1D profiles in Figs. \ref{['fig: swb']} and \ref{['fig: snr']}.
• Figure 2: Detailed one-dimensional profiles of ISM parameters through the center of the SWB at $(R, Z) \approx (7, 0)$ kpc. (a) Total magnetic field $B_0$ ($\mu$G), showing slight compression at the shell but minimal interior perturbation. (b) Number densities ($\text{cm}^{-3}$) for electrons ($n_{e0}$, solid green), ions ($n_{i0}$, dotted blue), and positrons ($n_{p0}$, dashed red). Note the deep central cavity where ion density drops significantly. (c) Electron temperature $T_{e0}$ (K), exhibiting the characteristic "flat-top" plateau of the hot shocked wind ($T \sim 10^6$ K).
• Figure 3: Detailed one-dimensional profiles of ISM parameters through the center of the SNR at $(R, Z) \approx (10, 0)$ kpc. (a) Total magnetic field $B_0$ ($\mu$G), highlighting the twin peaks of the shock-compressed shell and the central spike corresponding to the PWN. (b) Multi-species densities ($\text{cm}^{-3}$): The ion density ($n_{i0}$, dotted blue) forms a shell but vanishes in the core, while the positron density ($n_{p0}$, dashed red) rises centrally, creating a pair-dominated interior. (c) Electron temperature $T_{e0}$ (K), showing intense shock heating at the shell.
• Figure 4: Two-dimensional Galactic maps of KA soliton (a) amplitude $\psi_0$ and (b) width W, at $\kappa_e=150, \kappa_p=100$, $\theta=85^\circ$, $U=-0.06$ and $\sigma=1$. The parameters for the background (WIM) and the analytical profiles of embedded structures are given in Table \ref{['tab:param_background']} and Table \ref{['tab:profiles_parameters']}.
• Figure 5: One-dimensional profiles of plasma $\beta$ in (a), amplitude in (b) and width in (c) for $Z=0\,$kpc at SWB for same parameters as in Fig. \ref{['fig: full_map_amp_width']}.