Relativistic electron acceleration on the bow shock of Jupiter and past

Data

Observations for this research are from the Juno spacecraft of NASA³². The energetic particle knowledge are supplied by the JEDI⁴⁷, which measures ions and electrons from about 30 keV to 1 MeV with an vitality decision of round 20%. JEDI consists of three similar sensor heads (JEDI90, JEDI180 and JEDI270) distributed across the spacecraft to optimize pitch angle protection over a 160° × 12° subject of view with an angular decision of about 18°. The first two vitality bins of JEDI used on this research are contaminated and aren’t included within the evaluation, leading to 4 vitality bins protecting roughly 100 keV to 1 MeV, as proven in Fig. 2. Lower-energy ion and electron observations are obtained from the JADE⁴⁸. JADE consists of two electron sensors (JADE-E) and an ion sensor (JADE-I), each measuring ions with vitality per cost from 10 eV q⁻¹ to 46.5 keV q⁻¹ throughout 64 vitality channels and electrons with vitality per cost from 30 eV q⁻¹ to 32 keV q⁻¹, with a time decision that’s mode dependent and corresponds to about 2 min within the offered occasion. Magnetic subject vector knowledge are sourced from the Magnetic Field Investigation (MAG) instrument⁴⁹, which makes use of two fluxgate magnetometers to offer measurements with a temporal decision of 1 s. All knowledge are offered within the JSO coordinate system, a Jupiter-centred body during which the x-axis factors to the Sun, the y-axis is within the anti-direction of the orbital movement of Jupiter and the z-axis completes the right-handed system⁵⁰.

Data post-processing and density calculations

The uncooked instrument knowledge have been processed to generate the merchandise used on this evaluation. The JEDI energy-time spectrograms (Fig. 1c) have been created by averaging knowledge from all three sensors and all look instructions. During the commentary interval, the instrument operated in a low-resolution mode, binning counts into six logarithmically spaced vitality channels from about 30 keV to 1 MeV and into 300 s time bins. The rely charges related to the transient occasion, starting from about 20 to 60 counts per second, are thought of statistically vital. The electron vitality effectivity correction detailed in ref. ⁵¹ was utilized, though its impact is minimal within the low-radiation surroundings close to the magnetospheric boundary of Jupiter. For JADE, proton densities have been derived from JADE-I knowledge utilizing a numerical integration technique on SPECIES=3 knowledge⁵². Although JADE-I isn’t optimized for photo voltaic wind measurements⁵³, this technique has been proven to be per forward-modelled Maxwellian suits for comparable occasions²⁸. The omnidirectional differential quantity intensities for JADE-E have been calculated by averaging the noticed intensities over 48 look instructions, that are binned onboard within the low charge science mode of the instrument⁴⁸.

Bow shock and foreshock transient characterization

In Extended Data Fig. 1, a magnified timeseries of the foreshock transient interval (11:30–13:30 UTC) is proven. Energetic particle intensification and plasma density depletion start at about 12:30 UTC, with a localized compression marking the trailing fringe of the construction at roughly 12:50 UTC, typical options of foreshock transients^4,5,22,54,55.

To higher characterize the worldwide surroundings throughout this encounter, we use the native magnetic subject circumstances and the shock regular vector estimated in ref. ⁵⁶. Using this, we acquire a standard vector of [0.77, 0.45, −0.44], per the duskward Juno location. The orientation of the magnetic subject with respect to this regular is proven in Extended Data Fig. 1e, suggesting that the shock orientation transitions from an indirect or quasi-parallel to a quasi-perpendicular one. Specifically, in the course of the formation and commentary of the transient itself, the orientation turns into much more quasi-parallel. This shock geometry (with θ_Bn ≲ 60°) is anticipated to supply substantial populations of foreshock suprathermal particles^57,58,59 related to the formation of foreshock transients^4,5,22,54.

Particle knowledge additional help this interpretation. The presence of diffuse, isotropic suprathermal ions and electrons signifies that the spacecraft is residing throughout the foreshock area⁶⁰. Specifically, the pitch angle distributions (PADs) of ions and particularly electrons present a transparent isotropic inhabitants of accelerated particles. These PADs reveal that particles are distributed throughout all pitch angles, a signature of well-scattered populations throughout the foreshock. This is in settlement with traits of accelerated electrons noticed throughout foreshock transients at Earth^5,18,19. An illustration of the surroundings and related transient is proven in Extended Data Fig. 2.

Focusing on the foreshock transient (12:30–12:50 UT), the electron PAD signature reveals a development because the transient passes by means of the spacecraft. This signature means that particles are accelerated within the approaching area, peaking throughout the transient and ceasing because the spacecraft exits the construction and the sphere rotates to a quasi-perpendicular regime after 12:50 UT. This strongly helps a neighborhood acceleration mechanism as a result of if the supply was exterior, energetic electrons could be observable over wider intervals. Instead, their strict localization to the transient construction implies they’re generated in situ quite than being remote-sensed. Regarding the broader spatial context, based mostly on spacecraft velocity (about 4 km s⁻¹) and the interval period, we estimate that Juno was residing roughly 1R_J upstream of the bow shock (Fig. 2, inset). This serves as an approximate estimate, because the bow shock at planetary flanks can change location quickly. This estimate is per observations at Earth, during which transients are noticed at round 1−4R_E (refs. ^19,61,62).

To decide the precise geometry and scale of the noticed foreshock transient, we first established its orientation utilizing minimal variance evaluation (MVA) on the magnetic subject vector knowledge within the JSO coordinate system⁶³. This method identifies the principal axes of the variance of the magnetic subject by discovering the eigenvalues (λ_max ≥ λ_int ≥ λ_min) and the corresponding eigenvectors of the covariance matrix of the magnetic subject over the interval containing the transient crossing. The eigenvector related to the minimal eigenvalue (n_MVA) is interpreted as the conventional path to the boundary of the transient, assuming a quasi-planar construction. The validity of this regular was confirmed by guaranteeing a big ratio of the intermediate to minimal eigenvalues (λ_int/λ_min ≫ 1). With the boundary regular established, we then estimated the dimensions measurement of the transient, L, alongside this path utilizing the single-spacecraft timing technique. The scale is calculated as L = |v_sw ⋅ n_MVA| × Δt, the place v_sw is the upstream photo voltaic wind velocity and Δt is the measured period of the passage of the spacecraft by means of the construction. Finally, the convection electrical subject −V × B factors in the direction of the transient sheet, which permits particles to pay attention and kind the noticed transient. The total methodological strategy we adopted is a standardized course of usually achieved when single spacecraft in situ observations can be found^4,5,18,28. Specifically, for our case, we used a typical upstream photo voltaic wind velocity of v_sw = [400, 0, 0] km s⁻¹ in JSO coordinates, which is in settlement with estimations of velocity throughout that interval, and calculated the dimensions as L = |v_sw ⋅ n_MVA| × Δt, the place Δt was taken as a 15-min period of the passage of the spacecraft in the course of the transient occasion. It ought to be famous that this 15-min interval, whereas comparatively conservative, supplies a practical vary of values for the spatial scale evaluation (described under). The consequence of this evaluation is supplied in Extended Data Table 1.

Finally, a power-law fitted to the energetic tail of the JADE and JEDI knowledge (≥10 keV) in the course of the foreshock transient outcomes (Fig. 2b) in a spectral index of P ≈ −1.85 ± 0.2 with the precise 95% confidence interval being decided utilizing a non-parametric bootstrap evaluation with 1,000 iterations⁶⁴. This worth suggests an acceleration course of with a signature and effectivity much like that of DSA. The obtained index is well-bounded by the canonical DSA restrict of P = −1.5, a function of environment friendly acceleration, and can also be per the anticipated spectral softening from −1.5 (non-relativistic) in the direction of −2 for electrons at relativistic energies at ≳1 MeV, as they’re above the electron relaxation mass vitality⁶⁵.

Spatial scales and vitality limits

The most vitality attainable by a charged particle is essentially constrained by the bodily properties of the accelerating surroundings. This restrict is named the Hillas criterion³, which relates the utmost particle vitality to the accessible potential drop throughout the system. For a attribute magnetic subject B and movement velocity V, the induced motional electrical subject creates a possible distinction throughout a scale L that in the end limits the utmost attainable vitality of a particle, with cost q, to E_max = qBLV.

In the particular context of diffusive shock acceleration, this restrict will be expressed as a confinement situation requiring that the upstream diffusion size of the particle, L_d, stays similar to or smaller than the system measurement L (ref. ³⁴). The diffusion size (L_d) will be estimated by means of the expression L_d ≈ D/V_sh, the place D is the spatial diffusion coefficient for the utmost vitality, and V_sh is the rate of the shock. For sturdy scattering, as is commonly assumed within the foreshock areas of planetary bow shocks^2,5, the scattering approaches the Bohm restrict, at which the diffusion coefficient is (Dapprox frac{1}{3}v{r}_{{rm{g}}}), with v being the relativistic velocity of the particle and r_g its gyroradius. By equating the diffusion size to the system measurement (L ≈ L_d), we recuperate the rate dependence inherent within the Hillas criterion: (Lapprox frac{1}{3}frac{v}{{V}_{{rm{sh}}}}{r}_{{rm{g}}}).

To derive a quantitative expression for the utmost vitality from this relationship, we specific the rate v of the particle and gyroradius r_g = p/(qB) when it comes to its whole vitality E_whole = E_max + mc², the place p is the momentum of the particle, q is the particle cost and B is the magnetic subject magnitude. The relativistic relations p = γmv and ({E}_{,{rm{whole}}}^{2}={(laptop)}^{2}+{(m{c}^{2})}^{2}) suggest (p=sqrt{{E}_{,{rm{whole}}}^{2}-{(m{c}^{2})}^{2}}/c) and v = pc²/E_whole. Substituting these into the confinement situation provides

which ends up in the quadratic equation for the particle’s whole vitality,

Solving this equation for the optimistic vitality root supplies the precise resolution for the utmost whole vitality {that a} size-limited shock acceleration area can produce:

the place the time period A = 3qBLV_sh describes the properties of the accelerator. The most kinetic vitality is then discovered by subtracting the remainder mass vitality of the particle, E_max = E_max,whole − mc². In the ultrarelativistic restrict (v → c and E_whole ≫ mc²), this relation reduces to the straightforward linear kind E_max ≃ αL, with α ≃ 3qBV_sh.

The calculation of vitality limits in equation (3) requires three parameters: the native upstream magnetic subject energy (B), the shock velocity (V_sh) and the attribute system measurement (L). In this research, we use the upstream magnetic subject for B and the relative velocity between the transient and the first shock for V_sh. The acceleration area measurement (L) is described by the spatial extent of the transient, which is sufficiently equated to the precursor or foreshock area.

To mannequin the connection between the acceleration area measurement (L) and the attribute system measurement (S) as proven in Fig. 3, we assumed a power-law dependency of the shape L = okay ⋅ S^m (ref. ²⁸). Two separate power-law fashions based mostly on planetary observations are developed: a ‘typical’ mannequin and an ‘extreme’ mannequin. For the standard mannequin, we used the standoff distance of every planet because the system measurement (S) and its typical noticed acceleration area measurement because the related worth for L_typ. For the intense mannequin, we used the identical standoff distances (S) however paired them with the utmost noticed acceleration area sizes (L_ext). We carried out a linear match for every mannequin within the log–log area utilizing extraordinary least squares as applied by the statsmodels (v.0.14.4) of Python library. The ensuing suits are proven in Fig. 3a.

The extension from planetary to astrophysical scales follows three factors of reasoning: (1) We set up empirically that the acceleration area measurement (L) scales systematically with the worldwide shock measurement (S) throughout planetary environments, the place in situ observations verify that large-scale foreshock transients are the first acceleration websites. (2) We apply the Hillas criterion to narrate this acceleration scale (L) to the utmost particle vitality (E_max), a relationship validated by means of the offered Juno observations, earlier analysis and not too long ago proven to be relevant at planetary scales^5,34. (3) We mix these scaling relations below the premise that foreshock transient-driven acceleration operates universally throughout collisionless shocks, impartial of whether or not the system is a planetary bow shock, protostellar jet or supernova remnant. This last step represents an extrapolation past what’s testable with in situ experimentation and depends on theoretical expectations relating to the physics of collisionless shocks in numerous astrophysical contexts, additional mentioned under. To mannequin the utmost particle energies for the planetary methods proven in Fig. 3(b), we set the acceleration measurement (L) to the utmost noticed worth (L_ext). For the astrophysical objects, the place direct observations of L are unavailable, we used our excessive match mannequin to estimate a worth. Specifically, we calculated the higher boundary of the 95% prediction interval for a brand new commentary on the object’s system measurement (S). This interval accounts for each the uncertainty within the fitted mannequin and the inherent variability within the knowledge. The ensuing vitality vary for every object, subsequently, displays the uncertainty propagated from its native shock velocity and magnetic subject parameters⁶⁶. This strategy supplies an estimate for the utmost achievable vitality that’s constrained by direct commentary for the planets and by a sturdy statistical extrapolation for the astrophysical shocks.

Physical justification for astrophysical applicability

The extension of the planetary-derived scaling to astrophysical environments essentially entails an extrapolation past verifiable limits, as in situ measurements of foreshock localized processes at kiloparsec distances aren’t possible for the foreseeable future. However, a number of bodily arguments help this extrapolation as an affordable theoretical expectation, grounded in our understanding of collisionless shock physics and the properties of astrophysical environments.

First, the astrophysical shocks that we contemplate (the protostellar jet HH 211 and the supernova remnants SN 1987A and SN 1006) characterize curved, collisionless shocks, identified from numerical simulations and oblique observational proof to exhibit notable foreshock areas of their quasi-parallel geometries^2,8. The curved geometry of those shocks ensures that quasi-parallel configurations are naturally current throughout substantial parts of the shock floor, much like planetary bow shocks.

Second, though the expected acceleration area scales of 10⁸–10¹⁰ km for these astrophysical methods can’t be instantly resolved, our estimates characterize conservative decrease bounds. These values are per impartial estimates of foreshock sizes^67,68. The consistency between our extrapolated predictions and these observational constraints supplies help for the validity of our scaling strategy, even within the absence of direct spatial decision of the acceleration areas themselves.

Third, the formation of large-scale foreshock transient can come up spontaneously by means of plasma instabilities and wave-particle interactions, whereas transient formation can also be facilitated by large-scale discontinuities within the upstream medium²². These variable upstream media are an anticipated function of astrophysical quasi-parallel shock environments. The turbulent interstellar medium comprises magnetized filaments and buildings extending to parsec scales^69,70,71. Furthermore, supernova shock waves are anticipated to work together with advanced circumstellar buildings created by the progenitor star itself. These embrace stellar wind-blown bubbles and beforehand ejected dense shells of fabric, which may lengthen to parsec scales^72,73,74. When the increasing supernova shock encounters these pre-existing buildings, it might generate large-scale magnetic subject discontinuities. These encounters present excellent circumstances for seeding transient phenomena and driving particle acceleration by means of transient processes. Finally, the very excessive shock speeds noticed at supernova remnants permit Alfvénic Mach numbers (M_A) to achieve values of the order of 10²–10³, inflicting significantly sturdy foreshock areas upstream of their quasi-parallel geometries which might be extremely beneficial for the formation and persistence of large-scale foreshock buildings. We observe, nevertheless, that the utmost extent of those buildings could also be regulated by the coherence size of the upstream magnetic subject (L_B). For remnants within the galactic disk (for instance, SN 1572), during which the shock radius can exceed L_B, the acceleration area may be affected by the coherence scale, presumably making a plateau within the acceleration scale size. This is in distinction to high-latitude methods resembling SN 1006 or younger SRNs resembling SN7D21987A, during which L_B > S. Overall, the interaction between tangled fields and foreshock improvement within the regime during which S > L_B stays an open query for future investigation.

Although these arguments set up the bodily plausibility of extending our planetary framework to astrophysical scales, we emphasize that this extrapolation stays tentative within the absence of direct observational affirmation. Nevertheless, the interior consistency of our framework, validated at planetary scales and yielding predictions for SN 1006 that match noticed most energies of about 100 TeV, supplies encouraging help for the underlying bodily image.

Parameters of Solar System planets

To apply our acceleration mannequin, we first outlined the bodily parameter area for every surroundings. This area consists of the attribute measurement (S), the acceleration area scale measurement (L), the shock velocity (V_sh) and the upstream magnetic subject energy (B) of the system. For the planetary environments, the system measurement S is the standoff distance of the bow shock, and L is the standard scale of a foreshock transient. These values have been adopted from established literature^28,75, with the dataset for Jupiter supplemented by the occasion offered on this work. Minor refinements to the values have been made on the idea of the latest standoff distance statistics and observations from every planetary surroundings^{76,77,78,79,80}. A full vary for the shock velocity, V_sh, was chosen for all environments to account for the number of dynamic circumstances noticed at interplanetary shocks. This vary is predicated on the standard relative speeds between the photo voltaic wind plasma and transient compressive buildings, as documented in a number of earlier works, with the higher restrict primarily describing typical photo voltaic wind velocity^4,5,28,81. Finally, the vary for the upstream magnetic subject, B, at every planet was estimated utilizing the Parker spiral mannequin, which describes the evolution of the heliospheric magnetic subject energy with rising distance from the Sun⁸².

Parameters for HH 211

HH 211 lies at roughly 1,000 gentle years away and drives a slim bipolar jet with an related bow shock. James Webb Space Telescope measurements confirmed speeds of about 80–100 km s⁻¹ (determine 5 in ref. ⁷). Polarimetric ALMA/SMA observations have proven envelope magnetic fields of 40–100 μG (4–10 nT) inside about 700 AU (refs. ^83,84). The projected separation from the protostar to the bow shock is about 1,000 AU (determine 1 in ref. ⁷), which we take because the efficient standoff distance. To account for variability, we took 700 AU because the system measurement (S), and a shock velocity between 50 km s⁻¹ and 150 km s⁻¹ based mostly on the noticed outer shock. The comparatively chilly surroundings upstream of those propagating jets permits the native Alfvén Mach quantity (M_A) to be comparatively excessive, subsequently inflicting the formation of a robust foreshock precursor.

Parameters of SN 1987A

SN 1987A within the Large Magellanic Cloud has been extensively monitored by means of radio and X-ray observations. During the interplay with the dense equatorial ring, the ahead shock decelerated to round 2,300 km s⁻¹, later re-accelerating to three,500–3,600 km s⁻¹ (refs. ^85,86,87); for our evaluation, we used a conservative vary between 2,000 km s⁻¹ and 4,000 km s⁻¹. For the upstream interstellar magnetic subject, we count on values of the order of 1–5 μG (0.1–0.5 nT)^88,89. Finally, for an equivalence of a standoff distance of the increasing shock (or arc), we used the radius of the blast wave of 0.1 laptop, leading to an roughly 200,000 AU system measurement (S).

Parameters of SN 1006

SN 1006 is a younger, shell-type supernova remnant that gives compelling proof for environment friendly particle acceleration to very excessive energies, much like different historic remnants, resembling Tycho⁹⁰, Cas A⁹¹ and RX J1713.7–3946⁹². Its shock velocity has been well-constrained by observations of its growth, with estimates usually putting it within the vary of three,500–4,500 km s⁻¹. The energy of the upstream ambient magnetic subject values is taken to be between 1 μG and a pair of μG (0.1–0.2 nT). The system measurement of the remnant is roughly 9–10 parsecs, which interprets to a system measurement (S) of roughly 2  × 10⁶ AU. A novel function of SN 1006 is the sturdy observational proof, derived from its X-ray synchrotron and TeV emission, that particles are being accelerated to energies of the order of 100 TeV inside its shocks^8,36,38. This makes SN 1006 a great take a look at case for evaluating acceleration frameworks and supplies a direct benchmark for the predictions of our mannequin.

The very excessive speeds noticed at SN shocks permit Alfvénic Mach numbers (M_A) to achieve values of the order of 10³−10⁴ inflicting a very sturdy foreshock upstream of their quasi-parallel geometry. This parameter vary estimation is troublesome to constrain as increasing SN shocks change throughout their lifetime. However, the parameter area above supplies an affordable order of magnitude estimate and uncertainty vary for the utmost obtainable particle vitality. The parameters for the planetary environments that have been used for the becoming are proven in Extended Data Table 2, and all parameters used for our generalized most vitality mannequin are summarized in Extended Data Table 3. It ought to be famous that the obtained shock velocity vary for every planetary surroundings is predominantly pushed by the relative velocity of the photo voltaic wind with respect to the shock, whereas for astrophysical objects, it’s primarily dictated by the outward shock growth to the comparatively secure interstellar medium. Finally, though the upstream magnetic subject values used on this research could also be elevated by compression results throughout the foreshock or precursor area, the acceleration areas themselves are characterised by native magnetic depressions. Consequently, the magnetic subject inside these transient buildings is anticipated to be notably weaker than within the surrounding surroundings.