A standard expectation of primordial cosmological inflation is that it dilutes all relics created before its onset to unobservable levels. We present a counterexample to this expectation by demonstrating that a network of cosmic strings diluted by inflation can regrow to a level that is potentially observable today in gravitational waves~(GWs). In contrast to undiluted cosmic strings, whose primary GW signals are typically in the form of a stochastic GW background, the leading signal from a diluted cosmic string network can be distinctive bursts of GWs within the sensitivity reach of current and future GW observatories.

Cosmic strings are generically predicted in many extensions of the Standard Model of particle physics. We propose a new avenue for detecting cosmic strings through their effect on the filamentary structure in the cosmic web. Using cosmological simulations of the density wake from a cosmic string, we examine a variety of filament structure probes. We show that the largest effect of the cosmic string is an overdensity in the filament distribution around the string wake. The signal from the overdensity is stronger at higher redshift, and more robust with a wider field. We analyze the spatial distribution of filaments from a publicly available catalog of filaments built from SDSS galaxies. With existing data, we find no evidence for the presence of a cosmic string wake with string tension parameter $G\mu$ above $\sim 5\times 10^{-6}$. However, we project WFIRST will be able to detect a signal from such a wake at the $99\%$ confidence level at redshift $z=2$, with significantly higher confidence and the possibility of probing lower tensions ($G\mu \sim 10^{-6}$), at $z=10$. The sensitivity of this method is not competitive with constraints derived from the CMB. However, it provides an independent discovery channel at low redshift, which could be a smoking-gun in scenarios where the CMB bound can be weakened.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

We propose a new mechanism where asymmetric dark matter (ADM) and the baryon asymmetry are both generated in the same decay chain of a metastable weakly interacting massive particle (WIMP) after its thermal freeze-out. Dark matter and baryons are connected by a generalized baryon number that is conserved, while the DM asymmetry and baryon asymmetry compensate each other. This unified framework addresses the DM-baryon coincidence while inheriting the merit of the conventional WIMP miracle in predicting relic abundances of matter. Examples of renormalizable models realizing this scenario are presented. These models generically predict ADM with sub-GeV to GeV-scale mass that interacts with Standard Model quarks or leptons, thus rendering potential signatures at direct detection experiments sensitive to low mass DM. Other interesting phenomenological predictions are also discussed, including: LHC signatures of new intermediate particles with color or electroweak charge and DM induced nucleon decay; the long-lived WIMP may be within reach of future high energy collider experiments.

We report observations of a transient source \fermi\ from radio to \grs. \fermi\ was discovered by the {\it Fermi-LAT} in May 2017. Follow-up {\it Swift-XRT} observations revealed three flaring episodes through March 2018, and the peak X-ray flux is about $10^3$ higher than the {\it ROSAT all-sky survey (RASS)} flux upper limit. Optical spectral measurements taken by the {\it Magellan 6.5-m telescope} and the {\it Lick-Shane telescope} both show a largely featureless spectrum, strengthening the BL Lac interpretation first proposed by \citet{Bruni18}. The optical and mid-infrared (MIR) emission goes to a higher state in 2018, when the flux in high energies goes down to a lower level. Our {\it RATAN-600m} measurements at 4.8~GHz and 8.2~GHz do not indicate any significant radio flux variation over the monitoring seasons in 2017 and 2018, nor deviate from the archival {\it NVSS} flux level. During GeV flaring times, the spectrum is very hard ($\Gamma_\gamma\sim$1.7) in the GeV band and at times also very hard (($\Gamma_{\rm X}\lesssim2$) in the X-rays, similar to a high-synchrotron-peak (or even an extreme) BL Lac object, making \fermi\ a good target for ground-based {\it Cherenkov telescopes}.

We investigate the nature of the unidentified very-high-energy (VHE) gamma-ray object, HESS J1832-093, in a multi-wavelength context. Based on X-ray variability and spectral index ($\Gamma_X\sim\,1.5$), and its broad-band spectrum (which was remarkably similar to HESS J0632+057, a confirmed gamma-ray binary), HESS J1832-093 has been considered to be a strong gamma-ray binary candidate in previous works. In this work, we provide further evidence for this scenario. We obtained a spectrum of its IR counterpart using Gemini/Flamingo, finding absorption lines that are usually seen in massive stars, in particular O stars. We also obtained a rather steep ATCA spectrum ($\alpha=-1.18^{+1.04}_{-0.88}$) which prefers a gamma-ray binary over an AGN scenario. Based on spatial-spectral analysis and variability search, we found that 4FGL J1832.9-0913 is possible to be associated with SNR G22.7-0.2 rather than with HESS J1832-093 only.

Boosted dark matter (BDM) is a well-motivated class of dark matter candidates in which a small subcomponent is relativistic at the present time. In this work, we lay the foundation for BDM searches via its hadronic interactions at DUNE. We investigate BDM-nucleus scattering in detail by developing new dedicated event generation techniques with a parameterized detector simulation. We demonstrate that the excellent event reconstruction capability of DUNE enables a significantly broader discovery reach.

A zero initial velocity of the axion field is assumed in the conventional misalignment mechanism. We propose an alternative scenario where the initial velocity is nonzero, which may arise from an explicit breaking of the PQ symmetry in the early Universe. We demonstrate that, depending on the specifics about the initial velocity and the time order of the PQ symmetry breaking vs. inflation, this new scenario can alter the conventional prediction for the axion relic abundance in different, potentially significant ways. As a result, new viable parameter regions for axion dark matter may open up.

We study the kinematics of the M87 jet using the first year data of the KVN and VERA Array (KaVA) large program, which has densely monitored the jet at 22 and 43 GHz since 2016. We find that the apparent jet speeds generally increase from $\approx0.3c$ at $\approx0.5$ mas from the jet base to $\approx2.7c$ at $\approx20$ mas, indicating that the jet is accelerated from subluminal to superluminal speeds on these scales. We perform a complementary jet kinematic analysis by using archival Very Long Baseline Array monitoring data observed in $2005-2009$ at 1.7 GHz and find that the jet is moving at relativistic speeds up to $\approx5.8c$ at distances of $200-410$ mas. We combine the two kinematic results and find that the jet is gradually accelerated over a broad distance range that coincides with the jet collimation zone, implying that conversion of Poynting flux to kinetic energy flux takes place. If the jet emission consists of a single streamline, the observed trend of jet acceleration ($\Gamma\propto z^{0.16\pm0.01}$) is relatively slow compared to models of a highly magnetized jet. This indicates that Poynting flux conversion through the differential collimation of poloidal magnetic fields may be less efficient than expected. However, we find a non-negligible dispersion in the observed speeds for a given jet distance, making it difficult to describe the jet velocity field with a single power-law acceleration function. We discuss the possibility that the jet emission consists of multiple streamlines following different acceleration profiles, resulting in jet velocity stratification.

The supernova remnant (SNR) HESS J1731-347 is a young SNR which displays a non-thermal X-ray and TeV shell structure. A molecular cloud at a distance of 3.2 kpc is spatially coincident with the western part of the SNR, and it is likely hit by the SNR. The X-ray emission from this part of the shell is much lower than from the rest of the SNR. Moreover, a compact GeV emission region coincident with the cloud has been detected with a soft spectrum. These observations seem to imply a shock-cloud collision scenario at this area, where the stalled shock can no longer accelerate super-TeV electrons or maintain strong magnetic turbulence downstream, while the GeV cosmic rays (CRs) are released through this stalled shock. To test this hypothesis, we have performed a detailed Fermi-LAT reanalysis of the HESS J1731-347 region with over 9 years of data. We find that the compact GeV emission region displays a spectral power-law index of -2.4, whereas the GeV emission from the rest of the SNR (excluding the cloud region) has an index of -1.8. A hadronic model involving a shock-cloud collision scenario is built to explain the -ray emission from this area. It consists of three CR sources: run-away super-TeV CRs that have escaped from the fast shock, leaked GeV CRs from the stalled shock, and the local CR sea. The X-ray and -ray emission of the SNR excluding the shock-cloud interaction region is explained in a one-zone leptonic model. Our shock-cloud collision model explains well the GeV-TeV observations from both cloud regions around HESS J1731-347, i.e. from the cloud in contact with the SNR and from the more distant cloud which is coincident with the nearby TeV source HESS J1729-345. We find however that the leaked GeV CRs from the shock-cloud collision do not necessarily dominate the GeV emission from the clouds, due to a comparable contribution from the local CR sea.

We propose in this White Paper a concept for a space experiment using cold atoms to search for ultra-light dark matter, and to detect gravitational waves in the frequency range between the most sensitive ranges of LISA and the terrestrial LIGO/Virgo/KAGRA/INDIGO experiments. This interdisciplinary experiment, called Atomic Experiment for Dark Matter and Gravity Exploration (AEDGE), will also complement other planned searches for dark matter, and exploit synergies with other gravitational wave detectors. We give examples of the extended range of sensitivity to ultra-light dark matter offered by AEDGE, and how its gravitational-wave measurements could explore the assembly of super-massive black holes, first-order phase transitions in the early universe and cosmic strings. AEDGE will be based upon technologies now being developed for terrestrial experiments using cold atoms, and will benefit from the space experience obtained with, e.g., LISA and cold atom experiments in microgravity. This paper is based on a submission (v1) in response to the Call for White Papers for the Voyage 2050 long-term plan in the ESA Science Programme. ESA limited the number of White Paper authors to 30. However, in this version (v2) we have welcomed as supporting authors participants in the Workshop on Atomic Experiments for Dark Matter and Gravity Exploration held at CERN: ({\tt https://indico.cern.ch/event/830432/}), as well as other interested scientists, and have incorporated additional material.

Global cosmic strings are generically predicted in particle physics beyond the Standard Model, e.g., a post-inflationary global $U(1)$ symmetry breaking which may associate with axion-like dark matter. We demonstrate that although subdominant to Goldstone emission, gravitational waves (GWs) radiated from global strings can be observable with current/future GW detectors. The frequency spectrum of such GWs is also shown to be a powerful tool to probe the Hubble expansion rate of the Universe at times prior to the Big Bang nucleosynthesis where the standard cosmology has yet to be tested.

We report X-ray data analysis results obtained from Chandra, XMM-Newton, NuSTAR and Swift observations of PSR J2032+4127 taken before, during, and after the periastron on 2017 November 13. We found the first clear evidence of a change in the X-ray spectral index over the passage period, thanks to a broad and sensitive spectral coverage by XMM-Newton and NuSTAR. We analysed the joint XMM-Newton and NuSTAR observation epochs with power-law and broken power-law model. We have obtained change in spectral parameters before and after the periastron passage for both models. The spectra get softened after the passage. The evolution of the spectral index and break energy before and after the periastron may indicate a change in the physical state of shock-accelerated electrons.

We study the cosmological signatures of Invisibly Annihilating Dark Matter (IAnDM), where DM annihilates into dark radiation particles that are decoupled from the Standard Model (SM). In the simple benchmark model we consider here, such invisible annihilation determines the relic abundance of DM via dark thermal freeze-out. We demonstrate that IAnDM may reveal itself through observable, novel signatures that are correlated: scale-dependent $\Delta N_\text{eff}$ (number of extra effective neutrinos) in the Cosmic Microwave Background (CMB) spectrum due to DM residual annihilation, while the phase of acoustic peaks shifts towards the opposite direction relative to that due to SM neutrinos, resembling the effect due to scattering (fluidlike) thermal dark radiation; in addition, IAnDM induces modifications to the matter power spectrum that resemble those due to warm dark matter. Current data are sensitive to IAnDM with masses up to $\sim200$ keV, while future observations will improve this reach, especially if the late-time DM annihilation cross section is enhanced relative to the standard thermal value, which can be realized in a variety of scenarios. This work also demonstrates a new possibility of realizing thermal sub-MeV DM with observable signals.

We present the jet kinematics of the flat spectrum radio quasar (FSRQ) 4C +21.35 using time-resolved KaVA very long baseline interferometry array radio maps obtained from September 2014 to July 2016. During two out of three observing campaigns, observations were performed bi-weekly at 22 and 43 GHz quasi-simultaneously. At 22 GHz, we identified three jet components near the core with apparent speeds up to (14.4+/-2.1)c. The timing of the ejection of a new component detected in 2016 is consistent with a gamma-ray flare in November 2014. At 43 GHz, we found four inner jet (<3 mas) components with speeds from (3.5+/-1.4)c to (6.8+/-1.5)c. Jet component speeds tend to be higher with increasing distances from the core. We compared our data with archival Very Long Baseline Array (VLBA) data from the Boston University (BU) 43 GHz and the Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE) 15.4 GHz monitoring programs. Whereas MOJAVE data and our data are in good agreement, jet speeds obtained from the BU Program data in the same time period are about twice as high as the ones we obtain from the KaVA data. The discrepancy at 43 GHz indicates that radio arrays with different angular resolution identify and trace different jet features even when the data are obtained at the same frequency and at the same time. The flux densities of jet components decay exponentially, in agreement with a synchrotron cooling time scale of about 1 year. Using known electron Lorentz factor values (about 9,000), we estimate the magnetic field strength to be around 1-3 micro-Tesla. When adopting a jet viewing angle of 5 degrees, the intrinsic jet speed is of order 0.99c.

The KVN(Korean VLBI Network)-style simultaneous multi-frequency receiving mode is demonstrated to be promising for mm-VLBI observations. Recently, other Very long baseline interferometry (VLBI) facilities all over the globe start to implement compatible optics systems. Simultaneous dual/multi-frequency VLBI observations at mm wavelengths with international baselines are thus possible. In this paper, we present the results from the first successful simultaneous 22/43 GHz dual-frequency observation with KaVA(KVN and VERA array), including images and astrometric results. Our analysis shows that the newly implemented simultaneous receiving system has brought a significant extension of the coherence time of the 43 GHz visibility phases along the international baselines. The astrometric results obtained with KaVA are consistent with those obtained with the independent analysis of the KVN data. Our results thus confirm the good performance of the simultaneous receiving systems for the non-KVN stations. Future simultaneous observations with more global stations bring even higher sensitivity and micro-arcsecond level astrometric measurements of the targets.

Many motivated extensions of the Standard Model predict the existence of cosmic strings. Gravitational waves originating from the dynamics of the resulting cosmic string network have the ability to probe many otherwise inaccessible properties of the early universe. In this study we show how the spectrum of gravitational waves from a cosmic string network can be used to test the equation of state of the early universe prior to Big Bang Nucleosynthesis (BBN). We also demonstrate that current and planned gravitational wave detectors such as LIGO, LISA, DECIGO/BBO, and ET/CE have the potential to detect signals of a non-standard pre-BBN equation of state and evolution of the early universe (e.g., early non-standard matter domination or kination domination) or new degrees of freedom active in the early universe beyond the sensitivity of terrestrial collider experiments and cosmic microwave background measurements.

GeV flares from PSR B1259-63/LS 2883 were seen starting around 30 days after the two periastron passages in 2010 and 2014. The flares are clearly delayed compared to the occurrence of the X-ray and TeV flux peaks during the post-periastron disk crossing. While several attempts have been put forward to explain this phenomenon, the origin of these GeV flares remains a puzzle. Here we present a detailed analysis of the observational data taken by the Fermi and Swift observatories over the 2017 September periastron passage. For the first time, we find short-lived but powerful GeV flares on time scales of down to three hours. The onset of the GeV flaring period in 2017 is also delayed compared to those seen in 2011 and 2014. Supplemented by a re-analysis of previous data, we compare the Fermi/LAT, Swift/XRT and Swift/UVOT light curves in 2017 with those taken over the 2010 and 2014 periastrons, and difference in UVOT light curves are noted.

Supernova remnant (SNR)\,W28 is well known for its classic hadronic scenario, in which the TeV cosmic rays (CRs) released at the early stage of this mid-aged SNR are illuminating nearby molecular clouds (MCs). Overwhelming evidences have shown that the northeast of the SNR (W28-North) has already encountered with the MC clumps. Through this broken shell -- W28-North, we believe the CRs with energy down to $<$1\,GeV to be able to be injected into nearby MCs. To further testify this hadronic scenario, we first analyse the 9 years Fermi-LAT data in/around W28 with energy down to 0.3\,GeV. Our Fermi-LAT analysis display a 10-200 GeV skymap which spatially matches well with the known TeV sources -- HESS\,J1801-233 (W28-North), HESS\,J1800-240\,A,\,B\,\&\,C (240\,A\,B\,\&\,C). At low energy band, we has discovered a 0.5-1\,GeV blob located to the south of 240\,B\,\&\,C, and a low flux of 0.3-1\,GeV at 240\,A. A hadronic model is build to explain our analysis results and previous multi-wavelength observations of W28. Our model consists of three CR sources: The run-away CRs escaped from a strong shock; The leaked GeV CRs from the broken shell -- W28-North; And the local CR sea. Through modelling the SNR evolution, CR acceleration\,\&\,releasing, we have explained the GeV-TeV emission in/around SNR\,W28 (except for 240\,A) in one model. Both the damping of the magnetic waves by the neutrals and the decreased acceleration efficiency are taken into account in our model due to the mid-age of SNR W28.

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.

We consider the generic possibility that the Universe's energy budget includes some form of relativistic or semi-relativistic dark radiation (DR) with non-gravitational interactions with Standard Model (SM) particles. Such dark radiation may consist of SM singlets or a non-thermal, energetic component of neutrinos. If such DR is created at a relatively recent epoch, it can carry sufficient energy to leave a detectable imprint in experiments designed to search for very weakly interacting particles: dark matter and underground neutrino experiments. We analyze this possibility in some generality, assuming that the interactive dark radiation is sourced by late decays of an unstable particle, potentially a component of dark matter, and considering a variety of possible interactions between the dark radiation and SM particles. Concentrating on the sub-GeV energy region, we derive constraints on different forms of DR using the results of the most sensitive neutrino and dark matter direct detection experiments. In particular, for interacting dark radiation carrying a typical momentum of $\sim30$~MeV$/c$, both types of experiments provide competitive constraints. This study also demonstrates that non-standard sources of neutrino emission (e.g. via dark matter decay) are capable of creating a "neutrino floor" for dark matter direct detection that is closer to current bounds than is expected from standard neutrino sources.

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.

We study UV complete theories where the Standard Model (SM) gauge group is extended with a new abelian $U(1)$, and the field content is augmented by an arbitrary number of scalar and fermion SM singlets, potentially including dark matter (DM) candidates. Considerations such as classical and quantum gauge invariance of the full theory and S-matrix unitarity, not applicable within a simplified model approach, are shown to have significant phenomenological consequences. The lack of gauge anomalies leads to compact relations among the $U(1)$ fermion charges, and puts a lower bound on the number of dark fermions. Contrary to naive expectations, the DM annihilation to Zh is found to be p-wave suppressed, as hinted by perturbative unitarity of S-matrix, with dramatic implications for DM thermal relic density and indirect searches. Within this framework, the interplay between dark matter, new vector boson and Higgs physics is rather natural and generic.

Cosmic strings are generic cosmological predictions of many extensions of the Standard Model of particle physics, such as a $U(1)^\prime$ symmetry breaking phase transition in the early universe or remnants of superstring theory. Unlike other topological defects, cosmic strings can reach a scaling regime that maintains a small fixed fraction of the total energy density of the universe from a very early epoch until today. If present, they will oscillate and generate gravitational waves with a frequency spectrum that imprints the dominant sources of total cosmic energy density throughout the history of the universe. We demonstrate that current and future gravitational wave detectors, such as LIGO and LISA, could be capable of measuring the frequency spectrum of gravitational waves from cosmic strings and discerning the energy composition of the universe at times well before primordial nucleosynthesis and the cosmic microwave background where standard cosmology has yet to be tested. This work establishes a benchmark case that gravitational waves may provide an unprecedented, powerful tool for probing the evolutionary history of the very early universe.

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

The standard paradigm of collisionless cold dark matter is in tension with measurements on large scales. In particular, the best fit values of the Hubble rate $H_0$ and the matter density perturbation $\sigma_8$ inferred from the cosmic microwave background seem inconsistent with the results from direct measurements. We show that both problems can be solved in a framework in which dark matter consists of two distinct components, a dominant component and a subdominant component. The primary component is cold and collisionless. The secondary component is also cold, but interacts strongly with dark radiation, which itself forms a tightly coupled fluid. The growth of density perturbations in the subdominant component is inhibited by dark acoustic oscillations due to its coupling to the dark radiation, solving the $\sigma_8$ problem, while the presence of tightly coupled dark radiation ameliorates the $H_0$ problem. The subdominant component of dark matter and dark radiation continue to remain in thermal equilibrium until late times, inhibiting the formation of a dark disk. We present an example of a simple model that naturally realizes this scenario in which both constituents of dark matter are thermal WIMPs. Our scenario can be tested by future stage-IV experiments designed to probe the CMB and large scale structure.

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.

We present here the results of the first part of the VLBI Ecliptic Plane Survey (VEPS) program. The goal of the program is to find all compact sources within $7.5^\circ$ of the ecliptic plane which are suitable as calibrators for anticipated phase referencing observations of spacecraft and determine their positions with accuracy at the 1.5~nrad level. We run the program in two modes: the search mode and the refining mode. In the search mode, a complete sample of all sources brighter than 50 mJy at 5 GHz listed in the Parkes-MIT-NRAO (PMN) and Green Bank 6~cm (GB6) catalogs, except those previously detected with VLBI, is observed. In the refining mode, the positions of all ecliptic plane sources, including those found in the search mode, are improved. By October 2016, thirteen 24-hr sessions that targeted all sources brighter than 100~mJy have been observed and analyzed. Among 3320 observed target sources, 555 objects have been detected. We also conducted a number of follow-up VLBI experiments in the refining mode and improved the positions of 249 ecliptic plane sources.

Runaway stars form bow shocks by ploughing through the interstellar medium at supersonic speeds and are promising sources of non-thermal emission of photons. One of these objects has been found to emit non-thermal radiation in the radio band. This triggered the development of theoretical models predicting non-thermal photons from radio up to very-high-energy (VHE, E $\geq 0.1$TeV) gamma rays. Subsequently, one bow shock was also detected in X-ray observations. However, the data did not allow discrimination between a hot thermal and a non-thermal origin. Further observations of different candidates at X-ray energies showed no evidence for emission at the position of the bow shocks either. A systematic search in the \textit{Fermi}-LAT energy regime resulted in flux upper limits for 27 candidates listed in the E-BOSS catalogue. Here we perform the first systematic search for VHE gamma-ray emission from bow shocks of runaway stars. Using all available archival H.E.S.S. data we search for very-high-energy gamma-ray emission at the positions of bow shock candidates listed in the second E-BOSS catalogue release. Out of the 73 bow shock candidates in this catalogue, 32 have been observed with H.E.S.S. None of the observed 32 bow shock candidates in this population study show significant emission in the H.E.S.S. energy range. Therefore, flux upper limits are calculated in five energy bins and the fraction of the kinetic wind power that is converted into VHE gamma rays is constrained. Emission from stellar bow shocks is not detected in the energy range between 0.14 and 18 TeV. The resulting upper limits constrain the level of VHE gamma-ray emission from these objects down to 0.1-1 % of the kinetic wind energy.

Cosmological observations are becoming increasingly sensitive to the effects of light particles in the form of dark radiation (DR) at the time of recombination. The conventional observable of effective neutrino number, $N_{\rm eff}$, is insufficient for probing generic, interacting models of DR. In this work, we perform likelihood analyses which allow both free-streaming effective neutrinos (parametrized by $N_{\rm eff}$) and interacting effective neutrinos (parametrized by $N_{\rm fld}$). We motivate an alternative parametrization of DR in terms of $N_{\rm tot}$ (total effective number of neutrinos) and $f_{\rm fs}$ (the fraction of effective neutrinos which are free-streaming), which is less degenerate than using $N_{\rm eff}$ and $N_{\rm fld}$. Using the Planck 2015 likelihoods in conjunction with measurements of baryon acoustic oscillations (BAO), we find constraints on the total amount of beyond the Standard Model effective neutrinos (both free-streaming and interacting) of $\Delta N_{\rm tot} < 0.39$ at 2$\sigma$. In addition, we consider the possibility that this scenario alleviates the tensions between early-time and late-time cosmological observations, in particular the measurements of $\sigma_8$ (the amplitude of matter power fluctuations at 8$h^{-1}$ Mpc), finding a mild preference for interactions among light species. We further forecast the sensitivities of a variety of future experiments, including Advanced ACTPol (a representative CMB Stage-III experiment), CMB Stage-IV, and the Euclid satellite. This study is relevant for probing non-standard neutrino physics as well as a wide variety of new particle physics models beyond the Standard Model that involve dark radiation.

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.

Supernova remnants exhibit shock fronts (shells) that can accelerate charged particles up to very high energies. In the past decade, measurements of a handful of shell-type supernova remnants in very-high-energy gamma rays have provided unique insights into the acceleration process. Among those objects, RX$~$J1713.7-3946 (also known as G347.3-0.5) has the largest surface brightness, allowing us in the past to perform the most comprehensive study of morphology and spatially resolved spectra of any such very-high-energy gamma-ray source. Here we present extensive new H.E.S.S. measurements of RX$~$J1713.7-3946, almost doubling the observation time compared to our previous publication. Combined with new improved analysis tools, the previous sensitivity is more than doubled. The H.E.S.S. angular resolution of $0.048^\circ$ ($0.036^\circ$ above 2 TeV) is unprecedented in gamma-ray astronomy and probes physical scales of 0.8 (0.6) parsec at the remnant's location. The new H.E.S.S. image of RX$~$J1713.7-3946 allows us to reveal clear morphological differences between X-rays and gamma rays. In particular, for the outer edge of the brightest shell region, we find the first ever indication for particles in the process of leaving the acceleration shock region. By studying the broadband energy spectrum, we furthermore extract properties of the parent particle populations, providing new input to the discussion of the leptonic or hadronic nature of the gamma-ray emission mechanism.

The addition of a 28 m Cherenkov telescope (CT5) to the H.E.S.S. array extended the experiment's sensitivity to lower energies. The lowest energy threshold is obtained using monoscopic analysis of data taken with CT5, providing access to gamma-ray energies below 100 GeV. Such an extension of the instrument's energy range is particularly beneficial for studies of Active Galactic Nuclei with soft spectra, as expected for those at a redshift > 0.5. The high-frequency peaked BL Lac objects PKS 2155-304 (z = 0.116) and PG 1553+113 (0.43 < z < 0.58) are among the brightest objects in the gamma-ray sky, both showing clear signatures of gamma-ray absorption at E > 100 GeV interpreted as being due to interactions with the extragalactic background light (EBL). Multiple observational campaigns of PKS 2155-304 and PG 1553+113 were conducted during 2013 and 2014 using the full H.E.S.S. II instrument. A monoscopic analysis of the data taken with the new CT5 telescope was developed along with an investigation into the systematic uncertainties on the spectral parameters. The energy spectra of PKS 2155-304 and PG 1553+113 were reconstructed down to energies of 80 GeV for PKS 2155-304, which transits near zenith, and 110 GeV for the more northern PG 1553+113. The measured spectra, well fitted in both cases by a log-parabola spectral model (with a 5.0 sigma statistical preference for non-zero curvature for PKS 2155-304 and 4.5 sigma for PG 1553+113), were found consistent with spectra derived from contemporaneous Fermi-LAT data, indicating a sharp break in the observed spectra of both sources at E ~ 100 GeV. When corrected for EBL absorption, the intrinsic H.E.S.S. II mono and Fermi-LAT spectrum of PKS 2155-304 was found to show significant curvature. For PG 1553+113, however, no significant detection of curvature in the intrinsic spectrum could be found within statistical and systematic uncertainties.

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

A robust mechanism was recently proposed in which thermal freeze-out of WIMPs can provide a unified origin of dark matter and baryon abundances in our universe. We point out that this WIMP-triggered baryogenesis mechanism can exhibit a rich collider phenomenology and be tested at the current and near-future experiments at LHC, even in the case where the WIMPs are completely devoid of SM gauge and higgs portal interactions, as may be motivated by the persistent null results of WIMP dark matter searches. We catalogue a rich array of LHC signatures robustly present in such a scenario. In particular, the simplest such implementation can already offer a very clean signal of a TeV-scale resonance that decays to diphotons with a cross section that can easily be within the reach of the current and near-future LHC runs in the region of parameter space that leads to a successful baryogenesis. Other characteristic signatures include the production of multi-bottom and/or multi-top quarks, promptly or displaced. An even more exotic possibility is the production of two separate sets of isolated emerging jets connected by a charged track, which may require new dedicated studies. Finally, di-nucleon decay can also provide a powerful probe of the mechanism.

Aims: Following the detection of the fast radio burst FRB150418 by the SUPERB project at the Parkes radio telescope, we aim to search for very-high energy gamma-ray afterglow emission. Methods: Follow-up observations in the very-high energy gamma-ray domain were obtained with the H.E.S.S. imaging atmospheric Cherenkov telescope system within 14.5 hours of the radio burst. Results: The obtained 1.4 hours of gamma-ray observations are presented and discussed. At the 99 % C.L. we obtained an integral upper limit on the gamma-ray flux of (E>350 GeV) < 1.33 x 10^-8 m^-2s^-1. Differential flux upper limits as function of the photon energy were derived and used to constrain the intrinsic high-energy afterglow emission of FRB 150418. Conclusions: No hints for high-energy afterglow emission of FRB 150418 were found. Taking absorption on the extragalactic background light into account and assuming a distance of z = 0.492 based on radio and optical counterpart studies and consistent with the FRB dispersion, we constrain the gamma-ray luminosity at 1 TeV to L < 5.1 x 10^47 erg/s at 99% C.L.

Most of the gamma-ray emitting active galactic nuclei (AGN) are blazars, although there is still a small fraction of non-blazar AGN in the Fermi/LAT catalog. Among these misaligned gamma-ray-emitting AGN, a few can be classified as Compact Symmetric Objects (CSOs). In contrast to blazars in which gamma-ray emission is generally thought to originate from highly beamed relativistic jets, the source of gamma-ray emission in unbeamed CSOs remains an open question. The rarity of the gamma-ray emitting CSOs is a mystery as well. Here we present the radio properties of the gamma-ray CSO candidate 2234+282.

The quasar 3C~286 is one of two compact steep spectrum sources detected by the {\it Fermi}/LAT. Here, we investigate the radio properties of the parsec(pc)-scale jet and its (possible) association with the $\gamma$-ray emission in 3C~286. The Very Long Baseline Interferometry (VLBI) images at various frequencies reveal a one-sided core--jet structure extending to the southwest at a projected distance of $\sim$1 kpc. The component at the jet base showing an inverted spectrum is identified as the core, with a mean brightness temperature of $2.8\times 10^{9}$~K. The jet bends at about 600 pc (in projection) away from the core, from a position angle of $-135^\circ$ to $-115^\circ$. Based on the available VLBI data, we inferred the proper motion speed of the inner jet as $0.013 \pm 0.011$ mas yr$^{-1}$ ($\beta_{\rm app} = 0.6 \pm 0.5$), corresponding to a jet speed of about $0.5\,c$ at an inclination angle of $48^\circ$ between the jet and the line of sight of the observer. The brightness temperature, jet speed and Lorentz factor are much lower than those of $\gamma$-ray-emitting blazars, implying that the pc-scale jet in 3C~286 is mildly relativistic. Unlike blazars in which $\gamma$-ray emission is in general thought to originate from the beamed innermost jet, the location and mechanism of $\gamma$-ray emission in 3C~286 may be different as indicated by the current radio data. Multi-band spectrum fitting may offer a complementary diagnostic clue of the $\gamma$-ray production mechanism in this source.

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

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

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 DUNE/LBNF program aims to address key questions in neutrino physics and astroparticle physics. Realizing DUNE's potential to reconstruct low-energy particles in the 10-100 MeV energy range will bring significant benefits for all DUNE's science goals. In neutrino physics, low-energy sensitivity will improve neutrino energy reconstruction in the GeV range relevant for the kinematics of DUNE's long-baseline oscillation program. In astroparticle physics, low-energy capabilities will make DUNE's far detectors the world's best apparatus for studying the electron-neutrino flux from a supernova. This will open a new window to unrivaled studies of the dynamics and neutronization of a star's central core in real time, the potential discovery of the neutrino mass hierarchy, provide new sensitivity to physics beyond the Standard Model, and evidence of neutrino quantum-coherence effects. The same capabilities will also provide new sensitivity to `boosted dark matter' models that are not observable in traditional direct dark matter detectors.

A bright optical flare was detected in the high-redshift ($z=2.133$) quasar CGRaBS J0809+5341 on 2014 April 13. The absolute magnitude of the object reached $-30.0$ during the flare, making it the brightest one (in flaring stage) among all known quasars so far. The 15 GHz flux density of CGRaBS J0809+5341 monitored in the period from 2008 to 2016 also reached its peak at the same time. To reveal any structural change possibly associated with the flare in the innermost radio structure of the quasar, we conducted a pilot very long baseline interferometry (VLBI) observation of CGRaBS J0809+5341 using the European VLBI Network (EVN) at 5 GHz on 2014 November 18, about seven months after the prominent optical flare. Three epochs of follow-up KaVA (Korean VLBI Network and VLBI Exploration of Radio Astrometry Array) observations were carried out at 22 and 43 GHz frequencies from 2015 February 25 to June 4, with the intention of exploring a possibly emerging new radio jet component associated with the optical flare. However, these high-resolution VLBI observations revealed only the milliarcsecond-scale compact "core" that was known in the quasar from earlier VLBI images, and showed no sign of any extended jet structure. Neither the size, nor the flux density of the "core" changed considerably after the flare according to our VLBI monitoring. The results suggest that any putative radio ejecta associated with the major optical and radio flare could not yet be separated from the "core" component, or the newly-born jet was short-lived.

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.

PKS 0202$+$149 is a low-power radio source with blazar-like $\gamma$-ray AGN characteristics and its properties and classification are investigated in relation to its $\gamma$-ray characteristics. It shows a hint of low frequency turnover at about 200 MHz. Radio imaging data of 0202$+$149 at different frequencies show differing morphologies on both kilo-parsec (kpc) and parsec (pc) scales. The overall source shows a triple structure of core and double lobes with a total projected size of $\sim$1.3 kpc. The compact source structure of 0202+149 is reminiscent of a compact steep spectrum (CSS) source. At pc scales a core-jet structure extends $\sim$25 pc (in projection) at a position angle perpendicular to the kpc-scale structure. The curved pc-scale structure with jet and inner lobe suggests that the CSS nuclear activity has recently re-started although its power has been decreasing, while the kpc-scale lobes are relics of earlier activity. A maximum apparent superluminal motion of $\sim16\,c$ is detected in the jet components, indicating a highly relativistic jet flow. The brightness temperature of the core is lower than the average value found for highly-beamed $\gamma$-ray AGN, indicating a lower radio power and a relatively lower Doppler boosting factor. The CSS radio classification indicates that blazar-like $\gamma$-ray properties can also be manifested in low-power CSS radio sources with appropriate jet and beaming properties.

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.

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.

In a broad class of theories, the relic abundance of dark matter is determined by interactions internal to a thermalized dark sector, with no direct involvement of the Standard Model (SM). We point out that these theories raise an immediate cosmological question: how was the dark sector initially populated in the early universe? Motivated in part by the difficulty of accommodating large amounts of entropy carried in dark radiation with cosmic microwave background measurements of the effective number of relativistic species at recombination, $N_{\mathrm{eff}}$, we aim to establish which admissible cosmological histories can populate a thermal dark sector that never reaches thermal equilibrium with the SM. The minimal cosmological origin for such a dark sector is asymmetric reheating, when the same mechanism that populates the SM in the early universe also populates the dark sector at a lower temperature. Here we demonstrate that the resulting inevitable inflaton-mediated scattering between the dark sector and the SM can wash out a would-be temperature asymmetry, and establish the regions of parameter space where temperature asymmetries can be generated in minimal reheating scenarios. Thus obtaining a temperature asymmetry of a given size either restricts possible inflaton masses and couplings or necessitates a non-minimal cosmology for one or both sectors. As a side benefit, we develop techniques for evaluating collision terms in the relativistic Boltzmann equation when the full dependence on Bose-Einstein or Fermi-Dirac phase space distributions must be retained, and present several new results on relativistic thermal averages in an appendix.

The Supernova Remnant (SNR) HESS J1731-347 displays strong non-thermal TeV gamma-ray and X-ray emission, thus the object is at present time accelerating particles to very high energies. A distinctive feature of this young SNR is the nearby (~30 pc in projection) extended source HESS J1729-345, which is currently unidentified but is in spatial projection coinciding with known molecular clouds (MC). We model the SNR evolution to explore if the TeV emission from HESS J1729-345 can be explained as emission from runaway hadronic cosmic rays (CRs) that are illuminating these MCs. The observational data of HESS J1729-345 and HESS J1731-347 can be reproduced using core-collapse SN models for HESS J1731-347. Starting with different progenitor stars and their pre-supernova environment, we model potential SNR evolution histories along with the CR acceleration in the SNR and the diffusion of the CRs. A simplified 3-dimensional structure of the MCs is introduced based on 12CO data, adopting a distance of 3.2 kpc to the source. A Monte Carlo-based diffusion model for the escaping CRs is developed to deal with the inhomogeneous environment. The fast SNR forward shock speed as implied from the X-ray data can easily be explained when employing scenarios with progenitor star masses between 20 and 25 solar masses, where the SNR shock is still expanding inside the main sequence (MS)-bubble at present time. The TeV spectrum of HESS J1729-345 is satisfactorily fitted by the emission from the highest-energy CRs that have escaped the SNR, using a standard galactic CR diffusion coefficient in the inter-clump medium. The TeV image of HESS J1729-345 can be explained with a reasonable 3-dimensional structure of MCs. The TeV emission from the SNR itself is dominated by leptonic emission in this model. We also explore scenarios where the shock is starting to encounter the dense MS progenitor wind bubble shell.

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.

We consider theories where dark matter is composed of a thermal relic of weak scale mass, whose couplings to the Standard Model (SM) are however too small to give rise to the observed abundance. Instead, the abundance is set by annihilation to light hidden sector states that carry no charges under the SM gauge interactions. In such a scenario the constraints from direct and indirect detection, and from collider searches for dark matter, can easily be satisfied. The masses of such light hidden states can be protected by symmetry if they are Nambu-Goldstone bosons, fermions, or gauge bosons. These states can then contribute to the cosmic energy density as dark radiation, leading to observable signals in the cosmic microwave background (CMB). Furthermore, depending on whether or not the light hidden sector states self-interact, the fraction of the total energy density that free-streams is either decreased or increased, leading to characteristic effects on both the scalar and tensor components of the CMB anisotropy that allows these two cases to be distinguished. The magnitude of these signals depends on the number of light degrees of freedom in the hidden sector, and on the temperature at which it kinetically decouples from the SM. We consider a simple model that realizes this scenario, based on a framework in which the SM and hidden sector are initially in thermal equilibrium through the Higgs portal, and show that the resulting signals are compatible with recent Planck results, while large enough to be detected in upcoming experiments such as CMBPol and CMB Stage-IV. Invisible decays of the Higgs into hidden sector states at colliders can offer a complementary probe of this model.

We present the particle-in-cell (PIC) simulation results of the interaction of a high-energy lepton plasma flow with background electron-proton plasma and focus on the acceleration processes of the protons. It is found that the acceleration follows a two-stage processes. In the first stage, protons are accelerated transversely (perpendicular to the lepton flow) by the turbulent magnetic field "islands" generated via the strong Weibel-type instabilities. The accelerated protons shows a perfect inverse-power energy spectrum. As the interaction continues, a shockwave structure forms and the protons in front of the shockwave are reflected at twice of the shock speed, resulting in a quasi-monoenergetic peak located near 200MeV under the simulation parameters.

We present a new multi-component dark matter model with a novel experimental signature that mimics neutral current interactions at neutrino detectors. In our model, the dark matter is composed of two particles, a heavier dominant component that annihilates to produce a boosted lighter component that we refer to as boosted dark matter. The lighter component is relativistic and scatters off electrons in neutrino experiments to produce Cherenkov light. This model combines the indirect detection of the dominant component with the direct detection of the boosted dark matter. Directionality can be used to distinguish the dark matter signal from the atmospheric neutrino background. We discuss the viable region of parameter space in current and future experiments.

It was recently proposed that weakly interacting massive particles (WIMP) may provide new ways of generating the observed baryon asymmetry in the early universe, as well as addressing the cosmic coincidence between dark matter and baryon abundances. This suggests a new possible connection between weak scale new particle physics and modern cosmology. This review summarizes the general ideas and simple model examples of the two recently proposed WIMP baryogenesis mechanisms: baryogenesis from WIMP dark matter annihilation during thermal freezeout, and baryogenesis from metastable WIMP decay after thermal freezeout. This letter also discusses the interesting phenomenology of these models, in particular the experimental signals that can be probed in the intensity frontier experiments and the Large Hadron Collider (LHC) experiments.

The generation of the asymmetric cosmic baryon abundance requires a departure from thermal equilibrium in the early universe. In a large class of baryogenesis models, the baryon asymmetry results from the out-of-equilibrium decay of a new, massive particle. We highlight that in the interesting scenario where this particle has a weak scale mass, this out-of-equilibrium condition requires a proper decay length larger than O(1) mm. Such new fields are within reach of the LHC, at which they can be pair produced leaving a distinctive, displaced-vertex signature. This scenario is realized in the recently proposed mechanism of baryogenesis where the baryon asymmetry is produced through the freeze-out and subsequent decay of a meta-stable weakly interacting massive particle ("WIMP baryogenesis"). In analogy to missing energy searches for WIMP dark matter, the LHC is an excellent probe of these new long-lived particles responsible for baryogenesis via the low-background displaced vertex channel. In our paper, we estimate the limits on simplified models inspired by WIMP baryogenesis from two of the most sensitive collider searches by CMS and ATLAS with 8 TeV LHC data. We also estimate the LHC reach at 13 TeV using current strategies, and demonstrate that O(10) improvements in cross-section limits can be achieved by requiring two displaced vertices while lowering kinematic thresholds. For meta-stable WIMPs produced through electroweak interactions, the high luminosity LHC is sensitive to masses up to 2.5 TeV for lifetimes around 1 cm, while for singlets pair-produced through the off-shell-Higgs portal, the LHC is sensitive to production cross sections of O(50) ab for benchmark masses around 150 GeV. Our analysis and proposals also generally apply to displaced vertex signatures from other new physics such as hidden valley models, twin Higgs models and displaced supersymmetry.

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

Puppis A is an interesting ~4 kyr-old supernova remnant (SNR) that shows strong evidence of interaction between the forward shock and a molecular cloud. It has been studied in detail from radio frequencies to high-energy (HE, 0.1-100 GeV) gamma-rays. An analysis of the Fermi-LAT data has shown an extended HE gamma-ray emission with a 0.2-100 GeV spectrum exhibiting no significant deviation from a power law, unlike most of the GeV-emitting SNRs known to be interacting with molecular clouds. This makes it a promising target for imaging atmospheric Cherenkov telescopes (IACTs) to probe the gamma-ray emission above 100 GeV. Very-high-energy (VHE, E >= 0.1 TeV) gamma-ray emission from Puppis A is for the first time searched for with the High Energy Stereoscopic System (H.E.S.S.). The analysis of the H.E.S.S. data does not reveal any significant emission towards Puppis A. The derived upper limits on the differential photon flux imply that its broadband gamma-ray spectrum must exhibit a spectral break or cutoff. By combining Fermi-LAT and H.E.S.S. measurements, the 99% confidence level upper limits on such a cutoff are found to be 450 and 280 GeV, assuming a power law with a simple exponential and a sub-exponential cutoff, respectively. It is concluded that none of the standard limitations (age, size, radiative losses) on the particle acceleration mechanism, assumed to be still on-going at present, can explain the lack of VHE signal. The scenario in which particle acceleration has ceased some time ago is considered as an alternative explanation. The HE/VHE spectrum of Puppis A could then exhibit a break of non-radiative origin, (as observed in several other interacting SNRs, albeit at somewhat higher energies) owing to the interaction with dense and neutral material in particular towards the northeastern region.

An annihilation signal of dark matter is searched for from the central region of the Milky Way. Data acquired in dedicated ON/OFF observations of the Galactic center region with H.E.S.S. are analyzed for this purpose. No significant signal is found in a total of $\sim 9$ h of ON/OFF observations. Upper limits on the velocity averaged cross section, $<\sigma v >$, for the annihilation of dark matter particles with masses in the range of $\sim 300$ GeV to $\sim 10$ TeV are derived. In contrast to previous constraints derived from observations of the Galactic center region, the constraints that are derived here apply also under the assumption of a central core of constant dark matter density around the center of the Galaxy. Values of $<\sigma v >$ that are larger than $3\cdot 10^{-24}\:\mathrm{cm^3/s}$ are excluded for dark matter particles with masses between $\sim 1$ and $\sim 4$ TeV at 95% CL if the radius of the central dark matter density core does not exceed $500$ pc. This is the strongest constraint that is derived on $<\sigma v>$ for annihilating TeV mass dark matter without the assumption of a centrally cusped dark matter density distribution in the search region.

The Large Magellanic Cloud, a satellite galaxy of the Milky Way, has been observed with the High Energy Stereoscopic System (H.E.S.S.) above an energy of 100 billion electron volts for a deep exposure of 210 hours. Three sources of different types were detected: the pulsar wind nebula of the most energetic pulsar known N 157B, the radio-loud supernova remnant N 132D and the largest non-thermal X-ray shell - the superbubble 30 Dor C. The unique object SN 1987A is, surprisingly, not detected, which constrains the theoretical framework of particle acceleration in very young supernova remnants. These detections reveal the most energetic tip of a gamma-ray source population in an external galaxy, and provide via 30 Dor C the unambiguous detection of gamma-ray emission from a superbubble.

Very high energy (VHE, $E>$100 GeV) $\gamma$-ray flaring activity of the high-frequency peaked BL Lac object \pg\ has been detected by the \hess\ telescopes. The flux of the source increased by a factor of 3 during the nights of 2012 April 26 and 27 with respect to the archival measurements with hint of intra-night variability. No counterpart of this event has been detected in the \fla\ data. This pattern is consistent with VHE $\gamma$ ray flaring being caused by the injection of ultrarelativistic particles, emitting $\gamma$ rays at the highest energies. The dataset offers a unique opportunity to constrain the redshift of this source at \bestz\ using a novel method based on Bayesian statistics. The indication of intra-night variability is used to introduce a novel method to probe for a possible Lorentz Invariance Violation (LIV), and to set limits on the energy scale at which Quantum Gravity (QG) effects causing LIV may arise. For the subluminal case, the derived limits are $\textrm{E}_{\rm QG,1}>4.10\times 10^{17}$ GeV and $\textrm{E}_{\rm QG,2}>2.10\times 10^{10}$ GeV for linear and quadratic LIV effects, respectively.

We initiate the study of novel thermal dark matter (DM) scenarios where present-day annihilation of DM in the galactic center produces boosted stable particles in the dark sector. These stable particles are typically a subdominant DM component, but because they are produced with a large Lorentz boost in this process, they can be detected in large volume terrestrial experiments via neutral-current-like interactions with electrons or nuclei. This novel DM signal thus combines the production mechanism associated with indirect detection experiments (i.e. galactic DM annihilation) with the detection mechanism associated with direct detection experiments (i.e. DM scattering off terrestrial targets). Such processes are generically present in multi-component DM scenarios or those with non-minimal DM stabilization symmetries. As a proof of concept, we present a model of two-component thermal relic DM, where the dominant heavy DM species has no tree-level interactions with the standard model and thus largely evades direct and indirect DM bounds. Instead, its thermal relic abundance is set by annihilation into a subdominant lighter DM species, and the latter can be detected in the boosted channel via the same annihilation process occurring today. Especially for dark sector masses in the 10 MeV-10 GeV range, the most promising signals are electron scattering events pointing toward the galactic center. These can be detected in experiments designed for neutrino physics or proton decay, in particular Super-K and its upgrade Hyper-K, as well as the PINGU/MICA extensions of IceCube. This boosted DM phenomenon highlights the distinctive signatures possible from non-minimal dark sectors.