The Origin of Cosmic Rays from Supernova Remnants

The origin of cosmic rays is one of the key questions in high-energy astrophysics. Supernovae have been always considered as the dominant sources of cosmic rays below the energy spectrum knee. Multi-wavelength observations indeed show that supernova remnants are capable for accelerating particles into sub-PeV (10${}^{15}$ eV) energies. Diffusive shock acceleration is considered as one of the most efficient acceleration mechanisms of astrophysical high-energy particles, which may just operate effectively in the large-scale shocks of supernova remnants. Recently, a series of high-precision ground and space experiments have greatly promoted the study of cosmic rays and supernova remnants. New observational features challenge the classical acceleration model by diffusive shock and the application to the scenario of supernova remnants for the origin of Galactic cosmic rays, and have deepened our understanding to the cosmic high-energy phenomena. In combination with the time evolution of radiation energy spectrum of supernova remnants, a time-dependent particle acceleration model is established, which can not only explain the anomalies in cosmic-ray distributions around 200 GV, but also naturally form the cosmic-ray spectrum knee, even extend the contribution of supernova particle acceleration to cosmic ray flux up to the spectrum ankle. This model predicts that the high-energy particle transport behavior is dominated by the turbulent convection, which needs to be verified by future observations and plasma numerical simulations relevant to the particle transport.     

����λ��α�෨�͵����в̽����޸�����

?????λ??α????????в????????????????????????к????????????????????λ??α????????????????????в????ж???????????????λ??α????????????????????????????????????????????????????????к??????????     

Three-dimensional hydrodynamical simulations of the large-scale structure of W50−SS433

We present 3D hydrodynamical simulations of a precessing jet propagating inside a supernova remnant (SNR) shell, particularly applied to the W50−SS433 system in a search for the origin of its peculiar elongated morphology. Several runs were carried out with different values for the mass-loss rate of the jet, the initial radius of the SNR, and the opening angle of the precession cone. We found that our models successfully reproduce the scale and morphology of W50 when the opening angle of the jets is set to 10° or if this angle linearly varies with time. For these models, more realistic runs were made considering that the remnant is expanding into an interstellar medium with an exponential density profile (as H  i observations suggest). Taking into account all these ingredients, the large-scale morphology of the W50−SS433 system, including the asymmetry between the lobes (formed by the jet–SNR interaction), is well reproduced.     

How rapidly do neutron stars spin at birth? Constraints from archival X-ray observations of extragalactic supernovae

Traditionally, studies aimed at inferring the distribution of birth periods of neutron stars are based on radio surveys. Here we propose an independent method to constrain the pulsar spin periods at birth based on their X-ray luminosities. In particular, the observed luminosity distribution of supernovae (SNe) poses a constraint on the initial rotational energy of the embedded pulsars, via the     correlation found for radio pulsars, and under the assumption that this relation continues to hold beyond the observed range. We have extracted X-ray luminosities (or limits) for a large sample of historical SNe observed with Chandra , XMM and Swift , which have been firmly classified as core-collapse SNe. We have then compared these observational limits with the results of Monte Carlo simulations of the pulsar X-ray luminosity distribution for a range of values of the birth parameters. We find that a pulsar population dominated by millisecond periods at birth is ruled out by the data.     

Non-thermal emission from old supernova remnants

We study the non-thermal emission from old shell-type supernova remnants (SNRs) on the frame of a time-dependent model. In this model, the time-dependent non-thermal spectra of both primary electrons and protons as well as secondary electron/positron (e^±) pairs can be calculated numerically by taking into account the evolution of the secondary e^± pairs produced from proton–proton (p–p) interactions as accelerated protons collide with the ambient matter in an SNR. The multiwavelength photon spectrum for a given SNR can be produced through leptonic processes such as electron/positron synchrotron radiation, bremsstrahlung and inverse Compton scattering as well as hadronic interaction. Our results indicate that the non-thermal emission of the secondary e^± pairs is becoming more and more prominent when the SNR ages in the radiative phase because the source of the primary electrons has been cut off and the electron synchrotron energy loss is significant for a radiative SNR, whereas the secondary e^± pairs can be produced continuously for a long time in the phase due to the large energy-loss time for the p–p interaction. We apply the model to two old SNRs, G8.7−0.1 and G23.3−0.3, and the predicted results can explain the observed multiwavelength photon spectra for the two sources.     

Star formation triggered by supernova explosions in young galaxies

We study the evolution of supernova remnants in a low-metallicity medium   Z /Z_⊙= 10⁻⁴ to 10⁻²  in the early universe, using one-dimensional hydrodynamics with non-equilibrium chemistry. Once a post-shock layer is able to cool radiatively, a dense shell forms behind the shock. If this shell becomes gravitationally unstable and fragments into pieces, next-generation stars are expected to form from these fragments. To explore the possibility of this triggered star formation, we apply a linear perturbation analysis of an expanding shell to our results and constrain the parameter range of ambient density, explosion energy and metallicity where fragmentation of the shell occurs. For the explosion energy of  10⁵¹ erg (10⁵² erg)  , the shell fragmentation occurs for ambient densities higher than  ≳10² cm⁻³ (10 cm⁻³  ), respectively. This condition depends little on the metallicity in the ranges we examined. We find that the mode of star formation triggered occurs only in massive  (≳10⁸ M_⊙)  haloes.     

The multi-band non-thermal emission from the supernova remnant RX J1713.7−3946

Non-thermal X-rays and very high energy (VHE) γ-rays have been detected from the supernova remnant (SNR) RX J1713.7−3946, and the recent observations with the Suzaku satellite clearly reveal a spectral cut-off in the X-ray spectrum, which directly relates to the cut-off of the energy spectrum of the parent electrons. However, whether the origin of the VHE γ-rays from the SNR is hadronic or leptonic is still in debate. We studied the multi-band non-thermal emission from RX J1713.7−3946 based on a semi-analytical approach towards the non-linear shock acceleration process by including the contribution of the accelerated electrons to the non-thermal radiation. The results show that the multi-band observations on RX J1713.7−3946 can be well explained in the model with appropriate parameters, and the TeV γ-rays have hadronic origin, i.e. they are produced via proton–proton interactions as the relativistic protons accelerated by the shock collide with the ambient matter.     

Saturated magnetic field amplification at supernova shocks

Cosmic ray streaming instabilities at supernova shocks are discussed in the quasi-linear diffusion formalism which takes into account the feedback effect of wave growth on the cosmic ray streaming motion. In particular, the non-resonant instability that leads to magnetic field amplification in the short wavelength regime is considered. The linear growth rate is calculated using kinetic theory for a streaming distribution. We show that the non-resonant instability is actually driven by a compensating current in the background plasma. The non-resonant instability can develop into a non-linear regime generating turbulence. The saturation of the amplified magnetic fields due to particle diffusion in the turbulence is derived analytically. It is shown that the evolution of parallel and perpendicular cosmic ray pressures is predominantly determined by non-resonant diffusion. However, the saturation is determined by resonant diffusion which tends to reduce the streaming motion through pitch angle scattering. The saturated level can exceed the mean background magnetic field.     

Dynamical feedback of self-generated magnetic fields in cosmic ray modified shocks

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 notably reduces the modification of the shock precursor. The total compression factors which are obtained for parameters typical of supernova remnants are   R _tot∼ 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 is 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.