skip to main content


Title: Particle-in-cell simulations of Alfvén wave parametric decay in a low-beta plasma

We study the parametric decay instability of parallel-propagating Alfvén waves in a low-beta plasma using one-dimensional fully kinetic simulations. We focus for the first time on the conversion of the energy stored in the initial Alfvén wave into particle internal energy, and on its partition between particle species. We show that compressible fluctuations generated by the decay of the pump wave into a secondary ion-acoustic mode and a reflected Alfvén wave contribute to the gain of internal energy via two distinct mechanisms. First, the ion-acoustic mode leads nonlinearly to proton trapping and proton phase-space mixing, in agreement with previous work based on hybrid simulations. Second, during the nonlinear stage, a compressible front of the fast type develops at the steepened edge of the backward Alfvén wave leading to a field-aligned proton beam propagating backwards at the Alfvén speed. We find that parametric decay heats preferentially protons, which gain approximately 50 % of the pump wave energy in the form of internal energy. However, we find that electrons are also energized and that they contribute to the total energy balance by gaining 10 % of the pump wave energy. By investigating energy partition and particle heating during parametric decay, our results contribute to the determination of the role of compressible and kinetic effects in wave-driven models of the solar wind.

 
more » « less
Award ID(s):
2141564
NSF-PAR ID:
10476844
Author(s) / Creator(s):
; ;
Publisher / Repository:
Journal of Plasma Physics
Date Published:
Journal Name:
Journal of Plasma Physics
Volume:
89
Issue:
2
ISSN:
0022-3778
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We present the results of numerical studies of the whistler wave parametric decay instability in the system with the suppressed Landau damping of ion acoustic waves (IAWs) based on the self-consistent Darwin particle-in-cell (PIC) model. It has been demonstrated that a monochromatic whistler wave launched along the background magnetic field couples to a counter-propagating whistler mode and co-propagating ion acoustic mode. The coupling of the electromagnetic mode to the electrostatic mode is guided by a ponderomotive force that forms spatio-temporal beat patterns in the longitudinal electric field generated by the counter-propagating whistler and the pump whistler wave. The threshold amplitude for the instability is determined to be δB w / B 0 = 0.028 and agrees with a prediction for the ion decay instability: δB w / B 0 = 0.042 based on the linear kinetic damping rates, and δB w / B 0 = 0.030 based on the simulation derived damping rates. Increasing the amplitude of the pump whistler wave, the secondary and tertiary decay thresholds are reached, and cascading parametric decay from the daughter whistler modes is observed. At the largest amplitude ( δB w / B 0 ∼ 0.1) the primary IAW evolves into a short-lived and highly nonlinear structure. The observed dependence of the IAW growth rate on the pump wave amplitude agrees with the expected trend; however, quantitatively, the growth rate of the IAW is larger than expected from theoretical predictions. We discuss the relevant space regimes where the instability could be observed and extensions to the parametric coupling of whistler waves with the electron acoustic wave (EAW). 
    more » « less
  2. With the support of hybrid-kinetic simulations and analytic theory, we describe the nonlinear behaviour of long-wavelength non-propagating (NP) modes and fast magnetosonic waves in high- $\beta$ collisionless plasmas, with particular attention to their excitation of and reaction to kinetic micro-instabilities. The perpendicularly pressure balanced polarization of NP modes produces an excess of perpendicular pressure over parallel pressure in regions where the plasma $\beta$ is increased. For mode amplitudes $|\delta B/B_0| \gtrsim 0.3$ , this excess excites the mirror instability. Particle scattering off these micro-scale mirrors frustrates the nonlinear saturation of transit-time damping, ensuring that large-amplitude NP modes continue their decay to small amplitudes. At asymptotically large wavelengths, we predict that the mirror-induced scattering will be large enough to interrupt transit-time damping entirely, isotropizing the pressure perturbations and morphing the collisionless NP mode into the magnetohydrodynamic (MHD) entropy mode. In fast waves, a fluctuating pressure anisotropy drives both mirror and firehose instabilities when the wave amplitude satisfies $|\delta B/B_0| \gtrsim 2\beta ^{-1}$ . The induced particle scattering leads to delayed shock formation and MHD-like wave dynamics. Taken alongside prior work on self-interrupting Alfvén waves and self-sustaining ion-acoustic waves, our results establish a foundation for new theories of electromagnetic turbulence in low-collisionality, high- $\beta$ plasmas such as the intracluster medium, radiatively inefficient accretion flows and the near-Earth solar wind. 
    more » « less
  3. Abstract Parker Solar Probe (PSP) observed predominately Alfvénic fluctuations in the solar wind near the Sun where the magnetic field tends to be radially aligned. In this paper, two magnetic-field-aligned solar wind flow intervals during PSP’s first two orbits are analyzed. Observations of these intervals indicate strong signatures of parallel/antiparallel-propagating waves. We utilize multiple analysis techniques to extract the properties of the observed waves in both magnetohydrodynamic (MHD) and kinetic scales. At the MHD scale, outward-propagating Alfvén waves dominate both intervals, and outward-propagating fast magnetosonic waves present the second-largest contribution in the spectral energy density. At kinetic scales, we identify the circularly polarized plasma waves propagating near the proton gyrofrequency in both intervals. However, the sense of magnetic polarization in the spacecraft frame is observed to be opposite in the two intervals, although they both possess a sunward background magnetic field. The ion-scale plasma wave observed in the first interval can be either an inward-propagating ion cyclotron wave (ICW) or an outward-propagating fast-mode/whistler wave in the plasma frame, while in the second interval it can be explained as an outward ICW or inward fast-mode/whistler wave. The identification of the exact kinetic wave mode is more difficult to confirm owing to the limited plasma data resolution. The presence of ion-scale waves near the Sun suggests that ion cyclotron resonance may be one of the ubiquitous kinetic physical processes associated with small-scale magnetic fluctuations and kinetic instabilities in the inner heliosphere. 
    more » « less
  4. Abstract

    We investigate the local proton energization at magnetic discontinuities/intermittent structures and the corresponding kinetic signatures in velocity phase space in Alfvénic (high cross helicity) and non-Alfvénic (low cross helicity) wind streams observed by Parker Solar Probe. By means of the partial variance of increments method, we find that the hottest proton populations are localized around compressible, coherent magnetic structures in both types of wind. Analysis of parallel and perpendicular temperature distributions suggest that the Alfvénic wind undergoes preferential enhancements ofTat such structures, whereas the non-Alfvénic wind experiences preferentialTenhancements. Although proton beams are present in both types of wind, the proton velocity distribution function displays distinct features. Hot beams, i.e., beams with beam-to-core perpendicular temperatureT⊥,b/T⊥,cup to three times larger than the total distribution anisotropy, are found in the non-Alfvénic wind, whereas colder beams are in the Alfvénic wind. Our data analysis is complemented by 2.5D hybrid simulations in different geometrical setups, which support the idea that proton beams in Alfvénic and non-Alfvénic wind have different kinetic properties and different origins. The development of a perpendicular nonlinear cascade, favored in balanced turbulence, allows a preferential relative enhancement of the perpendicular plasma temperature and the formation of hot beams. Cold field-aligned beams are instead favored by Alfvén wave steepening. Non-Maxwellian distribution functions are found near discontinuities and intermittent structures, pointing to the fact that the nonlinear formation of small-scale structures is intrinsically related to the development of highly nonthermal features in collisionless plasmas. Our results contribute to understanding the role of different coherent structures in proton energization and their implication in collisionless energy dissipation processes in space plasmas.

     
    more » « less
  5. ABSTRACT

    Cosmic rays (CRs) are thought to be an important feedback mechanism in star-forming galaxies. They can provide an important source of pressure support and possibly drive outflows. We perform multidimensional CR magnetohydrodynamic simulations including transport by streaming and diffusion to investigate wind launching from an initially hydrostatic atmosphere by CRs. We estimate a characteristic Eddington limit on the CR flux for which the CR force exceeds gravity and compare it to simulated systems. Scaling our results to conditions in star-forming galaxies, we find that CRs are likely to contribute to driving outflows for a broad range of star formation environments. We quantify the momentum and energy transfer between CRs and gas, along with the associated mass outflow rates under different assumptions about the relative importance of streaming and diffusion for transport. In simulations with streaming, we observe the growth and saturation of the CR acoustic instability, but the CRs and gas remain well coupled, with CR momentum transferred efficiently to the gas even when this instability is present. Higher CR fluxes transfer more energy to the gas and drive stronger outflows. When streaming is present, most of the transferred energy takes the form of Alfvén wave heating of the gas, raising its pressure and internal energy, with a lower fractional contribution to the kinetic energy of the outflow. We also consider runs with radiative cooling, which modifies gas temperature and pressure profiles but does not seem to have a large impact on the mass outflow for super-Eddington CR fluxes.

     
    more » « less