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.

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.

The H.E.S.S. Galactic Plane Survey (HGPS) represents one of the most sensitive surveys of the Galactic Plane at very high energies (VHE, 0.1-100 TeV). However the source detection algorithm of the HGPS pipeline is not well-suited for complex regions, including sources with shell-like morphologies. As an alternative and complementary approach, we have investigated blind search methods for VHE gamma-ray source detection based on well-known and widely used image processing and pattern recognition techniques. Our goal is to build in a short amount of computational time a list of potentially valuable objects without prior case-specific morphological assumptions. We aim to classify and rank the detected objects in order to identify only the most promising source candidates for further multi-wavelength-association searches, dedicated analyses, or deeper observations. In the approach proposed, we extract sparse and pertinent structural information from the significance maps using a edge detection operator. We then apply a Hough circle transform and detect a collection of objects as local maxima in the Hough space. On the basis of morphological parameters we can characterize different object classes. Classification can be used to identify valuable source candidates sharing the characteristics of well-known sources. We show that using these pattern recognition techniques we can detect objects with partial circular symmetry irrespective of a morphological template (e.g. point-like, Gaussian-like, or shell-like). All the shell-type supernova remnants (SNRs) catalogued in the HGPS are associated with at least one detected object. Catalogue cross-matches indicate that several detected objects not catalogued in the HGPS are spatially coincident with multi-wavelength counterparts. Further investigation on the most promising candidates will be conducted in dedicated follow-up analyses.

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$.

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.

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.

Gamma-ray bursts (GRBs) are brief flashes of gamma rays, considered to be the most energetic explosive phenomena in the Universe. The emission from GRBs comprises a short (typically tens of seconds) and bright prompt emission, followed by a much longer afterglow phase. During the afterglow phase, the shocked outflow -- produced by the interaction between the ejected matter and the circumburst medium -- slows down, and a gradual decrease in brightness is observed. GRBs typically emit most of their energy via gamma-rays with energies in the kiloelectronvolt-to-megaelectronvolt range, but a few photons with energies of tens of gigaelectronvolts have been detected by space-based instruments. However, the origins of such high-energy (above one gigaelectronvolt) photons and the presence of very-high-energy (more than 100 gigaelectronvolts) emission have remained elussive. Here we report observations of very-high-energy emission in the bright GRB 180720B deep in the GRB afterglow -ten hours after the end of the prompt emission phase, when the X-ray flux had already decayed by four orders of magnitude. Two possible explanations exist for the observed radiation: inverse Compton emission and synchrotron emission of ultrarelativistic electrons. Our observations show that the energy fluxes in the X-ray and gamma-ray range and their photon indices remain comparable to each other throughout the afterglow. This discovery places distinct constraints on the GRB environment for both emission mechanisms, with the inverse Compton explanation alleviating the particle energy requirements for the emission observed at late times. The late timing of this detection has consequences for the future observations of GRBs at the highest energies.

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.

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.

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.

We have explored the capabilities of dust extinction and $\gamma$ rays to probe the properties of the interstellar medium in the nearby anti-centre region. We have jointly modelled the $\gamma$-ray intensity and the stellar reddening, E(B-V) as a combination of H$_{\rm I}$-bright, CO-bright, and ionised gas components. The complementary information from dust reddening and $\gamma$ rays is used to reveal the dark gas not seen, or poorly traced, by H$_{\rm I}$, free-free, and $^{12}$CO emissions. We compare the total gas column densities, $N_{\rm{H}}$, derived from the $\gamma$ rays and stellar reddening with those inferred from a similar analysis (Remy et al. 2017) of $\gamma$ rays and of the optical depth of the thermal dust emission, $\tau_{353}$, at 353 GHz. We can therefore compare environmental variations in specific dust reddening, E(B-V)/$N_{\rm H}$, and in dust emission opacity (dust optical depth per gas nucleon), $\tau_{353}/N_{\rm{H}}$. Over the whole anti-centre region, we find an average E(B-V)/$N_{\rm H}$ ratio of $(2.02\pm0.48)\times$ $10^{-22}$~mag~cm$^2$, with maximum local variations of about $\pm30\%$ at variance with the two to six fold coincident increase seen in emission opacity as the gas column density increases. In the diffuse medium, the small variations in specific reddening, E(B-V)/$N_{\rm H}$ implies a rather uniform dust-to-gas mass ratio in the diffuse parts of the anti-centre clouds. The small amplitude of the E(B-V)/$N_{\rm H}$ variations with increasing $N_{\rm{H}}$ column density confirms that the large opacity $\tau_{353}/N_{\rm{H}}$ rise seen toward dense CO clouds is primarily due to changes in dust emissivity. The environmental changes are qualitatively compatible with model predictions based on mantle accretion on the grains and the formation of grain aggregates.

H$_{\rm I}$ 21-cm and $^{12}$CO 2.6-mm line emissions respectively trace the atomic and molecular gas phases, but they miss most of the opaque H$_{\rm I}$ and diffuse H$_2$ present in the Dark Neutral Medium (DNM) at the transition between the H$_{\rm I}$-bright and CO-bright regions. Jointly probing H$_{\rm I}$, CO, and DNM gas, we aim to constrain the threshold of the H$_{\rm I}$-H$_2$ transition. We also aim to measure gas mass fractions in the different phases and to test their relation to cloud properties. We have used dust optical depth measurements at 353 GHz, $\gamma$-ray maps at GeV energies, and H$_{\rm I}$ and CO line data to trace the gas column densities and map the DNM in nearby clouds toward the Galactic anticentre and Chamaeleon regions. The onset of detectable CO intensities varies by only a factor of 4 from cloud to cloud, between 0.6$ \times 10^{21}$ cm$^{-2}$ and 2.5$ \times 10^{21}$ cm$^{-2}$ in total gas column density. We observe larger H$_2$ column densities than linearly inferred from the CO intensities at $A_{\rm{V}}> 3$ magnitudes because of the large CO optical thickness: the additional H$_2$ mass in this regime represents on average 20$\%$ of the CO-inferred molecular mass. In the DNM envelopes, we find that the fraction of diffuse CO-dark H$_2$ in the molecular column densities decreases with increasing $A_{\rm{V}}$ in a cloud. For a half molecular DNM, the fraction decreases from more than 80$\%$ at 0.4 magnitude to less than 20$\%$ beyond 2 magnitudes. In mass, the DNM fraction varies with the cloud properties. The mass present in the DNM envelopes appears to scale with the molecular mass seen in CO as $M_{\rm H}^{\rm DNM} = (62 \pm 7) {M_{\rm H_2}^{\rm CO}}^{\,0.51 \pm 0.02}$ across two decades in mass. The phase transitions in these clouds show both common trends and environmental differences.

We aim to explore the capabilities of dust emission and rays for probing the properties of the interstellar medium in the nearby anti-centre region, using gamma-ray observations with the Fermi Large Area Telescope (LAT), and the thermal dust optical depth inferred from Planck and IRAS observations. In particular, we aim at quantifying potential variations in cosmic-ray density and dust properties per gas nucleon across the different gas phases and different clouds, and at measuring the CO-to-H2 conversion factor, X$_{CO}$ , in different environments. We have separated six nearby anti-centre clouds that are coherent in velocities and distances, from the Galactic-disc background in HI 21-cm and $^{12}$CO 2.6-mm line emission. We have jointly modelled the gamma-ray intensity recorded between 0.4 and 100 GeV, and the dust optical depth at 353 GHz as a combination of HI-bright, CO-bright, and ionised gas components. The complementary information from dust emission and gamma rays was used to reveal the gas not seen, or poorly traced, by HI , free-free, and $^{12}$CO emissions, namely (i) the opaque HI and diffuse H$_2$ present in the Dark Neutral Medium at the atomic-molecular transition, and (ii) the dense H$_2$ to be added where $^{12}$CO lines saturate. The measured interstellar gamma-ray spectra support a uniform penetration of the cosmic rays with energies above a few GeV through the clouds. We find a gradual increase in grain opacity as the gas becomes more dense. The increase reaches a factor of four to six in the cold molecular regions that are well shielded from stellar radiation. Consequently, the X$_{CO}$ factor derived from dust is systematically larger by 30% to 130% than the gamma-ray estimate. We also evaluate the average gamma-ray X$_{CO}$ factorfor each cloud, and find that X$_{CO}$ tends to decrease from diffuse to more compact molecular clouds, as expected from theory.

High resolution stellar spectroscopic surveys provide massive amounts of diffuse interstellar bands (DIBs) measurements. Data can be used to study the distribution of the DIB carriers and those environmental conditions that favor their formation. In parallel, recent studies have also proved that DIBs extracted from stellar spectra constitute new tools for building the 3D structure of the Galactic Interstellar Medium (ISM). The amount of details on the structure depends directly on the quantity of available lines of sight (LOS). Therefore there is a need to construct databases of high-quality DIB measurements as large as possible. We aim at providing the community with a catalog of high-quality measurements of the 1.5273 micron DIB towards a large fraction of the Apache Point Observatory Galactic Evolution Experiment (APOGEE) hot stars observed to correct for the telluric absorption and not used for ISM studies so far. This catalog would complement the extensive database recently extracted from the APOGEE observations and used for 3D ISM mapping. We devised a method to fit the stellar continuum of the hot calibration stars and extracted the DIB from the normalized spectrum. Severe selection criteria based on the absorption characteristics are applied to the results. In particular limiting constraints on the DIB widths and Doppler shifts are deduced from the HI 21 cm measurements, following a new technique of decomposition of the emission spectra. From ~16 000 available hot telluric spectra we have extracted ~ 6700 DIB measurements and their associated uncertainties. The statistical properties of the extracted absorptions are examined and our selection criteria are shown to provide a robust dataset. The resulting catalog contains the DIB total equivalent widths, central wavelengths and widths. We briefly illustrate its potential use for the stellar and interstellar communities.

Most of the celestial gamma rays detected by the Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope originate from the interstellar medium when energetic cosmic rays interact with interstellar nucleons and photons. Conventional point and extended source studies rely on the modeling of this diffuse emission for accurate characterization. We describe here the development of the Galactic Interstellar Emission Model (GIEM) that is the standard adopted by the LAT Collaboration and is publicly available. The model is based on a linear combination of maps for interstellar gas column density in Galactocentric annuli and for the inverse Compton emission produced in the Galaxy. We also include in the GIEM large-scale structures like Loop I and the Fermi bubbles. The measured gas emissivity spectra confirm that the cosmic-ray proton density decreases with Galactocentric distance beyond 5 kpc from the Galactic Center. The measurements also suggest a softening of the proton spectrum with Galactocentric distance. We observe that the Fermi bubbles have boundaries with a shape similar to a catenary at latitudes below 20 degrees and we observe an enhanced emission toward their base extending in the North and South Galactic direction and located within 4 degrees of the Galactic Center.