Search

Hybrid simulations are carried out to study the instabilities driven by the outer heliosheath pickup ions of ring-beam velocity distributions. Following the linear instability analysis in the companion paper, the electromagnetic instabilities at various propagation directions (0° ≤ θ ≤ 180°) are studied for pickup ion ring-beam distributions of different pickup angles (0° ≤ α ≤ 90°). Previous studies of the ring-beam pickup ion distributions have been limited to the parallel- and antiparallel-propagating modes (θ = 0° and 180°). In this paper, several one-dimensional hybrid simulations are performed to investigate the most unstable modes of different propagation angles predicted by the linear instability analysis. Moreover, two-dimensional hybrid simulations are conducted to examine the simultaneous excitation of the parallel, antiparallel, and oblique unstable modes and their role in the scattering of the outer heliosheath pickup ions. Our recent one-dimensional hybrid and particle-in-cell simulations with simulation domains along the background magnetic field demonstrated that at small pickup angles, the outer heliosheath pickup ions experience an incomplete scattering limited to the hemisphere of positive parallel velocities with respect to the background magnetic field. In contrast, the present two-dimensional hybrid simulation results show that when both parallel/antiparallel and oblique unstable modes are acting together, the pickup ions can be scattered into the hemisphere of negative parallel velocities and approach a near-isotropic velocity distribution. This is in conflict with the spatial retention scenario of the energetic neutral atom ribbon observed by the Interstellar Boundary Explorer, which requires pickup ions at small pickup angles to stream away from the ribbon directions. [ABSTRACT FROM AUTHOR]

The energetic neutral atom ribbon observed by the Interstellar Boundary Explorer spacecraft is believed to originate from the pickup ions in the outer heliosheath. The outer heliosheath pickup ions generally have a ring-beam velocity distribution at a certain pickup angle, α, the angle at which these ions are picked up by the interstellar magnetic field. The pickup ion ring-beam distributions can drive unstable waves of different propagation angles with respect to the background interstellar magnetic field, θ. Previous studies of the outer heliosheath pickup ion dynamics were mainly focused on ring-like pickup ion distributions with α ≈ 90° and/or the parallel- and antiparallel-propagating unstable waves (θ = 0° and 180°). The present study carries out linear kinetic instability analysis to investigate both the parallel and oblique unstable modes (0° ≤ θ ≤ 180°) driven by ring-beam pickup ion distributions of different pickup angles between 0° and 90°. Our linear instability analysis reveals that ring-beam pickup ion distributions can excite oblique mirror waves as well as parallel/quasi-parallel and oblique right- and left-helicity waves. The maximum growth rate among all the instabilities belongs to the parallel-propagating left-helicity waves at most pickup angles. Furthermore, the evolution of the unstable mirror waves by varying pickup angle indicates that as the pickup angle increases, the maximum growth rate of the mirror modes increases, while its propagation angle decreases. [ABSTRACT FROM AUTHOR]

The electromagnetic instabilities driven by the pickup ions in the outer heliosheath play a crucial role in interpreting the energetic neutral atom (ENA) ribbon observed by the Interstellar Boundary EXplorer (IBEX). Previous studies on the stability of the outer heliosheath pickup ions focused on ring-like velocity distributions with pickup angles equal or close to 90°. In this study, we investigate the instabilities driven by a ring-beam velocity distribution of the pickup ions in the outer heliosheath using linear instability analysis. The modes propagating along (both parallel and antiparallel to) the background magnetic field are examined, and the emphasis is on how the instabilities change with the pickup angle, α, varying from zero to 90°. The parallel thermal spread of the pickup ions is chosen to make the plasma parameters lie in the Alfvén cyclotron 'stability gap' to exclude the well-studied Alfvén cyclotron instability. Our linear instability analysis reveals that the unstable modes all have left helicity but may propagate parallel or antiparallel to the background magnetic field. At small α, the parallel-propagating modes are unstable in two separate wavenumber (and frequency) ranges. The two unstable ranges merge into one when α increases to 82°. In addition, the antiparallel-propagating modes become unstable when α is above 5°, but the parallel-propagating modes have larger growth rates until α becomes larger than αc = 82.8°. The growth rates of the parallel-propagating modes quickly drop below zero when α goes beyond αc, while the antiparallel-propagating modes remain unstable until α approaches 90°. [ABSTRACT FROM AUTHOR]

Two‐dimensional hybrid particle‐in‐cell (PIC) simulations are carried out on a constant L‐shell (or drift shell) surface of the dipole magnetic field to investigate the generation process of near‐equatorial fast magnetosonic waves (a.k.a equatorial noise; MSWs hereafter) in the inner magnetosphere. The simulation domain on a constant L‐shell surface adopted here allows wave propagation and growth in the azimuthal direction (as well as along the field line) and is motivated by the observations that MSWs propagate preferentially in the azimuthal direction in the source region. Furthermore, the equatorial ring‐like proton distribution used to drive MSWs in the present study is (realistically) weakly anisotropic. Consequently, the ring‐like velocity distribution projected along the field line by Liouville's theorem extends to rather high latitude, and linear instability analysis using the local plasma conditions predicts substantial MSW growth up to ±27° latitude. In the simulations, however, the MSW intensity maximizes near the equator and decreases quasi‐exponentially with latitude. Further analysis reveals that the stronger equatorward refraction at higher latitude due to the larger gradient of the dipole magnetic field strength prevents off‐equatorial MSWs from growing continuously, whereas MSWs of equatorial origin experience little refraction and can fully grow. Furthermore, the simulated MSWs exhibit a rather complex wave field structure varying with latitude, and the scattering of energetic ring‐like protons in response to MSW excitation occurs faster than the bounce period of those protons so that they do not necessarily follow Liouville's theorem during MSW excitation. Key Points: The 2‐D hybrid PIC simulations are carried out on a constant L‐shell surface to simulate the excitation of MSWs for the first timeDespite the extended unstable region in latitude, MSWs do not grow well if they get latitudinally out of resonanceScattering of ring‐like protons during MSW excitation is local so that those protons do not necessarily follow Liouville's theorem [ABSTRACT FROM AUTHOR]

Rising‐tone chorus waves have already been successfully produced in a mirror magnetic field with the use of one‐ and two‐dimensional particle‐in‐cell (PIC) simulations. However, in reality, the background magnetic field in the inner Earth's magnetosphere is a dipole magnetic field, unlike symmetric mirror fields. In this paper, with the two‐dimensional (2‐D) general curvilinear PIC (gcPIC) code, we investigate the generation of rising‐tone chorus waves in the dipole magnetic field configuration. The plasma consists of three components: immobile ions, cold background, and hot electrons. In order to save computational resource, the topology of the magnetic field is roughly equal to that at L = 0.6 RE, although the plasma parameters corresponding to those at L = 6 RE (RE is the Earth's radius) are used. Whistler mode waves are first excited around the magnetic equator by the hot electrons with a temperature anisotropy. The excited whistler mode waves propagate almost parallel and antiparallel to the background magnetic field in their source region, which is limited at ∣λ ∣  ≤ 3° (where λ is the magnetic latitude). When the waves leave from the source region and propagate toward high latitudes, both their amplitude and wave normal angle become larger. However, the group velocity of the waves is directed toward high latitudes almost along the magnetic field. During such a process, the waves have a frequency chirping, as shown by a rising tone in the frequency‐time spectrogram. To our best knowledge, it is for the first time that rising‐tone chorus are generated in a dipole magnetic field with a PIC simulation. Key Points: The rising‐tone chorus waves are generated in a dipole magnetic field with a 2‐D gcPIC simulation for the first timeWhistler waves are excited around the magnetic equator and propagate almost along the background magnetic fieldWhen whistler waves propagate toward high latitudes, their wave normal angles increase with the growth of the latitude [ABSTRACT FROM AUTHOR]

The whistler anisotropy instability is studied in a magnetized, homogeneous, collisionless plasma model. The electrons (denoted by subscript e) are represented initially with a single bi-Maxwellian velocity distribution with a temperature anisotropy T⊥e/T∥e>1, where ⊥ and ∥ denote directions perpendicular and parallel to the background magnetic field Bo, respectively. Kinetic linear dispersion theory predicts that, if the ratio of the electron plasma frequency ωe to the electron cyclotron frequency Ωe is greater than unity and β∥e≥0.025, the maximum growth rate of this instability is at parallel propagation, where the fluctuating fields are strictly electromagnetic. At smaller values of β∥e, however, the maximum growth rate shifts to propagation oblique to Bo and the fluctuating electric fields become predominantly electrostatic. Linear theory and two-dimensional particle-in-cell simulations are used to examine the consequences of this transition. Three simulations are carried out, with initial β∥e=0.10, 0.03, and 0.01. The fluctuating fields of the β∥e=0.10 run are predominantly electromagnetic, with nonlinear consequences similar to those of simulations already described in the literature. In contrast, the growth of fluctuations at oblique propagation in the low electron β runs leads to a significant δE∥, which heats the electrons leading to the formation of a substantial suprathermal component in the electron parallel velocity distribution. [ABSTRACT FROM AUTHOR]

This paper describes a theory of the velocity space scattering of charged particles in a static magnetic field composed of a uniform background field and a sum of transverse, circularly polarized, magnetic waves. When that sum has many terms the autocorrelation time required for particle orbits to become effectively randomized is small compared with the time required for the particle velocity distribution to change significantly. In this regime the deterministic equations of motion can be transformed into stochastic differential equations of motion. The resulting stochastic velocity space scattering is described, in part, by a pitch angle diffusion rate that is a function of initial pitch angle and properties of the wave spectrum. Numerical solutions of the deterministic equations of motion agree with the theory at all pitch angles, for wave energy densities up to and above the energy density of the uniform field, and for different wave spectral shapes. [ABSTRACT FROM AUTHOR]