We present an analysis of X-ray spectra of the high mass X-ray binary 4U 0115+634 as observed with Suzaku and RXTE in 2011 July, during the fading phase of a giant X-ray outburst. We used a continuum model consisting of an absorbed cutoff power-law and an ad-hoc Gaussian emission feature centered around 8.5 keV, which we discuss to be due to cyclotron emission. Our results are consistent with a fundamental cyclotron absorption line centered at ${\sim}10.2$ keV for all observed flux ranges. At the same time we rule out significant influence of the 8.5 kev Gaussian on the CRSF parameters, which are not consistent with the cyclotron line energies and depths of previously reported flux-dependent descriptions. We also show that some continuum models can lead to artificial line-like residuals in the analyzed spectra, which are then misinterpreted as unphysically strong cyclotron lines. Specifically, our results do not support the existence of a previously claimed additional cyclotron feature at ${\sim}15$ keV. Apart from these features, we find for the first time evidence for a He-like Fe XXV emission line at ${\sim}6.7$ keV and weak H-like Fe XXVI emission close to ${\sim}7.0$ keV.

Cyclotron lines, also called cyclotron resonant scattering features (CRSF) are spectral features, generally appearing in absorption, in the X-ray spectra of objects containing highly magnetized neutron stars, allowing the direct measurement of the magnetic field strength in these objects. Cyclotron features are thought to be due to resonant scattering of photons by electrons in the strong magnetic fields. The main content of this contribution focusses on electron cyclotron lines as found in accreting X-ray binary pulsars (XRBP) with magnetic fields on the order of several 1012 Gauss. Also, possible proton cyclotron lines from single neutron stars with even stronger magnetic fields are briefly discussed.With regard to electron cyclotron lines, we present an updated list of XRBPs that show evidence of such absorption lines. The first such line was discovered in a 1976 balloon observation of the accreting binary pulsar Hercules X-1, it is considered to be the first direct measurement of the magnetic field of a neutron star. As of today (mid 2018), we list 36 XRBPs showing evidence of one ore more electron cyclotron absorption line(s). A few have been measured only once and must be confirmed (several more objects are listed as candidates). In addition to the Tables of objects, we summarize the evidence of variability of the cyclotron line as a function of various parameters (especially pulse phase, luminosity and time), and add a discussion of the different observed phenomena and associated attempts of theoretical modeling. We also discuss our understanding of the underlying physics of accretion onto highly magnetized neutron stars. For proton cyclotron lines, we present tables with seven neutron stars and discuss their nature and the physics in these objects.

Centaurus A (Cen A) is the nearest radio galaxy discovered as a very-high-energy (VHE; 100 GeV-100 TeV) $\gamma$-ray source by the High Energy Stereoscopic System (H.E.S.S.). It is a faint VHE $\gamma$-ray emitter, though its VHE flux exceeds both the extrapolation from early Fermi-LAT observations as well as expectations from a (misaligned) single-zone synchrotron-self Compton (SSC) description. The latter satisfactorily reproduces the emission from Cen A at lower energies up to a few GeV. New observations with H.E.S.S., comparable in exposure time to those previously reported, were performed and eight years of Fermi-LAT data were accumulated to clarify the spectral characteristics of the $\gamma$-ray emission from the core of Cen A. The results allow us for the first time to achieve the goal of constructing a representative, contemporaneous $\gamma$-ray core spectrum of Cen A over almost five orders of magnitude in energy. Advanced analysis methods, including the template fitting method, allow detection in the VHE range of the core with a statistical significance of 12$\sigma$ on the basis of 213 hours of total exposure time. The spectrum in the energy range of 250 GeV-6 TeV is compatible with a power-law function with a photon index $\Gamma=2.52\pm0.13_{\mathrm{stat}}\pm0.20_{\mathrm{sys}}$. An updated Fermi-LAT analysis provides evidence for spectral hardening by $\Delta\Gamma\simeq0.4\pm0.1$ at $\gamma$-ray energies above $2.8^{+1.0}_{-0.6}$ GeV at a level of $4.0\sigma$. The fact that the spectrum hardens at GeV energies and extends into the VHE regime disfavour a single-zone SSC interpretation for the overall spectral energy distribution (SED) of the core and is suggestive of a new $\gamma$-ray emitting component connecting the high-energy emission above the break energy to the one observed at VHE energies.

G15.9+0.2 is a Galactic shell-type supernova remnant (SNR), which was detected in radio and has been confirmed in X-rays based on Chandra observations. An X-ray point source CXOUJ181852.0-150213 has been detected and suggested to be an associated neutron star. In a recent study, we have confirmed the source to be a central compact object (CCO). We have studied the SNR using high-resolution X-ray data taken with Chandra in combination with infrared (IR) data in order to understand its emission and to derive its physical parameters. This will also help to constrain, e.g., the age of the CCO and the environment in which it was born. The spectral analysis of the X-ray emission using the new Chandra data and the comparison to the IR data have shown that the SNR is relatively young with an age of a few thousand years and that its emission is dominated by that of shocked interstellar medium (ISM). However, the analysis of the spectrum of the bright eastern shell shows that there is an additional emission component with enhanced abundances of {\alpha} elements and Fe, suggesting ejecta emission. The multi-wavelength emission is consistent with SNR G15.9+0.2 expanding in an ISM with a density gradient, while there is also colder material located in front of the SNR, which absorbs its thermal X-ray emission in the softer bands.

Spectral lines are among the most powerful signatures for dark matter (DM) annihilation searches in very-high-energy $\gamma$-rays. The central region of the Milky Way halo is one of the most promising targets given its large amount of DM and proximity to Earth. We report on a search for a monoenergetic spectral line from self-annihilations of DM particles in the energy range from 300 GeV to 70 TeV using a two-dimensional maximum likelihood method taking advantage of both the spectral and spatial features of signal versus background. The analysis makes use of Galactic Center (GC) observations accumulated over ten years (2004 - 2014) with the H.E.S.S. array of ground-based Cherenkov telescopes. No significant $\gamma$-ray excess above the background is found. We derive upper limits on the annihilation cross section $\langle\sigma v\rangle$ for monoenergetic DM lines at the level of $\sim4\times10^{-28}$ cm$^{3}$s$^{-1}$ at 1 TeV, assuming an Einasto DM profile for the Milky Way halo. For a DM mass of 1 TeV, they improve over the previous ones by a factor of six. The present constraints are the strongest obtained so far for DM particles in the mass range 300 GeV - 70 TeV. Ground-based $\gamma$-ray observations have reached sufficient sensitivity to explore relevant velocity-averaged cross sections for DM annihilation into two $\gamma$-ray photons at the level expected from the thermal relic density for TeV DM particles.

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.

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 \pm 0.50_{stat} \pm 0.10_{syst}$ and normalization (at $E_0$=1.0 TeV) $N_0(E_0)=(0.58 \pm 0.22_{stat} \pm 0.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.

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.

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.

Recent observations of Her X-1 with NuSTAR and INTEGRAL in 2016 have provided evidence that the 20-year decay of the cyclotron line energy found between 1996 and 2015 has ended and that an inversion with a new increase, possibly similar to the one observed around 1990-1993, has started. We consider this a strong motivation for further observations and for enhanced efforts to significantly improve our theoretical understanding of the accretion process in binary X-ray pulsars. We speculate about the physics behind the long-term decay and its inversion, a possible cyclic behavior, and correlations with other variable observables.

Index of H.E.S.S. conference proceedings to the 35th ICRC, Busan, Korea

We present a comprehensive spectral analysis of the BeXRB GRO J1008-57 over a luminosity range of three orders of magnitude using NuSTAR, Suzaku and RXTE data. We find significant evolution of the spectral parameters with luminosity. In particular the photon index hardens with increasing luminosity at intermediate luminosities between $10^{36}$ $-$ $10^{37}$ erg s$^{-1}$. This evolution is stable and repeatedly observed over different outbursts. However, at the extreme ends of the observed luminosity range, we find that the correlation breaks down, with a significance level of at least $3.7\sigma$. We conclude that these changes indicate transitions to different accretion regimes, which are characterized by different deceleration processes, such as Coulomb or radiation breaking. We compare our observed luminosity levels of these transitions to theoretical predications and discuss the variation of those theoretical luminosity values with fundamental neutron star parameters. Finally, we present detailed spectroscopy of the unique "triple peaked" outburst in 2014/15 which does not fit in the general parameter evolution with luminosity. The pulse profile on the other hand is consistent with what is expected at this luminosity level, arguing against a change in accretion geometry. In summary, GRO J1008-57 is an ideal target to study different accretion regimes due to the well constrained evolution of its broad-band spectral continuum over several orders of magnitude in luminosity.

We report on the analysis of XMM-Newton observations of the non-thermal shell-type supernova remnant HESS J1731-347 (G353.6-0.7). For the first time the complete remnant shell has been covered in X-rays, which allowed direct comparison with radio and TeV observations. We carried out a spatially resolved spectral analysis of XMM-Newton data and confirmed the previously reported non-thermal power-law X-ray spectrum of the source with negligible variations of spectral index across the shell. On the other hand, the X-ray absorption column is strongly variable and correlates with the CO emission thus confirming that the absorbing material must be in the foreground and reinforcing the previously suggested lower limit on distance. Finally, we find that the X-ray emission of the remnant is suppressed towards the Galactic plane, which points to lower shock velocities in this region, likely due to the interaction of the shock with the nearby molecular cloud.

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 large jet kinetic power and non-thermal processes occurring in the microquasar SS 433 make this source a good candidate for a very high-energy (VHE) gamma-ray emitter. Gamma-ray fluxes have been predicted for both the central binary and the interaction regions between jets and surrounding nebula. Also, non-thermal emission at lower energies has been previously reported. We explore the capability of SS 433 to emit VHE gamma rays during periods in which the expected flux attenuation due to periodic eclipses and precession of the circumstellar disk periodically covering the central binary system is expected to be at its minimum. The eastern and western SS433/W50 interaction regions are also examined. We aim to constrain some theoretical models previously developed for this system. We made use of dedicated observations from MAGIC and H.E.S.S. from 2006 to 2011 which were combined for the first time and accounted for a total effective observation time of 16.5 h. Gamma-ray attenuation does not affect the jet/medium interaction regions. The analysis of a larger data set amounting to 40-80 h, depending on the region, was employed. No evidence of VHE gamma-ray emission was found. Upper limits were computed for the combined data set. We place constraints on the particle acceleration fraction at the inner jet regions and on the physics of the jet/medium interactions. Our findings suggest that the fraction of the jet kinetic power transferred to relativistic protons must be relatively small to explain the lack of TeV and neutrino emission from the central system. At the SS433/W50 interface, the presence of magnetic fields greater 10$\mu$G is derived assuming a synchrotron origin for the observed X-ray emission. This also implies the presence of high-energy electrons with energies up to 50 TeV, preventing an efficient production of gamma-ray fluxes in these interaction regions.

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.

Be/X-ray binaries are the most populous class of High Mass X-ray Binaries. Their X-ray duty cycle is tightly related to the optical companion wind activity, which in turn can be studied through optical spectroscopical dedicated observations. We study optical spectral features of the Be circumstellar disk to test their long-term variability and their relation with the X-ray activity. Special attention has been given to the H$\alpha$ emission line, that is one of the best tracers of the disk conditions. We obtained optical broadband medium resolution spectra from a dedicated campaign with the Anglo-Australian Telescope and the Southern African Large Telescope in 2014-2015. Data span over one entire binary orbit, and cover both X-ray quiescent and moderately active periods. We used Balmer emission lines to follow the evolution of the circumstellar disk. We observe prominent spectral features, like double-peaked H$\alpha$ and H$\beta$ emission lines. The H$\alpha$ V/R ratio significantly changes over a time scale of about one year. Our observations are consistent with a system observed at a large inclination angle ($i\gtrsim60^{\circ}$). The derived circumstellar disk size shows that the disk evolves from a configuration that prevents accretion onto the neutron star, to one that allows only moderate accretion. This is in agreement with the contemporary observed X-ray activity. Our results are interpreted within the context of inefficient tidal truncation of the circumstellar disk, as expected for this source's binary configuration. We derived the H$\beta$-emitting region size, which results about half of the corresponding H$\alpha$-emitting disk, and constrain the luminosity class of V850 Cen as III-V, consistent with the previously proposed class.}

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.

X-ray spectra of accreting pulsars are generally observed to vary with their X-ray luminosity. In particular, the hardness of the X-ray continuum is found to depend on luminosity. In a few sources, the correlation between the energy of the cyclotron resonance scattering feature (CRSF) and the luminosity is clear. Different types (signs) of the correlation are believed to reflect different accretion modes. We analyse two NuSTAR observations of the transient accreting pulsar Cep X-4 during its 2014 outburst. Our analysis is focused on a detailed investigation of the dependence of the CRSF energy and of the spectral hardness on X-ray luminosity, especially on short timescales. To investigate the spectral changes as a function of luminosity within each of the two observations, we used the intrinsic variability of the source on the timescale of individual pulse cycles (tens of seconds), the so-called pulse-to-pulse variability. We find that the NuSTAR spectrum of Cep X-4 contains two CRSFs: the fundamental line at ~30 keV and its harmonic at ~55 keV. We find for the first time that the energy of the fundamental CRSF increases and the continuum becomes harder with increasing X-ray luminosity not only between the two observations, that is, on the long timescale, but also within an individual observation, on the timescale of a few tens of seconds. We investigate these dependencies in detail including their non-linearity. We discuss a possible physical interpretation of the observed behaviour in the frame of a simple one-dimensional model of the polar emitting region with a collisionless shock formed in the infalling plasma near the neutron star surface. With this model, we are able to reproduce the observed variations of the continuum hardness ratio and of the CRSF energy with luminosity.

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.

Cyclotron resonant scattering features (CRSFs) are formed by scattering of X-ray photons off quantized plasma electrons in the strong magnetic field (of the order 10^12 G) close to the surface of an accreting X-ray pulsar. The line profiles of CRSFs cannot be described by an analytic expression. Numerical methods such as Monte Carlo (MC) simulations of the scattering processes are required in order to predict precise line shapes for a given physical setup, which can be compared to observations to gain information about the underlying physics in these systems. A versatile simulation code is needed for the generation of synthetic cyclotron lines. Sophisticated geometries should be investigatable by making their simulation possible for the first time. The simulation utilizes the mean free path tables described in the first paper of this series for the fast interpolation of propagation lengths. The code is parallelized to make the very time consuming simulations possible on convenient time scales. Furthermore, it can generate responses to mono-energetic photon injections, producing Green's functions, which can be used later to generate spectra for arbitrary continua. We develop a new simulation code to generate synthetic cyclotron lines for complex scenarios, allowing for unprecedented physical interpretation of the observed data. An associated XSPEC model implementation is used to fit synthetic line profiles to NuSTAR data of Cep X-4. The code has been developed with the main goal of overcoming previous geometrical constraints in MC simulations of CRSFs. By applying this code also to more simple, classic geometries used in previous works, we furthermore address issues of code verification and cross-comparison of various models. The XSPEC model and the Green's function tables are available online at http://www.sternwarte.uni-erlangen.de/research/cyclo .

The X-ray spectra of the neutron stars located in the centers of supernova remnants Cas A and HESS J1731-347 are well fit with carbon atmosphere models. These fits yield plausible neutron star sizes for the known or estimated distances to these supernova remnants. The evidence in favor of the presence of a pure carbon envelope at the neutron star surface is rather indirect and is based on the assumption that the emission is generated uniformly by the entire stellar surface. Although this assumption is supported by the absence of pulsations, the observational upper limit on the pulsed fraction is not very stringent. In an attempt to quantify this evidence, we investigate the possibility that the observed spectrum of the neutron star in HESS J1731-347 is a combination of the spectra produced in a hydrogen atmosphere of the hotspots and of the cooler remaining part of the neutron star surface. The lack of pulsations in this case has to be explained either by a sufficiently small angle between the neutron star spin axis and the line of sight, or by a sufficiently small angular distance between the hotspots and the neutron star rotation poles. As the observed flux from a non-uniformly emitting neutron star depends on the angular distribution of the radiation emerging from the atmosphere, we have computed two new grids of pure carbon and pure hydrogen atmosphere model spectra accounting for Compton scattering. Using new hydrogen models, we have evaluated the probability of a geometry that leads to a pulsed fraction below the observed upper limit to be about 8.2 %. Such a geometry thus seems to be rather improbable but cannot be excluded at this stage.

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.

HESS J1731-347 is a shell-type supernova remnant emitting both TeV gamma rays and non-thermal X-ray photons, spatially coincident with the radio SNR G353.6-0.7. Hadronic and leptonic scenarios (or a blend of both) are discussed in the literature to explain the TeV emission from the object. In 2011, a $\gamma$-ray excess was also found in the neighborhood of the source (HESS J1729-345). Here we present results of an updated analysis obtained with the meanwhile available additional H.E.S.S. data. Beyond HESS J1731-347, the analysis reveals the morphology of the emission of the adjacent TeV source HESS J1729-345 and the emission in between the two sources in greater detail. The results permit us to correlate the TeV emission outside of the SNR with molecular gas tracers, and to confront the data with scenarios in which the TeV emission outside the SNR is produced by escaping cosmic rays.

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.

The RXTE observed four outbursts of the accreting X-ray binary transient source, GX 304-1 in 2010 and 2011. We present results of detailed 3-100 keV spectral analysis of 69 separate observations, and report a positive correlation between cyclotron line parameters, as well as other spectral parameters, with power law flux. The cyclotron line energy, width and depth versus flux, and thus luminosity, correlations show a flattening of the relationships with increasing luminosity, which are well described by quasi-spherical or disk accretion that yield the surface magnetic field to be ~60 keV. Since HEXTE cluster A was fixed aligned with the PCA field of view and cluster B was fixed viewing a background region 1.5 degrees off of the source direction during these observations near the end of the RXTE mission, the cluster A background was estimated from cluster B events using HEXTEBACKEST. This made possible the detection of the ~55 keV cyclotron line and an accurate measurement of the continuum. Correlations of all spectral parameters with the primary 2-10 keV power law flux reveal it to be the primary driver of the spectral shape. The accretion is found to be in the collisionless shock braking regime.

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 < E < 300 GeV) and very high energy (VHE, E > 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.

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.

Electron cyclotron resonant scattering features (CRSFs) are observed as absorption-like lines in the spectra of X-ray pulsars. A significant fraction of the computing time for Monte Carlo simulations of these quantum mechanical features is spent on the calculation of the mean free path for each individual photon before scattering, since it involves a complex numerical integration over the scattering cross section and the (thermal) velocity distribution of the scattering electrons. We aim to numerically calculate interpolation tables which can be used in CRSF simulations to sample the mean free path of the scattering photon and the momentum of the scattering electron. The tables also contain all the information required for sampling the scattering electron's final spin. The tables were calculated using an adaptive Simpson integration scheme. The energy and angle grids were refined until a prescribed accuracy is reached. The tables are used by our simulation code to produce artificial CRSF spectra. The electron momenta sampled during these simulations were analyzed and justified using theoretically determined boundaries. We present a complete set of tables suited for mean free path calculations of Monte Carlo simulations of the cyclotron scattering process for conditions expected in typical X-ray pulsar accretion columns (0.01<B/B_{crit}<=0.12, where B_{crit}=4.413x10^{13} G and 3keV<=kT<15keV). The sampling of the tables is chosen such that the results have an estimated relative error of at most 1/15 for all points in the grid. The tables are available online at http://www.sternwarte.uni-erlangen.de/research/cyclo.

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.

The inner region of the Milky Way halo harbors a large amount of dark matter (DM). Given its proximity, it is one of the most promising targets to look for DM. We report on a search for the annihilations of DM particles using $\gamma$-ray observations towards the inner 300 parsecs of the Milky Way, with the H.E.S.S. array of ground-based Cherenkov telescopes. The analysis is based on a 2D maximum likelihood method using Galactic center (GC) data accumulated by H.E.S.S. over the last 10 years (2004-2014), and does not show any significant $\gamma$-ray signal above background. Assuming Einasto and Navarro-Frenk-White DM density profiles at the GC, we derive upper limits on the annihilation cross section $\langle \sigma v\rangle$. These constraints are the strongest obtained so far in the TeV DM mass range and improve upon previous limits by a factor 5. For the Einasto profile, the constraints reach $\langle \sigma v\rangle$ values of $\rm 6\times10^{-26} cm^3s^{-1}$ in the $W^+W^-$ channel for a DM particle mass of 1.5 TeV, and $\rm 2\times10^{-26} cm^3s^{-1}$ in the $\tau^+\tau^-$ channel for 1 TeV mass. For the first time, ground-based $\gamma$-ray observations have reached sufficient sensitivity to probe $\langle \sigma v\rangle$ values expected from the thermal relic density for TeV DM particles.

We present our study of the central point source CXOU J181852.0-150213 in the young Galactic supernova remnant (SNR) G15.9+0.2 based on the recent ~90 ks Chandra observations. The point source was discovered in 2005 in shorter Chandra observations and was hypothesized to be a neutron star associated with the SNR. Our X-ray spectral analysis strongly supports the hypothesis of a thermally emitting neutron star associated with G15.9+0.2. We conclude that the object belongs to the class of young cooling low-magnetized neutron stars referred to as central compact objects (CCOs). We modeled the spectrum of the neutron star with a blackbody spectral function and with our hydrogen and carbon neutron star atmosphere models, assuming that the radiation is uniformly emitted by the entire stellar surface. Under this assumption, only the carbon atmosphere models yield a distance that is compatible with a source located in the Galaxy. In this respect, CXOU J181852.0-150213 is similar to two other well-studied CCOs, the neutron stars in Cas A and in HESS J1731-347, for which carbon atmosphere models were used to reconcile their emission with the known or estimated distances.

Microquasars are potential $\gamma$-ray emitters. Indications of transient episodes of $\gamma$-ray emission were recently reported in at least two systems: Cyg X-1 and Cyg X-3. The identification of additional $\gamma$-ray-emitting microquasars is required to better understand how $\gamma$-ray emission can be produced in these systems. Theoretical models have predicted very high-energy (VHE) $\gamma$-ray emission from microquasars during periods of transient outburst. Observations reported herein were undertaken with the objective of observing a broadband flaring event in the $\gamma$-ray and X-ray bands. Contemporaneous observations of three microquasars, GRS 1915+105, Circinus X-1, and V4641 Sgr, were obtained using the High Energy Spectroscopic System (H.E.S.S.) telescope array and the Rossi X-ray Timing Explorer (RXTE) satellite. X-ray analyses for each microquasar were performed and VHE $\gamma$-ray upper limits from contemporaneous H.E.S.S. observations were derived. No significant $\gamma$-ray signal has been detected in any of the three systems. The integral $\gamma$-ray photon flux at the observational epochs is constrained to be $I(> 560 {\rm\ GeV}) < 7.3\times10^{-13}$ cm$^{-2}$ s$^{-1}$, $I(> 560 {\rm\ GeV}) < 1.2\times10^{-12}$ cm$^{-2}$ s$^{-1}$, and $I(> 240 {\rm\ GeV}) < 4.5\times10^{-12}$ cm$^{-2}$ s$^{-1}$ for GRS 1915+105, Circinus X-1, and V4641 Sgr, respectively. The $\gamma$-ray upper limits obtained using H.E.S.S. are examined in the context of previous Cherenkov telescope observations of microquasars. The effect of intrinsic absorption is modelled for each target and found to have negligible impact on the flux of escaping $\gamma$-rays. When combined with the X-ray behaviour observed using RXTE, the derived results indicate that if detectable VHE $\gamma$-ray emission from microquasars is commonplace, then it is likely to be highly transient.

We took advantage of the large collecting area and good timing capabilities of the EPIC cameras on-board XMM-Newton to investigate the accretion geometry onto the magnetized neutron star hosted in the high mass X-ray binary EXO 2030+375 during the rise of a source Type-I outburst in 2014. We carried out a timing and spectral analysis of the XMM-Newton observation as function of the neutron star spin phase. We used a phenomenological spectral continuum model comprising the required fluorescence emission lines. Two neutral absorption components are present: one covering fully the source and one only partially. The same analysis was also carried out on two Suzaku observations of the source performed during outbursts in 2007 and 2012, to search for possible spectral variations at different luminosities. The XMM-Newton data caught the source at an X-ray luminosity of $2\times10^{36}$ erg s$^{-1}$ and revealed the presence of a narrow dip-like feature in its pulse profile that was never reported before. The width of this feature corresponds to about one hundredth of the neutron star spin period. From the results of the phase-resolved spectral analysis we suggest that this feature can be ascribed to the self-obscuration of the accretion stream passing in front of the observer line of sight. We inferred from the Suzaku observation carried out in 2007 that the self-obscuration of the accretion stream might produce a significantly wider feature in the neutron star pulsed profile at higher luminosities ($\gtrsim$$2\times10^{37}$ erg s$^{-1}$).

Using the High Energy Spectroscopic System (H.E.S.S.) telescopes we have discovered a steady and extended very high-energy (VHE) gamma-ray source towards the luminous blue variable candidate LBV1806-20, massive stellar cluster Cl*1806-20, and magnetar SGR1806-20. The new VHE source, HESSJ1808-204, was detected at a statistical significance of >6sigma (post-trial) with a photon flux normalisation $(2.9 \pm 0.4_{\rm stat} \pm 0.5_{\rm sys})\times 10^{-13}$\,ph\,cm$^{-2}$\,s$^{-1}$\,TeV$^{-1}$ at 1TeV and a power-law photon index of $2.3\pm0.2_{\rm stat}\pm 0.3_{\rm sys}$. The luminosity of this source (0.2 to 10 TeV; scaled to distance d=8.7kpc) is $L_{\rm VHE}\sim1.6 \times 10^{34}(d/{8.7kpc})^2$\,erg\,s$^{-1}$. The VHE gamma-ray emission is extended and is well fit by a single Gaussian with statistical standard deviation of $0.095^\circ \pm 0.015^\circ$. This extension is similar to that of the synchrotron radio nebula G10.0-0.3, which is thought to be powered by LBV1806-20. The VHE gamma-ray luminosity could be provided by the stellar wind luminosity of LBV1806-20 by itself and/or the massive star members of Cl*1806-20. Alternatively, magnetic dissipation (e.g. via reconnection) from SGR1806-20 can potentially account for the VHE luminosity. The origin and hadronic and/or leptonic nature of the accelerated particles responsible for HESSJ1808-204 is not yet clear. If associated with SGR1806-20, the potentially young age of the magnetar (650yr) can be used to infer the transport limits of these particles to match the VHE source size. This discovery provides new interest in the potential for high-energy particle acceleration from magnetars, massive stars, and/or stellar clusters.

Parallel computing and steadily increasing computation speed have led to a new tool for analyzing multiple datasets and datatypes: fitting several datasets simultaneously. With this technique, physically connected parameters of individual data can be treated as a single parameter by implementing this connection into the fit directly. We discuss the terminology, implementation, and possible issues of simultaneous fits based on the X-ray data analysis tool Interactive Spectral Interpretation System (ISIS). While all data modeling tools in X-ray astronomy allow in principle fitting data from multiple data sets individually, the syntax used in these tools is not often well suited for this task. Applying simultaneous fits to the transient X-ray binary GRO J1008-57, we find that the spectral shape is only dependent on X-ray flux. We determine time independent parameters such as, e.g., the folding energy E_fold, with unprecedented precision.

In a previous work, we introduced a tool for analyzing multiple datasets simultaneously, which has been implemented into ISIS. This tool was used to fit many spectra of X-ray binaries. However, the large number of degrees of freedom and individual datasets raise an issue about a good measure for a simultaneous fit quality. We present three ways to check the goodness of these fits: we investigate the goodness of each fit in all datasets, we define a combined goodness exploiting the logical structure of a simultaneous fit, and we stack the fit residuals of all datasets to detect weak features. These tools are applied to all RXTE-spectra from GRO 1008-57, revealing calibration features that are not detected significantly in any single spectrum. Stacking the residuals from the best-fit model for the Vela X-1 and XTE J1859+083 data evidences fluorescent emission lines that would have gone undetected otherwise.

The centroid energy Ecyc of the cyclotron line in the spectrum of the binary X-ray pulsar Her X-1 had been found to decrease with time on a time scale of a few tens of years - surprisingly short astrophysically. This was found for the pulse phase averaged line centroid energy using observational data from various X-ray satellites over the time frame 1996 to 2012, establishing a reduction of about 4 keV. Here we report on the result of a new observation by NuSTAR performed in August 2015. The earlier results are confirmed and strengthened, with respect to both, the dependence of Ecyc on flux (it is still present after 2006) and the dependence on time: the long-term decay continued with the same rate, corresponding to a reduction of about 5 keV in 20 years.

Galactic cosmic rays reach energies of at least a few Peta-electronvolts (1 PeV =$10^\mathbf{15}$ electron volts). This implies our Galaxy contains PeV accelerators (PeVatrons), but all proposed models of Galactic cosmic-ray accelerators encounter non-trivial difficulties at exactly these energies. Tens of Galactic accelerators capable of accelerating particle to tens of TeV (1 TeV =$10^\mathbf{12}$ electron volts) energies were inferred from recent gamma-ray observations. None of the currently known accelerators, however, not even the handful of shell-type supernova remnants commonly believed to supply most Galactic cosmic rays, have shown the characteristic tracers of PeV particles: power-law spectra of gamma rays extending without a cutoff or a spectral break to tens of TeV. Here we report deep gamma-ray observations with arcminute angular resolution of the Galactic Centre regions, which show the expected tracer of the presence of PeV particles within the central 10~parsec of the Galaxy. We argue that the supermassive black hole Sagittarius A* is linked to this PeVatron. Sagittarius A* went through active phases in the past, as demonstrated by X-ray outbursts and an outflow from the Galactic Centre. Although its current rate of particle acceleration is not sufficient to provide a substantial contribution to Galactic cosmic rays, Sagittarius A* could have plausibly been more active over the last $\gtrsim 10^{6-7}$ years, and therefore should be considered as a viable alternative to supernova remnants as a source of PeV Galactic cosmic rays.

We have performed a full time- and luminosity-resolved spectral analysis of the high-mass X-ray binary 4U 1538-522 using the available RXTE, INTEGRAL, and Suzaku data, examining both phase-averaged and pulse-phase-constrained datasets and focusing on the behavior of the cyclotron resonance scattering feature (CRSF). No statistically significant trend between the energy of the CRSF and luminosity is observed in the combined dataset. However, the CRSF energy appears to have increased by ~1.5 keV in the ~8.5 years between the RXTE and Suzaku measurements, with Monte Carlo simulations finding the Suzaku measurement 4.6$\sigma$ above the RXTE points. Interestingly, the increased Suzaku CRSF energy is much more significant and robust in the pulse-phase-constrained spectra from the peak of the main pulse, suggesting a change that is limited to a single magnetic pole. The 7 years of RXTE measurements do not show any strongly-significant evolution with time on their own. We discuss the significance of the CRSF's behavior with respect to luminosity and time in the context of historical observations of this source as well as recent observational and theoretical work concerning the neutron star accretion column, and suggest some mechanisms by which the observed change over time could occur.

Central compact objects (CCOs) are thought to be young thermally emitting isolated neutron stars that were born during the preceding core-collapse supernova explosion. Here we present evidence that at least in one case the CCO could have been formed within a binary system. We show that the highly reddened optical source IRAS~17287$-$3443, located $25^{\prime \prime}$ away from the CCO candidate XMMUJ173203.3$-$344518 and classified previously as a post asymptotic giant branch star, is indeed surrounded by a dust shell. This shell is heated by the central star to temperatures of $\sim90$\,K and observed as extended infrared emission in 8-160\,$\mu$m band. The dust temperature also increases in the vicinity of the CCO which implies that it likely resides within the shell. We estimate the total dust mass to be $\sim0.4-1.5\,M_\odot$ which significantly exceeds expected dust yields by normal stars and thus likely condensed from supernova ejecta. Taking into account that both the age of the supernova remnant and the duration of active mass loss phase by the optical star are much shorter than the total lifetime of either object, the supernova and the onset of the active mass loss phase of the companion have likely occurred approximately simultaneously. This is most easily explained if the evolution of both objects is interconnected. We conclude, therefore, that both stars were likely members of the same binary system disrupted by a supernova.

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.

Model spectra of neutron star atmospheres are nowadays widely used to fit the observed thermal X-ray spectra of neutron stars. This fitting is the key element in the method of the neutronstar radius determination. Here, we present the basic assumptions used for the neutron star atmosphere modeling as well as the main qualitative features of the stellar atmospheres leading to the deviations of the emergent model spectrum from blackbody. We describe the properties of two of our model atmosphere grids: (i) pure carbon atmospheres for relatively cool neutron stars (1--4 MK) and (ii) hot atmospheres with Compton scattering taken into account. The results obtained by applying these grids to model the X-ray spectra of the central compact object in supernova remnant HESS 1731-347, and two X-ray bursting neutron stars in low-mass X-ray binaries, 4U 1724-307 and 4U 1608-52, are presented. Possible systematic uncertainties associated with the obtained neutron star radii are discussed.

We present a timing and spectral analysis of the X-ray pulsar XTE J1946+274 observed with Suzaku during an outburst decline in 2010 October and compare with previous results. XTE J1946+274 is a transient X-ray binary consisting of a Be-type star and a neutron star with a 15.75 s pulse period in a 172 d orbit with 2-3 outbursts per orbit during phases of activity. We improve the orbital solution using data from multiple instruments. The X-ray spectrum can be described by an absorbed Fermi-Dirac cutoff power law model along with a narrow Fe K line at 6.4 keV and a weak Cyclotron Resonance Scattering Feature (CRSF) at ~35 keV. The Suzaku data are consistent with the previously observed continuum flux versus iron line flux correlation expected from fluorescence emission along the line of sight. However, the observed iron line flux is slightly higher, indicating the possibility of a higher iron abundance or the presence of non-uniform material. We argue that the source most likely has only been observed in the subcritical (non-radiation dominated) state since its pulse profile is stable over all observed luminosities and the energy of the CRSF is approximately the same at the highest (~$5 \times 10^{37} $erg s$^{-1}$) and lowest (~$5 \times 10^{36} $erg s$^{-1}$) observed 3-60 keV luminosities.

The central compact object in the supernova remnant HESS J1731-347 appears to be the hottest observed isolated cooling neutron star. The cooling theory of neutron stars enables one to explain observations of this star by assuming the presence of strong proton superfluidity in the stellar core and the existence of the surface heat blanketing envelope which almost fully consists of carbon. The cooling model of this star is elaborated to take proper account of the neutrino emission due to neutron-neutron collisions which is not suppressed by proton superfluidity. Using the results of spectral fits of observed thermal spectra for the distance of 3.2 kpc and the cooling theory for the neutron star of age 27 kyr, new constraints on the stellar mass and radius are obtained which are more stringent than those derived from the spectral fits alone.

Context: Be/X-ray binaries show outbursts with peak luminosities up to a few times $10^{37}\,$erg/s, during which they can be observed and studied in detail. Most (if not all) Be/X-ray binaries harbour accreting pulsars, whose X-ray spectra in many cases contain cyclotron resonant scattering features related to the magnetic field of the sources. Spectral variations as a function of luminosity and of the rotational phase of the neutron star are observed in many accreting pulsars. Aims: We explore X-ray spectral and timing properties of the Be/X-ray binary GX 304-1 during an outburst episode. Specifically, we investigate the behavior of the cyclotron resonant scattering feature, the continuum spectral parameters, the pulse period, and the energy- and luminosity-resolved pulse profiles. We combine the luminosity-resolved spectral and timing analysis to probe the accretion geometry and the beaming patterns of the rotating neutron star. Methods: We analyze the INTEGRAL data from the two JEM-X modules, ISGRI and SPI, covering the January-February 2012 outburst, divided in six observations. We obtain pulse profiles in two energy bands, phase-averaged and phase-resolved spectra for each observation. Results: We confirm the positive luminosity-dependence of the cyclotron line energy in GX 304-1, and report a dependence of the photon index on luminosity. Using a pulse-phase connection technique, we find a pulse period solution valid for the entire outburst. Our pulse-phase resolved analysis shows, that the centroid energy of the cyclotron line is varying only slightly with pulse phase, while other spectral parameters show more pronounced variations. Our results are consistent with a scenario in which, as the pulsar rotates, we are exploring only a small portion of its beam pattern.

Using RXTE/ASM archival data, we investigate the behaviour of the spectral hardness ratio as a function of X-ray luminosity in a sample of six transient X-ray pulsars (EXO 2030+375, GX 304-1, 4U 0115+63, V 0332+63, A 0535+26 and MXB 0656-072). In all sources we find that the spectral hardness ratio defined as $F_{5-12\mathrm{keV}}/ F_{1.33-3\mathrm{keV}}$ increases with the ASM flux (1.33--12 keV) at low luminosities and then saturates or even slightly decreases above some critical X-ray luminosity falling into the range $\sim(3-7)\times10^{37}$~erg~s$^{-1}$. Two-dimensional structure of accretion columns in the radiation-diffusion limit is calculated for two possible geometries (filled and hollow cylinder) for mass accretion rates $\dot M$ ranging from $10^{17}$ to 1.2$\times 10^{18}$~g s$^{-1}$. The observed spectral behaviour in the transient X-ray pulsars with increasing $\dot M$ can be reproduced by a Compton saturated sidewall emission from optically thick magnetized accretion columns with taking into account the emission reflected from the neutron star atmosphere. At $\dot M$ above some critical value $\dot M_{cr}\sim (6-8)\times 10^{17}$~g~s$^{-1}$, the hight of the column becomes such that the contribution of the reflected component to the total emission starts decreasing, which leads to the saturation and even slight decrease of the spectral hardness. Hollow-cylinder columns have a smaller height than the filled-cylinder ones, and the contribution of the reflected component in the total emission does not virtually change with $\dot M$ (and hence the hardness of the continuum monotonically increases) up to higher mass accretion rates than $\dot M_{cr}$ for the filled columns.

Context: The magnetic field is a crucial ingredient of neutron stars. It governs the physics of accretion and of the resulting high-energy emission in accreting pulsars. Studies of the cyclotron resonant scattering features (CRSFs) seen as absorption lines in the X-ray spectra of the pulsars permit direct measuremets of the field strength. Aims: From an analysis of a number of pointed observations with different instruments, the energy of CRSF, Ecyc, has recently been found to decay in Her X-1, which is one of the best-studied accreting pulsars. We present our analysis of a homogeneous and almost uninterrupted monitoring of the line energy with Swift/BAT. Methods: We analyzed the archival Swift/BAT observations of Her X-1 from 2005 to 2014. The data were used to measure the CRSF energy averaged over several months. Results: The analysis confirms the long-term decay of the line energy. The downward trend is highly significant and consistent with the trend measured with the pointed observations: dEcyc/dt ~-0.3 keV per year. Conclusions: The decay of Ecyc either indicates a local evolution of the magnetic field structure in the polar regions of the neutron star or a geometrical displacement of the line-forming region due to long-term changes in the structure of the X-ray emitting region. The shortness of the observed timescale of the decay, -Ecyc/(dEcyc/dt) ~ 100 yr, suggests that trend reversals and/or jumps of the line energy might be observed in the future.

Re-observations with the H.E.S.S. telescope array of the very-high-energy (VHE) source HESS J1018-589 A coincident with the Fermi-LAT $\gamma$-ray binary 1FGL J1018.6-5856 have resulted in a source detection significance of more than 9$\sigma$, and the detection of variability ($\chi^2$/$\nu$ of 238.3/155) in the emitted $\gamma$-ray flux. This variability confirms the association of HESS J1018-589 A with the high-energy $\gamma$-ray binary detected by Fermi-LAT, and also confirms the point-like source as a new very-high-energy binary system. The spectrum of HESS J1018-589 A is best fit with a power-law function with photon index $\Gamma = 2.20 \pm 0.14_{\rm stat} \pm 0.2_{\rm sys}$. Emission is detected up to ~20 TeV. The mean differential flux level is $(2.9 \pm 0.4)\times10^{-13}$ TeV$^{-1}$ cm$^{-2}$ s$^{-1}$ at 1 TeV, equivalent to ~1% of the flux from the Crab Nebula at the same energy. Variability is clearly detected in the night-by-night lightcurve. When folded on the orbital period of 16.58 days, the rebinned lightcurve peaks in phase with the observed X-ray and high-energy phaseograms. The fit of the H.E.S.S. phaseogram to a constant flux provides evidence of periodicity at the level of 3$\sigma$. The shape of the VHE phaseogram and measured spectrum suggest a low inclination, low eccentricity system with a modest impact from VHE $\gamma$-ray absorption due to pair production ($\tau$ < 1 at 300 GeV).