-
A search for new supernova remnants (SNRs) has been conducted using TeV
gamma-ray data from the H.E.S.S. Galactic plane survey. As an identification
criterion, shell morphologies that are characteristic for known resolved TeV
SNRs have been used. Three new SNR candidates were identified in the H.E.S.S.
data set with this method. Extensive multiwavelength searches for counterparts
were conducted. A radio SNR candidate has been identified to be a counterpart
to HESS J1534-571. The TeV source is therefore classified as a SNR. For the
other two sources, HESS J1614-518 and HESS J1912+101, no identifying
counterparts have been found, thus they remain SNR candidates for the time
being. TeV-emitting SNRs are key objects in the context of identifying the
accelerators of Galactic cosmic rays. The TeV emission of the relativistic
particles in the new sources is examined in view of possible leptonic and
hadronic emission scenarios, taking the current multiwavelength knowledge into
account.
-
Aims. We study gamma-ray emission from the shell-type supernova remnant (SNR)
RX J0852.0$-$4622 to better characterize its spectral properties and its
distribution over the SNR.
Methods. The analysis of an extended High Energy Spectroscopic System
(H.E.S.S.) data set at very high energies (E > 100 GeV) permits detailed
studies, as well as spatially resolved spectroscopy, of the morphology and
spectrum of the whole RX J0852.0$-$4622 region. The H.E.S.S. data are combined
with archival data from other wavebands and interpreted in the framework of
leptonic and hadronic models. The joint Fermi-LAT-H.E.S.S. spectrum allows the
direct determination of the spectral characteristics of the parent particle
population in leptonic and hadronic scenarios using only GeV-TeV data.
Results. An updated analysis of the H.E.S.S. data shows that the spectrum of
the entire SNR connects smoothly to the high-energy spectrum measured by
Fermi-LAT. The increased data set makes it possible to demonstrate that the
H.E.S.S. spectrum deviates significantly from a power law and is well described
by both a curved power law and a power law with an exponential cutoff at an
energy of $E_\mathrm{cut} = (6.7 \pm 1.2_\mathrm{stat} \pm 1.2_\mathrm{syst})$
TeV. The joint Fermi-LAT-H.E.S.S. spectrum allows the unambiguous
identification of the spectral shape as a power law with an exponential cutoff.
No significant evidence is found for a variation of the spectral parameters
across the SNR, suggesting similar conditions of particle acceleration across
the remnant. A simple modeling using one particle population to model the SNR
emission demonstrates that both leptonic and hadronic emission scenarios remain
plausible. It is also shown that at least a part of the shell emission is
likely due to the presence of a pulsar wind nebula around PSR J0855$-$4644.
-
Recent reports of possible correlations between high energy neutrinos
observed by IceCube and Active Galactic Nuclei (AGN) activity sparked a burst
of publications that attempt to predict the neutrino flux of these sources.
However, often rather crude estimates are used to derive the neutrino rate from
the observed photon spectra. In this work neutrino fluxes were computed in a
wide parameter space. The starting point of the model was a representation of
the full spectral energy density (SED) of \textit{3C 279}. The time-dependent
hybrid model that was used for this study takes into account the full $p\gamma$
reaction chain as well as proton synchrotron, electron-positron-pair cascades
and the full SSC scheme. We compare our results to estimates frequently used in
the literature. This allows to identify regions in the parameter space for
which such estimates are still valid and those in which they can produce
significant errors. Furthermore, if estimates for the Doppler factor, magnetic
field, proton and electron densities of a source exist, the expected IceCube
detection rate is readily available.
-
PKS 0625-354 (z=0.055) was observed with the four H.E.S.S. telescopes in 2012
during 5.5 hours. The source was detected above an energy threshold of 200 GeV
at a significance level of 6.1$\sigma$. No significant variability is found in
these observations. The source is well described with a power-law spectrum with
photon index $\Gamma=2.84\pm0.50_{stat}\pm0.10_{syst}$ and normalization (at
$E_0$=1.0 TeV)
$N_0(E_0)=(0.58\pm0.22_{stat}\pm0.12_{syst})\times10^{-12}$\,TeV$^{-1}$cm$^{-2}$s$^{-1}$.
Multi-wavelength data collected with Fermi-LAT, Swift-XRT, Swift-UVOT, ATOM and
WISE are also analysed. Significant variability is observed only in the
Fermi-LAT $\gamma$-ray and Swift-XRT X-ray energy bands. Having a good
multi-wavelength coverage from radio to very high energy, we performed a
broadband modelling from two types of emission scenarios. The results from a
one zone lepto-hadronic, and a multi-zone leptonic models are compared and
discussed. On the grounds of energetics, our analysis favours a leptonic
multi-zone model. Models associated to the X-ray variability constraint
supports previous results suggesting a BL Lac nature of PKS 0625-354, with,
however, a large-scale jet structure typical of a radio galaxy.
-
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
$\mathrm{E} \approx \mathrm{3}\times \mathrm{10}^\mathrm{15}$ 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
$\gamma$-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\leq 7~\textrm{cm}^\textrm{-3}$ and electron-to-proton energy fractions above
10~TeV to $\epsilon_\textrm{ep} \leq 5\times 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
SNRs.
-
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.
-
The H.E.S.S. collaboration has discovered a new very high energy (VHE, E $>$
0.1 TeV) $\gamma$-ray source, HESS J1741-302, located in the Galactic plane.
Despite several attempts to constrain its nature, no plausible counterpart has
been found so far at X-ray and MeV/GeV $\gamma$-ray energies, and the source
remains unidentified. An analysis of 145-hour of observations of HESS J1741-302
at VHEs has revealed a steady and relatively weak TeV source ($\sim$1$\%$ of
the Crab Nebula flux), with a spectral index of $\Gamma$ = 2.3 $\pm$
0.2$_{\text{stat}}$ $\pm$ 0.2$_{\text{sys}}$, extending to energies up to 10
TeV without any clear signature of a cut-off. In a hadronic scenario, such a
spectrum implies an object with particle acceleration up to energies of several
hundred TeV. Contrary to most H.E.S.S. unidentified sources, the angular size
of HESS J1741-302 is compatible with the H.E.S.S. point spread function at
VHEs, with an extension constrained to be below 0.068$^{\circ}$ at a 99$\%$
confidence level. The $\gamma$-ray emission detected by H.E.S.S. can be
explained both within a hadronic scenario, due to collisions of protons with
energies of hundreds of TeV with dense molecular clouds, and in a leptonic
scenario, as a relic pulsar wind nebula, possibly powered by the middle-aged
(20 kyr) pulsar PSR B1737-30. A binary scenario, related to the compact radio
source 1LC 358.266+0.038 found to be spatially coincident with the best fit
position of HESS J1741-302, is also envisaged.
-
We search for high-energy gamma-ray emission from the binary neutron star
merger GW170817 with the H.E.S.S. Imaging Air Cherenkov Telescopes. The
observations presented here have been obtained starting only 5.3h after
GW170817. The H.E.S.S. target selection identified regions of high probability
to find a counterpart of the gravitational wave event. The first of these
regions contained the counterpart SSS17a that has been identified in the
optical range several hours after our observations. We can therefore present
the first data obtained by a ground-based pointing instrument on this object. A
subsequent monitoring campaign with the H.E.S.S. telescopes extended over
several days, covering timescales from 0.22 to 5.2 days and energy ranges
between $270\,\mathrm{GeV}$ to $8.55\,\mathrm{TeV}$. No significant gamma-ray
emission has been found. The derived upper limits on the very-high-energy
gamma-ray flux for the first time constrain non-thermal, high-energy emission
following the merger of a confirmed binary neutron star system.
-
Index of H.E.S.S. conference proceedings to the 35th ICRC, Busan, Korea
-
Understanding the physics of beam-plasma systems is of great importance for a
large number of astrophysical systems. Recently the validity of
particle-in-cell simulations of these systems has been questioned by Shalaby et
al. (2017) We analyse one specific beam-plasma and its numerical simulations
Kempf et al. (2016) that has been critized and prove that in fact the
simulations performed fulfill the new validity criteria introduced by Shalaby
et al. (2016).
-
Very high-energy gamma-rays (VHE, E>100 GeV) propagating over cosmological
distances can interact with the low-energy photons of the extragalactic
background light (EBL) and produce electron-positron pairs. The transparency of
the universe to VHE gamma-rays is then directly related to the spectral energy
distribution (SED) of the EBL. The observation of features in the VHE energy
spectra of extragalactic sources allows the EBL to be measured, which otherwise
is very difficult to determine. An EBL-model independent measurement of the EBL
SED with the H.E.S.S. array of Cherenkov telescopes is presented. It is
obtained by extracting the EBL absorption signal from the reanalysis of
high-quality spectra of blazars. From H.E.S.S. data alone the EBL signature is
detected at a significance of 9.5 sigma, and the intensity of the EBL obtained
in different spectral bands is presented together with the associated gamma-ray
horizon.
-
The diffuse very high-energy (VHE, > 100 GeV) gamma-ray emission observed in
the central 200 pc of the Milky Way by H.E.S.S. was found to follow the dense
matter distribution in the Central Molecular Zone (CMZ) up to a longitudinal
distance of about 130 pc to the Galactic Centre (GC), where the flux rapidly
decreases. Recent sensitive H.E.S.S. analysis revealed that the cosmic-ray (CR)
density profile drops with the distance to the centre, making data compatible
with a steady cosmic PeVatron at the GC. In this paper, we extend this analysis
to obtain for the first time a detailed characterisation of the correlation
with matter and to search for additional features and individual gamma-ray
sources in the inner 200 pc. Taking advantage of 250 hours of H.E.S.S. data and
improved analysis techniques we perform a detailed morphology study of the
diffuse VHE emission observed from the GC ridge and reconstruct its total
spectrum. We show that the emission correlated with dense matter covers the
full CMZ and that its flux is about half the total diffuse emission flux. We
also detect some emission at higher latitude likely produced by hadronic
collisions of CRs in less dense regions of the GC interstellar medium. We
detect an additional emission component centred on the GC and extending over
about 15 pc that is consistent with the existence of a strong CR density
gradient and confirms the presence of a CR accelerator at the very centre of
our Galaxy. We show that the spectrum of the full ridge diffuse emission is
compatible with the one previously derived from the central regions, suggesting
that a single population of particles fills the entire CMZ. Finally, we report
the discovery of a VHE gamma-ray source near the GC radio arc and argue that it
is produced by the pulsar wind nebula candidate G0.13-0.11.
-
Runaway stars form bow shocks by ploughing through the interstellar medium at
supersonic speeds and are promising sources of non-thermal emission of photons.
One of these objects has been found to emit non-thermal radiation in the radio
band. This triggered the development of theoretical models predicting
non-thermal photons from radio up to very-high-energy (VHE, E $\geq 0.1$TeV)
gamma rays. Subsequently, one bow shock was also detected in X-ray
observations. However, the data did not allow discrimination between a hot
thermal and a non-thermal origin. Further observations of different candidates
at X-ray energies showed no evidence for emission at the position of the bow
shocks either. A systematic search in the \textit{Fermi}-LAT energy regime
resulted in flux upper limits for 27 candidates listed in the E-BOSS catalogue.
Here we perform the first systematic search for VHE gamma-ray emission from bow
shocks of runaway stars. Using all available archival H.E.S.S. data we search
for very-high-energy gamma-ray emission at the positions of bow shock
candidates listed in the second E-BOSS catalogue release. Out of the 73 bow
shock candidates in this catalogue, 32 have been observed with H.E.S.S. None of
the observed 32 bow shock candidates in this population study show significant
emission in the H.E.S.S. energy range. Therefore, flux upper limits are
calculated in five energy bins and the fraction of the kinetic wind power that
is converted into VHE gamma rays is constrained. Emission from stellar bow
shocks is not detected in the energy range between 0.14 and 18 TeV. The
resulting upper limits constrain the level of VHE gamma-ray emission from these
objects down to 0.1-1 % of the kinetic wind energy.
-
The nonlinear fluid theory developed by Schamel suggests a modified KdV
equation to describe the temporal evolution of ion acoustic (IA) solitons in
the presence of trapped electrons. The validity of this theory is studied here
by verifying solitons main characteristic, i.e., stability against successive
mutual collisions. We have employed a kinetic model as a more comprehensive
theory than the fluid one, and utilized a fully kinetic simulation approach
(both ions and electrons are treated based on the Vlasov equation). In the
simulation approach, these solitons are excited self-consistently by employing
the nonlinear process of IA solitons formation from an initial density
perturbation (IDP). The effect of the size of IDPs on the chain formation is
proved by the simulation code as a benchmark test. It is shown that the IA
solitons, in presence of trapped electrons, can retain their features (both in
spatial and velocity direction) after successive mutual collisions. The
collisions here include encounters of IA solitons with the same trapping
parameter, while differing in size. Kinetic simulation results reveal a
complicated behavior during a collision between IA solitons in contrast to the
fluid theory predictions and simulations. In the range of parameters considered
here two oppositely propagating solitons rotate around their collective center
in the phase space during a collision, independent of their trapping
parameters. Furthermore, they exchange some portions of their trapped
populations.
-
The nine-year H.E.S.S. Galactic Plane Survey (HGPS) 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 $\dot{E}$. This seems to be caused both
by an increase of extension with decreasing $\dot{E}$, and hence with time,
compatible with a power law $R_\mathrm{PWN}(\dot{E}) \sim \dot{E}^{-0.65 \pm
0.20}$, and by a mild decrease of TeV gamma-ray luminosity with decreasing
$\dot{E}$, compatible with $L_{1-10\,\mathrm{TeV}} \sim \dot{E}^{0.59 \pm
0.21}$. We also find that the offsets of pulsars with respect to the wind
nebula centres 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
$L_{1-10\,\mathrm{TeV}}/\dot{E}$. In addition to 14 HGPS sources considered as
firmly identified pulsar wind nebulae and 5 additional pulsar wind nebulae
taken from literature, we find 10 HGPS sources that form 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 theory.
-
The overtaking collisions of ion-acoustic solitons (IASs) in presence of
trapping effects of electrons are studied based on a fully kinetic simulation
approach. The method is able to provide all the kinetic details of the process
alongside the fluid-level quantities self consistently. Solitons are produced
naturally by utilizing the chain formation phenomenon, then are arranged in a
new simulation box to test different scenarios of overtaking collisions. Three
achievements are reported here. Firstly, simulations prove the long-time life
span of the ion-acoustic solitons in the presence of trapping effect of
electrons (kinetic effects), which serves as the benchmark of the simulation
code. Secondly, their stability against overtaking mutual collisions is
established by creating collisions between solitons with different number and
shapes of trapped electrons, i.e. different trapping parameter. Finally,
details of solitons during collisions for both ions and electrons are provided
on both fluid and kinetic levels. These results show that on the kinetic level,
trapped electron population accompanying each of the solitons are exchanged
between the solitons during the collision. Furthermore, the behavior of
electron holes accompanying solitons contradicts the theory about the electron
holes interaction developed based on kinetic theory. They also show behaviors
much different from other electron holes witnessed in processes such as
nonlinear Landau damping (Bernstein-Greene-Kruskal -BGK- modes) or beam-plasma
interaction (like two-beam instability).
-
Supernova remnants exhibit shock fronts (shells) that can accelerate charged
particles up to very high energies. In the past decade, measurements of a
handful of shell-type supernova remnants in very-high-energy gamma rays have
provided unique insights into the acceleration process. Among those objects,
RX$~$J1713.7-3946 (also known as G347.3-0.5) has the largest surface
brightness, allowing us in the past to perform the most comprehensive study of
morphology and spatially resolved spectra of any such very-high-energy
gamma-ray source. Here we present extensive new H.E.S.S. measurements of
RX$~$J1713.7-3946, almost doubling the observation time compared to our
previous publication. Combined with new improved analysis tools, the previous
sensitivity is more than doubled. The H.E.S.S. angular resolution of
$0.048^\circ$ ($0.036^\circ$ above 2 TeV) is unprecedented in gamma-ray
astronomy and probes physical scales of 0.8 (0.6) parsec at the remnant's
location.
The new H.E.S.S. image of RX$~$J1713.7-3946 allows us to reveal clear
morphological differences between X-rays and gamma rays. In particular, for the
outer edge of the brightest shell region, we find the first ever indication for
particles in the process of leaving the acceleration shock region. By studying
the broadband energy spectrum, we furthermore extract properties of the parent
particle populations, providing new input to the discussion of the leptonic or
hadronic nature of the gamma-ray emission mechanism.
-
The propagation of ionic perturbations in a dusty plasma is considered
through a three-species kinetic simulation approach, in which the temporal
evolution of all three elements i.e. electrons, ions and dust particles are
followed based on the Vlasov equation coupled with the Poisson equation. Two
cases are focused upon: firstly a fully electron depleted dusty plasma, i.e., a
plasma consisting of ions and dust-particles. The second case includes dusty
plasmas with large electron-to-ion temperature ratios. The main features of the
ionic waves in these two settings including the dispersion relation and the
Landau damping rate are studied. It is shown that the dispersion relation of
the ionic waves perfectly matches the dispersion relation of Langmuir waves and
hence are called Langmuir-like ionic waves and can be considered as
ion-Langmuir waves. These waves can be theoretically predicted by the
dispersion relation of the dust-ion-acoustic waves. The transition of ionic
waves from dust-ion-acoustic to Langmuir-like waves are shown to be
sharp/smooth in first/second case. The Landau damping rates based on simulation
results are presented and compared with theoretical predictions wherever
possible.
-
The addition of a 28 m Cherenkov telescope (CT5) to the H.E.S.S. array
extended the experiment's sensitivity to lower energies. The lowest energy
threshold is obtained using monoscopic analysis of data taken with CT5,
providing access to gamma-ray energies below 100 GeV. Such an extension of the
instrument's energy range is particularly beneficial for studies of Active
Galactic Nuclei with soft spectra, as expected for those at a redshift > 0.5.
The high-frequency peaked BL Lac objects PKS 2155-304 (z = 0.116) and PG
1553+113 (0.43 < z < 0.58) are among the brightest objects in the gamma-ray
sky, both showing clear signatures of gamma-ray absorption at E > 100 GeV
interpreted as being due to interactions with the extragalactic background
light (EBL). Multiple observational campaigns of PKS 2155-304 and PG 1553+113
were conducted during 2013 and 2014 using the full H.E.S.S. II instrument. A
monoscopic analysis of the data taken with the new CT5 telescope was developed
along with an investigation into the systematic uncertainties on the spectral
parameters. The energy spectra of PKS 2155-304 and PG 1553+113 were
reconstructed down to energies of 80 GeV for PKS 2155-304, which transits near
zenith, and 110 GeV for the more northern PG 1553+113. The measured spectra,
well fitted in both cases by a log-parabola spectral model (with a 5.0 sigma
statistical preference for non-zero curvature for PKS 2155-304 and 4.5 sigma
for PG 1553+113), were found consistent with spectra derived from
contemporaneous Fermi-LAT data, indicating a sharp break in the observed
spectra of both sources at E ~ 100 GeV. When corrected for EBL absorption, the
intrinsic H.E.S.S. II mono and Fermi-LAT spectrum of PKS 2155-304 was found to
show significant curvature. For PG 1553+113, however, no significant detection
of curvature in the intrinsic spectrum could be found within statistical and
systematic uncertainties.
-
Aims: Following the detection of the fast radio burst FRB150418 by the SUPERB
project at the Parkes radio telescope, we aim to search for very-high energy
gamma-ray afterglow emission. Methods: Follow-up observations in the very-high
energy gamma-ray domain were obtained with the H.E.S.S. imaging atmospheric
Cherenkov telescope system within 14.5 hours of the radio burst. Results: The
obtained 1.4 hours of gamma-ray observations are presented and discussed. At
the 99 % C.L. we obtained an integral upper limit on the gamma-ray flux of
(E>350 GeV) < 1.33 x 10^-8 m^-2s^-1. Differential flux upper limits as function
of the photon energy were derived and used to constrain the intrinsic
high-energy afterglow emission of FRB 150418. Conclusions: No hints for
high-energy afterglow emission of FRB 150418 were found. Taking absorption on
the extragalactic background light into account and assuming a distance of z =
0.492 based on radio and optical counterpart studies and consistent with the
FRB dispersion, we constrain the gamma-ray luminosity at 1 TeV to L < 5.1 x
10^47 erg/s at 99% C.L.
-
List of contributions from the Cherenkov Telescope Array (CTA) Consortium
presented at the 6th International Symposium on High-Energy Gamma-Ray Astronomy
(Gamma 2016), July 11-15, 2016, in Heidelberg, Germany.
-
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_{\mu\mu}$ playing a
major role. Since the diffusion coefficient $D_{\mu\mu}$ 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.
-
A fully kinetic simulation approach, treating each plasma component based on
the Vlasov equation, is adopted to study the disintegration of an initial
density perturbation (IDP) into a number of ion-acoustic solitary waves (IASWs)
in the presence of the trapping effect of electrons. The non-linear fluid
theory developed by Schamel has identified three separate regimes of
ion-acoustic solitary waves based on the trapping parameter. Here, the
disintegration process and the resulting self-consistent IASWs are studied in a
wide range of trapping parameters covering all the three regimes continuously.
The dependency of features such as the time of disintegration, the number,
speed and size of IASWs on the trapping parameter are focused upon. It is shown
that an increase in this parameter slows down the propagation of IASWs while
decreases their sizes in the phase space. These features of IASWs tend to
saturate for large value of trapping parameters. The disintegration time shows
a more complicated behavior than what was predicted by the theoretical
approach. Also for the case of trapping parameters bigger than one, propagation
of IASWs is observed in contrast with the theoretical predictions. The kinetic
simulation results unveil a smooth and well-defined dependency of solitary
waves' features on the trapping parameter, showing the possibility of bridging
all the three regimes. Finally, it is shown that for beta around zero, the
electron phase space structure of the accompanying vortex stays symmetric. The
effect of the electron-to-ion temperature ratio on the disintegration and the
propagation of IASWs are considered as a benchmarking test of the simulation
code (in the nonlinear regime).
-
A search for dark matter line-like signals was performed in the vicinity of
the Galactic Centre by the H.E.S.S. experiment on observational data taken in
2014. An unbinned likelihood analysis was developed to improve the sensitivity
to line-like signals. The upgraded analysis along with newer data extend the
energy coverage of the previous measurement down to 100 GeV. The 18 h of data
collected with the H.E.S.S. array allow one to rule out at 95% CL the presence
of a 130 GeV line (at $l = -1.5^{\circ}, b = 0^{\circ}$ and for a dark matter
profile centered at this location) previously reported in Fermi-LAT data. This
new analysis overlaps significantly in energy with previous Fermi-LAT and
H.E.S.S. results. No significant excess associated with dark matter
annihilations was found in the energy range 100 GeV to 2 TeV and upper limits
on the gamma-ray flux and the velocity weighted annihilation cross-section are
derived adopting an Einasto dark matter halo profile. Expected limits for
present and future large statistics H.E.S.S. observations are also given.
-
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. The detection of a $\gamma$-ray
source coincident with W49B at very high energies (VHE; E > 100 GeV) with the
H.E.S.S. Cherenkov telescopes is reported together with a study of the source
with 5 years of Fermi-LAT high energy $\gamma$-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\pm20$ MeV and
$8.4_{-2.5}^{+2.2}$ GeV, the latter being 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
$\gamma$-ray emission produced through neutral-pion decay.
A search for new supernova remnants (SNRs) has been conducted using TeV gamma-ray data from the H.E.S.S. Galactic plane survey. As an identification criterion, shell morphologies that are characteristic for known resolved TeV SNRs have been used. Three new SNR candidates were identified in the H.E.S.S. data set with this method. Extensive multiwavelength searches for counterparts were conducted. A radio SNR candidate has been identified to be a counterpart to HESS J1534-571. The TeV source is therefore classified as a SNR. For the other two sources, HESS J1614-518 and HESS J1912+101, no identifying counterparts have been found, thus they remain SNR candidates for the time being. TeV-emitting SNRs are key objects in the context of identifying the accelerators of Galactic cosmic rays. The TeV emission of the relativistic particles in the new sources is examined in view of possible leptonic and hadronic emission scenarios, taking the current multiwavelength knowledge into account.
Aims. We study gamma-ray emission from the shell-type supernova remnant (SNR) RX J0852.0$-$4622 to better characterize its spectral properties and its distribution over the SNR. Methods. The analysis of an extended High Energy Spectroscopic System (H.E.S.S.) data set at very high energies (E > 100 GeV) permits detailed studies, as well as spatially resolved spectroscopy, of the morphology and spectrum of the whole RX J0852.0$-$4622 region. The H.E.S.S. data are combined with archival data from other wavebands and interpreted in the framework of leptonic and hadronic models. The joint Fermi-LAT-H.E.S.S. spectrum allows the direct determination of the spectral characteristics of the parent particle population in leptonic and hadronic scenarios using only GeV-TeV data. Results. An updated analysis of the H.E.S.S. data shows that the spectrum of the entire SNR connects smoothly to the high-energy spectrum measured by Fermi-LAT. The increased data set makes it possible to demonstrate that the H.E.S.S. spectrum deviates significantly from a power law and is well described by both a curved power law and a power law with an exponential cutoff at an energy of $E_\mathrm{cut} = (6.7 \pm 1.2_\mathrm{stat} \pm 1.2_\mathrm{syst})$ TeV. The joint Fermi-LAT-H.E.S.S. spectrum allows the unambiguous identification of the spectral shape as a power law with an exponential cutoff. No significant evidence is found for a variation of the spectral parameters across the SNR, suggesting similar conditions of particle acceleration across the remnant. A simple modeling using one particle population to model the SNR emission demonstrates that both leptonic and hadronic emission scenarios remain plausible. It is also shown that at least a part of the shell emission is likely due to the presence of a pulsar wind nebula around PSR J0855$-$4644.
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.
Understanding the physics of beam-plasma systems is of great importance for a large number of astrophysical systems. Recently the validity of particle-in-cell simulations of these systems has been questioned by Shalaby et al. (2017) We analyse one specific beam-plasma and its numerical simulations Kempf et al. (2016) that has been critized and prove that in fact the simulations performed fulfill the new validity criteria introduced by Shalaby et al. (2016).
Very high-energy gamma-rays (VHE, E>100 GeV) propagating over cosmological distances can interact with the low-energy photons of the extragalactic background light (EBL) and produce electron-positron pairs. The transparency of the universe to VHE gamma-rays is then directly related to the spectral energy distribution (SED) of the EBL. The observation of features in the VHE energy spectra of extragalactic sources allows the EBL to be measured, which otherwise is very difficult to determine. An EBL-model independent measurement of the EBL SED with the H.E.S.S. array of Cherenkov telescopes is presented. It is obtained by extracting the EBL absorption signal from the reanalysis of high-quality spectra of blazars. From H.E.S.S. data alone the EBL signature is detected at a significance of 9.5 sigma, and the intensity of the EBL obtained in different spectral bands is presented together with the associated gamma-ray horizon.
The diffuse very high-energy (VHE, > 100 GeV) gamma-ray emission observed in the central 200 pc of the Milky Way by H.E.S.S. was found to follow the dense matter distribution in the Central Molecular Zone (CMZ) up to a longitudinal distance of about 130 pc to the Galactic Centre (GC), where the flux rapidly decreases. Recent sensitive H.E.S.S. analysis revealed that the cosmic-ray (CR) density profile drops with the distance to the centre, making data compatible with a steady cosmic PeVatron at the GC. In this paper, we extend this analysis to obtain for the first time a detailed characterisation of the correlation with matter and to search for additional features and individual gamma-ray sources in the inner 200 pc. Taking advantage of 250 hours of H.E.S.S. data and improved analysis techniques we perform a detailed morphology study of the diffuse VHE emission observed from the GC ridge and reconstruct its total spectrum. We show that the emission correlated with dense matter covers the full CMZ and that its flux is about half the total diffuse emission flux. We also detect some emission at higher latitude likely produced by hadronic collisions of CRs in less dense regions of the GC interstellar medium. We detect an additional emission component centred on the GC and extending over about 15 pc that is consistent with the existence of a strong CR density gradient and confirms the presence of a CR accelerator at the very centre of our Galaxy. We show that the spectrum of the full ridge diffuse emission is compatible with the one previously derived from the central regions, suggesting that a single population of particles fills the entire CMZ. Finally, we report the discovery of a VHE gamma-ray source near the GC radio arc and argue that it is produced by the pulsar wind nebula candidate G0.13-0.11.
Runaway stars form bow shocks by ploughing through the interstellar medium at supersonic speeds and are promising sources of non-thermal emission of photons. One of these objects has been found to emit non-thermal radiation in the radio band. This triggered the development of theoretical models predicting non-thermal photons from radio up to very-high-energy (VHE, E $\geq 0.1$TeV) gamma rays. Subsequently, one bow shock was also detected in X-ray observations. However, the data did not allow discrimination between a hot thermal and a non-thermal origin. Further observations of different candidates at X-ray energies showed no evidence for emission at the position of the bow shocks either. A systematic search in the \textit{Fermi}-LAT energy regime resulted in flux upper limits for 27 candidates listed in the E-BOSS catalogue. Here we perform the first systematic search for VHE gamma-ray emission from bow shocks of runaway stars. Using all available archival H.E.S.S. data we search for very-high-energy gamma-ray emission at the positions of bow shock candidates listed in the second E-BOSS catalogue release. Out of the 73 bow shock candidates in this catalogue, 32 have been observed with H.E.S.S. None of the observed 32 bow shock candidates in this population study show significant emission in the H.E.S.S. energy range. Therefore, flux upper limits are calculated in five energy bins and the fraction of the kinetic wind power that is converted into VHE gamma rays is constrained. Emission from stellar bow shocks is not detected in the energy range between 0.14 and 18 TeV. The resulting upper limits constrain the level of VHE gamma-ray emission from these objects down to 0.1-1 % of the kinetic wind energy.
The nonlinear fluid theory developed by Schamel suggests a modified KdV equation to describe the temporal evolution of ion acoustic (IA) solitons in the presence of trapped electrons. The validity of this theory is studied here by verifying solitons main characteristic, i.e., stability against successive mutual collisions. We have employed a kinetic model as a more comprehensive theory than the fluid one, and utilized a fully kinetic simulation approach (both ions and electrons are treated based on the Vlasov equation). In the simulation approach, these solitons are excited self-consistently by employing the nonlinear process of IA solitons formation from an initial density perturbation (IDP). The effect of the size of IDPs on the chain formation is proved by the simulation code as a benchmark test. It is shown that the IA solitons, in presence of trapped electrons, can retain their features (both in spatial and velocity direction) after successive mutual collisions. The collisions here include encounters of IA solitons with the same trapping parameter, while differing in size. Kinetic simulation results reveal a complicated behavior during a collision between IA solitons in contrast to the fluid theory predictions and simulations. In the range of parameters considered here two oppositely propagating solitons rotate around their collective center in the phase space during a collision, independent of their trapping parameters. Furthermore, they exchange some portions of their trapped populations.
The nine-year H.E.S.S. Galactic Plane Survey (HGPS) 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 $\dot{E}$. This seems to be caused both by an increase of extension with decreasing $\dot{E}$, and hence with time, compatible with a power law $R_\mathrm{PWN}(\dot{E}) \sim \dot{E}^{-0.65 \pm 0.20}$, and by a mild decrease of TeV gamma-ray luminosity with decreasing $\dot{E}$, compatible with $L_{1-10\,\mathrm{TeV}} \sim \dot{E}^{0.59 \pm 0.21}$. We also find that the offsets of pulsars with respect to the wind nebula centres 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 $L_{1-10\,\mathrm{TeV}}/\dot{E}$. In addition to 14 HGPS sources considered as firmly identified pulsar wind nebulae and 5 additional pulsar wind nebulae taken from literature, we find 10 HGPS sources that form 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 theory.
The overtaking collisions of ion-acoustic solitons (IASs) in presence of trapping effects of electrons are studied based on a fully kinetic simulation approach. The method is able to provide all the kinetic details of the process alongside the fluid-level quantities self consistently. Solitons are produced naturally by utilizing the chain formation phenomenon, then are arranged in a new simulation box to test different scenarios of overtaking collisions. Three achievements are reported here. Firstly, simulations prove the long-time life span of the ion-acoustic solitons in the presence of trapping effect of electrons (kinetic effects), which serves as the benchmark of the simulation code. Secondly, their stability against overtaking mutual collisions is established by creating collisions between solitons with different number and shapes of trapped electrons, i.e. different trapping parameter. Finally, details of solitons during collisions for both ions and electrons are provided on both fluid and kinetic levels. These results show that on the kinetic level, trapped electron population accompanying each of the solitons are exchanged between the solitons during the collision. Furthermore, the behavior of electron holes accompanying solitons contradicts the theory about the electron holes interaction developed based on kinetic theory. They also show behaviors much different from other electron holes witnessed in processes such as nonlinear Landau damping (Bernstein-Greene-Kruskal -BGK- modes) or beam-plasma interaction (like two-beam instability).
Supernova remnants exhibit shock fronts (shells) that can accelerate charged particles up to very high energies. In the past decade, measurements of a handful of shell-type supernova remnants in very-high-energy gamma rays have provided unique insights into the acceleration process. Among those objects, RX$~$J1713.7-3946 (also known as G347.3-0.5) has the largest surface brightness, allowing us in the past to perform the most comprehensive study of morphology and spatially resolved spectra of any such very-high-energy gamma-ray source. Here we present extensive new H.E.S.S. measurements of RX$~$J1713.7-3946, almost doubling the observation time compared to our previous publication. Combined with new improved analysis tools, the previous sensitivity is more than doubled. The H.E.S.S. angular resolution of $0.048^\circ$ ($0.036^\circ$ above 2 TeV) is unprecedented in gamma-ray astronomy and probes physical scales of 0.8 (0.6) parsec at the remnant's location. The new H.E.S.S. image of RX$~$J1713.7-3946 allows us to reveal clear morphological differences between X-rays and gamma rays. In particular, for the outer edge of the brightest shell region, we find the first ever indication for particles in the process of leaving the acceleration shock region. By studying the broadband energy spectrum, we furthermore extract properties of the parent particle populations, providing new input to the discussion of the leptonic or hadronic nature of the gamma-ray emission mechanism.
The propagation of ionic perturbations in a dusty plasma is considered through a three-species kinetic simulation approach, in which the temporal evolution of all three elements i.e. electrons, ions and dust particles are followed based on the Vlasov equation coupled with the Poisson equation. Two cases are focused upon: firstly a fully electron depleted dusty plasma, i.e., a plasma consisting of ions and dust-particles. The second case includes dusty plasmas with large electron-to-ion temperature ratios. The main features of the ionic waves in these two settings including the dispersion relation and the Landau damping rate are studied. It is shown that the dispersion relation of the ionic waves perfectly matches the dispersion relation of Langmuir waves and hence are called Langmuir-like ionic waves and can be considered as ion-Langmuir waves. These waves can be theoretically predicted by the dispersion relation of the dust-ion-acoustic waves. The transition of ionic waves from dust-ion-acoustic to Langmuir-like waves are shown to be sharp/smooth in first/second case. The Landau damping rates based on simulation results are presented and compared with theoretical predictions wherever possible.
The addition of a 28 m Cherenkov telescope (CT5) to the H.E.S.S. array extended the experiment's sensitivity to lower energies. The lowest energy threshold is obtained using monoscopic analysis of data taken with CT5, providing access to gamma-ray energies below 100 GeV. Such an extension of the instrument's energy range is particularly beneficial for studies of Active Galactic Nuclei with soft spectra, as expected for those at a redshift > 0.5. The high-frequency peaked BL Lac objects PKS 2155-304 (z = 0.116) and PG 1553+113 (0.43 < z < 0.58) are among the brightest objects in the gamma-ray sky, both showing clear signatures of gamma-ray absorption at E > 100 GeV interpreted as being due to interactions with the extragalactic background light (EBL). Multiple observational campaigns of PKS 2155-304 and PG 1553+113 were conducted during 2013 and 2014 using the full H.E.S.S. II instrument. A monoscopic analysis of the data taken with the new CT5 telescope was developed along with an investigation into the systematic uncertainties on the spectral parameters. The energy spectra of PKS 2155-304 and PG 1553+113 were reconstructed down to energies of 80 GeV for PKS 2155-304, which transits near zenith, and 110 GeV for the more northern PG 1553+113. The measured spectra, well fitted in both cases by a log-parabola spectral model (with a 5.0 sigma statistical preference for non-zero curvature for PKS 2155-304 and 4.5 sigma for PG 1553+113), were found consistent with spectra derived from contemporaneous Fermi-LAT data, indicating a sharp break in the observed spectra of both sources at E ~ 100 GeV. When corrected for EBL absorption, the intrinsic H.E.S.S. II mono and Fermi-LAT spectrum of PKS 2155-304 was found to show significant curvature. For PG 1553+113, however, no significant detection of curvature in the intrinsic spectrum could be found within statistical and systematic uncertainties.
Aims: Following the detection of the fast radio burst FRB150418 by the SUPERB project at the Parkes radio telescope, we aim to search for very-high energy gamma-ray afterglow emission. Methods: Follow-up observations in the very-high energy gamma-ray domain were obtained with the H.E.S.S. imaging atmospheric Cherenkov telescope system within 14.5 hours of the radio burst. Results: The obtained 1.4 hours of gamma-ray observations are presented and discussed. At the 99 % C.L. we obtained an integral upper limit on the gamma-ray flux of (E>350 GeV) < 1.33 x 10^-8 m^-2s^-1. Differential flux upper limits as function of the photon energy were derived and used to constrain the intrinsic high-energy afterglow emission of FRB 150418. Conclusions: No hints for high-energy afterglow emission of FRB 150418 were found. Taking absorption on the extragalactic background light into account and assuming a distance of z = 0.492 based on radio and optical counterpart studies and consistent with the FRB dispersion, we constrain the gamma-ray luminosity at 1 TeV to L < 5.1 x 10^47 erg/s at 99% C.L.
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_{\mu\mu}$ playing a major role. Since the diffusion coefficient $D_{\mu\mu}$ 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.
A fully kinetic simulation approach, treating each plasma component based on the Vlasov equation, is adopted to study the disintegration of an initial density perturbation (IDP) into a number of ion-acoustic solitary waves (IASWs) in the presence of the trapping effect of electrons. The non-linear fluid theory developed by Schamel has identified three separate regimes of ion-acoustic solitary waves based on the trapping parameter. Here, the disintegration process and the resulting self-consistent IASWs are studied in a wide range of trapping parameters covering all the three regimes continuously. The dependency of features such as the time of disintegration, the number, speed and size of IASWs on the trapping parameter are focused upon. It is shown that an increase in this parameter slows down the propagation of IASWs while decreases their sizes in the phase space. These features of IASWs tend to saturate for large value of trapping parameters. The disintegration time shows a more complicated behavior than what was predicted by the theoretical approach. Also for the case of trapping parameters bigger than one, propagation of IASWs is observed in contrast with the theoretical predictions. The kinetic simulation results unveil a smooth and well-defined dependency of solitary waves' features on the trapping parameter, showing the possibility of bridging all the three regimes. Finally, it is shown that for beta around zero, the electron phase space structure of the accompanying vortex stays symmetric. The effect of the electron-to-ion temperature ratio on the disintegration and the propagation of IASWs are considered as a benchmarking test of the simulation code (in the nonlinear regime).
A search for dark matter line-like signals was performed in the vicinity of the Galactic Centre by the H.E.S.S. experiment on observational data taken in 2014. An unbinned likelihood analysis was developed to improve the sensitivity to line-like signals. The upgraded analysis along with newer data extend the energy coverage of the previous measurement down to 100 GeV. The 18 h of data collected with the H.E.S.S. array allow one to rule out at 95% CL the presence of a 130 GeV line (at $l = -1.5^{\circ}, b = 0^{\circ}$ and for a dark matter profile centered at this location) previously reported in Fermi-LAT data. This new analysis overlaps significantly in energy with previous Fermi-LAT and H.E.S.S. results. No significant excess associated with dark matter annihilations was found in the energy range 100 GeV to 2 TeV and upper limits on the gamma-ray flux and the velocity weighted annihilation cross-section are derived adopting an Einasto dark matter halo profile. Expected limits for present and future large statistics H.E.S.S. observations are also given.