Giovanni Morlino

Gamma ray emission from SNR RX J1713.7-3946 and the origin of galactic cosmic rays

We calculate the flux of non-thermal radiations from the supernova remnant (SNR) RX J1713.7-3946 in the context of the non-linear theory of particle acceleration at shocks, which allows us to take into account self-consistently the dynamical reaction of the accelerated particles, the generation of magnetic fields in the shock proximity and the dynamical reaction of the magnetic field on the plasma. When the fraction of particles which get accelerated is of the order of ∼10 -4 , we find that the strength of the magnetic field obtained as a result of streaming instability induced by cosmic rays is compatible with the interpretation of the X-ray emitting filaments being produced by strong synchrotron losses in ∼100 μG magnetic fields. The maximum energy of accelerated protons is ≥10 5 GeV. If the X-ray filaments are explained in alternative ways, the constraint on the magnetic field downstream of the shock disappears and the HESS data can be marginally fitted with ICS of relativistic electrons off a complex population of photons, tailored to comprise cosmic microwave background and ambient infrared/optical photons. The fit, typically poor at the highest energies, requires a large density of target photons within the remnant; only a fraction of the order of ∼10 -6 of the background particles gets accelerated; the local magnetic field is of the order of ∼20 μG and the maximum energy of protons is much lower than the knee energy. Current HESS gamma-ray observations combined with recent X-ray observations by Suzaku do not allow as yet to draw a definitive conclusion on whether RX J1713.7-3946 is an efficient cosmic ray accelerator, although at the present time a hadronic interpretation of HESS data seems more likely. We discuss the implications of our results for the GLAST gamma-ray telescope, which should be able to discriminate the two scenarios discussed above, by observing the shape of the gamma-ray spectrum at lower energies.

The microphysics of collisionless shock waves

Collisionless shocks, that is shocks mediated by electromagnetic processes, are customary in space physics and in astrophysics. They are to be found in a great variety of objects and environments: magnetospheric and heliospheric shocks, supernova remnants, pulsar winds and their nebulæ, active galactic nuclei, gamma-ray bursts and clusters of galaxies shock waves. Collisionless shock microphysics enters at different stages of shock formation, shock dynamics and particle energization and/or acceleration. It turns out that the shock phenomenon is a multi-scale non-linear problem in time and space. It is complexified by the impact due to high-energy cosmic rays in astrophysical environments. This review adresses the physics of shock formation, shock dynamics and particle acceleration based on a close examination of available multi-wavelength or in situ observations, analytical and numerical developments. A particular emphasis is made on the different instabilities triggered during the shock formation and in association with particle acceleration processes with regards to the properties of the background upstream medium. It appears that among the most important parameters the background magnetic field through the magnetization and its obliquity is the dominant one. The shock velocity that can reach relativistic speeds has also a strong impact over the development of the micro-instabilities and the fate of particle acceleration. Recent developments of laboratory shock experiments has started to bring some new insights in the physics of space plasma and astrophysical shock waves. A special section is dedicated to new laser plasma experiments probing shock physics.

Cosmic ray driven Galactic winds

The escape of cosmic rays from the Galaxy leads to a gradient in the cosmic ray pressure that acts as a force on the background plasma, in the direction opposite to the gravitational pull. If this force is large enough to win against gravity, a wind can be launched that removes gas from the Galaxy, thereby regulating several physical processes, including star formation. The dynamics of these cosmic ray driven winds is intrinsically non-linear in that the spectrum of cosmic rays determines the characteristics of the wind (velocity, pressure, magnetic field) and in turn the wind dynamics affects the cosmic ray spectrum. Moreover, the gradient of the cosmic ray distribution function causes excitation of Alfven waves, that in turn determine the scattering properties of cosmic rays, namely their diffusive transport. These effects all feed into each other so that what we see at the Earth is the result of these non-linear effects. Here we investigate the launch and evolution of such winds, and we determine the implications for the spectrum of cosmic rays by solving together the hydrodynamical equations for the wind and the transport equation for cosmic rays under the action of self-generated diffusion and advection with the wind and the self-excited Alfven waves.

On the radial distribution of Galactic cosmic rays

The spectrum and morphology of the diffuse Galactic gamma-ray emission carries valuable information on cosmic ray (CR) propagation. Recent results obtained by analyzing Fermi-LAT data accumulated over seven years of observation show a substantial variation of the CR spectrum as a function of the distance from the Galactic Center. The spatial distribution of the CR density in the outer Galaxy appears to be weakly dependent upon the galactocentric distance, as found in previous studies as well, while the density in the central region of the Galaxy was found to exceed the value measured in the outer Galaxy. At the same time, Fermi-LAT data suggest a gradual spectral softening while moving outward from the center of the Galaxy to its outskirts. These findings represent a challenge for standard calculations of CR propagation based on assuming a uniform diffusion coefficient within the Galactic volume. Here we present a model of non-linear CR propagation in which transport is due to particle scattering and advection off self-generated turbulence. We find that for a realistic distribution of CR sources following the spatial distribution of supernova remnants and the space dependence of the magnetic field on galactocentric distance, both the spatial profile of CR density and the spectral softening can easily be accounted for.

Broad Balmer line emission and cosmic ray acceleration efficiency in supernova remnant shocks

Balmer emission may be a powerful diagnostic tool to test the paradigm of cosmic ray (CR) acceleration in young supernova remnant (SNR) shocks. The width of the broad Balmer line is a direct indicator of the downstream plasma temperature. In case of efficient particle acceleration an appreciable fraction of the total kinetic energy of the plasma is channeled into CRs, therefore the downstream temperature decreases and so does the broad Balmer line width. This width also depends on the level of thermal equilibration between ions and neutral hydrogen atoms in the downstream. Since in general in young SNR shocks only a few charge exchange (CE) reactions occur before ionization, equilibration between ions and neutrals is not reached, and a kinetic description of the neutrals is required in order to properly compute Balmer emission. We provide a method for the calculation of Balmer emission using a self-consistent description of the shock structure in the presence of neutrals and CRs. We use a recently developed semi-analytical approach, where neutral particles, ionized plasma, accelerated particles and magnetic fields are all coupled together through the mass, momentum and energy flux conservation equations. The distribution of neutrals is obtained from the full Boltzmann equation in velocity space, coupled to Maxwellian ions through ionization and CE processes. The computation is also improved with respect to previous work thanks to a better approximation for the atomic interaction rates. We find that for shock speeds >2500km/s the distribution of broad neutrals never approaches a Maxwellian and its moments differ from those of the ionized component. These differences reflect into a smaller FWHM than predicted in previous calculations, where thermalization was assumed. The method presented here provides a realistic estimate of particle acceleration efficiency in Balmer dominated shocks.

Non-linear diffusion of cosmic rays escaping from supernova remnants – I. The effect of neutrals

Supernova remnants are believed to be the main sources of galactic cosmic rays (CR). Within this framework, particles are accelerated at supernova remnant shocks and then released in the interstellar medium. The mechanism through which CRs are released and the way in which they propagate still remain open issues. The main difficulty is the high non-linearity of the problem: CRs themselves excite the magnetic turbulence that confines them close to their sources. We solve numerically the coupled differential equations describing the evolution in space and time of the escaping particles and of the waves generated through the CR streaming instability. The warm ionized and warm neutral phases of the interstellar medium are considered. These phases occupy the largest fraction of the disc volume, where most supernovae explode, and are characterized by the significant presence of neutral particles. The friction between those neutrals and ions results in a very effective wave damping mechanism. It is found that streaming instability affects the propagation of CRs even in the presence of ion-neutral friction. The diffusion coefficient can be suppressed by more than a factor of ∼2 over a region of few tens of pc around the remnant. The suppression increases for smaller distances. The propagation of ≈10 GeV particles is affected for several tens of kiloyears after escape, while ≈1 TeV particles are affected for few kiloyears. This might have a great impact on the interpretation of gamma-ray observations of molecular clouds located in the vicinity of supernova remnants.

Cosmic ray acceleration and Balmer emission from SNR 0509-67.5

Context. Observation of Balmer lines from the region around the forward shock of supernova remnants may provide precious information on the shock dynamics and on the efficiency of particle acceleration at the shock. Aims. We calculate the Balmer line emission and the shape of the broad Balmer line for parameter values suitable for SNR 0509-67.5, as a function of the cosmic ray acceleration efficiency and of the level of thermal equilibrium between electrons and protons behind the shock. This calculation aims to use the width of the broad Balmer line emission to infer the cosmic ray acceleration efficiency in this remnant. Methods. We use the recently developed nonlinear theory of diffusive shock acceleration in the presence of neutrals. The semi-analytical approach that we developed includes a description of magnetic field amplification as due to resonant streaming instability, the dynamical reaction of both accelerated particles and turbulent magnetic field on the shock, and all channels of interaction between neutral atoms and background plasma that change the shock dynamics. Results. We achieve a quantitative assessment of the CR acceleration efficiency in SNR 0509-67.5 as a function of the shock velocity and different levels of electron-proton thermalization in the shock region. If the shock moves faster than ~4500 km s -1 , one can conclude that particle acceleration must be taking place with an efficiency of several tens of percent. For lower shock velocity the evidence of particle acceleration becomes less clear because of the uncertainty in the electron-ion equilibrium downstream. We also discuss the role of future measurements of the narrow Balmer line.

Spectra of accelerated particles at supernova shocks in the presence of neutral hydrogen: the case of Tycho

The presence of neutral hydrogen in the shock proximity changes the structure of the shock and affects the spectra of particles accelerated through the first-order Fermi mechanism. This phenomenon has profound implications for the interpretation of the multifrequency spectra of radiation from supernova remnants. Neutrals that undergo charge exchange with hot ions downstream of the shock may result in fast neutrals moving towards the upstream gas, where they can suffer additional charge exchange or ionisation reactions, thereby depositing energy and momentum upstream. Here we discuss the implications of this neutral return flux, which was already predicted in our previous work on neutral mediated supernova shocks, and show how the spectra of accelerated particles turn out to be appreciably steeper than

Gamma rays and neutrinos from SNR RX J1713.7–3946

p^{-4}

The role of ionization in the shock acceleration theory

, thereby affecting the gamma ray spectra from supernova remnants in general and from Tycho specifically. The theory that describes non-linear diffusive shock acceleration in the presence of neutral hydrogen has been developed in recent years. Here we use a semi-analytical theory developed in previous work and specialise our predictions to the case of the Tycho supernova shock, where there is evidence that the spectrum of the accelerated cosmic rays is steeper than expected from the traditional theory of diffusive shock acceleration. We show that, if the fraction of neutral hydrogen in the vicinity of the Tycho supernova shock is, as suggested by observations, ~70-90, then spectra of accelerated protons steeper than