Supernova remnants (SNRs) are believed to be the main sources of galactic cosmic rays. Discovery of the non-thermal component in X-ray spectrum of SN 1006 in 1995 and detection of a number of SNRs by H.E.S.S. strengthen the investigation of SNRs. SN 1006 remains to be one of the most interesting objects for high-energy astrophysics. Electrons accelerated by the shock are the source of the non-thermal radiation in radio, X-rays (via synchrotron emission) and \gamma-rays (via inverse-Compton process). Experimental images of SN 1006 are known in all these bands, including the very-high energy \gamma-ray range. An important task is therefore to model the distribution of the surface brightness in SNRs. We develop a method for synthesis of SNR maps due to the non-thermal radiation of electrons in radio, X-rays and \gamma-rays. In particular, the method takes into account the injection of particles and the behaviour of magnetic field at shocks with different obliquities as well as the radiation losses of electrons downstream of the shock. The method is used to model images of SN 1006. The images well correlate with the observations in radio and X-ray but not in \gamma-ray range. For coincide in \gamma-ray range it is needed to increases parameters \eta and E_{f}.
Deep Dive into Non-thermal images of SN 1006: from radio to gamma-rays.
Supernova remnants (SNRs) are believed to be the main sources of galactic cosmic rays. Discovery of the non-thermal component in X-ray spectrum of SN 1006 in 1995 and detection of a number of SNRs by H.E.S.S. strengthen the investigation of SNRs. SN 1006 remains to be one of the most interesting objects for high-energy astrophysics. Electrons accelerated by the shock are the source of the non-thermal radiation in radio, X-rays (via synchrotron emission) and \gamma-rays (via inverse-Compton process). Experimental images of SN 1006 are known in all these bands, including the very-high energy \gamma-ray range. An important task is therefore to model the distribution of the surface brightness in SNRs. We develop a method for synthesis of SNR maps due to the non-thermal radiation of electrons in radio, X-rays and \gamma-rays. In particular, the method takes into account the injection of particles and the behaviour of magnetic field at shocks with different obliquities as well as the radiati
It is common believe that SNRs are the main sources of cosmic rays at least up to the knee at 3 × 10 15 eV and even beyond (Berezhko & V ölk, 2006), (Blandford & Eichler, 1987). Electrons accelerated at the shock of SNRs radiate in wide range from radio to γ-rays. The SN1006 is the prime example of such an object. This object is a type of Ia SNR located at a distance of 2.2 kpc (Winkler et al., 2003).
The remnant of SN 1006 was first identified in the radio (Gurdner & Miln, 1965). In 1982, the first paper went describing were published about describe the first radio images of SN 1006 and Tycho (Reynolds & Chavalier, 1982).
Synchrotron emission from accelerated electrons is detected up to the X-rays in SN 1006 (Rothenflug et al., 2004). In 1995, the first papers went out about domination of non-thermal X-ray component in some regions of SN 1006 (Koyama et al., 1995).
Investigation of the high energy γ-rays was the next step. The first observations in this range are done by the CANGAROO telescope (Tanimori et al., 1998). The first images of SN 1006 as a TeV source was reported by H.E.S.S. team in 2008.
It is suitable to model this object because it has the simple symmetrical distribution of surface brightness in all ranges. Almost spherical form of SN 1006 provide evidences that it develops in a uniform interstellar medium (ISM). Therefore, we model SN 1006 and synthesize it is images in different bands as SNR in uniform ISM and uniform interstellar magnetic field (ISMF).
SN 1006 is on the adiabatic stage of evolution. Therefore, for description of hydrodynamics we use the Sedov self-similar solutions (Sedov, 1959). The calculations of the magnetic field in the SNR interier and on the shock, are bathed on the model of Reynolds (1998). We use the classic 1 Beshley.vasyl@gmail.com, Institute for Applied Problems in Mechanics and Mathematics model of kinetics of the charged particles. The spectrum of relativistic electrons is power-law with the exponential cut-off:
where s is a spectral index and E max the maximum energy of electrons. Distribution of spectrum of electrons is downstream of the shock is described in paper (Reynolds, 1998) and generalized in (Petruk & Beshley, 2008).
It is also needed to take into account the azimuth distribution of the maximal energy. It is not isotropic as it is given by observations (Rothenflug et al., 2004). The Reynolds developed three models of maximum energy: 1)by the electron radiative losses (due to synchrotron and inverse Compton, loss-limited), 2) by the limited time of the acceleration (time-limited) and 3) by properties of microphysics (escape-limited). For description of distribution of maximum energy we use the second model (timelimited) with parameter η, defined as λ = ηr g , where λ is the mean free path and r g -gyroradius and η the “gyrofactor”.
Reynolds (1998) considered thee models for injection: quasi-parallel (qpar), isotropic (iso), quasiperpendicular (qper). Morphology of SNR under quasiperpendicular injection is similar to morphology for isotropic case (Orlando, 2007). Experimental data reject quasi-perpendicular model (Petruk, 2008), if the SNR is assumed to be in uniform ISM. Therefore we consider only isotropic injection.
Spectral index. With the use of radio observation it is possible to find the value of spectral index. Spectral index in this case is s = 2.2 (Rothenflug et al., 2004).
Radiative losses. For the account of radiation losses of electrons we use the parameter
is the reduced fiducial energy, where t is the age of SNR, B the magnetic field strength after parallel shock. This parameter is a measure of importable of radiative losses in modification of the electron spectrum (Reynolds, 1998). If 0 < E f < 1 then the radiative losses are very important if E f > 1 radiative losses may be neglected. In our calculations
Other parameter is orientations of ISMF in respect to the observer (φ o ). Aspect angel equals 70 o (Petruk, 2008). In the present paper we consider two mechanisms of radiation of electrons. This is synchrotron (Petruk & Beshley , 2007), (Petruk & Beshley, 2008) end inverse Compton(IC) (Petruk, 2009).
The surface brightness at an arbitrary point of SNR is calculated as integral along the line of sight, taking into account distributions of the magnetic field, density of electrons etc.
Using results of observation we impose limitation on the parameters of model.
Maximum energy. Distribution of the maximum energy is taken from the observations of satellite XMM-Newton. For different azimuthal angle, E max is different. Observations agree well with time-limited model of Reynolds (1998) with η = 1.
Using the classic model of kinetics of electrons the maps of surface brightness distribution were built for three ranges: radio, X-ray and γ-ray. The model and observational images of SN 1006 are shown on Fig. 1. The limbs of brightness in the simulated images coincide in radio and X-ray ranges. However we have no coincidence in
…(Full text truncated)…
This content is AI-processed based on ArXiv data.