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

The observed accelerated expansion of the Universe may be explained by dark energy or the breakdown of general relativity (GR) on cosmological scales. When the latter case, a modified gravity scenario, is considered, it is often assumed that the background evolution is the same as the $\Lambda$CDM model but the density perturbation evolves differently. In this paper, we investigate more general classes of modified gravity, where both the background and perturbation evolutions are deviated from those in the $\Lambda$CDM model. We introduce two phase diagrams, $\alpha{\rm-}f\sigma _8$ and $H{\rm-}f\sigma _8$ diagrams; $H$ is the expansion rate, $f\sigma_8$ is a combination of the growth rate of the Universe and the normalization of the density fluctuation which is directly constrained by redshift-space distortions, and $\alpha$ is a parameter which characterizes the deviation of gravity from GR and can be probed by gravitational lensing. We consider several specific examples of Horndeski's theory, which is a general scalar-tensor theory, and demonstrate how deviations from the $\Lambda$CDM model appears in the $\alpha{\rm-}f\sigma _8$ and $H{\rm-}f\sigma _8$ diagrams. The predicted deviations will be useful for future large-scale structure observations to exclude some of the modified gravity models.

We apply the independent component analysis (ICA) to the real data from a gravitational wave detector for the first time. Specifically we use the iKAGRA data taken in April 2016, and calculate the correlations between the gravitational wave strain channel and 35 physical environmental channels. Using a couple of seismic channels which are found to be strongly correlated with the strain, we perform ICA. Injecting a sinusoidal continuous signal in the strain channel, we find that ICA recovers correct parameters with enhanced signal-to-noise ratio, which demonstrates usefulness of this method. Among the two implementations of ICA used here, we find the correlation method yields the optimal result for the case environmental noises act on the strain channel linearly.

We report the first evidence for high-mass star formation triggered by collisions of molecular clouds in M33. Using the Atacama Large Millimeter/submillimeter Array, we spatially resolved filamentary structures of giant molecular cloud 37 in M33 using $^{12}$CO($J$ = 2-1), $^{13}$CO($J$ = 2-1), and C$^{18}$O($J$ = 2-1) line emission at a spatial resolution of $\sim$2 pc. There are two individual molecular clouds with a systematic velocity difference of $\sim$6 km s$^{-1}$. Three continuum sources representing up to $\sim$10 high-mass stars with the spectral types of B0V-O7.5V are embedded within the densest parts of molecular clouds bright in the C$^{18}$O($J$ = 2-1) line emission. The two molecular clouds show a complementary spatial distribution with a spatial displacement of $\sim$6.2 pc, and show a V-shaped structure in the position-velocity diagram. These observational features traced by CO and its isotopes are consistent with those in high-mass star-forming regions created by cloud-cloud collisions in the Galactic and Magellanic Cloud HII regions. Our new finding in M33 indicates that the cloud-cloud collision is a promising process to trigger high-mass star formation in the Local Group.

We study the stochastic gravitational wave (GW) background induced by the primordial scalar perturbation with the spectrum having a lognormal peak of width $\Delta$ at $k=k_*$. We derive an analytical formula for the GW spectrum $\Omega_\text{GW}$ for both narrow ($\Delta\ll1$) and broad ($\Delta\gtrsim1$) peaks. In the narrow-peak case, the spectrum has a double peak feature with the sharper peak at $k= 2k_*/\sqrt{3}$. On the infrared (IR) side of the spectrum, we find power-law behavior with a break at $k=k_b$ in the power-law index where it chages from $k^3$ on the far IR side to $k^2$ on the near IR side. We find the ratio of the break frequency to the peak frequency is determined by $\Delta$ as $f_b/f_p\approx\sqrt{3}\Delta$, where $f_b$ and $f_p$ are the break and peak frequencies, respectively. In the broad-peak case, we find the GW spectrum also has a lognormal peak at $k=k_*$ but with a smaller width of $\Delta/\sqrt2$. Using these derived analytic formulae, we also present expressions for the maximum values of $\Omega_\text{GW}$ for both narrow and broad cases. Our results will provide a useful tool in searching for the induced GW signals in the coming decades.

The scalar perturbation induced gravitational waves are a probe of the primordial density perturbation spectrum on small scales. In this paper, we show that they can also probe the thermal history of the universe. We assume the universe underwent a stage with a constant equation of state parameter $w$, followed by the radiation-dominated stage of the conventional bigbang universe. We find that the infrared slope of the power spectrum of the induced stochastic gravitational wave background for decelerating cosmologies is related to the equation of state of the universe. Furthermore, the induced gravitational wave spectrum has in general a broken power-law shape around the scale of reheating. Interestingly, below the threshold $w=0$ of the equation of state parameter, the broken power-law presents a peak for a dirac delta peak in the scalar spectrum. For a finite width peak, the threshold changes to $w=-1/15$ depending on the value of the width. In some cases, such a broken power-law gravitational wave spectrum may degenerate to the spectrum from other sources like phase transitions or global cosmic strings.

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

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

Aims: We use new deep XMM-Newton observations of the northern disk of M 31 to trace the hot interstellar medium (ISM) in unprecedented detail and to characterise the physical properties of the X-ray emitting plasmas. Methods: We used all XMM-Newton data up to and including our new observations to produce the most detailed image yet of the hot ISM plasma in a grand design spiral galaxy such as our own. We compared the X-ray morphology to multi-wavelength studies in the literature to set it in the context of the multi-phase ISM. We performed spectral analyses on the extended emission using our new observations as they offer sufficient depth and count statistics to constrain the plasma properties. Data from the Panchromatic Hubble Andromeda Treasury were used to estimate the energy injected by massive stars and their supernovae. We compared these results to the hot gas properties. Results: The brightest emission regions were found to be correlated with populations of massive stars, notably in the 10 kpc star-forming ring. The plasma temperatures in the ring regions are ~0.2 keV up to ~0.6 keV. We suggest this emission is hot ISM heated in massive stellar clusters and superbubbles. We derived X-ray luminosities, densities, and pressures for the gas in each region. We also found large extended emission filling low density gaps in the dust morphology of the northern disk, notably between the 5 kpc and 10 kpc star-forming rings. We propose that the hot gas was heated and expelled into the gaps by the populations of massive stars in the rings. Conclusions. It is clear that the massive stellar populations are responsible for heating the ISM to X-ray emitting temperatures, filling their surroundings, and possibly driving the hot gas into the low density regions. Overall, the morphology and spectra of the hot gas in the northern disk of M 31 is similar to other galaxy disks.

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

Primordial black holes (PBHs) are a viable candidate for dark matter if the PBH masses are in the currently unconstrained "sublunar" mass range. We revisit the possibility that PBHs were produced by nucleation of false vacuum bubbles during inflation. We show that this scenario can produce a population of PBHs that simultaneously accounts for all dark matter, explains the candidate event in Subaru Hyper Suprime-Cam (HSC) data, and contains both heavy black holes as observed by LIGO and very heavy seeds of super-massive black holes. We demonstrate with numerical studies that future observations of HSC, as well as other optical surveys, such as LSST, will be able to provide a definitive test for this generic PBH formation mechanism if it is the dominant source of dark matter.

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

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

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

Here we present a radio continuum study based on new and archival data from the Australia Telescope Compact Array towards N 103B, a young (<=1000 yrs) spectroscopically confirmed type Ia SNR in the Large Magellanic Cloud (LMC) and proposed to have originated from a single degenerate progenitor. The radio morphology of this SNR is asymmetrical with two bright regions towards the north-west and south-west of the central location as defined by radio emission. N 103B identified features include: a radio spectral index of -0.75+-0.01 (consistent with other young type Ia SNRs in the Galaxy); a bulk SNR expansion rate as in X-rays; morphology and polarised electrical field vector measurements where we note radial polarisation peak towards the north-west of the remnant at both 5500 and 9000 MHz. The spectrum is concave-up and the most likely reason is the non-linear diffusive shock acceleration effects or presence of two different populations of ultra-relativistic electrons. We also note unpolarized clumps near the south-west region which is in agreement with this above scenario. We derive a typical magnetic field strength for N 103B, of 16.4 microG for an average rotation measurement of 200 rad m^-2. However, we estimate the equipartition field to be of the order of ~235 microG with an estimated minimum energy of Emin=6.3*10^48 erg. The close (~0.5 degree) proximity of N 103B to the LMC mid-plane indicates that an early encounter with dense interstellar medium may have set an important constrain on SNR evolution. Finally, we compare features of N 103B, to six other young type Ia SNRs in the LMC and Galaxy, with a range of proposed degeneracy scenarios to highlight potential differences due to a different models. We suggest that the single degenerate scenario might point to morphologically asymmetric type Ia supernova explosions.

Primordial black hole (PBH) formation is a more generic phenomenon than was once thought. The dynamics of a scalar field in inflationary universe can produce PBHs under mild assumptions regarding the scalar potential. In the early universe, light scalar fields develop large expectation values during inflation and subsequently relax to the minimum of the effective potential at a later time. During the relaxation process, an initially homogeneous scalar condensate can fragment into lumps via an instability similar to the gravitational (Jeans) instability, where the scalar self-interactions, rather than gravity, play the leading role. The fragmentation of the scalar field into lumps (e.g. Q-balls or oscillons) creates matter composed of relatively few heavy "particles", whose distribution is subject to significant fluctuations unconstrained by comic microwave background (CMB) observations and unrelated to the large-scale structure. If this matter component comes to temporarily dominate the energy density before the scalar lumps decay, PBHs can be efficiently produced during the temporary matter-dominated era. We develop a general analytic framework for description of PBH formation in this class of models. We highlight the differences between the scalar fragmentation scenario and other commonly considered PBH formation models. Given the existence of the Higgs field and the preponderance of scalar fields within supersymmetric and other models of new physics, PBHs constitute an appealing and plausible candidate for dark matter.

We report the discovery of a very young high-mass X-ray binary (HMXB) system associated with the supernova remnant (SNR) MCSNRJ0513-6724 in the Large Magellanic Cloud (LMC), using XMM-Newton X-ray observations. The HMXB is located at the geometrical centre of extended soft X-ray emission, which we confirm as an SNR. The HMXB spectrum is consistent with an absorbed power law with spectral index ~1.6 and a luminosity of 7x10^{33} ergs/s (0.2--12 keV). Tentative X-ray pulsations are observed with a periodicity of 4.4 s and the OGLE I-band light curve of the optical counterpart from more than 17.5 years reveals a period of 2.2324\pm0.0003 d, which we interpret as the orbital period of the binary system. The X-ray spectrum of the SNR is consistent with non-equilibrium shock models as expected for young/less evolved SNRs. From the derived ionisation time scale we estimate the age of the SNR to be <6 kyr. The association of the HMXB with the SNR makes it the youngest HMXB, in the earliest evolutionary stage known to date. A HMXB as young as this can switch on as an accreting pulsar only when the spin period has reached a critical value. Under this assumption, we obtain an upper limit to the magnetic field of < 5x10^{11} G. This implies several interesting possibilities including magnetic field burial, possibly by an episode of post-supernova hyper-critical accretion. Since these fields are expected to diffuse out on a timescale of 10^{3}-10^{4} years, the discovery of a very young HMXB can provide us the unique opportunity to observe the evolution of the observable magnetic field for the first time in X-ray binaries.

We study the general infrared behavior of the power spectrum of a stochastic gravitational wave background produced by stress tensor in the form bilinear in certain dynamical degrees of freedom. We find $\Omega_{\text{GW}}\propto k^3$ for a very wide class of the sources which satisfy a set of reasonable conditions. Namely, the $k^3$ scaling is universally valid when the source term is bounded in both frequency and time, is effective in a radiation-dominated stage, and for $k$ smaller than all the physical scales associated with the source, like the peak frequency, peak width, and time duration, etc.. We also discuss possible violations of these conditions and their physical implications.

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

We present two new radio continuum images from the Australian Square Kilometre Array Pathfinder (ASKAP) survey in the direction of the Small Magellanic Cloud (SMC). These images are part of the Evolutionary Map of the Universe (EMU) Early Science Project (ESP) survey of the Small and Large Magellanic Clouds. The two new source lists produced from these images contain radio continuum sources observed at 960 MHz (4489 sources) and 1320 MHz (5954 sources) with a bandwidth of 192 MHz and beam sizes of 30.0"x30.0" and 16.3"x15.1", respectively. The median Root Mean Squared (RMS) noise values are 186$\mu$Jy beam$^{-1}$ (960 MHz) and 165$\mu$Jy beam$^{-1}$ (1320 MHz). To create point source catalogues, we use these two source lists, together with the previously published Molonglo Observatory Synthesis Telescope (MOST) and the Australia Telescope Compact Array (ATCA) point source catalogues to estimate spectral indices for the whole population of radio point sources found in the survey region. Combining our ASKAP catalogues with these radio continuum surveys, we found 7736 point-like sources in common over an area of 30 deg$^2$. In addition, we report the detection of two new, low surface brightness supernova remnant candidates in the SMC. The high sensitivity of the new ASKAP ESP survey also enabled us to detect the bright end of the SMC planetary nebula sample, with 22 out of 102 optically known planetary nebulae showing point-like radio continuum emission. Lastly, we present several morphologically interesting background radio galaxies.

Aims: We present a comprehensive study of the supernova remnant (SNR) population of the Small Magellanic Cloud (SMC). We measure multiwavelength properties of the SMC SNRs and compare them to those of the Large Magellanic Cloud (LMC) population. Methods: This study combines the large dataset of XMM-Newton observations of the SMC, archival and recent radio continuum observations, an optical line emission survey, and new optical spectroscopic observations. We can thus build a complete and clean sample of 19 confirmed and 4 candidate SNRs. The homogeneous X-ray spectral analysis allows to search for SN ejecta and Fe K line emission, and to measure interstellar medium (ISM) abundances. We estimate the ratio of core-collapse to type Ia supernova rates of the SMC based on the X-ray properties and the local stellar environment of each SNR. Results : After the removal of unconfirmed or misclassified objects, and the addition of two newly confirmed SNRs based on multi-wavelength features, we present a final list of 21 confirmed SNRs and 2 candidates. While no Fe K line is detected even for the brightest and youngest SNR, we find X-ray evidence of SN ejecta in 11 SNRs. We estimate a ratio of 4.7$_{-1.9} ^{+0.6}$ core-collapse supernova to every type Ia SN, three times higher than in the LMC. The difference can be ascribed to the absence of the enhanced star formation episode in the SMC, which occurred in the LMC 0.5-1.5 Gyr ago. The hot-gas abundances of O, Ne, Mg, and Fe are 0.1-0.2 times solar. Their ratios with respect to SMC stellar abundances reflect the effects of dust depletion and partial dust destruction in SNR shocks. We find evidence that the ambient medium probed by SMC SNRs is less disturbed and less dense on average than in the LMC, consistent with the different morphologies of the two galaxies.

We examine the new Galactic supernova remnant (SNR) candidate, G23.11+0.18, as seen by the Murchison Widefield Array (MWA) radio telescope. We describe the morphology of the candidate and find a spectral index of -0.63+/-0.05 in the 70-170MHz domain. A coincident TeV gamma-ray detection in High-Energy Stereoscopic System (HESS) data supports the SNR nature of G23.11+0.18 and suggests that G23.11+0.18 is accelerating particles beyond TeV energies, thus making this object a promising new cosmic ray hadron source candidate. The remnant cannot be seen in current optical, infrared and X-ray data-sets. We do find, however, a dip in CO-traced molecular gas at a line-of-sight velocity of ~85 km/s, suggesting the existence of a G23.11+0.18 progenitor wind-blown bubble. Furthermore, the discovery of molecular gas clumps at a neighbouring velocity towards HESS J1832-085 adheres to the notion that a hadronic gamma-ray production mechanism is plausible towards the north of the remnant. Based on these morphological arguments, we propose an interstellar medium association for G23.11+0.18 at a kinematic distance of 4.6+/-0.8 kpc.

RX J0046.5$-$7308 is a shell-type supernova remnant (SNR) in the Small Magellanic Cloud (SMC). We carried out new $^{12}$CO($J$ = 1-0, 3-2) observations toward the SNR using Mopra and the Atacama Submillimeter Telescope Experiment. We found eight molecular clouds (A-H) along the X-ray shell of the SNR. The typical cloud size and mass are $\sim$10-15 pc and $\sim$1000-3000 $M_{\odot}$, respectively. The X-ray shell is slightly deformed and has the brightest peak in the southwestern shell where two molecular clouds A and B are located. The four molecular clouds A, B, F, and G have high intensity ratios of $^{12}$CO($J$ = 3-2) / $^{12}$CO($J$ = 1-0) $> 1.2$, which are not attributable to any identified internal infrared sources or high-mass stars. The HI cavity and its expanding motion are found toward the SNR, which are likely created by strong stellar winds from a massive progenitor. We suggest that the molecular clouds A-D, F, and G and HI clouds within the wind-blown cavity at $V_\mathrm{LSR} = 117.1$-122.5 km s$^{-1}$ are to be associated with the SNR. The X-ray spectroscopy reveals the dynamical age of $26000^{+1000}_{-2000}$ yr and the progenitor mass of $\gtrsim 30$ $M_{\odot}$, which is also consistent with the proposed scenario. We determine physical conditions of the giant molecular cloud LIRS 36A using the large velocity gradient analysis with archival datasets of the Atacama Large Millimeter/submillimeter Array; the kinematic temperature is $72^{+50}_{-37}$ K and the number density of molecular hydrogen is $1500^{+600}_{-300}$ cm$^{-3}$. The next generation of $\gamma$-ray observations will allow us to study the pion-decay $\gamma$-rays from the molecular clouds in the SMC SNR.

We study the primordial tensor perturbation produced from the double inflationary scenario with an intermediate break stage. Because of the transitions, the power spectrum deviates from the vacuum one and there will appear oscillatory behavior. In the case of a scalar-type curvature perturbation, it is known that the amplitude of these oscillations may be enhanced to result in the power spectrum larger than the one for the vacuum case. One might expect the similar enhancement for the tensor perturbation as well. Unfortunately, it is found that when the equation of state (EOS) parameter $w=p/\rho$ of the break stage is a constant with $w>-1/3$, the amplitude of oscillations is never large enough to enhance the power spectrum. On the contrary, the power spectrum is found to be suppressed even on those scales that leave the horizon at the first inflationary stage and remain superhorizon throughout the entire stage. We identify the cause of this suppression with the correction terms in additional to the leading order constant solution on superhorizon scales. We argue that our result is general in the sense that any intermediate break stage during inflation cannot yield an enhancement of the tensor spectrum as long as the Hubble expansion rate is non-increasing in time.

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

A previous study of the X-ray luminosity function of the X-ray sources in the field of the Draco dwarf spheroidal (dSph) galaxy indicated the presence of a population of unknown X-ray sources in the soft energy range of 0.5-2 keV. In 2015, Draco dSph was observed again in twenty-six deep XMM-Newton, observations providing an opportunity for a new study of the yet unclassified sources. We apply the classification criteria presented in our previous multi-wavelength study of the X-ray sources of the Draco dSph to the sources detected in the combined 2009 and 2015 XMM-Newton data set. These criteria are based on X-ray studies and properties of the optical, near-infrared, and mid-infrared counterparts and allows us to distinguish background active galactic nuclei~(AGNs) and galaxies from other types of X-ray sources. We present the classification of X-ray sources, for which the counterpart is identified as a stellar object based on our criteria from multi-wavelength data. We identify three new symbiotic stars in the Draco dSph with X-ray luminosities between $\sim$3.5$\times10^{34}$ erg s$^{-1}$ and 5.5$\times10^{34}$ erg s$^{-1}$. The X-ray spectral analysis shows that two of the classified symbiotic stars are $\beta$-type. This is the first identification of this class of symbiotic stars in a nearby galaxy. Eight sources are classified as Galactic M dwarfs in the field of the Draco dSph. The distances of these M dwarfs are between$\sim$140-800 pc, their X-ray luminosities are between $10^{28}-10^{29}$ erg s$^{-1}$ and the logarithmic ratio of X-ray to bolometric luminosity, log$(\frac{L_\text{X}}{L_\text{bol}})$, is between $-3.4$ to $-2.1$. The multiple observations allowed us to investigate flare activity of the M dwarfs. Moreover, we classified three foreground sources, located at distances of the order of $\sim$1-3 kpc in the field of the Draco dSph.

We study gravitational waves (GWs) induced by non-Gaussian curvature perturbations. We calculate the density parameter per logarithmic frequency interval, $\Omega_\text{GW}(k)$, given that the power spectrum of the curvature perturbation $\mathcal{P}_\mathcal{R}(k)$ has a narrow peak at some small scale $k_*$, with a local-type non-Gaussianity, and constrain the nonlinear parameter $f_\text{NL}$ with the future LISA sensitivity curve as well as with constraints from the abundance of the primordial black holes (PBHs). We find that the non-Gaussian contribution to $\Omega_\text{GW}$ increases as $k^3$, peaks at $k/k_*=4/\sqrt{3}$, and has a sharp cutoff at $k=4k_*$. The non-Gaussian part can exceed the Gaussian part if $\mathcal{P}_\mathcal{R}(k)f_\text{NL}^2\gtrsim1$. If both a slope $\Omega_\text{GW}(k)\propto k^\beta$ with $\beta\sim3$ and the multiple-peak structure around a cutoff are observed, it can be recognized as a smoking gun of the primordial non-Gaussianity. We also find that if PBHs with masses of $10^{20}\text{g}$ to $10^{22}\text{g}$ are identified as cold dark matter of the Universe, the corresponding GWs must be detectable by LISA-like detectors, irrespective of the value of $\mathcal{P}_\mathcal{R}$ or $f_\text{NL}$.

In this paper, we present the analysis and results of a direct measurement of the cosmic-ray proton spectrum with the CALET instrument onboard the International Space Station, including the detailed assessment of systematic uncertainties. The observation period used in this analysis is from October 13, 2015 to August 31, 2018 (1054 days). We have achieved the very wide energy range necessary to carry out measurements of the spectrum from 50 GeV to 10 TeV covering, for the first time in space, with a single instrument the whole energy interval previously investigated in most cases in separate subranges by magnetic spectrometers (BESS-TeV, PAMELA, and AMS-02) and calorimetric instruments (ATIC, CREAM, and NUCLEON). The observed spectrum is consistent with AMS-02 but extends to nearly an order of magnitude higher energy, showing a very smooth transition of the power-law spectral index from -2.81 +- 0.03 (50--500 GeV) neglecting solar modulation effects (or -2.87 +- 0.06 including solar modulation effects in the lower energy region) to -2.56 +- 0.04 (1--10 TeV), thereby confirming the existence of spectral hardening and providing evidence of a deviation from a single power law by more than 3 sigma.

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

Supernova remnants serve as nearby laboratories relevant to many areas in Astrophysics, from stellar and galaxy evolution to extreme astrophysics and the formation of the heavy elements in the Universe. The Chandra X-ray mission has enabled a giant leap forward in studying both SNRs and their compact stellar remnants on sub-arcsecond scale. However, such high-resolution imaging studies have been mostly limited to the nearby and/or relatively bright objects. There is no question that we are missing a large population, especially in external galaxies. Within our own Galaxy, we are presented with new fundamental questions related to neutron stars' diversity, kicks, relativistic winds and the way these objects interact with, and impact, their host environments. In this white paper, we highlight some of the breakthroughs to be achieved with future X-ray missions (such as the proposed AXIS probe) equipped with sub-arcsecond imaging resolution and an order of magnitude improvement in sensitivity.

We report the discovery of a new Small Magellanic Cloud Pulsar Wind Nebula (PWN) at the edge of the Supernova Remnant (SNR)-DEM S5. The pulsar powered object has a cometary morphology similar to the Galactic PWN analogs PSR B1951+32 and 'the mouse'. It is travelling supersonically through the interstellar medium. We estimate the Pulsar kick velocity to be in the range of 700-2000 km/s for an age between 28-10 kyr. The radio spectral index for this SNR PWN pulsar system is flat (-0.29 $\pm$ 0.01) consistent with other similar objects. We infer that the putative pulsar has a radio spectral index of -1.8, which is typical for Galactic pulsars. We searched for dispersion measures (DMs) up to 1000 cm/pc^3 but found no convincing candidates with a S/N greater than 8. We produce a polarisation map for this PWN at 5500 MHz and find a mean fractional polarisation of P $\sim 23$ percent. The X-ray power-law spectrum (Gamma $\sim 2$) is indicative of non-thermal synchrotron emission as is expected from PWN-pulsar system. Finally, we detect DEM S5 in Infrared (IR) bands. Our IR photometric measurements strongly indicate the presence of shocked gas which is expected for SNRs. However, it is unusual to detect such IR emission in a SNR with a supersonic bow-shock PWN. We also find a low-velocity HI cloud of $\sim 107$ km/s which is possibly interacting with DEM S5. SNR DEM S5 is the first confirmed detection of a pulsar-powered bow shock nebula found outside the Galaxy.

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

The recent detection of 60Fe in the cosmic rays provides conclusive evidence that there is a recently synthesized component (few MY) in the GCRs (Binns et al. 2016). In addition, these nuclei must have been synthesized and accelerated in supernovae near the solar system, probably in the Sco-Cen OB association subgroups, which are about 100 pc distant from the Sun. Recent theoretical work on the production of r-process nuclei appears to indicate that it is difficult for SNe to produce the solar system abundances relative to iron of r-process elements with high atomic number (Z), including the actinides (Th, U, Np, Pu, and Cm). Instead, it is believed by many that the heaviest r-process nuclei, or perhaps even all r-process nuclei, are produced in binary neutron star mergers. Since we now know that there is at least a component of the GCRs that has been recently synthesized and accelerated, models of r-process production by SNe and BNSM can be tested by measuring the relative abundances of these ultra-heavy r-process nuclei, and especially the actinides, since they are radioactive and provide clocks that give the time interval from nucleosynthesis to detection at Earth. Since BNSM are believed to be much less frequent in our galaxy than SNe (roughly 1000 times less frequent, the ratios of the actinides, each with their own half-life, will enable a clear determination of whether the heaviest r-process nuclei are synthesized in SNe or in BNSM. In addition, the r-process nuclei for the charge range from 34 to 82 can be used to constrain models of r-process production in BNSM and SNe. Thus, GCRs become a multi-messenger component in the study of BNSM and SNe.

We show that there is a class of spacetime curvature singularities which can be resolved with metric and matter field transformations. As an example, we consider an anisotropic power-law inflation model with gauge and scalar fields in which a space-like curvature singularity exists at the beginning of time. First, we provide a transformation of the metric to the flat geometry. The transformation is regular in the whole region of spacetime except for the singularity. In general, matter fields are still singular after such a metric transformation. However, we explicitly show that there is a case in which the singular behavior of the matter fields can be completely removed by a field re-definition. Since the action is invariant under any metric and matter field transformations, the regularity of the action at the original singularity is a necessary condition for the complete removal of a singularity.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

The CALorimetric Electron Telescope (CALET) is a high-energy astroparticle physics space experiment installed on the International Space Station (ISS), developed and operated by Japan in collaboration with Italy and the United States. The CALET mission goals include the investigation of possible nearby sources of high-energy electrons, of the details of galactic particle acceleration and propagation, and of potential signatures of dark matter. CALET measures the cosmic-ray electron + positron flux up to 20 TeV, gamma-rays up to 10 TeV, and nuclei with Z=1 to 40 up to 1,000 TeV for the more abundant elements during a long-term observation aboard the ISS. Starting science operation in mid-October 2015, CALET performed continuous observation without major interruption with close to 20 million triggered events over 10 GeV per month. Based on the data taken during the first two-years, we present an overview of CALET observations: uses w/o major interruption 1) Electron + positron energy spectrum, 2) Nuclei analysis, 3) Gamma-ray observation including a characterization of on-orbit performance. Results of the electromagnetic counterpart search for LIGO/Virgo gravitational wave events are discussed as well.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

We carried out new $^{12}$CO($J$ = 1-0, 3-2) observations of a N63A supernova remnant (SNR) from the LMC using ALMA and ASTE. We find three giant molecular clouds toward the northeast, east, and near the center of the SNR. Using the ALMA data, we spatially resolved clumpy molecular clouds embedded within the optical nebulae in both the shock-ionized and photoionized lobes discovered by previous H$\alpha$ and [S II] observations. The total mass of the molecular clouds is $\sim$$800$ $M_{\odot}$ for the shock-ionized region and $\sim$$1700$ $M_{\odot}$ for the photoionized region. Spatially resolved X-ray spectroscopy reveals that the absorbing column densities toward the molecular clouds are $\sim$$1.5$-$6.0\times10^{21}$ cm$^{-2}$, which are $\sim$$1.5$-$15$ times less than the averaged interstellar proton column densities for each region. This means that the X-rays are produced not only behind the molecular clouds, but also in front of them. We conclude that the dense molecular clouds have been completely engulfed by the shock waves, but have still survived erosion owing to their high-density and short interacting time. The X-ray spectrum toward the gas clumps is well explained by an absorbed power-law or high-temperature plasma models in addition to the thermal plasma components, implying that the shock-cloud interaction is efficiently working for both the cases through the shock ionization and magnetic field amplification. If the hadronic gamma-ray is dominant in the GeV band, the total energy of cosmic-ray protons is calculated to be $\sim$$0.3$-$1.4\times10^{49}$ erg with the estimated ISM proton density of $\sim$$190\pm90$ cm$^{-3}$, containing both the shock-ionized gas and neutral atomic hydrogen.

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

KAGRA is a second-generation interferometric gravitational-wave detector with 3-km arms constructed at Kamioka, Gifu in Japan. It is now in its final installation phase, which we call bKAGRA (baseline KAGRA), with scientific observations expected to begin in late 2019. One of the advantages of KAGRA is its underground location of at least 200 m below the ground surface, which brings small seismic motion at low frequencies and high stability of the detector. Another advantage is that it cools down the sapphire test mass mirrors to cryogenic temperatures to reduce thermal noise. In April-May 2018, we have operated a 3-km Michelson interferometer with a cryogenic test mass for 10 days, which was the first time that km-scale interferometer was operated at cryogenic temperatures. In this article, we report the results of this "bKAGRA Phase 1" operation. We have demonstrated the feasibility of 3-km interferometer alignment and control with cryogenic mirrors.

We have compiled the most complete census of High-Mass X-ray Binaries (HMXBs) in the Small Magellanic Cloud with the aim to investigate the formation efficiency of young accreting binaries in its low metallicity environment. In total, we use 127 X-ray sources with detections in our \chandra X-ray Visionary Program (XVP), supplemented by 14 additional (likely and confirmed) HMXBs identified by \cite{2016A&A...586A..81H} that fall within the XVP area, but are not either detected in our survey (9 sources) or matched with any XVP source that has at least one OB counterpart in the OGLE-III catalog (5 sources). Specifically, we examine the number ratio of the HMXBs [N(HMXBs)] to {\it (a)} the number of OB stars, {\it (b)} the local star-formation rate (SFR), and {\it (c)} the stellar mass produced during the specific star-formation burst, all as a function of the age of their parent stellar populations. Each of these indicators serves a different role, but in all cases we find that the HMXB formation efficiency increases as a function of time (following a burst of star formation) up to $\sim$40--60\,Myr, and then gradually decreases. The peak formation efficiency N(HMXB)/SFR is (49 $\pm$ 14) $[10^{-5}~{\rm M_{\odot}/yr}]^{-1}$, in good agreement with previous estimates of the average formation efficiency in the broad $\sim$20--60\,Myr age range. The frequency of HMXBs is a factor of 8$\times$ higher than at $\sim$10\,Myr, and 4$\times$ higher than at $\sim$260\,Myr, i.e. at earlier and later epochs, respectively.

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

The recent detections of gravitational waves (GWs) reported by LIGO/Virgo collaborations have made significant impact on physics and astronomy. A global network of GW detectors will play a key role to solve the unknown nature of the sources in coordinated observations with astronomical telescopes and detectors. Here we introduce KAGRA (former name LCGT; Large-scale Cryogenic Gravitational wave Telescope), a new GW detector with two 3-km baseline arms arranged in the shape of an "L", located inside the Mt. Ikenoyama, Kamioka, Gifu, Japan. KAGRA's design is similar to those of the second generations such as Advanced LIGO/Virgo, but it will be operating at the cryogenic temperature with sapphire mirrors. This low temperature feature is advantageous for improving the sensitivity around 100 Hz and is considered as an important feature for the third generation GW detector concept (e.g. Einstein Telescope of Europe or Cosmic Explorer of USA). Hence, KAGRA is often called as a 2.5 generation GW detector based on laser interferometry. The installation and commissioning of KAGRA is underway and its cryogenic systems have been successfully tested in May, 2018. KAGRA's first observation run is scheduled in late 2019, aiming to join the third observation run (O3) of the advanced LIGO/Virgo network. In this work, we describe a brief history of KAGRA and highlights of main feature. We also discuss the prospects of GW observation with KAGRA in the era of O3. When operating along with the existing GW detectors, KAGRA will be helpful to locate a GW source more accurately and to determine the source parameters with higher precision, providing information for follow-up observations of a GW trigger candidate.

The 30 Dor C superbubble is unique for its synchrotron X-ray shell, as well as being the first superbubble to be detected in TeV $\gamma$-rays, though the dominant TeV emission mechanism, i.e., leptonic or hadronic, is still unclear. We aim to use new Chandra observations of 30 Dor C to resolve the synchrotron shell in unprecedented detail and to estimate the magnetic ($B$) field in the postshock region, a key discriminator between TeV $\gamma$-ray emission mechanisms. We extracted radial profiles in the 1.5-8 keV range from various sectors around the synchrotron shell and fitted these with a projected and point spread function convolved postshock volumetric emissivity model to determine the filament widths. We then calculated the postshock magnetic field strength from these widths. We found that most of the sectors were well fitted with our postshock model and the determined $B$-field values were low, all with best fits $\lesssim 20~\mu$G. Upper limits on the confidence intervals of three sectors reached $\gtrsim 30~\mu$G though these were poorly constrained. The generally low $B$-field values suggests a leptonic-dominated origin for the TeV $\gamma$-rays. Our postshock model did not provide adequate fits to two sectors. We found that one sector simply did not provide a clean enough radial profile, while the other could be fitted with a modified model where the projected profile falls off abruptly below $\sim0.8$ times the shell radius, yielding a postshock $B$-field of 4.8 (3.7-11.8) $\mu$G which is again consistent with the leptonic TeV $\gamma$-ray mechanism. Alternatively, the observed profiles in these sectors could result from enhancements around a shock-cloud interaction as suggested in previous works. The average $B$-field determined around the X-ray synchrotron shell of 30 Dor C suggests the leptonic scenario as the dominant emission mechanism for the TeV $\gamma$-rays.

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

We found that an intermediate velocity cloud (IVC) IVC 86-36 in HI 21 cm emission shows a head-tail distribution toward the Galactic plane with marked parallel filamentary streamers, which is extended over 40 degrees in the sky. The distance of IVC 86-36 is constrained to be less than ~3 kpc from absorption of a background star as determined from opticalspectroscopy. There is a bridge feature in velocity between the IVC and the local ISM with velocity separation of ~50 km s-1, which may indicate dynamical interaction of the IVC with the disk. If the interaction is correct, the distance estimate d of the IVC ranges from 200 pc to 3 kpc, and the mass of the IVC head is estimated to be 7X103(d/1kpc)2Msol. The IVC shares similar properties to the Smith cloud located at 12 kpc, including the head-tail distribution, streamers, and bridge feature, while the mass of the IVC is less than ~0.1 of the Smith cloud. A comparison between the Hi and the Planck/IRAS dust emission indicates that the dust emission of IVC 86-36 is not detectable in spite of its HI column density of 2X10^20 cm-2, indicating low metalicity of IVC 86-36 by a factor of ~< 0.2 as compared with the solar neighbor. We conclude that IVC 86-38 is an infalling cloud which likely originated in the low-metallicity environment of the Galactic halo or the Magellanic system.

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

We carried out new observations of two fields in the northern ring of M31 with XMM-Newton with two exposures of 100 ks each and obtained a complete list of X-ray sources down to a sensitivity limit of ~7 x 10^34 erg s^-1 (0.5 - 2.0 keV). The major objective of the observing programme was the study of the hot phase of the ISM in M31. The analysis of the diffuse emission and the study of the ISM is presented in a separate paper. We analysed the spectral properties of all detected sources using hardness ratios and spectra if the statistics were high enough. We also checked for variability. We cross-correlated the source list with the source catalogue of a new survey of the northern disc of M31 carried out with Chandra and Hubble (Panchromatic Hubble Andromeda Treasury, PHAT) as well as with other existing catalogues. We detected a total of 389 sources, including 43 foreground stars and candidates and 50 background sources. Based on the comparison to the Chandra/PHAT survey, we classify 24 hard X-ray sources as new candidates for X-ray binaries (XRBs). In total, we identified 34 XRBs and candidates and 18 supernova remnants (SNRs) and candidates. Three of the four brightest SNRs show emission mainly below 2 keV, consistent with shocked ISM. The spectra of two of them also require an additional component with a higher temperature. The SNR [SPH11] 1535 has a harder spectrum and might suggest that there is a pulsar-wind nebula inside the SNR. We find five new sources showing clear time variability. We also studied the spectral properties of the transient source SWIFT J004420.1+413702, which shows significant variation in flux over a period of seven months (June 2015 to January 2016) and associated change in absorption. Based on the likely optical counterpart detected in the Chandra/PHAT survey, the source is classified as a low-mass X-ray binary.

We study a model where quintessence potential $e^{-\xi\phi}$ coupled to Higgs potential. We calculate the evolution of the quintessence, and track the running of the effective Higgs self-coupling. We find it slightly larger than that of the standard model in the past. Requiring the electroweak vacuum to be absolutely stable in inflationary era, we find a lower bound $\xi> 0.35\pm 0.05$, where the uncertainty is mainly from the measurement of the top quark mass. This lower bound, together with the upper bound from the observation for dark energy $\xi\lesssim0.6$, narrows down the parameter space and makes it possible to test this model in the near future. Interestingly, the bound on $\xi$, if actually shown to be the case by observation, supports the recently proposed Swampland Conjecture.

The X-ray source populations within galaxies are typically difficult to identify and classify from X-ray data alone. We are able to break through this barrier by combining deep new Chandra ACIS-I observations with extensive Hubble Space Telescope imaging from the PHAT survey of the M31 disk. We detect 373 X-ray sources down to 0.35-8.0 keV flux of 10$^{-15}$ erg cm$^{-2}$ s$^{-1}$ over 0.4 square degrees, 170 of which are reported for the first time. We identify optical counterpart candidates for 188 of the 373 sources, after using the HST data to correct the absolute astrometry of our Chandra imaging to 0.1$"$. While 58 of these 188 are associated with point sources potentially in M31, over half (107) of the counterpart candidates are extended background galaxies, 5 are star clusters, 12 are foreground stars, and 6 are supernova remnants. Sources with no clear counterpart candidate are most likely to be undetected background galaxies and low-mass X-ray binaries in M31. The 58 point sources that are not consistent with foreground stars are bright enough that they could be high mass stars in M31; however, all but 8 have optical colors inconsistent with single stars, suggesting that many could be background galaxies or binary counterparts. For point-like counterparts, we examine the star formation history of the surrounding stellar populations to look for a young component that could be associated with a high mass X-ray binary. For the 40 point-like counterpart candidates associated with young populations, we find that their age distribution has two peaks at 15-20 Myr and 40-50 Myr. If we only consider the 8 counterpart candidates with typical high-mass main sequence optical star colors, their age distribution peaks mimic those of the sample of 40. Finally, we find that intrinsic faintness, and not extinction, is the main limitation for finding further counterpart candidates.

We point out that the gravitational-wave event GW150914 observed by the LIGO detectors can be explained by the coalescence of primordial black holes (PBHs). It is found that the expected PBH merger rate would exceed the rate estimated by the LIGO scientific Collaboration and Virgo Collaboration if PBHs were the dominant component of dark matter, while it can be made compatible if PBHs constitute a fraction of dark matter. Intriguingly, the abundance of PBHs required to explain the suggested lower bound on the event rate, $> 2$ events ${\rm Gpc}^{-3} {\rm yr}^{-1}$, roughly coincides with the existing upper limit set by the nondetection of the cosmic microwave background spectral distortion. This implies that the proposed PBH scenario may be tested in the not-too-distant future.

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

The Advanced Telescope for High Energy Astrophysics (Athena) is the X-ray observatory large mission selected by the European Space Agency (ESA), within its Cosmic Vision 2015-2025 programme, to address the "Hot and Energetic Universe" scientific theme (Nandra et al. 2013), and it is provisionally due for launch in the early 2030s. The Square Kilometer Array (SKA) is the next generation radio observatory and consists of two telescopes, one comprised of dishes operating at mid frequencies (SKA1-MID) and located in South Africa, and the other comprised of Log-Periodic antennas operating at low radio frequencies (SKA1-LOW), which will be located in Australia (Braun et al. 2017). The scientific commissioning of the radio telescope is planned to begin in 2021-2022. The SKA-Athena Synergy Team (SAST) has been tasked to single out the potential scientific synergies between Athena and SKA. The astrophysical community was involved in this exercise primarily through a dedicated SKA-Athena Synergy Workshop, which took place on April 24-25, 2017 at SKAO, Jodrell Bank, Manchester. The final result of the synergy exercise, this White Paper, describes in detail a number of scientific opportunities that will be opened up by the combination of Athena and SKA, these include: 1. the Cosmic Dawn; 2. the Evolution of black holes and galaxies; 3. Active galaxy feedback in galaxy clusters; 4. Non-thermal phenomena in galaxy clusters; 5. Detecting the cosmic web; 6. Black-hole accretion physics and astrophysical transients; 7. Galactic astronomy: stars, planets, pulsars and supernovae.

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.

We present results on searches for gamma-ray counterparts of the LIGO/Virgo gravitational-wave events using CALorimetric Electron Telescope ({\sl CALET}) observations. The main instrument of {\sl CALET}, CALorimeter (CAL), observes gamma-rays from $\sim1$ GeV up to 10 TeV with a field of view of nearly 2 sr. In addition, the {\sl CALET} gamma-ray burst monitor (CGBM) views $\sim$3 sr and $\sim2\pi$ sr of the sky in the 7 keV -- 1 MeV and the 40 keV -- 20 MeV bands, respectively, by using two different crystal scintillators. The {\sl CALET} observations on the International Space Station started in October 2015, and here we report analyses of events associated with the following gravitational wave events: GW151226, GW170104, GW170608, GW170814 and GW170817. Although only upper limits on gamma-ray emission are obtained, they correspond to a luminosity of $10^{49}\sim10^{53}$ erg s$^{-1}$ in the GeV energy band depending on the distance and the assumed time duration of each event, which is approximately the order of luminosity of typical short gamma-ray bursts. This implies there will be a favorable opportunity to detect high-energy gamma-ray emission in further observations if additional gravitational wave events with favorable geometry will occur within our field-of-view. We also show the sensitivity of {\sl CALET} for gamma-ray transient events which is the order of $10^{-7}$~erg\,cm$^{-2}$\,s$^{-1}$ for an observation of 100~s duration.

G15.9+0.2 is a Galactic shell-type supernova remnant (SNR), which was detected in radio and has been confirmed in X-rays based on Chandra observations. An X-ray point source CXOUJ181852.0-150213 has been detected and suggested to be an associated neutron star. In a recent study, we have confirmed the source to be a central compact object (CCO). We have studied the SNR using high-resolution X-ray data taken with Chandra in combination with infrared (IR) data in order to understand its emission and to derive its physical parameters. This will also help to constrain, e.g., the age of the CCO and the environment in which it was born. The spectral analysis of the X-ray emission using the new Chandra data and the comparison to the IR data have shown that the SNR is relatively young with an age of a few thousand years and that its emission is dominated by that of shocked interstellar medium (ISM). However, the analysis of the spectrum of the bright eastern shell shows that there is an additional emission component with enhanced abundances of {\alpha} elements and Fe, suggesting ejecta emission. The multi-wavelength emission is consistent with SNR G15.9+0.2 expanding in an ISM with a density gradient, while there is also colder material located in front of the SNR, which absorbs its thermal X-ray emission in the softer bands.

Extended results on the cosmic-ray electron + positron spectrum from 11 GeV to 4.8 TeV are presented based on observations with the Calorimetric Electron Telescope (CALET) on the International Space Station utilizing the data up to November 2017. The analysis uses the full detector acceptance at high energies, approximately doubling the statistics compared to the previous result. CALET is an all-calorimetric instrument with a total thickness of 30 $X_0$ at normal incidence and fine imaging capability, designed to achieve large proton rejection and excellent energy resolution well into the TeV energy region. The observed energy spectrum in the region below 1 TeV shows good agreement with Alpha Magnetic Spectrometer (AMS-02) data. In the energy region below $\sim$300 GeV, CALET's spectral index is found to be consistent with the AMS-02, Fermi Large Area Telescope (Fermi-LAT) and Dark Matter Particle Explorer (DAMPE), while from 300 to 600 GeV the spectrum is significantly softer than the spectra from the latter two experiments. The absolute flux of CALET is consistent with other experiments at around a few tens of GeV. However, it is lower than those of DAMPE and Fermi-LAT with the difference increasing up to several hundred GeV. The observed energy spectrum above $\sim$1 TeV suggests a flux suppression consistent within the errors with the results of DAMPE, while CALET does not observe any significant evidence for a narrow spectral feature in the energy region around 1.4 TeV. Our measured all-electron flux, including statistical errors and a detailed breakdown of the systematic errors, is tabulated in the Supplemental Material in order to allow more refined spectral analyses based on our data.

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

We study the degeneracy of the primordial spectrum and bispectrum of the curvature perturbation in single field inflationary models with a class of features in the inflaton potential. The feature we consider is a discontinuous change in the shape of the potential and is controlled by a couple of parameters that describe the strength of the discontinuity and the change in the potential shape. This feature produces oscillations of the spectrum and bispectrum around the comoving scale $k=k_0$ that exits the horizon when the inflaton passes the discontinuity. We find that the effects on the spectrum and almost all configurations of the bispectrum including the squeezed limit depend on a single quantity which is a function of the two parameters defining the feature. This leads to a degeneracy, i.e. different features of the inflaton potential can produce the same observational effects. However, we find that the degeneracy in the bispectrum is removed at the equilateral limit around $k=k_0$. This can be used to discriminate different models which give the same spectrum.

Given observations of B-mode polarization power spectrum of the cosmic microwave background (CMB), we can reconstruct power spectra of primordial tensor modes from the early Universe without assuming their functional form such as a power-law spectrum. Shape of the reconstructed spectra can then be used to probe the origin of tensor modes in a model-independent manner. We use the Fisher matrix to calculate the covariance matrix of tensor power spectra reconstructed in bins. We find that the power spectra are best reconstructed at wavenumbers in the vicinity of $k\approx 6\times 10^{-4}$ and $5\times 10^{-3}~{\rm Mpc}^{-1}$, which correspond to the "reionization bump" at $\ell\lesssim 6$ and "recombination bump" at $\ell\approx 80$ of the CMB B-mode power spectrum, respectively. The error bar between these two wavenumbers is larger because of lack of the signal between the reionization and recombination bumps. The error bars increase sharply towards smaller (larger) wavenumbers because of the cosmic variance (CMB lensing and instrumental noise). To demonstrate utility of the reconstructed power spectra we investigate whether we can distinguish between various sources of tensor modes including those from the vacuum metric fluctuation and SU(2) gauge fields during single-field slow-roll inflation, open inflation and massive gravity inflation. The results depend on the model parameters, but we find that future CMB experiments are sensitive to differences in these models. We make our calculation tool available on-line.

According to the equivalence principal, the long wavelength perturbations must not have any dynamical effect on the short scale physics up to ${\cal O} (k_L^2/k_s^2)$. Their effect can be always absorbed to a coordinate transformation locally. So any physical effect of such a perturbation appears only on scales larger than the scale of the perturbation. The bispectrum in the squeezed limit of the curvature perturbation in single-field slow-roll inflation is a good example, where the long wavelength effect is encoded in the spectral index through Maldacena's consistency relation. This implies that one should be able to derive the bispectrum in the squeezed limit without resorting to the in-in formalism in which one computes perturbative corrections field-theoretically. In this short paper, we show that the $\delta N$ formalism as it is, or more generically the separate universe approach, when applied carefully can indeed lead to the correct result for the bispectrum in the squeezed limit. Hence despite the common belief that the $\delta N$ formalism is incapable of recovering the consistency relation within itself, it is in fact self-contained and consistent.

We study a model of inflation in which a scalar field $\chi$ is non-minimally coupled to Starobinsky's $R^2$ gravity. After transforming it to the Einstein frame, a new scalar field, the scalaron $\phi$, will appear and couple to $\chi$ with a nontrivial field metric, while $\chi$ acquires a positive mass via the non-minimal coupling. Initially inflation occurs along the $\phi$ direction with $\chi$ trapped near its origin by this induced mass. After $\phi$ crosses a critical value, it starts rolling down rapidly and proceeds todamped oscillations around an effective local minimum determined by the value of $\chi$, while inflation still continues, driven by the $\chi$ field at this second stage where the effect of the non-minimal coupling becomes negligible. The presence of the damped oscillations during the transition from the first to second stage of inflation causes enhancement and oscillation features in the power spectrum of the curvature perturbation. Assuming that the oscillations may be treated perturbatively, we calculate these features by using the $\delta N$ formalism, and discuss its observational implications to large scale CMB anomalies or primordial black hole formation, depending on the scale of the features.

Spectral lines are among the most powerful signatures for dark matter (DM) annihilation searches in very-high-energy $\gamma$-rays. The central region of the Milky Way halo is one of the most promising targets given its large amount of DM and proximity to Earth. We report on a search for a monoenergetic spectral line from self-annihilations of DM particles in the energy range from 300 GeV to 70 TeV using a two-dimensional maximum likelihood method taking advantage of both the spectral and spatial features of signal versus background. The analysis makes use of Galactic Center (GC) observations accumulated over ten years (2004 - 2014) with the H.E.S.S. array of ground-based Cherenkov telescopes. No significant $\gamma$-ray excess above the background is found. We derive upper limits on the annihilation cross section $\langle\sigma v\rangle$ for monoenergetic DM lines at the level of $\sim4\times10^{-28}$ cm$^{3}$s$^{-1}$ at 1 TeV, assuming an Einasto DM profile for the Milky Way halo. For a DM mass of 1 TeV, they improve over the previous ones by a factor of six. The present constraints are the strongest obtained so far for DM particles in the mass range 300 GeV - 70 TeV. Ground-based $\gamma$-ray observations have reached sufficient sensitivity to explore relevant velocity-averaged cross sections for DM annihilation into two $\gamma$-ray photons at the level expected from the thermal relic density for TeV DM particles.

A search for new supernova remnants (SNRs) has been conducted using TeV gamma-ray data from the H.E.S.S. Galactic plane survey. As an identification criterion, shell morphologies that are characteristic for known resolved TeV SNRs have been used. Three new SNR candidates were identified in the H.E.S.S. data set with this method. Extensive multiwavelength searches for counterparts were conducted. A radio SNR candidate has been identified to be a counterpart to HESS J1534-571. The TeV source is therefore classified as a SNR. For the other two sources, HESS J1614-518 and HESS J1912+101, no identifying counterparts have been found, thus they remain SNR candidates for the time being. TeV-emitting SNRs are key objects in the context of identifying the accelerators of Galactic cosmic rays. The TeV emission of the relativistic particles in the new sources is examined in view of possible leptonic and hadronic emission scenarios, taking the current multiwavelength knowledge into account.

Aims. We study gamma-ray emission from the shell-type supernova remnant (SNR) RX J0852.0$-$4622 to better characterize its spectral properties and its distribution over the SNR. Methods. The analysis of an extended High Energy Spectroscopic System (H.E.S.S.) data set at very high energies (E > 100 GeV) permits detailed studies, as well as spatially resolved spectroscopy, of the morphology and spectrum of the whole RX J0852.0$-$4622 region. The H.E.S.S. data are combined with archival data from other wavebands and interpreted in the framework of leptonic and hadronic models. The joint Fermi-LAT-H.E.S.S. spectrum allows the direct determination of the spectral characteristics of the parent particle population in leptonic and hadronic scenarios using only GeV-TeV data. Results. An updated analysis of the H.E.S.S. data shows that the spectrum of the entire SNR connects smoothly to the high-energy spectrum measured by Fermi-LAT. The increased data set makes it possible to demonstrate that the H.E.S.S. spectrum deviates significantly from a power law and is well described by both a curved power law and a power law with an exponential cutoff at an energy of $E_\mathrm{cut} = (6.7 \pm 1.2_\mathrm{stat} \pm 1.2_\mathrm{syst})$ TeV. The joint Fermi-LAT-H.E.S.S. spectrum allows the unambiguous identification of the spectral shape as a power law with an exponential cutoff. No significant evidence is found for a variation of the spectral parameters across the SNR, suggesting similar conditions of particle acceleration across the remnant. A simple modeling using one particle population to model the SNR emission demonstrates that both leptonic and hadronic emission scenarios remain plausible. It is also shown that at least a part of the shell emission is likely due to the presence of a pulsar wind nebula around PSR J0855$-$4644.

The CALorimetric Electron Telescope (CALET), launched for installation on the International Space Station (ISS) in August, 2015, has been accumulating scientific data since October, 2015. CALET is intended to perform long-duration observations of high-energy cosmic rays onboard the ISS. CALET directly measures the cosmic-ray electron spectrum in the energy range of 1 GeV to 20 TeV with a 2% energy resolution above 30 GeV. In addition, the instrument can measure the spectrum of gamma rays well into the TeV range, and the spectra of protons and nuclei up to a PeV. In order to operate the CALET onboard ISS, JAXA Ground Support Equipment (JAXA-GSE) and the Waseda CALET Operations Center (WCOC) have been established. Scientific operations using CALET are planned at WCOC, taking into account orbital variations of geomagnetic rigidity cutoff. Scheduled command sequences are used to control the CALET observation modes on orbit. Calibration data acquisition by, for example, recording pedestal and penetrating particle events, a low-energy electron trigger mode operating at high geomagnetic latitude, a low-energy gamma-ray trigger mode operating at low geomagnetic latitude, and an ultra heavy trigger mode, are scheduled around the ISS orbit while maintaining maximum exposure to high-energy electrons and other high-energy shower events by always having the high-energy trigger mode active. The WCOC also prepares and distributes CALET flight data to collaborators in Italy and the United States. As of August 31, 2017, the total observation time is 689 days with a live time fraction of the total time of approximately 84%. Nearly 450 million events are collected with a high-energy (E>10 GeV) trigger. By combining all operation modes with the excellent-quality on-orbit data collected thus far, it is expected that a five-year observation period will provide a wealth of new and interesting results.

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.