Search for: All records

Abstract Collisionless plasma shocks are a common feature of many space and astrophysical systems. They are sources of high-energy particles and nonthermal emission, channeling as much as 20% of the shock’s energy into nonthermal particles. The generation and acceleration of these nonthermal particles have been previously studied and shown to affect shock hydrodynamics to the zeroth order. In this work, we use self-consistent hybrid particle-in-cell simulations to examine the effect of self-generated nonthermal ion populations on the nature of collisionless, quasi-parallel shocks. Accelerated nonthermal particles downstream of the shock diffuse into the upstream region, taking energy away from the shock, which increases the compression ratio, slows the shock down, and flattens the nonthermal population’s spectral index for lower-Mach-number shocks. We show that this enhances shock compressibility when the heat flux is included in the Rankine–Hugoniot jump conditions, results that are roughly consistent with previous theories of “cosmic-ray-modified shocks.” Additionally, the simulation data show that heat flux and enthalpy flux cancels out in the upstream region, yielding a relatively simple, alternative closure for the jump conditions which accurately predict for the shock speed and compression ratio. The results have the potential to explain discrepancies between predictions and observations in a wide range of systems, such as inaccuracies in predictions of the arrival times of coronal mass ejections and the conflicting radio and X-ray observations of intracluster shocks. These effects will likely need to be included in fluid modeling to predict shock evolution accurately. 

Large-amplitude electrostatic fluctuations are routinely observed by spacecraft upon traversal of collisionless shocks in the heliosphere. Kinetic simulations of shocks have struggled to reproduce the amplitude of such fluctuations, complicating efforts to un- derstand their influence on energy dissipation and shock structure. In this paper, 1D particle-in-cell simulations with realistic proton-to-electron mass ratio are used to show that in cases with upstream electron temperature Te exceeding the ion temperature Ti, the magnitude of the fluctuations increases with the electron plasma-to-cyclotron frequency ratio ωpe/Ωce, reaching realistic values at ωpe/Ωce ≳ 30. The large-amplitude fluctuations in the simulations are shown to be associated with electrostatic solitary structures, such as ion phase-space holes. In the cases where upstream temperature ratio is reversed, the magnitude of the fluctuations remains small. 

Abstract Collisionless low-Mach-number shocks are abundant in astrophysical and space plasma environments, exhibiting complex wave activity and wave–particle interactions. In this paper, we present 2D Particle-in-Cell (PIC) simulations of quasi-perpendicular nonrelativistic (v_sh≈ (5500–22000) km s⁻¹) low-Mach-number shocks, with a specific focus on studying electrostatic waves in the shock ramp and precursor regions. In these shocks, an ion-scale oblique whistler wave creates a configuration with two hot counterstreaming electron beams, which drive unstable electron acoustic waves (EAWs) that can turn into electrostatic solitary waves (ESWs) at the late stage of their evolution. By conducting simulations with periodic boundaries, we show that the EAW properties agree with linear dispersion analysis. The characteristics of ESWs in shock simulations, including their wavelength and amplitude, depend on the shock velocity. When extrapolated to shocks with realistic velocities (v_sh≈ 300 km s⁻¹), the ESW wavelength is reduced to one-tenth of the electron skin depth and the ESW amplitude is anticipated to surpass that of the quasi-static electric field by more than a factor of 100. These theoretical predictions may explain a discrepancy, between PIC and satellite measurements, in the relative amplitude of high- and low-frequency electric field fluctuations. 

The shock-drift acceleration of ions at quasi-perpendicular shocks is a well-known kinetic mechanism for the acceleration of a small fraction of incoming ions to high energy. Here, we use a suite of sixteen hybrid simulations of quasi-perpendicular collisionless shocks over the range of Alfvén Mach number 4.3≤MA≤15.8 (corresponding to a range of fast magnetosonic Mach numbers 2.6≤Mf≤9.4) and shock-normal angle 45°≤θBn≤90° to identify the velocity-space signature of shock-drift acceleration using the field-particle correlation technique. We show that the features of the ion velocity distribution in the shock foot and ramp regions can be clearly interpreted by analysis of the single-particle trajectory of a reflected ion through the full 3D-3V phase space. The characteristic features of the velocity-space signature of shock-drift acceleration remain qualitatively robust over the full parameter range of our simulations, providing a potential means for its identification using single-point spacecraft measurements. At higher Alfvén Mach numbers MA≳8 (Mf≳5), kinetic instabilities generate fluctuations of the electromagnetic fields within the shock transition region, leading to fluctuations in and smearing out of the resulting velocity-space signatures, but the signature remains generally robust and identifiable. The results on the shock-drift acceleration of ions presented here represent a novel means to determine more completely the partitioning of upstream bulk flow kinetic energy into plasma heating, particle acceleration, and electromagnetic fields in collisionless shocks. 

Abstract Shocks in collisionless astrophysical plasmas redistribute some of the incident flow energy into both thermal and non‐thermal energy. Quantifying the partition of that energy amongst various particle species or their sub‐populations, and electromagnetic energy, represents a fundamental goal of shock physics. It embodies the role of the equation of state for the system. Here we apply a framework to assess all the incident and downstream energy fluxes at a crossing of Earth's bow shock for which the upstream magnetic field was roughly aligned with the shock normal direction. Such quasi‐parallel shocks are known to be non‐steady and to produce significant populations of suprathermal particles. We quantify the evolution of all the important carriers of energy flux through the shock region. We sub‐divide the proton population into thermal, suprathermal, and energetic components in order to investigate the shock's efficiency in energizing the nonthermal particles. While the largest energy fluxes are found in the incident proton ram energy and downstream proton thermal enthalpy fluxes, a significant suprathermal population pervades the regions both up‐ and downstream. We also evaluate the energy fluxes attributable to fluctuations in the fluid and field parameters. 

The current state of the art thermal particle measurements in the solar wind are insufficient to address many long standing, fundamental physical processes. The solar wind is a weakly collisional ionized gas experiencing collective effects due to long-range electromagnetic forces. Unlike a collisionally mediated fluid like Earth’s atmosphere, the solar wind is not in thermodynamic or thermal equilibrium. For that reason, the solar wind exhibits multiple particle populations for each particle species. We can mostly resolve the three major electron populations (e.g., core, halo, strahl, and superhalo) in the solar wind. For the ions, we can sometimes separate the proton core from a secondary proton beam and heavier ion species like alpha-particles. However, as the solar wind becomes cold or hot, our ability to separate these becomes more difficult. Instrumental limitations have prevented us from properly resolving features within each ion population. This destroys our ability to properly examine energy budgets across transient, discontinuous phenomena (e.g., shock waves) and the evolution of the velocity distribution functions. Herein we illustrate both the limitations of current instrumentation and why higher resolutions are necessary to properly address the fundamental kinetic physics of the solar wind. This is accomplished by directly comparing to some current solar wind observations with calculations of velocity moments to illustrate the inaccuracy and incompleteness of poor resolution data. 

A large number of heliophysicists from across career levels, institution types, and job titles came together to support a poster at Heliophysics 2050 and the position papers for the 2024 Heliophysics decadal survey titled “Cultivating a Culture of Inclusivity in Heliophysics,” “The Importance of Policies: It’s not just a pipeline problem,” and “Mentorship within Heliophysics.” While writing these position papers, the number of people who privately shareddisturbing stories and experiences of bullying and harassmentwas shocking. The number of people who privately expressed howburned outthey were was staggering. The number of people who privately spoke about how theyconsidered leaving the field for their and their family’s healthwas astounding. And for as much good there is in our community, it is still atoxic environmentfor many. If we fail to do something now, our field will continue to suffer. While acknowledging the ongoing growth that we as individuals must work toward, we call on our colleagues to join us in working on organizational, group, and personal levels toward a truly inclusive culture, for the wellbeing of our colleagues and the success of our field. This work includes policies, processes, and commitments to promote:accountabilityfor bad actors;financial securitythrough removing the constant anxiety about funding;prioritizationof mental health and community through removing constant deadlines and constant last-minute requests;a collaborative culturerather than a hyper-competitive one; anda community where people can thrive as whole personsand do not have to give up a healthy or well-rounded life to succeed. 

Using the field–particle correlation technique, we examine the particle energization in a three-dimensional (one spatial dimension and two velocity dimensions; 1D-2V) continuum Vlasov–Maxwell simulation of a perpendicular magnetized collisionless shock. The combination of the field–particle correlation technique with the high-fidelity representation of the particle distribution function provided by a direct discretization of the Vlasov equation allows us to ascertain the details of the exchange of energy between the electromagnetic fields and the particles in phase space. We identify the velocity-space signatures of shock-drift acceleration of the ions and adiabatic heating of the electrons arising from the perpendicular collisionless shock by constructing a simplified model with the minimum ingredients necessary to produce the observed energization signatures in the self-consistent Vlasov–Maxwell simulation. We are thus able to completely characterize the energy transfer in the perpendicular collisionless shock considered here and provide predictions for the application of the field–particle correlation technique to spacecraft measurements of collisionless shocks.