On the Nature and Origin of Bipolar Electrostatic Structures in the Earth's Bow Shock

Abstract

We present a statistical analysis of large-amplitude bipolar electrostatic structures measured by Magnetospheric Multiscale spacecraft in the Earth's bow shock. The analysis is based on 371 large-amplitude bipolar structures collected in nine supercritical quasi-perpendicular Earth's bow shock crossings. We find that 361 of the bipolar structures have negative electrostatic potentials, and only 10 structures (< 3%) have positive potentials. The bipolar structures with negative potentials are interpreted in terms of ion phase space holes produced by ion streaming instabilities, particularly the two-stream instability between incoming and reflected ions. We obtain an upper estimate for the amplitudes of the ion phase space holes that is in agreement with the measurements. The bipolar structures with positive potentials could be electron phase space holes produced by electron two-stream instabilities. We argue that the negligible number of electron phase space holes among large-amplitude bipolar structures is due to the electron hole transverse instability, the criterion for which is highly restrictive at ω_pe/ω_ce ≫ 1, a parameter range typical of collisionless shocks in the heliosphere and various astrophysical environments. Our analysis indicates that the original mechanism of electron surfing acceleration involving electron phase space holes is not likely to be efficient in realistic collisionless shocks.

1. Introduction

The Earth's bow shock is a natural laboratory for probing plasma processes in supercritical collisionless shocks, because the Alfvén Mach number of the solar wind flow typically exceeds the second critical value [1]. Of vital importance in the physics of collisionless shocks are the origin and consequences of various wave activities produced in the shock transition region under realistic background plasma parameters [2, 3]. In particular, there is currently a lack of a thorough understanding of the origin and consequences of electrostatic fluctuations measured in the Earth's bow shock [4–6], which are likely present in astrophysical shocks as well [7–9]. Numerical simulations have demonstrated that electrostatic fluctuations and coherent electrostatic structures produced by various instabilities in a shock transition region can result in efficient electron acceleration (surfing and stochastic surfing mechanisms) at high Mach numbers typical of various astrophysical shocks [8–12] and at low Mach numbers typical of the Earth's bow shock [13]. However, the efficiency of those acceleration processes in realistic collisionless shocks, where typical background plasma parameters are hardly achievable in modern numerical simulations, remains poorly understood. Among the critical background plasma parameters is ω_pe/ω_ce, the ratio of the electron plasma frequency to the electron cyclotron frequency [8], which is typically large in space and astrophysical plasmas [14] but much smaller in numerical simulations of collisionless shocks [3, 15]. The analysis of spacecraft measurements in the Earth's bow shock characterized by ω_pe/ω_ce ~ 100 should advance our understanding of the origin and consequences of electrostatic fluctuations in collisionless shocks under realistic background plasma parameters.

Spacecraft measurements in the Earth's bow shock showed that electric field fluctuations above a few hundred hertz tend to be electrostatic [5, 16, 17]. Analyses of waveform measurements demonstrated that the electrostatic wave activity is produced by ion-acoustic waves [4, 18–21], bipolar electrostatic structures [22–27], electron cyclotron harmonics [6, 28], and low-hybrid waves [29]. The ion-acoustic waves have electric field oriented generally oblique to the magnetic field, wavelengths on the order of a few tens of Debye lengths, and amplitudes of up to a few hundred mV/m [4, 20, 21, 30]. The bipolar structures have similar properties, except that their spatial scales are of only a few Debye lengths [23, 25–27].

The bipolar structures originally measured aboard the Wind spacecraft were interpreted in terms of electron phase space holes [23, 24], which are coherent structures with positive electrostatic potentials appearing in a nonlinear stage of various electron streaming instabilities [31–33]. However, detailed analyses of bipolar structures measured in several crossings of the Earth's bow shock by Cluster spacecraft [25] and Magnetospheric Multiscale spacecraft (MMS) [26, 27] showed that the bipolar structures had negative electrostatic potentials, which is inconsistent with the interpretation in terms of electron phase space holes. Wang et al. [27] have recently considered more than one hundred bipolar structures measured in a particular crossing of the Earth's bow shock and provided strong arguments that the bipolar structures are ion phase space holes produced by the two-stream instability between incoming and reflected ions. The ion two-stream instability was hypothesized to be a source of electrostatic fluctuations in the Earth's bow shock early on by Formisano and Torbert [18], Akimoto and Winske [34], and Fuselier et al. [35], but there have been no strong arguments in favor of that hypothesis until the recent MMS measurements [27].

The analysis of bipolar structures based on MMS measurements is currently restricted to a few of the Earth's bow shock crossings [26, 27]. In this paper we present results of a statistical analysis of bipolar structures measured by MMS in nine crossings of the Earth's bow shock. The analysis provides valuable information on instabilities operating in the Earth's bow shock and reveals the presence, albeit in negligible numbers, of bipolar structures with positive potentials, which could be electron phase space holes. The rarity of electron phase space holes among large-amplitude bipolar structures is most likely due to the electron hole transverse instability, the criterion for which is strongly dependent on ω_pe/ω_ce [36, 37]. We discuss the implications of our results for the surfing acceleration mechanism, which suggested that electron phase space holes might be involved in efficient electron acceleration in high-Mach-number collisionless shocks [10, 11].

2. MMS Observations

To collect a statistically representative dataset of bipolar electrostatic structures, we considered nine crossings of the Earth's bow shock by MMS. The selection criteria for the bow shock crossings were the presence of continuous measurements in burst mode in the shock transition region, a quasi-perpendicular configuration, and a supercritical regime of the shock, that is, an Alfvén Mach number M_A greater than 3. The selected Earth's bow shock crossings are listed in Table 1. We note that shocks #6 and #8 were previously considered by Vasko et al. [26] and Wang et al. [27], but the analysis of Vasko et al. [26] was restricted to a few bipolar structures, while Wang et al. [27] excluded from consideration a far-downstream region of the shock. We used measurements taken by the following instruments aboard the MMS: DC-coupled magnetic field at 128 S/s (samples per second) provided by Digital and Analogue Fluxgate Magnetometers [38], AC-coupled electric fields at 8,192 S/s provided by Axial Double Probe [39], and Spin-Plane Double Probe [40], and electron moments at 0.03 s cadence and ion moments at 0.15 s cadence provided by the Fast Plasma Investigation instrument [41]. The electric field was measured by four voltage-sensitive spherical probes on 60 m antennas in the spacecraft spin plane (almost in the ecliptic plane) and two probes on roughly 15 m axial antennas along the spin axis (almost perpendicular to the ecliptic plane). The voltages of the opposing probes measured with respect to the spacecraft are used to compute the electric field and estimate the velocity of propagation and other parameters of bipolar electrostatic structures [26, 42]. Because the spatial scales of bipolar structures are generally comparable to the antenna lengths, corresponding correction factors are included to accurately estimate the electric fields (see Vasko et al. [26] for methodological details).

Figure 1 presents an overview of the selected Earth's bow shock crossings. The panels display the magnetic field magnitude and parallel electron temperature measured aboard MMS1. The perpendicular electron temperature is almost identical to the parallel temperature (not shown here), because in supercritical shocks the macroscopic electron heating is essentially isotropic [43, 44]. The other MMS located within a few tens of kilometers of MMS1 yield essentially identical overviews of the shocks. The dates and time intervals corresponding to the shocks are shown in the panels, while Table 1 presents estimates of the Alfvén Mach number M_A and angle θ_Bn between the upstream magnetic field and the shock normal. The selected shocks are supercritical and quasi-perpendicular, with 3 ≲ M_A ≲ 13 and . The shock normals were assumed to be directed toward the upstream region and estimated using the procedure of Vinas and Scudder [45], which is based on considering the Rankine-Hugoniot conditions. Table 1 also presents estimates of the ion and electron beta parameters β_i and β_e, the ion temperature T_i, and the electron-to-ion temperature ratio T_e/T_i in the upstream region. Because MMS measurements of the ion temperature in the solar wind are not accurate, we used ion temperature measurements by the Wind spacecraft¹ located at the L1 Earth-Sun Lagrange point, which is about two hundred of Earth's radius upstream of the Earth's bow shock. An ion temperature estimate T_i shown in Table 1 for a particular shock represents an averaged value over 20 min around a time instant of the shock crossing aboard MMS1. Ion temperature variations around an averaged value over these 20 min are also shown in Table 1 and demonstrate that the presented estimates of T_i and hence β_i are rather reliable.

Figure 1

Figure 2 presents the electric field measured by MMS1 at 8,192 S/s in shock #8. The electric field magnitude shown in Figure 2A demonstrates that the electric field fluctuations reach amplitudes of up to a few hundred mV/m. Figure 2B displays the electric field in a local magnetic field-aligned coordinate system measured over the 35 ms interval highlighted in Figure 2A. We can see that some of the intense electric field fluctuations correspond to bipolar electrostatic structures. The bipolar structures are noticeable in the electric field components both parallel and perpendicular to the local magnetic field, thereby indicating that the electric fields of the bipolar structures are oriented oblique to the magnetic field, typically at a few tens of degrees [26]. We looked through the electric fields measured aboard the four MMS in the selected Earth's bow shock crossings and collected a dataset of 371 bipolar structures with electric field amplitudes exceeding 50 mV/m. The selection was restricted to large-amplitude bipolar structures to make a careful analysis of these structures feasible. The selection procedure, based on visual inspection, is subjective and probably not exhaustive, but it is adequate enough for collecting a representative dataset of bipolar structures. In the resulting dataset the ratio between the maximum and minimum values of any bipolar electric field does not exceed 2, and for more than 95% of the structures the ratio was below 1.5, so that the selected structures may indeed be referred to as bipolar structures. In addition to bipolar structures, large-amplitude electric field fluctuations in the Earth's bow shock are produced by electrostatic wave-packets and tripolar structures (Figure 2B). The electrostatic wave-packets and the bipolar and tripolar structures are highly likely to be from a common origin (produced by the same instability or instabilities), but in this study we concentrate purely on bipolar structures. The occurrence of the selected bipolar structures in shock transition regions is demonstrated in Figure 1, while Table 1 gives the number of bipolar structures selected in each shock. The number of bipolar structures can be rather different in various shocks, but none of the parameters M_A, θ_Bn, β_e, β_i, and T_e/T_i in Table 1 can explain the observed variation in the number of bipolar structures from one shock to another. In what follows, we explain the methodology and present results of the statistical analysis of the selected large-amplitude bipolar structures.

Figure 2

Figure 3 presents the analysis of a particular bipolar structure measured in shock #8. We consider voltage signals induced on voltage-sensitive probes by the electric field of the bipolar structure. Figures 3a,b plot the voltage signals V₁ vs. −V₂ and V₃ vs. −V₄ measured by two pairs of opposing probes on 60 m antennas in the spacecraft spin plane, while Figure 3c plots the voltage signals V₅ vs. −V₆ measured by two opposing probes on roughly 15 m antennas parallel to the spin axis. The noticeable correlations between the voltage signals of the opposing probes reflect the high quality of the electric field measurements [42]. Figure 3d presents components E_ij of the electric field along the antenna directions computed using the voltage signals of the opposing probes, E_ij ∝ (V_i − V_j)/2l_ij, with correction factors for the antenna frequency response and finite antenna length (see Vasko et al. [26] for the methodology), where l₁₂ = l₃₄ = 60 m and l₅₆ = 14.6 m are antenna lengths. The bipolar profiles of the electric field components E_ij indicate that the bipolar structure is essentially one-dimensional, while Figure 3e presents the dominant electric field E_l determined using minimum variance analysis [46]. The profiles of the electric field components E_ij and the time delays between voltage signals measured by the opposing probes allow us to determine the sign of the electrostatic potential of the bipolar structure. The bipolar structure propagates from probe 1 to probe 2, while the electric field E₁₂ is first negative and then positive. This implies that the electric field E₁₂ has a convergent spatial configuration and, hence, the bipolar structure has to have a negative electrostatic potential. Similar analysis based on the other opposing probes consistently shows that the bipolar structure has a negative electrostatic potential.

Figure 3

The time delays Δt_ij between the voltage signals of the opposing probes are indicated in Figures 3a–c. We use the time delays to estimate the velocity V_s of the bipolar structure in the spacecraft frame, , and the direction k of propagation of the bipolar structure, k_ij = −V_sΔt_ij/l_ij [the minus sign was missed in Vasko et al. [26] but did not affect their results]. We have found that the bipolar structure propagates with velocity V_s ≈ 50 km/s along vector k ≈ (0.43, −0.21, 0.88), which is within 8° of the electric field direction. The small angle between the electric field and the propagation direction is a strong indication that the bipolar structure is approximately a one-dimensional electrostatic structure. The estimated velocity allows us to translate temporal profiles into spatial profiles. Figure 3e shows that the spatial scale of the bipolar structure, determined as half the distance between minimum and maximum values of the electric field E_l, is l ≈ 17 m. Figure 3e presents the electrostatic potential of the bipolar structure, computed as Φ = ∫E_lV_sdt, and demonstrates that the amplitude of the electrostatic potential is Φ₀ ≈ −3.8 V. In terms of the local Debye length λ_D and electron temperature T_e, we have l ≈ 2λ_D and eΦ₀ ≈ −0.15 T_e, where e is the electron charge. This bipolar structure is similar to previously reported bipolar structures with negative electrostatic potentials in the Earth's bow shock [25–27].

The extensive dataset of bipolar structures that we collected allows us to detect bipolar structures with positive electrostatic potentials not reported previously. Figure 4 presents the analysis of one of these bipolar structures measured in shock #4. The noticeable correlations between the voltage signals of the opposing probes in Figures 4a–c indicate the high quality of the electric field measurements. Similar bipolar profiles of the electric field components E_ij suggest that the bipolar structure should be approximately one-dimensional. The bipolar structure propagates from probe 2 to probe 1, while the electric field E₁₂ is first negative and then positive. This implies that the electric field E₁₂ has a divergent spatial configuration and, hence, the bipolar structure has to have a positive electrostatic potential. Analysis based on the other opposing probes consistently shows that the bipolar structure has a positive electrostatic potential. Using the time delays Δt_ij indicated in Figures 4a–c, we have found that the bipolar structure propagates with velocity V_s ≈ 45 km/s along direction k = (−0.37, 0.34, −0.86). The propagation direction is within 7° of the electric field direction, which is a strong indication that the bipolar structure is indeed a one-dimensional structure. Figure 4e shows that the spatial scale of the bipolar structure is l ≈ 25 m and the amplitude of the electrostatic potential is Φ₀ ≈ 7 V. In terms of the local Debye length and electron temperature, we have l ≈ 4λ_D and eΦ₀ ≈ 0.23 T_e.

Figure 4

The bipolar structures presented in Figures 3, 4 can be considered perfect in that all time delays between the voltage signals of the opposing probes could be well-determined. This is usually not the case for structures in the collected dataset, as demonstrated by a particular bipolar structure presented in Figure 5. Figures 5a–c show that the time delays can be well determined between the opposing probes in the spin plane, but not between the opposing probes on the spin axis. In addition, there are many bipolar structures in our dataset with time delays that are well determined only for one pair of the opposing probes. For bipolar structures such that at least one of the time delays is not well determined, we need to use a method of analysis different from the one described for bipolar structures in Figures 3, 4. Figures 3d, 4d shows that although the time delay between V₅ and −V₆ cannot be determined, the profile of the electric field E₅₆ is bipolar and similar to the profiles of E₁₂ and E₃₄. In fact, the electric field components E_ij always have similar bipolar profiles for the selected bipolar structures, which is a strong indication that the bipolar structures are approximately one-dimensional structures. Because the bipolar structures are electrostatic, we assume the propagation direction k to be parallel to the electric field direction determined by minimum variance analysis. We note that minimum variance analysis yields the electric field direction and, hence, the propagation direction k up to a sign chosen such that the velocity V_s is positive, and so the bipolar structure propagates parallel to the electric field direction. The velocity of propagation is estimated as V_s = −k_ijl_ij/Δt_ij, where ij is a pair of opposing probes with the best determined time delay Δt_ij. Whether a bipolar structure has a negative or a positive electrostatic potential is also determined using that pair of probes. In cases where time delays are well determined on two antennas, those antennas always provided consistent results on the sign of the electrostatic potential of a bipolar structure and similar estimates of the velocity V_s. In particular, the bipolar structure presented in Figure 5 has a positive electrostatic potential as inferred from both pairs of opposing probes in the spin plane and propagates with velocity V_s ≈ 85 km/s. Figure 5e shows that the bipolar structure has spatial scale l ≈ 20 m and amplitude of the electrostatic potential Φ₀ ≈ 3 V. In terms of the local Debye length and electron temperature, we have l ≈ 2λ_D and eΦ₀ ≈ 0.12 T_e.

Figure 5

We have estimated the parameters of the selected bipolar structures using the interferometry method based on the antenna with the best determined time delay between voltage signals of the opposing probes. We have found that 361 of the bipolar structures have negative electrostatic potentials, while the remaining 10 have positive electrostatic potentials. In Figure 1, the occurrence of the bipolar structures is indicated along with the types of the bipolar structures, while Table 1 presents the numbers of bipolar structures with negative and positive potentials observed in each shock. The bipolar structures with positive potentials are observed in shocks #2, #3, #4, and #8; four of these structures are observed around shock ramps, while the other six occur in far-downstream regions. None of parameters M_A, θ_Bn, β_{e, i}, and T_e/T_i shown in Table 1 could be identified as being critical for the appearance of bipolar structures with positive electrostatic potentials.

Figure 6 presents a summary of the parameters of the bipolar structures. Figure 6a shows that the bipolar structures, whether they have positive and negative potentials, propagate with velocities from a few tens of km/s up to a few hundred km/s in the spacecraft frame, which is much smaller than the typical electron thermal velocity in the Earth's bow shock. Typical amplitudes |Φ₀| of the bipolar structures are a few volts, though they can be up to 30 V. Figure 6b shows that in terms of local electron temperature, the bipolar structures typically have amplitudes of |eΦ₀| ≲ 0.1 T_e, though these can be as large as 0.5 T_e. Typical spatial scales of the bipolar structures are a few λ_D, though they can be up to about 10λ_D. There is a noticeable positive correlation between |eΦ₀|/T_e and l/λ_D, whose physical nature will be clarified in the next section. Figures 6a,b demonstrate that velocities, spatial scales and absolute amplitudes of the bipolar structures with positive and negative potentials are basically similar.

Figure 6

3. Theoretical Interpretation

The bipolar electrostatic structures with negative potentials are most likely ion phase space holes, which are coherent structures formed in a nonlinear stage of various ion streaming instabilities [31, 48–50]. Interestingly, very early numerical simulations showed that ion streaming instabilities can potentially produce Debye-scale electrostatic structures with negative potentials in the Earth's bow shock [51, 52]. In the linear stage of an ion streaming instability, electrostatic fluctuations grow at the expense of the energy of resonant ions. In a nonlinear stage, the amplitude of the electrostatic fluctuations becomes sufficiently large to trap a fraction of resonant ions into potential wells of the electrostatic fluctuations, resulting in the formation of vortices in the ion phase space. The merging of these vortices leads to formation of electrostatic solitary waves called ion phase space holes [48, 49].

Although the details of ion velocity distribution functions, whose instability results in formation of ion phase space holes in the Earth's bow shock, are not known, we can estimate the growth rates of that instability. The instability saturates when the bounce frequency of ions trapped within electrostatic fluctuations becomes comparable to the initial increment of the instability [53]. Therefore, the initial increment γ should be of the order of the bounce frequency of ions trapped within the observed bipolar structures, . This condition can be written as where ω_pi is the ion plasma frequency. Adopting typical parameters of the bipolar structures (Figure 6), we find initial increments on the order of a fraction of the ion plasma frequency, γ ~ 0.1 ω_pi. In supercritical shocks, increments of that order can in principle be provided by the two-stream instability between incoming and reflected ions. The maximum growth rate of ion-acoustic waves driven by that instability is reached when both incoming and reflected ions are cold [34, 54]: , where α_ref is the fraction of incoming ions reflected to the upstream region and has a typical value of 0.1 in the Earth's bow shock [47]. Equation (1) can be written as which explains the positive correlation between the amplitude and spatial scale of the bipolar structures noticeable in Figure 6b. Because at finite temperatures of incoming and reflected ions the increment of the two-stream instability γ is smaller than γ_max [55], we can give the following upper estimate for the amplitudes of the bipolar structures:

Figure 6b demonstrates that the amplitudes of the bipolar structures are indeed below the upper estimate for typical fractions α_ref of reflected ions.

Direct identification of the ion two-stream instability using measurements of the ion velocity distribution function in the Earth's bow shock is expected to be complicated. The observed parameters of the bipolar structures indicate that the inverse bounce frequency of trapped ions is of the order of one millisecond. The relaxation of ion beams seeding the formation of ion phase space holes is expected to occur on a time scale of a few tens of inverse bounce periods [53], while ion velocity distribution functions aboard MMS are available at a 0.15 s cadence. Therefore, the measured ion velocity distribution functions are expected to correspond to a marginally stable relaxed state. The fast relaxation of the ion velocity distribution function explains the results of the stability analysis of Goodrich et al. [56], who demonstrated that measured ion velocity distribution functions are stable with respect to simultaneously measured ion-acoustic waves. An alternative approach to identifying the nature of an instability seeding the formation of ion phase space holes is based on analysis and comparison of the properties of these bipolar structures to predictions of a particular instability. Wang et al. [27] have suggested that approach and provided strong arguments in favor of the ion two-stream instability as the source of bipolar structures with negative potentials.

The analysis in section 2 showed that in addition to bipolar structures with negative potentials, there are bipolar structures with positive potentials, albeit in negligible numbers. The bipolar structures with positive potentials can be electron phase space holes [31, 32] produced by various electron streaming instabilities, which are expected to operate in the transition region of the Earth's bow shock [57, 58]. In terms of parameters, these bipolar structures are similar to slow electron phase space holes observed in reconnection current sheets [59, 60] and the plasma sheet boundary layer [61] in the Earth's magnetosphere. The bump-on-tail instability cannot be a source of these bipolar structures, because in that case the bipolar structures would propagate with velocities that are a fraction of the electron thermal velocity [57, 62]. Instabilities capable of producing slow electron phase space holes include the electron two-stream instability [58, 62, 63] and the Buneman instability [64, 65], though the latter is less likely to operate under conditions typical of the Earth's bow shock. Whichever instability potentially produces the bipolar structures with positive potentials, the initial increment should be on the order of the bounce frequency of electrons trapped within the bipolar electric field: where ω_pe is the electron plasma frequency. Adopting typical parameters of the bipolar structures (Figure 6), we find initial increments that are a fraction of the electron plasma frequency, γ ~ 0.1 ω_pe. Increments of this magnitude can in principle be provided by the electron two-stream instability [58]. Although the bipolar structures with positive potentials can be electron phase space holes, we note the alternative interpretation of these bipolar structures in terms of ion-acoustic solitons arising from ion-acoustic fluctuations due to electron and ion fluid nonlinearities [66, 67]. The detailed analysis of this alternative interpretation is left for a separate study.

The critical result of the statistical analysis in section 2 is that only 10 out of 371 bipolar structures, accounting for <3% of the structures, have positive potentials and can be interpreted in terms of electron phase space holes. Does the scarcity of bipolar structures with positive potentials imply that electron streaming instabilities rarely operate in the Earth's bow shock, or that the formation of these structures is hindered in conditions typical of the Earth's bow shock? We lean toward the second scenario. According to existing simulation and theoretical studies [36, 37], one-dimensional electron phase space holes can be stable with respect to the transverse instability only under the condition that the bounce frequency of trapped electrons should be smaller than the electron cyclotron frequency, ω_b ≲ ω_ce. This stability criterion can be written in the form of an upper estimate for electron hole amplitudes: With realistic background parameters, ω_pe/ω_ce ~ 100, an electron phase space hole with the spatial scale of a few Debye lengths can be stable with respect to the transverse instability, provided that , which implies Φ₀ ≲ 0.01 V, where we have assumed that by the order of magnitude T_e ~ 10 eV. Taking into account that a few Debye lengths is a few tens of meters, we find that electron phase space holes can be stable with respect to the transverse instability provided that the electric field amplitude is rather small, E ≲ 1 mV/m. Electron phase space holes with electric field amplitudes larger than about 1 mV/m should be rarely observed in the Earth's bow shock because of the transverse instability. Thus, the transverse instability may explain the observed scarcity of bipolar structures with positive potentials among large-amplitude bipolar structures in the Earth's bow shock.

4. Discussion and Conclusion

We have presented statistical analysis of large-amplitude bipolar electrostatic structures measured in the Earth's bow shock by Magnetospheric Multiscale spacecraft. We have found that over 97% of the bipolar structures in the Earth's bow shock have negative potentials. We have interpreted these bipolar structures in terms of ion phase space holes and obtained an upper estimate for the amplitude of these structures. The copious amount of ion phase space holes in the Earth's bow shock is a strong indication that ion streaming instabilities, and in particular the two-stream instability between incoming and reflected ions, drive electrostatic fluctuations in the Earth's bow shock. The statistical analysis has also shown that <3% of the bipolar structures have positive potentials and that these can in principle be interpreted in terms of electron phase space holes. The observed rarity of bipolar structures with positive potentials may be due to the transverse instability of electron phase space holes, the criterion for which is strongly dependent on ω_pe/ω_ce.

The results presented here imply that the efficiency of the electron surfing acceleration mechanism originally demonstrated in 1D Particle-In-Cell simulations of high-Mach-number shocks [10, 11] will be restricted. In the surfing acceleration mechanism, electron phase space holes appearing in a nonlinear stage of the Buneman instability operating in the shock transition region facilitate prolonged trapping of electrons around the shock ramp, resulting in efficient electron acceleration by the motional electric field. We note that the transverse instability of electron phase space holes is suppressed in one-dimensional simulations, because that instability is a multi-dimensional effect [36, 37]. In one-dimensional simulations, electron phase space holes are generally stable over an indefinitely long time at any ω_pe/ω_ce [62, 63]. In realistic three-dimensional shocks and under realistic background plasma parameters, ω_pe/ω_ce ≫ 1, the transverse instability of electron phase space holes will strongly reduce the efficiency of electron acceleration via the original surfing acceleration mechanism. Thus, our analysis of large-amplitude bipolar structures in the Earth's bow shock has demonstrated that the surfing acceleration mechanism involving electron phase space holes is highly unlikely to be efficient in realistic collisionless shocks. The stochastic surfing acceleration mechanism is a prospective mechanism, though, because it does not require the presence of coherent electrostatic structures [9, 12].

References

• 1.

  KennelCFEdmistonJPHadaT. A quarter century of collisionless shock research. In: StoneRGTsurutaniBT, editors. Collisionless Shocks in the Heliosphere. A Tutorial Review.Washington, DC: American Geophysical Union (1985). p. 1–36.

  • Google Scholar

• 2.

  PapadopoulosK. Microinstabilities and anomalous transport. In: StoneRGTsurutaniBT, editors. Collisionless Shocks in the Heliosphere. A Tutorial Review.Washington, DC: American Geophysical Union (1985). p. 59–90.

  • Google Scholar

• 3.

  KrasnoselskikhVBalikhinMWalkerSNSchwartzSSundkvistDLobzinVet al. The dynamic quasiperpendicular shock: cluster discoveries. Space Sci Rev. (2013) 178:535–98. 10.1007/s11214-013-9972-y

  • CrossRef
  • Google Scholar

• 4.

  BaleSDMozerFS. Measurement of large parallel and perpendicular electric fields on electron spatial scales in the terrestrial bow shock. Phys Rev Lett. (2007) 98:205001. 10.1103/PhysRevLett.98.205001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5.

  MozerFSSundkvistD. Electron demagnetization and heating in quasi-perpendicular shocks. J Geophys Res. (2013) 118:5415–20. 10.1002/jgra.50534

  • CrossRef
  • Google Scholar

• 6.

  WilsonLBSibeckDGBrenemanAWLe ContelOCullyCTurnerDLet al. Quantified energy dissipation rates in the terrestrial bow shock: 2. Waves and dissipation. J Geophys Res. (2014) 119:6475–95. 10.1002/2014JA019930

  • CrossRef
  • Google Scholar

• 7.

  PapadopoulosK. Electron heating in superhigh mach number shocks. Astrophys Space Sci. (1988) 144:535–47. 10.1007/BF00793203

  • CrossRef
  • Google Scholar

• 8.

  ShimadaNHoshinoM. Electron heating and acceleration in the shock transition region: background plasma parameter dependence. Phys Plasmas. (2004) 11:1840–9. 10.1063/1.1652060

  • CrossRef
  • Google Scholar

• 9.

  MatsumotoYAmanoTHoshinoM. Electron accelerations at high mach number shocks: two-dimensional particle-in-cell simulations in various parameter regimes. Astphys J. (2012) 755:109. 10.1088/0004-637X/755/2/109

  • CrossRef
  • Google Scholar

• 10.

  HoshinoMShimadaN. Nonthermal electrons at high mach number shocks: electron shock surfing acceleration. Astphys J. (2002) 572:880–7. 10.1086/340454

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11.

  SchmitzHChapmanSCDendyRO. Electron preacceleration mechanisms in the foot region of high Alfvénic mach number shocks. Astphys J. (2002) 579:327–36. 10.1086/341733

  • CrossRef
  • Google Scholar

• 12.

  AmanoTHoshinoM. Electron shock surfing acceleration in multidimensions: two-dimensional particle-in-cell simulation of collisionless perpendicular shock. Astphys J. (2009) 690:244–51. 10.1088/0004-637X/690/1/244

  • CrossRef
  • Google Scholar

• 13.

  UmedaTYamaoMYamazakiR. Electron acceleration at a low mach number perpendicular collisionless shock. Astphys J. (2009) 695:574–9. 10.1088/0004-637X/695/1/574

  • CrossRef
  • Google Scholar

• 14.

  TrolandTHHeilesC. Interstellar magnetic field strengths and gas densities: observational and theoretical perspectives. Astphys J. (1986) 301:339. 10.1086/163904

  • CrossRef
  • Google Scholar

• 15.

  ComişelH. On electron adiabaticity in collisionless shocks. Front Phys. (2016) 4:29. 10.3389/fphy.2016.00029

  • CrossRef
  • Google Scholar

• 16.

  FredricksRWCoronitiFVKennelCFScarfFL. Fast time-resolved spectra of electrostatic turbulence in the earth's bow shock. Phys Rev Lett. (1970) 24:994–8. 10.1103/PhysRevLett.24.994

  • CrossRef
  • Google Scholar

• 17.

  RodriguezPGurnettDA. Correlation of bow shock plasma wave turbulence with solar wind parameters. J Geophys Res. (1976) 81:2871. 10.1029/JA081i016p02871

  • CrossRef
  • Google Scholar

• 18.

  FormisanoVTorbertR. Ion acoustic wave forms generated by ion-ion streams at the Earth's bow shock. Geophys Res Lett. (1982) 9:207–10. 10.1029/GL009i003p00207

  • CrossRef
  • Google Scholar

• 19.

  BalikhinMWalkerSTreumannRAlleyneHKrasnoselskikhVGedalinMet al. Ion sound wave packets at the quasiperpendicular shock front. Geophys Res Lett. (2005) 32:L24106. 10.1029/2005GL024660

  • CrossRef
  • Google Scholar

• 20.

  HullAJLarsonDEWilberMScudderJDMozerFSRussellCTet al. Large-amplitude electrostatic waves associated with magnetic ramp substructure at Earth's bow shock. Geophys Res Lett. (2006) 33:L15104. 10.1029/2005GL025564

  • CrossRef
  • Google Scholar

• 21.

  GoodrichKAErgunRSchwartzSJWilsonLBNewmanDWilderFDet al. MMS observations of electrostatic waves in an oblique shock crossing. J Geophys Res. (2018) 123:9430–42. 10.1029/2018JA025830

  • CrossRef
  • Google Scholar

• 22.

  MatsumotoHKojimaHKasabaYMiyakeTAndersonRMukaiT. Plasma waves in the upstream and bow shock regions observed by geotail. Adv Space Res. (1997) 20:683–93. 10.1016/S0273-1177(97)00456-0

  • CrossRef
  • Google Scholar

• 23.

  BaleSDKelloggPJLarsenDELinRPGoetzKLeppingRP. Bipolar electrostatic structures in the shock transition region: evidence of electron phase space holes. Geophys Res Lett. (1998) 25:2929–32. 10.1029/98GL02111

  • CrossRef
  • Google Scholar

• 24.

  BaleSDHullALarsonDELinRPMuschiettiLKelloggPJet al. Electrostatic turbulence and debye-scale structures associated with electron thermalization at collisionless shocks. Astphys J Lett. (2002) 575:L25–8. 10.1086/342609

  • CrossRef
  • Google Scholar

• 25.

  HobaraYWalkerSNBalikhinMPokhotelovOAGedalinMKrasnoselskikhVet al. Cluster observations of electrostatic solitary waves near the Earth's bow shock. J Geophys Res. (2008) 113:A05211. 10.1029/2007JA012789

  • CrossRef
  • Google Scholar

• 26.

  VaskoIYMozerFSKrasnoselskikhVVArtemyevAVAgapitovOVBaleSDet al. Solitary waves across supercritical quasi-perpendicular shocks. Geophys Res Lett. (2018) 45:5809–17. 10.1029/2018GL077835

  • CrossRef
  • Google Scholar

• 27.

  WangRVaskoIYMozerFSBaleSDArtemyevAVBonnellJWet al. Electrostatic turbulence and Debye-scale structures in collisionless shocks. Astphys J Lett. (2020) 889:L9. 10.3847/2041-8213/ab6582

  • CrossRef
  • Google Scholar

• 28.

  BrenemanAWCattellCAKerstenKParadiseASchreinerSKelloggPJet al. STEREO and wind observations of intense cyclotron harmonic waves at the Earth's bow shock and inside the magnetosheath. J Geophys Res. (2013) 118:7654–64. 10.1002/2013JA019372

  • CrossRef
  • Google Scholar

• 29.

  WalkerSNBalikhinMAAlleyneHSCKHobaraYAndréMDunlopMW. Lower hybrid waves at the shock front: a reassessment. Ann Geophys. (2008) 26:699–707. 10.5194/angeo-26-699-2008

  • CrossRef
  • Google Scholar

• 30.

  FuselierSAGurnettDA. Short wavelength ion waves upstream of the earth's bow shock. J Geophys Res. (1984) 89:91–104. 10.1029/JA089iA01p00091

  • CrossRef
  • Google Scholar

• 31.

  SchamelH. Electron holes, ion holes and double layers. Electrostatic phase space structures in theory and experiment. Phys Rep. (1986) 140:161–91. 10.1016/0370-1573(86)90043-8

  • CrossRef
  • Google Scholar

• 32.

  HutchinsonIH. Electron holes in phase space: what they are and why they matter. Phys Plasmas. (2017) 24:055601. 10.1063/1.4976854

  • CrossRef
  • Google Scholar

• 33.

  MozerFSAgapitovOVGilesBVaskoI. Direct observation of electron distributions inside millisecond duration electron holes. Phys Rev Lett. (2018) 121:135102. 10.1103/PhysRevLett.121.135102

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34.

  AkimotoKWinskeD. Ion-acoustic-like waves excited by the reflected ions at the earth's bow shock. J Geophys Res. (1985) 90:12095–103. 10.1029/JA090iA12p12095

  • CrossRef
  • Google Scholar

• 35.

  FuselierSAGaryPSThomsenMFBameSJGurnettDA. Ion beams and the ion/ion acoustic instability upstream from the Earth's bow shock. J Geophys Res. (1987) 92:4740–4. 10.1029/JA092iA05p04740

  • CrossRef
  • Google Scholar

• 36.

  MuschiettiLRothICarlsonCWErgunRE. Transverse instability of magnetized electron holes. Phys Rev Lett. (2000) 85:94–7. 10.1103/PhysRevLett.85.94

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37.

  HutchinsonIH. Transverse instability magnetic field thresholds of electron phase-space holes. Phys Rev E. (2019) 99:053209. 10.1103/PhysRevE.99.053209

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38.

  RussellCTAndersonBJBaumjohannWBromundKRDearbornDFischerDet al. The magnetospheric multiscale magnetometers. Space Sci Rev. (2016) 199:189–256. 10.1007/s11214-014-0057-3

  • CrossRef
  • Google Scholar

• 39.

  ErgunRETuckerSWestfallJGoodrichKAMalaspinaDMSummersDet al. The axial double probe and fields signal processing for the MMS mission. Space Sci Rev. (2016) 199:167–88. 10.1007/s11214-014-0115-x

  • CrossRef
  • Google Scholar

• 40.

  LindqvistPAOlssonGTorbertRBKingBGranoffMRauDet al. The spin-plane double probe electric field instrument for MMS. Space Sci Rev. (2016) 199:137–65. 10.1007/s11214-014-0116-9

  • CrossRef
  • Google Scholar

• 41.

  PollockCMooreTJacquesABurchJGlieseUSaitoYet al. Fast plasma investigation for magnetospheric multiscale. Space Sci Rev. (2016) 199:331–406. 10.1007/s11214-016-0245-4

  • CrossRef
  • Google Scholar

• 42.

  MozerFS. DC and low-frequency double probe electric field measurements in space. J Geophys Res. (2016) 121:10. 10.1002/2016JA022952

  • CrossRef
  • Google Scholar

• 43.

  ScudderJD. A review of the physics of electron heating at collisionless shocks. Adv Space Res. (1995) 15:181–223. 10.1016/0273-1177(94)00101-6

  • CrossRef
  • Google Scholar

• 44.

  SchwartzSJHenleyEMitchellJKrasnoselskikhV. Electron temperature gradient scale at collisionless shocks. Phys Rev Lett. (2011) 107:215002. 10.1103/PhysRevLett.107.215002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45.

  VinasAFScudderJD. Fast and optimal solution to the ‘Rankine-Hugoniot problem’. J Geophys Res. (1986) 91:39–58. 10.1029/JA091iA01p00039

  • CrossRef
  • Google Scholar

• 46.

  SonnerupBUÖScheibleM. Minimum and maximum variance analysis. In: PaschmannGDalP editors. Analysis Methods for Multi-Spacecraft Data.ISSI Scientific Reports Series. Noordwijk: ESA Publications Division (1998). p. 185–220.

  • Google Scholar

• 47.

  LeroyMMWinskeDGoodrichCCWuCSPapadopoulosK. The structure of perpendicular bow shocks. J Geophys Res. (1982) 87:5081–94. 10.1029/JA087iA07p05081

  • CrossRef
  • Google Scholar

• 48.

  Kofoed-HansenOPecseliHLTrulsenJ. Coherent structures in numerically simulated plasma turbulence. Phys Scripta. (1989) 40:280–94. 10.1088/0031-8949/40/3/004

  • CrossRef
  • Google Scholar

• 49.

  BørveSPécseliHLTrulsenJ. Ion phase-space vortices in 2.5-dimensional simulations. J Plasma Phys. (2001) 65:107–29. 10.1017/S0022377801008947

  • CrossRef
  • Google Scholar

• 50.

  GoldmanMVNewmanDLErgunRE. Phase-space holes due to electron and ion beams accelerated by a current-driven potential ramp. Nonlin Process Geophys. (2003) 10:37–44. 10.5194/npg-10-37-2003

  • CrossRef
  • Google Scholar

• 51.

  PapadopoulosK. Ion thermalization in the Earth's bow shock. J Geophys Res. (1971) 76:3806. 10.1029/JA076i016p03806

  • CrossRef
  • Google Scholar

• 52.

  BiskampDWelterH. Structure of the Earth's bow shock. J Geophys Res. (1972) 77:6052. 10.1029/JA077i031p06052

  • CrossRef
  • Google Scholar

• 53.

  SagdeevRZGaleevAA. Nonlinear Plasma Theory.New York: W. A. Benjamin (1969).

  • Google Scholar

• 54.

  OhiraYTakaharaF. Oblique ion two-stream instability in the foot region of a collisionless shock. Astphys J. (2008) 688:320–6. 10.1086/592182

  • CrossRef
  • Google Scholar

• 55.

  GarySPOmidiN. The ion-ion acoustic instability. J Plasma Phys. (1987) 37:45–61. 10.1017/S0022377800011983

  • CrossRef
  • Google Scholar

• 56.

  GoodrichKAErgunRSchwartzSJWilsonLBJohlanderANewmanDet al. Impulsively reflected ions: a plausible mechanism for ion acoustic wave growth in collisionless shocks. J Geophys Res. (2019) 124:1855–65. 10.1029/2018JA026436

  • CrossRef
  • Google Scholar

• 57.

  ThomsenMFBarrHCGarySPFeldmanWCColeTE. Stability of electron distributions within the earth's bow shock. J Geophys Res. (1983) 88:3035–45. 10.1029/JA088iA04p03035

  • CrossRef
  • Google Scholar

• 58.

  GedalinM. Two-stream instability of electrons in the shock front. Geophys Res Lett. (1999) 26:1239–42. 10.1029/1999GL900239

  • CrossRef
  • Google Scholar

• 59.

  KhotyaintsevYVVaivadsAAndréMFujimotoMRetinòAOwenCJ. Observations of slow electron holes at a magnetic reconnection site. Phys Rev Lett. (2010) 105:165002. 10.1103/PhysRevLett.105.165002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 60.

  GrahamDBKhotyaintsevYVVaivadsAAndréM. Electrostatic solitary waves and electrostatic waves at the magnetopause. J Geophys Res. (2016) 121:3069–92. 10.1002/2015JA021527

  • CrossRef
  • Google Scholar

• 61.

  NorgrenCAndréMVaivadsAKhotyaintsevYV. Slow electron phase space holes: magnetotail observations. Geophys Res Lett. (2015) 42:1654–61. 10.1002/2015GL063218

  • CrossRef
  • Google Scholar

• 62.

  OmuraYMatsumotoHMiyakeTKojimaH. Electron beam instabilities as generation mechanism of electrostatic solitary waves in the magnetotail. J Geophys Res. (1996) 101:2685–98. 10.1029/95JA03145

  • CrossRef
  • Google Scholar

• 63.

  MorseRLNielsonCW. One-, two-, and three-dimensional numerical simulation of two-beam plasmas. Phys Rev Lett. (1969) 23:1087–90. 10.1103/PhysRevLett.23.1087

  • CrossRef
  • Google Scholar

• 64.

  BunemanO. Instability, turbulence, and conductivity in current-carrying plasma. Phys Rev Lett. (1958) 1:8–9. 10.1103/PhysRevLett.1.8

  • CrossRef
  • Google Scholar

• 65.

  NorgrenCAndréMGrahamD. B.KhotyaintsevYuV.VaivadsA.Slow electron holes in multicomponent plasmas. Geophys Res Lett. (2015) 42:7264–72. 10.1002/2015GL065390

  • CrossRef
  • Google Scholar

• 66.

  ZakharovVEKuznetsovEA. Three-dimensional solitons. Soviet J Exp Theor Phys. (1974) 39:285.

  • Google Scholar

• 67.

  TranM. Q. Ion acoustic solitons in a plasma: a review of their experimental properties and related theories. Phys Scripta. (1979) 20:317–27. 10.1088/0031-8949/20/3-4/004

  • CrossRef
  • Google Scholar