Source author record

Damiano Caprioli

Damiano Caprioli appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

27works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

27 published item(s)

preprint2022arXiv

A Kinetic Study of the Saturation of the Bell Instability

The nonresonant cosmic ray instability, predicted by Bell (2004), is thought to play an important role in the acceleration and confinement of cosmic rays (CR) 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 hybrid simulations (with kinetic ions and fluid electrons), we study the non-linear evolution of the Bell instability as a function of the parameters of the CR population. We find that saturation is achieved when the magnetic pressure in the amplified field is comparable to the initial CR momentum flux.

preprint2022arXiv

Cosmic-ray generated bubbles around their sources

Cosmic rays are thought to escape their sources streaming along the local magnetic field lines. We show that this phenomenon generally leads to the excitation of both resonant and non-resonant streaming instabilities. The self-generated magnetic fluctuations induce particle diffusion in extended regions around the source, so that cosmic rays build up a large pressure gradient. By means of two-dimensional (2D) and three-dimensional (3D) hybrid particle-in-cell simulations, we show that such a pressure gradient excavates a cavity around the source and leads to the formation of a cosmic-ray dominated bubble, inside which diffusivity is strongly suppressed. Based on the trends extracted from self-consistent simulations, we estimate that, in the absence of severe damping of the self-generated magnetic fields, the bubble should keep expanding until pressure balance with the surrounding medium is reached, corresponding to a radius of $\sim 10-50$ pc. The implications of the formation of these regions of low diffusivity for sources of Galactic cosmic rays are discussed. Special care is devoted to estimating the self-generated diffusion coefficient and the grammage that cosmic rays might accumulate in the bubbles before moving into the interstellar medium. Based on the results of 3D simulations, general considerations on the morphology of the $γ$-ray and synchrotron emission from these extended regions also are outlined.

preprint2022arXiv

Relativistic Asymmetric 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 non-thermal 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.

preprint2022arXiv

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.

preprint2021arXiv

Lepton-driven Non-resonant Streaming Instability

A strong super-Alfvénic drift of energetic particles (or cosmic rays, CRs) in a magnetized plasma can amplify the magnetic field significantly through non-resonant streaming instability (NRSI). While the traditional analysis is done for an ion current, here we use kinetic particle-in-cell simulations to study how the NRSI behaves when it is driven by electrons or by a mixture of electrons and positrons. In particular, we characterize growth rate, spectrum, and helicity of the unstable modes, as well the level of magnetic field at saturation. Our results are potentially relevant for several space/astrophysical environments (e.g, electron strahl in the solar wind, at oblique non-relativistic shocks, around pulsar wind nebulae) and also in laboratory experiments.

preprint2020arXiv

Electron acceleration in non-relativistic quasi-perpendicular collisionless shocks

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é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.

preprint2020arXiv

Electron energy partition across interplanetary shocks: III. Analysis

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.

preprint2020arXiv

Nonsecondary origin of cosmic ray positrons

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.

preprint2019arXiv

dHybridR: a Hybrid--Particle-in-Cell Code Including Relativistic Ion Dynamics

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.

preprint2019arXiv

Electron energy partition across interplanetary shocks: II. Statistics

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$, $β{\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.

preprint2016arXiv

Supernova Remnants in the Local Group I: A model for the radio luminosity function and visibility times of supernova remnants

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, $ε_{\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.

preprint2015arXiv

Cosmic-ray Acceleration and Propagation

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 $γ$-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 $γ$-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.

preprint2015arXiv

Magnetohydrodynamic-Particle-in-Cell Method for Coupling Cosmic Rays with a Thermal Plasma: Application to Non-relativistic Shocks

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én 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.

preprint2015arXiv

Novae as Tevatrons: Prospects for CTA and IceCube

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.

preprint2015arXiv

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 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.

preprint2014arXiv

Hybrid Simulations of Particle Acceleration at 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.

preprint2014arXiv

Simulations and Theory of Ion Injection at Non-relativistic Collisionless 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.

preprint2014arXiv

Simulations of Ion Acceleration at Non-relativistic Shocks. I. Acceleration Efficiency

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énic 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.

preprint2014arXiv

Simulations of Ion Acceleration at Non-relativistic Shocks. II. Magnetic Field Amplification

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énic 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.

preprint2014arXiv

Simulations of Ion Acceleration at Non-relativistic Shocks. III. Particle Diffusion

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.

preprint2013arXiv

From E. Fermi to Fermi-LAT: watching particle acceleration in supernova remnants

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 - γ-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 γ-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 γ-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

preprint2013arXiv

New insights on hadron acceleration at supernova remnant shocks

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.

preprint2012arXiv

Cosmic-ray acceleration in supernova remnants: non-linear theory revised

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}$.

preprint2012arXiv

Fermi acceleration at supernova remnant shocks

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.

preprint2011arXiv

Supernova remnants as cosmic ray factories

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.

preprint2011arXiv

Understanding hadronic $γ$-ray emission from supernova remnants

We aim to test the plausibility of a theoretical framework in which the $γ-$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 $γ-$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 $γ-$ray experiments is finally outlined and critically discussed.

preprint2011arXiv

Understanding hadronic gamma-ray emission from supernova remnants

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.