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 Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

The $\rho$ Oph molecular cloud is one of the best examples of spinning dust emission, first detected by the Cosmic Background Imager (CBI). Here we present 4.5 arcmin observations with CBI 2 that confirm 31 GHz emission from $\rho$ Oph W, the PDR exposed to B-type star HD 147889, and highlight the absence of signal from S1, the brightest IR nebula in the complex. In order to quantify an association with dust-related emission mechanisms, we calculated correlations at different angular resolutions between the 31 GHz map and proxies for the column density of IR emitters, dust radiance and optical depth templates. We found that the 31 GHz emission correlates best with the PAH column density tracers, while the correlation with the dust radiance improves when considering emission that is more extended (from the shorter baselines), suggesting that the angular resolution of the observations affects the correlation results. A proxy for the spinning dust emissivity reveals large variations within the complex, with a dynamic range of 25 at 3$\sigma$ and a variation by a factor of at least 23, at 3$\sigma$, between the peak in $\rho$ Oph W and the location of S1, which means that environmental factors are responsible for boosting spinning dust emissivities locally.

Forthcoming space-based observations will require high-quality point-spread function (PSF) models for weak gravitational lensing measurements. One approach to generating these models is using a wavefront model based on the known telescope optics. We present an empirical framework for validating such models to confirm that they match the actual PSF to within requirements by comparing the models to the observed light distributions of isolated stars. We apply this framework to Tiny Tim, the standard tool for generating model PSFs for the Hubble Space Telescope (HST), testing its models against images taken by HST's Advanced Camera for Surveys in the Wide Field Channel. We show that Tiny Tim's models, in the default configuration, differ significantly from the observed PSFs, most notably in their sizes. We find that the quality of Tiny Tim PSFs can be improved through fitting the full set of Zernike polynomial coefficients which characterise the optics, to the point where the practical significance of the difference between model and observed PSFs is negligible for most use cases, resulting in additive and multiplicative biases both of order approximately 4e-4. We also show that most of this improvement can be retained through using an updated set of Zernike coefficients, which we provide.

We report the discovery of a remnant radio AGN J1615+5452 in the field of ELAIS-N1. GMRT continuum observations at 150, 325 and 610 MHz combined with archival data from the 1.4 GHz NVSS survey were used to derive the radio spectrum of the source. At a redshift $z \sim$ 0.33, J1615+5452 has a linear size of $\sim$ 100 kpc and spectral indices ranging between $\alpha^{1400}_{610} < -1.5$ and $\alpha^{325}_{150} = -0.61 \pm 0.12$. While the source has a diffuse radio emission at low frequencies, we do not find evidence of core, jets or hotspots in the 1.4\,GHz VLA data of $\sim 5$ arcsec angular resolution. Such morphological properties coupled with a curved radio spectrum suggest that the AGN fueling mechanisms undergo a shortage of energy supply which is typical of a dying radio AGN. This is consistent with the observed steep curvature in the spectrum $\Delta\alpha \approx -1$, the estimated synchrotron age of $t_{\rm s}=76.0^{+7.4}_{-8.7}$ Myr and a $t_{\rm off}/t_{\rm s}$ ratio of $\sim 0.3$.

We present a point-source detection algorithm that employs the second order Spherical Mexican Hat wavelet filter (SMHW2), and use it on C-BASS northern intensity data to produce a catalogue of point-sources. This catalogue allows us to cross-check the C-BASS flux-density scale against existing source surveys, and provides the basis for a source mask which will be used in subsequent C-BASS and cosmic microwave background (CMB) analyses. The SMHW2 allows us to filter the entire sky at once, avoiding complications from edge effects arising when filtering small sky patches. The algorithm is validated against a set of Monte Carlo simulations, consisting of diffuse emission, instrumental noise, and various point-source populations. The simulated source populations are successfully recovered. The SMHW2 detection algorithm is used to produce a $4.76\,\mathrm{GHz}$ northern sky source catalogue in total intensity, containing 1784 sources and covering declinations $\delta\geq-10^{\circ}$. The C-BASS catalogue is matched with the Green Bank 6\,cm (GB6) and Parkes-MIT-NRAO (PMN) catalogues over their areas of common sky coverage. From this we estimate the $90$ per cent completeness level to be approximately $610\,\mathrm{mJy}$, with a corresponding reliability of $98$ per cent, when masking the brightest $30$ per cent of the diffuse emission in the C-BASS northern sky map. We find the C-BASS and GB6 flux-density scales to be consistent with one another to within approximately $4$ per cent.

The capacity of Stage IV lensing surveys to measure the neutrino mass sum and differentiate between the normal and inverted mass hierarchies depends on the impact of nuisance parameters describing small-scale baryonic astrophysics and intrinsic alignments. For a Euclid-like survey, we perform the first combined weak lensing and galaxy clustering Fisher analysis with baryons, intrinsic alignments, and massive neutrinos for both hierarchies. We use a matter power spectrum generated from a halo model that captures the impact of baryonic feedback and adiabatic contraction. For weak lensing, we find that baryons cause severe degradation to forecasts of the neutrino mass sum, $\Sigma$, approximately doubling $\sigma_{\Sigma}$. We show that including galaxy clustering constraints from Euclid and BOSS, and cosmic microwave background (CMB) Planck priors, can reduce this degradation to $\sigma_{\Sigma}$ to 9% and 16% for the normal and inverted hierarchies respectively. The combined forecasts, $\sigma_{\Sigma_{\rm{NH}}}=0.034\, \rm{eV}$ and $\sigma_{\Sigma_{\rm{IH}}}=0.034\, \rm{eV}$, preclude a meaningful distinction of the hierarchies but could be improved upon with future CMB priors on $n_s$ and information from neutrinoless double beta decay to achieve a 2$\sigma$ distinction. The effect of intrinsic alignments on forecasts is shown to be minimal, with $\sigma_{\Sigma}$ even experiencing mild improvements due to information from the intrinsic alignment signal. We find that while adiabatic contraction and intrinsic alignments will require careful calibration to prevent significant biasing of $\Sigma$, there is less risk presented by feedback from energetic events like AGN and supernovae.

The Square Kilometre Array (SKA) will answer fundamental questions about the origin, evolution, properties, and influence of magnetic fields throughout the Universe. Magnetic fields can illuminate and influence phenomena as diverse as star formation, galactic dynamics, fast radio bursts, active galactic nuclei, large-scale structure, and Dark Matter annihilation. Preparations for the SKA are swiftly continuing worldwide, and the community is making tremendous observational progress in the field of cosmic magnetism using data from a powerful international suite of SKA pathfinder and precursor telescopes. In this contribution, we revisit community plans for magnetism research using the SKA, in the light of these recent rapid developments. We focus in particular on the impact that new radio telescope instrumentation is generating, thus advancing our understanding of key SKA magnetism science areas, as well as the new techniques that are required for processing and interpreting the data. We discuss these recent developments in the context of the ultimate scientific goals for the SKA era.

We update the constraints on the location of the nearest UHECR source. By analyzing recent data from the Pierre Auger Observatory using start-of-the-art CR propagation models, we reaffirm the need of local sources with a distance less than 25-100~Mpc, depending on mass composition. A new fast semi-analytical method for the propagation of UHECR in environments with turbulent magnetic fields is developed. The onset of an enhancement and a low-energy magnetic horizon of cosmic rays from sources located within a particular distance range is demonstrated. We investigate the distance to the nearest source, taking into account these magnetic field effects. The results obtained highlight the robustness of our constrained distances to the nearest source.

An analysis is presented of the optical polarimetric and multicolour photometric ($BVRJ$) behaviour of the blazar PKS 2155$-$304 during an outburst in 2010. This flare develops over roughly 117 days, with a flux doubling time $\tau \sim 11$ days that increases from blue to red wavelengths. The polarization angle is initially aligned with the jet axis but rotates by roughly $90^\circ$ as the flare grows. Two distinct states are evident at low and high fluxes. Below 18 mJy, the polarization angle takes on a wide range of values, without any clear relation to the flux. In contrast, there is a positive correlation between the polarization angle and flux above 18 mJy. The polarization degree does not display a clear correlation with the flux. We find that the photopolarimetric behaviour for the high flux state can be attributed to a variable component with a steady power-law spectral energy distribution and high optical polarization degree (13.3%). These properties are interpreted within the shock-in-jet model, which shows that the observed variability can be explained by a shock that is seen nearly edge-on. Some parameters derived for the relativistic jet within the shock-in-jet model are: $B=0.06$ G for the magnetic field, $\delta=22.3$ for the Doppler factor and $\Phi=2.6^\circ$ for the viewing angle.

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

The Large Scale Polarization Explorer (LSPE) is a cosmology program for the measurement of large scale curl-like features (B-modes) in the polarization of the Cosmic Microwave Background. Its goal is to constrain the background of inflationary gravity waves traveling through the universe at the time of matter-radiation decoupling. The two instruments of LSPE are meant to synergically operate by covering a large portion of the northern microwave sky. LSPE/STRIP is a coherent array of receivers planned to be operated from the Teide Observatory in Tenerife, for the control and characterization of the low-frequency polarized signals of galactic origin; LSPE/SWIPE is a balloon-borne bolometric polarimeter based on 330 large throughput multi-moded detectors, designed to measure the CMB polarization at 150 GHz and to monitor the polarized emission by galactic dust above 200 GHz. The combined performance and the expected level of systematics mitigation will allow LSPE to constrain primordial B-modes down to a tensor/scalar ratio of $10^{-2}$. We here report the status of the STRIP pre-commissioning phase and the progress in the characterization of the key subsystems of the SWIPE payload (namely the cryogenic polarization modulation unit and the multi-moded TES pixels) prior to receiver integration.

The data from the Euclid mission will enable the measurement of the photometric redshifts, angular positions, and weak lensing shapes for over a billion galaxies. This large dataset will allow for cosmological analyses using the angular clustering of galaxies and cosmic shear. The cross-correlation (XC) between these probes can tighten constraints and it is therefore important to quantify their impact for Euclid. In this study we carefully quantify the impact of XC not only on the final parameter constraints for different cosmological models, but also on the nuisance parameters. In particular, we aim at understanding the amount of additional information that XC can provide for parameters encoding systematic effects, such as galaxy bias or intrinsic alignments (IA). We follow the formalism presented in Euclid Collaboration: Blanchard et al. (2019) and make use of the codes validated therein. We show that XC improves the dark energy Figure of Merit (FoM) by a factor $\sim 5$, whilst it also reduces the uncertainties on galaxy bias by $\sim 17\%$ and the uncertainties on IA by a factor $\sim 4$. We observe that the role of XC on the final parameter constraints is qualitatively the same irrespective of the galaxy bias model used. We also show that XC can help in distinguishing between different IA models, and that if IA terms are neglected then this can lead to significant biases on the cosmological parameters. We find that the XC terms are necessary to extract the full information content from the data in future analyses. They help in better constraining the cosmological model, and lead to a better understanding of the systematic effects that contaminate these probes. Furthermore, we find that XC helps in constraining the mean of the photometric-redshift distributions, but it requires a more precise knowledge of this mean in order not to degrade the final FoM. [Abridged]

Context. Future weak lensing surveys, such as the Euclid mission, will attempt to measure the shapes of billions of galaxies in order to derive cosmological information. These surveys will attain very low levels of statistical error, and systematic errors must be extremely well controlled. In particular, the point spread function (PSF) must be estimated using stars in the field, and recovered with high accuracy. Aims. The aims of this paper are twofold. Firstly, we took steps toward a nonparametric method to address the issue of recovering the PSF field, namely that of finding the correct PSF at the position of any galaxy in the field, applicable to Euclid. Our approach relies solely on the data, as opposed to parametric methods that make use of our knowledge of the instrument. Secondly, we studied the impact of imperfect PSF models on the shape measurement of galaxies themselves, and whether common assumptions about this impact hold true in an Euclid scenario. Methods. We extended the recently proposed resolved components analysis approach, which performs super-resolution on a field of under-sampled observations of a spatially varying, image-valued function. We added a spatial interpolation component to the method, making it a true 2-dimensional PSF model. We compared our approach to PSFEx, then quantified the impact of PSF recovery errors on galaxy shape measurements through image simulations. Results. Our approach yields an improvement over PSFEx in terms of the PSF model and on observed galaxy shape errors, though it is at present far from reaching the required Euclid accuracy. We also find that the usual formalism used for the propagation of PSF model errors to weak lensing quantities no longer holds in the case of an Euclid-like PSF. In particular, different shape measurement approaches can react differently to the same PSF modeling errors.

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

We present the luminosity function (LF) for ultraluminous Ly$\alpha$ emitting galaxies (LAEs) at z = 6.6. We define ultraluminous LAEs (ULLAEs) as galaxies with logL(Ly$\alpha$) > 43.5 erg s$^{-1}$. We select our main sample using the g', r', i', z', and NB921 observations of a wide-area (30 deg$^2$) Hyper Suprime-Cam survey of the North Ecliptic Pole (NEP) field. We select candidates with g', r', i' > 26, NB921 $\leq$ 23.5, and NB921 - z' $\leq$ 1.3. Using the DEIMOS spectrograph on Keck II, we confirm 9 of our 14 candidates as ULLAEs at z = 6.6 and the remaining 5 as an AGN at z = 6.6, two [OIII]$\lambda$5007 emitting galaxies at z = 0.84 and z = 0.85, and two non-detections. This emphasizes the need for full spectroscopic follow-up to determine accurate LFs. In constructing the ULLAE LF at z = 6.6, we combine our 9 NEP ULLAEs with two previously discovered and confirmed ULLAEs in the COSMOS field: CR7 and COLA1. We apply rigorous corrections for incompleteness based on simulations. We compare our ULLAE LF at z = 6.6 with LFs at z = 5.7 and z = 6.6 from the literature. Our data reject some previous LF normalizations and power law indices, but they are broadly consistent with others. Indeed, a comparative analysis of the different literature LFs suggests that none is fully consistent with any of the others, making it critical to determine the evolution from z = 5.7 to z = 6.6 using LFs constructed in exactly the same way at both redshifts.

The Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a flattened shape, discrete geological units, and noticeable albedo heterogeneity. However, there is little surface color and compositional heterogeneity. No evidence for satellites, ring or dust structures, gas coma, or solar wind interactions was detected. By origin MU69 appears consistent with pebble cloud collapse followed by a low velocity merger of its two lobes.

Stage IV weak lensing experiments will offer more than an order of magnitude leap in precision. We must therefore ensure that our analyses remain accurate in this new era. Accordingly, previously ignored systematic effects must be addressed. In this work, we evaluate the impact of the reduced shear approximation and magnification bias, on the information obtained from the angular power spectrum. To first-order, the statistics of reduced shear, a combination of shear and convergence, are taken to be equal to those of shear. However, this approximation can induce a bias in the cosmological parameters that can no longer be neglected. A separate bias arises from the statistics of shear being altered by the preferential selection of galaxies and the dilution of their surface densities, in high-magnification regions. The corrections for these systematic effects take similar forms, allowing them to be treated together. We calculated the impact of neglecting these effects on the cosmological parameters that would be determined from Euclid, using cosmic shear tomography. To do so, we employed the Fisher matrix formalism, and included the impact of the super-sample covariance. We also demonstrate how the reduced shear correction can be calculated using a lognormal field forward modelling approach. These effects cause significant biases in Omega_m, sigma_8, n_s, Omega_DE, w_0, and w_a of -0.53 sigma, 0.43 sigma, -0.34 sigma, 1.36 sigma, -0.68 sigma, and 1.21 sigma, respectively. We then show that these lensing biases interact with another systematic: the intrinsic alignment of galaxies. Accordingly, we develop the formalism for an intrinsic alignment-enhanced lensing bias correction. Applying this to Euclid, we find that the additional terms introduced by this correction are sub-dominant.

The Euclid mission will observe well over a billion galaxies out to $z\sim6$ and beyond. This will offer an unrivalled opportunity to investigate several key questions for understanding galaxy formation and evolution. The first step for many of these studies will be the selection of a sample of quiescent and star-forming galaxies, as is often done in the literature by using well known colour techniques such as the `UVJ' diagram. However, given the limited number of filters available for the Euclid telescope, the recovery of such rest-frame colours will be challenging. We therefore investigate the use of observed Euclid colours, on their own and together with ground-based u-band observations, for selecting quiescent and star-forming galaxies. The most efficient colour combination, among the ones tested in this work, consists of the (u-VIS) and (VIS-J) colours. We find that this combination allows users to select a sample of quiescent galaxies complete to above $\sim70\%$ and with less than 15$\%$ contamination at redshifts in the range $0.75<z<1$. For galaxies at high-z or without the u-band complementary observations, the (VIS-Y) and (J-H) colours represent a valid alternative, with $>65\%$ completeness level and contamination below 20$\%$ at $1<z<2$ for finding quiescent galaxies. In comparison, the sample of quiescent galaxies selected with the traditional UVJ technique is only $\sim20\%$ complete at $z<3$, when recovering the rest-frame colours using mock Euclid observations. This shows that our new methodology is the most suitable one when only Euclid bands, along with u-band imaging, are available.

The synchrotron cooling of relativistic electrons is one of the most effective radiation mechanisms in astrophysics. It not only accompanies the process of particle acceleration but also has feedback on the formation of the energy distribution of the parent electrons. The radiative cooling time of electrons decreases with energy as $t_{\rm syn} \propto 1/E$; correspondingly the overall radiation efficiency increases with energy. On the other hand, this effect strictly limits the maximum energy of individual photons. Even in the so-called extreme accelerators, where the acceleration proceeds at the highest possible rate, $t_{\rm acc} = E/eBc$, the synchrotron radiation cannot extend well beyond the characteristic energy determined by the electron mass and the fine-structure constant: $h \nu_{\rm max} \sim m_e c^2/\alpha \sim 100 \rm\,MeV$. In a highly turbulent medium with a broad magnetic field distribution, this limit is shifted towards higher energies. However, for any realistic distribution of the magnetic field, it cannot significantly exceed $ h \nu_{\rm max}$. In this paper, we propose a model in which the formation of synchrotron radiation takes place in compact magnetic blobs located inside the particle accelerator. We develop a formalism for calculations of synchrotron radiation emerging from such systems. We demonstrate that for certain combinations of parameters characterizing the accelerator and the magnetic blobs, the synchrotron radiation can extend into the very high energy (VHE) domain, $h \nu_{\rm max} \geq 100 \ \rm GeV$. We briefly outline the possible astrophysical implications of the results in the context of VHE radiation of blazars and gamma-ray bursts.

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.

The most powerful persistent accelerators in the Universe are jetted active galaxies. Galaxies whose jets are directed towards Earth, so called blazars, dominate the extragalactic gamma-ray sky. Still, most of the highest-energy accelerators likely elude detection. These extreme blazars, whose radiated energy can peak beyond 10 TeV, are ideal targets to study particle acceleration and radiative processes, and may provide links to cosmic rays and astrophysical neutrinos. The growing number of extreme blazars observed at TeV energies has been critical for the emergence of gamma-ray cosmology, including measurements of the extragalactic background light, tight bounds on the intergalactic magnetic field, and constraints on exotic physics at energies inaccessible with human-made accelerators. Tremendous progress is expected in the decade to come, particularly with the deployment of the Cherenkov Telescope Array.

Anomalous Microwave Emission (AME) is a significant component of Galactic diffuse emission in the frequency range $10$-$60\,$GHz and a new window into the properties of sub-nanometre-sized grains in the interstellar medium. We investigate the morphology of AME in the $\approx10^{\circ}$ diameter $\lambda$ Orionis ring by combining intensity data from the QUIJOTE experiment at $11$, $13$, $17$ and $19\,$GHz and the C-Band All Sky Survey (C-BASS) at $4.76\,$GHz, together with 19 ancillary datasets between $1.42$ and $3000\,$GHz. Maps of physical parameters at $1^{\circ}$ resolution are produced through Markov Chain Monte Carlo (MCMC) fits of spectral energy distributions (SEDs), approximating the AME component with a log-normal distribution. AME is detected in excess of $20\,\sigma$ at degree-scales around the entirety of the ring along photodissociation regions (PDRs), with three primary bright regions containing dark clouds. A radial decrease is observed in the AME peak frequency from $\approx35\,$GHz near the free-free region to $\approx21\,$GHz in the outer regions of the ring, which is the first detection of AME spectral variations across a single region. A strong correlation between AME peak frequency, emission measure and dust temperature is an indication for the dependence of the AME peak frequency on the local radiation field. The AME amplitude normalised by the optical depth is also strongly correlated with the radiation field, giving an overall picture consistent with spinning dust where the local radiation field plays a key role.

Recent developments of transition-edge sensors (TESs), based on extensive experience in ground-based experiments, have been making the sensor techniques mature enough for their application on future satellite CMB polarization experiments. LiteBIRD is in the most advanced phase among such future satellites, targeting its launch in Japanese Fiscal Year 2027 (2027FY) with JAXA's H3 rocket. It will accommodate more than 4000 TESs in focal planes of reflective low-frequency and refractive medium-and-high-frequency telescopes in order to detect a signature imprinted on the cosmic microwave background (CMB) by the primordial gravitational waves predicted in cosmic inflation. The total wide frequency coverage between 34GHz and 448GHz enables us to extract such weak spiral polarization patterns through the precise subtraction of our Galaxy's foreground emission by using spectral differences among CMB and foreground signals. Telescopes are cooled down to 5Kelvin for suppressing thermal noise and contain polarization modulators with transmissive half-wave plates at individual apertures for separating sky polarization signals from artificial polarization and for mitigating from instrumental 1/f noise. Passive cooling by using V-grooves supports active cooling with mechanical coolers as well as adiabatic demagnetization refrigerators. Sky observations from the second Sun-Earth Lagrangian point, L2, are planned for three years. An international collaboration between Japan, USA, Canada, and Europe is sharing various roles. In May 2019, the Institute of Space and Astronautical Science (ISAS), JAXA selected LiteBIRD as the strategic large mission No. 2.

The space-borne gravitational wave interferometer, LISA, is expected to detect signals from numerous binary white dwarfs. At small orbital separation, rapid rotation and large tidal bulges may allow for the stellar internal structure to be probed through such observations. Finite-size effects are encoded in quantities like the moment of inertia ($I$), tidal Love number (Love), and quadrupole moment ($Q$). The universal relations among them (I-Love-Q relations) can be used to reduce the number of parameters in the gravitational-wave templates. We here study I-Love-Q relations for more realistic white dwarf models than used in previous studies. In particular, we extend previous works by including (i) differential rotation and (ii) internal temperature profiles taken from detailed stellar evolution calculations. We use the publicly available stellar evolution code MESA to generate cooling models of both low- and high-mass white dwarfs. We show that differential rotation causes the I-Q relation (and similarly the Love-Q relation) to deviate from that of constant rotation. We also find that the introduction of finite temperatures causes the white dwarf to move along the zero-temperature mass sequence of I-Q values, moving towards values that suggest a lower mass. We further find that after only a few Myrs, high-mass white dwarfs are well-described by the zero-temperature model, suggesting that the relations with zero-temperature may be good enough in most practical cases. Low-mass, He-core white dwarfs with thick hydrogen envelopes may undergo long periods of H burning which sustain the stellar temperature and allow deviations from the I-Love-Q relations for longer times.

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.

We study the properties of star-forming galaxies selected at 610 MHz with the GMRT in a survey covering $\sim$1.86 deg$^2$ down to a noise of $\sim$7.1\,$\mu$Jy / beam. These were identified by combining multiple classification diagnostics: optical, X-ray, infrared and radio data. Of the 1685 SFGs from the GMRT sample, 496 have spectroscopic redshifts whereas 1189 have photometric redshifts. We find that the IRRC of star-forming galaxies, quantified by the infrared-to-1.4 GHz radio luminosity ratio $\rm{q_{IR}}$, decreases with increasing redshift: $\rm{q_{IR}\,=\,2.86\pm0.04(1\,+\,z)^{-0.20\pm0.02}}$ out to $z \sim 1.8$. We use the $\rm{V/V_{max}}$ statistic to quantify the evolution of the co-moving space density of the SFG sample. Averaged over luminosity our results indicate $\rm{\langle V/V_{max} \rangle}$ to be $\rm{0.51\,\pm\, 0.06}$, which is consistent with no evolution in overall space density. However we find $\rm V/V_{max}$ to be a function of radio luminosity, indicating strong luminosity evolution with redshift. We explore the evolution of the SFGs radio luminosity function by separating the source into five redshift bins and comparing to theoretical model predictions. We find a strong redshift trend that can be fitted with a pure luminosity evolution of the form $\rm{L_{610\,MHz}\,\propto\,(\,1+\,z)^{(2.95\pm0.19)-(0.50\pm0.15)z}}$. We calculate the cosmic SFR density since $\rm{z \sim 1.5}$ by integrating the parametric fits of the evolved 610\,MHz luminosity function. Our sample reproduces the expected steep decline in the star formation rate density since $\rm{z\,\sim\,1}$.

We assess the ionising effect of low energy protostellar cosmic rays in protoplanetary disks around a young solar mass star for a wide range of disk parameters. We assume a source of low energy cosmic rays located close to the young star which travel diffusively through the protoplanetary disk. We use observationally inferred values from nearby star-forming regions for the total disk mass and the radial density profile. We investigate the influence of varying the disk mass within the observed scatter for a solar mass star. We find that for a large range of disk masses and density profiles that protoplanetary disks are "optically thin" to low energy ($\sim$3 GeV) cosmic rays. At $R\sim10$au, for all of the disks that we consider ($M_\mathrm{disk}=6.0\times10^{-4} - 2.4\times 10^{-2}M_\odot$), the ionisation rate due to low energy stellar cosmic rays is larger than that expected from unmodulated galactic cosmic rays. This is in contrast to our previous results which assumed a much denser disk which may be appropriate for a more embedded source. At $R\sim70$au, the ionisation rate due to stellar cosmic rays dominates in $\sim$50% of the disks. These are the less massive disks with less steep density profiles. At this radius there is at least an order of magnitude difference in the ionisation rate between the least and most massive disk that we consider. Our results indicate, for a wide range of disk masses, that low energy stellar cosmic rays provide an important source of ionisation at the disk midplane at large radii ($\sim$70au).

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.

We provide predictions of the yield of $7<z<9$ quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; $\Phi$) with redshift, $\Phi\propto10^{k(z-6)}$, $k=-0.72$, and a further steeper rate of decline, $k=-0.92$; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we use an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to $J_{AB}\sim23$. Quasars at $z>8$ may be selected from Euclid $OYJH$ photometry alone, but selection over the redshift interval $7<z<8$ is greatly improved by the addition of $z$-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond $z=6$. For the case that the decline of the QLF accelerates beyond $z=6$, with $k=-0.92$, Euclid should nevertheless find over 100 quasars with $7.0<z<7.5$, and $\sim25$ quasars beyond the current record of $z=7.5$, including $\sim8$ beyond $z=8.0$. The first Euclid quasars at $z>7.5$ should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, $7<z<8$, $M_{1450}<-25$, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at $J_{AB}\sim23$. The precision with which $k$ can be determined over $7<z<8$ depends on the value of $k$, but assuming $k=-0.72$ it can be measured to a 1 sigma uncertainty of 0.07.

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

This is the first of a series of papers based on sensitive 610 MHz observations of the ELAIS N1 field using the Giant Metrewave Radio Telescope. We describe the observations, processing and source catalogue extraction from a deep image with area of 1.86\,deg$^{2}$ and minimum noise of $\sim$7.1\,$\mu$Jy/beam. We compile a catalogue of 4290 sources with flux densities in the range 28.9\,$\mu$Jy - 0.503\,Jy, and derive the Euclidean-normalized differential source counts for sources with flux densities brighter than $\rm{35.5\,\mu\,Jy}$. Our counts show a flattening at 610 MHz flux densities below 1 mJy. Below the break the counts are higher than previous observations at this frequency, but generally consistent with recent models of the low-frequency source population. The radio catalogue is cross-matched against multi-wavelength data leading to identifications for 92\% and reliable redshifts for 72\% of our sample, with 19\% of the redshifts based on spectroscopy. For the sources with redshifts we use radio and X-ray luminosity, optical spectroscopy and mid-infrared colours to search for evidence of the presence of an Active Galactic Nucleus (AGN). We compare our identifications to predictions of the flux density distributions of star forming galaxies (SFGs) and AGN, and find a good agreement assuming the majority of the sources without redshifts are SFGs. We derive spectral index distributions for a sub-sample. The majority of the sources are steep spectra, with a median spectral index that steepens with frequency; $\mathrm{\alpha^{325}_{610}\,=\,-0.80\,\pm\,0.29}$, $\mathrm{\alpha^{610}_{1400}\,=\,-0.83\,\pm\,0.31}$ and $\mathrm{\alpha^{610}_{5000}\,=\,-1.12\,\pm\,0.15}$.

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.

In this work we use a multi-messenger approach to determine if the high energy diffuse neutrino flux observed by the IceCube Observatory can originate from $\gamma$-ray sources powered by Cosmic Rays interactions with gas. Typical representatives of such sources are Starburst and Ultra-Luminous Infrared Galaxies. Using the three most recent calculations of the non-blazar contribution to the extragalactic $\gamma$-ray background measured by the Fermi-LAT collaboration, we find that a hard power-law spectrum with spectral index $\alpha \leq 2.12$ is compatible with all the estimations for the allowed contribution from non-blazar sources, within 1$\sigma$. Using such a spectrum we are able to interpret the IceCube results, showing that various classes of hadronically powered $\gamma$-ray galaxies can provide the dominant contribution to the astrophysical signal above 100 TeV and about half of the contribution to the energy flux between 10-100 TeV. With the addition of neutrinos from the Galactic plane, it is possible to saturate the IceCube signal at high energy. Our result shows that these sources are still well motivated candidates.

Black holes are unique among astrophysical sources: they are the simplest macroscopic objects in the Universe, and they are extraordinary in terms of their ability to convert energy into electromagnetic and gravitational radiation. Our capacity to probe their nature is limited by the sensitivity of our detectors. The LIGO/Virgo interferometers are the gravitational-wave equivalent of Galileo's telescope. The first few detections represent the beginning of a long journey of exploration. At the current pace of technological progress, it is reasonable to expect that the gravitational-wave detectors available in the 2035-2050s will be formidable tools to explore these fascinating objects in the cosmos, and space-based detectors with peak sensitivities in the mHz band represent one class of such tools. These detectors have a staggering discovery potential, and they will address fundamental open questions in physics and astronomy. Are astrophysical black holes adequately described by general relativity? Do we have empirical evidence for event horizons? Can black holes provide a glimpse into quantum gravity, or reveal a classical breakdown of Einstein's gravity? How and when did black holes form, and how do they grow? Are there new long-range interactions or fields in our universe, potentially related to dark matter and dark energy or a more fundamental description of gravitation? Precision tests of black hole spacetimes with mHz-band gravitational-wave detectors will probe general relativity and fundamental physics in previously inaccessible regimes, and allow us to address some of these fundamental issues in our current understanding of nature.

The recent milestones in multi-messenger astronomy have opened new ways to study the Unverse. The unprecedented gravitational wave (GW) follow-up campaigns established the power that the combination of different messengers has to identify and study the nature and evolution of astrophysical phenomena. Here we focus on the search for high-energy gamma ray emission as electromagnetic counterpart of compact binary coalescences with the H.E.S.S. Imaging Air Cherenkov Telescopes (IACTs). In this contribution, the optimized strategies developed specifically for the prompt follow-up of gravitational wave events with H.E.S.S are presented. As illustration, the successful observation campaigns up to this time will be described, including the ones during Observation Run O2 on the binary black hole (BH-BH) merger GW170814 and the binary neutron star (NS-NS) merger GW170817, and an update on recent events occurring during O3. Results of these searches are presented and the constraints that prompt observations can put on very-high-energy, non-thermal emission, are briefly discussed. Finally, an outlook on further improvements for the gravitational waves follow-up program with H.E.S.S. will be provided.

Stabilizing the gain of a radio astronomy receiver is of great importance for sensitive radio intensity mapping. In this paper we discuss a stabilization method using a continuous-wave reference signal injected into the signal chain and tracked in a single channel of the spectrometer to correct for the gain variations of the receiver. This method depends on the fact that gain fluctuations of the receiver are strongly correlated across the frequency band, which we can show is the case for our experimental setup. This method is especially suited for receivers with a digital back-end with high spectral resolution and moderate dynamic range. The sensitivity of the receiver is unaltered except for one lost frequency channel. We present experimental results using a new 4-8.5 GHz receiver with a digital back-end that shows substantial reduction of the 1/ f noise and the 1/ f knee frequency.

We summarise the results found following a study of the effects that Galactic magnetic fields can have on the propagation of cosmic rays from local extragalactic sources. This study focuses on the coherent structures in the Jansson-Farrar Galactic magnetic field model, namely the: disk field, a toroidal field, and an x-field components. The phenomena of Galactic magnetic field shadowing, source deflection, and Galactic magnetic field tunnel vision are all noted. Attention throughout this study is placed on particles with rigidity around $10^{18.5}$ EV, believed to dominate the cosmic ray spectrum above the ankle. The Galactic magnetic field component predominantly responsible for giving rise to each of these different effects was determined.

The cosmic microwave background $B$-mode signal is potentially weaker than the diffuse Galactic foregrounds over most of the sky at any frequency. A common method of separating the CMB from these foregrounds is via pixel-based parametric-model fitting. There are not currently enough all-sky maps to fit anything more than the most simple models of the sky. By simulating the emission in seven representative pixels, we demonstrate that the inclusion of a 5 GHz data point allows for more complex models of low-frequency foregrounds to be fitted than at present. It is shown that the inclusion of the CBASS data will significantly reduce the uncertainties in a number of key parameters in the modelling of both the galactic foregrounds and the CMB. The extra data allow estimates of the synchrotron spectral index to be constrained much more strongly than is presently possible, with corresponding improvements in the accuracy of the recovery of the CMB amplitude. However, we show that to place good limits on models of the synchrotron spectral curvature will require additional low-frequency data.

In modern weak-lensing surveys, the common approach to correct for residual systematic biases in the shear is to calibrate shape measurement algorithms using simulations. These simulations must fully capture the complexity of the observations to avoid introducing any additional bias. In this paper we study the importance of faint galaxies below the observational detection limit of a survey. We simulate simplified Euclid VIS images including and excluding this faint population, and measure the shift in the multiplicative shear bias between the two sets of simulations. We measure the shear with three different algorithms: a moment-based approach, model fitting, and machine learning. We find that for all methods, a spatially uniform random distribution of faint galaxies introduces a shear multiplicative bias of the order of a few times $10^{-3}$. This value increases to the order of $10^{-2}$ when including the clustering of the faint galaxies, as measured in the Hubble Space Telescope Ultra-Deep Field. The magnification of the faint background galaxies due to the brighter galaxies along the line of sight is found to have a negligible impact on the multiplicative bias. We conclude that the undetected galaxies must be included in the calibration simulations with proper clustering properties down to magnitude 28 in order to reach a residual uncertainty on the multiplicative shear bias calibration of a few times $10^{-4}$, in line with the $2\times10^{-3}$ total accuracy budget required by the scientific objectives of the Euclid survey. We propose two complementary methods for including faint galaxy clustering in the calibration simulations.

Much data on the Galactic polarized radio emission has been gathered in the last five decades. All-sky surveys have been made, but only in narrow, widely spaced frequency bands, and the data are inadequate for the characterization of Faraday rotation, the main determinant of the appearance of the polarized radio sky at decimetre wavelengths. We describe a survey of the polarized radio emission from the Southern sky, aiming to characterize the magneto-ionic medium, particularly the strength and configuration of the magnetic field. This work is part of the Global Magneto-Ionic Medium Survey (GMIMS). We have designed and built a feed and receiver covering the band 300 to 900 MHz for the CSIRO Parkes 64-m Telescope. We have surveyed the entire sky between declinations -90 and +20 degrees. We present data covering 300 to 480 MHz with angular resolution 81' to 45'. The survey intensity scale is absolutely calibrated, based on measurements of resistors at known temperatures and on an assumed flux density and spectral index for Taurus A. Data are presented as brightness temperatures. We have applied Rotation Measure Synthesis to the data to obtain a Faraday depth cube of resolution 5.9 radians per metre squared, sensitivity of 60 mK of polarized intensity, and angular resolution 1.35 degrees. The data presented in this paper are available at the Canadian Astronomy Data Centre.

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

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

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

We present the first joint inference of standard sirens and gravitational wave weak lensing by filtering of the same dataset. We imagine a post-LISA scenario emerging around the late 2030s when LISA will have accumulated a number of detections at high redshift; LIGO-VIRGO will have finished observing at low redshift, and Einstein Telescope will have started making new observations out to redshifts possibly overlapping with LISA. Euclid and other cosmological probes will have provided constraints at percent level by then, but mostly exhausted their ability to improve any further. We derive forecasts assuming $\sim1\,\text{deg}^{-2}$ detected sources, in conjunction with a spectroscopic follow-up (e.g. by Euclid, DESI, or ATHENA). Thanks to the statistical power of standard sirens as a geometry probe -- lifting key degeneracies in the gravitational wave weak lensing -- and no external priors assumed, the constraints on dark matter and its clustering, namely $\Omega_m$ and $\sigma_8$, could be achieved to $2\%$ and $3\%$. The Hubble constant could be constrained to better than $1\%$ in all cases; the dark energy density, $\Omega_\Lambda$, to $2\%$, and curvature, $\Omega_K$, to $0.02$; the amplitude and spectral tilt of the scalar fluctuations, $\ln(10^{10}A_s)$ and $n_s$, to $2\%$ and $7\%$. As a completely independent cosmological probe, with less calibration requirements, the joint inference of standard sirens and gravitational wave weak lensing might help solve the tensions currently observed between other cosmological probes, such as CMB, galaxy lensing, and Type Ia SNs, and distinguish between residual systematics and new physics.

Sub-Gauss magnetic fields are crucial for many physical processes from Galactic $kpc$ scale to stellar $Au$ scale, but are nontrivial to detect. All present magnetic diagnoses trace only one component of a magnetic field at best. Here we report the first observational results that unveil 3D topology as well as the strength of a sub-milliGauss magnetic field beyond our solar system. We found that two weak neutral iron absorption lines from the ground state in the atmosphere of $89~Her$ produced counterintuitive high-amplitude polarizations and consistent polarization angles, exactly in line with the theoretical prediction from ground state alignment (GSA). Our analysis reveals the first sub-AU scale magnetic field on $89~Her$, that is $1.3kpc$ away, has a 3D orientation aligned to the stellar outflow axis and a strength of $70\mu{G} \lesssim B \lesssim150\mu{G}$, thus substantially improving the accuracy by five orders of magnitude compared to the previous $10G$ upper limit set by non-detection of Zeeman effect. This long-awaited discovery is anticipated to usher in a new era of precise measurement of sub-gauss magnetic field in the Universe.

We review open questions and prospects for progress in ultrahigh-energy cosmic ray (UHECR) research, based on a series of discussions that took place during the `The High-Energy Universe: Gamma-Ray, Neutrino, and Cosmic-ray Astronomy' MIAPP workshop in 2018. Specifically, we overview open questions on the origin of the bulk of UHECRs, the UHECR mass composition, the origin of the end of the cosmic-ray spectrum, the transition from Galactic to extragalactic cosmic rays, the effect of magnetic fields on the trajectories of UHECRs, anisotropy expectations for specific astrophysical scenarios, hadronic interactions, and prospects for discovering neutral particles as well as new physics at ultrahigh energies. We also briefly present upcoming and proposed UHECR experiments and discuss their projected science reach.

We have previously presented a reconstruction of Horndeski scalar-tensor theories from linear cosmological observables. It includes free nonlinear terms which can be added onto the reconstructed covariant theory without affecting the background and linear dynamics. After discussing the uniqueness of these correction terms, we apply this nonlinear freedom to a range of different applications. First we demonstrate how the correction terms can be configured to endow the reconstructed models with screening mechanisms such as the chameleon, k-mouflage and Vainshtein effects. A further implication is the existence of classes of Horndeski models that are degenerate with standard cosmology to an arbitrary level in the cosmological perturbations. Particularly interesting examples are kinetically self-accelerating models that mimic the dynamics of the cosmological constant to an arbitrary degree in perturbations. Finally, we develop the reconstruction method further to the level of higher-order effective field theory, which under the restriction to a luminal propagation speed of gravitational waves introduces two new free functions per order. These functions determine the corresponding correction terms in the linearly reconstructed action at the same order. Our results enable the connection of linear cosmological constraints on generalised modifications of gravity and dark energy with the nonlinear regime and astrophysical probes for a more global interpretation of the wealth of forthcoming cosmological survey data.

We outline the science motivation for SGSO, the Southern Gamma-Ray Survey Observatory. SGSO will be a next-generation wide field-of-view gamma-ray survey instrument, sensitive to gamma-rays in the energy range from 100 GeV to hundreds of TeV. Its science topics include unveiling galactic and extragalactic particle accelerators, monitoring the transient sky at very high energies, probing particle physics beyond the Standard Model, and the characterization of the cosmic ray flux. SGSO will consist of an air shower detector array, located in South America. Due to its location and large field of view, SGSO will be complementary to other current and planned gamma-ray observatories such as HAWC, LHAASO, and CTA.

Although many high-energy neutrinos detected by the IceCube telescope are believed to have anextraterrestrial origin, their astrophysical sources remain a mystery. Recently, an unprecedenteddiscovery of a high-energy muon neutrino event coincident with a multiwavelength flare from ablazar, TXS 0506+056, shed some light on the origin of the neutrinos. It is usually believed that ablazar is produced by a relativistic jet launched from an accreting supermassive black hole (SMBH).Here we show that the high-energy neutrino event can be interpreted by the inelastic hadronuclearinteractions between the accelerated cosmic-ray protons in the relativistic jet and the dense gasclouds in the vicinity of the SMBH. Such a scenario only requires a moderate proton power in thejet, which could be much smaller than that required in the conventional hadronic model whichinstead calls upon the photomeson process. Meanwhile, the flux of the multiwavelength flare fromthe optical to gamma-ray band can be well explained by invoking a second radiation zone in thejet at a larger distance to the SMBH. In our model, the neutrino emission lasts a shorter time thanthe multiwavelength flare so the neutrino event is not necessarily correlated with the flare but it is probably accompanied by a spectrum hardening above a few GeV.

The C-Band All-Sky Survey C-BASS is a high-sensitivity all-sky radio survey at an angular resolution of 45 arcmin and a frequency of 4.7 GHz. We present a total intensity 4.7 GHz map of the North Celestial Pole (NCP) region of sky, above declination +80 deg, which is limited by source confusion at a level of ~0.6 mK rms. We apply the template-fitting (cross-correlation) technique to WMAP and Planck data, using the C-BASS map as the synchrotron template, to investigate the contribution of diffuse foreground emission at frequencies ~20-40 GHz. We quantify the anomalous microwave emission (AME) that is correlated with far-infrared dust emission. The AME amplitude does not change significantly (<10%) when using the higher frequency C-BASS 4.7 GHz template instead of the traditional Haslam 408 MHz map as a tracer of synchrotron radiation. We measure template coefficients of $9.93\pm0.35$ and $9.52\pm0.34$ K per unit $\tau_{353}$ when using the Haslam and C-BASS synchrotron templates, respectively. The AME contributes $55\pm2\,\mu$K rms at 22.8 GHz and accounts for ~60% of the total foreground emission. Our results suggest that a harder (flatter spectrum) component of synchrotron emission is not dominant at frequencies >5 GHz; the best-fitting synchrotron temperature spectral index is $\beta=-2.91\pm0.04$ from 4.7 to 22.8 GHz and $\beta=-2.85\pm0.14$ from 22.8 to 44.1 GHz. Free-free emission is weak, contributing ~$7\,\mu$K rms (~7%) at 22.8 GHz. The best explanation for the AME is still electric dipole emission from small spinning dust grains.

Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3\pi$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$\sigma$ point source limiting sensitivities in the stacked 3$\pi$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ \Delta ra, \Delta dec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI.

The C-Band All-Sky Survey (C-BASS) is an all-sky full-polarization survey at a frequency of 5 GHz, designed to provide data complementary to the all-sky surveys of WMAP and Planck and future CMB B-mode polarization imaging surveys. We describe the design and performance of the digital backend used for the northern part of the survey. In particular we describe the features that efficiently implement the demodulation and filtering required to suppress contaminating signals in the time-ordered data, and the capability for real-time correction of detector non-linearity and receiver balance.

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.

Flux distribution is an important tool to understand the variability processes in active galactic nuclei. We now have available a great deal of observational evidences pointing towards the presence of log-normal components in the high energy light curves, and different models have been proposed to explain these data. Here, we collect some of the recent developments on this topic using the well-known blazar Mrk 501 as example of complex and interesting aspects coming from its flux distribution in different energy ranges and at different timescales. The observational data we refer to are those collected in a complementary manner by Fermi-LAT over multiple years, and by the First G-APD Cherenkov Telescope (FACT) telescope and the H.E.S.S. array in correspondence of the bright flare of June 2014.

Accurate and precise covariance matrices will be important in enabling planned cosmological surveys to detect new physics. Standard methods imply either the need for many N-body simulations in order to obtain an accurate estimate, or a precise theoretical model. We combine these approaches by constructing a likelihood function conditioned on simulated and theoretical covariances, consistently propagating noise from the finite number of simulations and uncertainty in the theoretical model itself using an informative Inverse-Wishart prior. Unlike standard methods, our approach allows the required number of simulations to be less than the number of summary statistics. We recover the linear 'shrinkage' covariance estimator in the context of a Bayesian data model, and test our marginal likelihood on simulated mock power spectrum estimates. We conduct a thorough investigation into the impact of prior confidence in different choices of covariance models on the quality of model fits and parameter variances. In a simplified setting we find that the number of simulations required can be reduced if one is willing to accept a mild degradation in the quality of model fits, finding that even weakly informative priors can help to reduce the simulation requirements. We identify the correlation matrix of the summary statistics as a key quantity requiring careful modelling. Our approach can be easily generalized to any covariance model or set of summary statistics, and elucidates the role of hybrid estimators in cosmological inference.

We study the cosmic evolution of the magnetic fields of a large sample of spiral galaxies in a cosmologically representative volume by employing a semi-analytic galaxy formation model and numerical dynamo solver in tandem. We start by deriving time- and radius-dependent galaxy properties using the galform galaxy formation model, which are then fed into the nonlinear mean-field dynamo equations. These are solved to give the large-scale (mean) field as a function of time and galactocentric radius for a thin disc, assuming axial symmetry. A simple prescription for the evolution of the small-scale (random) magnetic field component is also adopted. We find that, while most massive galaxies are predicted to have large-scale magnetic fields at redshift z=0, a significant fraction of them are expected to contain negligible large-scale field. Our model indicates that, for most of the galaxies containing large-scale magnetic fields today, the mean-field dynamo becomes active at z<3. We compute the radial profiles of pitch angle, and find broad agreement with observational data for nearby galaxies.

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.

We investigate non-thermal electron and nuclei energy losses within the binary neutron star merger remnant produced by the event GW170817. The lack of a cooling feature within the detected synchrotron emission from the source is used to constrain the magnetic field at the mG level, assuming that this emission is electron synchrotron in origin, and that the accelerated spectrum in the electrons follows the form $dN/dE_e \propto E_e^{-2}$. The level of subsequent gamma-ray emission from the source is demonstrated to provide a further constraint on the source magnetic field strength. We also put forward alternative strong ($\sim$G) magnetic field scenarios able to support this emission. For such stronger fields, the photo-disintegration of non-thermal nuclei within the source is considered, and a bottleneck period of $\sim$5-30 days is found when this process peaks. We find that this class of source is in principle able to support the population of cosmic rays detected at Earth below the "ankle".

We argue that pitch angles of the azimuthally averaged large-scale or mean magnetic fields in nearby spiral galaxies inferred from observations can tentatively be explained with simple galactic dynamo models. Agreement is not perfect, but is reasonable considering the uncertainty in dynamo parameters.

The growing wealth of cosmological observations places increasingly more stringent constraints on dark energy and alternative gravity models. Particularly successful in efficiently probing the vast model space has been the effective field theory of dark energy and modified gravity, providing a unified framework for generalised cosmological predictions. However, the optimal parametrisation of the free time-dependent functions inherent to the formalism is still unresolved. It should respect a multitude of requirements, ranging from simplicity, generality, and representativity of known theories to computational efficiency. But in particular, for theoretical soundness, the parameter space should adhere to strict stability requirements. We have recently proposed an inherently stable effective field theory with physical basis of Planck mass evolution, sound speed of the scalar field fluctuation, kinetic coefficient, and background expansion, which covers Horndeski models with luminal speed of gravity. Here we devise a parametrisation of these basis functions that can straightforwardly be configured to evade theoretical pathologies such as ghost or gradient instabilities or to accommodate further theoretical priors such as (sub)luminal scalar sound speed. The parametrisation is simple yet general, conveniently represents a broad range of known dark energy and gravitational theories, and with a simple additional approximation can be rendered numerically highly efficient. Finally, by operating in our new basis, we show that there are no general limitations from stability requirements on the current values that can be assumed by the phenomenological modification of the Poisson equation and the gravitational slip besides the exclusion of anti-gravity. The inherently stable effective field theory is ready for implementation in parameter estimation analyses employing linear cosmological observations.

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 investigate the use of estimators of weak lensing power spectra based on a flat-sky implementation of the Pseudo-Cl (PCl) technique, where the masked shear field is transformed without regard for masked regions of sky. This masking mixes power, and E-convergence and B-modes. To study the accuracy of forward-modelling and full-sky power spectrum recovery we consider both large-area survey geometries, and small-scale masking due to stars and a checkerboard model for field-of-view gaps. The power spectrum for the large-area survey geometry is sparsely-sampled and highly oscillatory, which makes modelling problematic. Instead, we derive an overall calibration for large-area mask bias using simulated fields. The effects of small-area star masks can be accurately corrected for, while the checkerboard mask has oscillatory and spiky behaviour which leads to percent biases. Apodisation of the masked fields leads to increased biases and a loss of information. We find that we can construct an unbiased forward-model of the raw PCls, and recover the full-sky convergence power to within a few percent accuracy for both Gaussian and lognormal-distributed shear fields. Propagating this through to cosmological parameters using a Fisher-Matrix formalism, we find we can make unbiased estimates of parameters for surveys up to 1,200 deg$^2$ with 30 galaxies per arcmin$^2$, beyond which the percent biases become larger than the statistical accuracy. This implies a flat-sky PCl analysis is accurate for current surveys but a Euclid-like survey will require higher accuracy.

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.

DG Tau A, a class-II young stellar object (YSO) displays both thermal, and non-thermal, radio emission associated with its bipolar jet. To investigate the nature of this emission, we present sensitive ($\sigma\sim2\,\mu {\rm Jy \,beam^{-1}}$), Karl G.\ Jansky Very Large Array (VLA) $6$ and $10\,{\rm GHz}$ observations. Over $3.81\,{\rm yr}$, no proper motion is observed towards the non-thermal radio knot C, previously thought to be a bowshock. Its quasi-static nature, spatially-resolved variability and offset from the central jet axis supports a scenario whereby it is instead a stationary shock driven into the surrounding medium by the jet. Towards the internal working surface, knot A, we derive an inclination-corrected, absolute velocity of $258\pm23\, {\rm km\,s^{-1}}$. DG Tau A's receding counterjet displays a spatially-resolved increase in flux density, indicating a variable mass loss event, the first time such an event has been observed in the counterjet. For this ejection, we measure an ionised mass loss rate of $(3.7\pm1.0) \times 10^{-8}\, {\rm M_\odot\, yr^{-1}}$ during the event. A contemporaneous ejection in the approaching jet isn't seen, showing it to be an asymmetric process. Finally, using radiative transfer modelling, we find that the extent of the radio emission can only be explained with the presence of shocks, and therefore reionisation, in the flow. Our modelling highlights the need to consider the relative angular size of optically thick, and thin, radio emission from a jet, to the synthesised beam, when deriving its physical conditions from its spectral index.

Forthcoming large-scale surveys will soon attempt to measure cosmic shear to an unprecedented level of accuracy, requiring a similarly high level of accuracy in the shear measurements of galaxies. Factors such as pixelisation, imperfect point-spread function (PSF) correction, and pixel noise can all directly or indirectly lead to biases in shear measurements, and so it can be necessary for shear measurement methods to be calibrated against internal, external, or simulated data to minimize bias. It is thus important to understand the nature of this calibration. In this paper, we show that a typical calibration procedure will on average leave no residual additive bias, but will leave a residual multiplicative bias. Additionally, the errors on the post-calibration bias parameters will be changed, and on average increased, from the errors on the pre-calibration measurements of these parameters, but that this is generally worth the benefit in decreasing the expected value of the multiplicative bias. We find that in most typical cases, it is worthwhile to apply a first-order bias correction, while a higher-order bias correction is only worthwhile for methods with intrinsically high multiplicative bias ($>10$ per cent) or when the simulation size is very small ($<10^6$ simulated galaxies).

Baryonic effects on large-scale structure, if not accounted for, can significantly bias dark energy constraints. As the detailed physics of the baryons is not yet well-understood, correcting for baryon effects introduces additional parameters which must be marginalized over, increasing the uncertainties on the inferred cosmological parameters. Forthcoming weak lensing surveys are aiming for percent-level precision on the dark energy equation of state, so the problem must be thoroughly examined. We use a halo model with analytic modifications which capture the impact of adiabatic contraction of baryons and feedback on the matter power spectrum, and generalize the Navarro-Frenk-White profile to account for a possible inner core. A Fisher analysis predicts degradations of 40% in the $w_0$-$w_a$ Figure of Merit for a Euclid-like survey, and up to 80% for other cosmological parameters. We forecast potential inner core constraints of a few $\mathrm{kpc}$, while for a fixed inner core, adiabatic concentration and feedback parameters are constrained to a few percent. We explore the scales where baryons and dark energy contribute most to the Fisher information, finding that probing to increasingly non-linear scales does little to reduce degradation. Including external baryon information improves our forecasts, but limiting degradation to 1% requires strong priors. Adding Planck cosmic microwave background priors improves the Figure of Merit by a factor of 2.7 and almost completely recovers the individual marginalized errors on $w_0$ and $w_a$. We also quantify the calibration of baryon modelling required to reduce biases of dark energy forecasts to acceptable levels for forthcoming lensing surveys.

Upcoming surveys with new radio observatories such as the Square Kilometer Array will generate a wealth of imaging data containing large numbers of radio galaxies. Different classes of radio galaxies can be used as tracers of the cosmic environment, including the dark matter density field, to address key cosmological questions. Classifying these galaxies based on morphology is thus an important step toward achieving the science goals of next generation radio surveys. Radio galaxies have been traditionally been classified as Fanaroff-Riley (FR) I and II, although some exhibit more complex 'bent' morphologies arising from environmental factors or intrinsic properties. In this work we present the FIRST Classifier, an on-line system for automated classification of Compact and Extended radio sources. We developed the FIRST Classifier based on a trained Deep Convolutional Neural Network Model to automate the morphological classification of compact and extended radio sources observed in the FIRST radio survey. Our model achieved an overall accuracy of 97% and a recall of 98%, 100%, 98% and 93% for Compact, BENT, FRI and FRII galaxies respectively. The current version of the FIRST classifier is able to predict the morphological class for a single source or for a list of sources as Compact or Extended (FRI, FRII and BENT).

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 C-Band All-Sky Survey (C-BASS) is an all-sky full-polarisation survey at a frequency of 5 GHz, designed to provide complementary data to the all-sky surveys of WMAP and Planck, and future CMB B-mode polarization imaging surveys. The observing frequency has been chosen to provide a signal that is dominated by Galactic synchrotron emission, but suffers little from Faraday rotation, so that the measured polarization directions provide a good template for higher frequency observations, and carry direct information about the Galactic magnetic field. Telescopes in both northern and southern hemispheres with matched optical performance are used to provide all-sky coverage from a ground-based experiment. A continuous-comparison radiometer and a correlation polarimeter on each telescope provide stable imaging properties such that all angular scales from the instrument resolution of 45 arcmin up to full sky are accurately measured. The northern instrument has completed its survey and the southern instrument has started observing. We expect that C-BASS data will significantly improve the component separation analysis of Planck and other CMB data, and will provide important constraints on the properties of anomalous Galactic dust and the Galactic magnetic field.

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.

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

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.

The C-Band All-Sky Survey (C-BASS) is an experiment to image the whole sky in intensity and polarization at 5 GHz. The primary aim of C-BASS is to provide low-frequency all-sky maps of the Galactic emission which will enable accurate component separation analysis of both existing and future CMB intensity and polarization imaging surveys. Here we present an overview of the experiment and an update on the current status of observations. We present simulation results showing the expected improvement in the recovery of CMB and foreground signals when including C-BASS data as an additional low-frequency channel, both for intensity and polarization. We also present preliminary results from the northern part of the sky survey.

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.

Recently we have derived a set of mapping relations that enables the reconstruction of the family of Horndeski scalar-tensor theories which reproduce the background dynamics and linear perturbations of a given set of effective field theory of dark energy coefficients. In this paper we present a number of applications of this reconstruction. We examine the form of the underlying theories behind different phenomenological parameterizations of modified gravity and dark energy used in the literature, as well as examine theories that exhibit weak gravity, linear shielding, and minimal self-acceleration. Finally, we propose a new inherently stable parametrization basis for modified gravity and dark energy models.

Recent advances in the $\gamma$-ray observations of solar flares by the \textit{Fermi} satellite, demand revisions in the hadronic $\gamma$-ray flux computation below 1 GeV. In this work we utilize recently updated pion production cross sections, along with an accurate description of low energy nuclear interactions. Applying these new interaction descriptions to model the \textit{Fermi} Large Area Telescope (LAT) solar flare data, we infer the primary particle spectral parameters. Application of this new cross section description leads to significantly different spectral parameters compared to those obtained previously. Furthermore, the inclusion of nuclei in these calculations leads to a primary spectrum that is generally harder than that required from proton only considerations. Lastly, the flare data at lower MeV energies, detected by the \textit{Fermi} Gamma-ray Burst Monitor (GBM), is demonstrated to provide additional low-energy spectral information.

Anomalous Microwave Emission (AME) is a component of diffuse Galactic radiation observed at frequencies in the range $\approx 10$-60 GHz. AME was first detected in 1996 and recognised as an additional component of emission in 1997. Since then, AME has been observed by a range of experiments and in a variety of environments. AME is spatially correlated with far-IR thermal dust emission but cannot be explained by synchrotron or free-free emission mechanisms, and is far in excess of the emission contributed by thermal dust emission with the power-law opacity consistent with the observed emission at sub-mm wavelengths. Polarization observations have shown that AME is very weakly polarized ($\lesssim 1$%). The most natural explanation for AME is rotational emission from ultra-small dust grains ("spinning dust"), first postulated in 1957. Magnetic dipole radiation from thermal fluctuations in the magnetization of magnetic grain materials may also be contributing to the AME, particularly at higher frequencies ($\gtrsim 50$ GHz). AME is also an important foreground for Cosmic Microwave Background analyses. This paper presents a review and the current state-of-play in AME research, which was discussed in an AME workshop held at ESTEC, The Netherlands, June 2016.