Acceleration of Electrons and Ions by an “Almost” Astrophysical Shock in the Heliosphere

Collisionless shock waves, ubiquitous in the Universe, are crucial for particle acceleration in various astrophysical systems. Currently, the heliosphere is the only natural environment available for their in situ study. In this work, we showcase the collective acceleration of electrons and ions by one of the fastest in situ shocks ever recorded, observed by the pioneering Parker Solar Probe at only 34.5 million km from the Sun. Our analysis of this unprecedented, near-parallel shock shows electron acceleration up to 6 MeV amidst intense multiscale electromagnetic wave emissions. We also present evidence of a variable shock structure capable of injecting and accelerating ions from the solar wind to high energies through a self-consistent process. The exceptional capability of the probe’s instruments to measure electromagnetic fields in a shock traveling at 1% the speed of light has enabled us, for the first time, to confirm that the structure of a strong heliospheric shock aligns with theoretical models of strong shocks observed in astrophysical environments. This alignment offers viable avenues for understanding astrophysical shock processes and the self-consistent acceleration of charged particles.


INTRODUCTION
Collisionless shock waves (CSWs), resulting from converging flows in tenuous plasma, are a fundamental phenomenon in plasma physics (Sagdeev 1966;Galeev 1976;Lembege et al. 2004) capable of heating the plasma and accelerating charged particles (Kennel et al. 1985;Blandford & Eichler 1987).While ubiquitous in diverse environments, from planetary bow shocks to supernova remnants (SNRs), our comprehension of particle energization (heating and acceleration) at these shocks remains incomplete (Malkov & Drury 2001;Lembege et al. 2004).It is mostly the radiation generated by energized electrons that make the detection of CSWs, like those in SNRs, possible (Helder et al. 2012).Direct in situ examination of astrophysical shocks is currently not possible, leading to a limited understanding of the exact mechanisms behind particle acceleration and ensuing radiation.At present, heliospheric shocks are the only natural CSWs accessible for direct in situ study, making them essential for comprehending shock structure, evolution, Corresponding author: Immanuel Christopher Jebaraj immanuel.c.jebaraj@gmail.comand particle energization mechanisms.In particular, the shock transition, which marks the finite region between the upstream and downstream flows, plays a crucial role in understanding the redistribution of kinetic energy (Balikhin et al. 1993;Krasnosselskikh et al. 1995;Agapitov et al. 2023).Consequently, heliospheric shocks, which may in many aspects be similar to young SNR shocks, have become a key focus of research (Terasawa 2003;Gedalin 2023).
In magnetized collisionless shock waves (CSWs), such as those found in supernova remnants (SNRs) and the heliosphere, particle motion is determined by the angle θ Bn between the upstream magnetic field and the shock normal.CSWs are classified into quasi-parallel (θ Bn < 45 • ) and quasi-perpendicular (θ Bn > 45 • ) based on this angle.Observational surveys have often found the brightest radio and X-ray synchrotron emissions in regions of SNRs where the shock is quasi-parallel (e.g.Giuffrida et al. 2022;Vink et al. 2022).They are also thought to be the most efficient particle accelerators (Vink 2020) and are of particular interest for in situ exploration.Quasi-parallel shocks have been understood for nearly 80 years (Moiseev & Sagdeev 1963;Kennel & Sagdeev 1967;Quest 1988), but observational efforts have not yet yielded conclusive results (Burgess et al. 2005;Balikhin et al. 2023).Another critical parameter in the study of CSWs is the Alfvén Mach number (M A ), which is directly related to the proportion of flow kinetic energy dissipated at the shock.It is defined as the ratio of the upstream flow speed to the characteristic wave speed in the plasma, known as the Alfvén speed (v A ).Typically, heliospheric shocks exhibit M A values less than 102 (Masters et al. 2013a), while SNR shocks can reach M A values around 10 3 (Vink 2020).A final important parameter is the plasma-β, which is the ratio of kinetic plasma pressure to magnetic field pressure.It is generally believed to control the evolutionary characteristics of shocks.Those developing in a low plasma-β regime tend to exhibit large amplitude overshoot magnetic fields (Russell et al. 1982).
SNR shocks such as Tycho and Vela Jr. expand with velocities around 3000 km s −1 , and are established accelerators of cosmic rays > 10 14 eV (100 TeV) by channeling up to ∼10% of their kinetic energy into particle acceleration (Helder et al. 2009).A key factor in the acceleration of high-energy particles is the shock's ability to confine particles between converging scattering centers, such as electromagnetic waves, on either side of the CSW.These scattering centers, fundamental to the theory of diffusive shock acceleration (DSA), enable particles to repeatedly cross the shock and gain energy in the process (Krymskii 1977;Axford et al. 1977;Blandford & Ostriker 1978;Bell 1978).This energy acquisition process persists until the particles can escape the shock's influence, conventionally dubbed, the duration of particle confinement (Hillas 1984).Planetary bow shocks, the strongest heliospheric shocks, are incapable of confining particles long enough to accelerate relativistic particles due to their small size (Hoppe & Russell 1982;Terasawa 2011).
CSWs driven by solar eruptive events such as coronal mass ejections (CMEs), can occasionally accelerate particles to energies of ≥ 10 9 eV (GeV) (Reames 1999).The most significant acceleration occurs in the relatively denser solar atmosphere when these shocks are fast but not very large, highlighting the importance of their early evolution (Afanasiev et al. 2018).This process is somewhat analogous to that of young SNR shocks, which are capable of accelerating particles to 10 15 eV (PeV) during their rapid evolution within a dense medium (Aharonian et al. 2019).A recent statistical study of the strongest CME-driven shocks from the past decade suggests that high-energy electrons are likely accelerated by the shock through a mechanism similar to that affecting ions (Dresing et al. 2022).This draws attention to the precise process of shock-induced particle injection and acceleration through DSA, processes that are actively studied (Malkov & Völk 1998;Malkov & Drury 2001).The complexity of this process is exacerbated by the significant mass disparity between ions and electrons, with protons being 1836 times heavier.This separation of scales necessitates distinct processes for the energization of protons and electrons at their respective scales (Lembege et al. 2004).
On March 13, 2023, during its fifteenth perihelion, NASA's Parker Solar Probe (PSP) mission (Fox et al. 2016) encountered an extraordinary interplanetary (IP) shock at merely 0.24 astronomical units (a.u.) from the Sun.This CME-driven shock, one of the fastest IP shocks ever recorded in situ, was traveling at 1% the speed of light and exhibited a near-parallel geometry.In this study, we detail the distinct features observed in the transition region of this shock and its impact on particle behavior, enhancing our understanding of strong shocks in both heliospheric and astrophysical contexts.This investigation was enabled by PSP overcoming engineering challenges of entering the Sun's atmosphere, coupled with the exceptional capabilities of its high fidelity instrument suites to measure such phenomena in great detail.Information about PSP's instrument suites used in this study can be found in the Appendix A. In Fig. 1(A,B), we present magnetic field and ion measurements during a 9-hour period (03:00 to 12:00 UT) on March 13, 2023.The coordinate system used throughout the manuscript is the inertial RTN (Radial-Tangential-Normal) system, where the radial component R is oriented along the Sun-spacecraft line, the transverse component T is defined to be orthogonal to the rotational axis of the Sun and the radial component, i.e., T = Ω ⊙ × R, while the normal component N completes the orthogonal right-handed triad and, in this case, is aligned with the normal of the ecliptic plane.The two panels provide details on the pre-event conditions and the in situ arrival of the shock driven by a large magnetic structure, the CME.The upstream bulk flow speed, measured by SPAN-i, remains consistent throughout the period, averaging 410 ± 40 km s −1 in the 15 minutes prior to the shock's arrival.However, this consistency does not extend to the number density, which is underestimated, as detailed in Appendix B.2. Similarly, the upstream magnetic field is averaged over the same period, while the downstream magnetic field is averaged over a subsequent 10-minute period; both are listed in Table 1.
PSP observed the shock propagating almost exactly radially outward, as indicated by the estimated shock normal (n RTN ), which is also listed in the table.Figure 1(C-E) provides a zoomed-in view of the shock's arrival at 07:13 UT, marked by a vertical line, and covers the 75 minutes leading up to it.During the shock's passage, the downstream plasma flow exceeded SPAN-i's measurement capabilities, as depicted in Fig. 1(B,D,E).By combining data from SPAN-i and EPI-Lo and utilizing the principle of mass flux conservation across the shock (detailed in Appendix C.3), we calculate the shock speed to be ∼ 2800 ± 300 km s −1 in the spacecraft frame.This classifies it as one of the fastest shocks ever recorded in situ, placing it in the same category as the one described in Russell et al. (2013).
Additionally, the occurrence of such a strong shock (M A ∼ 9.1 ± 1.35) in a low-β plasma environment (∼ 0.16 ± 0.031), along with its near-parallel geometry (8 • ± 4 • ), makes it an exceptionally rare in situ observation.A list of the most relevant plasma and shock parameters, along with their uncertainties is presented in Table 1.The uncertainty for each parameter is subsequently propagated when estimating related parameters, following the general fractional error propagation.For instance, the shock normal is derived as a average from minimum variance analysis (MVA) and the magnetic coplanarity technique (MCT), detailed in Appendix C.1.Additionally, unorthodox methodologies employed to estimate plasma density contribute to the uncertainty of the shock parameters, as discussed in Appendix B.2.A notable characteristic of quasi-parallel shocks is the patchwork of complex magnetic structures making up an extended transition region, typically in the order of 100s of ion-inertial lengths (λ i , the characteristic distance over which ions respond to electromagnetic changes.Burgess et al. 2005).This contrasts with quasi-perpendicular shocks, known for their sharp transition in the form of a ramp and overshoot (Krasnoselskikh et al. 2013).This complexity of quasi-parallel shocks is illustrated in Fig. 1(F-H), which display the magnetic field, its fluctuations, and the ion temperature over a 2-minute interval (07:12-07:14 UT).As we approach the shock transition, marked by the vertical line, there is a noticeable increase in the amplitude of magnetic field fluctuations (δB/B), indicative of intense upstream wave activity.The data also show that the magnitude of the magnetic field intensifies before the transition, likely as a result of shock mediation by energetic particles.This amplification of the magnetic field and the intensification of fluctuations is also predicted by theoretical models of strictly parallel shocks (e.g.Wang et al. 2022).The shock transition is characterized by a δB/B ≥ 0.5, preceding the ramp-overshoot region, which exhibits a δB/B ≥ 2. Figure 1H displays a sudden rise in ion temperature at the transition, consistent with the expected behavior for magnetohydrodynamic (MHD) shocks (Sagdeev 1966).While it is evident that the temperature sharply increases at the transition, the SPAN-i measurements do not accurately reflect downstream conditions.

Electron acceleration to ultra-relativistic energies
A particularly notable and unprecedented observation concerning this strong quasi-parallel shock is the local acceleration of electrons, with energies ranging from tens of keV (sub-relativistic) to 6 MeV (ultra-relativistic). Observations of the local acceleration of relativistic electrons (300 keV to 3 MeV) is exceptionally rare in IP shocks (Talebpour Sheshvan et al. 2023;Jebaraj et al. 2023) and is even less frequent at planetary bow shocks, regardless of their strength (Masters et al. 2013b).The occurrence of ultra-relativistic electrons at these shocks has never been documented before.The rarity of local electron acceleration at IP shocks can be attributed to a combination of scale-related challenges and a lack of high-fidelity instrumentation.Physical challenges, such as the energization of electrons from a thermal core, require specific processes that affect electrons at scales close to their gyrofrequency (ω ce ).These processes are related to fundamental issues of energy redistribution in CSWs (Balikhin et al. 1993;Schwartz et al. 2011).Additionally, their acceleration also necessitates the presence of oblique waves across a wide range of frequencies.Figure 2(A,B) presents energetic electron observations during a 90-minute period surrounding the shock's arrival.Figure 2A displays the counts per second (s −1 ) of relativistic electrons (>1 MeV) throughout this interval, while Fig. 2B provides differential energy fluxes of lower-energy electrons, ranging from 50 to 350 keV.Notably, there is a significant increase in electron fluxes and counts across a wide range of energy when approaching the shock (between 06:40 and 07:13 UT), thus illustrating its impact on electron acceleration.The observed range of electron energies, spanning sub-relativistic to ultra-relativistic, implies the presence of multiscale electromagnetic structures in and around the shock transition, affecting electrons across all these energies.Fig. 2B shows that while sub-relativistic electrons are significantly affected by the shock transition, relativistic and ultrarelativistic electrons seem largely unaffected.This distinction arises from the increasing radius of gyration (Larmor radius) of electrons with higher kinetic energy.The upstream data reveal fluctuations in electron fluxes at distinct times, specifically at 06:00 UT and 06:40 UT, suggesting localized upstream phenomena.Nonetheless, the limited resolution of electron observations restricts a more detailed analysis of these variations.

Self-consistent injection and acceleration of solar wind ions
In Fig. 2C, we showcase the ion fluxes measured by EPI-Lo's sunward-pointing detector (W3) at the shock, and up to 90 minutes prior, to investigate the motion of energetic ions.Measurements from all EPI-Lo wedges are presented in Figs. 8 and 9 of Appendix D. A striking observation is the intense population of ions streaming ahead of the shock, persisting for the entire 90-minute duration, peaking between 06:40 and 07:13 UT. Figure 2D illustrates the fluctuations in the magnetic field (δB/B) on a characteristic time scale of τ = 1 minute, corresponding to a frequency of approximately 0.02 Hz.This frequency roughly matches the resonance frequency for the lowest energy of the streaming ions further from the shock, namely ∼1 MeV.The fluctuations are categorized into transversal (δB ⊥ /B) and parallel (δB ∥ /B) components, with the former being significantly larger than the latter within the same timeframe.The technical details of these analyses can be found in Appendix E.
Another intriguing feature is the observed absence of flux at energies below ∼1 MeV, before 07:05 UT.This absence of flux extends to higher energies further from the shock, especially prior to 06:40 UT.In Fig. 2E, the omni-directional intensity spectra of protons are shown for three distinct time periods: at the shock (07:13-07:15 UT), near-upstream (07:10 UT), and far-upstream (06:40 UT).These spectra highlight a marked absence of flux near 1 MeV, followed by a significant flux increase towards 10 MeV, forming what can be described as a spectral roll-over at low energies.The significant absence of low-energy particles away from the shock may be due to the fact that only high-energy particles diffuse farther from the shock.Approaching the shock however, the spectrum evolves into a single power law with a spectral index of E −1 beyond ∼250 keV.Below approximately 250 keV, the spectral index becomes steeper, approaching E −2 .However, we refrain from interpreting this finding, as it is potentially erroneous to perform a spectral fit over a range that does not extend across at least an order of magnitude.
Figure 2F presents the Fourier power spectra for B ∥ and B ⊥ during a 33-minute interval upstream (06:40-07:13 UT) and compares it to the pre-event background (01:30-02:07 UT).For a detailed description of the spectral analysis, please refer to Appendix E.2.Generally, upstream power for both B ∥ and B ⊥ exceed pre-event levels by more than an order of magnitude at frequencies above the resonance frequency for particles of energy ∼1 MeV (f 1MeV res ).Below this frequency, we observe little difference between the upstream and pre-event power of B ⊥ .In contrast, there is an increase in B ∥ relative to the pre-event power above f 5MeV res .However, it is important to note that frequencies at or below 10 −2 Hz are subject to significant uncertainties due to various factors, as detailed in the Appendix E.2.
Zooming into the shock ramp, we note an uptick in both the amplitude and intensity of these fluctuations.In Fig. 2(G,H), we show the energetic protons and magnetic field fluctuations within a 2-minute window (07:12-07:14 UT) around the ramp.Here, we detect an intense proton population with energies between 100 and ∼250 keV which correspond to the energy range where the spectra becomes steeper in Fig. 2E.This occurs concurrently with highintensity, large-amplitude δB ⊥ /B ≥ 0.5 fluctuations.These fluctuations are estimated at τ = 5 seconds, corresponding to resonant frequencies above 0.2 Hz, which match those of protons in the 100s of keV range.The δB ∥ /B also show an increase as we approach the ramp, especially after 07:13 UT.These observations concur with prior studies on the challenges of particle confinement at quasi-parallel shocks leading to an extensive wave foreshock (Kennel et al. 1986).It also supports theoretical predictions that wave power declines with increasing distance from the shock owing to the reduced particle fluxes (Bell 1978;Vainio & Laitinen 2007).
Finally, the continuous detection of ions with energies less than 100 keV to 10 MeV near the shock is noteworthy.Figure 2C highlights data from the W3 detector of EPI-Lo, but all EPI-Lo detectors registered ions above 1 MeV for about 90 minutes before the shock arrival.In contrast, ≤300 keV protons were primarily observed near the shock transition.It is likely that these are either reflected or gyrating ions as previously suggested by Kennel (1981) and thus deviate from the E −1 law below ∼ 250 keV.Distinguishing between the two components would require higher-cadence particle measurements.However, in the case of near-specular reflection, the energy of the ions would primarily peak at twice the speed (four times the energy) of the upstream bulk flow (in the shock rest frame), as most of the ions are reflected only once (Gedalin 2016).Considering the quasi-parallel geometry of the shock, these ions can diffuse ahead of the shock and participate in DSA, thereby gaining significantly higher energies.An E −1 scaling law supports this process occurring in a lossless manner, meaning that the protons (ions) are confined in the vicinity of the shock without being lost.Conversely, this could also suggest that they are replenished at the same rate as they are lost, underscoring the efficiency of the shock in particle acceleration.

Structure of the quasi-parallel shock transition
In contrast to the infinitesimal boundary between two flows marking a shock in MHD, the shock transition is where processes crucial for regulating the distribution of kinetic energy occur, such as dispersion, dissipation, and non-linear steepening (Krasnoselskikh 1985;Krasnoselskikh et al. 2002).In Fig. 3, we present a comprehensive analysis of the shock transition, with further details available in Appendix F. Figure 3A illustrates the magnetic field components and magnitude over a 1.75-second interval (07:13:10.25-07:13:12.00),complemented by Fig. 3B, which depicts the electric field.
Figure 3(C,D) depict the dynamic spectra of magnetic and electric field perturbations, respectively, derived from the Morlet wavelet transform.These panels emphasize the shock ramp and its associated precursor wave train, notably within the 20-30 Hz frequency range (indicated by arrows) which is higher than the ion gyrofrequency (ω ci ∼ 1.4 Hz).The wave train observed here consists of circularly polarized waves, as seen in Fig. 3G, the plane containing the intermediate and maximum variance axis.Consequently, they also exhibit high planarity in the plane containing the maximum and minimum variance axis as shown in Fig. 3H.Such characteristics indicate that these are electromagnetic whistler mode waves, which are occasionally observed ahead of the shock ramp and are associated with fundamental shock processes (Sundkvist et al. 2012).These waves are right-hand circularly polarized in the spacecraft frame and occur between the ion and electron gyrofrequencies (ω ci < ω < ω ce ).As the wave number (k) increases and they approach ω ce , these waves can also increase in amplitude and may steepen.
The wave normal angle, depicted in Fig. 3E (θ kB , the angle between the wave normal and the background magnetic field direction), indicates that both the precursors and the ramp are quasi-perpendicular to the magnetic field.Their oblique propagation is an essential factor to be a standing wave in the shock frame.Figure 3F illustrates the radial component of the Poynting flux, estimated in the spacecraft frame, indicating an anti-sunward propagation direction of the observed wave activity near the shock.The technical material related to the analysis of the precursor whistler waves can be found in Appendix F.2. From this analysis, the estimated phase speed (v ph ) is approximately 400 km s −1 .Given that θ kB is approximately 75 • , the radial component is calculated to be around 2600±200 km s −1 .This analysis of the precursor whistler wave train lends credibility to the estimated shock speed.
More importantly, we demonstrate that the ramp and precursors display a quasi-perpendicular geometry.This finding is in line with theoretical predictions (Gedalin 1998;Gedalin et al. 2015), which suggest that a quasi-parallel shock becomes oblique, even quasi-perpendicular, due to the increase in transversal magnetic field amplitude in the ramp-overshoot region.However, prior to our study, observational support for this was limited (Balikhin et al. 2023).We have clearly demonstrated that the enhancement of the transversal magnetic field at the ramp-overshoot leads to a quasi-perpendicular geometry.The degree of enhancement is particularly noteworthy given that the shock develops in a low plasma-β environment, where large ramp-overshoot structures have previously not even been observed in either oblique or quasi-perpendicular obstacle shocks (Russell et al. 1982).

DISCUSSION
Strong collisionless shock waves, spanning a wide range of Mach numbers in the heliosphere, have been extensively studied to enhance our understanding of particle acceleration in astrophysical shocks, such as those in SNRs.However, in contrast to SNR shocks, even the strongest heliospheric shocks (planetary bow shocks) are ineffective accelerators of relativistic particles (Hoppe & Russell 1982).This inefficiency is largely attributed to their smaller size and the weak magnetization of the ambient medium, resulting in inadequate particle confinement (Terasawa 2011).Conversely, in- terplanetary traveling shocks are significantly larger but rarely as strong as planetary bow shocks.Remote observations have occasionally identified CME-driven shocks possessing strengths comparable to planetary bow shocks and speeds similar to some SNR shocks (Yurchyshyn et al. 2005).Particle acceleration by these shocks is anticipated to be highly efficient (Afanasiev et al. 2018).Measurements from PSP enable us to investigate these shocks in situ, unprecedentedly close to the Sun.This proximity significantly narrows the gap between our understanding of heliospheric shocks and SNRs, allowing for an in-depth study of efficient acceleration mechanisms.
In this study, we have presented extraordinary observations of a strong (M A ∼ 9.1 ± 1.35), fast (v sh ∼ 2800 ± 300 km s −1 ), and near-parallel (θ Bn ∼ 8 • ± 4 • ) traveling interplanetary shock.This shock, one of the fastest IP traveling shocks observed in situ to date, exhibits a magnetic structure that aligns with theoretical models.Figure 4A illustrates the clear rotation of the transverse magnetic field components at the shock transition, not obscured by nonlinear structures formed in the foreshock.This alignment with theoretical predictions is essential for making reasonable comparisons with SNR shocks.
As expected from such a prominent shock, it proved to be an efficient accelerator of particles.A notable observation is the local acceleration of electrons across a wide range of energies by such a near-parallel shock, often disregarded as potential accelerators of energetic electrons (Masters et al. 2013b).Typically, in the MHD context of a near-parallel shock, the magnetic field amplification is minimal and smooth, lacking fine structures, which means the electrons are barely affected by it.Electron acceleration presents a significant challenge due to the substantial scale difference from ions, as highlighted by the mass ratio.In the case of quasi-perpendicular shocks, electron acceleration is believed to occur via simple adiabatic reflection due to the large magnetic gradient.However, only a small number of particles participate due to the steep energy requirements for injection (Jebaraj et al. 2023).Conversely, in quasi-parallel shocks, where there is no single, large magnetic gradient, electron acceleration may require the presence of highfrequency oblique wave activity.These waves can then induce either adiabatic or non-adiabatic behavior in electrons.The presence of such waves is illustrated in the schematic of the shock transition presented in Fig. 4B.
The schematic, based on the magnetic field data in Fig. 4C, details various spatial scales and directions for clarity.It demonstrates that at scales where the shock transition is considered a simple MHD discontinuity, it is quasiparallel.At smaller scales, within the shock's transition region, we identify whistler precursors and the ramp-overshoot structure.This area is marked by quasi-perpendicular electromagnetic waves, instrumental in the injection, trapping, and eventual acceleration of electrons to relativistic energies.Such processes have been supported by numerical simulations (Riquelme & Spitkovsky 2011).Once electrons attain sufficient energy and their Larmor radii increase adequately, they can interact with ion-driven waves.Thus, accurately quantifying these distinct spatial scales is critical, as it influences the redistribution of flow kinetic energy into processes like electron heating and acceleration (Balikhin et al. 1993;Schwartz et al. 2011).
Finally, we have demonstrated for the first time that the quasi-perpendicular ramp-overshoot region reflects ions, whereas the quasi-parallel upstream region allows their escape.The latter inability of the quasi-parallel shock to effectively confine particles at the ramp, results in the generation of ion-scale waves in the upstream region.These phenomena collectively create a continuous process at the shock transition, essential for particle injection and sustained acceleration, resulting in an uninterrupted ion population.We have observed this continual presence of ions with energies ranging from less than 100 keV to 10 MeV.We found the concurrent occurrence of >1 MeV streaming ions and intense transverse waves up to 90 minutes prior to shock arrival, suggesting active diffusive shock acceleration.This is supported by the low-energy roll-over observed in the upstream energy spectra, which is unstable leading to the growth of ion-driven waves (Shapiro et al. 1998).Conversely, the confinement of low-energy ions near the shock leads to a consistent E −1 energy spectrum, in line with DSA predictions (with no losses) for a shock with r gas = 4.
In this study, we used in situ evidence to demonstrate that a strong quasi-parallel shock can self-consistently accelerate both ions and electrons from the background medium to high energies.Such as self-consistent mechanism is enabled by the amplification of both the magnitude of the mean magnetic field and the level of multi-scale fluctuations relative to it.This is compelling even from an elementary prediction that a shock must either be very large or highly magnetized to effectively confine high energy particles (Hillas 1984;Terasawa 2011).Our findings naturally lead to the conclusion that strong shocks close to the Sun are highly effective at particle confinement by virtue of being highly magnetized.This makes them comparable to much larger shocks with lesser magnetization.These rare observations were made possible by the pioneering capabilities of the Parker Solar Probe and its unprecedented proximity to the Sun.If the physics of cosmic shocks, like those in SNR blast waves, adhere to the fundamental principles of CSWs, then understanding the evolution of the shock transition is crucial for facilitating the acceleration of charged particles.
experiment was developed and is operated under NASA contract NNN06AA01C.Our study primarily utilizes full electromagnetic fields measured by the FIELDS instrument suite on board the PSP spacecraft (Bale et al. 2016).The electric field measurements are made using the electric fields instrument (EFI) consisting of two pairs of dipole electric field antennas oriented in the TN plane and extending beyond the PSP heat shield, and a fifth antenna located behind the heat shield on the instrument boom; the location of antenna V5 in the wake of PSP means the R component is susceptible to detrimental interference by the wake electric field and cannot be reliably interpreted (Bale et al. 2016).Two three-component flux-gate magnetometers (MAG) measure the magnetic field from DC to approximately 60 Hz during aphelion and up to 293 vector measurements per second during 2 to 4 days around perihelion.The latter is used in the present study.
Additionally, we incorporate high-frequency electric field measurements from both the high-frequency and lowfrequency receivers (HFR & LFR) of the Radio Frequency Spectrometer (RFS; Pulupa et al. 2017).The RFS includes four electric antennas and measures over a wide frequency range, spanning from 1 kHz to 20 MHz at a 3.5 second cadence during the science phase.A particularly relevant application of these measurements in our study is in determining the electron plasma frequency, ω pe .

A.2. Solar Wind Electrons, Alphas, and Protons (SWEAP)
Considering that shock waves manifest as discontinuities in both electromagnetic fields and plasma, our study employed the Solar Wind Electrons, Alphas, and Protons (SWEAP) instrument suite (Kasper et al. 2016), specifically using the ion electrostatic analyzer referred to as the Solar Probe ANalyzer for Ions (SPAN-i) (Livi et al. 2022).SPAN-i measures the 3D velocity distributions of solar wind ions with an additional time-of-flight component to distinguish between protons and heavier ion species.For this event, we used the L3 sf00 data product corresponding to uncontaminated protons at a 3.5s cadence in the energy range 20 eV to ∼20 keV.However, it is important to note a limitation: due to the placement of the instrument on the ram side of PSP, its field of view (FOV) is obstructed by the Thermal Protection System (TPS), which means only partial plasma moments are available.Significant deflections in the solar wind flow (characterized by the angle between the magnetic field vector B and the plasma flow velocity vector v, denoted as θ Bv ), with respect to the spacecraft frame, may result in parts of the bulk distribution being undetected by the instrument.When the proton velocity distribution function (VDF) sufficiently departs from the FOV of SPAN-i, the temperature measurement is largely overestimated since the instrument will only capture the wings of the distribution.
An example of this deflection is evident at 05:30 UT in Fig. 1B, where the VDF remains unrecorded, as deflections in the +T direction are most susceptible to measurement reliability.Additionally, during the shock crossing period, the proton VDF not only sufficiently deflected from the FOV, but also exceeded the energy range of SPAN-i.Consequently, the downstream plasma measurements become entirely artificial and fail to reflect the actual conditions.Regarding the temperature at the time of the shock's arrival, the proton VDF's tails are exceptionally broad, resulting in an anomalously high measured temperature.During the period of time when the upstream velocity was obtained, however, the VDF was sufficiently in the FOV, as verified by plotting the VDFs from the L2 sf00 data product and confirming that the majority of the bulk flow is unobstructed (not shown).However, parts of the VDF other than the core may be unaccounted for when estimating the number density resulting in discrepancies.

A.3. Integrated Science Investigation of the Sun (IS⊙IS)
To analyze the energetic particle populations, in particular electrons and protons, we used the Integrated Science Investigation of the Sun, IS⊙IS instrument suite (McComas et al. 2016).It measures energetic particles from ∼20 keV to over 100 MeV/nuc with two Energetic Particle Instruments (EPI), EPI-Lo (Hill et al. 2017) and EPI-Hi (Wiedenbeck et al. 2017).EPI-Lo is a time-of-flight mass spectrometer with 80 separate apertures to provide a ∼2π-wide FOV.We use a few different data products from EPI-Lo, namely, the ChanP (triple coincidence protons), ChanR (high time cadence triple coincidence proton measurements), ChanT (time-of-flight-only ion measurements), and ChanE (used primarily to measure electrons) channels.EPI-Hi includes two low energy telescopes (LETs) measuring ions from ∼1 to 20 MeV/nuc; one double-ended with a sunward facing aperture, LETA, and anti-sunward aperture, LETB, and a singleended telescope, LETC, pointing orthogonally to the LETA instrument axis.The higher energies (>10 MeV/nuc) are captured in EPI-Hi by a double-ended high energy telescope (HET) with one side, HETA, pointed sunward and the other, HETB, anti-sunward.The sunward and anti-sunward designations are based on the assumption of the nominal magnetic field line.The downstream distribution of the bulk plasma flow is beyond the measurement capabilities of the SPAN-i instrument for E15, which corresponds to approximately 20 keV or 2000 km s −1 .Consequently, it was not feasible to directly obtain plasma parameters such as density and velocity at the moment of the shock crossing, parameters that are crucial for understanding the shock wave's characteristics.
To gauge the extent to which the downstream bulk flow was shifted, we utilized data from the EPI-Lo ChanT, which extends down to 30 keV.This is illustrated in Fig. 1D of the main text, where part of the distribution is observable at 30 keV.This observation strengthens the assertion that the bulk flow's kinetic energy exceeded SPANi's measurement range and was not merely deflected away from its sensor.However, it is a complex task to match a thermal distribution across two different instruments.Therefore, we infer that the center of the distribution lies somewhere between 20 and 30 keV.Since these measurements represent the kinetic energy (E ke ) of the bulk flow, we estimated the downstream velocity (v d ) using the formula v d = 2E ke /m p , where m p denotes the proton mass.This calculation yields a downstream velocity of approximately 2200 ± 200 km s −1 .The ±200 km s −1 uncertainty is due to the large energy bin width (20 to 30 keV) that was considered.Comparatively, the upstream bulk flow velocity (v u ), measured by the spacecraft, was around 410 ± 40 km s −1 in the 2-minute interval preceding the shock's arrival.While it is possible to get an approximation for the magnitude, v d of the downstream flow vector, it is not possible to obtain the direction of the flow.

B.2. Estimation of electron density
To verify the plasma density measurements from SPAN-i, we employed the electron plasma frequency measured by the FIELDS/RFS electric field antennas as a proxy.Under usual conditions, the RFS antennas determine electron density from quasi-thermal noise (Meyer-Vernet et al. 2017).However, this approach was not feasible during the observed period because the Debye length (λ De ) was larger than the effective length of the antennas (Moncuquet et al. 2020).In such cases, electrostatic waves, and especially electron beam-driven Langmuir waves typically resonating near the electron plasma frequency (ω pe ) can be used.In cases where electron sound waves are present, a result of electron cyclotron instabilities, the frequency can be downshifted to approximately 30-60% of ω pe (i.e., 0.4 − 0.7ω pe ) (Lobzin et al. 2005).Therefore, the highest frequency waves detected in our observations are presumed to be Langmuir waves.
When observed with a dipole antenna pair like the FIELDS/RFS (namely V 1 − V 2 and V 3 − V 4 ), electrostatic waves should exhibit high phase coherence (Krasnoselskikh et al. 2011).We constructed a phase coherence spectrum to identify such coherent signals both upstream and downstream of the shock.Coherence can be estimated as, .
In the above equations, X0 X * 0 is the auto spectrum from the (Channel 0) ⃗ X * 1 is the cross spectrum between the two channels.The top panel of Fig. 5 displays this coherence spectrum, clearly distinguishing upstream electrostatic waves near the 60 ± 10 kHz frequency channels.We focused exclusively on highly coherent signals, enabling distinct identification of the waves at the shock transition and the upstream.Non-stationary processes and general noise typically exhibit near-zero phase coherency.Intriguingly, two coherent downstream signals were also detected, marked by arrows, at frequencies of 125 ± 10 kHz.A notable capability of the FIELDS/RFS receivers is their measurement of full Stokes parameters, providing robust verification of the wave types observed.
Langmuir waves, being highly field-aligned, are strongly linearly polarized.Linear polarization is quantified by the Stokes Q parameter, normalized with the total intensity (Stokes I), expressed as Q/I.Both Stokes Q and I are estimated as, The second panel of Fig. 5 presents these strongly linearly polarized electrostatic waves upstream, with the coherent downstream signals also showing a relatively high degree of linear polarization.
During the upstream interval from 07:06 to 07:12 UT, the highest frequency recorded was 65 ± 10 kHz.This uncertainty is due to the frequency channel's 4.5% bandwidth and the variation in the highest frequency wave observed during this period.For the downstream interval, we utilized the frequencies of the two strongly coherent and linearly polarized signals, 125 ± 10 kHz.In situ shock arrival (07:13:10 UT) To estimate the electron plasma density using Langmuir waves, we can use the linear dispersion relation for Langmuir waves in the spacecraft frame: where k is the wavenumber, v th = 2k B T e /m e represents the electron thermal speed, and θ kB is the angle between the solar wind flow and the wave vector k.The equation includes a final term accounting for the Doppler shift caused by the plasma flow past the spacecraft.These two factors may cause ω to deviate from ω pe : the increase in frequency due to the thermal speed and the Doppler shift.It is noteworthy to mention that the Doppler shift may be both positive and negative and accounts for the largest contribution to the observed frequency.The range of k was estimated by Graham et al. (2021), who found 0.11 ≤ kλ De ≤ 0.25 for typical parameters found in the solar wind.For instance, they considered an electron temperature T e ∼ 15 ± 5 eV, similar to the values obtained by Liu et al. (2023) using the PSP measurements, T e ∼ 20 ± 5 eV .The Doppler shift accounts for a solar wind speed V sw ∼ 400 km/s, which is similar to the v u = V sw measured upstream of the shock.Graham et al. (2021) also considered linearly polarized Langmuir waves, which are field-aligned.By using these parameters in the linear dispersion relation, we estimated the Langmuir wave frequency to range from 0.95ω pe to 1.15ω pe .As expected, the variation in ω is largely attributable to the Doppler effect.It is worth noting that this estimation might change when considering a more accurate but complex dispersion relation, such as the quasi-electrostatic slow extraordinary mode wave presented in Larosa et al. (2022).The longitudinal component of the wave, indicated by the ratio , becomes prominent when k is on the order of this ratio.However, for the conditions presented here, where The ω pe is of particular interest, as it is directly related to the electron plasma density (n e ) which is given by the formula: where, e is the elementary charge, ϵ 0 is the vacuum permittivity, and m e is the electron mass.The n e can be estimated by rearranging the formula for ω pe : Given the range of ω with respect to ω pe , we can directly estimate a range of values for ω pe .We estimate the ω pe corresponding to ω = 65 ± 10 kHz is ω pe = 61.5 ± 10 kHz, and for ω ∼ 125 ± 10 kHz we estimate ω pe = 121 ± 12 kHz.Next, by applying the measured values of ω pe , we can estimate n e .The upstream density using ω pe = 61.5 ± 10 kHz results in an estimated average electron plasma density of 46.5 ± 15 cm −3 .However, data from SPAN-i, as shown in the last panel of Fig. 5, indicate an average proton density of approximately 22 cm −3 , which is a factor of two smaller than our estimate.While, ω pe = 121 ± 12 kHz corresponds to an average electron density of 184 ± 35 cm −3 .We may consider the SPAN-i measurement of number density to be a lower limit in our estimations.

C.1. Shock normal
Identifying the shock wave's normal direction necessitates pinpointing its transition.To achieve this, we implemented a moving average technique with a 30-second sliding window to filter out high-frequency signals and isolate the magnetic gradient.We observed an increase in the mean field 30 seconds prior to the transition, indicated by the principal shock jump at 07:13 UT.For the estimation of the normal, we employed two methods: minimum variance analysis (MVA; Sonnerup & Scheible 1998), and the magnetic coplanarity theorem (MCT; Colburn & Sonett 1966).Due to uncertainties and the absence of necessary plasma measurements such as velocity, we could not apply other methods such as velocity coplanarity and mixed-mode techniques, usually regarded as more robust (Paschmann & Schwartz 2000).
For the MVA technique, applied when a spacecraft encounters a transition layer like a shock front, the relevant equations are: Here, M µν equals ⟨B µ B ν ⟩ − ⟨B µ ⟩⟨B ν ⟩, and nMVA aligns with the smallest eigenvalue, λ min .We assess reliability using the ratio λ int /λ max and B n /B, with λ int /λ max ≥ 2 and B n /B ≤ 0.3 indicating a well-defined shock normal.
We selected the time window for MVA to span from 07:03 to 07:23 UT, i.e. centered on the shock ramp at 07:13 UT.Within this interval, the eigenvalue ratio between the intermediate and minimum variance directions, λ int /λ min , was approximately 7.5, suggesting a good distinction-the three eigenvalues are λ max = 8298, λ int = 3070, and λ min = 408.The normal, indicated by the direction of the minimum variance eigenvector, was determined to be [0.976,−0.052, −0.216] in the inertial-RTN coordinate system.However, the magnetic field ratio B n /B ≈ 0.66 suggested that the normal direction was not well-defined, even after accounting for the high-frequency components of the field fluctuations.For a well-defined normal direction, it is expected that B n /B < 0.1.However, in the case of quasi-parallel shocks, where B n is on the order of B, the ratio of eigenvalues (λ int /λ min ) serves as a far better indicator of quality (Paschmann & Schwartz 2000).
To corroborate these findings, we applied the MCT, the only other feasible method using solely magnetic field data.The MCT calculates the shock normal nMCT , assuming B u , B d , and n are coplanar: We set the upstream and downstream windows for this analysis from 07:04 to 07:12 UT and from 07:15 to 07:23 UT, respectively.The choice of the interval is motivated by previous statistical studies of in situ shocks (Kilpua et al. 2015).The estimated shock normal n was found to be [0.995,−0.101, −0.026] in the RTN coordinate system.Despite the limitations of both methods, their yielding of similar results lends credibility to the estimated shock normal.

C.2. Shock geometry
We can determine the shock angle, θ, which is the angle between the shock normal n and the upstream magnetic field B u , using the following definition: Our calculations yielded shock angles of approximately θ MVA Bn ∼ 12.5 • for the MVA method and θ MCT Bn ∼ 3.9 • for the MCT method.Hence, both methods suggest that the shock was quasi-parallel.The shock was measured in situ at a distance of approximately 0.23 astronomical units (a.u.), roughly 3 hours and 55 minutes following the solar eruption.The FIELDS/RFS measurements presented in Fig. 6, the only remote sensing observation in this period, served as the basis for estimating the transit speed.The solar eruption is marked by the dispersionless signal in the FIELDS/RFS receiver when the probe was at 0.24 a.u.This timing suggests a ballistic speed for the eruption-driven shock of approximately 2450 km s −1 .While this offers a basic first-order estimate of the expected speed of the piston (CME), it does not provide an accurate measure of the local shock speed.Knowing the shock normal vector n is essential for estimating the local speed of the shock, a critical parameter for analyzing different shock rest frames.We utilize the upstream velocity (v u = 410 ± 40 km s −1 ) and the downstream velocity (v d = 2200 ± 200 km s −1 , along with the shock normal and the expected density compression ratio of a strong shock (r gas ∼ 4), to calculate the shock speed.The caveat of not having proper plasma measurements downstream is that, it is not possible to obtain the vector direction therefore, only the magnitude of v d is used.To estimate the density compression ratio, we use the proxy measurements of density from the FIELDS/RFS: n d = 184 ± 35 cm −3 for downstream density and n u = 46.5 ± 15 cm −3 for upstream density.From these values, we derive the density compression ratio r gas = n d /n u , which approximates to 4, aligning with the characteristics of a strong shock.We then use the conservation of mass flux across the parallel shock to estimate the shock speed (v sh ) as, gas yielding an estimated speed of v sh ∼ 2800 ± 300 km s −1 which seems to align with the ballistic estimation.The direction of the vector dot product is completely artificial as only the magnitude of the downstream flow (v d ) is used in the calculation.Then we estimate the Alfvénic Mach number (M A ) as the ratio of the upstream velocity component normal to the shock (V u ) to the Alfvén speed (v A ).The Alfvén speed is calculated as v A = B u / √ µ 0 n u m p , where B u is the upstream magnetic field magnitude, and n u is the upstream proton number density.Here, we use the electron number density obtained from proxies.We then compute Consequently, for this range of V u , M A is estimated to be approximately 9.1 ± 1.35.This estimation can be further validated using the proxies established in (Gedalin et al. 2021).We can estimate the value of M A using the expression: Here, B max = 330 nT and B u = 90 ± 5 nT are the overshoot and the upstream magnetic field, respectively.The normalized potential jump (s) is calculated as s = 2ϕ NIF /m p v 2 u , where ϕ NIF is the cross-shock electrostatic potential, m p is the proton mass, and v u is the upstream velocity.The value for s = 0.15 is selected as a reasonable adjustment of the median cross-shock potential to the magnetic overshoot.This selection is justified based on Dimmock et al. (2012), who demonstrate that for M A ∼ 9, s is approximately 0.25, and this value decreases with decreasing θ Bn .Therefore, a value of s = 0.15 for a nearly parallel shock is justified.Inserting these values into the equation yields M A ∼ 9, which closely aligns with the value estimated from observations.
In a magnetized medium with β ≪ 1, the critical Mach number (M c ) for near-parallel shocks can be as low as ∼ 1.5 Kennel et al. (1985).The ratio M A /M c ∼ 6.07 ± 0.9 clearly indicates that the shock is significantly supercritical.However, energetic particles and associated waves may lead to a change in M c (Laming 2022)

D. ENERGETIC PARTICLE OBSERVATIONS
During the event, and particularly at the arrival of the interplanetary (IP) shock, particle intensities increased dramatically.Due to these high intensities, EPI-Hi transitioned to dynamic threshold (DT) modes, which involve raising the energy thresholds on various detector segments to limit the effective instrumental geometry factor for incident protons, helium, and electrons.EPI-Lo, which was less affected by the high intensities, was therefore the primary observing instrument.
In response to these conditions, we utilized all available observations from EPI-Lo.Initially, we employed EPI-Lo's full time-of-flight data (ChanT), as shown in Fig. 1D of the paper, to demonstrate that the bulk flow reached energies within EPI-Lo's range.As ChanT only utilizes particle time-of-flight through the instrument, this channel cannot distinguish ion species.While Fig. 1D displays measurements from only the sunward detector, we supplemented this data with measurements from all wedges in Fig. 7.
Furthermore, we used the energy resolution of EPI-Lo ChanP and the temporal resolution of ChanR to analyze the various proton populations at the shock.Measurements from the sunward pointing detector (W3) are presented in Fig. 2C of the main text.Complementary measurements from all wedges, supporting the discussions in the main text, are provided in Figs. 8 and 9.

E.2. Fourier analysis
In Fig. 2F, we present power spectral densities of the magnetic field in the window 06:40-07:13 UT and in a background window 01:30-02:07.We note that different periods prior to the solar eruption yielded similar background spectra.The background window was chosen to be before the solar eruption at 03:16 UT.This choice is because of the probe's proximity to the Sun, meaning that the fastest particles arrive almost instantaneously in its vicinity.To avoid any potential modifications of the magnetic field influencing our measurements, we selected the pre-event background.Figure 10 shows magnetic field observation during these periods.We performed a Fourier transform (applying a Hann window) of each magnetic field component B R , B T , and B N .The magnetic field in both intervals is strongly dominated by B R .Therefore, we show the B R power spectral density, which estimates the power parallel to the background magnetic field, and the sum of B T and B N power spectral densities, which estimates the total transverse power.The presented spectra have been smoothed by averaging over adjacent frequencies.
There are two important considerations regarding the data analysis conducted here.The first concerns the window sizes, set at approximately 35 minutes, which correspond to a maximum possible wave period of 2100 seconds (approximately 5 × 10 −4 Hz).A wave with a period of 300 seconds (approximately 3 × 10 −3 Hz) would be sampled around seven times.Therefore, for frequencies below 10 −2 Hz or wave periods shorter than 100 seconds (sampled about 21 times), the statistical representation might not accurately reflect reality.The second consideration pertains to the time series analysis itself, which represents a random process-in this case, the magnetic field.The estimations made here are under the assumption of ergodicity, and long time averaging together with the azimuthal motion of PSP may break this behavior.This limitation becomes more significant when PSP is closer to the Sun, affecting the range of wave periods from the same plasma that can be analyzed using Fourier analysis.

E.3. Estimation of characteristic frequencies
In order to understand the effect of the upstream magnetic fluctuations on the ions, we estimate the frequencies of Alfvén waves that resonate with 1 MeV and 5 MeV protons.When considering waves moving parallel to the background magnetic field, wave-particle cyclotron resonance is governed by the Doppler condition (Vainio 2000), Time (UT) where ω is the resonant wave frequency, k ∥ is the wave number, v ∥ is the particle speed along the magnetic field, Ω is the proton cyclotron frequency, and n is an integer.We consider particles with pitch-angle cosine µ = 1, and thus v ∥ = vµ = v, where v is the speed of the particle.We take n = 1 corresponding to the fundamental resonance and consider outward-propagating Alfvén waves with the dispersion relation ω = k ∥ (v u + v A ) (in the spacecraft frame), where v u is the solar wind speed and v A is the Alfvén speed.Here, we ignore the motion of the spacecraft and approximate that the solar wind flows along the magnetic field.Finally, the resonant frequency becomes F. WAVE ANALYSIS OF THE SHOCK RAMP

F.1. Morlet wavelet analysis
To determine the properties of the electromagnetic waves in our observations, we applied a Morlet wavelet transform on the time series of both magnetic and electric fields.Wavelets are generally preferred for non-stationary signals such as whistlers and therefore it allows the estimation of power distribution as a function of time and frequency.In short, it reveals the temporal evolution of spectral parameters of wave activity (Bendat & Piersol 2011).Here, we have followed the same wavelet parameters and methodology as used by Sundkvist et al. (2012) and Karbashewski et al. (2023).

F.2. Whistler wave analysis
An important aspect of any quasi-parallel shock is the number of large amplitude magnetic structures throughout the magnetic gradient of the shock.It is widely understood that the general structure of the quasi-parallel shock is highly non-stationary and therefore time-dependent.At a shock propagating at high speeds such as the one studied here, we would require a close to 1 kHz resolution for both magnetic fields and plasma properties to truly understand the evolution in detail.In our observations, the FIELDS instrument suite was the only one capable of capturing phenomena in the scale of several ion-inertial lengths (λ i = c/ω pi ).We performed a spectral and phase wave analysis of the data which revealed several features in the vicinity the the shock.The large amplitude waves in the frequency range of 0.1-10 Hz are located close to the ramp, both upstream and downstream.The observed local peak of the electric field indicates the crossing of the ramp at 07:13:11.45.The ion-inertial length estimation is based on electron density processing using the FIELDS plasma waves measurement.In the upstream, we found the n e to be ∼ 46.5 ± 15 cm −3 .Assuming that the plasma is quasi-neutral, the ion density may not vary far from the electron density, and as such would correspond to an ion-inertial length, λ i ∼ 32.5 ± 3 km.
The distinct whistler precursor is observed in the upstream region of the quasi-perpendicular ramp.Several wave periods are seen at frequencies around 30 Hz.The waves are right-hand circularly polarized (Fig. 3H) and have amplitudes up to 35 nT.The wave normal angles (WNA, θ kB ) are derived by making use of the singular value decomposition (SVD) technique (Karbashewski et al. 2023).The WNA for the whistler wave was found to be around 73 ± 3 • (Fig. 3E).The oblique propagation is the essential factor for the wave precursor to be a standing wave in the shock frame.
The Poynting flux can only reliably be determined along the Sun-spacecraft R-axis due to the unreliable R-component electric field measurement.The Poynting flux in the spacecraft frame is calculated using the complex ⃗ E w and ⃗ B w spectra: S z = Re(E wx • B * wy − E wy • B * wx ) from Fig. 3(C,D) and subsequently mapped to a bi-symmetric logarithmic scale shown in Fig. 3F.The wave phase speed (v ph ) is estimated using the flow speed and the WNA.

F.2.1. Schematic of the Shock transition
The schematic in Fig. 4B presents a 2D simplified reconstruction of the shock structure based on the values and direction of the magnetic field recorded by PSP during the shock crossing presented in Fig. 4C.The changing velocity of the plasma flow is reflected in the nonlinear (but isotropic in each particular point of the schematic) spatial scales.The physical scales are indicated by the squares of 50λ i × 50λ i .The color scheme represents the magnetic field magnitude.

Figure 1 .
Figure 1.Magnetic field and plasma measurements during the shock encounter.Panels A and B provide a 9-hour overview of the magnetic field and the one-dimensional proton energy flux, capturing the pre-event solar wind, the in situ shock passage, the following sheath region, and the initial portion of the associated CME ejecta.The SPAN-i data downstream of the shock do not accurately represent reality, as the flow is largely deflected away from the detector's field of view.Panels C to E offer an 80-minute close-up leading to the shock arrival, showing variations in the bulk ion distribution as captured by the SPAN-i and EPI-Lo instruments.Finally, Panels F to H focus on a 2-minute snapshot close to the shock transition, emphasizing fluctuations in the total magnetic field as well as the change in proton temperature.The downstream proton temperature is marked with hatching due to significant instrumental uncertainties.

Figure 2 .
Figure 2. Energetic particles and the evolution of the wave foreshock.Panels A-D show the time-energy plots for high-energy electron counts detected by EPI-Hi, low-energy electron flux from EPI-Lo, proton flux from EPI-Lo, and the amplitude of transversal (⊥) and parallel (∥) magnetic fluctuations.In panel E, we present the proton intensity spectra for three selected times (06:40, 07:10, and at the shock, averaged as 07:13-07:15 UT) with a E −1 power-law shown for reference, which indicates the changes in energetic particle distribution.Panel F shows the Fourier power spectra of BR (estimating B ∥ power) and the sum of BT and BN (estimating B ⊥ power) in the foreshock and the background.Variations in the proton flux and magnetic fluctuations over a two-minute period surrounding the shock transition are shown in panels G and H.

Figure 3 .
Figure 3. Analysis of the waves at the shock transition.Panels A and B show the magnetic and electric field components.The duration of the precursor whistlers (between 07:13:11.10-07:13:11.45) is annotated on top.Panel A also indicates the approximate spatial scale of the upstream in units of λi estimated in the plasma frame.Panels C and D show the Morlet wavelet spectrum of the magnetic and electric fields.The θ kB angle between the wave vector k and the magnetic field B is shown in panel E. In panel F, we show the Poynting flux estimated in the spacecraft frame S sc .Finally, the hodographs of precursor whistler waves in the maximum vs intermediate variance plane and the maximum vs minimum variance planes are presented in panels G and H, respectively.In both panels, the blue arrow denotes the mean field B0.

Figure 4 .
Figure 4. Transition region of the quasi-parallel shock.Panel A displays the time evolution (depicted by the color blue to red) of the transversal magnetic field components.The dotted line through the middle represents the mean field, B0.Panel B illustrates a two-dimensional schematic of the shock transition in the plasma frame.The spatial scales are in λi (represented with 50λi×50λi squares) and are estimated for the upstream and downstream separated at the edge of the quasi-parallel shock (Q ∥ ).The grayscale represents the change in B across the shock.The direction of B and the shock normal are also marked by arrows.The dark stripes in the middle represent the quasi-perpendicular structure at the ramp and the whistler precursors.Panel C presents the magnetic field magnitude on which the schematic is based upon.
Estimation of the downstream flow speed

Figure 5 .
Figure 5.In situ arrival of the shock as observed through high-frequency electric field data from the FIELDS/RFS.The top panel displays the coherence spectra, while the middle panel presents the normalized linear polarization, Stokes Q/I.The estimated upstream and downstream plasma frequencies (ωpe) are indicated by horizontal black lines, with arrows highlighting the observed electrostatic waves.The bottom panel juxtaposes the estimated electron density (ne) -with the range of uncertainty depicted as a red shaded area -against the ion density measured by SPAN-i.

CFigure 6 .
Figure 6.Dynamic radio spectra constructed from the FIELDS/RFS receivers.The Solar flare onset, and the in situ shock arrival are indicated on top of the panel.
Figure 7. Ion Energy Fluxes Measured by EPI-Lo ChanT, the full time-of-flight data.The vertical line indicates the moment of the in situ shock's arrival.Grey-shaded areas in the graph represent intervals where data was not available.

Figure 10 .
Figure 10.Window selection for Fourier analysis.The pink shaded areas indicate the FFT (Fast Fourier Transform) windows utilized for estimating the power spectral density shown in Fig. 2F.

Table 1 .
List of relevant plasma and shock parameters associated with the March 13, 2023 event.