The discovery of direct evidences for the acceleration of high energetic particles at the shell supernova remnant RXJ1713.7-3946 underlined the need to calculate the cosmic ray (CR) distribution in ...the Galaxy on a spatial grid fine enough to resolve the changes in the CR density due to these kind of objects. It was shown before by Büsching et al. Büsching, I., Kopp, A., Pohl, M., Schlickeiser, R., Perrot, C., Grenier, I. Cosmic-ray propagation properties for an origin in supernova remnants. ApJ 619, 314–326, 2005 that the discrete nature (both in space and time) of super novae (SN) as sources of Galactic CR leads to CR spectra changing in space and time, resulting in a range of possible CR spectra at a given location in the Galaxy. As the most frequent SN types Ib and II are found within spiral arms, one can expect a significant difference of the range of possible spectra in and outside spiral arms. We investigate the variability of the local interstellar CR proton spectrum during the motion of the Sun in and out of spiral arms in its journey around the Galactic center. Using the code described by Büsching et al. Büsching, I., Kopp, A., Pohl, M., Schlickeiser, R., Perrot, C., Grenier, I. Cosmic-ray propagation properties for an origin in supernova remnants. ApJ 619, 314–326, 2005, the proton CR density in the Galaxy is calculated with high spatial and temporal resolution (75
pc in galactocentric radius
r and azimuth
φ at the position of the Sun and 20
pc in
z (perpendicular to the galactic plane), with a time step of 1
kyr), assuming stochastically distributed point sources with a probability distribution that resembles the spiral structure of our Galaxy. We find the averaged CR flux outside the spiral arm to be about 50% of that inside a spiral arm. We further find spatial and temporal variations of the CR flux inside spiral arms.
We present calculations of the propagation times and energy losses of cosmic rays as they are transported through the heliosphere. By calculating these quantities for a spatially 1D scenario, we ...benchmark our numerical model, which uses stochastic differential equations to solve the relevant transport equation, with known analytical solutions. The comparison is successful and serves as a vindication of the modeling approach. A spatially 3D version of the modulation model is subsequently used to calculate the propagation times and energy losses of galactic electrons and protons in different drift cycles. We find that the propagation times of electrons are longer than those of the protons at the same energy. Furthermore, the propagation times are longer in the drift cycle when the particles reach the Earth by drifting inward along the heliospheric current sheet. The calculated energy losses follow this same general trend. The energy losses suffered by the electrons are comparable to those of the protons, which is in contrast to the generally held perception that electrons experience little energy losses during their propagation through the heliosphere.
Key Points
Heliospheric propagation times of cosmic ray protons & electrons are calculated
Exact calculations of energy losses for protons and electrons are presented
The SDE approach lends itself to the calculation of these quantities
PSR J0437-4715 is one of the closest millisecond pulsars (MSPs) to Earth, lying at a distance of ∼140
pc. This pulsar has a characteristic age of 4.9
Gyr and a relatively low spindown power of ∼10
34
...ergs/s. During its rather long lifetime, a large fraction of the energy output has been in the form of multi-TeV electrons. In this paper, we investigate the possible contribution of this nearby MSP to the local interstellar electron spectrum (LIS). The old age of the system justifies a steady-state evaluation of the contribution from this pulsar to the LIS. We calculate the electron spectrum at the light cylinder in the framework of a General Relativistic polar cap (PC) model, and use this as an injection spectrum in a diffusion model. The younger Geminga pulsar is also very close to Earth and warrants investigation. A steady-state approach is however no longer justified, so we use an impulsive injection model. We will present results of a study of the contribution from these pulsars to the cosmic ray (CR) LIS. Our calculations show that pulsars like Geminga can make a non-negligible contribution to the LIS.
The measurement of an excess in the cosmic-ray electron spectrum between 300 and 800 GeV by the ATIC experiment has – together with the PAMELA detection of a rise in the positron fraction up to ...≈100 GeV – motivated many interpretations in terms of dark matter scenarios; alternative explanations assume a nearby electron source like a pulsar or supernova remnant. Here we present a measurement of the cosmic-ray electron spectrum with H.E.S.S. starting at 340 GeV. While the overall electron flux measured by H.E.S.S. is consistent with the ATIC data within statistical and systematic errors, the H.E.S.S. data exclude a pronounced peak in the electron spectrum as suggested for interpretation by ATIC. The H.E.S.S. data follow a power-law spectrum with spectral index of 3.0±0.1(stat.)± 0.3(syst.), which steepens at about 1 TeV.
In this study, we investigate how restrictive the γ-ray emission from the Galactic center region, as seen by HESS and other Cherenkov air shower arrays, is against various models for cosmic ray ...injection. We derive diffusion coefficients which fit the observed spatial scales of diffuse γ-ray emission from the extended emission associated with the molecular clouds SgrA, B and C. Using these diffusion coefficients, we then obtain a limit for time scale of assumed recent proton acceleration near the SMBH, as the spatial size of SgrA
∗ in VHE γ-rays has to be consistent with the observed unresolved HESS point source size at this position. The signal from this hadronic component may be mixed with the expected VHE inverse Compton emission from the nearby unresolved pulsar wind nebula.
The detection of fast variations of the tera-electron volt (TeV) (10¹² eV) γ-ray flux, on time scales of days, from the nearby radio galaxy M87 is reported. These variations are about 10 times as ...fast as those observed in any other wave band and imply a very compact emission region with a dimension similar to the Schwarzschild radius of the central black hole. We thus can exclude several other sites and processes of the γ-ray production. The observations confirm that TeV γ rays are emitted by extragalactic sources other than blazars, where jets are not relativistically beamed toward the observer.
We present a newly developed numerical modulation model to study the transport of galactic and Jovian electrons in the heliosphere. The model employs stochastic differential equations (SDEs) to solve ...the corresponding transport equation in five dimensions (time, energy, and three spatial dimensions) which is difficult to accomplish with the numerical schemes used in finite difference models. Modeled energy spectra for galactic electrons are compared for the two drift cycles to observations at Earth. Energy spectra and radial intensity profiles of galactic and Jovian electrons are compared successfully to results from previous studies. In line with general drift considerations, it is found that most 100 MeV electrons observed at Earth enter the heliosphere near the equatorial regions in the A > 0 cycle, while they enter mainly over the polar regions in the A < 0 cycle. Our results indicate that 100 MeV electrons observed at Earth originate at the heliopause with ~600 MeV undergoing adiabatic cooling during their transport to Earth. The mean propagation time of these particles varies between ~180 and 300 days, depending on the drift cycle. For 10 MeV Jovian electrons observed at Earth, a mean propagation time of ~40 days is obtained. During this time, the azimuthal position of the Jovian magnetosphere varies by ~1?. At a 50 AU observational point, the mean propagation time of these electrons increases to ~370 days with an azimuthal position change of Jupiter of ~20?. The SDE approach is very effective in calculating these propagation times.
The recent discovery of diffuse, very high energy (VHE) g-radiation from the Galactic center ridge by the HESS telescope allows for the first time the direct determination of the parameters of ...Galactic cosmic-ray propagation models. Whereas this discovery showed that the diffuse g-radiation can be explained by the interaction of VHE cosmic-ray (CR) protons with the interstellar gas located in several giant molecular clouds, we show in this paper that the associated diffusion coefficient for the protons depends on the epoch of activity of the central source of protons: Assuming that the supernova remnant (SNR) Sgr A East was responsible for the particle acceleration, we infer a diffusion coefficient for the Galactic center region of =1-5 kpc super(2) Myr super(-1) for a mean proton energy of 63 TeV. More specifically, for impulsive injection in a 5-10 kyr SNR, we infer a value of =1-2 kpc super(2) Myr super(-1), whereas for source activity timescales equal to the age of the SNR, the diffusion coefficient would increase to g 5 kpc super(2) Myr super(-1). These values are smaller than those inferred from local CR abundances. Finally, the above-mentioned values of for impulsive injection are equally valid if the required transient source of protons was due to an earlier epoch of stellar infall into the black hole Sgr A*.
Aims. To provide a significantly improved probability distribution for the H-test for periodicity in X-ray and γ-ray arrival times, which is already extensively used by the γ-ray pulsar community. ...Also, to obtain an analytical probability distribution for stacked test statistics in the case of a search for pulsed emission from an ensemble of pulsars where the significance per pulsar is relatively low, making individual detections insignificant on their own. This information is timely given the recent rapid discovery of new pulsars with the Fermi-LAT t γ-ray telescope. Methods. Approximately 1014 realisations of the H-statistic (H) for random (white) noise is calculated from a random number generator for which the repetition cycle is ≫1014. From these numbers the probability distribution P(>H) is calculated. Results. The distribution of H is found to be exponential with parameter λ = 0.4 so that the cumulative probability distribution P(>H)= $\exp$(-λ H). If we stack independent values for H, the sum of K such values would follow the Erlang-K distribution with parameter λ for which the cumulative probability distribution is also a simple analytical expression. Conclusions. Searches for weak pulsars with unknown pulse profile shapes in the Fermi-LAT, Agile or other X-ray data bases should benefit from the H-test since it is known to be powerful against a broad range of pulse profiles, which introduces only a single statistical trial if only the H-test is used. The new probability distribution presented here favours the detection of weaker pulsars in terms of an improved sensitivity relative to the previously known distribution.
Observations of the Crab nebula with HESS Aharonian, F.; Akhperjanian, A. G.; Bazer-Bachi, A. R. ...
Astronomy & astrophysics,
10/2006, Volume:
457, Issue:
3
Journal Article
Peer reviewed
Open access
Context.The Crab nebula was observed with the HESS stereoscopic Cherenkov-telescope array between October 2003 and January 2005 for a total of 22.9 h (after data quality selection). This period of ...time partly overlapped with the commissioning phase of the experiment; observations were made with three operational telescopes in late 2003 and with the complete 4 telescope array in January–February 2004 and October 2004–January 2005. Aims.Observations of the Crab nebula are discussed and used as an example to detail the flux and spectral analysis procedures of HESS. The results are used to evaluate the systematic uncertainties in HESS flux measurements. Methods.The Crab nebula data are analysed using standard HESS analysis procedures, which are described in detail. The flux and spectrum of γ-rays from the source are calculated on run-by-run and monthly time-scales, and a correction is applied for long-term variations in the detector sensitivity. Comparisons of the measured flux and spectrum over the observation period, along with the results from a number of different analysis procedures are used to estimate systematic uncertainties in the measurements. Results.The data, taken at a range of zenith angles between $45^{\circ}$ and $65^{\circ}$, show a clear signal with over 7500 excess events. The energy spectrum is found to follow a power law with an exponential cutoff, with photon index $\Gamma = 2.39$ ± $0.03_{{\rm {stat}}}$ and cutoff energy $E_{\rm c} = (14.3 \pm 2.1_{{\rm {stat}}})~{\rm TeV}$ between 440 GeV and 40 TeV. The observed integral flux above 1 TeV is $(2.26 \pm 0.08_{{\rm {stat}}}) \times 10^{-11}~{{\rm cm}^{-2}~{\rm s}^{-1}}$. The estimated systematic error on the flux measurement is estimated to be 20%, while the estimated systematic error on the spectral slope is 0.1.