We study diffusive shock acceleration (DSA) of electrons in non-relativistic quasi-perpendicular shocks using self-consistent one-dimensional particle-in-cell (PIC) simulations. By exploring the parameter space of sonic and Alfv\'{e}nic Mach numbers we find that high Mach number quasi-perpendicular shocks can efficiently accelerate electrons to power-law downstream spectra with slopes consistent with DSA prediction. Electrons are reflected by magnetic mirroring at the shock and drive non-resonant waves in the upstream. Reflected electrons are trapped between the shock front and upstream waves and undergo multiple cycles of shock drift acceleration before the injection into DSA. Strong current-driven waves also temporarily change the shock obliquity and cause mild proton pre-acceleration even in quasi-perpendicular shocks, which otherwise do not accelerate protons. These results can be used to understand nonthermal emission in supernova remnants and intracluster medium in galaxy clusters.

We present X-ray and radio observations of the Fast Blue Optical Transient (FBOT) CRTS-CSS161010 J045834-081803 (CSS161010 hereafter) at t=69-531 days. CSS161010 shows luminous X-ray ($L_x\sim5\times 10^{39}\,\rm{erg\,s^{-1}}$) and radio ($L_{\nu}\sim10^{29}\,\rm{erg\,s^{-1}Hz^{-1}}$) emission. The radio emission peaked at ~100 days post transient explosion and rapidly decayed. We interpret these observations in the context of synchrotron emission from an expanding blastwave. CSS161010 launched a mildly relativistic outflow with velocity $\Gamma\beta c\ge0.55c$ at ~100 days. This is faster than the non-relativistic AT2018cow ($\Gamma\beta c\sim0.1c$) and closer to ZTF18abvkwla ($\Gamma\beta c\ge0.3c$ at 63 days). The inferred initial kinetic energy of CSS161010 ($E_k\gtrsim10^{51}$ erg) is comparable to that of long Gamma Ray Bursts (GRBs), but the ejecta mass that is coupled to the mildly relativistic outflow is significantly larger ($\sim0.01-0.1\,\rm{M_{\odot}}$). This is consistent with the lack of observed gamma-rays. The luminous X-rays were produced by a different emission component to the synchrotron radio emission. CSS161010 is located at ~150 Mpc in a dwarf galaxy with stellar mass $M_{*}\sim10^{7}\,\rm{M_{\odot}}$ and specific star formation rate $sSFR\sim 0.3\,\rm{Gyr^{-1}}$. This mass is among the lowest inferred for host-galaxies of explosive transients from massive stars. Our observations of CSS161010 are consistent with an engine-driven aspherical explosion from a rare evolutionary path of a H-rich stellar progenitor, but we cannot rule out a stellar tidal disruption event on a centrally-located intermediate mass black hole. Regardless of the physical mechanism, CSS161010 establishes the existence of a new class of rare (rate $<0.4\%$ of the local core-collapse supernova rate) H-rich transients that can launch mildly relativistic outflows.

We investigate the possibility of secondary production as the source of cosmic ray (CR) positrons, taking into account observed electron and antiproton spectra and allowing that the reported steepening in the AMS-02 positron spectrum may be due to radiative losses. Using a simple argument based on the antiproton spectrum, we show that positrons cannot be purely secondary. Namely, the antiproton to positron ratio requires that positrons be uncooled up to a few hundred GeV. This requirement implies that CR protons escape from the Galaxy too quickly to reproduce the observed antiproton flux. We further show that this result rules out also more complex secondary production scenarios such as reacceleration and the Nested Leaky Box. Thus, we conclude that while antiprotons can be produced exclusively in secondary interactions, a primary source of CR positrons -- be it dark matter or pulsars -- is required to explain their spectrum.

We present the results of a detailed investigation into the physical conditions in interstellar material interacting with the supernova remnant IC 443. Our analysis is based on a comprehensive examination of high-resolution far-ultraviolet spectra obtained with the Space Telescope Imaging Spectrograph onboard the Hubble Space Telescope of two stars behind IC 443. One of our targets (HD 43582) probes gas along the entire line of sight through the supernova remnant, while the other (HD 254755) samples material located ahead of the primary supernova shock front. We identify low velocity quiescent gas in both directions and find that the densities and temperatures in these components are typical of diffuse atomic and molecular clouds. Numerous high velocity components are observed in the absorption profiles of neutral and singly-ionized atomic species toward HD 43582. These components exhibit a combination of greatly enhanced thermal pressures and significantly reduced dust-grain depletions. We interpret this material as cooling gas in a recombination zone far downstream from shocks driven into neutral gas clumps. The pressures derived for a group of ionized gas components at high positive velocity toward HD 43582 are lower than those of the other shocked components, pointing to pressure inhomogeneities across the remnant. A strong very high velocity component near -620 km/s is seen in the absorption profiles of highly-ionized species toward HD 43582. The velocity of this material is consistent with the range of shock velocities implied by observations of soft thermal X-ray emission from IC 443. Moderately high-velocity gas toward HD 254755 may represent shocked material from a separate foreground supernova remnant.

Analysis of model fit results of 15,210 electron velocity distribution functions (VDFs), observed within $\pm$2 hours of 52 interplanetary (IP) shocks by the Wind spacecraft near 1 AU, is presented as the third and final part on electron VDFs near IP shocks. The core electrons and protons dominate in the magnitude and change in the partial-to-total thermal pressure ratio, with the core electrons often gaining as much or more than the protons. Only a moderate positive correlation is observed between the electron temperature and the kinetic energy change across the shock, while weaker, if any, correlations were found with any other macroscopic shock parameter. No VDF parameter correlated with the shock normal angle. The electron VDF evolves from a narrowly peaked core with flaring suprathermal tails in the upstream to either a slightly hotter core with steeper tails or much hotter flattop core with even steeper tails downstream of the weaker and strongest shocks, respectively. Both quasi-static and fluctuating fields are examined as possible mechanisms modifying the VDF but neither is sufficient alone. For instance, flattop VDFs can be generated by nonlinear ion acoustic wave stochastic acceleration (i.e., inelastic collisions) while other work suggested they result from the combination of quasi-static and fluctuating fields. This three-part study shows that not only are these systems not thermodynamic in nature, even kinetic models may require modification to include things like inelastic collision operators to properly model electron VDF evolution across shocks or in the solar wind.

It has been proposed that Ultra-High-Energy Cosmic Rays (UHECRs) up to $10^{20}$eV could be produced in the relativistic jets of powerful Active Galactic Nuclei (AGNs) via a one-shot reacceleration of lower-energy CRs produced in supernova remnants (the $\textit{espresso}$ mechanism). We test this theory by propagating particles in realistic 3D magneto-hydrodynamic simulations of ultra-relativistic jets and find that about $10\%$ of the CRs entering the jet are boosted by at least a factor of $\sim\Gamma^2$ in energy, where $\Gamma$ is the jet's effective Lorentz factor, in agreement with the analytical expectations. Furthermore, about $0.1\%$ of the CRs undergo two or more shots and achieve boosts well in excess of $\Gamma^2$. Particles are typically accelerated up to the Hillas limit, suggesting that the $\textit{espresso}$ mechanism may promote galactic-like CRs to UHECRs even in AGN jets with moderate Lorentz factors, and not in powerful blazars only. Finally, we find that the sign of the toroidal magnetic field in the jet and in the cocoon controls the angular distribution of the reaccelerated particles, leading to a UHECR emission that may be either quasi-isotropic or beamed along the jet axis. These findings strongly support the idea that $\textit{espresso}$ acceleration in AGN jets can account for the UHECR spectra, chemical composition, and arrival directions measured by Auger and Telescope Array.

We present the first plasma simulations obtained with the code dHybridR, a hybrid particle-in-cell code with fluid electrons and both thermal and energetic ions that retain relativistic dynamics. dHybridR is constructed to study astrophysical and space-physics problems where a few energetic non-thermal particles (i.e., cosmic rays, CRs) affect the overall dynamics of a non-relativistic plasma, such as CR-driven instabilities, collisionless shocks, magnetic reconnection, turbulence, etc.In this method paper we provide some applications to linear (resonant/non-resonant CR streaming instability) and strongly non-linear (parallel shocks) problems that show the capabilities of the code. In particular, we provide the first self-consistent hybrid runs that show the acceleration of relativistic ions at non-relativistic shocks; CRs develop a power-law in momentum, which translates to a broken power law in energy that exhibits a steepening around the ion rest mass, as predicted by the theory of diffusive shock acceleration. We present examples of 2D dHybridR runs relevant for fast shocks in radio supernovae, whose evolution can be followed in real time, and 3D runs of low-Mach-number heliospheric shocks, which can be compared with in-situ spacecraft observations.

A statistical analysis of 15,210 electron velocity distribution function (VDF) fits, observed within $\pm$2 hours of 52 interplanetary (IP) shocks by the $Wind$ spacecraft near 1 AU, is presented. This is the second in a three-part series on electron VDFs near IP shocks. The electron velocity moment statistics for the dense, low energy core, tenuous, hot halo, and field-aligned beam/strahl are a statistically significant list of values illustrated with both histograms and tabular lists for reference and baselines in future work. The beam/strahl fit results in the upstream are currently the closest thing to a proper parameterization of the beam/strahl electron velocity moments in the ambient solar wind. This work will also serve as a 1 AU baseline and reference for missions like $Parker \ Solar \ Probe$ and $Solar \ Orbiter$. The median density, temperature, beta, and temperature anisotropy values for the core(halo)[beam/strahl] components, with subscripts $ec$($eh$)[$eb$], of all fit results respectively are $n{\scriptstyle_{ec(h)[b]}}$ $\sim$ 11.3(0.36)[0.17] $cm^{-3}$, $T{\scriptstyle_{ec(h)[b], tot}}$ $\sim$ 14.6(48.4)[40.2] $eV$, $\beta{\scriptstyle_{ec(h)[b], tot}}$ $\sim$ 0.93(0.11)[0.05], and $\mathcal{A}{\scriptstyle_{ec(h)[b]}}$ $\sim$ 0.98(1.03)[0.93]. The nuanced details of the fitting method and data product description were published in Paper I and the detailed analysis of the results will be shown in Paper III.

Using hybrid simulations (kinetic ions--fluid electrons), we test the linear theory predictions of the cosmic ray (CR) streaming instability. We consider two types of CR distribution functions: a "hot" distribution where CRs are represented by a drifting power law in momentum and an anisotropic "beam" of monochromatic particles. Additionally, for each CR distribution we scan over different CR densities to transition from triggering the resonant to the non-resonant (Bell) streaming instability. We determine the growth rates of these instabilities in simulations by fitting an exponential curve during the linear stage, and we show that they agree well with the theoretical predictions as a function of wave number agree. We also examine the magnetic helicity as a function of time and wave number, finding a general good agreement with the predictions, as well as some unexpected non-linear features to the instability development.

The nonresonant cosmic ray instability, predicted by Bell (2004), is thought to play an important role in the acceleration and confinement of cosmic rays (CRs) close to supernova remnants. Despite its importance, the exact mechanism responsible for the saturation of the instability has not been determined, and there is no first-principle prediction for the amplitude of the saturated magnetic field. Using a survey of self-consistent kinetic hybrid simulations (with kinetic ions and fluid electrons), we study the saturation of the non-resonant streaming instability as a function of the parameters of both the thermal background plasma and the CR population. The strength of the saturated magnetic field has important implications for both CR acceleration in supernova remnants and CR diffusion in the Galaxy.

Espresso is a novel acceleration model for Ultra-High-Energy Cosmic Rays (UHECRs), where lower-energy CRs produced in supernova remnants experience a one-shot reacceleration in the relativistic jets of powerful Active Galactic Nuclei (AGNs) to reach energies up to $10^{20}$ eV. To test the espresso framework, we follow UHECR acceleration bottom-up from injection to the highest energies by propagating 100,000 particles in realistic 3D magneto-hydrodynamic (MHD) simulations of ultra-relativistic jets. We find that simulations agree well with analytical expectations in terms of trajectories of individual particles. We also quantify that $\sim 10\%$ of CR seeds gain a factor of $\sim\Gamma^2$ in energy, where $\Gamma$ is the jet's effective Lorentz factor; moreover, about $0.1\%$ of the particles undergo two or more shots to achieve gains in excess of $\Gamma^2$. Particles are generally accelerated up to the jet's Hillas limit, indicating that the espresso mechanism should boost galactic CRs to UHECRs in typical AGN jets. Finally, we discuss how espresso acceleration in AGN jets is consistent with UHECR spectra and chemical composition, and also with the UHECR arrival directions measured by Auger and Telescope Array.

We discuss the recent developments in the theory of diffusive shock acceleration (DSA) by using both first-principle kinetic plasma simulations and analytical theory based on the solution of the convection/diffusion equation. In particular, we show how simulations reveal that the spectra of accelerated particles are significantly steeper than the $E^{-2}$ predicted by the standard theory of DSA for strong shocks, in agreement with several observational pieces of evidence. We single out which standard assumptions of test-particle and non-linear DSA are violated in the presence of strong (self-generated) magnetic turbulence and put forward a novel theory in better agreement with the particle spectra inferred with multi-wavelength observations of young SN remnants, radio-supernovae, and Galactic cosmic rays.

A multi-institutional, multi-national science team will soon submit a NASA proposal to build a constellation of spacecraft to fly into the near-Earth solar wind in a swarm spanning a multitude of scales in order to obtain critically needed measurements that will reveal the underlying dynamics of magnetized turbulence. This white paper, submitted to the Plasma 2020 Decadal Survey Committee, provides a brief overview of turbulent systems that constitute an area of compelling plasma physics research, including why this mission is needed, and how this mission will achieve the goal of revealing how energy is transferred across scales and boundaries in plasmas throughout the universe.

Using a semi-analytic model of non-linear diffusive shock acceleration, we model the 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."

We present comprehensive observations and analysis of the energetic H-stripped SN 2016coi (a.k.a. ASASSN-16fp), spanning the $\gamma$-ray through optical and radio wavelengths, acquired within the first hours to $\sim$420 days post explosion. Our campaign confirms the identification of He in the SN ejecta, which we interpret to be caused by a larger mixing of Ni into the outer ejecta layers. From the modeling of the broad bolometric light curve we derive a large ejecta mass to kinetic energy ratio ($M_{\rm{ej}}\sim 4-7\,\rm{M_{\odot}}$, $E_{\rm{k}}\sim 7-8\times 10^{51}\,\rm{erg}$). The small [\ion{Ca}{ii}] \lam\lam7291,7324 to [\ion{O}{i}] \lam\lam6300,6364 ratio ($\sim$0.2) observed in our late-time optical spectra is suggestive of a large progenitor core mass at the time of collapse. We find that SN 2016coi is a luminous source of X-rays ($L_{X}>10^{39}\,\rm{erg\,s^{-1}}$ in the first $\sim100$ days post explosion) and radio emission ($L_{8.5\,GHz}\sim7\times 10^{27}\,\rm{erg\,s^{-1}Hz^{-1}}$ at peak). These values are in line with those of relativistic SNe (2009bb, 2012ap). However, for SN 2016coi we infer substantial pre-explosion progenitor mass-loss with rate $\dot M \sim (1-2)\times 10^{-4}\,\rm{M_{\odot}yr^{-1}}$ and a sub-relativistic shock velocity $v_{sh}\sim0.15c$, in stark contrast with relativistic SNe and similar to normal SNe. Finally, we find no evidence for a SN-associated shock breakout $\gamma$-ray pulse with energy $E_{\gamma}>2\times 10^{46}\,\rm{erg}$. While we cannot exclude the presence of a companion in a binary system, taken together, our findings are consistent with a massive single star progenitor that experienced large mass loss in the years leading up to core-collapse, but was unable to achieve complete stripping of its outer layers before explosion.

Analysis of 15314 electron velocity distribution functions (VDFs) within $\pm$2 hours of 52 interplanetary (IP) shocks observed by the \emph{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 velocity distribution function, while both the halo and beam/strahl components were best fit to bi-kappa velocity distribution function. This is the first statistical study to show that the core electron distribution is better fit to a self-similar velocity distribution function 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 range of values defined by the lower and upper quartiles for the kappa exponents are $\kappa{\scriptstyle_{ec}}$ $\sim$ 5.40--10.2 for the core, $\kappa{\scriptstyle_{eh}}$ $\sim$ 3.58--5.34 for the halo, and $\kappa{\scriptstyle_{eb}}$ $\sim$ 3.40--5.16 for the beam/strahl. The lower-to-upper quartile range of symmetric bi-self-similar core exponents are $s{\scriptstyle_{ec}}$ $\sim$ 2.00--2.04, and asymmetric bi-self-similar core exponents are $p{\scriptstyle_{ec}}$ $\sim$ 2.20--4.00 for the parallel exponent, and $q{\scriptstyle_{ec}}$ $\sim$ 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.

Galaxy clusters are the largest and most massive bound objects resulting from cosmic hierarchical structure formation. Baryons account for somewhat more than 10% of that mass, with roughly 90% of the baryonic matter distributed throughout the clusters as hot ($T>1$ keV), high-$\beta$, very weakly collisional plasma; the so-called "intracluster medium" (ICM). Cluster mergers, close gravitational encounters and accretion, along with violent feedback from galaxies and relativistic jets from active galactic nuclei, drive winds, gravity waves, turbulence and shocks within the ICM. Those dynamics, in turn, generate cluster-scale magnetic fields and accelerate and mediate the transport of high-energy charged particles. Kinetic-scale, collective plasma processes define the basic character and fundamental signatures of these ICM phenomena, which are observed primarily by X-ray and radio astronomers.

This paper briefly reviews a number of fundamental measurements that need to be made in order to characterize turbulence in space plasmas such as the solar wind. It has long been known that many of these quantities require simultaneous multipoint measurements to attain a proper characterization that would reveal the fundamental physics of plasma turbulence. The solar wind is an ideal plasma for such an investigation, and it now appears to be technologically feasible to carry out such an investigation, following the pioneering Cluster and MMS missions. Quantities that need to be measured using multipoint measurements include the two-point, two-time second correlation function of velocity, magnetic field and density, and higher order statistical objects such as third and fourth order structure functions. Some details of these requirements are given here, with a eye towards achieving closure on fundamental questions regarding the cascade rate, spectral anisotropy, characteristic coherent structures, intermittency, and dissipation mechanisms that describe plasma turbuelence, as well as its variability with plasma parameters in the solar wind. The motivation for this discussion is the current planning for a proposed Helioswarm mission that would be designed to make these measurements,leading to breakthrough understanding of the physics of space and astrophysical turbulence.

We use fully kinetic particle-in-cell simulations with unprecedentedly large transverse box sizes to study particle acceleration in weakly-magnetized mildly relativistic shocks traveling at a velocity $\approx 0.75c$ and a Mach number of 15. We examine both subluminal (quasi-parallel) and superluminal (quasi-perpendicular) magnetic field orientations. We find that quasi-parallel shocks are mediated by a filamentary non-resonant (Bell) instability driven by non-thermal ions, producing magnetic fluctuations on scales comparable to the ion gyro-radius. In quasi-parallel shocks, both electrons and ions are accelerated into non-thermal power-laws whose maximum energy grows linearly with time. The upstream heating of electrons is small, and the two species enter the shock front in rough thermal equilibrium. The shock's structure is complex; the current of reflected non-thermal ions evacuates cavities in the upstream which form filaments of amplified magnetic fields once advected downstream. At late times, $10\%$ of the shock's energy goes into non-thermal protons and $\gtrsim10\%$ into magnetic fields. We find that properly capturing the magnetic turbulence driven by the non-thermal ions is important for properly measuring the energy fraction of non-thermal electrons, $\epsilon_e$. We find $\epsilon_e\sim 5\times10^{-4}$ for quasi-parallel shocks with $v=0.75c$, slightly larger than what was measured in simulations of non-relativistic shocks. In quasi-perpendicular shocks, no non-thermal power-law develops in ions or electrons. The ion acceleration efficiency in quasi-parallel shocks suggests that astrophysical objects that could host mildly relativistic quasi-parallel shocks -- for example, the jets of active galactic nuclei or microquasars -- may be important sources of cosmic rays and their secondaries, such as gamma-rays and neutrinos.

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 hydro simulations across the adiabatic and radiative stages for the gas-only case, and it predicts that typical CR acceleration efficiencies ($\approx 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.

We have performed 2D hybrid simulations of non-relativistic collisionless shocks in the presence of pre-existing energetic particles ("seeds"); such a study applies, for instance, to the re-acceleration of Galactic cosmic rays (CRs) in supernova remnant (SNR) shocks and solar wind energetic particles in heliospheric shocks. Energetic particles can be effectively reflected and accelerated regardless of shock inclination via a process that we call \emph{diffusive shock re-acceleration}. We find that reaccelerated seeds can drive the streaming instability in the shock upstream and produce effective magnetic field amplification. This can eventually trigger the injection of thermal protons even at oblique shocks that ordinarily cannot inject thermal particles. We characterize the current in reflected seeds, finding that it tends to a universal value $J\simeq e n_{cr} v_{sh}$, where $e n_{cr}$ is the seed charge density and $v_{sh}$ is the shock velocity. When applying our results to SNRs, we find that the re-acceleration of Galactic CRs can excite the Bell instability to non-linear levels in less than $\sim 10$ yr, thereby providing a minimum level of magnetic field amplification for any SNR shock. Finally, we discuss the relevance of diffusive shock re-acceleration also for other environments, such as heliospheric shocks, Galactic superbubbles, and clusters of galaxies.

We suggest that ultra-high-energy (UHE) cosmic rays (CRs) may be accelerated in ultra-relativistic flows via a one-shot mechanism, the "espresso" acceleration, in which already-energetic particles are generally boosted by a factor of $\sim\Gamma^2$ in energy, where $\Gamma$ is the flow Lorentz factor. More precisely, we consider blazar-like jets with $\Gamma\gtrsim 30$ propagating into a halo of "seed" CRs produced in supernova remnants, which can accelerate UHECRs up to $10^{20}$\,eV. Such a re-acceleration process naturally accounts for the chemical composition measured by the Pierre Auger Collaboration, which resembles the one around and above the knee in the CR spectrum, and is consistent with the distribution of potential sources in the local universe, particularly intriguing is the coincidence of the powerful blazar Mrk 421 with the hotspot reported by the Telescope Array Collaboration.

We study the thermalization, injection, and acceleration of ions with different mass/charge ratios, $A/Z$, in non-relativistic collisionless shocks via hybrid (kinetic ions-fluid electrons) simulations. In general, ions thermalize to a post-shock temperature proportional to $A$. When diffusive shock acceleration is efficient, ions develop a non-thermal 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.

We propose that ultra-high-energy (UHE) cosmic rays (CRs) above $10^{18}$eV are produced in relativistic jets of powerful active galactic nuclei via an original mechanism, which we dub "espresso" acceleration: "seed" galactic CRs with energies $\lesssim 10^{17}$eV that penetrate the jet sideways receive a "one-shot'" boost of a factor of $\sim\Gamma^2$ in energy, where $\Gamma$ is the Lorentz factor of the relativistic flow. For typical jet parameters, a few percent of the CRs in the host galaxy can undergo this process, and powerful blazars with $\Gamma\gtrsim 30$ may accelerate UHECRs up to more than $10^{20}$eV. The chemical composition of espresso-accelerated UHECRs is determined by that at the Galactic CR knee and is expected to be proton-dominated at $10^{18}$eV and increasingly heavy at higher energies, in agreement with recent observations made at the Pierre Auger Observatory.

Supernova remnants (SNRs) in Local Group galaxies offer unique insights into the origin of different types of supernovae. In order to take full advantage of these insights, one must understand the intrinsic and environmental diversity of SNRs in the context of their host galaxies. We introduce a semi-analytic model that reproduces the statistical properties of a radio continuum-selected SNR population, taking into account the detection limits of radio surveys, the range of SN kinetic energies, the measured ISM and stellar mass distribution in the host galaxy from multi-wavelength images and the current understanding of electron acceleration and field amplification in SNR shocks from first-principle kinetic simulations. Applying our model to the SNR population in M33, we reproduce the SNR radio luminosity function with a median SN rate of $\sim 3.1 \times 10^{-3}$ per year and an electron acceleration efficiency, $\epsilon_{\rm{e}} \sim 4.2 \times 10^{-3}$. We predict that the radio visibility times of $\sim 70\%$ of M33 SNRs will be determined by their Sedov-Taylor lifetimes, and correlated with the measured ISM column density, $N_H$ ($t_{\rm{vis}} \propto N_H^{-a}$, with $a \sim 0.33$) while the remaining will have visibility times determined by the detection limit of the radio survey. These observational constraints on the visibility time of SNRs will allow us to use SNR catalogs as `SN surveys' to calculate SN rates and delay time distributions in the Local Group.

Foreshock disturbances -- large-scale (~1000 km to >30,000 km), transient (~5-10 per day - lasting ~10s of seconds to several minutes) structures [1,2] - generated by suprathermal (>100 eV to 100s of keV) ions [3,4] arise upstream of Earth's bow shock formed by the solar wind colliding with the Earth's magnetosphere. They have recently been found to accelerate ions to energies of several keV [5,6]. Although electrons in Saturn's high Mach number (M > 40) bow shock can be accelerated to relativistic energies (nearly 1000 keV) [7], it has hitherto been thought impossible to accelerate electrons at the much weaker (M < 20) Earth's bow shock beyond a few 10s of keV [8]. Here we report observations of electrons energized by foreshock disturbances to energies up to at least ~300 keV. Although such energetic electrons have been previously reported, their presence has been attributed to escaping magnetospheric particles [9,10] or solar events [11]. These relativistic electrons are not associated with any solar activity nor are they of magnetospheric origin. Further, current theories of ion acceleration in foreshock disturbances cannot account for electrons accelerated to the observed relativistic energies [12-17]. These electrons are clearly coming from the disturbances, leaving us with no explanation as to their origin.

The discovery of novae as sources of ~GeV gamma-rays highlights the key role of shocks and relativistic particle acceleration in these transient systems. Although there is evidence for a spectral cut-off above energies ~1-100 GeV at particular epochs in some novae, the maximum particle energy achieved in these accelerators has remained an open question. The high densities of the nova ejecta (~10 orders of magnitude larger than in supernova remnants) render the gas far upstream of the shock neutral and shielded from ionizing radiation. The amplification of the magnetic field needed for diffusive shock acceleration requires ionized gas, thus confining the acceleration process to a narrow photo-ionized layer immediately ahead of the shock. Based on the growth rate in this layer of the hybrid non-resonant cosmic ray current-driven instability (considering also ion-neutral damping), we quantify the maximum particle energy, Emax, across the range of shock velocities and upstream densities of interest. We find values of Emax ~ 10 GeV - 10 TeV, which are broadly consistent with the inferred spectral cut-offs, but which could also in principle lead to emission extending to higher energies >100 GeV accessible to atmosphere Cherenkov telescopes, such as the planned Cherenkov Telescope Array (CTA). Detecting TeV neutrinos with IceCube in hadronic scenarios appears to be more challenging, although the prospects are improved if the shock power during the earliest, densest phases of the nova outburst is higher than is implied by the observed GeV light curves, due to downscattering of the gamma-rays by electrons within the ejecta. Novae provide ideal nearby laboratories to study magnetic field amplification and the onset of cosmic ray acceleration, because other time-dependent sources (e.g. radio supernovae) typically occur too distant to detect as gamma-ray sources.

The origin of cosmic rays (CRs) has puzzled scientists since the pioneering discovery by Victor Hess in 1912. In the last decade, however, modern supercomputers have opened a new window on the processes regulating astrophysical collisionless plasmas, allowing the study of CR acceleration via first-principles kinetic simulations. At the same time, a new-generation of X-ray and $\gamma$-ray telescopes has been collecting evidence that Galactic CRs are accelerated in the blast waves of supernova remnants (SNRs). I present state-of-the-art particle-in-cells simulations of non-relativistic shocks, in which ion and electron acceleration efficiency and magnetic field amplification are studied in detail as a function of the shock parameters. I then discuss the theoretical and observational counterparts of these findings, comparing them with predictions of diffusive shock acceleration theory and with multi-wavelength observations of young SNRs. I especially outline some major open questions, such as the possible causes of the steep CR spectra inferred from $\gamma$-ray observations of SNRs and the origin of the knee in the Galactic CR spectrum. Finally, I put such a theoretical understanding in relation with CR propagation in the Galaxy in order to bridge the gap between acceleration in sources and measurements of CRs at Earth.

We formulate a magnetohydrodynamic-particle-in-cell (MHD-PIC) method for describing the interaction between collisionless cosmic ray (CR) particles and a thermal plasma. The thermal plasma is treated as a fluid, obeying equations of ideal MHD, while CRs are treated as relativistic Lagrangian particles subject to the Lorentz force. Backreaction from CRs to the gas is included in the form of momentum and energy feedback. In addition, we include the electromagnetic feedback due to CR-induced Hall effect that becomes important when the electron-ion drift velocity of the background plasma induced by CRs approaches the Alfv\'en velocity. Our method is applicable on scales much larger than the ion inertial length, bypassing the microscopic scales that must be resolved in conventional PIC methods, while retaining the full kinetic nature of the CRs. We have implemented and tested this method in the Athena MHD code, where the overall scheme is second-order accurate and fully conservative. As a first application, we describe a numerical experiment to study particle acceleration in non-relativistic shocks. Using a simplified prescription for ion injection, we reproduce the shock structure and the CR energy spectra obtained with more self-consistent hybrid-PIC simulations, but at substantially reduced computational cost. We also show that the CR-induced Hall effect reduces the growth rate of the Bell instability and affects the gas dynamics in the vicinity of the shock front. As a step forward, we are able to capture the transition of particle acceleration from non relativistic to relativistic regimes, with momentum spectrum $f(p)\sim p^{-4}$ connecting smoothly through the transition, as expected from the theory of Fermi acceleration.

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 non-resonant hybrid (Bell) instability. Protons are injected into DSA after a few gyro-cycles of shock drift acceleration (SDA), while electrons are first pre-heated 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 at Earth and the multi-wavelength emission of astrophysical objects such as supernova remnants, radio supernovae, and galaxy clusters. We find the normalization of electron power-law is $\lesssim 10^{-2}$ that of the protons for strong nonrelativistic shocks.

We use kinetic hybrid simulations (kinetic ions - fluid electrons) to characterize the fraction of ions that are accelerated to non-thermal energies at non-relativistic collisionless shocks. We investigate the properties of the shock discontinuity and show that shocks propagating almost along the background magnetic field (quasi-parallel shocks) reform quasi-periodically on ion cyclotron scales. Ions that impinge on the shock when the discontinuity is the steepest are specularly reflected. This is a necessary condition for being injected, but it is not sufficient. Also by following the trajectories of reflected ions, we calculate the minimum energy needed for injection into diffusive shock acceleration, as a function of the shock inclination. We construct a minimal model that accounts for the ion reflection from quasi-periodic shock barrier, for the fraction of injected ions, and for the ion spectrum throughout the transition from thermal to non-thermal energies. This model captures the physics relevant for ion injection at non-relativistic astrophysical shocks with arbitrary strengths and magnetic inclinations, and represents a crucial ingredient for understanding the diffusive shock acceleration of cosmic rays.

We use 2D and 3D hybrid (kinetic ions - fluid electrons) simulations to investigate particle acceleration and magnetic field amplification at non-relativistic astrophysical shocks. We show that diffusive shock acceleration operates for quasi-parallel configurations (i.e., when the background magnetic field is almost aligned with the shock normal) and, for large sonic and Alfv\'enic Mach numbers, produces universal power-law spectra proportional to p^(-4), where p is the particle momentum. The maximum energy of accelerated ions increases with time, and it is only limited by finite box size and run time. Acceleration is mainly efficient for parallel and quasi-parallel strong shocks, where 10-20% of the bulk kinetic energy can be converted to energetic particles, and becomes ineffective for quasi-perpendicular shocks. Also, the generation of magnetic turbulence correlates with efficient ion acceleration, and vanishes for quasi-perpendicular configurations. At very oblique shocks, ions can be accelerated via shock drift acceleration, but they only gain a factor of a few in momentum, and their maximum energy does not increase with time. These findings are consistent with the degree of polarization and the morphology of the radio and X-ray synchrotron emission observed, for instance, in the remnant of SN 1006. We also discuss the transition from thermal to non-thermal particles in the ion spectrum (supra-thermal region), and we identify two dynamical signatures peculiar of efficient particle acceleration, namely the formation of an upstream precursor and the alteration of standard shock jump conditions.

We use large hybrid (kinetic protons-fluid electrons) simulations to investigate the transport of energetic particles in self-consistent electromagnetic configurations of collisionless shocks. In previous papers of this series, we showed that ion acceleration may be very efficient (up to $10-20\%$ in energy), and outlined how the streaming of energetic particles amplifies the upstream magnetic field. Here, we measure particle diffusion around shocks with different strengths, finding that the mean free path for pitch-angle scattering of energetic ions is comparable with their gyroradii calculated in the self-generated turbulence. For moderately-strong shocks, magnetic field amplification proceeds in the quasi-linear regime, and particles diffuse according to the self-generated diffusion coefficient, i.e., the scattering rate depends only on the amount of energy in modes with wavelengths comparable with the particle gyroradius. For very strong shocks, instead, the magnetic field is amplified up to non-linear levels, with most of the energy in modes with wavelengths comparable to the gyroradii of highest-energy ions, and energetic particles experience Bohm-like diffusion in the amplified field. We also show how enhanced diffusion facilitates the return of energetic particles to the shock, thereby determining the maximum energy that can be achieved in a given time via diffusive shock acceleration. The parametrization of the diffusion coefficient that we derive can be used to introduce self-consistent microphysics into large-scale models of cosmic ray acceleration in astrophysical sources, such as supernova remnants and clusters of galaxies.

We use large hybrid simulations to study ion acceleration and generation of magnetic turbulence due to the streaming of particles that are self-consistently accelerated at non-relativistic shocks. When acceleration is efficient, we find that the upstream magnetic field is significantly amplified. The total amplification factor is larger than 10 for shocks with Alfv\'enic Mach number $M=100$, and scales with the square root of $M$. The spectral energy density of excited magnetic turbulence is determined by the energy distribution of accelerated particles, and for moderately-strong shocks ($M\lesssim30$) agrees well with the prediction of resonant streaming instability, in the framework of quasilinear theory of diffusive shock acceleration. For $M\gtrsim30$, instead, Bell's non-resonant hybrid (NRH) instability is predicted and found to grow faster than resonant instability. NRH modes are excited far upstream by escaping particles, and initially grow without disrupting the current, their typical wavelengths being much shorter than the current ions' gyroradii. Then, in the nonlinear stage, most unstable modes migrate to larger and larger wavelengths, eventually becoming resonant in wavelength with the driving ions, which start diffuse. Ahead of strong shocks we distinguish two regions, separated by the free-escape boundary: the far upstream, where field amplification is provided by the current of escaping ions via NRH instability, and the shock precursor, where energetic particles are effectively magnetized, and field amplification is provided by the current in diffusing ions. The presented scalings of magnetic field amplification enable the inclusion of self-consistent microphysics into phenomenological models of ion acceleration at non-relativistic shocks.

We present the results of large hybrid (kinetic ions - fluid electrons) simulations of particle acceleration at non-relativistic collisionless shocks. Ion acceleration efficiency and magnetic field amplification are investigated in detail as a function of shock inclination and strength, and compared with predictions of diffusive shock acceleration theory, for shocks with Mach number up to 100. Moreover, we discuss the relative importance of resonant and Bell's instability in the shock precursor, and show that diffusion in the self-generated turbulence can be effectively parametrized as Bohm diffusion in the amplified magnetic field.

In this paper we present the first formulation of the theory of non-linear particle acceleration in collisionless shocks in the presence of neutral hydrogen in the acceleration region. The dynamical reaction of the accelerated particles, the magnetic field amplification and the magnetic dynamical effects on the shock are also included. The main new aspect consists however in accounting for charge exchange and ionization of neutral hydrogen, which profoundly change the structure of the shock, as discussed in our previous work. This important dynamical effect of neutrals is mainly associated to the so-called neutral return flux, namely the return of hot neutrals from the downstream region to the upstream, where they deposit energy and momentum through charge exchange and ionization. We also present the self-consistent calculation of Balmer line emission from the shock region and discuss how to use measurements of the anomalous width of the different components of the Balmer line to infer the cosmic ray acceleration efficiency in supernova remnants showing Balmer emission: the broad Balmer line, which is due to charge exchange of hydrogen atoms with hot ions downstream of the shock, is shown to become narrower as a result of the energy drainage into cosmic rays, while the narrow Balmer line, due to charge exchange in the cosmic-ray--induced precursor, is shown to become broader. In addition to these two well-known components, the neutral return flux leads to the formation of a third component with intermediate width: this too contains information on ongoing processes at the shock.

Supernova remnants (SNRs) have been regarded for many decades as the sources of Galactic cosmic rays (CRs) up to a few PeV. However, only with the advent of Fermi-LAT it has been possible to detect - at least in some SNRs - \gamma-rays whose origin is unequivocally hadronic, namely due to the decay of neutral pions produced by collisions between relativistic nuclei and the background plasma. When coupled with observations in other bands (from radio to TeV \gamma-rays), Fermi-LAT data present evidence for CR spectra significantly steeper than the standard prediction of diffusive shock acceleration, forcing us to rethink our theoretical understanding of efficient particle energization at strong shocks. We outline how, by including the effects of CR-triggered magnetic field amplification, it is possible to reconcile non-linear models of diffusive shock acceleration with \gamma-ray observations, in particular providing a successful application of such a theory to Tycho's SNR. Finally, we show how kinetic simulations can investigate the microphysics of the non-linear coupling of accelerated particles and magnetic fields, probing from first principles the efficiency of the Fermi mechanism at strong shocks

We used unprecedentedly large 2D and 3D hybrid (kinetic ions - fluid electrons) simulations of non-relativistic collisionless strong shocks in order to investigate the effects of self-consistently accelerated ions on the overall shock dynamics. The current driven by suprathermal particles streaming ahead of the shock excites modes transverse to the background magnetic field. The Lorentz force induced by these self-amplified fields tends to excavate tubular, underdense, magnetic-field-depleted cavities that are advected with the fluid and perturb the shock surface, triggering downstream turbulent motions. These motions further amplify the magnetic field, up to factors of 50-100 in knot-like structures. Once downstream, the cavities tend to be filled by hot plasma plumes that compress and stretch the magnetic fields in elongated filaments; this effect is particularly evident if the shock propagates parallel to the background field. Highly-magnetized knots and filaments may provide explanations for the rapid X-ray variability observed in RX J1713.7-3946 and for the regular pattern of X-ray bright stripes detected in Tycho's supernova remnant.

We outline the main features of nuclei acceleration at supernova remnant forward shocks, stressing the crucial role played by self-amplified magnetic fields in determining the energy spectrum observed in this class of sources. In particular, we show how the standard predictions of the non-linear theory of diffusive shock acceleration has to be completed with an additional ingredient, which we propose to be the enhanced velocity of the magnetic irregularities particles scatter against, to reconcile the theory of efficient particle acceleration with recent observations of gamma-ray bright supernova remnants.

We investigate the physics of particle acceleration at non-relativistic shocks exploiting two different and complementary approaches, namely a semi-analytic modeling of cosmic-ray modified shocks and large hybrid (kinetic protons/fluid electrons) simulations. The former technique allows us to extract some information from the multi-wavelength observations of supernova remnants, especially in the gamma-ray band, while the latter returns fundamental insights into the details of particle injection and magnetic field amplification via plasma instabilities. In particular, we present the results of large hybrid simulations of non-relativistic shocks, discussing the properties of the transition from the thermal to the non-thermal component, the spectrum of which turns out to be the power-law predicted by first-order Fermi acceleration. Along with a rather effective magnetic field amplification, we find that more than 20% of the bulk energy is converted in non-thermal particles, altering significantly the dynamics of the shock and leading to the formation of a precursor.

We apply the non-linear diffusive shock acceleration theory in order to describe the properties of two supernova remnants, SN 1572 (Tycho) and SN 1604 (Kepler). By analyzing the multi-wavelength spectra, we infer that both Tycho's and Kepler's forward shocks are accelerating protons up to ~500 TeV, channeling into cosmic rays more than 10 per cent of their kinetic energy. We find that the streaming instability induced by cosmic rays is consistent with the X-ray morphology of the remnants, indicating a very efficient magnetic field amplification (up to ~300 microG). In the case of Tycho we explain the gamma-ray spectrum from the GeV up to the TeV band as due to pion decay produced in nuclear collisions by accelerated nuclei scattering against the background gas. On the other hand, due to the larger distance, the gamma-ray emission from Kepler is not detected, being below the sensitivity of the present detectors, but it should be detectable by the Cerenkov Telescope Array.

A collisionless shock may be strongly modified by the presence of neutral atoms through the processes of charge exchange between ions and neutrals and ionization of the latter. These two processes lead to exchange of energy and momentum between charged and neutral particles both upstream and downstream of the shock. In particular, neutrals that suffer a charge exchange downstream with shock-heated ions generate high velocity neutrals that have a finite probability of returning upstream. These neutrals might then deposit heat in the upstream plasma through ionization and charge exchange, thereby reducing the fluid Mach number. A consequence of this phenomenon, that we refer to as the "neutral return flux", is a reduction of the shock compression factor and the formation of a shock precursor upstream. The scale length of the precursor is determined by the ionization and charge exchange interaction lengths of fast neutrals moving towards upstream infinity. In the case of a shock propagating in the interstellar medium, the effects of ion-neutral interactions are especially important for shock velocities < 3000 km/s. Such propagation velocities are common among shocks associated with supernova remnants, the primary candidate sources for the acceleration of Galactic cosmic rays. We then investigate the effects of the return flux of neutrals on the spectrum of test-particles accelerated at the shock. We find that, for shocks slower than ~3000 km/s, the particle energy spectrum steepens appreciably with respect to the naive expectation for a strong shock, namely E^-2.

A rapidly growing amount of evidences, mostly coming from the recent gamma-ray observations of Galactic supernova remnants (SNRs), is seriously challenging our understanding of how particles are accelerated at fast shocks. The cosmic-ray (CR) spectra required to account for the observed phenomenology are in fact as steep as $E^{-2.2}--E^{-2.4}$, i.e., steeper than the test-particle prediction of first-order Fermi acceleration, and significantly steeper than what expected in a more refined non-linear theory of diffusive shock acceleration. By accounting for the dynamical back-reaction of the non-thermal particles, such a theory in fact predicts that the more efficient the particle acceleration, the flatter the CR spectrum. In this work we put forward a self-consistent scenario in which the account for the magnetic field amplification induced by CR streaming produces the conditions for reversing such a trend, allowing --- at the same time --- for rather steep spectra and CR acceleration efficiencies (about 20%) consistent with the hypothesis that SNRs are the sources of Galactic CRs. In particular, we quantitatively work out the details of instantaneous and cumulative CR spectra during the evolution of a typical SNR, also stressing the implications of the observed levels of magnetization on both the expected maximum energy and the predicted CR acceleration efficiency. The latter naturally turns out to saturate around 10-30%, almost independently of the fraction of particles injected into the acceleration process as long as this fraction is larger than about $10^{-4}$.

We aim to test the plausibility of a theoretical framework in which the $\gamma-$ray emission detected from supernova remnants is of hadronic origin, i.e., due to the decay of neutral pions produced in nuclear collisions involving relativistic nuclei. In particular, we investigate how the nature of the circumstellar medium affects the evolution of a remnant and of its $\gamma-$ray emission, stressing the role of magnetic field amplification in the prediction of expected particle spectra. A phenomenological scenario consistent with both the underlying Physics and the larger and larger amount of observational data provided by the present generation of $\gamma-$ray experiments is finally outlined and critically discussed.

Very recent gamma-ray observations of G120.1+1.4 (Tycho's) supernova remnant (SNR) by Fermi-LAT and VERITAS provided new fundamental pieces of information for understanding particle acceleration and non-thermal emission in SNRs. We want to outline a coherent description of Tycho's properties in terms of SNR evolution, shock hydrodynamics and multi-wavelength emission by accounting for particle acceleration at the forward shock via first order Fermi mechanism. We adopt here a quick and reliable semi-analytical approach to non-linear diffusive shock acceleration which includes magnetic field amplification due to resonant streaming instability and the dynamical backreaction on the shock of both cosmic rays (CRs) and self-generated magnetic turbulence. We find that Tycho's forward shock is accelerating protons up to at least 500 TeV, channelling into CRs about the 10 per cent of its kinetic energy. Moreover, the CR-induced streaming instability is consistent with all the observational evidences indicating a very efficient magnetic field amplification (up to ~300 micro Gauss). In such a strong magnetic field the velocity of the Alfv\'en waves scattering CRs in the upstream is expected to be enhanced and to make accelerated particles feel an effective compression factor lower than 4, in turn leading to an energy spectrum steeper than the standard prediction {\propto} E^-2. This latter effect is crucial to explain the GeV-to-TeV gamma-ray spectrum as due to the decay of neutral pions produced in nuclear collisions between accelerated nuclei and the background gas. The self-consistency of such an hadronic scenario, along with the fact that the concurrent leptonic mechanism cannot reproduce both the shape and the normalization of the detected the gamma-ray emission, represents the first clear and direct radiative evidence that hadron acceleration occurs efficiently in young Galactic SNRs.

We aim to test the plausibility of a theoretical framework in which the gamma-ray emission detected from supernova remnants may be of hadronic origin, i.e., due to the decay of neutral pions produced in nuclear collisions involving relativistic nuclei. In particular, we investigate the effects induced by magnetic field amplification on the expected particle spectra, outlining a phenomenological scenario consistent with both the underlying Physics and the larger and larger amount of observational data provided by the present generation of gamma experiments, which seem to indicate rather steep spectra for the accelerated particles. In addition, in order to study to study how pre-supernova winds might affect the expected emission in this class of sources, the time-dependent gamma-ray luminosity of a remnant with a massive progenitor is worked out. Solid points and limitations of the proposed scenario are finally discussed in a critical way.

In this work we investigate particle acceleration in supernova remnant shocks within a semi-analytical formalism which self-consistently accounts for particle acceleration, amplification of the magnetic field via streaming instability and back-reaction of both accelerated particles and magnetic turbulence on the shock dynamics. In particular, we study the interplay between particle injection and magnetic field amplification, showing how a phenomenological but reasonable saturation of the standard streaming instability may lead to quite steep spectra for the accelerated particles. We comment on the implications that such a scenario may have on the comprehension of the diffuse spectrum of Galactic cosmic rays and of gamma-ray observations of single remnants.

We present the non-linear theory of shock acceleration applied to SNRs expanding into partially neutral plasma. Using this theory we show how the Balmer lines detected from young SNRs can be used to test the efficiency of shocks in the production of cosmic rays. In particular we investigate the effect of charge-exchange between protons and neutral hydrogen occurring in the precursor formed ahead of the shock. In this precursor the CR pressure accelerate the ionized component of the plasma and a relative velocity between protons and neutral hydrogen is established. On the other hand the charge-exchange process tends to equilibrate ions and neutrals resulting in the heating of both components. We show that even when the shock converts only a few per cent of the total bulk kinetic energy into CRs, the heating is efficient enough to produce a detectable broadening of the narrow Balmer lines emitted by the neutral hydrogen.

We investigate the theoretical and observational implications of the acceleration of protons and heavier nuclei in supernova remnants (SNRs). By adopting a semi-analytical technique, we study the non-linear interplay among particle acceleration, magnetic field generation and shock dynamics, outlining a self-consistent scenario for the origin of the spectrum of Galactic cosmic rays as produced in this class of sources. Moreover, the inferred chemical abundances suggest nuclei heavier than Hydrogen to be relevant not only in the shock dynamics but also in the calculation of the gamma-ray emission from SNRs due to the decay of neutral pions produced in nuclear interactions.

We describe a semi-analytical approach to non-linear diffusive shock acceleration in the case in which nuclei other than protons are also accelerated. The structure of the shock is determined by the complex interplay of all nuclei, and in turn this shock structure determines the spectra of all components. The magnetic field amplification upstream is described as due to streaming instability of all nuclear species. The amplified magnetic field is then taken into account for its dynamical feedback on the shock structure as well as in terms of the induced modification of the velocity of the scattering centers that enters the particle transport equation. The spectra of accelerated particles are steep enough to be compared with observed cosmic ray spectra only if the magnetic field is sufficiently amplified and the scattering centers have high speed in the frame of the background plasma. We discuss the implications of this generalized approach on the structure of the knee in the all-particle cosmic ray spectrum, which we interpret as due to an increasingly heavier chemical composition above $10^{15}$eV. The effects of a non trivial chemical composition at the sources on the gamma ray emission from a supernova remnant when gamma rays are of hadronic origin are also discussed.

The supernova paradigm for the origin of galactic cosmic rays has been deeply affected by the development of the non-linear theory of particle acceleration at shock waves. Here we discuss the implications of applying such theory to the calculation of the spectrum of cosmic rays at Earth as accelerated in supernova remnants and propagating in the Galaxy. The spectrum is calculated taking into account the dynamical reaction of the accelerated particles on the shock, the generation of magnetic turbulence which enhances the scattering near the shock, and the dynamical reaction of the amplified field on the plasma. Most important, the spectrum of cosmic rays at Earth is calculated taking into account the flux of particles escaping from upstream during the Sedov-Taylor phase and the adiabatically decompressed particles confined in the expanding shell and escaping at later times. We show how the spectrum obtained in this way is well described by a power law in momentum with spectral index close to -4, despite the concave shape of the instantaneous spectra of accelerated particles. On the other hand we also show how the shape of the spectrum is sensible to details of the acceleration process and environment which are and will probably remain very poorly known.

We provide a both qualitative and quantitative comparison among different approaches aimed to solve the problem of non-linear diffusive acceleration of particles at shocks. In particular, we show that state-of-the-art models (numerical, Monte Carlo and semi-analytical), even if based on different physical assumptions and implementations, for typical environmental parameters lead to very consistent results in terms of shock hydrodynamics, cosmic ray spectrum and also escaping flux spectrum and anisotropy. Strong points and limits of each approach are also discussed, as a function of the problem one wants to study.

We present here a semi-analytical solution of the problem of particle acceleration at non-linear shock waves with a free escape boundary at some location upstream. This solution, besides allowing us to determine the spectrum of particles accelerated at the shock front, including the shape of the cutoff at some maximum momentum, also allows us to determine the spectrum of particles escaping the system from upstream. This latter aspect of the problem is crucial for establishing a connection between the accelerated particles in astrophysical sources, such as supernova remnants, and the cosmic rays observed at the Earth. An excellent approximate solution, which leads to a computationally fast calculation of the structure of shocks with an arbitrary level of cosmic ray modification, is also obtained.

The theory of Non-Linear Diffusive Shock Acceleration (NLDSA) predicts the formation of a precursor upstream of the shock, where accelerated particles diffuse and induce magnetic field amplification through streaming instability. The non detection of this precursor in X-rays in {\it Chandra} observations of the north-eastern region of SN 1006 (G329.6+14.6) led to impose an upper limit to the X-ray emission generated by accelerated electrons diffusing in this precursor, at an emissivity level of <1.5 per cent of the emission from the downstream region (Long et al. 2003). This has been used as an argument against Fermi acceleration at this shock. Here we calculate the spectrum and spatial distribution of accelerated particles in SN 1006 and show that Chandra results (including more recent data) are in perfect agreement with the predictions of NLDSA suggesting efficient particle acceleration and magnetic field amplification upstream of the shock by a factor ~10.

Stationary solutions to the problem of particle acceleration at shock waves in the non-linear regime, when the dynamical reaction of the accelerated particles on the shock cannot be neglected, are known to show a prominent energy flux escaping from the shock towards upstream infinity. On physical grounds, the escape of particles from the upstream region of a shock has to be expected in all those situations in which the maximum momentum of accelerated particles, $p_{max}$, decreases with time, as is the case for the Sedov-Taylor phase of expansion of a shell Supernova Remnant, when both the shock velocity and the cosmic ray induced magnetization decrease. In this situation, at each time $t$, particles with momenta larger than $p_{max}(t)$ leave the system from upstream, carrying away a large fraction of the energy if the shock is strongly modified by the presence of cosmic rays. This phenomenon is of crucial importance for explaining the cosmic ray spectrum detected at Earth. In this paper we discuss how this escape flux appears in the different approaches to non-linear diffusive shock acceleration, and especially in the quasi-stationary semi-analytical kinetic ones. We apply our calculations to the Sedov-Taylor phase of a typical supernova remnant, including in a self-consistent way particle acceleration, magnetic field amplification and the dynamical reaction on the shock structure of both particles and fields. Within this framework we calculate the temporal evolution of the maximum energy reached by the accelerated particles and of the escape flux towards upstream infinity. The latter quantity is directly related to the cosmic ray spectrum detected at Earth.

We present a semi-analytical kinetic calculation of the process of non-linear diffusive shock acceleration (NLDSA) which includes the magnetic field amplification due to cosmic ray induced streaming instability, the dynamical reaction of the amplified magnetic field and the possible effects of turbulent heating. The approach is specialized to parallel shock waves and the parameters we chose are the ones appropriate to forward shocks in Supernova Remnants. Our calculation allows us to show that the net effect of the amplified magnetic field is to enhance the maximum momentum of accelerated particles while reducing the concavity of the spectra, with respect to the standard predictions of NLDSA. This is mainly due to the dynamical reaction of the amplified field on the shock, which noticeably reduces the modification of the shock precursor. The total compression factors which are obtained for parameters typical of supernova remnants are $R_{tot}\sim 7-10$, in good agreement with the values inferred from observations. The strength of the magnetic field produced through excitation of streaming instability is found in good agreement with the values inferred for several remnants if the thickness of the X-ray rims are interpreted as due to severe synchrotron losses of high energy electrons. We also discuss the relative role of turbulent heating and magnetic dynamical reaction in driving the reduction of the precursor modification.

Recent observations of greatly amplified magnetic fields ($\delta B/B\sim 100$) around supernova shocks are consistent with the predictions of the non-linear theory of particle acceleration (NLT), if the field is generated upstream of the shock by cosmic ray induced streaming instability. The high acceleration efficiencies and large shock modifications predicted by NLT need however to be mitigated to confront observations, and this is usually assumed to be accomplished by some form of turbulent heating. We show here that magnetic fields with the strength inferred from observations have an important dynamical role on the shock, and imply a shock modification substantially reduced with respect to the naive unmagnetized case. The effect appears as soon as the pressure in the turbulent magnetic field becomes comparable with the pressure of the thermal gas. The relative importance of this unavoidable effect and of the poorly known turbulent heating is assessed. More specifically we conclude that even in the cases in which turbulent heating may be of some importance, the dynamical reaction of the field cannot be neglected, as instead is usually done in most current calculations.

Particle acceleration at non-relativistic shocks can be very efficient, leading to the appearance of non-linear effects due to the dynamical reaction of the accelerated particles on the shock structure and to the non-linear amplification of the magnetic field in the shock vicinity. The value of the maximum momentum $p_{max}$ in these circumstances cannot be estimated using the classical results obtained within the framework of test particle approaches. We provide here the first attempt at estimating $p_{max}$ in the cosmic ray modified regime, taking into account the non-linear effects mentioned above.