Fast Particle Acceleration in 3D Hybrid Simulations of Quasiperpendicular Shocks.

Understanding the conditions conducive to particle acceleration at collisionless, nonrelativistic shocks is important for the origin of cosmic rays. We use hybrid (kinetic ions-fluid electrons) kinetic simulations to investigate particle acceleration and magnetic field amplification at nonrelativistic, weakly magnetized, quasiperpendicular shocks. So far, no self-consistent kinetic simulation has reported nonthermal tails at quasiperpendicular shocks. Unlike 2D simulations, 3D runs show that protons develop a nonthermal tail spontaneously (i.e., from the thermal bath and without preexisting magnetic turbulence). They are rapidly accelerated via shock drift acceleration up to a maximum energy determined by their escape upstream. We discuss the implications of our results for the phenomenology of heliospheric shocks, supernova remnants, and radio supernovae.

Epidemic Management via Imperfect Testing: A Multi-criterial Perspective.

Diagnostic testing may represent a key component in response to an ongoing epidemic, especially if coupled with containment measures, such as mandatory self-isolation, aimed to prevent infectious individuals from furthering onward transmission while allowing non-infected individuals to go about their lives. However, by its own nature as an imperfect binary classifier, testing can produce false negative or false positive results. Both types of misclassification are problematic: while the former may exacerbate the spread of disease, the latter may result in unnecessary isolation mandates and socioeconomic burden. As clearly shown by the COVID-19 pandemic, achieving adequate protection for both people and society is a crucial, yet highly challenging task that needs to be addressed in managing large-scale epidemic transmission. To explore the trade-offs imposed by diagnostic testing and mandatory isolation as tools for epidemic containment, here we present an extension of the classical Susceptible-Infected-Recovered model that accounts for an additional stratification of the population based on the results of diagnostic testing. We show that, under suitable epidemiological conditions, a careful assessment of testing and isolation protocols can contribute to epidemic containment, even in the presence of false negative/positive results. Also, using a multi-criterial framework, we identify simple, yet Pareto-efficient testing and isolation scenarios that can minimize case count, isolation time, or seek a trade-off solution for these often contrasting epidemic management objectives.

Relativistic Asymmetric Magnetic Reconnection.

We derive basic scaling equations for relativistic magnetic reconnection in the general case of asymmetric inflow conditions and obtain predictions for the outflow Lorentz factor and the reconnection rate. Kinetic particle-in-cell simulations show that the outflow speeds as well as the nonthermal spectral index are constrained by the inflowing plasma with the weaker magnetic energy per particle, in agreement with the scaling predictions. These results are significant for understanding nonthermal emission from reconnection in magnetically dominated astrophysical systems, many of which may be asymmetric in nature. The results provide a quantitative approach for including asymmetry on reconnection in the relativistic regime.

Spectrum of Electrons Accelerated in Supernova Remnants.

Using a semianalytic model of nonlinear diffusive shock acceleration, we model the total spectrum of cosmic ray (CR) electrons accelerated by supernova remnants (SNRs). Because electrons experience synchrotron losses in the amplified magnetic fields characteristic of SNRs, they exhibit substantially steeper spectra than protons. In particular, we find that the difference between the electron and proton spectral index (power law slope) ranges from 0.1 to 0.4. Our findings must be reckoned with theories of Galactic CR transport, which often assume that electrons and protons are injected with the same slope, and may especially have implications for the observed "positron excess".

Effect of Cosmic Rays on the Evolution and Momentum Deposition of Supernova Remnants.

Using a semianalytical approach based on the thin-shell approximation, we calculate the long-term evolution of supernova remnants (SNRs) while also accounting for the cosmic rays (CRs) accelerated at their blast waves. Our solution reproduces the results of state-of-the-art fluid simulations across the adiabatic and radiative stages for the gas-only case, and it predicts that typical CR acceleration efficiencies (≈10%) can boost SNR momentum deposition by a factor of 2-3. This enhancement can become as large as an order of magnitude in environments in which the gas experiences more severe radiative losses. This result may have a crucial impact on modeling the effect of supernova feedback on star formation and galaxy evolution.

Chemical Enhancements in Shock-Accelerated Particles: Ab initio Simulations.

We study the thermalization, injection, and acceleration of ions with different mass/charge ratios, A/Z, in nonrelativistic collisionless shocks via hybrid (kinetic ions-fluid electrons) simulations. In general, ions thermalize to a postshock temperature proportional to A. When diffusive shock acceleration is efficient, ions develop a nonthermal tail whose extent scales with Z and whose normalization is enhanced as (A/Z)^{2} so that incompletely ionized heavy ions are preferentially accelerated. We discuss how these findings can explain observed heavy-ion enhancements in Galactic cosmic rays.

Simultaneous acceleration of protons and electrons at nonrelativistic quasiparallel collisionless shocks.

We study diffusive shock acceleration (DSA) of protons and electrons at nonrelativistic, high Mach number, quasiparallel, collisionless shocks by means of self-consistent 1D particle-in-cell simulations. For the first time, both species are found to develop power-law distributions with the universal spectral index -4 in momentum space, in agreement with the prediction of DSA. We find that scattering of both protons and electrons is mediated by right-handed circularly polarized waves excited by the current of energetic protons via nonresonant hybrid (Bell) instability. Protons are injected into DSA after a few gyrocycles of shock drift acceleration (SDA), while electrons are first preheated via SDA, then energized via a hybrid acceleration process that involves both SDA and Fermi-like acceleration mediated by Bell waves, before eventual injection into DSA. Using the simulations we can measure the electron-proton ratio in accelerated particles, which is of paramount importance for explaining the cosmic ray fluxes measured on Earth and the multiwavelength emission of astrophysical objects such as supernova remnants, radio supernovae, and galaxy clusters. We find the normalization of the electron power law is ≲10^{-2} of the protons for strong nonrelativistic shocks.

The supernova remnant SN 1006 as a Galactic particle accelerator.

The origin of cosmic rays is a pivotal open issue of high-energy astrophysics. Supernova remnants are strong candidates to be the Galactic factory of cosmic rays, their blast waves being powerful particle accelerators. However, supernova remnants can power the observed flux of cosmic rays only if they transfer a significant fraction of their kinetic energy to the accelerated particles, but conclusive evidence for such efficient acceleration is still lacking. In this scenario, the shock energy channeled to cosmic rays should induce a higher post-shock density than that predicted by standard shock conditions. Here we show this effect, and probe its dependence on the orientation of the ambient magnetic field, by analyzing deep X-ray observations of the Galactic remnant of SN 1006. By comparing our results with state-of-the-art models, we conclude that SN 1006 is an efficient source of cosmic rays and obtain an observational support for the quasi-parallel acceleration mechanism.

Electron Energy Partition across Interplanetary Shocks. I. Methodology and Data Product.

Analyses of 15,314 electron velocity distribution functions (VDFs) within ±2 hr of 52 interplanetary (IP) shocks observed by the Wind spacecraft near 1 au are introduced. The electron VDFs are fit to the sum of three model functions for the cold dense core, hot tenuous halo, and field-aligned beam/strahl component. The best results were found by modeling the core as either a bi-kappa or a symmetric (or asymmetric) bi-self-similar VDF, while both the halo and beam/strahl components were best fit to bi-kappa VDF. This is the first statistical study to show that the core electron distribution is better fit to a self-similar VDF than a bi-Maxwellian under all conditions. The self-similar distribution deviation from a Maxwellian is a measure of inelasticity in particle scattering from waves and/or turbulence. The ranges of values defined by the lower and upper quartiles for the kappa exponents are κ ec ~ 5.40-10.2 for the core, κ eh ~ 3.58-5.34 for the halo, and κ eb ~ 3.40-5.16 for the beam/strahl. The lower-to-upper quartile range of symmetric bi-self-similar core exponents is s ec ~ 2.00-2.04, and those of asymmetric bi-self-similar core exponents are p ec ~ 2.20-4.00 for the parallel exponent and q ec ~ 2.00-2.46 for the perpendicular exponent. The nuanced details of the fit procedure and description of resulting data product are also presented. The statistics and detailed analysis of the results are presented in Paper II and Paper III of this three-part study.