We present the analysis of 9.7 years Fermi-LAT data of the middle-aged supernova remnant W44 and the massive molecular gas complex that surrounds it. We derived a high-quality spectral energy distribution of gamma-radiation of the shell over three decades. The very hard spectrum below 1 GeV supports the earlier claims regarding the hadronic origin of radiation. We also confirm the presence of two extended $\gamma$-ray structures located at two opposite edges of the remnant along its major axis. Based on the high-resolution gas maps, we demonstrate that the gamma-ray structures are caused by the enhanced cosmic-ray density rather than the gradient of the gas distribution. We argue that the revealed cosmic-ray "clouds" suggest an anisotropic character of the escape of high-energy particles from the shell along the magnetic field of the remnant.

The nearby radio galaxy Centaurus A belongs to a class of Active Galaxies that are very luminous at radio wavelengths. The majority of these galaxies show collimated relativistic outflows known as jets, that extend over hundreds of thousands of parsecs for the most powerful sources. Accretion of matter onto the central super-massive black hole is believed to fuel these jets and power their emission, with the radio emission being related to the synchrotron radiation of relativistic electrons in magnetic fields. The origin of the extended X-ray emission seen in the kiloparsec-scale jets from these sources is still a matter of debate, although Cen A's X-ray emission has been suggested to originate in electron synchrotron processes. The other possible explanation is Inverse Compton (IC) scattering with CMB soft photons. Synchrotron radiation needs ultra-relativistic electrons ($\sim50$ TeV), and given their short cooling times, requires some continuous re-acceleration mechanism to be active. IC scattering, on the other hand, does not require very energetic electrons, but requires jets that stay highly relativistic on large scales ($\geq$1 Mpc) and that remain well-aligned with the line of sight. Some recent evidence disfavours inverse Compton-CMB models, although other evidence seems to be compatible with them. In principle, the detection of extended gamma-ray emission, directly probing the presence of ultra-relativistic electrons, could distinguish between these options, but instruments have hitherto been unable to resolve the relevant structures. At GeV energies there is also an unusual spectral hardening in Cen A, whose explanation is unclear. Here we report observations of Cen A at TeV energies that resolve its large-scale jet. We interpret the data as evidence for the acceleration of ultra-relativistic electrons in the jet, and favour the synchrotron explanation for the X-rays.

We report the detection of high-energy gamma-ray signal towards the young star-forming region, W40. Using 10-year Pass 8 data from the Fermi Large Area Telescope (Fermi-LAT), we extracted an extended gamma-ray excess region with a significance of about 18sigma. The radiation has a spectrum with a photon index of 2.49 +/- 0.01. The spatial correlation with the ionized gas content favors the hadronic origin of the gamma-ray emission. The total cosmic-ray (CR) proton energy in the gamma-ray production region is estimated to be the order of 10^47 erg. However, this could be a small fraction of the total energy released in cosmic rays (CRs) by local accelerators, presumably by massive stars, over the lifetime of the system. If so, W40, together with earlier detections of gamma-rays from Cygnus cocoon, Westerlund 1, Westerlund 2, NGC 3603, and 30 Dor C, supports the hypothesis that young star clusters are effective CR factories. The unique aspect of this result is that the gamma-ray emission is detected, for the first time, from a stellar cluster itself, rather than from the surrounding "cocoons".

The detection of the first electromagnetic counterpart to the binary neutron star (BNS) merger remnant GW170817 established the connection between short $\gamma$-ray bursts and BNS mergers. It also confirmed the forging of heavy elements in the ejecta (a so-called kilonova) via the r-process nucleosynthesis. The appearance of non-thermal radio and X-ray emission, as well as the brightening, which lasted more than 100 days, were somewhat unexpected. Current theoretical models attempt to explain this temporal behavior as either originating from a relativistic off-axis jet or a kilonova-like outflow. In either scenario, there is some ambiguity regarding how much energy is transported in the non-thermal electrons versus the magnetic field of the emission region. Combining the VLA (radio) and Chandra (X-ray) measurements with observations in the GeV-TeV domain can help break this ambiguity, almost independently of the assumed origin of the emission. Here we report for the first time on deep H.E.S.S. observations of GW170817 / GRB 170817A between 124 and 272 days after the BNS merger with the full H.E.S.S. array of telescopes, as well as on an updated analysis of the prompt (<5 days) observations with the upgraded H.E.S.S. phase-I telescopes. We discuss implications of the H.E.S.S. measurement for the magnetic field in the context of different source scenarios.

We report on the detection of very-high-energy (VHE; $E > 100$ GeV) $\gamma$-ray emission from the BL Lac objects KUV 00311-1938 and PKS 1440-389 with the High Energy Stereoscopic System (H.E.S.S.). H.E.S.S. observations were accompanied or preceded by multi-wavelength observations with Fermi/LAT, XRT and UVOT on board the Swift satellite, and ATOM. Based on an extrapolation of the Fermi/LAT spectrum towards the VHE $\gamma$-ray regime, we deduce a 95% confidence level upper limit on the unknown redshift of KUV 00311-1938 of z < 0.98, and of PKS 1440-389 of z < 0.53. When combined with previous spectroscopy results the redshift of KUV 00311-1938 is constrained to $0.51 \leq z < 0.98$ and for PKS 1440-389 to $0.14 \lessapprox z < 0.53$.

Analyzing Chandra data of Tycho's supernova remnant (SNR) taken in 2000, 2003, 2007, 2009, and 2015, we search for time variable features of synchrotron X-rays in the southwestern part of the SNR, where stripe structures of hard X-ray emission were previous found. By comparing X-ray images obtained at each epoch, we discover a knot-like structure in the northernmost part of the stripe region became brighter particularly in 2015. We also find a bright filamentary structure gradually became fainter and narrower as it moved outward. Our spectral analysis reveal that not only the nonthermal X-ray flux but also the photon indices of the knot-like structure change from year to year. During the period from 2000 to 2015, the small knot shows brightening of $\sim 70\%$ and hardening of $\Delta \Gamma \sim 0.45$. The time variability can be explained if the magnetic field is amplified to $\sim 100~\mathrm{\mu G}$ and/or if magnetic turbulence significantly changes with time.

In nuclei of starburst galaxies, the combination of an enhanced rate of supernova explosions and a high gas density suggests that cosmic rays can be efficiently produced, and that most of them lose their energy before escaping these regions, resulting in a large flux of secondary products, including neutrinos. Although the flux inferred from an individual starburst region is expected to be well below the sensitivity of current neutrino telescopes, such sources may provide a substantial contribution to the diffuse neutrino flux measured by IceCube. Here we compute the gamma-ray and neutrino flux due to starburst galaxies based on a physical model of cosmic ray transport in a starburst nucleus, and accounting for the redshift evolution of the number density of starburst sources as inferred from recent measurements of the star formation rate. The model accounts for gamma-ray absorption both inside the sources and in the intergalactic medium. The latter process is responsible for electromagnetic cascades, which also contribute to the diffuse gamma-ray background at lower energies. The conditions for acceleration of cosmic ray protons up to energies exceeding $ \sim 10 \, \rm PeV$ in starburst regions, necessary for the production of PeV neutrinos, are investigated in a critical way. We show that starburst nuclei can account for the diffuse neutrino flux above $\sim 200 \, \rm TeV$, thereby producing $\lesssim 40 \%$ of the extragalactic diffuse gamma-ray background. Below $\sim 200 \, \rm TeV$, the flux from starburst appears to be somewhat lower than the observed one, where both the Galactic contribution and the flux of atmospheric neutrinos may account for the difference.

The synchrotron cooling of relativistic electrons is one of the most effective radiation mechanisms in astrophysics. It not only accompanies the process of particle acceleration but also has feedback on the formation of the energy distribution of the parent electrons. The radiative cooling time of electrons decreases with energy as $t_{\rm syn} \propto 1/E$; correspondingly the overall radiation efficiency increases with energy. On the other hand, this effect strictly limits the maximum energy of individual photons. Even in the so-called extreme accelerators, where the acceleration proceeds at the highest possible rate, $t_{\rm acc} = E/eBc$, the synchrotron radiation cannot extend well beyond the characteristic energy determined by the electron mass and the fine-structure constant: $h \nu_{\rm max} \sim m_e c^2/\alpha \sim 100 \rm\,MeV$. In a highly turbulent medium with a broad magnetic field distribution, this limit is shifted towards higher energies. However, for any realistic distribution of the magnetic field, it cannot significantly exceed $ h \nu_{\rm max}$. In this paper, we propose a model in which the formation of synchrotron radiation takes place in compact magnetic blobs located inside the particle accelerator. We develop a formalism for calculations of synchrotron radiation emerging from such systems. We demonstrate that for certain combinations of parameters characterizing the accelerator and the magnetic blobs, the synchrotron radiation can extend into the very high energy (VHE) domain, $h \nu_{\rm max} \geq 100 \ \rm GeV$. We briefly outline the possible astrophysical implications of the results in the context of VHE radiation of blazars and gamma-ray bursts.

We demonstrate that the exact solution for the spectrum of synchrotron radiation from an isotropic population of mono-energetic electrons in turbulent magnetic field with Gaussian distribution of local field strengths can be expressed in the simple analytic form: $\left( \frac{{\rm d} \dot{N}}{{\rm d} \omega} \right)_t = \frac{\alpha}{3} \frac{1}{\gamma^2} \left( 1 + \frac{1}{x^{2/3}} \right) \exp \left( - 2 x^{2/3} \right)$, where $x = \frac{\omega}{\omega_0}\, ; \omega_0 = \frac{4}{3} \gamma^2 \frac{eB_0}{m_e c}\, .$ We use this expression to find approximate synchrotron spectra for power-law electron distributions with $\propto \exp\left( -\left[ \gamma/\gamma_0 \right]^\beta\right)$ type high-energy cut-off; the resulting synchrotron spectrum has the exponential cut-off factor with frequency raised to $2\beta/(3\beta+4)$ power in the exponent. For the power-law electron distribution without high-energy cut-off, we find the coefficient $a_m$ as a function of the power-law index, which results in exact expression for the synchrotron spectrum when using monochromatic (i.e., each electron radiates at frequency $\omega_m = a_m \gamma^2 \, \frac{e B_0}{m_e c}$) approximation.

Aims. Colliding wind binary systems have long been suspected to be high-energy (HE; 100 MeV < E < 100 GeV) {\gamma}-ray emitters. {\eta} Car is the most prominent member of this object class and is confirmed to emit phase-locked HE {\gamma} rays from hundreds of MeV to ~100 GeV energies. This work aims to search for and characterise the very-high-energy (VHE; E >100 GeV) {\gamma}-ray emission from {\eta} Car around the last periastron passage in 2014 with the ground-based High Energy Stereoscopic System (H.E.S.S.). Methods. The region around {\eta} Car was observed with H.E.S.S. between orbital phase p = 0.78 - 1.10, with a closer sampling at p {\approx} 0.95 and p {\approx} 1.10 (assuming a period of 2023 days). Optimised hardware settings as well as adjustments to the data reduction, reconstruction, and signal selection were needed to suppress and take into account the strong, extended, and inhomogeneous night sky background (NSB) in the {\eta} Car field of view. Tailored run-wise Monte-Carlo simulations (RWS) were required to accurately treat the additional noise from NSB photons in the instrument response functions. Results. H.E.S.S. detected VHE {\gamma}-ray emission from the direction of {\eta} Car shortly before and after the minimum in the X-ray light-curve close to periastron. Using the point spread function provided by RWS, the reconstructed signal is point-like and the spectrum is best described by a power law. The overall flux and spectral index in VHE {\gamma} rays agree within statistical and systematic errors before and after periastron. The {\gamma}-ray spectrum extends up to at least ~400 GeV. This implies a maximum magnetic field in a leptonic scenario in the emission region of 0.5 Gauss. No indication for phase-locked flux variations is detected in the H.E.S.S. data.

We present results of spectrometric studies based on the observations of very strong 3C 279 flares in high and very-high-energy bands and discuss their implications regarding the origin of radiation mechanisms. The FSRQ 3C 279 (z=0.536) is one of the most luminous gamma-ray emitting AGN. It shows variability on time scales down to minutes during strong flares detected by the Fermi-LAT. We have analyzed LAT and Swift-XRT data for the flaring periods in June 2015 and January 2018, when the flux above 100 MeV during hourly time intervals could exceed $3 \times 10^{-5}$ ph cm$^{-2}$ s$^{-1}$. The X-ray flux derived from Swift X-ray observations described by a very hard spectrum with photon index ~1.5 is typically explained by Inverse Compton scattering of low energy electrons and protons. Here we consider an alternative interpretation which describes the entire band from X-ray to very high energy as a results of synchrotron emission by ultra relativistic electrons or protons.

PSR B1259-63/LS 2883 is a gamma-ray binary system consisting of a pulsar in an eccentric orbit around a bright Oe stellar-type companion star that features a dense circumstellar disc. The high- and very-high-energy (HE, VHE) gamma-ray emission from PSR B1259-63/LS 2883 around the times of its periastron passage are characterised, in particular, at the time of the HE gamma-ray flares reported to have occurred in 2011, 2014, and 2017. Spectra and light curves were derived from observations conducted with the H.E.S.S.-II array in 2014 and 2017. A local double-peak profile with asymmetric peaks in the VHE light curve is measured, with a flux minimum at the time of periastron $t_p$ and two peaks coinciding with the times at which the neutron star crosses the companion's circumstellar disc ($\sim t_p \pm 16$ d). A high VHE gamma-ray flux is also observed at the times of the HE gamma-ray flares ($\sim t_p + 30$ d) and at phases before the first disc crossing ($\sim t_p - 35$ d). PSR B1259-63/LS 2883 displays periodic flux variability at VHE gamma-rays without clear signatures of super-orbital modulation in the time span covered by H.E.S.S. observations. In contrast, the photon index of the measured power-law spectra remains unchanged within uncertainties for about 200 d around periastron. Lower limits on exponential cut-off energies up to $\sim 40$ TeV are placed. At HE gamma-rays, PSR B1259-63/LS 2883 has now been detected also before and after periastron, close to the disc crossing times. Repetitive flares with distinct variability patterns are detected in this energy range. Such outbursts are not observed at VHEs, although a relatively high emission level is measured. The spectra obtained in both energy regimes displays a similar slope, although a common physical origin either in terms of a related particle population, emission mechanism, or emitter location is ruled out.

We present a numerical study of the properties of the flow produced by the collision of a magnetized anisotropic pulsar wind with the circumbinary environment. We focus on studying the impact of the high wind magnetization on the geometrical structure of the shocked flow. This work is an extension of our earlier studies that focused on a purely hydrodynamic interaction and weak wind magnetization. We consider the collision in the axisymmetric approximation, that is, the pulsar rotation axis is assumed to be oriented along the line between the pulsar and the optical star. The increase of the magnetization results in the expansion of the opening cone in which the shocked pulsar wind propagates. This effect is explained in the frameworks of the conventional theory of collimation of magnetized winds. This finding has a direct implication for scenarios that involve Doppler boosting as the primary mechanism behind the GeV flares detected with the Fermi/LAT from PSR B1259-63/LS2883. The maximum enhancement of the apparent emission is determined by the ratio of $4\pi$ to the solid in which the shocked pulsar wind propagates. Our simulations suggest that this enhancement factor is decreased by the impact of the magnetic field.

The Crab Nebula is an extreme particle accelerator boosting the energy of electrons up to a few PeV ($10^{15} \ \rm eV$), close to the maximum energy allowed by theory. The physical conditions in the acceleration site and the nature of the acceleration process itself remain highly uncertain. The key information about the highest energy accelerated particles is contained in the synchrotron and inverse Compton (IC) channels of radiation at energies above 1~MeV and 100~TeV, respectively. The recent report of detection of ultra-high energy gamma-ray signal from the Crab Nebula up to 300~TeV allows one to determine the energy distribution of the highest energy electrons and to derive the magnetic field strength in the acceleration region, $B\leq120\rm\,\mu G$, in a parameter-free way. This estimate brings new constraints on the properties of non-thermal particle distributions and puts important constraints on the MHD models for the Crab Nebula, in particular on the feasible magnetization and anisotropy of the pulsar wind. The calculations of synchrotron and IC emission show that future observations with instruments allowing detection of the Crab Nebula above 300~TeV and above 1~MeV will clarify the conditions allowing acceleration of electrons beyond PeV energies in the Crab Nebula. In particular, one will (1) verify the hypothetical multi-component composition of the electron energy distribution, and (2) determine the magnetic field strength in the regions responsible for the acceleration of PeV electrons.

Flat-spectrum radio-quasars (FSRQs) are rarely detected at very-high-energies (VHE; E>100 GeV) due to their low-frequency-peaked SEDs. At present, only 6 FSRQs are known to emit VHE photons, representing only 7% of the VHE extragalactic catalog. Following the detection of MeV-GeV gamma-ray flaring activity from the FSRQ PKS 0736+017 (z=0.189) with Fermi, the H.E.S.S. array of Cherenkov telescopes triggered ToO observations on February 18, 2015, with the goal of studying the gamma-ray emission in the VHE band. H.E.S.S. ToO observations were carried out during the nights of February 18, 19, 21, and 24, 2015. Together with Fermi-LAT, the multi-wavelength coverage of the flare includes Swift observations in soft-X-rays and optical/UV, and optical monitoring (photometry and spectro-polarimetry) by the Steward Observatory, the ATOM, the KAIT and the ASAS-SN telescope. VHE emission from PKS 0736+017 was detected with H.E.S.S. during the night of February 19, 2015, only. Fermi data indicate the presence of a gamma-ray flare, peaking at the time of the H.E.S.S. detection, with a flux doubling time-scale of around six hours. The gamma-ray flare was accompanied by at least a 1 mag brightening of the non-thermal optical continuum. No simultaneous observations at longer wavelengths are available for the night of the H.E.S.S. detection. The gamma-ray observations with H.E.S.S. and Fermi are used to put constraints on the location of the gamma-ray emitting region during the flare: it is constrained to be just outside the radius of the broad-line-region with a bulk Lorentz factor $\simeq 20$, or at the level of the radius of the dusty torus with Gamma > 60. PKS 0736+017 is the seventh FSRQ known to emit VHE photons and, at z=0.189, is the nearest so far. The location of the gamma-ray emitting region during the flare can be tightly constrained thanks to opacity, variability, and collimation arguments.

Cosmic ray (CR) spectra, both measured upon their arrival at the Earth's atmosphere and inferred from the emission in supernova remnants (SNRs), appear to be significantly steeper than the "standard" diffusive shock acceleration (DSA) theory predicts. Although the reconstruction of the primary spectra introduces an additional steepening due to propagation effects, there is a growing consensus in the CR community that these corrections fall short to explain the newest high-precision data. Using 2D hybrid simulations, we investigate a new mechanism that may steepen the spectrum during the acceleration in SNR shocks. Most of the DSA treatments are limited to homogeneous shock environments. To investigate whether inhomogeneity effects can produce the necessary extra steepening, we assume that the magnetic field changes its angle along the shock front. The rationale behind this approach is the strong dependence of the DSA efficiency upon the field angle, $\theta_\mathrm{Bn}$. Our results show that the variation of shock obliquity along its face results in a noticeable steepening of the DSA spectrum. Compared to simulations of quasi-parallel shocks, we observe an increase of the spectral index by $\Delta q=0.1-0.15$. Possible extrapolation of the limited simulation results to more realistic SNR conditions are briefly considered.

High energy $\gamma$ rays from Giant Molecular Clouds (GMCs) carry direct information about the spatial and energy distributions of Galactic Cosmic Rays (CRs). The recently released catalogs of GMCs contain sufficiently massive clouds to be used as barometers for probing, through their $\gamma$-ray emission, the density of CRs throughout the Galactic Disk. Based on the data of \fermi{}, we report the discovery of $\gamma$-ray signals from nineteen GMCs located at distances up to 12.5 kpc. The galactocentric radial distribution of the CR density derived from the $\gamma$-ray and CO observations of these objects, as well as from some nearby clouds that belong to the Gould Belt complex, unveil a homogeneous \textquotedblleft sea" of CRs with a constant density and spectral shape close to the flux of directly (locally) measured CRs. We found noticeable deviations from the \textquotedblleft sea level" only in some locations characterized by enhanced CR density in the galactocentric 4--6 kpc ring. Furthermore, we found a hint for fluctuations of the CR density in different locations within the same 4--6 kpc ring. The confirmation of this result with the next-generation $\gamma$-ray detectors based on the higher quality data and denser coverage of galactocentric distances, would have dramatic implications for the understanding of the origin of Galactic CRs.

The escape process of particles accelerated at supernova remnant (SNR) shocks is one of the poorly understood aspects of the shock acceleration theory. Here we adopt a phenomenological approach to study the particle escape and its impact on the gamma-ray spectrum resulting from hadronic collisions both inside and outside of a middle-aged SNR. Under the assumption that in the spatial region immediately outside of the remnant the diffusion coefficient is suppressed with respect to the average Galactic one, we show that a significant fraction of particles are still located inside the SNR long time after their nominal release from the acceleration region. This fact results into a gamma-ray spectrum that resembles a broken power law, similar to those observed in several middle-aged SNRs. Above the break, the spectral steepening is determined by the diffusion coefficient outside of the SNR and by the time dependence of maximum energy. Consequently, the comparison between the model prediction and actual data will contribute to determining these two quantities, the former being particularly relevant within the predictions of the gamma-ray emission from the halo of escaping particles around SNRs which could be detected with future Cherenkov telescope facilities. We also calculate the spectrum of run-away particles injected into the Galaxy by an individual remnant. Assuming that the acceleration stops before the SNR enters the snowplow phase, we show that the released spectrum can be a featureless power law only if the accelerated spectrum has a slope alpha > 4.

We argue that the anomalously high fluxes of positrons and antiprotons found in cosmic rays (CR) can be satisfactorily explained by introducing two additional elements to the current "standard" paradigm of Galactic CRs. First, we propose that the antiparticles are effectively produced in interactions of primary CRs with the surrounding gas not only in the interstellar medium (ISM) but also inside the accelerators. Secondly, we postulate the existence of two source populations injecting CRs into the ISM with different, (1) soft (close to $FI \propto E^{-2.3}$) and (2) hard ($FII \propto E^{-1.8}$ or harder), energy distributions. Assuming that CRs in the 2nd population of accelerators accumulate "grammage" of the order of $1 \ \rm g/cm^2$ before their leakage into ISM, we can explain the energy distributions and absolute fluxes of both positrons and antiprotons, as well as the fluxes of secondary nuclei of the (Li,Be,B) group. The superposition of contributions of two source populations also explains the reported hardening of the spectra of CR protons and nuclei above 200 GV. The 2nd source population accelerating CRs with a rate at the level below 10 percent of the power of the 1st source population, can be responsible for the highest energy protons and nuclei of Galactic CRs up to the "knee" around $10^{15} \ \rm eV$.

The Crab nebula is one of the most studied cosmic particle accelerators, shining brightly across the entire electromagnetic spectrum up to very high-energy gamma rays. It is known from radio to gamma-ray observations that the nebula is powered by a pulsar, which converts most of its rotational energy losses into a highly relativistic outflow. This outflow powers a pulsar wind nebula (PWN), a region of up to 10~light-years across, filled with relativistic electrons and positrons. These particles emit synchrotron photons in the ambient magnetic field and produce very high-energy gamma rays by Compton up-scattering of ambient low-energy photons. While the synchrotron morphology of the nebula is well established, it was up to now not known in which region the very high-energy gamma rays are emitted. Here we report that the Crab nebula has an angular extension at gamma-ray energies of 52 arcseconds (assuming a Gaussian source width), significantly larger than at X-ray energies. This result closes a gap in the multi-wavelength coverage of the nebula, revealing the emission region of the highest energy gamma rays. These gamma rays are a new probe of a previously inaccessible electron and positron energy range. We find that simulations of the electromagnetic emission reproduce our new measurement, providing a non-trivial test of our understanding of particle acceleration in the Crab nebula.

We present a detailed study of the gamma-ray emission from the direction of the star-forming region W30 based on a decade of the Fermi Large Area Telescope data in the 0.3-300 GeV photon energy range. The morphological and spectral analyses allow us to resolve the gamma-ray emission into two extended structures from different origins. One of them mostly overlaps with the supernova remnant (SNR) G8.7-0.1 and has a soft spectrum that resembles with the spectra of other middle-aged SNRs interacting with molecular clouds. The other shows remarkable spatial and spectral consistency with the TeV emission from HESS J1804-216, and its spectrum could be naturally explained by inverse Compton scattering of electrons like a number of TeV gamma-ray emitting pulsar wind nebulae. Thus we attribute this source to the nebula around the pulsar PSR J1803-2137.

The direct measurements of cosmic rays (CRs), after correction for the propagation effects in the interstellar medium, indicate that their source spectra are likely to be significantly steeper than the canonical $E^{-2}$ spectrum predicted by the standard Diffusive Shock Acceleration (DSA) mechanism. The DSA has long been held responsible for the production of galactic CRs in supernova remnant (SNR) shocks. The $\gamma$-ray 'probes' of the acceleration spectra of CRs on-the-spot, inside of the SNRs, lead to the same conclusion. We show that the steep acceleration spectrum can be attributed to the $combination$ of (i) spherical expansion, (ii) tilting of the magnetic field along the shock surface and (iii) shock deceleration. Because of (i) and (ii), the DSA is efficient only on two ``polar caps'' of a spherical shock where the local magnetic field is within $\simeq45^{\circ}$ to its normal. The shock-produced spectrum observed edge-on steepens with the particle energy because the number of freshly accelerated particles with lower energies continually adds up to a growing acceleration region. We demonstrate the steepening effect by obtaining an exact self-similar solution for the particle acceleration at expanding shock surface with an arbitrary energy dependence of particle diffusivity $\kappa$. We show that its increase toward higher energy steepens the spectrum, which deeply contrasts with the standard DSA spectrum where $\kappa$ cancels out.

The flat spectrum radio quasar 3C 279 is known to exhibit pronounced variability in the high-energy ($100\,$MeV$<E<100\,$GeV) $\gamma$-ray band, which is continuously monitored with Fermi-LAT. During two periods of high activity in April 2014 and June 2015 Target-of-Opportunity observations were undertaken with H.E.S.S. in the very-high-energy (VHE, $E>100\,$GeV) $\gamma$-ray domain. While the observation in 2014 provides an upper limit, the observation in 2015 results in a signal with $8.7\,\sigma$ significance above an energy threshold of $66\,$GeV. No VHE variability has been detected during the 2015 observations. The VHE photon spectrum is soft and described by a power-law index of $4.2\pm 0.3$. The H.E.S.S. data along with a detailed and contemporaneous multiwavelength data set provide constraints on the physical parameters of the emission region. The minimum distance of the emission region from the central black hole is estimated using two plausible geometries of the broad-line region and three potential intrinsic spectra. The emission region is confidently placed at $r\gtrsim 1.7\times10^{17}\,$cm from the black hole, i.e., beyond the assumed distance of the broad-line region. Time-dependent leptonic and lepto-hadronic one-zone models are used to describe the evolution of the 2015 flare. Neither model can fully reproduce the observations, despite testing various parameter sets. Furthermore, the H.E.S.S. data are used to derive constraints on Lorentz invariance violation given the large redshift of 3C 279.

Observations from the radio to the gamma-ray wavelengths indicate that supernova remnant (SNR) shocks are sites of effective particle acceleration. It has been proposed that the presence of dense clumps in the environment where supernovae explode might have a strong impact in shaping the hadronic gamma-ray spectrum. Here we present a detailed numerical study about the penetration of relativistic protons into clumps which are engulfed by a SNR shock, taking into account the magneto-hydrodynamical properties of the background plasma. We show that the spectrum of protons inside clumps is much harder than that in the diffuse inter-clump medium and we discuss the implications for the formation of the spectrum of hadronic gamma rays, which does not reflect anymore the acceleration spectrum of protons, resulting substantially modified inside the clumps due to propagation effects. For the Galactic SNR RX J1713.7-3946, we show that a hadronic scenario including dense clumps inside the remnant shell is able to reproduce the broadband gamma-ray spectrum from GeV to TeV energies. Moreover, we argue that small clumps crossed by the shock could provide a natural explanation to the non-thermal X-ray variability observed in some hot spots of RX J1713.7-3946. Finally we discuss the detectability of gamma-ray emission from clumps with the upcoming Cherenkov Telescope Array and the possible detection of the clumps themselves through molecular lines.

We investigate the interstellar medium (ISM) towards seven TeV gamma-ray sources thought to be pulsar wind nebulae (PWNe) using Mopra molecular line observations at 7mm [CS(1-0), SiO(1-0,v=0)], Nanten CO(1-0) data and the SGPS/GASS HI survey. We have discovered several dense molecular clouds co-located to these TeV gamma-ray sources , which allows us to search for cosmic-rays (CRs) coming from progenitor SNRs or, potentially, from PWNe. We notably found SiO(1-0,v=0) emission towards HESS J1809-193, highlighting possible interaction between the adjacent supernova remnant SNR G011.0-0.0 and the molecular cloud at d $\sim$ 3.7 kpc. Using morphological features, and comparative studies of our column densities with those obtained from X-ray measurements, we claim a distance d $\sim$ 8.6 - 9.7 kpc for SNR G292.2-00.5, d $\sim$ 3.5 - 5.6 kpc for PSR J1418-6058 and d $\sim$ 1.5 kpc for the new SNR candidate found towards HESS J1303-631. From our mass and density estimates of selected molecular clouds, we discuss signatures of hadronic/leptonic components from PWNe and their progenitor SNRs. Interestingly, the molecular gas, which overlaps HESS J1026-582 at d $\sim$ 5 kpc, may support a hadronic origin. We find however that this scenario requires an undetected cosmic-ray accelerator to be located at d $\lt$ 10 pc from the molecular cloud. For HESS J1809-193, the cosmic-rays which have escaped SNR G011.0-0.0 could contribute to the TeV gamma-ray emission. Finally, from the hypothesis that at most 20% the pulsar spin down power could be converted into CRs, we find that, among the studied PWNe, only those from PSR J1809-1917 could potentially contribute to the TeV emission.

The shock waves of supernova remnants (SNRs) are prominent candidates for the acceleration of the Galactic cosmic rays. SNR RX J1713.7$-$3946 is one well-studied particle accelerator in our Galaxy because of its strong non-thermal X-ray and gamma-ray radiation. We have performed NuSTAR (3-79 keV) observations of the northwest rim of RX J1713.7$-$3946, where is the brightest part in X-ray and the shock speed is about 4000 km s$^{-1}$. The spatially resolved X-ray emission from RX J1713.7$-$3946 is detected up to 20 keV for the first time. The hard X-ray image in 10-20 keV is broadly similar to the soft-band image in 3-10 keV. The typical spectrum is described by power-law model with exponential cutoff with the photon index $\Gamma$=2.15 and the cutoff energy $\varepsilon_c$=18.8 keV. Using a synchrotron radiation model from accelerated electrons in the loss-limited case, the cutoff energy parameter ranges 0.6-1.9 keV, varying from region to region. Combined with the previous measurement of the shock speed, the acceleration of electrons is close to the Bohm-limit regime in the outer edge, while the standard picture of accelerated particles limited by synchrotron radiation in SNR shock is not applicable in the inner edge and the filamentary structure.

The high rate of star formation and supernova explosions of starburst galaxies make them interesting sources of high energy radiation. Depending upon the level of turbulence present in their interstellar medium, the bulk of cosmic rays produced inside starburst galaxies may lose most of their energy before escaping, thereby making these sources behave as calorimeters, at least up to some maximum energy. Contrary to previous studies, here we investigate in detail the conditions under which cosmic ray confinement may be effective for electrons and nuclei and we study the implications of cosmic ray confinement in terms of multifrequency emission from starburst nuclei and production of high energy neutrinos. The general predictions are then specialized to three cases of active starbursts, namely M82, NGC253 and Arp220. Both primary and secondary electrons, as well as electron-positron pairs produced by gamma ray absorption inside starburst galaxies are taken into account. Electrons and positrons produced as secondary products of hadronic interactions are found to be responsible for most of the emission of leptonic origin. In particular, synchrotron emission of very high energy secondary electrons produces an extended emission of hard X-rays that represent a very interesting signature of hadronic process in starburst galaxies, potentially accessible to current and future observations in the X-ray band. A careful understanding of both the production and absorption of gamma rays in starburst galaxies is instrumental to the assessment of the role of these astrophysical sources as sources of high energy astrophysical neutrinos.

Young core-collapse supernovae with dense-wind progenitors may be able to accelerate cosmic-ray hadrons beyond the knee of the cosmic-ray spectrum, and this may result in measurable gamma-ray emission. We searched for gamma-ray emission from ten supernovae observed with the High Energy Stereoscopic System (H.E.S.S.) within a year of the supernova event. Nine supernovae were observed serendipitously in the H.E.S.S. data collected between December 2003 and December 2014, with exposure times ranging from 1.4 hours to 53 hours. In addition we observed SN 2016adj as a target of opportunity in February 2016 for 13 hours. No significant gamma-ray emission has been detected for any of the objects, and upper limits on the $>1$ TeV gamma-ray flux of the order of $\sim$10$^{-13}$ cm$^{-2}$s$^{-1}$ are established, corresponding to upper limits on the luminosities in the range $\sim$2 $\times$ 10$^{39}$ erg s$^{-1}$ to $\sim$1 $\times$ 10$^{42}$ erg s$^{-1}$. These values are used to place model-dependent constraints on the mass-loss rates of the progenitor stars, implying upper limits between $\sim$2 $\times 10^{-5}$ and $\sim$2 $\times 10^{-3}$M$_{\odot}$yr$^{-1}$ under reasonable assumptions on the particle acceleration parameters.

PKS 1830-211 is a known macrolensed quasar located at a redshift of z=2.5. Its high-energy gamma-ray emission has been detected with the Fermi-LAT instrument and evidence for lensing was obtained by several authors from its high-energy data. Observations of PKS 1830-211 were taken with the H.E.S.S. array of Imaging Atmospheric Cherenkov Telescopes in August 2014, following a flare alert by the Fermi- LAT collaboration. The H.E.S.S observations were aimed at detecting a gamma-ray flare delayed by 20-27 days from the alert flare, as expected from observations at other wavelengths. More than twelve hours of good quality data were taken with an analysis threshold of $\sim67$ GeV. The significance of a potential signal is computed as a function of the date as well as the average significance over the whole period. Data are compared to simultaneous observations by Fermi-LAT. No photon excess or significant signal is detected. An upper limit on PKS 1830-211 flux above 67 GeV is computed and compared to the extrapolation of the Fermi-LAT flare spectrum.

Optically thin accretion plasmas can reach ion temperatures $T_{\rm i} \geq 10^{10}$K and thus trigger nuclear reactions. Using a large nuclear interactions network, we studied the radial evolution of the chemical composition of the accretion flow toward the black hole and computed the emissivity in nuclear $\gamma$-ray lines. In the advection dominated accretion flow (ADAF) regime, CNO and heavier nuclei are destroyed before reaching the last stable orbit. The overall luminosity in the de-excitation lines for a solar composition of plasma can be as high as few times $10^{-5}$ the accretion luminosity ($\dot{M}c^2$) and can be increased for heavier compositions up to $10^{-3}$. The efficiency of transformation of the kinetic energy of the outflow into high energy ($\geq 100$~MeV) $\gamma$-rays through the production and decay of $\pi^0$-mesons can be higher, up to $10^{-2}$ of the accretion luminosity. We show that in the ADAF model up to 15 percent of the mass of accretion matter can `evaporate' in the form of neutrons.

Very hot plasmas with ion temperature exceeding $10^{10}$ K can be formed in certain astrophysical environments. The distinct radiation signature of such plasmas is the $\gamma$-ray emission dominated by the prompt de-excitation nuclear lines and $\pi^0$-decay $\gamma$-rays. Using a large nuclear reaction network, we compute the time evolution of the chemical composition of such hot plasmas and their $\gamma$-ray line emissivity. At higher energies, we provide simple but accurate analytical presentations for the $\pi^0$-meson production rate and the corresponding $\pi^0\to2\gamma$ emissivity derived for the Maxwellian distribution of protons. We discuss the impact of the possible deviation of the high energy tail of the particle distribution function from the "nominal" Maxwellian distribution on the plasma $\gamma$-ray emissivity.

Sources of X-rays such as active galactic nuclei and X-ray binaries are often variable by orders of magnitude in luminosity over timescales of years. During and after these flares the surrounding gas is out of chemical and thermal equilibrium. We introduce a new implementation of X-ray radiative transfer coupled to a time-dependent chemical network for use in 3D magnetohydrodynamical simulations. A static fractal molecular cloud is irradiated with X-rays of different intensity, and the chemical and thermal evolution of the cloud are studied. For a simulated $10^5$ M$_\odot$ fractal cloud an X-ray flux $<0.01$ erg cm$^{-2}$ s$^{-1}$ allows the cloud to remain molecular, whereas most of the CO and H$_2$ are destroyed for a flux of $>1$ erg cm$^{-2}$ s$^{-1}$. The effects of an X-ray flare, which suddenly increases the X-ray flux by $10^5 \times$ are then studied. A cloud exposed to a bright flare has 99% of its CO destroyed in 10-20 years, whereas it takes $>10^3$ years for 99% of the H$_2$ to be destroyed. CO is primarily destroyed by locally generated far-UV emission from collisions between non-thermal electrons and H$_2$; He$^+$ only becomes an important destruction agent when the CO abundance is already very small. After the flare is over, CO re-forms and approaches its equilibrium abundance after $10^3-10^5$ years. This implies that molecular clouds close to Sgr A$^*$ in the Galactic Centre may still be out of chemical equilibrium, and we predict the existence of clouds near flaring X-ray sources in which CO has been mostly destroyed but H is fully molecular.

The identification of major contributors to the locally observed fluxes of Cosmic Rays (CRs) is a prime objective towards the resolution of the long-standing enigma of CRs. We report on a compelling similarity of the energy and radial distributions of multi-TeV CRs extracted from observations of very high energy (VHE) $\gamma$-rays towards the Galactic Center (GC) and two prominent clusters of young massive stars, Cyg~OB2 and Westerlund~1. This resemblance we interpret as a hint that CRs responsible for the diffuse VHE $\gamma$-ray emission from the GC are accelerated by the ultracompact stellar clusters located in the heart of GC. The derived $1/r$ decrement of the CR density with the distance from a star cluster is a distinct signature of continuous, over a few million years, CR injection into the interstellar medium. The lack of brightening of the $\gamma$-ray images toward the stellar clusters excludes the leptonic origin of $\gamma$-radiation. The hard, $\propto E^{-2.3}$ type power-law energy spectra of parent protons continues up to $\sim$ 1 PeV. The efficiency of conversion of kinetic energy of stellar winds to CRs can be as high as 10 percent implying that the young massive stars may operate as proton PeVatrons with a dominant contribution to the flux of highest energy galactic CRs.

We investigate the nature of the accelerated particles responsible for the production of the gamma-ray emission observed from the middle-aged supernova remnant (SNR) HB 21. The analysis of more than nine years of Fermi LAT data leads to the observation of an extended emission positionally in agreement with the SNR HB 21. The bulk of this gamma-ray emission is detected from the remnant; photons up to $\sim$10\,GeV show clear evidence of curvature at the lower energies. The remnant is characterized by an extension of $0^{\circ}.83$, that is, 30\% smaller than claimed in previous studies. The increased statistics allows us also to resolve a point-like source at the edge of the remnant, in proximity to a molecular cloud of the Cyg OB7 complex. In the southern part of the remnant, a hint of an additional gamma-ray excess in correspondence to shocked molecular clouds is observed. The spectral energy distribution of the SNR shows evidence of a break around 400 MeV, which can be properly fitted within both the hadronic and leptonic scenario. The pion-decay mechanism reproduces well the gamma rays, postulating a proton spectrum with a slope $\sim 2.5$ and with a steepening around tens of GeV, which could be explained by the energy-dependent escape of particles from the remnant. In the leptonic scenario the electron spectrum within the SNR matches closely the locally measured spectrum. This remarkable and novel result shows that SNR HB 21 could be a direct contributor to the population of Galactic electrons. In the leptonic scenario, we find that the local electron spectrum with a break around 2 GeV, closely evokes the best-fitting parental spectrum within this SNR. If such a scenario is confirmed, this would indicate that the SNR might be a source of Galactic background electrons.

The detection of the pulsed $\sim 1 $~TeV gamma-ray emission from the Crab pulsar reported by MAGIC and VERITAS collaborations demands a substantial revision of existing models of particle acceleration in the pulsar magnetosphere. In this regard model independent restrictions on the possible production site of the VHE photons become an important issue. In this paper, we consider limitations imposed by the process of conversion of VHE gamma rays into $e^{\pm}$ pairs in the magnetic field of the pulsar magnetosphere. Photons with energies exceeding 1~TeV are effectively absorbed even at large distances from the surface of the neutron star. Our calculations of magnetic absorption in the force-free magnetosphere show that the twisting of the magnetic field due to the pulsar rotation makes the magnetosphere more transparent compared to the dipole magnetosphere. The gamma-ray absorption appears stronger for photons emitted in the direction of rotation than in the opposite direction. There is a small angular cone inside which the magnetosphere is relatively transparent and photons with energy $1.5$~TeV can escape from distances beyond $0.1$~light cylinder radius ($R_{\rm{lc}}$). The emission surface from where photons can be emitted in the observer's direction further restricts the sites of VHE gamma-ray production. For the observation angle $57^{\circ}$ relative to the Crab pulsar axis of rotation and the orthogonal rotation, the emission surface in the open field line region is located as close as $0.4\,R_{\rm{lc}}$ from the stellar surface for a dipole magnetic field, and $0.1\,R_{\rm{lc}}$ for a force-free magnetic field.

The blazar Mrk 501 (z=0.034) was observed at very-high-energy (VHE, $E\gtrsim 100$~GeV) gamma-ray wavelengths during a bright flare on the night of 2014 June 23-24 (MJD 56832) with the H.E.S.S. phase-II array of Cherenkov telescopes. Data taken that night by H.E.S.S. at large zenith angle reveal an exceptional number of gamma-ray photons at multi-TeV energies, with rapid flux variability and an energy coverage extending significantly up to 20 TeV. This data set is used to constrain Lorentz invariance violation (LIV) using two independent channels: a temporal approach considers the possibility of an energy dependence in the arrival time of gamma rays, whereas a spectral approach considers the possibility of modifications to the interaction of VHE gamma rays with extragalactic background light (EBL) photons. The non-detection of energy-dependent time delays and the non-observation of deviations between the measured spectrum and that of a supposed power-law intrinsic spectrum with standard EBL attenuation are used independently to derive strong constraints on the energy scale of LIV ($E_{\rm{QG}}$) in the subluminal scenario for linear and quadratic perturbations in the dispersion relation of photons. For the case of linear perturbations, the 95% confidence level limits obtained are $E_{\rm{QG},1} > 3.6 \times 10^{17} \ \rm{GeV} $ using the temporal approach and $E_{\rm{QG},1} > 2.6 \times 10^{19} \ \rm{GeV}$ using the spectral approach. For the case of quadratic perturbations, the limits obtained are $E_{\rm{QG},2} > 8.5 \times 10^{10} \ \rm{GeV} $ using the temporal approach and $E_{\rm{QG},2} > 7.8 \times 10^{11} \rm{ GeV}$ using the spectral approach.

Aims: We present a detailed view of the pulsar wind nebula (PWN) HESS J1825-137. We aim to constrain the mechanisms dominating the particle transport within the nebula, accounting for its anomalously large size and spectral characteristics. Methods: The nebula is studied using a deep exposure from over 12 years of H.E.S.S. I operation, together with data from H.E.S.S. II improving the low energy sensitivity. Enhanced energy-dependent morphological and spatially-resolved spectral analyses probe the Very High Energy (VHE, E > 0.1 TeV) gamma-ray properties of the nebula. Results: The nebula emission is revealed to extend out to 1.5 degrees from the pulsar, ~1.5 times further than previously seen, making HESS J1825--137, with an intrinsic diameter of ~100 pc, potentially the largest gamma-ray PWN currently known. Characterisation of the nebula's strongly energy-dependent morphology enables the particle transport mechanisms to be constrained. A dependence of the nebula extent with energy of R $\propto$ E^\alpha with \alpha = -0.29 +/- 0.04 (stat) +/- 0.05 (sys) disfavours a pure diffusion scenario for particle transport within the nebula. The total gamma-ray flux of the nebula above 1~TeV is found to be (1.12 +/- 0.03 (stat) +/- 0.25 (sys)) $\times 10^{-11}$ cm$^{-2}$ s$^{-1}$, corresponding to ~64% of the flux of the Crab Nebula. Conclusions: HESS J1825-137 is a PWN with clear energy-dependent morphology at VHE gamma-ray energies. This source is used as a laboratory to investigate particle transport within middle-aged PWNe. Deep observations of this highly spatially-extended PWN enable a spectral map of the region to be produced, providing insights into the spectral variation within the nebula.

The total cosmic ray electron spectrum (electrons plus positrons) exhibits a break at a particle energy of $\sim 1\rm~TeV$ and extends without any attenuation up to $\rm \sim 20~ TeV $. Synchrotron and inverse Compton energy losses strongly constrain both the age and the distance of the potential sources of TeV and multi-TeV electrons to $\rm\approx 10^5~yr$ and $\rm \approx 100-500~pc$, depending on both the absolute value and energy dependence of the cosmic ray diffusion coefficient. This suggests that only a few, or just one nearby discrete source may explain the observed spectrum of high energy electrons. On the other hand the measured positron fraction, after initially increasing with particle energy, saturates at a level well below 0.5 and likely drops above $\sim 400-500$ GeV. This means that the local source(s) of TeV electrons should not produce positrons in equal amount, ruling out scenarios involving pulsars/pulsar winds as the main sources of high energy leptons. In this paper we show that a single, local, and fading source can naturally account for the entire spectrum of cosmic ray electrons in the TeV domain. Even though the nature of such source remains unclear, we discuss known cosmic ray accelerators, such as supernova remnant and stellar wind shocks, which are believed to accelerate preferentially electrons rather than positrons.

We report on NuSTAR observations of the mixed morphology supernova remnant (SNR) W49B, focusing on its nonthermal emission. Whereas radio observations as well as recent gamma-ray observations evidenced particle acceleration in this SNR, nonthermal X-ray emission has not been reported so far. With the unprecedented sensitivity of NuSTAR in the hard X-ray band, we detect a significant power-law-like component extending up to $\sim 20~{\rm keV}$, most probably of nonthermal origin. The newly discovered component has a photon index of $\Gamma =1.4^{+1.0}_{-1.1}$ with an energy flux between 10 and 20 keV of $(3.3 \pm 0.7) \times 10^{-13}~{\rm erg}~{\rm cm}^{-2}~{\rm s}^{-1}$. The emission mechanism is discussed based on the NuSTAR data combined with those in other wavelengths in the literature. The NuSTAR data, in terms both of the spectral slope and of the flux, are best interpreted as nonthermal electron bremsstrahlung. If this scenario is the case, then the NuSTAR emission provides a new probe to sub-relativistic particles accelerated in the SNR.

A hotspot at a position compatible with the BL Lac object 1ES 2322-409 was serendipitously detected with H.E.S.S. during observations performed in 2004 and 2006 on the blazar PKS 2316-423. Additional data on 1ES 2322-409 were taken in 2011 and 2012, leading to a total live-time of 22.3h. Point-like very-high-energy (VHE; E>100GeV) $\gamma$-ray emission is detected from a source centred on the 1ES 2322-409 position, with an excess of 116.7 events at a significance of 6.0$\sigma$. The average VHE $\gamma$-ray spectrum is well described with a power law with a photon index $\Gamma=3.40\pm0.66_{\text{stat}}\pm0.20_{\text{sys}}$ and an integral flux $\Phi(E>200GeV) = (3.11\pm0.71_{\rm stat}\pm0.62_{\rm sys})\times10^{-12} cm^{-2} s^{-1}$, which corresponds to 1.1$\%$ of the Crab nebula flux above 200 GeV. Multi-wavelength data obtained with Fermi LAT, Swift XRT and UVOT, RXTE PCA, ATOM, and additional data from WISE, GROND and Catalina, are also used to characterise the broad-band non-thermal emission of 1ES 2322-409. The multi-wavelength behaviour indicates day-scale variability. Swift UVOT and XRT data show strong variability at longer scales. A spectral energy distribution (SED) is built from contemporaneous observations obtained around a high state identified in Swift data. A modelling of the SED is performed with a stationary homogeneous one-zone synchrotron-self-Compton (SSC) leptonic model. The redshift of the source being unknown, two plausible values were tested for the modelling. A systematic scan of the model parameters space is performed, resulting in a well-constrained combination of values providing a good description of the broad-band behaviour of 1ES 2322-409.

Dwarf spheroidal galaxies are among the most promising targets for detecting signals of Dark Matter (DM) annihilations. The H.E.S.S. experiment has observed five of these systems for a total of about 130 hours. The data are re-analyzed here, and, in the absence of any detected signals, are interpreted in terms of limits on the DM annihilation cross section. Two scenarios are considered: i) DM annihilation into mono-energetic gamma-rays and ii) DM in the form of pure WIMP multiplets that, annihilating into all electroweak bosons, produce a distinctive gamma-ray spectral shape with a high-energy peak at the DM mass and a lower-energy continuum. For case i), upper limits at 95\% confidence level of about $\langle \sigma v \rangle \lesssim 3 \times 10^{-25}$ cm$^3$ s$^{-1}$ are obtained in the mass range of 400 GeV to 1 TeV. For case ii), the full spectral shape of the models is used and several excluded regions are identified, but the thermal masses of the candidates are not robustly ruled out.

We present the results from the Hitomi Soft Gamma-ray Detector (SGD) observation of the Crab nebula. The main part of SGD is a Compton camera, which in addition to being a spectrometer, is capable of measuring polarization of gamma-ray photons. The Crab nebula is one of the brightest X-ray / gamma-ray sources on the sky, and, the only source from which polarized X-ray photons have been detected. SGD observed the Crab nebula during the initial test observation phase of Hitomi. We performed the data analysis of the SGD observation, the SGD background estimation and the SGD Monte Carlo simulations, and, successfully detected polarized gamma-ray emission from the Crab nebula with only about 5 ks exposure time. The obtained polarization fraction of the phase-integrated Crab emission (sum of pulsar and nebula emissions) is (22.1 $\pm$ 10.6)% and, the polarization angle is 110.7$^o$ + 13.2 / $-$13.0$^o$ in the energy range of 60--160 keV (The errors correspond to the 1 sigma deviation). The confidence level of the polarization detection was 99.3%. The polarization angle measured by SGD is about one sigma deviation with the projected spin axis of the pulsar, 124.0$^o$ $\pm$0.1$^o$.

We report on the measurement and investigation of pulsed high-energy $\gamma$-ray emission from the Vela pulsar, PSR B0833$-$45, based on 40.3 hours of observations with the largest telescope of H.E.S.S., CT5, in monoscopic mode, and on 8 years of data obtained with the Fermi-LAT. A dedicated very-low-threshold event reconstruction and analysis pipeline was developed and, together with the CT5 telescope response model, was validated using the Fermi-LAT data as reference. A pulsed $\gamma$-ray signal at a significance level of more than $15\sigma$ is detected from the P2 peak of the Vela pulsar light curve. Of a total of 15835 events, more than 6000 lie at an energy below 20 GeV, implying a significant overlap between H.E.S.S. II-CT5 and the Fermi-LAT. While the investigation of the pulsar light curve with the LAT confirms characteristics previously known up to 20 GeV, in the tens of GeV energy range, CT5 data show a change in the pulse morphology of P2, i.e., an extreme sharpening of its trailing edge, together with the possible onset of a new component at 3.4$\sigma$ significance level. Assuming a power-law model for the P2 spectrum, an excellent agreement is found for the photon indices ($\Gamma \simeq$ 4.1) obtained with the two telescopes above 10 GeV and an upper bound of 8% is derived on the relative offset between their energy scales. Using both instruments data, it is however shown that the spectrum of P2 in the 10-100 GeV has a pronounced curvature, i.e. a confirmation of the sub-exponential cutoff form found at lower energies with the LAT. This is further supported by the weak evidence for an emission above 100 GeV obtained with CT5. In contrast, converging indications are found from both CT5 and LAT data for the emergence of a hard component above 50 GeV in the leading wing (LW2) of P2, which possibly extends beyond 100 GeV.

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.

(Abridged) Context. NGC 253 is one of only two starburst galaxies found to emit $\gamma$-rays from hundreds of MeV to multi-TeV energies. Accurate measurements of the very-high-energy (VHE) (E $>$ 100 GeV) and high-energy (HE) (E $>$ 60 MeV) spectra are crucial to study the underlying particle accelerators and cosmic-ray interaction and transport. Aims. The measurement of the VHE $\gamma$-ray emission of NGC 253 published in 2012 by H.E.S.S. was limited by large systematic uncertainties. Here, a measurement of the $\gamma$-ray spectrum of NGC 253 is investigated in both HE and VHE $\gamma$-rays. Methods. The data of H.E.S.S. observations are reanalysed using an updated calibration and analysis chain. The $Fermi$-LAT analysis employs more than 8 years of data processed using pass 8. The cosmic-ray particle population is evaluated from the combined HE--VHE $\gamma$-ray spectrum using NAIMA. Results. The VHE $\gamma$-ray energy spectrum is best fit by a power-law with a flux normalisation of $(1.34\,\pm\,0.14^{\mathrm{stat}}\,\pm\,0.27^{\mathrm{sys}}) \times 10^{-13} \mathrm{cm^{-2} s^{-1} TeV^{-1}}$ at 1 TeV -- about 40 \% above, but compatible with the value obtained in Abramowski et al. (2012). The spectral index $\Gamma = 2.39 \pm 0.14^{\mathrm{stat}} \pm 0.25^{\mathrm{sys}}$ is slightly softer than but consistent with the previous measurement. At energies above $\sim$3 GeV the HE spectrum is consistent with a power-law ranging into the VHE part of the spectrum measured by H.E.S.S. Conclusions. Two scenarios for the starburst nucleus are tested, in which the gas in the starburst nucleus acts as a target for hadronic cosmic rays. In these two models, the level to which NGC\,253 acts as a calorimeter is estimated to a range of $f_{\rm cal} = 0.1$ to $1$ while accounting for the measurement uncertainties.

The very high energy (VHE) gamma-ray emission reported from a number of pulsar wind nebulae (PWNe) is naturally explained by the inverse Compton scattering of multi-TeV electrons. However, the physical dimensions of some gamma-ray-emitting PWNe significantly exceed the scales anticipated by the standard hydrodynamical paradigm of PWN formation. The most "disturbing" case in this regard is HESS J1825-137, which extends to distances $r\approx70\rm\,pc$ from the central pulsar PSR J1826-1334. If the gamma-ray emission is indeed produced inside the PWN, but not by electrons that escaped the nebula and diffuse in the interstellar medium (ISM), the formation of such an anomalously extended plerion could be realized, in a diluted environment with the hydrogen number density $n_{ism}\le10^{-2}\rm\,cm^{-3}$. In this paper, we explore an alternative scenario assuming that the pulsar responsible for the formation of the nebula initially had a very short rotation period. In this case, the sizes of both the PWN and the surrounding supernova remnant depend on the initial pulsar period, the braking index, and the ISM density. To check the feasibility of this scenario, we study the parameter space that would reproduce the size of HESS J1825-137. We show that this demand can be achieved if the braking index is small, $n\leq2$ and the pulsar birth period is short, $P_{\rm b}\simeq1\rm\,ms$. This scenario can reproduce the wind termination position, which is expected at $R_{ts}\simeq0.03\rm\,pc$, only in a dense environment with $n_{ism}\geq\rm1\,cm^{-3}$. The requirement of the dense surrounding gas is supported by the presence of molecular clouds found in the source vicinity.

In a strong magnetic field, a high-energy photon can be absorbed and then produce an electron-positron pair. The produced electron/positron will in turn radiate a high-energy photon via synchrotron radiation, which then initiates a cascade. We built a one-dimensional Monte-Carlo code to study the development of the cascade especially after it reaches the saturated status, when almost all the energy of the primary particles transfers to the photons. The photon spectrum in this status has a cut-off due to the absorption by magnetic fields, which is much sharper than the exponential one. Below the cut-off, the spectral energy distribution (SED) manifest itself as a broken power-law with a spectral index of $0.5$ and $0.125$, respectively, below and above the broken energy. The SED can be fitted by a simple analytical function, which is solely determined by the product of the cascade scale $R$ and the magnetic field perpendicular to the motion of the particle B_{\perp}, with an accuracy better than 96\%. The similarity of the spectrum to that from the cascade in an isotropic black-body photon field is also studied.

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.

The "pion-decay" bump is a distinct signature of the differential energy spectrum of $\gamma$-rays between 100 MeV and 1 GeV produced in hadronic interactions of accelerated particles (cosmic rays) with the ambient gas. We use the recent parametrisations of relevant cross-sections to study the formation of the "pion-decay" bump. The $\gamma$-ray spectrum below the maximum of this spectral feature can be distorted because of contributions of additional radiation components, in particular, due to the bremsstrahlung of secondary electrons and positrons, the products of decays of $\pi^\pm$-mesons, accompanying the $\pi^0$-production. At energies below 100 MeV, a non-negligible fraction of $\gamma$-ray flux could originate from interactions of sub-relativistic heavy ions. We study the impact of these radiation channels on the formation of the overall $\gamma$-ray spectrum based on a time-dependent treatment of evolution of energy distributions of the primary and secondary particles in the $\gamma$-ray production region.

We discuss the discovery potential of extended very-high-energy (VHE) neutrino sources by the future KM3 Neutrino Telescope (KM3NeT) in the context of the constraining power of the Cherenkov Telescope Array (CTA), designed for deep surveys of the sky in VHE gamma rays. The study is based on a comparative analysis of sensitivities of KM3NeT and CTA. We show that a minimum gamma-ray energy flux of E^2{\phi}_{\gamma}(10 TeV) > 1x10^{-12} TeV cm^{-2} s^{-1} is required to identify a possible neutrino counterpart with a 3{\sigma} significance and 10 years of KM3NeT observations with upgoing muons, if the source has an angular size of R_{src} = 0.1 deg and emits gamma rays with an E^{-2} energy spectrum through a full hadronic mechanism. This minimum gamma-ray flux is increased to the level of E^2{\phi}_{\gamma}(10 TeV) > 2x10^{-11} TeV cm^{-2} s^{-1} in case of sources with radial extension of R_{src} = 2.0 deg. The analysis methods are applied to the supernova remnant RX J1713.7-3946 and the Galactic Center Ridge, as well as to the recent HAWC catalog of multi-TeV gamma-ray sources.

Recent advances in the $\gamma$-ray observations of solar flares by the \textit{Fermi} satellite, demand revisions in the hadronic $\gamma$-ray flux computation below 1 GeV. In this work we utilize recently updated pion production cross sections, along with an accurate description of low energy nuclear interactions. Applying these new interaction descriptions to model the \textit{Fermi} Large Area Telescope (LAT) solar flare data, we infer the primary particle spectral parameters. Application of this new cross section description leads to significantly different spectral parameters compared to those obtained previously. Furthermore, the inclusion of nuclei in these calculations leads to a primary spectrum that is generally harder than that required from proton only considerations. Lastly, the flare data at lower MeV energies, detected by the \textit{Fermi} Gamma-ray Burst Monitor (GBM), is demonstrated to provide additional low-energy spectral information.

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.

The X-ray emission from the jets in Active Galactic Nuclei (AGN) carries important information on the distributions of relativistic electrons and magnetic fields on large scales. We reanalyze archival Chandra observations on the jet of M87 from 2000 to 2016 with a total exposure of 1460 kiloseconds to explore the X-ray emission characteristics along the jet. We investigate the variability behaviours of the nucleus and the inner jet component HST-1, and confirm indications for day-scale X-ray variability in the nucleus contemporaneous to the 2010 high TeV gamma-ray state. HST-1 shows a general decline in X-ray flux over the last few years consistent with its synchrotron interpretation. We extract the X-ray spectra for the nucleus and all knots in the jet, showing that they are compatible with a single power-law within the X-ray band. There are indications of the resultant X-ray photon index to exhibit a trend, with slight but significant index variations ranging from $\simeq 2.2$ (e.g. in knot D) to $\simeq 2.4-2.6$ (in the outer knots F, A, and B). When viewed in a multi-wavelength context, a more complex situation is arising. Fitting the radio to X-ray spectral energy distributions (SEDs) assuming a synchrotron origin, we show that a broken power-law electron spectrum with break energy $E_b$ around $1~(300\mu G/B)^{1/2}$ TeV allows a satisfactorily description of the multi-band SEDs for most of the knots. However, in the case of knots B, C and D we find indications that an additional high energy component is needed to adequately reproduce the broadband SEDs. We discuss the implications and suggest that a stratified jet model may account for the differences.

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.

We present results from the Hitomi X-ray observation of a young composite-type supernova remnant (SNR) G21.5$-$0.9, whose emission is dominated by the pulsar wind nebula (PWN) contribution. The X-ray spectra in the 0.8-80 keV range obtained with the Soft X-ray Spectrometer (SXS), Soft X-ray Imager (SXI) and Hard X-ray Imager (HXI) show a significant break in the continuum as previously found with the NuSTAR observation. After taking into account all known emissions from the SNR other than the PWN itself, we find that the Hitomi spectra can be fitted with a broken power law with photon indices of $\Gamma_1=1.74\pm0.02$ and $\Gamma_2=2.14\pm0.01$ below and above the break at $7.1\pm0.3$ keV, which is significantly lower than the NuSTAR result ($\sim9.0$ keV). The spectral break cannot be reproduced by time-dependent particle injection one-zone spectral energy distribution models, which strongly indicates that a more complex emission model is needed, as suggested by recent theoretical models. We also search for narrow emission or absorption lines with the SXS, and perform a timing analysis of PSR J1833$-$1034 with the HXI and SGD. No significant pulsation is found from the pulsar. However, unexpectedly, narrow absorption line features are detected in the SXS data at 4.2345 keV and 9.296 keV with a significance of 3.65 $\sigma$. While the origin of these features is not understood, their mere detection opens up a new field of research and was only possible with the high resolution, sensitivity and ability to measure extended sources provided by an X-ray microcalorimeter.

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 present paper investigates the temperature structure of the X-ray emitting plasma in the core of the Perseus cluster using the 1.8--20.0 keV data obtained with the Soft X-ray Spectrometer (SXS) onboard the Hitomi Observatory. A series of four observations were carried out, with a total effective exposure time of 338 ks and covering a central region $\sim7'$ in diameter. The SXS was operated with an energy resolution of $\sim$5 eV (full width at half maximum) at 5.9 keV. Not only fine structures of K-shell lines in He-like ions but also transitions from higher principal quantum numbers are clearly resolved from Si through Fe. This enables us to perform temperature diagnostics using the line ratios of Si, S, Ar, Ca, and Fe, and to provide the first direct measurement of the excitation temperature and ionization temperature in the Perseus cluster. The observed spectrum is roughly reproduced by a single temperature thermal plasma model in collisional ionization equilibrium, but detailed line ratio diagnostics reveal slight deviations from this approximation. In particular, the data exhibit an apparent trend of increasing ionization temperature with increasing atomic mass, as well as small differences between the ionization and excitation temperatures for Fe, the only element for which both temperatures can be measured. The best-fit two-temperature models suggest a combination of 3 and 5 keV gas, which is consistent with the idea that the observed small deviations from a single temperature approximation are due to the effects of projection of the known radial temperature gradient in the cluster core along the line of sight. Comparison with the Chandra/ACIS and the XMM-Newton/RGS results on the other hand suggests that additional lower-temperature components are present in the ICM but not detectable by Hitomi SXS given its 1.8--20 keV energy band.

The Hitomi SXS spectrum of the Perseus cluster, with $\sim$5 eV resolution in the 2-9 keV band, offers an unprecedented benchmark of the atomic modeling and database for hot collisional plasmas. It reveals both successes and challenges of the current atomic codes. The latest versions of AtomDB/APEC (3.0.8), SPEX (3.03.00), and CHIANTI (8.0) all provide reasonable fits to the broad-band spectrum, and are in close agreement on best-fit temperature, emission measure, and abundances of a few elements such as Ni. For the Fe abundance, the APEC and SPEX measurements differ by 16%, which is 17 times higher than the statistical uncertainty. This is mostly attributed to the differences in adopted collisional excitation and dielectronic recombination rates of the strongest emission lines. We further investigate and compare the sensitivity of the derived physical parameters to the astrophysical source modeling and instrumental effects. The Hitomi results show that an accurate atomic code is as important as the astrophysical modeling and instrumental calibration aspects. Substantial updates of atomic databases and targeted laboratory measurements are needed to get the current codes ready for the data from the next Hitomi-level mission.

We present Hitomi observations of N132D, a young, X-ray bright, O-rich core-collapse supernova remnant in the Large Magellanic Cloud (LMC). Despite a very short observation of only 3.7 ks, the Soft X-ray Spectrometer (SXS) easily detects the line complexes of highly ionized S K and Fe K with 16-17 counts in each. The Fe feature is measured for the first time at high spectral resolution. Based on the plausible assumption that the Fe K emission is dominated by He-like ions, we find that the material responsible for this Fe emission is highly redshifted at ~800 km/s compared to the local LMC interstellar medium (ISM), with a 90% credible interval of 50-1500 km/s if a weakly informative prior is placed on possible line broadening. This indicates (1) that the Fe emission arises from the supernova ejecta, and (2) that these ejecta are highly asymmetric, since no blue-shifted component is found. The S K velocity is consistent with the local LMC ISM, and is likely from swept-up ISM material. These results are consistent with spatial mapping that shows the He-like Fe concentrated in the interior of the remnant and the S tracing the outer shell. The results also show that even with a very small number of counts, direct velocity measurements from Doppler-shifted lines detected in extended objects like supernova remnants are now possible. Thanks to the very low SXS background of ~1 event per spectral resolution element per 100 ks, such results are obtainable during short pointed or slew observations with similar instruments. This highlights the power of high-spectral-resolution imaging observations, and demonstrates the new window that has been opened with Hitomi and will be greatly widened with future missions such as the X-ray Astronomy Recovery Mission (XARM) and Athena.

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.

We report a Hitomi observation of IGR J16318-4848, a high-mass X-ray binary system with an extremely strong absorption of N_H~10^{24} cm^{-2}. Previous X-ray studies revealed that its spectrum is dominated by strong fluorescence lines of Fe as well as continuum emission. For physical and geometrical insight into the nature of the reprocessing material, we utilize the high spectroscopic resolving power of the X-ray microcalorimeter (the soft X-ray spectrometer; SXS) and the wide-band sensitivity by the soft and hard X-ray imager (SXI and HXI) aboard Hitomi. Even though photon counts are limited due to unintended off-axis pointing, the SXS spectrum resolves Fe K{\alpha_1} and K{\alpha_2} lines and puts strong constraints on the line centroid and width. The line width corresponds to the velocity of 160^{+300}_{-70} km s^{-1}. This represents the most accurate, and smallest, width measurement of this line made so far from any X-ray binary, much less than the Doppler broadening and shift expected from speeds which are characteristic of similar systems. Combined with the K-shell edge energy measured by the SXI and HXI spectra, the ionization state of Fe is estimated to be in the range of Fe I--IV. Considering the estimated ionization parameter and the distance between the X-ray source and the absorber, the density and thickness of the materials are estimated. The extraordinarily strong absorption and the absence of a Compton shoulder component is confirmed. These characteristics suggest reprocessing materials which are distributed in a narrow solid angle or scattering primarily with warm free electrons or neutral hydrogen.