Search

Nonlinear dynamics of runaway electron induced wave instabilities can significantly modify the runaway distribution critical to tokamak operations. Here we present the first-ever fully kinetic simulations of runaway-driven instabilities towards nonlinear saturation in a warm plasma as in tokamak start up. It is found that the slow-X modes grow an order of magnitude faster than the whistler modes, and they parametrically decay to produce whistlers much faster than those directly driven by runaways. These parent-daughter waves, as well as secondary and tertiary wave instabilities, initiate a chain of wave-particle resonances that strongly diffuse runaways to the backward direction. This reduces almost half of the current carried by high-energy runaways, over a time scale orders of magnitude faster than experimental shot duration. These results beyond quasilinear analysis may impact anisotropic energetic electrons broadly in laboratory, space and astrophysics.

Primordial magnetic fields (PMFs) play a pivotal role in influencing small-scale fluctuations within the primordial density field, thereby enhancing the matter power spectrum within the context of the $\Lambda$CDM model at small scales. These amplified fluctuations accelerate the early formation of galactic halos and stars, which can be observed through advanced high-redshift observational techniques. Therefore, Stellar Mass Density (SMD) observations, which provide significant opportunities for detailed studies of galaxies at small scales and high redshifts, offer a novel perspective on small-scale cosmic phenomena and constrain the characteristics of PMFs. In this study, we compile 14 SMD data points at redshifts $z > 6$ and derive stringent constraints on the parameters of PMFs, which include the amplitude of the magnetic field at a characteristic scale of $\lambda=1\,{\rm Mpc}$, denoted as $B_0$, and the spectral index of the magnetic field power spectrum, $n_{\rm B}$. At 95\% confidence level, we establish upper limits of $B_0 < 4.44$ nG and $n_{\rm B} < -2.24$, along with a star formation efficiency of approximately $f_*^0 \sim 0.1$. If we fix $n_{\rm B}$ at specific values, such as $-2.85$, $-2.9$, and $-2.95$, the 95\% upper limits for the amplitude of the magnetic field can be constrained to 1.33 nG, 2.21 nG, and 3.90 nG, respectively. Finally, we attempt to interpret recent early observations provided by James Webb Space Telescope (JWST) using the theory of PMFs, and find that by selecting appropriate PMF parameters, it is possible to explain these results without significantly increasing the star formation efficiency.

Understanding plasma dynamics and nonthermal particle acceleration in 3D magnetic reconnection has been a long-standing challenge. In this paper, we explore these problems by performing large-scale fully kinetic simulations of multi-xline plasmoid reconnection with various parameters in both the weak and strong guide field regimes. In each regime, we have identified its unique 3D dynamics that leads to field-line chaos and efficient acceleration, and we have achieved nonthermal acceleration of both electrons and protons into power-law spectra. The spectral indices agree well with a simple Fermi acceleration theory that includes guide field dependence. In the low-guide-field regime, the flux-rope kink instability governs the 3D dynamics for efficient acceleration. The weak dependence of the spectra on the ion-to-electron mass ratio and $\beta$ ($\ll1$) implies that the particles are sufficiently magnetized for Fermi acceleration in our simulations. While both electrons and protons are injected at reconnection exhausts, protons are primarily injected by perpendicular electric fields through Fermi reflections and electrons are injected by a combination of perpendicular and parallel electric fields. The magnetic power spectra agree with in-situ magnetotail observations, and the spectral index may reflect a reconnection-driven size distribution of plasmoids instead of Goldreich-Sridhar vortex cascade. As the guide field becomes stronger, the oblique flux ropes of large sizes capture the main 3D dynamics for efficient acceleration. Intriguingly, the oblique flux ropes can also run into flux-rope kink instability to drive extra 3D dynamics. This work has broad implications for 3D reconnection dynamics and particle acceleration in heliophysics and astrophysics., Comment: 15 pages, 11 figures