New Method for Determination of Diffusion Coefficients in Turbulent Plasmas

Peer reviewe

Detection of variable VHE gamma-ray emission from the extra-galactic gamma-ray binary LMC P3

Context. Recently, the high-energy (HE, 0.1-100 GeV) $\gamma$-ray emission from the object LMC P3 in the Large Magellanic Cloud (LMC) has been discovered to be modulated with a 10.3-day period, making it the first extra-galactic $\gamma$-ray binary. Aims. This work aims at the detection of very-high-energy (VHE, >100 GeV) $\gamma$-ray emission and the search for modulation of the VHE signal with the orbital period of the binary system. Methods. LMC P3 has been observed with the High Energy Stereoscopic System (H.E.S.S.); the acceptance-corrected exposure time is 100 h. The data set has been folded with the known orbital period of the system in order to test for variability of the emission. Energy spectra are obtained for the orbit-averaged data set, and for the orbital phase bin around the VHE maximum. Results. VHE $\gamma$-ray emission is detected with a statistical significance of 6.4 $\sigma$. The data clearly show variability which is phase-locked to the orbital period of the system. Periodicity cannot be deduced from the H.E.S.S. data set alone. The orbit-averaged luminosity in the $1-10$ TeV energy range is $(1.4 \pm 0.2) \times 10^{35}$ erg/s. A luminosity of $(5 \pm 1) \times 10^{35}$ erg/s is reached during 20% of the orbit. HE and VHE $\gamma$-ray emissions are anti-correlated. LMC P3 is the most luminous $\gamma$-ray binary known so far.Comment: 5 pages, 3 figures, 1 table, accepted for publication in A&

Characterizing the gamma-ray long-term variability of PKS 2155-304 with H.E.S.S. and Fermi-LAT

Studying the temporal variability of BL Lac objects at the highest energies provides unique insights into the extreme physical processes occurring in relativistic jets and in the vicinity of super-massive black holes. To this end, the long-term variability of the BL Lac object PKS 2155-304 is analyzed in the high (HE, 100 MeV 200 GeV) gamma-ray domain. Over the course of ~9 yr of H.E.S.S observations the VHE light curve in the quiescent state is consistent with a log-normal behavior. The VHE variability in this state is well described by flicker noise (power-spectral-density index {\ss}_VHE = 1.10 +0.10 -0.13) on time scales larger than one day. An analysis of 5.5 yr of HE Fermi LAT data gives consistent results ({\ss}_HE = 1.20 +0.21 -0.23, on time scales larger than 10 days) compatible with the VHE findings. The HE and VHE power spectral densities show a scale invariance across the probed time ranges. A direct linear correlation between the VHE and HE fluxes could neither be excluded nor firmly established. These long-term-variability properties are discussed and compared to the red noise behavior ({\ss} ~ 2) seen on shorter time scales during VHE-flaring states. The difference in power spectral noise behavior at VHE energies during quiescent and flaring states provides evidence that these states are influenced by different physical processes, while the compatibility of the HE and VHE long-term results is suggestive of a common physical link as it might be introduced by an underlying jet-disk connection.Comment: 11 pages, 16 figure

Detailed spectral and morphological analysis of the shell type SNR RCW 86

Aims: We aim for an understanding of the morphological and spectral properties of the supernova remnant RCW~86 and for insights into the production mechanism leading to the RCW~86 very high-energy gamma-ray emission. Methods: We analyzed High Energy Spectroscopic System data that had increased sensitivity compared to the observations presented in the RCW~86 H.E.S.S. discovery publication. Studies of the morphological correlation between the 0.5-1~keV X-ray band, the 2-5~keV X-ray band, radio, and gamma-ray emissions have been performed as well as broadband modeling of the spectral energy distribution with two different emission models. Results:We present the first conclusive evidence that the TeV gamma-ray emission region is shell-like based on our morphological studies. The comparison with 2-5~keV X-ray data reveals a correlation with the 0.4-50~TeV gamma-ray emission.The spectrum of RCW~86 is best described by a power law with an exponential cutoff at $E_{cut}=(3.5\pm 1.2_{stat})$ TeV and a spectral index of $\Gamma$~$1.6\pm 0.2$. A static leptonic one-zone model adequately describes the measured spectral energy distribution of RCW~86, with the resultant total kinetic energy of the electrons above 1 GeV being equivalent to $\sim$0.1\% of the initial kinetic energy of a Type I a supernova explosion. When using a hadronic model, a magnetic field of $B$~100$\mu$G is needed to represent the measured data. Although this is comparable to formerly published estimates, a standard E$^{-2}$ spectrum for the proton distribution cannot describe the gamma-ray data. Instead, a spectral index of $\Gamma_p$~1.7 would be required, which implies that ~$7\times 10^{49}/n_{cm^{-3}}$erg has been transferred into high-energy protons with the effective density $n_{cm^{-3}}=n/ 1$ cm^-3. This is about 10\% of the kinetic energy of a typical Type Ia supernova under the assumption of a density of 1~cm^-3.Comment: accepted for publication by A&

Determining pitch-angle diffusion coefficients from test particle simulations

Transport and acceleration of charged particles in turbulent media is a topic of great interest in space physics and interstellar astrophysics. These processes are dominated by the scattering of particles off magnetic irregularities. The scattering process itself is usually described by small-angle scattering with the pitch-angle coefficient D μμ   playing a major role. Since the diffusion coefficient D μμ   can be determined analytically only for the approximation of quasi-linear theory, the determination of this coefficient from numerical simulations has, therefore, become more important. So far these simulations yield particle tracks for small-scale scattering, which can then be interpreted using the running diffusion coefficients. This method has a limited range of validity. This paper presents two new methods that allow for the calculation of the pitch-angle diffusion coefficient from numerical simulations. These methods no longer analyse particle trajectories, but the change of particle distribution functions. It is shown that they provide better resolved results and allow for the analysis of strong turbulence. The application of these methods to Monte Carlo simulations of particle scattering and hybrid MHD-particle simulations is presented. Both analysis methods are able to recover the diffusion coefficients used as input for the Monte Carlo simulations and provide better results in MHD simulations especially for stronger turbulence. </p

Multi-messenger observations of a binary neutron star merger

On 2017 August 17 a binary neutron star coalescence candidate (later designated GW170817) with merger time 12:41:04 UTC was observed through gravitational waves by the Advanced LIGO and Advanced Virgo detectors. The Fermi Gamma-ray Burst Monitor independently detected a gamma-ray burst (GRB 170817A) with a time delay of ~1.7 s with respect to the merger time. From the gravitational-wave signal, the source was initially localized to a sky region of 31 deg2 at a luminosity distance of 40+8-8 Mpc and with component masses consistent with neutron stars. The component masses were later measured to be in the range 0.86 to 2.26 Mo. An extensive observing campaign was launched across the electromagnetic spectrum leading to the discovery of a bright optical transient (SSS17a, now with the IAU identification of AT 2017gfo) in NGC 4993 (at ~40 Mpc) less than 11 hours after the merger by the One- Meter, Two Hemisphere (1M2H) team using the 1 m Swope Telescope. The optical transient was independently detected by multiple teams within an hour. Subsequent observations targeted the object and its environment. Early ultraviolet observations revealed a blue transient that faded within 48 hours. Optical and infrared observations showed a redward evolution over ~10 days. Following early non-detections, X-ray and radio emission were discovered at the transient’s position ~9 and ~16 days, respectively, after the merger. Both the X-ray and radio emission likely arise from a physical process that is distinct from the one that generates the UV/optical/near-infrared emission. No ultra-high-energy gamma-rays and no neutrino candidates consistent with the source were found in follow-up searches. These observations support the hypothesis that GW170817 was produced by the merger of two neutron stars in NGC4993 followed by a short gamma-ray burst (GRB 170817A) and a kilonova/macronova powered by the radioactive decay of r-process nuclei synthesized in the ejecta

The 2012 flare of PG 1553+113 seen with H.E.S.S. and Fermi-LAT

Very high energy (VHE, E > 100 GeV) γ-ray flaring activity of the high-frequency peaked BL Lac object PG 1553+113 has been detected by the H.E.S.S. telescopes. The flux of the source increased by a factor of 3 during the nights of 2012 April 26 and 27 with respect to the archival measurements with a hint of intra-night variability. No counterpart of this event has been detected in the Fermi-Large Area Telescope data. This pattern is consistent with VHE γ-ray flaring being caused by the injection of ultrarelativistic particles, emitting γ-rays at the highest energies. The dataset offers a unique opportunity to constrain the redshift of this source at z = 0.49 ± 0.04 using a novel method based on Bayesian statistics. The indication of intra-night variability is used to introduce a novel method to probe for a possible Lorentz invariance violation (LIV), and to set limits on the energy scale at which Quantum Gravity (QG) effects causing LIV may arise. For the subluminal case, the derived limits are EQG,1 > 4.10 × 1017 GeV and EQG,2 > 2.10 × 1010 GeV for linear and quadratic LIV effects, respectivel

The population of TeV pulsar wind nebulae in the H.E.S.S. Galactic Plane Survey

The nine-year H.E.S.S. Galactic Plane Survey (HGPS) has yielded the most uniform observation scan of the inner Milky Way in the TeV gamma-ray band to date. The sky maps and source catalogue of the HGPS allow for a systematic study of the population of TeV pulsar wind nebulae found throughout the last decade. To investigate the nature and evolution of pulsar wind nebulae, for the first time we also present several upper limits for regions around pulsars without a detected TeV wind nebula. Our data exhibit a correlation of TeV surface brightness with pulsar spin-down power Ė. This seems to be caused both by an increase of extension with decreasing Ė, and hence with time, compatible with a power law RPWN(Ė) ~Ė−0.65±0.20, and by a mild decrease of TeV gamma-ray luminosity with decreasing Ė, compatible with L1−10 TeV ~Ė0.59±0.21. We also find that the offsets of pulsars with respect to the wind nebula centre with ages around 10 kyr are frequently larger than can be plausibly explained by pulsar proper motion and could be due to an asymmetric environment. In the present data, it seems that a large pulsar offset is correlated with a high apparent TeV efficiency L1−10 TeV∕Ė. In addition to 14 HGPS sources considered firmly identified pulsar wind nebulae and 5 additional pulsar wind nebulae taken from literature, we find 10 HGPS sources that are likely TeV pulsar wind nebula candidates. Using a model that subsumes the present common understanding of the very high-energy radiative evolution of pulsar wind nebulae, we find that the trends and variations of the TeV observables and limits can be reproduced to a good level, drawing a consistent picture of present-day TeV data and theor

Population study of Galactic supernova remnants at very high γ-ray energies with H.E.S.S.

Shell-type supernova remnants (SNRs) are considered prime candidates for the acceleration of Galactic cosmic rays (CRs) up to the knee of the CR spectrum at E ≈ 3 × 1015 eV. Our Milky Way galaxy hosts more than 350 SNRs discovered at radio wavelengths and at high energies, of which 220 fall into the H.E.S.S. Galactic Plane Survey (HGPS) region. Of those, only 50 SNRs are coincident with a H.E.S.S source and in 8 cases the very high-energy (VHE) emission is firmly identified as an SNR. The H.E.S.S. GPS provides us with a legacy for SNR population study in VHE γ-rays and we use this rich data set to extract VHE flux upper limits from all undetected SNRs. Overall, the derived flux upper limits are not in contradiction with the canonical CR paradigm. Assuming this paradigm holds true, we can constrain typical ambient density values around shell-type SNRs to n ≤ 7 cm−3 and electron-to-proton energy fractions above 10 TeV to ϵep ≤ 5 × 10−3. Furthermore, comparisons of VHE with radio luminosities in non-interacting SNRs reveal a behaviour that is in agreement with the theory of magnetic field amplification at shell-type SNR

The supernova remnant W49B as seen with H.E.S.S. and Fermi-LAT

The supernova remnant (SNR) W49B originated from a core-collapse supernova that occurred between one and four thousand years ago, and subsequently evolved into a mixed-morphology remnant, which is interacting with molecular clouds (MC). Gamma-ray observations of SNR-MC associations are a powerful tool to constrain the origin of Galactic cosmic rays, as they can probe the acceleration of hadrons through their interaction with the surrounding medium and subsequent emission of non-thermal photons. We report the detection of a γ-ray source coincident with W49B at very high energies (VHE; E > 100 GeV) with the H.E.S.S. Cherenkov telescopes together with a study of the source with five years of Fermi-LAT high-energy γ-ray (0.06–300 GeV) data. The smoothly connected, combined source spectrum, measured from 60 MeV to multi-TeV energies, shows two significant spectral breaks at 304 ± 20 MeV and 8.4−2.5+2.2 GeV; the latter is constrained by the joint fit from the two instruments. The detected spectral features are similar to those observed in several other SNR-MC associations and are found to be indicative of γ-ray emission produced through neutral-pion deca