The excess of electron recoil events seen by the XENON1T experiment has been interpreted as a potential signal of axion-like particles (ALPs), either produced in the Sun, or constituting part of the dark matter halo of the Milky Way. It has also been explained as a consequence of trace amounts of tritium in the experiment. We consider the evidence for the solar and dark-matter ALP hypotheses from the combination of XENON1T data and multiple astrophysical probes, including horizontal branch stars, red giants, and white dwarfs. We briefly address the influence of ALP decays and supernova cooling. While the different datasets are in clear tension for the case of solar ALPs, all measurements can be simultaneously accommodated for the case of a sub-dominant fraction of dark-matter ALPs. Nevertheless, this solution requires the tuning of several a priori unknown parameters, such that for our choices of priors a Bayesian analysis shows no strong preference for the ALP interpretation of the XENON1T excess over the background hypothesis.

Nested sampling (NS) is an invaluable tool in data analysis in modern astrophysics, cosmology, gravitational wave astronomy and particle physics. We identify a previously unused property of NS related to order statistics: the insertion indexes of new live points into the existing live points should be uniformly distributed. This observation enabled us to create a novel cross-check of single NS runs. The tests can detect when an NS run failed to sample new live points from the constrained prior and plateaus in the likelihood function, which break an assumption of NS and thus leads to unreliable results. We applied our cross-check to NS runs on toy functions with known analytic results in 2 - 50 dimensions, showing that our approach can detect problematic runs on a variety of likelihoods, settings and dimensions. As an example of a realistic application, we cross-checked NS runs performed in the context of cosmological model selection. Since the cross-check is simple, we recommend that it become a mandatory test for every applicable NS run.

We present BubbleProfiler, a C++ software package for finding field profiles in bubble walls and calculating the bounce action during phase transitions involving multiple scalar fields. Our code uses a recently proposed perturbative method for potentials with multiple fields and a shooting method for single field cases. BubbleProfiler is constructed with modularity, flexibility and practicality in mind. These principles extend from the input of an arbitrary potential with multiple scalar fields in various forms, through the code structure, to the testing suite. After reviewing the physics context, we describe how the methods are implemented in BubbleProfiler, provide an overview of the code structure and detail usage scenarios. We present a number of examples that serve as test cases of BubbleProfiler and comparisons to existing public codes with similar functionality. We also show a physics application of BubbleProfiler in the scalar singlet extension of the Standard Model of particle physics by calculating the action as a function of model parameters during the electroweak phase transition. BubbleProfiler completes an important link in the toolchain for studying the properties of the thermal phase transition driving baryogenesis and properties of gravitational waves in models with multiple scalar fields. The code can be obtained from: https://github.com/bubbleprofiler/bubbleprofiler

Searches for supersymmetric electroweakinos have entered a crucial phase, as the integrated luminosity of the Large Hadron Collider is now high enough to compensate for their weak production cross-sections. Working in a framework where the neutralinos and charginos are the only light sparticles in the Minimal Supersymmetric Standard Model, we use gambit to perform a detailed likelihood analysis of the electroweakino sector. We focus on the impacts of recent ATLAS and CMS searches with 36 fb$^{-1}$ of 13 TeV proton-proton collision data. We also include constraints from LEP and invisible decays of the $Z$ and Higgs bosons. Under the background-only hypothesis, we show that current LHC searches do not robustly exclude any range of neutralino or chargino masses. However, a pattern of excesses in several LHC analyses points towards a possible signal, with neutralino masses of $(m_{\tilde{\chi}_1^0}, m_{\tilde{\chi}_2^0}, m_{\tilde{\chi}_3^0}, m_{\tilde{\chi}_4^0})$ = (8-155, 103-260, 130-473, 219-502) GeV and chargino masses of $(m_{\tilde{\chi}_1^{\pm}}, m_{\tilde{\chi}_2^{\pm}})$ = (104-259, 224-507) GeV at the 95% confidence level. The lightest neutralino is mostly bino, with a possible modest Higgsino or wino component. We find that this excess has a combined local significance of $3.3\sigma$, subject to a number of cautions. If one includes LHC searches for charginos and neutralinos conducted with 8 TeV proton-proton collision data, the local significance is lowered to 2.9$\sigma$. We briefly consider the implications for dark matter, finding that the correct relic density can be obtained through the Higgs-funnel and $Z$-funnel mechanisms, even assuming that all other sparticles are decoupled. All samples, gambit input files and best-fit models from this study are available on Zenodo.

We present global analyses of effective Higgs portal dark matter models in the frequentist and Bayesian statistical frameworks. Complementing earlier studies of the scalar Higgs portal, we use GAMBIT to determine the preferred mass and coupling ranges for models with vector, Majorana and Dirac fermion dark matter. We also assess the relative plausibility of all four models using Bayesian model comparison. Our analysis includes up-to-date likelihood functions for the dark matter relic density, invisible Higgs decays, and direct and indirect searches for weakly-interacting dark matter including the latest XENON1T data. We also account for important uncertainties arising from the local density and velocity distribution of dark matter, nuclear matrix elements relevant to direct detection, and Standard Model masses and couplings. In all Higgs portal models, we find parameter regions that can explain all of dark matter and give a good fit to all data. The case of vector dark matter requires the most tuning and is therefore slightly disfavoured from a Bayesian point of view. In the case of fermionic dark matter, we find a strong preference for including a CP-violating phase that allows suppression of constraints from direct detection experiments, with odds in favour of CP violation of the order of 100:1. Finally, we present DDCalc 2.0.0, a tool for calculating direct detection observables and likelihoods for arbitrary non-relativistic effective operators.

We investigate the impact of uncertainty in the velocity distribution of dark matter on direct detection experiments. We construct an multinomial prior with a hyperparameter $\beta$ that describes the strength of our belief in an isotropic Maxwell-Boltzmann velocity distribution. By varying $\beta$, we interpolate between a halo-independent and halo-dependent analysis. We present a novel approximation for the marginalisation of this prior that is applicable to any counting experiment. With this formula, we investigate the impact of the uncertainty in limits from XENON1T. For dark matter masses greater than about 60 GeV, we find extremely mild sensitivity to the distribution. Below about 60 GeV, the limit weakens by less than an order of magnitude if we assume an isotropic distribution in the galactic frame. If we permit anisotropic distributions, the limit further weakens, but at most by about two orders of magnitude. Lastly, we check the impact of parametric uncertainties and discuss the possible inclusion and impact of our technique in global fits.

The Dark Matter Particle Explorer (DAMPE) recently released measurements of the electron spectrum with a hint of a narrow peak at about 1.4 TeV. We investigate dark matter (DM) models that could produce such a signal by annihilation in a nearby subhalo whilst simultaneously satisfying constraints from DM searches. In our model-independent approach, we consider all renormalizable interactions via a spin 0 or 1 mediator between spin 0 or 1/2 DM particles and the Standard Model leptons. We find that of the 20 combinations, 10 are ruled out by velocity or helicity suppression of the annihilation cross section to fermions. The remaining 10 models, though, evade constraints from the relic density, collider and direct detection searches, and include models of spin 0 and 1/2 DM coupling to a spin 0 or 1 mediator. We delineate the regions of mediator mass and couplings that could explain the DAMPE excess. In all cases the mediator is required to be heaver than about 2 TeV by LEP limits.

Using the DAMA/LIBRA anomaly as an example, we formalise the notion of halo-independence in the context of Bayesian statistics and quantified maximum entropy. We consider an infinite set of possible profiles, weighted by an entropic prior and constrained by a likelihood describing noisy measurements of modulated moments by DAMA/LIBRA. Assuming an isotropic dark matter (DM) profile in the galactic rest frame, we find the most plausible DM profiles and predictions for unmodulated signal rates at DAMA/LIBRA. The entropic prior contains an a priori unknown regularisation factor, $\beta$, that describes the strength of our conviction that the profile is approximately Maxwellian. By varying $\beta$, we smoothly interpolate between a halo-independent and a halo-dependent analysis, thus exploring the impact of prior information about the DM profile.

We present an application, Superplot, for calculating and plotting statistical quantities relevant to parameter inference from a "chain" of samples drawn from a parameter space, produced by e.g. MultiNest. A simple graphical interface allows one to browse a chain of many variables quickly, and make publication quality plots of, inter alia, one- and two-dimensional profile likelihood, posterior pdf (with kernel density estimation), confidence intervals and credible regions. In this short manual, we document installation and basic usage, and define all statistical quantities and conventions. The code is fully compatible with Linux, Windows and Mac OSX. Furthermore, if preferred, all functionality is available through the command line rather than a graphical interface.

We explore the MSSM with 9 free parameters (p9MSSM) that have been selected as a minimum set that allows an investigation of neutralino dark matter and collider signatures while maintaining consistency with several constraints. These include measurement of the dark matter relic density from PLANCK, main properties of the discovered Higgs boson, LHC direct SUSY searches, recent evidence for a Standard Model-like BR(Bs -> \mu+ \mu-), and the measurement of \delta(g-2), plus a number of other electroweak and flavor physics constraints. We perform a simulation of two LHC direct SUSY searches at sqrt(s)=8 TeV: the CMS inclusive \alpha_T search for squarks and gluinos and the CMS electroweak production search with 3l+E_T^miss in the final state. We use the latter to identify the regions of the parameter space, consistent at 2\sigma\ with \delta(g-2), that are not excluded by the direct limits from the electroweak production. We find that they correspond to a neutralino mass in the window 200 GeV<m_\chi<500 GeV. We also implement the likelihood for the XENON100 exclusion bound, in which we consider for the first time the impact of a recent determination of the \Sigma_{\pi N} term from CHAOS data, \Sigma_{\pi N}=43\pm12 MeV. We show that in light of this measurement, the present statistical impact of the XENON100 bound is greatly reduced, although future sensitivities of the LUX and XENON1T experiments will have decisive impact on the mixed bino/higgsino composition of the neutralino. We point out some tension between the constraints from \delta(g-2) and XENON100. Finally, we present prospects for various indirect searches of dark matter, namely \gamma-ray fluxes from dSphs and the Galactic Center at Fermi-LAT, and the positron flux at AMS02. We also show the 5-year sensitivity on the spin-dependent neutralino-proton cross section due to neutrino fluxes from the Sun at IceCube.