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.

On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. Preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope (ENGRAVE) collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical/near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS and VINROUGE projects also carried out a search on this event. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN) possibly generated by this NS-BH merger, and for the strategy of future searches. Altogether, our observations allow us to exclude a KN with large ejecta mass $M\gtrsim 0.1\,\mathrm{M_\odot}$ to a high ($>90\%$) confidence, and we can exclude much smaller masses in a subsample of our observations. This disfavours the tidal disruption of the neutron star during the merger. Despite the sensitive instruments involved in the campaign, given the distance of S190814bv we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundreds Mpc will be detected only by large facilities with both high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event.

Long-duration gamma-ray bursts (GRBs) originate from ultra-relativistic jets launched from the collapsing cores of dying massive stars. They are characterised by an initial phase of bright and highly variable radiation in the keV-MeV band that is likely produced within the jet and lasts from milliseconds to minutes, known as the prompt emission. Subsequently, the interaction of the jet with the external medium generates external shock waves, responsible for the afterglow emission, which lasts from days to months, and occurs over a broad energy range, from the radio to the GeV bands. The afterglow emission is generally well explained as synchrotron radiation by electrons accelerated at the external shock. Recently, an intense, long-lasting emission between 0.2 and 1 TeV was observed from the GRB 190114C. Here we present the results of our multi-frequency observational campaign of GRB~190114C, and study the evolution in time of the GRB emission across 17 orders of magnitude in energy, from $5\times10^{-6}$ up to $10^{12}$\,eV. We find that the broadband spectral energy distribution is double-peaked, with the TeV emission constituting a distinct spectral component that has power comparable to the synchrotron component. This component is associated with the afterglow, and is satisfactorily explained by inverse Compton upscattering of synchrotron photons by high-energy electrons. We find that the conditions required to account for the observed TeV component are not atypical, supporting the possibility that inverse Compton emission is commonly produced in GRBs.

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

We report the results of optical follow-up observations of 29 gravitational-wave triggers during the first half of the LIGO-Virgo Collaboration (LVC) O3 run with the Gravitational-wave Optical Transient Observer (GOTO) in its prototype 4-telescope configuration (GOTO-4). While no viable electromagnetic counterpart candidate was identified, we estimate our 3D (volumetric) coverage using test light curves of on- and off-axis gamma-ray bursts and kilonovae. In cases where the source region was observable immediately, GOTO-4 was able to respond to a GW alert in less than a minute. The average time of first observation was $8.79$ hours after receiving an alert ($9.90$ hours after trigger). A mean of $732.3$ square degrees were tiled per event, representing on average $45.3$ per cent of the LVC probability map, or $70.3$ per cent of the observable probability. This coverage will further improve as the facility scales up alongside the localisation performance of the evolving gravitational-wave detector network. Even in its 4-telescope prototype configuration, GOTO is capable of detecting AT2017gfo-like kilonovae beyond 200 Mpc in favourable observing conditions. We cannot currently place meaningful electromagnetic limits on the population of distant ($\hat{D}_L = 1.3$ Gpc) binary black hole mergers because our test models are too faint to recover at this distance. However, as GOTO is upgraded towards its full 32-telescope, 2 node (La Palma & Australia) configuration, it is expected to be sufficiently sensitive to cover the predicted O4 binary neutron star merger volume, and will be able to respond to both northern and southern triggers.

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 report on the observations of gamma-ray burst (GRB) 190114C by the Fermi Gamma-ray Space Telescope and the Neil Gehrels Swift Observatory. The early-time observations reveal multiple emission components that evolve independently, with a delayed power-law component that exhibits significant spectral attenuation above 40 MeV in the first few seconds of the burst. This power-law component transitions to a harder spectrum that is consistent with the afterglow emission observed at later times. This afterglow component is clearly identifiable in the GBM and BAT light curves as a slowly fading emission component on which the rest of the prompt emission is superimposed. As a result, we are able to constrain the transition from internal shock to external shock dominated emission. We find that the temporal and spectral evolution of the broadband afterglow emission can be well modeled as synchrotron emission from a forward shock propagating into a wind-like circumstellar environment and find that high-energy photons observed by Fermi LAT are in tension with the theoretical maximum energy that can be achieved through synchrotron emission from a shock. These violations of the maximum synchrotron energy are further compounded by the detection of very high energy (VHE) emission above 300 GeV by MAGIC concurrent with our observations. We conclude that the observations of VHE photons from GRB 190114C necessitates either an additional emission mechanism at very high energies that is hidden in the synchrotron component in the LAT energy range, an acceleration mechanism that imparts energy to the particles at a rate that is faster than the electron synchrotron energy loss rate, or revisions of the fundamental assumptions used in estimating the maximum photon energy attainable through the synchrotron process.

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.

The composition of comets in the Solar System come in multiple groups thought to encode information about their formation in different regions of the outer protosolar disk. The recent discovery of the second interstellar object, 2I/Borisov, allows for spectroscopic investigations into its gas content and a preliminary classification of it within the Solar System comet taxonomies to test the applicability of planetesimal formation models to other stellar systems. We present spectroscopic and imaging observations from 2019 September 20th to October 26th at the Bok, MMT, and LBT telescopes. We identify CN in the comet's spectrum and set precise upper limits on the abundance of C2 on all dates. We use a Haser model to convert our integrated fluxes to production rates and find Q(CN) = 5.0 +/- 2.0 * 10^24 mol/s on September 20th and Q(CN) = 1.1 - 1.9 * 10^24 mol/s on later dates, both consistent with contemporaneous observations. We set our lowest upper limit on a C2 production rate, Q(C2) < 1.6 * 10^23 mol/s, on October 10th. The measured ratio upper limit for that date, Q(C2)/Q(CN) < 0.095 indicates that 2I/Borisov is strongly in the (carbon chain) 'depleted' taxonomic group. The only comparable Solar System comets have detected ratios near this limit, making 2I/Borisov statistically likely to be more depleted than any known comet. Most 'depleted' comets are Jupiter Family Comets, perhaps indicating a similiarity in formation conditions between the most depleted of the JFCs and 2I/Borisov. More work is needed to understand the applicability of our knowledge of Solar System comet taxonomies onto interstellar objects, and we discuss future work that could help clarify the usefulness of the approach.

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 Neil Gehrels Swift Observatory carried out prompt searches for gravitational wave (GW) events detected by the LIGO/Virgo Collaboration (LVC) during the second observing run ("O2"). Swift performed extensive tiling of eight LVC triggers, two of which had very low false-alarm rates (GW 170814 and the epochal GW 170817), indicating a high confidence of being astrophysical in origin; the latter was the first GW event to have an electromagnetic counterpart detected. In this paper we describe the follow-up performed during O2 and the results of our searches. No GW electromagnetic counterparts were detected; this result is expected, as GW 170817 remained the only astrophysical event containing at least one neutron star after LVC's later retraction of some events. A number of X-ray sources were detected, with the majority of identified sources being active galactic nuclei. We discuss the detection rate of transient X-ray sources and their implications in the O2 tiling searches. Finally, we describe the lessons learned during O2, and how these are being used to improve the \swift\ follow-up of GW events. In particular, we simulate a population of GRB afterglows to evaluate our source ranking system's ability to differentiate them from unrelated and uncatalogued X-ray sources. We find that $\approx$60-70% of afterglows whose jets are oriented towards Earth will be given high rank (i.e., "interesting" designation) by the completion of our second follow-up phase (assuming their location in the sky was observed), but that this fraction can be increased to nearly 100% by performing a third follow-up observation of sources exhibiting fading behavior.

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 report our identification of the optical afterglow and host galaxy of the short-duration gamma-ray burst GRB 160821B. The spectroscopic redshift of the host is $z=0.162$, making it one of the lowest redshift sGRBs identified by Swift. Our intensive follow-up campaign using a range of ground-based facilities as well as HST, XMM and Swift, shows evidence for a late-time excess of optical and near-infrared emission in addition to a complex afterglow. The afterglow light-curve at X-ray frequencies reveals a narrow jet, $\theta_j\sim1.9^{+0.10}_{-0.03}$ deg, that is refreshed at $>1$ day post-burst by a slower outflow with significantly more energy than the initial outflow that produced the main GRB. Observations of the 5 GHz radio afterglow shows a reverse shock into a mildly magnetised shell. The optical and near-infrared excess is fainter than AT2017gfo associated with GW170817, and is well explained by a kilonova with dynamic ejecta mass $M_{\rm dyn}=(1.0\pm0.6)\times10^{-3}$ M$_{\odot}$ and a secular (postmerger) ejecta mass with $M_{\rm pm}=(1.0\pm0.6)\times10^{-2}$ M$_\odot$, consistent with a binary neutron star merger resulting in a short-lived massive neutron star. This optical and near-infrared dataset provides the best-sampled kilonova light-curve without a gravitational wave trigger to date.

One of the central scientific goals of the next-generation Cherenkov Telescope Array (CTA) is the detection and characterization of gamma-ray bursts (GRBs). CTA will be sensitive to gamma rays with energies from about 20 GeV, up to a few hundred TeV. The energy range below 1 TeV is particularly important for GRBs. CTA will allow exploration of this regime with a ground-based gamma-ray facility with unprecedented sensitivity. As such, it will be able to probe radiation and particle acceleration mechanisms at work in GRBs. In this contribution, we describe POSyTIVE, the POpulation Synthesis Theory Integrated project for very high-energy emission. The purpose of the project is to make realistic predictions for the detection rates of GRBs with CTA, to enable studies of individual simulated GRBs, and to perform preparatory studies for time-resolved spectral analyses. The mock GRB population used by POSyTIVE is calibrated using the entire 40-year dataset of multi-wavelength GRB observations. As part of this project we explore theoretical models for prompt and afterglow emission of long and short GRBs, and predict the expected radiative output. Subsequent analyses are performed in order to simulate the observations with CTA, using the publicly available ctools and Gammapy frameworks. We present preliminary results of the design and implementation of this project.

We present X-ray and optical observations of the short duration gamma-ray burst GRB 071227 and its host at $z=0.381$, obtained using \textit{Swift}, Gemini South and the Very Large Telescope. We identify a short-lived and moderately bright optical transient, with flux significantly in excess of that expected from a simple extrapolation of the X-ray spectrum at 0.2-0.3 days after burst. We fit the SED with afterglow models allowing for high extinction and thermal emission models that approximate a kilonova to assess the excess' origins. While some kilonova contribution is plausible, it is not favoured due to the low temperature and high luminosity required, implying superluminal expansion and a large ejecta mass of $\sim 0.1$ M$_{\odot}$. We find, instead, that the transient is broadly consistent with power-law spectra with additional dust extinction of $E(B-V)\sim0.4$ mag, although a possibly thermal excess remains in the \textit{z}-band. We investigate the host, a spiral galaxy with an edge-on orientation, resolving its spectrum along its major axis to construct the galaxy rotation curve and analyse the star formation and chemical properties. The integrated host emission shows evidence for high extinction, consistent with the afterglow findings. The metallicity and extinction are consistent with previous studies of this host and indicate the galaxy is a typical, but dusty, late-type SGRB host.

The Transient High-Energy Sky and Early Universe Surveyor (THESEUS) is a mission concept developed in the last years by a large European consortium and currently under study by the European Space Agency (ESA) as one of the three candidates for next M5 mission (launch in 2032). THESEUS aims at exploiting high-redshift GRBs for getting unique clues to the early Universe and, being an unprecedentedly powerful machine for the detection, accurate location (down to arcsec) and redshift determination of all types of GRBs (long, short, high-z, under-luminous, ultra-long) and many other classes of transient sources and phenomena, at providing a substantial contribution to multi-messenger time-domain astrophysics. Under these respects, THESEUS will show a strong synergy with the large observing facilities of the future, like E-ELT, TMT, SKA, CTA, ATHENA, in the electromagnetic domain, as well as with next-generation gravitational-waves and neutrino detectors, thus greatly enhancing their scientific return.

The number of long gamma-ray bursts (GRBs) known to have occurred in the distant Universe (z greater than 5) is small (approx 15), however these events provide a powerful way of probing star formation at the onset of galaxy evolution. In this paper, we present the case for GRB100205A being a largely overlooked high-redshift event. While initially noted as a high-z candidate, this event and its host galaxy have not been explored in detail. By combining optical and near-infrared Gemini afterglow imaging (at t less than 1.3 days since burst) with deep late-time limits on host emission from the Hubble Space Telescope, we show that the most likely scenario is that GRB100205A arose in the redshift range 4-8. GRB100205A is an example of a burst whose afterglow, even at 1 hour post-burst, could only be identified by 8m class IR observations, and suggests that such observations of all optically dark bursts may be necessary to significantly enhance the number of high-redshift GRBs known.

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.

We report the X-ray data analysis of two transient ultraluminous X-ray sources (ULXs, hereafter X1 and X2) located in the nearby galaxy NGC 7090. While they were not detected in the 2004 XMM-Newton and 2005 Chandra observations, their 0.3-10 keV X-ray luminosities reached $>3\times10^{39}\,\mathrm{erg\,s^{-1}}$ in later XMM-Newton or Swift observations, showing increases in flux by a factor of $>80$ and $>300$ for X1 and X2, respectively. X1 showed indications of spectral variability: at the highest luminosity, its X-ray spectra can be fitted with a powerlaw ($\Gamma=1.55\pm0.15$), or a multicolour disc model with $T_{\mathrm{in}}=2.07^{+0.30}_{-0.23}$ keV; the X-ray spectrum became softer ($\Gamma=2.67^{+0.69}_{-0.64}$), or cooler ($T_\mathrm{in}=0.64^{+0.28}_{-0.17}$ keV) at lower luminosity. No strong evidence for spectral variability was found for X2. Its X-ray spectra can be fitted with a simple powerlaw model ($\Gamma=1.61^{+0.55}_{-0.50}$), or a multicolour disc model ($1.69^{+1.17}_{-0.48}$ keV). A possible optical counterpart for X1 is revealed in HST imaging. No optical variability is found, indicating that the optical radiation may be dominated by the companion star. Future X-ray and optical observations are necessary to determine the true nature of the compact object.

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 present a study of 21 dark gamma-ray burst (GRB) host galaxies, predominantly using X-ray afterglows obtained with the Chandra X-Ray Observatory (CXO) to precisely locate the burst in deep Hubble Space Telescope (HST) imaging of the burst region. The host galaxies are well-detected in F160W in all but one case and in F606W imaging in approx 60 per cent of cases. We measure magnitudes and perform a morphological analysis of each galaxy. The asymmetry, concentration and ellipticity of the dark burst hosts are compared against the host galaxies of optically bright GRBs. In agreement with other studies, we find that dark GRB hosts are redder and more luminous than the bulk of the GRB host population. The distribution of projected spatial offsets for dark GRBs from their host galaxy centroids is comparable to that of optically-bright bursts. The dark GRB hosts are physically larger, more massive and redder, but are morphologically similar to the hosts of bright GRBs in terms of concentration and asymmetry. Our analysis constrains the fraction of high redshift (z greater than 5) GRBs in the sample to approx 14 per cent, implying an upper limit for the whole long-GRB population of less than 4.4 per cent. If dust is the primary cause of afterglow darkening amongst dark GRBs, the measured extinction may require a clumpy dust component in order to explain the observed offset and ellipticity distributions.

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.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

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.

During its first observing run, in late 2015, the advanced LIGO facility announced 3 gravitational wave (GW) triggers to electromagnetic follow-up partners. Two of these have since been confirmed as being of astrophysical origin: both are binary black hole mergers at ~500 Mpc; the other trigger was later found not to be astrophysical. In this paper we report on the Swift follow up observations of the second and third triggers, including details of 21 X-ray sources detected; none of which can be associated with the GW event. We also consider the challenges that the next GW observing run will bring as the sensitivity and hence typical distance of GW events will increase. We discuss how to effectively use galaxy catalogues to prioritise areas for follow up, especially in the presence of distance estimates from the GW data. We also consider two galaxy catalogues and suggest that the high completeness at larger distances of the 2MASS Photometric Redshift Catalogue (2MPZ) makes it very well suited to optimise Swift follow-up observations.

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

In this paper we describe the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to accretion flows in the strong field gravity regime around both stellar-mass and supermassive black-holes. eXTP has the unique capability of using advanced 'spectral-timing-polarimetry' techniques to analyze the rapid variations with three orthogonal diagnostics of the flow and its geometry, yielding unprecedented insight into the inner accreting regions, the effects of strong field gravity on the material within them and the powerful outflows which are driven by the accretion process.

The short gamma-ray burst (GRB) 170817A was the first GRB associated with a gravitational-wave event. Due to the exceptionally low luminosity of the prompt $\gamma$-ray and the afterglow emission, the origin of both radiation components is highly debated. The most discussed models for the burst and the afterglow include a regular GRB jet seen off-axis and the emission from the cocoon encompassing a "choked" jet. Here, we report low radio-frequency observations at 610 and 1390~MHz obtained with the Giant Metrewave Radio Telescope (GMRT). Our observations span a range of $\sim7$ to $\sim152$ days after the burst. The afterglow started to emerge at these low frequencies about 60~days after the burst. The $1390$~MHz light curve barely evolved between 60 and 150 days, but its evolution is also marginally consistent with a $F_\nu\propto t^{0.8}$ rise seen in higher frequencies. We model the radio data and archival X-ray, optical and high-frequency radio data with models of top-hat and Gaussian structured GRB jets. We performed a Markov Chain Monte Carlo analysis of the structured-jet parameter space. Though highly degenerate, useful bounds on the posterior probability distributions can be obtained. Our bounds of the viewing angle are consistent with that inferred from the gravitational wave signal. We estimate the energy budget in prompt emission to be an order of magnitude lower than that in the afterglow blast-wave.

We present SPLOT, a small scale pilot survey to test the potential of snapshot (single epoch) linear imaging polarimetry as a supplementary tool to traditional transient follow-up. Transients exist in a vast volume of observational parameter space and polarimetry has the potential to highlight sources of scientific interest and add value to near real-time transient survey streams. We observed a sample of $\sim 50$ randomly selected optical transients with the EFOSC2 and SofI instruments, on the 3.6m New Technology Telescope (NTT) to test the feasibility of the survey. Our sample contained a number of interesting individual sources: a variety of supernovae, X-ray binaries, a tidal disruption event, blazar outbursts, and, by design, numerous transients of unknown nature. We discuss the results, both for the individual sources and the survey in detail. We provide an overview on the success and limitations of SPLOT and also describe a novel calibration method for removing instrumental polarisation effects from Nasymth-mounted telescopes. We find that a SPLOT-like survey would be a benefit to the large scale future transient survey streams such as LSST. The polarimetric measurements have added scientific value to a significant number of the sources and, most importantly, has shown the potential to highlight unclassified transient sources of scientific interest for further study.

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.

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.

Gamma-ray bursts (GRBs) occurring in the local Universe constitute an interesting sub-class of the GRB family, since their luminosity is on average lower than that of their cosmological analogs. We aim to contribute to the study of local bursts by reporting the case of GRB 171205A. This source was discovered by Swift Burst Alert Telescope (BAT) on 2017, December 5 and soon associated with a low redshift host galaxy (z=0.037), and an emerging SN (SN 2017iuk). We analyzed the full Swift, dataset, comprising the UV-Optical Telescope (UVOT), X-ray Telescope (XRT) and BAT data. In addition, we employed the Konus-Wind high energy data as a valuable extension at gamma-ray energies. The photometric SN signature is clearly visible in the UVOT u, b and v filters. The maximum emission is reached at ~ 13 (rest frame) days, and the whole bump resembles that of SN 2006aj, but lower in magnitude and with a shift in time of +2 d. A prebump in the v-band is also clearly visible, and this is the first time that such a feature is not observed achromatically in GRB-SNe. Its physical origin cannot be easily explained. The X-ray spectrum shows an intrinsic Hydrogen column density N_H,int = 7.4(+4.1 -3.6) X 10^20 / cm^2$, which is at the low end of the N_H,int, even considering just low redshift GRBs. The spectrum also features a thermal component, which is quite common in GRBs associated with SNe, but whose origin is still a matter of debate. Finally, the isotropic energy in the gamma-ray band, E_iso = 2.18(+0.63 -0.50) X 10^49 erg, is lower than those of cosmological GRBs. Combining this value with the peak energy in the same band, E_p=125(+141 -37) keV, implies that GRB 171205A is an outlier of the Amati relation, as are some other low redshift GRBs, and its emission mechanism should be different from that of canonical, farther away GRBs.

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.

Simultaneous XMM-Newton, NuSTAR and HST observations, performed in March 2017, of the nearby ($z=0.184$) luminous quasar PDS 456 are presented. PDS 456 had a low X-ray flux compared to past observations, where the first of the two new XMM-Newton observations occurred during a pronounced dip in the X-ray lightcurve. The broad-band X-ray spectrum is highly absorbed, attenuated by a soft X-ray absorber of column density $N_{\rm H}=6\times10^{22}$ cm$^{-2}$. An increase in obscuration occurs during the dip, which may be due to an X-ray eclipse. In addition, the persistent, fast Fe K outflow is present, with velocity components of $-0.25c$ and $-0.4c$. The soft absorber is less ionized ($\log\xi=3$) compared to the iron K outflow ($\log\xi=5$) and is outflowing with a velocity of approximately $-0.2c$. A soft X-ray excess is present below 1 keV against the highly absorbed continuum and can be attributed to the re-emission from a wide angle wind. The complex X-ray absorption present in PDS 456 suggests that the wind is inhomogeneous, whereby the soft X-ray absorber originates from denser clumps or filaments which may form further out along the outflow. In contrast to the X-ray observations, the simultaneous UV spectrum of PDS 456 is largely unabsorbed, where only a very weak broad absorption trough is present bluewards of Ly$\alpha$, compared to a past observation in 2000 when the trough was significantly stronger. The relative weakness of the UV absorption may be due to the soft X-ray absorber being too highly ionized and almost transparent in the UV band.

The binary neutron star merger GW170817 was the first multi-messenger event observed in both gravitational and electromagnetic waves. The electromagnetic signal began approximately 2 seconds post-merger with a weak, short burst of gamma-rays, which was followed over the next hours and days by the ultraviolet, optical and near-infrared emission from a radioactively- powered kilonova. Later, non-thermal rising X-ray and radio emission was observed. The low luminosity of the gamma-rays and the rising non-thermal flux from the source at late times could indicate that we are outside the opening angle of the beamed relativistic jet. Alternatively, the emission could be arising from a cocoon of material formed from the interaction between a jet and the merger ejecta. Here we present late-time optical detections and deep near-infrared limits on the emission from GW170817 at 110 days post-merger. Our new observations are at odds with expectations of late-time emission from kilonova models, being too bright and blue. Instead, the emission arises from the interaction between the relativistic ejecta of GW170817 and the interstellar medium. We show that this emission matches the expectations of a Gaussian structured relativistic jet, which would have launched a high luminosity short GRB to an aligned observer. However, other jet structure or cocoon models can also match current data - the future evolution of the afterglow will directly distinguish the origin of the emission.

Asteroid (3200) Phaethon is a Near-Earth Apollo asteroid with an unusual orbit that brings it closer to the Sun than any other known asteroid. Its last close approach to the Earth was in mid-December 2017 and the next one will be on October 2026. Previous rotationally time-resolved spectroscopy of Phaethon showed that its spectral slope is slightly bluish, in agreement with its B/F taxonomic classification, but at some rotational phases, it changes to slightly reddish. Motivated by this result we performed time-resolved imaging polarimetry of Phaethon during its recent close approach to the Earth. Phaethon has a spin period of 3.604 hours and we found a variation of the linear polarisation with rotation. This seems to be a rare case in which such variation is unambiguously found, also a consequence of its fairly large amplitude. Combining this new information with the brightness and colour variation, as well as previously reported results from Arecibo radar observations, we conclude that there is no variation of the mineralogy across the surface of Phaeton. However, the observed change in the linear polarisation may be related to differences in the thickness of the surface regolith in different areas or local topographic features.

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.

Aim. The unusual transient ATLAS17aeu was serendipitously detected within the sky localisation of the gravitational wave trigger GW170104. The importance of a possible association with gravitational waves coming from a binary black hole merger led to an extensive follow-up campaign, with the aim of assessing a possible connection with GW170104. Methods. With several telescopes, we carried out both photometric and spectroscopic observations of ATLAS17aeu, for several epochs, between $\sim 3$ and $\sim 230$ days after the first detection. Results. We studied in detail the temporal and spectroscopic properties of ATLAS17aeu and its host galaxy. We detected spectral features similar to those of a broad lined supernova superposed to an otherwise typical long-GRB afterglow. Based on analysis of the optical light curve, spectrum and host galaxy SED, we conclude that the redshift of the source is probably $z \simeq 0.5 \pm 0.2$. Conclusions. While the redshift range we have determined is marginally compatible with that of the gravitational wave event, the presence of a supernova component and the consistency of this transient with the E$_{\rm p}$-E$_{\rm iso}$ correlation support the conclusion that ATLAS17aeu was associated with the long gamma-ray burst GRB170105A. This rules out the association of the GRB170105A/ATLAS17aeu transient with the gravitational wave event GW170104, which was due to a binary black hole merger.

The historic first joint detection of both gravitational wave and electromagnetic emission from a binary neutron star merger cemented the association between short gamma-ray bursts (SGRBs) and compact object mergers, as well as providing a well sampled multi-wavelength light curve of a radioactive kilonova (KN) for the first time. Here we compare the optical and near-infrared light curves of this KN, AT2017gfo, to the counterparts of a sample of nearby (z < 0.5) SGRBs to characterize their diversity in terms of their brightness distribution. Although at similar epochs AT2017gfo appears fainter than every SGRB-associated KN claimed so far, we find three bursts (GRBs 050509B, 061201 and 080905A) where, if the reported redshifts are correct, deep upper limits rule out the presence of a KN similar to AT2017gfo by several magnitudes. Combined with the properties of previously claimed KNe in SGRBs this suggests considerable diversity in the properties of KN drawn from compact object mergers, despite the similar physical conditions that are expected in many NS-NS mergers. We find that observer angle alone is not able to explain this diversity, which is likely a product of the merger type (NS-NS versus NS-BH) and the detailed properties of the binary (mass ratio, spins etc). Ultimately disentangling these properties should be possible through observations of SGRBs and gravitational wave sources, providing direct measurements of heavy element enrichment throughout the Universe.

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

Flares in the X-ray afterglow of gamma-ray bursts (GRBs) share more characteristics with the prompt emission than the afterglow, such as pulse profile and contained fluence. As a result, they are believed to originate from late-time activity of the central engine and can be used to constrain the overall energy budget. In this paper, we collect a sample of $19$ long GRBs observed by \emph{Swift}-XRT that contain giant flares in their X-ray afterglows. We fit this sample with a version of the magnetar propeller model, modified to include fallback accretion. This model has already successfully reproduced extended emission in short GRBs. Our best fits provide a reasonable morphological match to the light curves. However, $16$ out of $19$ of the fits require efficiencies for the propeller mechanism that approach $100\%$. The high efficiency parameters are a direct result of the high energy contained in the flares and the extreme duration of the dipole component, which forces either slow spin periods or low magnetic fields. We find that even with the inclusion of significant fallback accretion, in all but a few cases it is energetically challenging to produce prompt emission, afterglow and giant flares within the constraints of the rotational energy budget of a magnetar.

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.

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift $\sim$10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA).

The recent discovery of the electromagnetic counterpart of the gravitational wave source GW170817, has demonstrated the huge informative power of multi-messenger observations. During the next decade the nascent field of multi-messenger astronomy will mature significantly. Around 2030, third generation gravitational wave detectors will be roughly ten times more sensitive than the current ones. At the same time, neutrino detectors currently upgrading to multi km^3 telescopes, will include a 10 km^3 facility in the Southern hemisphere that is expected to be operational around 2030. In this review, we describe the most promising high frequency gravitational wave and neutrino sources that will be detected in the next two decades. In this context, we show the important role of the Transient High Energy Sky and Early Universe Surveyor (THESEUS), a mission concept proposed to ESA by a large international collaboration in response to the call for the Cosmic Vision Programme M5 missions. THESEUS aims at providing a substantial advancement in early Universe science as well as playing a fundamental role in multi-messenger and time-domain astrophysics, operating in strong synergy with future gravitational wave and neutrino detectors as well as major ground- and space-based telescopes. This review is an extension of the THESEUS white paper (Amati et al. 2017), also in light of the discovery of GW170817/GRB170817A that was announced on October 16th, 2017.

The observation of the electromagnetic counterpart of gravitational-wave (GW) transient GW170817 demonstrated the potential in extracting astrophysical information from multimessenger discoveries. The forthcoming deployment of the first telescopes of the Cherenkov Telescope Array (CTA) observatory will coincide with Advanced LIGO/Virgo's next observing run, O3, enabling the monitoring of gamma-ray emission at E > 20 GeV, and thus particle acceleration, from GW sources. CTA will not be greatly limited by the precision of GW localization as it will be be capable of rapidly covering the GW error region with sufficient sensitivity. We examine the current status of GW searches and their follow-up effort, as well as the status of CTA, in order to identify some of the general strategies that will enhance CTA's contribution to multimessenger discoveries.

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.

Recent advances in fitting prompt emission spectra in gamma-ray bursts (GRBs) are boosting our understanding of the still elusive origin of this radiation. These progresses have been possible thanks to a more detailed analysis of the low-energy part ($<$\,100\,keV) of the prompt spectrum, where the spectral shape is sometimes found to deviate from a simple power-law shape. This deviation is well described by a spectral break or, alternatively by the addition of a thermal component. Spectral data extending down to less than 1\,keV are extremely relevant for these studies, but presently they are available only for a small subsample of {\it Swift} GRBs observed by XRT (the X-ray telescope, 0.3-10\,keV) during the prompt emission. The space mission \th\ will allow a systematic study of prompt spectra from 0.3\,keV to several MeV. We show that observations performed by \th\ will allow us to discriminate between different models presently considered for GRB prompt studies, solving the long-standing open issue about the nature of the prompt radiation, with relevant consequences on the location of the emitting region, magnetic field strength and presence of thermal components.

The Infra-Red Telescope (IRT) on board the Transient High Energy Sky and Early Universe Surveyor (THESEUS) ESA M5 candidate mission will play a key role in identifying and characterizing moderate to high redshift Gamma-Ray Bursts afterglows. The IRT is the enabling instrument on board THESEUS for measuring autonomously the redshift of the several hundreds of GRBs detected per year by the Soft X-ray Imager (SXI) and the X- and Gamma-Ray Imaging Spectrometer (XGIS), and thus allowing the big ground based telescopes to be triggered on a redshift pre-selected sample, and finally fulfilling the cosmological goals of the mission. The IRT will be composed by a primary mirror of 0.7 m of diameter coupled to a single camera in a Cassegrain design. It will work in the 0.7-1.8 {\mu}m wavelength range, and will provide a 10x10 arc min imaging field of view with sub-arc second localization capabilities, and, at the same time, a 5x5 arc min field of view with moderate (R up to ~500) spectroscopic capabilities. Its sensitivity, mainly limited by the satellite jitter, is adapted to detect all the GRBs, localized by the SXI/XGIS, and to acquire spectra for the majority of them.

We summarize in this contribution the capabilities, design status, and the en- abling technologies of the Soft X-ray Imager (SXI) planned to be on-board the THESEUS mission. We describe its central role in making THESEUS a powerful machine to probe the physical conditions of the early Universe (close to the reionization era) and to explore the time-domain Universe.

The Transient High-Energy Sky and Early Universe Surveyor (THESEUS) is a mission concept developed in the last years by a large European consortium, with interest in prospective participation by research groups in USA and other non-European countries. As detailed in Amati et al. 2017 (arXiv:1710.04638) and Stratta et al. 2017 (arXiv:1712.08153), THESEUS aims at exploiting high-redshift Gamma-Ray Bursts for getting unique clues on the early Universe and, being an unprecedentedly powerful machine for the detection, accurate location and redshift determination of all types of GRBs (long, short, high-z, under-luminous, ultra-long) and many other classes of transient sources and phenomena, at providing a substantial contribution to multi-messenger astrophysics and time-domain astronomy. Under these respects, THESEUS will show a beautiful synergy with the large observing facilities of the future, like E-ELT, TMT, SKA, CTA, ATHENA, in the electromagnetic domain, as well as with next-generation gravitational-waves and neutrino detectors, thus enhancing importantly their scientific return. Moreover, it will also operate as a flexible IR and X-ray observatory, thus providing an even larger involvement of the scientific community, as is currently the case for the Swift mission. In order to further explore the magnificent prospective science of the mission, the THESEUS consortium organized a Workshop in Naples on October 5-6 2017. The programme included about 50 reviews and talks from worldwide recognized experts of the fields. The topics ranged from the description of the mission concept, instrumentation and technologies to the main, additional and observatory science, further showing the strong impact that THESEUS observations would have on several fields of astrophysics, cosmology and fundamental physics.

We will discuss the observing strategy of the Transient High Energy Sky and Early Universe Surveyor (THESEUS) mission proposed to ESA as a response to the M5 call for proposals. The description of THESEUS and its science goals can be found in the white paper by Amati et al. (2017).