The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was equipped with an interim camera (JPAS-Pathfinder), composed of one CCD with a 0.3deg$^2$ field-of-view and resolution of 0.23 arcsec pixel$^{-1}$. To demonstrate the scientific potential of J-PAS, with the JPAS-Pathfinder camera we carried out a survey on the AEGIS field (along the Extended Groth Strip), dubbed miniJPAS. We observed a total of $\sim 1$ deg$^2$, with the $56$ J-PAS filters, which include $54$ narrow band (NB, $\rm{FWHM} \sim 145$Angstrom) and two broader filters extending to the UV and the near-infrared, complemented by the $u,g,r,i$ SDSS broad band (BB) filters. In this paper we present the miniJPAS data set, the details of the catalogues and data access, and illustrate the scientific potential of our multi-band data. The data surpass the target depths originally planned for J-PAS, reaching $\rm{mag}_{\rm {AB}}$ between $\sim 22$ and $23.5$ for the NB filters and up to $24$ for the BB filters ($5\sigma$ in a $3$~arcsec aperture). The miniJPAS primary catalogue contains more than $64,000$ sources extracted in the $r$ detection band with forced photometry in all other bands. We estimate the catalogue to be complete up to $r=23.6$ for point-like sources and up to $r=22.7$ for extended sources. Photometric redshifts reach subpercent precision for all sources up to $r=22.5$, and a precision of $\sim 0.3$% for about half of the sample. (Abridged)

From the approximately $\sim$3,500 planetary nebulae (PNe) discovered in our Galaxy, only 14 are known to be members of the Galactic halo. Nevertheless, a systematic search for halo PNe has never been performed. In this study, we present new photometric diagnostic tools to identify compact PNe in the Galactic halo by making use of the novel 12-filter system projects, J-PLUS (Javalambre Photometric Local Universe Survey) and S-PLUS (Southern-Photometric Local Universe Survey). We reconstructed the IPHAS (Isaac Newton Telescope (INT) Photometric H${\alpha}$ Survey of the Northern Galactic Plane) diagnostic diagram and propose four new ones using i) the J-PLUS and S-PLUS synthetic photometry for a grid of photo-ionisation models of halo PNe, ii) several observed halo PNe, as well as iii) a number of other emission-line objects that resemble PNe. All colour-colour diagnostic diagrams are validated using two known halo PNe observed by J-PLUS during the scientific verification phase and the first data release (DR1) of S-PLUS and the DR1 of J-PLUS. By applying our criteria to the DR1s ($\sim$1,190 deg$^2$), we identified one PN candidate. However, optical follow-up spectroscopy proved it to be a H II region belonging to the UGC 5272 galaxy. Here, we also discuss the PN and two H II galaxies recovered by these selection criteria. Finally, the cross-matching with the most updated PNe catalogue (HASH) helped us to highlight the potential of these surveys, since we recover all the known PNe in the observed area. The tools here proposed to identify PNe and separate them from their emission-line contaminants proved to be very efficient thanks to the combination of many colours, even when applied -like in the present work- to an automatic photometric search that is limited to compact PNe.

The Southern Photometric Local Universe Survey (S-PLUS) is imaging ~9300 deg^2 of the celestial sphere in twelve optical bands using a dedicated 0.8 m robotic telescope, the T80-South, at the Cerro Tololo Inter-American Observatory, Chile. The telescope is equipped with a 9.2k by 9.2k e2v detector with 10 um pixels, resulting in a field-of-view of 2 deg^2 with a plate scale of 0.55"/pixel. The survey consists of four main subfields, which include two non-contiguous fields at high Galactic latitudes (8000 deg^2 at |b| > 30 deg) and two areas of the Galactic plane and bulge (for an additional 1300 deg^2). S-PLUS uses the Javalambre 12-band magnitude system, which includes the 5 u, g, r, i, z broad-band filters and 7 narrow-band filters centered on prominent stellar spectral features: the Balmer jump/[OII], Ca H+K, H-delta, G-band, Mg b triplet, H-alpha, and the Ca triplet. S-PLUS delivers accurate photometric redshifts (delta_z/(1+z) = 0.02 or better) for galaxies with r < 20 AB mag and redshift < 0.5, thus producing a 3D map of the local Universe over a volume of more than 1 (Gpc/h)^3. The final S-PLUS catalogue will also enable the study of star formation and stellar populations in and around the Milky Way and nearby galaxies, as well as searches for quasars, variable sources, and low-metallicity stars. In this paper we introduce the main characteristics of the survey, illustrated with science verification data highlighting the unique capabilities of S-PLUS. We also present the first public data release of ~336 deg^2 of the Stripe-82 area, which is available at http://datalab.noao.edu/splus.

We investigated 12 unclassified B[e] stars or candidates, 8 from the Galaxy, 2 from the Large Magellanic Cloud (LMC) and 2 from the Small Magellanic Cloud (SMC). Based on the analysis of high-resolution spectroscopic (FEROS) and photometric data, we confirmed the presence of the B[e] phenomenon for all objects of our sample, except for one (IRAS 07455-3143). We derived their effective temperature, spectral type, luminosity class, interstellar extinction and, using the distances from Gaia DR2, we obtained their bolometric magnitude, luminosity and radius. Modeling of the forbidden lines present in the FEROS spectra revealed information about the kinematics and geometry of the circumstellar medium of these objects. In addition, we analyzed the light curves of four stars, finding their most probable periods. The evolutionary stage of 11 stars of our sample is suggested from their position on the HR diagram, taking into account evolutionary tracks of stars with solar, LMC and SMC metallicities. As results, we identified B and B[e] supergiants, B[e] stars probably at the main sequence or close to its end, post-AGB and HAeB[e] candidates, and A[e] stars in the main sequence or in the pre-main sequence. However, our most remarkable results are the identification of the third A[e] supergiant (ARDB\,54, the first one in the LMC), and of an "LBV impostor" in the SMC (LHA 115-N82).

The tidal interaction of a Be star with a binary companion forms two spiral arms that cause orbital modulation of the Be disc structure. The aim of this work is to identify observables in which this modulation is apparent. The structure of a Be disc in a coplanar circular binary system is computed with a smoothed-particle hydrodynamics code, and a radiation transfer code calculates the spectral energy distribution. Line depolarisation was confirmed, with polarisation profiles nearly reverse to emission-line profiles. The continuum flux maximizes for pole-on discs, but photometric variability maximizes for edge-on discs. The linear polarisation exhibits one or two maxima per orbital cycle. While polarisation variability in visible passbands is important only at low inclinations, infrared bands may demonstrate high orbital variability even at large inclinations. More evident is the modulation in the polarisation angle (PA) for low inclinations. The latter can be used to track azimuthal asymmetries for pole-on discs, where the spectroscopic variability in the violet-to-red (V/R) emission-component ratio disappears. PA reversals coincide with phases where V/R=1, tracking lines of sight directed towards regions where the approaching and receding arms overlap. Continuum flux and polarisation are mostly in phase for neighbouring wavelength regions. It is suggested that studies of non-symmetric discs distorted by tidal forces from a secondary star may be used to study disc variabilities of other origins.

Using a sample of 31 main-sequence OB stars located between galactocentric distances 8.4 - 15.6 kpc, we aim to probe the present-day radial abundance gradients of the Galactic disk. The analysis is based on high-resolution spectra obtained with the MIKE spectrograph on the Magellan Clay 6.5-m telescope on Las Campanas. We used a non-NLTE analysis in a self-consistent semi-automatic routine based on TLUSTY and SYNSPEC to determine atmospheric parameters and chemical abundances. Stellar parameters (effective temperature, surface gravity, projected rotational velocity, microturbulence, and macroturbulence) and silicon and oxygen abundances are presented for 28 stars located beyond 9 kpc from the Galactic centre plus three stars in the solar neighborhood. The stars of our sample are mostly on the main-sequence, with effective temperatures between 20800 - 31300 K, and surface gravities between 3.23 - 4.45 dex. The radial oxygen and silicon abundance gradients are negative and have slopes of -0.07 dex/kpc and -0.09 dex/kpc, respectively, in the region $8.4 \leq R_G \leq 15.6$\,kpc. The obtained gradients are compatible with the present-day oxygen and silicon abundances measured in the solar neighborhood and are consistent with radial metallicity gradients predicted by chemodynamical models of Galaxy Evolution for a subsample of young stars located close to the Galactic plane.

We present a new methodology for the estimation of stellar atmospheric parameters from narrow- and intermediate-band photometry of the Javalambre Photometric Local Universe Survey (J-PLUS), and propose a method for target pre-selection of low-metallicity stars for follow-up spectroscopic studies. Photometric metallicity estimates for stars in the globular cluster M15 are determined using this method. By development of a neural-network-based photometry pipeline, we aim to produce estimates of effective temperature, $T_{\rm eff}$, and metallicity, [Fe/H], for a large subset of stars in the J-PLUS footprint. The Stellar Photometric Index Network Explorer, SPHINX, is developed to produce estimates of $T_{\rm eff}$ and [Fe/H], after training on a combination of J-PLUS photometric inputs and synthetic magnitudes computed for medium-resolution (R ~ 2000) spectra of the Sloan Digital Sky Survey. This methodology is applied to J-PLUS photometry of the globular cluster M15. Effective temperature estimates made with J-PLUS Early Data Release photometry exhibit low scatter, \sigma($T_{\rm eff}$) = 91 K, over the temperature range 4500 < $T_{\rm eff}$ (K) < 8500. For stars from the J-PLUS First Data Release with 4500 < $T_{\rm eff}$ (K) < 6200, 85 $\pm$ 3% of stars known to have [Fe/H] <-2.0 are recovered by SPHINX. A mean metallicity of [Fe/H]=-2.32 $\pm$ 0.01, with a residual spread of 0.3 dex, is determined for M15 using J-PLUS photometry of 664 likely cluster members. We confirm the performance of SPHINX within the ranges specified, and verify its utility as a stand-alone tool for photometric estimation of effective temperature and metallicity, and for pre-selection of metal-poor spectroscopic targets.

B[e] Supergiants are luminous evolved massive stars. The mass-loss during this phase creates a complex circumstellar environment with atomic, molecular, and dusty regions usually found in rings or disk-like structures. For a better comprehension of the mechanisms behind the formation of these rings, detailed knowledge about their structure and dynamics is essential. To address that, we obtained high-resolution optical and near-infrared spectra for 8 selected Galactic B[e] Supergiants, for which CO emission has been detected. Assuming Keplerian rotation for the disk, we combine the kinematics obtained from the CO bands in the near-IR with those obtained by fitting the forbidden emission [OI] $\lambda$5577, [OI] $\lambda\lambda$6300,6363, and [CaII] $\lambda\lambda$7291,7323 lines in the optical to probe the disk structure. We find that the emission originates from multiple ring structures around all B[e] Supergiants, with each one of them displaying a unique combination of rings regardless of whether the object is part of a binary system. The confirmed binaries display spectroscopic variations of their line intensities and profiles as well as photometric variability, whereas the ring structures around the single stars are stable.

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.

Prompted by peculiar spectroscopic variability observed in SDSS/APOGEE $H$-band spectra, we monitored the Be star HD 55606 using optical spectroscopy and found that it is an exotic double-lined spectroscopic binary (SB2) consisting of a Be star and a hot, compact companion that is probably an OB subdwarf (sdOB) star. Motion of the sdOB star is traced by its impact on the strong He~I lines, observed as radial velocity ($V_{r}$) variable, double-peaked emission profiles with narrow central absorption cores. Weak He II 4686 {\AA} absorption associated with the companion star is detected in most spectra. Use of the emission peaks of low-ionization emission lines to trace the Be star $V_{r}$ and the He I lines to trace the companion star $V_{r}$ yields a circular orbital solution with a 93.8-day period and masses of $M_{\rm Be}=6.2$ $M_{\rm \odot}$ and $M_{\rm sdOB}=0.9$ $M_{\rm \odot}$ in the case of $i=80^{\circ}$. HD 55606 exhibits a variety of phase-locked variability, including the development of shell lines twice per orbit. The shell phases coincide with variation in the double emission peak separations, and both forms of variability are likely caused by a two-armed spiral density perturbation in the Be disk. The intensity ratios of the double emission peaks are also phase-locked, possibly indicating heating by the sdOB star of the side of the Be disk facing it. HD 55606 is a new member of the growing sample of Be+sdOB binaries, in which the Be star's rapid rotation and ability to form a disk can be attributed to past mass transfer.

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.

J-PLUS is an ongoing 12-band photometric optical survey, observing thousands of square degrees of the Northern hemisphere from the dedicated JAST/T80 telescope at the Observatorio Astrof\'isico de Javalambre. T80Cam is a 2 sq.deg field-of-view camera mounted on this 83cm-diameter telescope, and is equipped with a unique system of filters spanning the entire optical range. This filter system is a combination of broad, medium and narrow-band filters, optimally designed to extract the rest-frame spectral features (the 3700-4000\AA\ Balmer break region, H$\delta$, Ca H+K, the G-band, the Mgb and Ca triplets) that are key to both characterize stellar types and to deliver a low-resolution photo-spectrum for each pixel of the sky observed. With a typical depth of AB $\sim 21.25$ mag per band, this filter set thus allows for an indiscriminate and accurate characterization of the stellar population in our Galaxy, it provides an unprecedented 2D photo-spectral information for all resolved galaxies in the local universe, as well as accurate photo-z estimates ($\Delta\,z\sim 0.01-0.03$) for moderately bright (up to $r\sim 20$ mag) extragalactic sources. While some narrow band filters are designed for the study of particular emission features ([OII]/$\lambda$3727, H$\alpha$/$\lambda$6563) up to $z < 0.015$, they also provide well-defined windows for the analysis of other emission lines at higher redshifts. As a result, J-PLUS has the potential to contribute to a wide range of fields in Astrophysics, both in the nearby universe (Milky Way, 2D IFU-like studies, stellar populations of nearby and moderate redshift galaxies, clusters of galaxies) and at high redshifts (ELGs at $z\approx 0.77, 2.2$ and $4.4$, QSOs, etc). With this paper, we release $\sim 36$ sq.deg of J-PLUS data, containing about $1.5\times 10^5$ stars and $10^5$ galaxies at $r<21$ mag.

PKS 0625-354 (z=0.055) was observed with the four H.E.S.S. telescopes in 2012 during 5.5 hours. The source was detected above an energy threshold of 200 GeV at a significance level of 6.1$\sigma$. No significant variability is found in these observations. The source is well described with a power-law spectrum with photon index $\Gamma =2.84 \pm 0.50_{stat} \pm 0.10_{syst}$ and normalization (at $E_0$=1.0 TeV) $N_0(E_0)=(0.58 \pm 0.22_{stat} \pm 0.12_{syst})\times10^{-12}$ TeV$^{-1}$cm$^{-2}$s$^{-1}$. Multi-wavelength data collected with Fermi-LAT, Swift-XRT, Swift-UVOT, ATOM and WISE are also analysed. Significant variability is observed only in the Fermi-LAT $\gamma$-ray and Swift-XRT X-ray energy bands. Having a good multi-wavelength coverage from radio to very high energy, we performed a broadband modelling from two types of emission scenarios. The results from a one zone lepto-hadronic, and a multi-zone leptonic models are compared and discussed. On the grounds of energetics, our analysis favours a leptonic multi-zone model. Models associated to the X-ray variability constraint supports previous results suggesting a BL Lac nature of PKS 0625-354, with, however, a large-scale jet structure typical of a radio galaxy.

Shell-type supernova remnants (SNRs) are considered prime candidates for the acceleration of Galactic cosmic rays (CRs) up to the knee of the CR spectrum at $\mathrm{E} \approx \mathrm{3}\times \mathrm{10}^\mathrm{15}$ eV. Our Milky Way galaxy hosts more than 350 SNRs discovered at radio wavelengths and at high energies, of which 220 fall into the H.E.S.S. Galactic Plane Survey (HGPS) region. Of those, only 50 SNRs are coincident with a H.E.S.S source and in 8 cases the very high-energy (VHE) emission is firmly identified as an SNR. The H.E.S.S. GPS provides us with a legacy for SNR population study in VHE $\gamma$-rays and we use this rich data set to extract VHE flux upper limits from all undetected SNRs. Overall, the derived flux upper limits are not in contradiction with the canonical CR paradigm. Assuming this paradigm holds true, we can constrain typical ambient density values around shell-type SNRs to $n\leq 7~\textrm{cm}^\textrm{-3}$ and electron-to-proton energy fractions above 10~TeV to $\epsilon_\textrm{ep} \leq 5\times 10^{-3}$. Furthermore, comparisons of VHE with radio luminosities in non-interacting SNRs reveal a behaviour that is in agreement with the theory of magnetic field amplification at shell-type SNRs.

In order to study the growth and evolution of circumstellar disks around classical Be stars, we analyze optical time-series photometry from the KELT survey with simultaneous infrared and visible spectroscopy from the APOGEE survey and BeSS database for a sample of 160 Galactic classical Be stars. The systems studied here show variability including transitions from a diskless to a disk-possessing state (and vice versa), and persistent disks that vary in strength, being replenished at either regularly or irregularly occurring intervals. We detect disk-building events (outbursts) in the light curves of 28\% of our sample. Outbursts are more commonly observed in early- (57\%), compared to mid- (27\%) and late-type (8\%) systems. A given system may show anywhere between 0 -- 40 individual outbursts in its light curve, with amplitudes ranging up to $\sim$0.5 mag and event durations between $\sim$2 -- 1000 days. We study how both the photometry and spectroscopy change together during active episodes of disk growth or dissipation, revealing details about the evolution of the circumstellar environment. We demonstrate that photometric activity is linked to changes in the inner disk, and show that, at least in some cases, the disk growth process is asymmetrical. Observational evidence of Be star disks both growing and clearing from the inside out is presented. The duration of disk buildup and dissipation phases are measured for 70 outbursts, and we find that the average outburst takes about twice as long to dissipate as it does to build up in optical photometry. Our analysis hints that dissipation of the inner disk occurs relatively slowly for late-type Be stars.

Context. Recently, the high-energy (HE, 0.1-100 GeV) $\gamma$-ray emission from the object LMC P3 in the Large Magellanic Cloud (LMC) has been discovered to be modulated with a 10.3-day period, making it the first extra-galactic $\gamma$-ray binary. Aims. This work aims at the detection of very-high-energy (VHE, >100 GeV) $\gamma$-ray emission and the search for modulation of the VHE signal with the orbital period of the binary system. Methods. LMC P3 has been observed with the High Energy Stereoscopic System (H.E.S.S.); the acceptance-corrected exposure time is 100 h. The data set has been folded with the known orbital period of the system in order to test for variability of the emission. Energy spectra are obtained for the orbit-averaged data set, and for the orbital phase bin around the VHE maximum. Results. VHE $\gamma$-ray emission is detected with a statistical significance of 6.4 $\sigma$. The data clearly show variability which is phase-locked to the orbital period of the system. Periodicity cannot be deduced from the H.E.S.S. data set alone. The orbit-averaged luminosity in the $1-10$ TeV energy range is $(1.4 \pm 0.2) \times 10^{35}$ erg/s. A luminosity of $(5 \pm 1) \times 10^{35}$ erg/s is reached during 20% of the orbit. HE and VHE $\gamma$-ray emissions are anti-correlated. LMC P3 is the most luminous $\gamma$-ray binary known so far.

B[e] supergiants (SGs) are massive post-main-sequence stars, surrounded by a complex circumstellar (CS) environment. The aim of this work is to investigate the structure and kinematics of the CS disc of the B[e] SG LHA 120-S 35. We used high-resolution optical spectra obtained in different years to model the forbidden emission lines and determine the kinematical properties of their line-forming regions, assuming Keplerian rotation. We also used low-resolution near-infrared (IR) spectra to explore the variability of molecular emission. LHA 120-S 35 displays spectral variability in both optical and IR regions. The P-Cygni line profiles of H I, as well as those of Fe II and O I, suggest the presence of a strong bipolar clumped wind. We distinguish density enhancements in the P-Cygni absorption component of the first Balmer lines, which show variations in both velocity and strength. The P-Cygni profile emission component is double-peaked, indicating the presence of a rotating CS disc. We also observe line-profile variations in the permitted and forbidden features of Fe II and O I. In the IR, we detect variations in the intensity of the H I emission lines as well as in the emission of the CO band-heads. Moreover, we find that the profiles of each [Ca II] and [O I] emission lines contain contributions from spatially different (complete or partial) rings. Globally, we find evidence of detached multi-ring structures, revealing density variations along the disc. We suggest that LHA 120-S 35 has passed through the red-supergiant (RSG) phase and evolves back bluewards in the Hertzsprung-Russell diagram. The formation of the complex CS structure could be the result of the wind-wind interactions of the post-RSG wind with the previously ejected material from the RSG. However, the presence of a binary companion can not be excluded. Finally, we find that LHA 120-S 35 belongs to a young stellar cluster.

We performed a spectroscopic and photometric analysis to study new eruptions in two luminous blue variables (LBVs) in the Magellanic Clouds. We detected a strong new eruption in the LBV R40 that reached $V \sim 9.2$ in 2016, which is around $1.3$ mag brighter than the minimum registered in 1985. During this new eruption, the star changed from an A-type to a late F-type spectrum. Based on photometric and spectroscopic empirical calibrations and synthetic spectral modeling, we determine that R\,40 reached $T_{\mathrm{eff}} = 5800-6300$~K during this new eruption. This object is thereby probably one of the coolest identified LBVs. We could also identify an enrichment of nitrogen and r- and s-process elements. We detected a weak eruption in the LBV R 110 with a maximum of $V \sim 9.9$ mag in 2011, that is, around $1.0$ mag brighter than in the quiescent phase. On the other hand, this new eruption is about $0.2$ mag fainter than the first eruption detected in 1990, but the temperature did not decrease below 8500 K. Spitzer spectra show indications of cool dust in the circumstellar environment of both stars, but no hot or warm dust was present, except by the probable presence of PAHs in R\,110. We also discuss a possible post-red supergiant nature for both stars.

We search for high-energy gamma-ray emission from the binary neutron star merger GW170817 with the H.E.S.S. Imaging Air Cherenkov Telescopes. The observations presented here have been obtained starting only 5.3h after GW170817. The H.E.S.S. target selection identified regions of high probability to find a counterpart of the gravitational wave event. The first of these regions contained the counterpart SSS17a that has been identified in the optical range several hours after our observations. We can therefore present the first data obtained by a ground-based pointing instrument on this object. A subsequent monitoring campaign with the H.E.S.S. telescopes extended over several days, covering timescales from 0.22 to 5.2 days and energy ranges between $270\,\mathrm{GeV}$ to $8.55\,\mathrm{TeV}$. No significant gamma-ray emission has been found. The derived upper limits on the very-high-energy gamma-ray flux for the first time constrain non-thermal, high-energy emission following the merger of a confirmed binary neutron star system.

The H.E.S.S. collaboration has discovered a new very high energy (VHE, E $>$ 0.1 TeV) $\gamma$-ray source, HESS J1741-302, located in the Galactic plane. Despite several attempts to constrain its nature, no plausible counterpart has been found so far at X-ray and MeV/GeV $\gamma$-ray energies, and the source remains unidentified. An analysis of 145-hour of observations of HESS J1741-302 at VHEs has revealed a steady and relatively weak TeV source ($\sim$1$\%$ of the Crab Nebula flux), with a spectral index of $\Gamma$ = 2.3 $\pm$ 0.2$_{\text{stat}}$ $\pm$ 0.2$_{\text{sys}}$, extending to energies up to 10 TeV without any clear signature of a cut-off. In a hadronic scenario, such a spectrum implies an object with particle acceleration up to energies of several hundred TeV. Contrary to most H.E.S.S. unidentified sources, the angular size of HESS J1741-302 is compatible with the H.E.S.S. point spread function at VHEs, with an extension constrained to be below 0.068$^{\circ}$ at a 99$\%$ confidence level. The $\gamma$-ray emission detected by H.E.S.S. can be explained both within a hadronic scenario, due to collisions of protons with energies of hundreds of TeV with dense molecular clouds, and in a leptonic scenario, as a relic pulsar wind nebula, possibly powered by the middle-aged (20 kyr) pulsar PSR B1737-30. A binary scenario, related to the compact radio source 1LC 358.266+0.038 found to be spatially coincident with the best fit position of HESS J1741-302, is also envisaged.

The first results of radiative transfer calculations on decretion discs of binary Be stars are presented. A smoothed particle hydrodynamics code computes the structure of Be discs in coplanar circular binary systems for a range of orbital and disc parameters. The resulting disc configuration consists of two spiral arms, and can be given as input into a Monte Carlo code, which calculates the radiative transfer along the line of sight for various observational coordinates. Making use of the property of steady disc structure in coplanar circular binaries, observables are computed as functions of the orbital phase. Orbital-phase series of line profiles are given for selected parameter sets under various viewing angles, to allow comparison with observations. Flat-topped profiles with and without superimposed multiple structures are reproduced, showing, for example, that triple-peaked profiles do not have to be necessarily associated with warped discs and misaligned binaries. It is demonstrated that binary tidal effects give rise to phase-locked variability of the violet-to-red (V/R) ratio of hydrogen emission lines. The V/R ratio exhibits two maxima per cycle; in certain cases those maxima are equal, leading to a clear new V/R cycle every half orbital period. This study opens a way in identifying binaries and in constraining the parameters of binary systems that exhibit phase-locked variations induced by tidal interaction with a companion star.

The Galactic object MWC 137 was suggested to belong to the group of B[e] supergiants. However, with its large-scale optical bipolar ring nebula and the high velocity jet and knots, it is a rather atypical representative of this class. We performed multi-wavelength observations spreading from the optical to the radio regime. Based on optical imaging and long-slit spectroscopic data we found that the northern parts of the large-scale nebula are predominantly blue-shifted, while the southern regions appear mostly red-shifted. We developed a geometrical model consisting of two double-cones. While various observational features can be approximated with such a scenario, the observed velocity pattern is more complex. Using near-infrared integral-field unit spectroscopy we studied the hot molecular gas in the close vicinity of the star. The emission from the hot CO gas arises in a small-scale disk revolving around the star on Keplerian orbits. While the disk itself cannot be spatially resolved, its emission is reflected by dust arranged in arc-like structures and clumps surrounding MWC 137 on small scales. In the radio regime we mapped the cold molecular gas in the outskirts of the optical nebula. We found that large amounts of cool molecular gas and warm dust embrace the optical nebula in the east, south and west. No cold gas or dust were detected in the north and north-western regions. Despite the new insights on the nebula kinematics gained from our studies, the real formation scenario of the large-scale nebula remains an open issue.

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

Very high-energy gamma-rays (VHE, E>100 GeV) propagating over cosmological distances can interact with the low-energy photons of the extragalactic background light (EBL) and produce electron-positron pairs. The transparency of the universe to VHE gamma-rays is then directly related to the spectral energy distribution (SED) of the EBL. The observation of features in the VHE energy spectra of extragalactic sources allows the EBL to be measured, which otherwise is very difficult to determine. An EBL-model independent measurement of the EBL SED with the H.E.S.S. array of Cherenkov telescopes is presented. It is obtained by extracting the EBL absorption signal from the reanalysis of high-quality spectra of blazars. From H.E.S.S. data alone the EBL signature is detected at a significance of 9.5 sigma, and the intensity of the EBL obtained in different spectral bands is presented together with the associated gamma-ray horizon.

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

The diffuse very high-energy (VHE, > 100 GeV) gamma-ray emission observed in the central 200 pc of the Milky Way by H.E.S.S. was found to follow the dense matter distribution in the Central Molecular Zone (CMZ) up to a longitudinal distance of about 130 pc to the Galactic Centre (GC), where the flux rapidly decreases. Recent sensitive H.E.S.S. analysis revealed that the cosmic-ray (CR) density profile drops with the distance to the centre, making data compatible with a steady cosmic PeVatron at the GC. In this paper, we extend this analysis to obtain for the first time a detailed characterisation of the correlation with matter and to search for additional features and individual gamma-ray sources in the inner 200 pc. Taking advantage of 250 hours of H.E.S.S. data and improved analysis techniques we perform a detailed morphology study of the diffuse VHE emission observed from the GC ridge and reconstruct its total spectrum. We show that the emission correlated with dense matter covers the full CMZ and that its flux is about half the total diffuse emission flux. We also detect some emission at higher latitude likely produced by hadronic collisions of CRs in less dense regions of the GC interstellar medium. We detect an additional emission component centred on the GC and extending over about 15 pc that is consistent with the existence of a strong CR density gradient and confirms the presence of a CR accelerator at the very centre of our Galaxy. We show that the spectrum of the full ridge diffuse emission is compatible with the one previously derived from the central regions, suggesting that a single population of particles fills the entire CMZ. Finally, we report the discovery of a VHE gamma-ray source near the GC radio arc and argue that it is produced by the pulsar wind nebula candidate G0.13-0.11.

High spatial resolution is the key for the understanding various astrophysical phenomena. But even with the future E-ELT, single dish instruments are limited to a spatial resolution of about 4 mas in the visible. For the closest objects within our Galaxy most of the stellar photosphere remains smaller than 1 mas. With the success of long baseline interferometry these limitations were soom overcome. Today low and high resolution interferometric instruments on the VLTI and CHARA offer an immense range of astrophysical studies. Combining more telescopes and moving to visible wavelengths broadens the science cases even more. With the idea of developing strong science cases for a future visible interferometer, we organized a science group around the following topics: pre-main sequence and main sequence stars, fundamental parameters, asteroseismology and classical pulsating stars, evolved stars, massive stars, active galactic nuclei (AGNs) and imaging techniques. A meeting was organized on the 15th and 16th of January, 2015 in Nice with the support of the Action Specific in Haute Resolution Angulaire (ASHRA), the Programme National en Physique Stellaire (PNPS), the Lagrange Laboratory and the Observatoire de la Cote d'Azur, in order to present these cases and to discuss them further for future visible interferometers. This White Paper presents the outcome of the exchanges. This book is dedicated to the memory of our colleague Olivier Chesneau who passed away at the age of 41.

The nine-year H.E.S.S. Galactic Plane Survey (HGPS) yielded the most uniform observation scan of the inner Milky Way in the TeV gamma-ray band to date. The sky maps and source catalogue of the HGPS allow for a systematic study of the population of TeV pulsar wind nebulae found throughout the last decade. To investigate the nature and evolution of pulsar wind nebulae, for the first time we also present several upper limits for regions around pulsars without a detected TeV wind nebula. Our data exhibit a correlation of TeV surface brightness with pulsar spin-down power $\dot{E}$. This seems to be caused both by an increase of extension with decreasing $\dot{E}$, and hence with time, compatible with a power law $R_\mathrm{PWN}(\dot{E}) \sim \dot{E}^{-0.65 \pm 0.20}$, and by a mild decrease of TeV gamma-ray luminosity with decreasing $\dot{E}$, compatible with $L_{1-10\,\mathrm{TeV}} \sim \dot{E}^{0.59 \pm 0.21}$. We also find that the offsets of pulsars with respect to the wind nebula centres with ages around 10 kyr are frequently larger than can be plausibly explained by pulsar proper motion and could be due to an asymmetric environment. In the present data, it seems that a large pulsar offset is correlated with a high apparent TeV efficiency $L_{1-10\,\mathrm{TeV}}/\dot{E}$. In addition to 14 HGPS sources considered as firmly identified pulsar wind nebulae and 5 additional pulsar wind nebulae taken from literature, we find 10 HGPS sources that form likely TeV pulsar wind nebula candidates. Using a model that subsumes the present common understanding of the very high-energy radiative evolution of pulsar wind nebulae, we find that the trends and variations of the TeV observables and limits can be reproduced to a good level, drawing a consistent picture of present-day TeV data and theory.

The peculiar emission-line star MWC 137 with its extended optical nebula was recently classified as B[e] supergiant. To study the spatial distribution of its circumstellar molecular gas on small and large scales, we obtained near-infrared and radio observations using SINFONI and APEX, respectively. We find that the hot CO gas is arranged in moving clumpy ring and shell structures close to the star, while a cold CO envelope is encircling the borders of the optical nebula from the south to the west.

Studying the temporal variability of BL Lac objects at the highest energies provides unique insights into the extreme physical processes occurring in relativistic jets and in the vicinity of super-massive black holes. To this end, the long-term variability of the BL Lac object PKS 2155-304 is analyzed in the high (HE, 100 MeV < E < 300 GeV) and very high energy (VHE, E > 200 GeV) gamma-ray domain. Over the course of ~9 yr of H.E.S.S observations the VHE light curve in the quiescent state is consistent with a log-normal behavior. The VHE variability in this state is well described by flicker noise (power-spectral-density index {\ss}_VHE = 1.10 +0.10 -0.13) on time scales larger than one day. An analysis of 5.5 yr of HE Fermi LAT data gives consistent results ({\ss}_HE = 1.20 +0.21 -0.23, on time scales larger than 10 days) compatible with the VHE findings. The HE and VHE power spectral densities show a scale invariance across the probed time ranges. A direct linear correlation between the VHE and HE fluxes could neither be excluded nor firmly established. These long-term-variability properties are discussed and compared to the red noise behavior ({\ss} ~ 2) seen on shorter time scales during VHE-flaring states. The difference in power spectral noise behavior at VHE energies during quiescent and flaring states provides evidence that these states are influenced by different physical processes, while the compatibility of the HE and VHE long-term results is suggestive of a common physical link as it might be introduced by an underlying jet-disk connection.

B[e] Supergiants are a phase in the evolution of some massive stars for which we have observational evidence but no predictions by any stellar evolution model. The mass-loss during this phase creates a complex circumstellar environment with atomic, molecular, and dust regions usually found in rings or disk-like structures. However, the detailed structure and the formation of the circumstellar environment are not well-understood, requiring further investigation. To address that we initiated an observing campaign to obtain a homogeneous set of high-resolution spectra in both the optical and NIR (using MPG-ESO/FEROS, GEMINI/Phoenix and VLT/CRIRES, respectively). We monitor a number of Galactic B[e] Supergiants, for which we examined the [OI] and [CaII] emission lines and the bandheads of the CO and SiO molecules to probe the structure and the kinematics of their formation regions. We find that the emission from each tracer forms either in a single or in multiple equatorial rings.

A search for dark matter line-like signals was performed in the vicinity of the Galactic Centre by the H.E.S.S. experiment on observational data taken in 2014. An unbinned likelihood analysis was developed to improve the sensitivity to line-like signals. The upgraded analysis along with newer data extend the energy coverage of the previous measurement down to 100 GeV. The 18 h of data collected with the H.E.S.S. array allow one to rule out at 95% CL the presence of a 130 GeV line (at $l = -1.5^{\circ}, b = 0^{\circ}$ and for a dark matter profile centered at this location) previously reported in Fermi-LAT data. This new analysis overlaps significantly in energy with previous Fermi-LAT and H.E.S.S. results. No significant excess associated with dark matter annihilations was found in the energy range 100 GeV to 2 TeV and upper limits on the gamma-ray flux and the velocity weighted annihilation cross-section are derived adopting an Einasto dark matter halo profile. Expected limits for present and future large statistics H.E.S.S. observations are also given.

The supernova remnant (SNR) W49B originated from a core-collapse supernova that occurred between one and four thousand years ago, and subsequently evolved into a mixed-morphology remnant, which is interacting with molecular clouds (MC). $\gamma$-ray observations of SNR/MC associations are a powerful tool to constrain the origin of Galactic cosmic-rays, as they can probe the acceleration of hadrons through their interaction with the surrounding medium and subsequent emission of non-thermal photons. The detection of a $\gamma$-ray source coincident with W49B at very high energies (VHE; E > 100 GeV) with the H.E.S.S. Cherenkov telescopes is reported together with a study of the source with 5 years of Fermi-LAT high energy $\gamma$-ray (0.06 - 300 GeV) data. The smoothly-connected combined source spectrum, measured from 60 MeV to multi-TeV energies, shows two significant spectral breaks at $304\pm20$ MeV and $8.4_{-2.5}^{+2.2}$ GeV, the latter being constrained by the joint fit from the two instruments. The detected spectral features are similar to those observed in several other SNR/MC associations and are found to be indicative of $\gamma$-ray emission produced through neutral-pion decay.

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

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

We aim to improve our knowledge on the structure and dynamics of the circumstellar disk of the LMC B[e] supergiant LHA 120-S 73. High-resolution optical and near-IR spectroscopic data were obtained over a period of 16 and 7 years, respectively. The spectra cover the diagnostic emission lines from [CaII] and [OI], as well as the CO bands. These features trace the disk at different distances from the star. We analyzed the kinematics of the individual emission regions by modeling their emission profiles. A low-resolution mid-infrared spectrum was obtained as well, which provides information on the composition of the dusty disk. All diagnostic emission features display double-peaked line profiles, which we interpret as due to Keplerian rotation. We find that LHA 120-S 73 is surrounded by at least four individual rings of material with alternating densities (or by a disk with strongly non-monotonic radial density distribution). Moreover, we find that the molecular ring must have gaps or at least strong density inhomogeneities, or in other words, a clumpy structure. The mid-infrared spectrum displays features of oxygen- and carbon-rich grain species, which indicates a long-lived, stable dusty disk. We cannot confirm the previously reported high value for the stellar rotation velocity. The line profile of HeI 5876 A is strongly variable in both width and shape and resembles of those seen in non-radially pulsating stars. A proper determination of the real underlying stellar rotation velocity is hence not possible. The existence of multiple stable and clumpy rings of alternating density recalls ring structures around planets. Although there is currently insufficient observational evidence, it is tempting to propose a scenario with one (or more) minor bodies or planets revolving around LHA 120-S 73 and stabilizing the ring system, in analogy to the shepherd moons in planetary systems.

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

The Pierre Auger Observatory has begun a major Upgrade of its already impressive capabilities, with an emphasis on improved mass composition determination using the surface detectors of the Observatory. Known as AugerPrime, the upgrade will include new 4 m$^2$ plastic scintillator detectors on top of all 1660 water-Cherenkov detectors, updated and more flexible surface detector electronics, a large array of buried muon detectors, and an extended duty cycle for operations of the fluorescence detectors. This Preliminary Design Report was produced by the Collaboration in April 2015 as an internal document and information for funding agencies. It outlines the scientific and technical case for AugerPrime. We now release it to the public via the arXiv server. We invite you to review the large number of fundamental results already achieved by the Observatory and our plans for the future.

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

Aims: We aim for an understanding of the morphological and spectral properties of the supernova remnant RCW~86 and for insights into the production mechanism leading to the RCW~86 very high-energy gamma-ray emission. Methods: We analyzed High Energy Spectroscopic System data that had increased sensitivity compared to the observations presented in the RCW~86 H.E.S.S. discovery publication. Studies of the morphological correlation between the 0.5-1~keV X-ray band, the 2-5~keV X-ray band, radio, and gamma-ray emissions have been performed as well as broadband modeling of the spectral energy distribution with two different emission models. Results:We present the first conclusive evidence that the TeV gamma-ray emission region is shell-like based on our morphological studies. The comparison with 2-5~keV X-ray data reveals a correlation with the 0.4-50~TeV gamma-ray emission.The spectrum of RCW~86 is best described by a power law with an exponential cutoff at $E_{cut}=(3.5\pm 1.2_{stat})$ TeV and a spectral index of $\Gamma$~$1.6\pm 0.2$. A static leptonic one-zone model adequately describes the measured spectral energy distribution of RCW~86, with the resultant total kinetic energy of the electrons above 1 GeV being equivalent to $\sim$0.1\% of the initial kinetic energy of a Type I a supernova explosion. When using a hadronic model, a magnetic field of $B$~100$\mu$G is needed to represent the measured data. Although this is comparable to formerly published estimates, a standard E$^{-2}$ spectrum for the proton distribution cannot describe the gamma-ray data. Instead, a spectral index of $\Gamma_p$~1.7 would be required, which implies that ~$7\times 10^{49}/n_{cm^{-3}}$erg has been transferred into high-energy protons with the effective density $n_{cm^{-3}}=n/ 1$ cm^-3. This is about 10\% of the kinetic energy of a typical Type Ia supernova under the assumption of a density of 1~cm^-3.

A measurement of the cosmic-ray spectrum for energies exceeding $4{\times}10^{18}$ eV is presented, which is based on the analysis of showers with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger Observatory between 1 January 2004 and 31 December 2013. The measured spectrum confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$ eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by a smooth suppression region. For the energy ($E_\text{s}$) at which the spectral flux has fallen to one-half of its extrapolated value in the absence of suppression, we find $E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$ eV.