
A new method providing general consistency constraints for
BeyondtheStandardModel (BSM) theories, using measurements at particle
colliders, is presented. The method, `Constraints On New Theories Using Rivet',
Contur, exploits the fact that particlelevel differential measurements made in
fiducial regions of phasespace have a high degree of modelindependence. These
measurements can therefore be compared to BSM physics implemented in Monte
Carlo generators in a very generic way, allowing a wider array of final states
to be considered than is typically the case. The Contur approach should be seen
as complementary to the discovery potential of direct searches, being designed
to eliminate inconsistent BSM proposals in a context where many (but perhaps
not all) measurements are consistent with the Standard Model. We demonstrate,
using a competitive simplified dark matter model, the power of this approach.
The Contur method is highly scaleable to other models and future measurements.

Cosmicray observations provide a powerful probe of dark matter annihilation
in the Galaxy. In this paper we derive constraints on heavy dark matter from
the recent precise AMS02 antiproton data. We consider all possible
annihilation channels into pairs of standard model particles. Furthermore, we
interpret our results in the context of minimal dark matter, including
higgsino, wino and quintuplet dark matter. We compare the cosmicray antiproton
limits to limits from $\gamma$ray observations of dwarf spheroidal galaxies
and to limits from $\gamma$ray and $\gamma$line observations towards the
Galactic center. While the latter limits are highly dependent on the dark
matter density distribution and only exclude a thermal wino for cuspy profiles,
the cosmicray limits are more robust, strongly disfavoring the thermal wino
dark matter scenario even for a conservative estimate of systematic
uncertainties.

Our Galactic Center, Sagittarius A* (Sgr A*), is believed to harbour a
supermassive black hole (BH), as suggested by observations tracking individual
orbiting stars. Upcoming submillimetre verylongbaselineinterferometry
(VLBI) images of Sgr A* carried out by the EventHorizonTelescope
Collaboration (EHTC) are expected to provide critical evidence for the
existence of this supermassive BH. We assess our present ability to use EHTC
images to determine if they correspond to a Kerr BH as predicted by Einstein's
theory of general relativity (GR) or to a BH in alternative theories of
gravity. To this end, we perform generalrelativistic magnetohydrodynamical
(GRMHD) simulations and use generalrelativistic radiative transfer (GRRT)
calculations to generate synthetic shadow images of a magnetised accretion flow
onto a Kerr BH. In addition, and for the first time, we perform GRMHD
simulations and GRRT calculations for a dilaton BH, which we take as a
representative solution of an alternative theory of gravity. Adopting the VLBI
configuration from the 2017 EHTC campaign, we find that it could be extremely
difficult to distinguish between BHs from different theories of gravity, thus
highlighting that great caution is needed when interpreting BH images as tests
of GR.

Observational efforts to image the immediate environment of a black hole at
the scale of the event horizon benefit from the development of efficient
imaging codes that are capable of producing synthetic data, which may be
compared with observational data. We aim to present RAPTOR, a new public code
that produces accurate images, animations, and spectra of relativistic plasmas
in strong gravity by numerically integrating the equations of motion of light
rays and performing timedependent radiative transfer calculations along the
rays. The code is compatible with any analytical or numerical spacetime. It is
hardwareagnostic and may be compiled and run both on GPUs and CPUs. We
describe the algorithms used in RAPTOR and test the code's performance. We have
performed a detailed comparison of RAPTOR output with that of other
radiativetransfer codes and demonstrate convergence of the results. We then
applied RAPTOR to study accretion models of supermassive black holes,
performing timedependent radiative transfer through general relativistic
magnetohydrodynamical (GRMHD) simulations and investigating the expected
observational differences between the socalled fastlight and slowlight
paradigms. Using RAPTOR to produce synthetic images and light curves of a GRMHD
model of an accreting black hole, we find that the relative difference between
fastlight and slowlight light curves is less than 5%. Using two distinct
radiativetransfer codes to process the same data, we find integrated flux
densities with a relative difference less than 0.01%. For twodimensional GRMHD
models, such as those examined in this paper, the fastlight approximation
suffices as long as errors of a few percent are acceptable. The convergence of
the results of two different codes demonstrates that they are, at a minimum,
consistent.

The first gravitational wave detections of mergers between black holes and
neutron stars represent a remarkable new regime of highenergy transient
astrophysics. The signals observed with LIGOVirgo detectors come from mergers
of extreme physical objects which are the end products of stellar evolution
within close binary systems. To better understand their origin and merger
rates, we have performed binary population syntheses at different metallicities
using the new stellar grid based populations synthesis code ComBinE. Starting
from newborn pairs of stars, we follow their evolution including mass loss,
mass transfer and accretion, common envelopes and supernova explosions. We
apply the binding energies of common envelopes based on dense grids of detailed
stellar structure models, make use of improved investigations of the subsequent
Case BB Rochelobe overflow and scale supernova kicks according to the
stripping of the exploding stars. We demonstrate that all the double black hole
mergers, GW150914, LVT151012, GW151226, GW170104, GW170608 and GW170814, as
well as the recently reported double neutron star merger GW170817, are
accounted for in our models in the appropriate metallicity regime. Our binary
interaction parameters are calibrated in particular to match the accurately
determined properties of Galactic double neutron star systems, and we discuss
their masses and types of supernova origin. Using our default values for the
input physics parameters, we find a double neutron star merger rate of about
3.0 Myr^1 for MilkyWay equivalent galaxies. Our absolute upper limit to the
mergerrate density of double neutron star systems is R=400 yr^1 Gpc^3 in the
local Universe (z=0).

We devise a new userfriendly tool interfaced with the Boltzmann code CLASS
to deal with any kind of exotic electromagnetic energy injection in the
universe and its impact on anisotropies of the Cosmic Microwave Background. It
makes use of the results from standard electromagnetic cascade calculations
develop in the context of WIMP annihilation, generalized to incorporate any
injection history. We first validate it on a specific WIMP scenario, the Higgs
Portal model, confirming that the standard effective onthespot treatment is
accurate enough. We then analyze the more involved example of evaporating
Primordial Black Holes (PBHs) with masses in the range
$[3\times10^{13},5\times10^{16}]$g, for which the standard approximations break
down. We derive robust CMB bounds on the relic density of evaporating PBHs,
ruling out the possibility for PBHs with a monochromatic distribution of masses
in the range $[3\times10^{13},2.5\times10^{16}]$g to represent all of the Dark
Matter in our Universe. Remarkably, we confirm with an accurate study that the
CMB bounds are several orders of magnitude stronger than those from the
galactic gammaray background in the range $[3\times10^{13},3\times10^{14}]$g.
A future CMB experiment like CORE+, or an experiment attempting at measuring
the 21 cm signal from the Dark Ages could greatly improve the sensitivity to
these models.

We apply the algorithm published by Liang et al. (2014) to describe the
Double Pulsar system J0737$$3039 and extract the sense of rotation of first
born recycled pulsar PSR J0737$$3039A. We find that this pulsar is rotating
prograde in its orbit. This is the first direct measurement of the sense of
rotation of a pulsar with respect to its orbit and a direct confirmation of the
rotating lighthouse model for pulsars. This result confirms that the spin
angular momentum vector is closely aligned with the orbital angular momentum,
suggesting that kick of the supernova producing the second born pulsar
J0737$$3039B was small.

Polarised radio emission from PSR J17452900 has already been used to
investigate the strength of the magnetic field in the Galactic Centre, close to
Sagittarius A*. Here we report how persistent radio emission from this
magnetar, for over four years since its discovery, has revealed large changes
in the observed Faraday rotation measure, by up to 3500 rad m$^{2}$ (a five
per cent fractional change). From simultaneous analysis of the dispersion
measure, we determine that these fluctuations are dominated by variations in
the projected magnetic field, rather than the integrated free electron density,
along the changing line of sight to the rapidly moving magnetar. From a
structure function analysis of rotation measure variations, and a recent epoch
of rapid change of rotation measure, we determine a minimum scale of magnetic
fluctuations of size ~ 2 au at the Galactic Centre distance, inferring PSR
J17452900 is just ~ 0.1 pc behind an additional scattering screen.

Observations of cosmic rays are a sensitive probe of dark matter annihilation
in our Galaxy. In this article we present an analysis of the AMS02 antiproton
data, reducing cosmicray propagation uncertainties by fitting at the same time
dark matter and propagation parameters. The result exhibits a possible hint for
dark matter pointing to an annihilation cross section close to the thermal
value. We investigate the compatibility of this signal with a dark matter
interpretation of the Galactic center excess seen in the FermiLAT gammaray
data and discuss implications for dark matter models.

Pulsar timing and gravitationalwave (GW) detectors are superb laboratories
to study gravity theories in the strongfield regime. Here we combine those
tools to test the monoscalartensor theory of Damour and EspositoFar{\`e}se
(DEF), which predicts nonperturbative scalarization phenomena for neutron stars
(NSs). First, applying Markovchain Monte Carlo techniques, we use the absence
of dipolar radiation in the pulsartiming observations of five binary systems
composed of a NS and a white dwarf, and eleven equations of state (EOSs) for
NSs, to derive the most stringent constraints on the two free parameters of the
DEF scalartensor theory. Since the binarypulsar bounds depend on the NS mass
and the EOS, we find that current pulsartiming observations leave
scalarization windows, i.e., regions of parameter space where scalarization can
still be prominent. Then, we investigate if these scalarization windows could
be closed and if pulsartiming constraints could be improved by
laserinterferometer GW detectors, when spontaneous (or dynamical)
scalarization sets in during the early (or late) stages of a binary NS (BNS)
evolution. For the early inspiral of a BNS carrying constant scalar charge, we
employ a Fisher matrix analysis to show that Advanced LIGO can improve
pulsartiming constraints for some EOSs, and nextgeneration detectors, such as
the Cosmic Explorer and Einstein Telescope, will be able to improve those
bounds for all eleven EOSs. Using the late inspiral of a BNS, we estimate that
for some of the EOSs under consideration the onset of dynamical scalarization
can happen early enough to improve the constraints on the DEF parameters
obtained by combining the five binary pulsars. Thus, in the near future the
complementarity of pulsar timing and direct observations of GWs on the ground
will be extremely valuable in probing gravity theories in the strongfield
regime.

We present first results from pilot observations using a phased array feed
(PAF) mounted on the Parkes 64m radio telescope. The observations presented
here cover a frequency range from 1150 to 1480 MHz and are used to show the
ability of PAFs to suppress standing wave problems by a factor of $\sim10$
which afflict normal feeds. We also compare our results with previous HIPASS
observations and with previous HI images of the Large Magellanic Cloud. Drift
scan observations of the GAMA G23 field resulted in direct HI detections at
$z=0.0043$ and $z=0.0055$ of HIPASS galaxies J224230 and J230930. Our new
measurements generally agree with archival data in spectral shape and flux
density, with small differences being due to differing beam patterns. We also
detect signal in the stacked HI data of 1094 individually undetected galaxies
in the GAMA G23 field in the redshift range $0.05 \leq z \leq 0.075$. Finally,
we use the low standing wave ripple and wide bandwidth of the PAF to set a
$3\sigma$ upper limit to any positronium recombination line emission from the
Galactic Centre of $<0.09$ K, corresponding to a recombination rate of
$<3.0\times10^{45}\,\mathrm{s}^{1}$.

We present a search for optical bursts from the repeating fast radio burst
FRB 121102 using simultaneous observations with the highspeed optical camera
ULTRASPEC on the 2.4m Thai National Telescope and radio observations with the
100m Effelsberg Radio Telescope. A total of 13 radio bursts were detected, but
we found no evidence for corresponding optical bursts in our 70.7ms frames.
The 5sigma upper limit to the optical flux density during our observations is
0.33 mJy at 767nm. This gives an upper limit for the optical burst fluence of
0.046 Jy ms, which constrains the broadband spectral index of the burst
emission to alpha < 0.2. Two of the radio pulses are separated by just 34 ms,
which may represent an upper limit on a possible underlying periodicity (a
rotation period typical of pulsars), or these pulses may have come from a
single emission window that is a small fraction of a possible period.

We report on the results of a 4year timing campaign of PSR~J2222$0137$, a
2.44day binary pulsar with a massive white dwarf (WD) companion, with the
Nan\c{c}ay, Effelsberg and Lovell radio telescopes. Using the Shapiro delay for
this system, we find a pulsar mass $m_{p}=1.76,\pm\,0.06,M_\odot$ and a WD mass
$m_{c}\,=\,1.293\,\pm\,0.025\, M_\odot$. We also measure the rate of advance of
periastron for this system, which is marginally consistent with the GR
prediction for these masses. The short lifetime of the massive WD progenitor
star led to a rapid Xray binary phase with little ($< \, 10^{2} \, M_\odot$)
mass accretion onto the neutron star (NS); hence, the current pulsar mass is,
within uncertainties, its birth mass; the largest measured to date. We discuss
the discrepancy with previous mass measurements for this system; we conclude
that the measurements presented here are likely to be more accurate. Finally,
we highlight the usefulness of this system for testing alternative theories of
gravity by tightly constraining the presence of dipolar radiation. This is of
particular importance for certain aspects of strongfield gravity, like
spontaneous scalarization, since the mass of PSR~J2222$0137$ puts that system
into a poorly tested parameter range.

We present the black hole accretion code (BHAC), a new multidimensional
generalrelativistic magnetohydrodynamics module for the MPIAMRVAC framework.
BHAC has been designed to solve the equations of ideal generalrelativistic
magnetohydrodynamics in arbitrary spacetimes and exploits adaptive mesh
refinement techniques with an efficient blockbased approach. Several
spacetimes have already been implemented and tested. We demonstrate the
validity of BHAC by means of various one, two, and threedimensional test
problems, as well as through a close comparison with the HARM3D code in the
case of a torus accreting onto a black hole. The convergence of a turbulent
accretion scenario is investigated with several diagnostics and we find
accretion rates and horizonpenetrating fluxes to be convergent to within a few
percent when the problem is run in three dimensions. Our analysis also involves
the study of the corresponding thermal synchrotron emission, which is performed
by means of a new generalrelativistic radiative transfer code, BHOSS. The
resulting synthetic intensity maps of accretion onto black holes are found to
be convergent with increasing resolution and are anticipated to play a crucial
role in the interpretation of horizonscale images resulting from upcoming
radio observations of the source at the Galactic Center.

We measure the effects of interstellar scattering on average pulse profiles
from 13 radio pulsars with simple pulse shapes. We use data from the LOFAR High
Band Antennas, at frequencies between 110 and 190~MHz. We apply a forward
fitting technique, and simultaneously determine the intrinsic pulse shape,
assuming single Gaussian component profiles. We find that the constant $\tau$,
associated with scattering by a single thin screen, has a powerlaw dependence
on frequency $\tau \propto \nu^{\alpha}$, with indices ranging from $\alpha =
1.50$ to $4.0$, despite simplest theoretical models predicting $\alpha = 4.0$
or $4.4$. Modelling the screen as an isotropic or extremely anisotropic
scatterer, we find anisotropic scattering fits lead to larger powerlaw
indices, often in better agreement with theoretically expected values. We
compare the scattering models based on the inferred, frequency dependent
parameters of the intrinsic pulse, and the resulting correction to the
dispersion measure (DM). We highlight the cases in which fits of extreme
anisotropic scattering are appealing, while stressing that the data do not
strictly favour either model for any of the 13 pulsars. The pulsars show
anomalous scattering properties that are consistent with finite scattering
screens and/or anisotropy, but these data alone do not provide the means for an
unambiguous characterization of the screens. We revisit the empirical $\tau$
versus DM relation and consider how our results support a frequency dependence
of $\alpha$. Very long baseline interferometry, and observations of the
scattering and scintillation properties of these sources at higher frequencies,
will provide further evidence.

Pulsars are natural cosmic clocks. On long timescales they rival the
precision of terrestrial atomic clocks. Using a technique called pulsar timing,
the exact measurement of pulse arrival times allows a number of applications,
ranging from testing theories of gravity to detecting gravitational waves. Also
an external reference system suitable for autonomous space navigation can be
defined by pulsars, using them as natural navigation beacons, not unlike the
use of GPS satellites for navigation on Earth. By comparing pulse arrival times
measured onboard a spacecraft with predicted pulse arrivals at a reference
location (e.g. the solar system barycenter), the spacecraft position can be
determined autonomously and with high accuracy everywhere in the solar system
and beyond. We describe the unique properties of pulsars that suggest that such
a navigation system will certainly have its application in future astronautics.
We also describe the ongoing experiments to use the clocklike nature of
pulsars to "construct" a galacticsized gravitational wave detector for
lowfrequency (f_GW ~1E9  1E7 Hz) gravitational waves. We present the
current status and provide an outlook for the future.

We evaluate dark matter (DM) limits from cosmicray antiproton observations
using the recent precise AMS02 measurements. We properly take into account
cosmicray propagation uncertainties fitting at the same time DM and
propagation parameters, and marginalizing over the latter. We find a
significant (~4.5 sigma) indication of a DM signal for DM masses near 80 GeV,
with a hadronic annihilation crosssection close to the thermal value, sigma v
~3e26 cm3s1. Intriguingly, this signal is compatible with the DM
interpretation of the Galactic center gammaray excess. Confirmation of the
signal will require a more accurate study of the systematic uncertainties,
i.e., the antiproton production crosssection, and modelling of the solar
modulation effect. Interpreting the AMS02 data in terms of upper limits on
hadronic DM annihilation, we obtain strong constraints excluding a thermal
annihilation crosssection for DM masses below about 50 GeV and in the range
between approximately 150 and 500 GeV, even for conservative propagation
scenarios. Except for the range around 80 GeV, our limits are a factor 4
stronger than the limits from gammaray observations of dwarf galaxies.

The annihilation of dark matter particles in the galactic halo of the Milky
Way may lead to cosmic ray signatures that can be probed by the AMS02
experiment, which has measured the composition and fluxes of charged cosmic
rays with unprecedented precision. Given the absence of characteristic spectral
features in the electron and positron fluxes measured by AMS02, we derive
upper limits on the dark matter annihilation cross section for leptophilic dark
matter models. Our limits are based on a new background model that describes
all recent measurements of the energy spectra of cosmic ray positrons and
electrons. For thermal dark matter relics, we can exclude dark matter masses
below about 100 GeV. We include the radiation of electroweak gauge bosons in
the dark matter annihilation process and compute the antiproton signal that can
be expected within leptophilic dark matter models.

A possible hint of dark matter annihilation has been found in Cuoco,
Korsmeier and Kr\"amer (2017) from an analysis of recent cosmicray antiproton
data from AMS02 and taking into account cosmicray propagation uncertainties
by fitting at the same time dark matter and propagation parameters. Here, we
extend this analysis to a wider class of annihilation channels. We find
consistent hints of a dark matter signal with an annihilation crosssection
close to the thermal value and with masses in range between 40 and 130 GeV
depending on the annihilation channel. Furthermore, we investigate in how far
the possible signal is compatible with the Galactic center gammaray excess and
recent observation of dwarf satellite galaxies by performing a joint global fit
including uncertainties in the dark matter density profile. As an example, we
interpret our results in the framework of the Higgs portal model.

SCYNet (SUSY Calculating Yield Net) is a tool for testing supersymmetric
models against LHC data. It uses neural network regression for a fast
evaluation of the profile likelihood ratio. Two neural network approaches have
been developed: one network has been trained using the parameters of the
11dimensional phenomenological Minimal Supersymmetric Standard Model
(pMSSM11) as an input and evaluates the corresponding profile likelihood ratio
within milliseconds. It can thus be used in global pMSSM11 fits without time
penalty. In the second approach, the neural network has been trained using
modelindependent signaturerelated objects, such as energies and particle
multiplicities, which were estimated from the parameters of a given new physics
model. While the calculation of the energies and particle multiplicities takes
up computation time, the corresponding neural network is more general and can
be used to predict the LHC profile likelihood ratio for a wider class of new
physics models.

The elucidation of the nature of new phenomena requires a multipronged
approach to understand the essential physics that underlies it. As an example,
we study the simplified model containing a new scalar singlet accompanied by
vectorlike quarks, as motivated by the recent diphoton excess at the LHC. To
be specific, we investigate three models with $SU(2)_L$doublet, vectorlike
quarks with Yukawa couplings to a new scalar singlet and which also couple
offdiagonally to corresponding Standard Model fermions of the first or third
generation through the usual Higgs boson. We demonstrate that three classes of
searches can play important and complementary roles in constraining this model.
In particular, we find that missing energy searches designed for superparticle
production, supply superior sensitivity for vectorlike quarks than the
dedicated new quark searches themselves.

We present the results of a closure phase analysis of 3 mm very long baseline
interferometry (VLBI) measurements performed on Sagittarius A* (Sgr A*). We
have analyzed observations made in May 2015 using the Very Long Baseline Array,
the Robert C. Byrd Green Bank Telescope and the Large Millimeter Telescope
Alfonso Serrano and obtained nonzero closure phase measurements on several
station triangles  indicative of a nonpointsymmetric source structure. The
data are fitted with an asymmetric source structure model in Sgr A*,
represented by a simple twocomponent model, which favours a fainter component
due East of the main source. This result is discussed in light of a scattering
screen with substructure or an intrinsically asymmetric source.

We present stateofthe art predictions for the production of supersymmetric
squarks and gluinos at the Large Hadron Collider (LHC), including softgluon
resummation up to nexttonexttoleading logarithmic (NNLL) accuracy, the
resummation of Coulomb corrections and the contribution from bound states. The
NNLL corrections enhance the crosssection predictions and reduce the scale
uncertainty to a level of 510%. The NNLL resummed crosssection predictions
can be obtained from the computer code NNLLfast, which also provides the scale
uncertainty and the pdf and alpha_s error.

Electroweak corrections can play a crucial role in dark matter annihilation.
The emission of gauge bosons, in particular, leads to a secondary flux
consisting of all Standard Model particles, and may be described by electroweak
fragmentation functions. To assess the quality of the fragmentation function
approximation to electroweak radiation in dark matter annihilation, we have
calculated the flux of secondary particles from gaugeboson emission in models
with Majorana fermion and vector dark matter, respectively. For both models, we
have compared cross sections and energy spectra of positrons and antiprotons
after propagation through the galactic halo in the fragmentation function
approximation and in the full calculation. Fragmentation functions fail to
describe the particle fluxes in the case of Majorana fermion annihilation into
light fermions: the helicity suppression of the lowestorder cross section in
such models cannot be lifted by the leading logarithmic contributions included
in the fragmentation function approach. However, for other classes of models
like vector dark matter, where the lowestorder cross section is not
suppressed, electroweak fragmentation functions provide a simple,
modelindependent and accurate description of secondary particle fluxes.

We analyse the excess in the $\gamma$ray emission from the center of our
galaxy observed by FermiLAT in terms of dark matter annihilation within the
scalar Higgs portal model. In particular, we include the astrophysical
uncertainties from the dark matter distribution and allow for unspecified
additional dark matter components. We demonstrate through a detailed numerical
fit that the strength and shape of the $\gamma$ray spectrum can indeed be
described by the model in various regions of dark matter masses and couplings.
Constraints from invisible Higgs decays, direct dark matter searches, indirect
searches in dwarf galaxies and for $\gamma$ray lines, and constraints from the
dark matter relic density reduce the parameter space to dark matter masses near
the Higgs resonance. We find two viable regions: one where the Higgsdark
matter coupling is of ${\cal O}(10^{2})$, and an additional dark matter
component beyond the scalar WIMP of our model is preferred, and one region
where the Higgsdark matter coupling may be significantly smaller, but where
the scalar WIMP constitutes a significant fraction or even all of dark matter.
Both viable regions are hard to probe in future direct detection and collider
experiments.