The Cherenkov Telescope Array potential for the study of young supernova remnants

Supernova remnants (SNRs) are among the most important targets for γ-ray observatories. Being prominent non-thermal sources, they are very likely responsible for the acceleration of the bulk of Galactic cosmic rays (CRs). To firmly establish the SNR paradigm for the origin of cosmic rays, it should be confirmed that protons are indeed accelerated in, and released from, SNRs with the appropriate flux and spectrum. This can be done by detailed theoretical models which account for microphysics of acceleration and various radiation processes of hadrons and leptons. The current generation of Cherenkov telescopes has insufficient sensitivity to constrain theoretical models. A new facility, the Cherenkov Telescope Array (CTA), will have superior capabilities and may finally resolve this long standing issue of high-energy astrophysics. We want to assess the capabilities of CTA to reveal the physics of various types of SNRs in the initial 2000 years of their evolution. During this time, the efficiency to accelerate cosmic rays is highest. We perform time-dependent simulations of the hydrodynamics, the magnetic fields, the cosmic-ray acceleration, and the non-thermal emission for type Ia, Ic and IIP SNRs. We calculate the CTA response to the γ-ray emission from these SNRs for various ages and distances, and we perform a realistic analysis of the simulated data. We derive distance limits for the detectability and resolvability of these SNR types at several ages. We test the ability of CTA to reconstruct their morphological and spectral parameters as a function of their distance. Finally, we estimate how well CTA data will constrain the theoretical models.     

Particle Acceleration at Shocks: Insights from Supernova Remnant Shocks

I review some basic properties of diffusive shock acceleration (DSA) in the context of young supernova remnants (SNRs). I also point out some key differences with cosmological, cluster-related shocks. DSA seems to be very efficient in strong, young SNR shocks. Provided the magnetic fields exceed some hundreds of μGauss (possibly amplified by CR related dynamics), these shocks can accelerate cosmic ray hadrons to PeV energies in the time available to them. Electron energies, limited by radiative losses, are likely limited to the TeV range. Injection of fresh particles at these shocks is poorly understood, but hadrons are much more easily injected than the more highly magnetized electrons. That seems supported by observational data, as well. So, while CR protons in young SNRs may play very major roles in the SNR evolution, the CR electron populations have minimal such impact, despite their observational importance.     

The origin of galactic cosmic rays

One century ago Viktor Hess carried out several balloon flights that led him to conclude that the penetrating radiation responsible for the discharge of electroscopes was of extraterrestrial origin. One century from the discovery of this phenomenon seems to be a good time to stop and think about what we have understood about Cosmic Rays. The aim of this review is to illustrate the ideas that have been and are being explored in order to account for the observable quantities related to cosmic rays and to summarize the numerous new pieces of observation that are becoming available. In fact, despite the possible impression that development in this field is somewhat slow, the rate of new discoveries in the last decade or so has been impressive, and mainly driven by beautiful pieces of observation. At the same time scientists in this field have been able to propose new, fascinating ways to investigate particle acceleration inside the sources, making use of multifrequency observations that range from the radio, to the optical, to X-rays and gamma rays. These ideas can now be confronted with data.I will mostly focus on supernova remnants as the most plausible sources of Galactic cosmic rays, and I will review the main aspects of the modern theory of diffusive particle acceleration at supernova remnant shocks, with special attention for the dynamical reaction of accelerated particles on the shock and the phenomenon of magnetic field amplification at the shock. Cosmic-ray escape from the sources is discussed as a necessary step to determine the spectrum of cosmic rays at the Earth. The discussion of these theoretical ideas will always proceed parallel to an account of the data being collected especially in X-ray and gamma-ray astronomy.In the end of this review I will also discuss the phenomenon of cosmic-ray acceleration at shocks propagating in partially ionized media and the implications of this phenomenon in terms of width of the Balmer line emission. This field of research has recently experienced a remarkable growth, in that Hα lines have been found to bear information on the cosmic-ray acceleration efficiency of supernova shocks.     

Particle acceleration test with Cas A multiwavelength emission

The investigation of supernova remnants (SNRs) across the electromagnetic spectrum from radio up to very high energy gamma-rays can serve as a test of the particle acceleration and touches on one of the unresolved problems of modern astrophysics, namely the origin of cosmic rays and the Galaxy's contribution to the overall cosmic ray spectrum. The multiwavelength observations of Cas A SNR demonstrated that structure and spectral features have clear signs of young SNRs and its overall properties make this object the best target to test a hypothesis of cosmic ray origin in SNRs. Studies of Cas A at very high energies by SHALON telescope showed the location of TeV gamma-ray emission region relative to the position of reveres shock. Also, the spectral energy distribution was obtained at high and very high energies. To describe the spectral and structural features of this SNR viewed in non-thermal emission, two approaches involving reverse and also both reverse and forward shocks to the mechanism of diffusive shock acceleration of cosmic rays in Cas A were applied. It is demonstrated that the observational properties of Cas A are well reproduced by the hadronic model with significant contribution of both the forward and reverse shocks in the generation of broadband emission. Calculation results suggest that the very high efficiency of particle acceleration in Cas A, which value is up to 25% of the supernova explosion energy with energy of accelerated particles not exceeding of $$ eV. Whereas, the forward shock model predicts the spectral characteristics of the TeV-gamma-emission corresponding to ones detected at 800 GeV–40 TeV that are the evidence of acceleration of the hadronic cosmic rays in shells of SNRs up to $$ eV     

Evidence for shock acceleration and intergalactic magnetic fields in a large-scale filament of galaxies ZwCl 2341.1+0000

We report the discovery of large-scale diffuse radio emission from what appears to be a large-scale filamentary network of galaxies in the region of cluster ZwCl 2341.1+0000, and stretching over an area of at least 6 h₅₀⁻¹ Mpc in diameter. Multicolour CCD observations yield photometric redshifts indicating that a significant fraction of the optical galaxies in this region is at a redshift of z=0.3. This is supported by spectroscopic measurements of 4 galaxies in the Sloan Digitized Sky Survey (SDSS) survey at a mean z=0.27. We present VLA images at λ=20 cm (NVSS) and 90 cm, showing the detailed radio structure of the filaments. Comparison with the high resolution FIRST radio survey shows that the diffuse emission is not due to known individual point sources. The diffuse radio-emission has a spectral index α≲−0.5, and is most likely synchrotron emission from relativistic charged particles in an inter-galactic magnetic field. Furthermore, this optical/radio structure is detected in X-rays by the ROSAT all-sky survey. It has a 0.1–2.4 keV luminosity of about 10⁴⁴ erg s⁻¹ and shows an extended highly non-relaxed morphology. These observations suggest that ZwCl 2341.1+0000 is possibly a proto-cluster of galaxies in which we are witnessing the process of structure formation. We show that the energetics of accretion shocks generated in forming large-scale structures are sufficient to produce enough high energy cosmic-ray (CR) electrons required to explain the observed radio emission, provided a magnetic field of strength B≳0.3 μG is present there. The latter is only a lower limit and the actual magnetic field is likely to be higher depending on the morphology of the emitting region. Finally, we show results from a numerical simulation of large-scale structure formation including acceleration of CR electrons at cosmological shocks and magnetic field evolution. Our results are in accord with the observed radio synchrotron and X-ray thermal bremsstrahlung fluxes. Thus we conclude that the reported radio detection is the first evidence of cosmic-ray particle acceleration taking place at cosmic shocks in a magnetized inter-galactic medium over scales of ≳5 h₅₀⁻¹ Mpc.     

Oscillatory cosmic-ray shock structures

A multiple scales analysis is used to derive a mixed Burgers-Korteweg-de Vries (BKdV) equation in the long wavelength regime for a two-fluid MHD model used to describe cosmic-ray acceleration by the first-order Fermi process in astrophysical shocks. The BKdV equation describes the time evolution of weak shocks in the theory of diffusive shock acceleration for all possible cosmic-ray pressures. Previous work on weak shocks in the cosmic-ray MHD model has assumed that dissipation alone is sufficient to balance nonlinearity, but, as cosmic-ray pressures become small, the weak shock becomes discontinous. By including Hall current effects into the MHD model, the low cosmic-ray pressure limit leads smoothly into solitary wave behaviour. For low cosmic-ray pressures, the shock has a downstream oscillatory precursor which is smoothed into the standard Taylor shock profile with increasing cosmic-ray pressure. As a by-product of the perturbation analysis, a dissipative KdV equation is derived. In conclusion, dispersive effects on Alfvén waves are discussed and a modulational stability analysis is presented.     

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.     

Energetic protons associated with a forward-reverse interplanetary shock pair at 1 A.U.

A forward-reverse interplanetary shock was observed on 25 March 1969 by the magnetometer and plasma detector on the HEOS-1 satellite. This relatively rare event was described by Chao et al (1972) who concluded that the shock pair was formed at a distance 0.10–0.13 A.U. upstream of the Earth as a result of the interaction between a fast and a slow solar wind streams. Simultaneous observations of 1 MeV solar proton fluxes were also performed on HEOS-1. A characteristic intensity peak was observed as the forward shock passed by the spacecraft. The evolution of the proton intensity, together with a detailed analysis of anisotropies and pitch angle distributions show a complex dynamic picture of the effect of the forward shock on the ambient proton population. Significant changes in particle fluxes are seen to be correlated with fluctuations in the magnetic field. It is suggested that simple geometrical models of shock-associated acceleration should be expanded to include the effect of magnetic fluctuations on particle fluxes. The interaction region limited by the forward and reverse shocks contained a large variety of magnetic fluctuations. Following the tangential discontinuity separating the fast solar wind stream from the preceding slow stream, a sunward flow was observed in the proton data, followed by a small but significant drop in intensity prior to the reverse shock.     

宇宙线的超新星遗迹起源

宇宙线的起源是高能天体物理的核心问题之一.一直以来,超新星爆发被认为是能谱膝区以下宇宙线的主要来源.多波段观测表明,超新星遗迹有能力加速带电粒子至亚PeV (10~(15)eV)能量.扩散激波加速被认为是最有效的天体高能粒子加速机制之一,而超新星遗迹的大尺度激波正好为这一机制提供平台.近年来,一系列较高精度的地面和空间实验极大地推动了对宇宙线以及超新星遗迹的研究.新的观测事实挑战着传统的扩散激波加速模型以及其在银河系宇宙线超新星遗迹起源学说上的应用,深化了人们对宇宙高能现象的认识.结合超新星遗迹辐射能谱的时间演化特性,构建的时间依赖的超新星遗迹粒子加速模型,不仅能够解释200 GV附近宇宙线的能谱反常,还自然地形成能谱膝区,甚至可以将超新星遗迹粒子加速对宇宙线能谱的贡献延伸至踝区.该模型预期超新星遗迹中粒子的输运行为表现为湍流扩散,这需要未来的观测以及与粒子输运相关的等离子体数值模拟工作来进一步验证.     

The physics of particle acceleration by collisionless shocks

Using analytic theory, test-particle simulations, and self-consistent hybrid simulations, we show that quasi-perpendicular shocks—those which propagate nearly perpendicular to the upstream magnetic field—accelerate particles directly out of the incident thermal population to energies much higher than the upstream ram energy of the plasma. It has already been established that quasi-parallel shocks—those which propagate nearly in the same direction as the upstream magnetic field—efficiently accelerate particles directly out of the incident thermal population; however, this has not yet been established for quasi-perpendicular shocks. Our results can be understood within the framework of the diffusive shock acceleration theory. We find that the accelerated-particle spectrum obtained from a more-general self-consistent hybrid plasma simulation are quantitatively consistent with a less-sophisticated test-particle simulation. The implications of this are discussed.     

On the Challenges of Cosmic-Ray Proton Shock Acceleration in the Intracluster Medium

Galaxy clusters host the largest particle accelerators in the Universe: Shock waves in the intracluster medium (ICM), a hot and ionised plasma, that accelerate particles to high energies. Radio observations pick up synchrotron emission in the ICM, proving the existence of accelerated cosmic-ray electrons. However, a sign of cosmic-ray protons, in form of $\gamma$-rays. remains undetected. This is know as the missing $\gamma$-ray problem and it directly challenges the shock acceleration mechanism at work in the ICM.Over the last decade, theoretical and numerical studies focused on improving our knowledge on the microphysics that govern the shock acceleration process in the ICM. These new models are able to predict a $\gamma$-ray signal, produced by shock accelerated cosmic-ray protons, below the detection limits set modern $\gamma$-ray observatories. In this review, we summarise the latest advances in solving the missing $\gamma$-ray problem.     

Critical Self-Organization of Astrophysical Shocks

There are two distinct regimes of the first-order Fermi acceleration of shocks. The first is a linear (test-particle) regime in which most of the shock energy goes into thermal and bulk motions of the plasma. The second is an efficient regime in which the shock energy goes into accelerated particles. Although the transition region between them is narrow, we identify the factors that drive the system toward a self-organized critical state between those two regimes. Using an analytic solution, we determine this critical state and calculate the spectra and maximum energy of accelerated particles.     

Cosmic Rays and Nonthermal Radiation in Middle-Aged Supernova Remnants

A nonlinear model of cosmic-ray acceleration at the shock fronts in the supernova remnants W28, W44, and IC433 is investigated. The hydrodynamic evolution of a supernova remnant, including the shock modification by the pressure of accelerated particles and the streaming instability of particles upstream of the shock propagating in a partially ionized interstellar gas, is modeled. The electromagnetic radiation generated by accelerated particles is calculated and compared with observations in a wide range of photon energies.     

A model of polarized X-ray emission from twinkling synchrotron supernova shells

Synchrotron X-ray emission components were recently detected in many young supernova remnants (SNRs). There is even an emerging class – SN 1006, RX J1713.72−3946, Vela Jr and others – that is dominated by non-thermal emission in X-rays, also probably of synchrotron origin. Such emission results from electrons/positrons accelerated well above TeV energies in the spectral cut-off regime. In the case of diffusive shock acceleration, which is the most promising acceleration mechanism in SNRs, very strong magnetic fluctuations with amplitudes well above the mean magnetic field must be present. Starting from such a fluctuating field, we have simulated images of polarized X-ray emission of SNR shells and show that these are highly clumpy with high polarizations up to 50 per cent. Another distinct characteristic of this emission is the strong intermittency, resulting from the fluctuating field amplifications. The details of this 'twinkling' polarized X-ray emission of SNRs depend strongly on the magnetic field fluctuation spectra, providing a potentially sensitive diagnostic tool. We demonstrate that the predicted characteristics can be studied with instruments that are currently being considered. These can give unique information on magnetic field characteristics and high-energy particle acceleration in SNRs.     

Cosmic Ray Production Rates in Supernova Remnants

Supernova Remnants (SNRs) are the most likely sources of the galactic cosmic rays up to energies of about 10¹⁵ eV/nuc. The large scale shock waves of SNRs are almost ideal sites to accelerate particles up to these highly non-thermal energies by a first order Fermi mechanism which operates through scattering of the particles at magnetic irregularities. In order to get an estimate on the total amount of the explosion energy E _SNconverted into high energy particles the evolution of a SNR has to be followed up to the final merging with the interstellar medium. This can only be done by numerical simulations since the non-linear modifications of the shock wave due to particle acceleration as well as radiative cooling processes at later SNR stages have to be considered in such investigations. Based on a large sample of numerical evolution calculations performed for different ambient densities n _ext, SN explosion energies, magnetic fields etc. we discuss the final ‘yields’ of cosmic rays at the final SNR stage where the Mach number of the shock waves drops below 2. At these times the cosmic rays start to diffuse out of the remnant. In the range of external densities of10^-2 ≤ n _ext/[cm^-3] ≤ 30 we find a the total acceleration efficiency of about 0.15 E _SN with an increase up to 0.24 E _SN at maximum for an external density of n _ext = 10 cm^-3. Since for the larger ambient densities radiative cooling can reduce significantly the total thermal energy content of the remnant dissipation of Alfvén waves can provide an important heating mechanism for the gas at these later stages. From the collisions of the cosmic rays with the thermal plasma neutral pions are generated which decay subsequently into observable γ-rays above 100 MeV. Hence, we calculate these γ-ray luminosities of SNRs and compare them with current upper limits of ground based γ-raytelescopes. The development of dense shells due to cooling of the thermal plasma increases the γ-ray luminosities and e.g. an external density of n _ext = 10 cm^-3 with E _SN = 10⁵¹ erg can lead to a γ-ray flux above 10^-6 ph cm^-2 s^-1 for a remnant located at a distance of 1 kpc. This revised version was published online in July 2006 with corrections to the Cover Date.     

Effect of Collisions and Magnetic Convergence on Electron Acceleration and Transport in Reconnecting Twisted Solar Flare Loops

We study a model of particle acceleration coupled with an MHD model of magnetic reconnection in unstable twisted coronal loops. The kink instability leads to the formation of helical currents with strong parallel electric fields resulting in electron acceleration. The motion of electrons in the electric and magnetic fields of the reconnecting loop is investigated using a test-particle approach taking into account collisional scattering. We discuss the effects of Coulomb collisions and magnetic convergence near loop footpoints on the spatial distribution and energy spectra of high-energy electron populations and possible implications on the hard X-ray emission in solar flares.     

Synchrotron X-ray diagnostics of cutoff shape of nonthermal electron spectrum at young supernova remnants

Synchrotron X-rays can be a useful tool to investigate electron acceleration at young supernova remnants(SNRs).At present,since the magnetic field configuration around the shocks of SNRs is uncertain,it is not clear whether electron acceleration is limited by SNR age,synchrotron cooling,or even escape from the acceleration region.We study whether the acceleration mechanism can be constrained by the cutoff shape of the electron spectrum around the maximum energy.We derive analytical formulae of the cutoff shape in each case where the maximum electron energy is determined by SNR age,synchrotron cooling and escape from the shock.They are related to the energy dependence of the electron diffusion coefficient.Next,we discuss whether information on the cutoff shape can be provided by observations in the near future which will simply give the photon indices and the flux ratios in the soft and hard X-ray bands.We find that if the power-law index of the electron spectrum is independently determined by other observations,then we can constrain the cutoff shape by comparing theoretical predictions of the photon indices and/or the flux ratios with observed data which will be measured by NuSTAR and/or ASTRO-H.Such study is helpful in understanding the acceleration mechanism.In particular,it will supply another independent constraint on the magnetic field strength around the shocks of SNRs.     

A possible explanation of photon emission from supernova remnants by jitter radiation

We investigate a possibility that non-thermal X-ray emission in a supernova remnant(SNR) is produced by jitter radiation, which is the analogue of synchrotron radiation in small-scale random magnetic fields. We can fit the multi-wavelength data of SNRs RX J1713.7-3946 (G347.3-0.5) and RX J0852.0-4622 (G266.6-1.2) by constructing pure jitter and inverse Compton (IC) emission models. We find that the physical fit parameters of random magnetic fields take values of several tens of μG strength and of the order of ∼10⁷ cm correlation length. These properties of random magnetic fields in collisionless shock of SNRs are discussed.      

The 26 December 2001 Solar Eruptive Event Responsible for GLE63: III. CME,Shock Waves,and Energetic Particles

The SOL2001-12-26 moderate solar eruptive event (GOES importance M7.1, microwaves up to 4000 sfu at 9.4 GHz, coronal mass ejection (CME) speed 1446 km?s^?1) produced strong fluxes of solar energetic particles and ground-level enhancement (GLE) of cosmic-ray intensity (GLE63). To find a possible reason for the atypically high proton outcome of this event, we study multi-wavelength images and dynamic radio spectra and quantitatively reconcile the findings with each other. An additional eruption probably occurred in the same active region about half an hour before the main eruption. The latter produced two blast-wave-like shocks during the impulsive phase. The two shock waves eventually merged around the radial direction into a single shock traced up to \(25~\mathrm{R}_{\odot}\) as a halo ahead of the expanding CME body, in agreement with an interplanetary Type II event recorded by the Radio and Plasma Wave Investigation (WAVES) experiment on the Wind spacecraft. The shape and kinematics of the halo indicate an intermediate regime of the shock between the blast wave and bow shock at these distances. The results show that i) the shock wave appeared during the flare rise and could accelerate particles earlier than usually assumed; ii) the particle event could be amplified by the preceding eruption, which stretched closed structures above the developing CME, facilitated its lift-off and escape of flare-accelerated particles, enabled a higher CME speed and stronger shock ahead; iii) escape of flare-accelerated particles could be additionally facilitated by reconnection of the flux rope, where they were trapped, with a large coronal hole; and iv) the first eruption supplied a rich seed population accelerated by a trailing shock wave.     

Nucleosynthesis and Mixing in Cassiopeia A

We present results from the first light observations of the Cassiopeia A supernova remnant (SNR) by the Chandra X-Ray Observatory. Based on representative spectra from four selected regions, we investigate the processes of nucleosynthesis and mixing in Cas A. We make the first unequivocal identification of iron-rich ejecta produced by explosive silicon burning in a young Galactic SNR. Elsewhere in the remnant, we see silicon-rich ejecta from explosive oxygen burning. The Fe-rich ejecta lie outside the Si-rich material, indicating that bulk motions were extensive and energetic enough in Cas A to cause a spatial inversion of a significant portion of the supernova core. It is likely that this inversion was caused by "Fe"-rich ejecta emerging in plumes from the rising bubbles in the neutrino-driven convection layer during the supernova explosion. In addition, the radioactive decay energy from 56Ni may have contributed to the subsequent evolution of the material. We have also discovered faint, well-defined filaments with featureless X-ray spectra that are possibly sites of cosmic-ray acceleration in Cas A.