Indirect detection of dark matter via its annihilation products is a key technique in the search for dark matter in the form of weakly interacting massive particles (WIMPs). Strong constraints exist on the annihilation of WIMPs to highly visible Standard Model final states such as photons or charged particles. In the case of s-wave annihilation, this typically eliminates thermal relic cross sections for dark matter of mass below $\cal{O}$(10) GeV. However, such limits typically neglect the possibility that dark matter may annihilate to assumed invisible or hard-to-detect final states, such as neutrinos. This is a difficult paradigm to probe due to the weak neutrino interaction cross section. Considering dark matter annihilation in the Galactic halo, we study the prospects for indirect detection using the Hyper-Kamiokande (HyperK) neutrino experiment, for dark matter of mass below 1 GeV. We undertake a dedicated simulation of the HyperK detector, which we benchmark against results from the similar Super-Kamiokande experiment and HyperK physics projections. We provide projections for the annihilation cross-sections that can be probed by HyperK for annihilation to muon or neutrino final states, and discuss uncertainties associated with the dark matter halo profile. For neutrino final states, we find that HyperK is sensitive to thermal annihilation cross-sections for dark matter with mass around 20 MeV, assuming an NFW halo profile. We also discuss the effects of Gadolinium enhancement, and prospects for improving the reach at low mass.

Despite considerable experimental progress large parts of the axion-like particle (ALP) parameter space remain difficult to probe in terrestrial experiments. In some cases, however, small-scale structure of the ALP dark matter (DM) distribution is strongly enhanced, offering opportunities for astrophysical tests. Such an enhancement can be produced by a period of pre-nucleosynthesis early matter domination (EMD). This cosmology arises in many ultraviolet completions and generates the correct relic abundance for weak coupling $f_a\sim 10^{16}$ GeV, ALP masses in the range $10^{-13}$ eV $<m_a < 1$ eV, and without fine-tuning of the initial misalignment angle. This range includes the QCD axion around $10^{-9}-10^{-8}$ eV. EMD enhances the growth of ALP small-scale structure, leading to the formation of dense ALP miniclusters. We study the interplay between the initial ALP oscillation, reheating temperature, and effective pressure to provide analytic estimates of the minicluster abundance and properties. ALP miniclusters in the EMD cosmology are denser and more abundant than in $\Lambda$CDM. While enhanced substructure generically reduces the prospects of direct detection experiments, we show that pulsar timing and lensing observations can discover these minihalos over a large range of ALP masses and reheating temperatures.

Many existing and proposed experiments targeting QCD axion dark matter (DM) can also search for a broad class of axion-like particles (ALPs). We analyze the experimental sensitivities to electromagnetically-coupled ALP DM in different cosmological scenarios with the relic abundance set by the misalignment mechanism. We obtain benchmark DM targets for the standard thermal cosmology, a pre-nucleosynthesis period of early matter domination, and a period of kination. These targets are theoretically simple and assume $\mathcal{O}(1)$ misalignment angles, avoiding fine-tuning of the initial conditions. We find that some experiments will have sensitivity to these ALP DM targets before they are sensitive to the QCD axion, and others can potentially reach interesting targets below the QCD band. The ALP DM abundance also depends on the origin of the ALP mass. Temperature-dependent masses that are generated by strong dynamics (as for the QCD axion) correspond to DM candidates with smaller decay constants, resulting in even better detection prospects.

We report the results of a global analysis of dark matter simplified models (DMSMs) with leptophobic mediator particles of spin one, considering the cases of both vector and axial-vector interactions with dark matter (DM) particles and quarks. We require the DMSMs to provide all the cosmological DM density indicated by Planck and other observations, and we impose the upper limits on spin-independent and -dependent scattering from direct DM search experiments. We also impose all relevant LHC constraints from searches for monojet events and measurements of the dijet mass spectrum. We model the likelihood functions for all the constraints and combine them within the MasterCode framework, and probe the full DMSM parameter spaces by scanning over the mediator and DM masses and couplings, not fixing any of the model parameters. We find, in general, two allowed regions of the parameter spaces: one in which the mediator couplings to Standard Model (SM) and DM particles may be comparable to those in the SM and the cosmological DM density is reached via resonant annihilation, and one in which the mediator couplings to quarks are $\lesssim 10^{-3}$ and DM annihilation is non-resonant. We find that the DM and mediator masses may well lie within the ranges accessible to LHC experiments. We also present predictions for spin-independent and -dependent DM scattering, and present specific results for ranges of the DM couplings that may be favoured in ultraviolet completions of the DMSMs.

Dark matter (DM) particles with mass in the sub-GeV range are an attractive alternative to heavier weakly-interacting massive particles, but direct detection of such light particles is challenging. If however DM-nucleus scattering leads to ionisation of the recoiling atom, the resulting electron may be detected even if the nuclear recoil is unobservable. We demonstrate that including this effect significantly enhances direct detection sensitivity to sub-GeV DM. Existing experiments set world-leading limits, and future experiments may probe the cross sections relevant for thermal freeze-out.

We use MasterCode to perform a frequentist analysis of the constraints on a phenomenological MSSM model with 11 parameters, the pMSSM11, including constraints from ~ 36/fb of LHC data at 13 TeV and PICO, XENON1T and PandaX-II searches for dark matter scattering, as well as previous accelerator and astrophysical measurements, presenting fits both with and without the $(g-2)_{\mu}$ constraint. The pMSSM11 is specified by the following parameters: 3 gaugino masses $M_{1,2,3}$, a common mass for the first-and second-generation squarks $m_{\tilde{q}}$ and a distinct third-generation squark mass $m_{\tilde{q}_3}$, a common mass for the first-and second-generation sleptons $m_{\tilde l}$ and a distinct third-generation slepton mass $m_{\tilde \tau}$, a common trilinear mixing parameter $A$, the Higgs mixing parameter $\mu$, the pseudoscalar Higgs mass $M_A$ and $\tan\beta$. In the fit including $(g-2)_{\mu}$, a Bino-like $\tilde\chi^0_1$ is preferred, whereas a Higgsino-like $\tilde \chi^0_1$ is favoured when the $(g-2)_{\mu}$ constraint is dropped. We identify the mechanisms that operate in different regions of the pMSSM11 parameter space to bring the relic density of the lightest neutralino, $\tilde\chi^0_1$, into the range indicated by cosmological data. In the fit including $(g-2)_{\mu}$, coannihilations with $\tilde \chi^0_2$ and the Wino-like $\tilde\chi^{\pm}_1$ or with nearly-degenerate first- and second-generation sleptons are favoured, whereas coannihilations with the $\tilde \chi^0_2$ and the Higgsino-like $\tilde\chi^{\pm}_1$ or with first- and second-generation squarks may be important when the $(g-2)_{\mu}$ constraint is dropped. Prospects remain for discovering strongly-interacting sparticles at the LHC as well as for discovering electroweakly-interacting sparticles at a future linear $e^+ e^-$ collider such as the ILC or CLIC.

We describe a likelihood analysis using MasterCode of variants of the MSSM in which the soft supersymmetry-breaking parameters are assumed to have universal values at some scale $M_{in}$ below the supersymmetric grand unification scale $M_{GUT}$, as can occur in mirage mediation and other models. In addition to $M_{in}$, such `sub-GUT' models have the 4 parameters of the CMSSM, namely a common gaugino mass $m_{1/2}$, a common soft supersymmetry-breaking scalar mass $m_0$, a common trilinear mixing parameter $A$ and the ratio of MSSM Higgs vevs $\tan\beta$, assuming that the Higgs mixing parameter $\mu > 0$. We take into account constraints on strongly- and electroweakly-interacting sparticles from $\sim 36$/fb of LHC data at 13 TeV and the LUX and 2017 PICO, XENON1T and PandaX-II searches for dark matter scattering, in addition to the previous LHC and dark matter constraints as well as full sets of flavour and electroweak constraints. We find a preference for $M_{in} \sim 10^5$ to $10^9$ GeV, with $M_{in} \sim M_{GUT}$ disfavoured by $\Delta \chi^2 \sim 3$ due to the ${\rm BR}(B_{s, d} \to \mu^+\mu^-)$ constraint. The lower limits on strongly-interacting sparticles are largely determined by LHC searches, and similar to those in the CMSSM. We find a preference for the LSP to be a Bino or Higgsino with $\tilde{\chi^0_1} \sim 1$ TeV, with annihilation via heavy Higgs bosons $H/A$ and stop coannihilation, or chargino coannihilation, bringing the cold dark matter density into the cosmological range. We find that spin-independent dark matter scattering is likely to be within reach of the planned LUX-Zeplin and XENONnT experiments. We probe the impact of the $(g-2)_\mu$ constraint, finding similar results whether or not it is included.

We perform a likelihood analysis of the constraints from accelerator experiments and astrophysical observations on supersymmetric (SUSY) models with SU(5) boundary conditions on soft SUSY-breaking parameters at the GUT scale. The parameter space of the models studied has 7 parameters: a universal gaugino mass $m_{1/2}$, distinct masses for the scalar partners of matter fermions in five- and ten-dimensional representations of SU(5), $m_5$ and $m_{10}$, and for the $\mathbf{5}$ and $\mathbf{\bar 5}$ Higgs representations $m_{H_u}$ and $m_{H_d}$, a universal trilinear soft SUSY-breaking parameter $A_0$, and the ratio of Higgs vevs $\tan \beta$. In addition to previous constraints from direct sparticle searches, low-energy and flavour observables, we incorporate constraints based on preliminary results from 13 TeV LHC searches for jets + MET events and long-lived particles, as well as the latest PandaX-II and LUX searches for direct Dark Matter detection. In addition to previously-identified mechanisms for bringing the supersymmetric relic density into the range allowed by cosmology, we identify a novel ${\tilde u_R}/{\tilde c_R} - \tilde{\chi}^0_1$ coannihilation mechanism that appears in the supersymmetric SU(5) GUT model and discuss the role of ${\tilde \nu_\tau}$ coannihilation. We find complementarity between the prospects for direct Dark Matter detection and SUSY searches at the LHC.

We perform a likelihood analysis of the minimal Anomaly-Mediated Supersymmetry Breaking (mAMSB) model using constraints from cosmology and accelerator experiments. We find that a wino-like or a Higgsino-like neutralino LSP, $m_{\tilde \chi^0_{1}}$, may provide the cold dark matter (DM) with similar likelihood. The upper limit on the DM density from Planck and other experiments enforces $m_{\tilde \chi^0_{1}} \lesssim 3~TeV$ after the inclusion of Sommerfeld enhancement in its annihilations. If most of the cold DM density is provided by the $\tilde \chi_0^1$, the measured value of the Higgs mass favours a limited range of $\tan \beta \sim 5$ (or for $\mu > 0$, $\tan \beta \sim 45$) but the scalar mass $m_0$ is poorly constrained. In the wino-LSP case, $m_{3/2}$ is constrained to about $900~TeV$ and ${m_{\tilde \chi^0_{1}}}$ to $2.9\pm0.1~TeV$, whereas in the Higgsino-LSP case $m_{3/2}$ has just a lower limit $\gtrsim 650TeV$ ($\gtrsim 480TeV$) and $m_{\tilde \chi^0_{1}}$ is constrained to $1.12 ~(1.13) \pm0.02~TeV$ in the $\mu>0$ ($\mu<0$) scenario. In neither case can the anomalous magnetic moment of the muon, ${(g-2)_\mu}$, be improved significantly relative to its Standard Model (SM) value, nor do flavour measurements constrain the model significantly, and there are poor prospects for discovering supersymmetric particles at the LHC, {though there} are some prospects for direct DM detection. On the other hand, if the ${m_{\tilde \chi^0_{1}}}$ contributes only a fraction of the cold DM density, {future LHC $E_T$-based searches for gluinos, squarks and heavier chargino and neutralino states as well as disappearing track searches in the wino-like LSP region will be relevant}, and interference effects enable ${\rm BR}(B_{s, d} \to \mu^+\mu^-)$ to agree with the data better than in the SM in the case of wino-like DM with $\mu > 0$.

We have constructed a series of non-rotating quasi-hydrostatic evolutionary models for the M2 Iab supergiant Betelgeuse ($\alpha~Orionis$). Our models are constrained by multiple observed values for the temperature, luminosity, surface composition and mass loss for this star, along with the parallax distance and high resolution imagery that determines its radius. We have then applied our best-fit models to analyze the observed variations in surface luminosity and the size of detected surface bright spots as the result of up-flowing convective material from regions of high temperature in the surface convective zone. We also attempt to explain the intermittently observed periodic variability in a simple radial linear adiabatic pulsation model. Based upon the best fit to all observed data, we suggest a best progenitor mass estimate of $ 20 ^{+5}_{-3} M_\odot$ and a current age from the start of the zero-age main sequence of $8.0 - 8.5$ Myr based upon the observed ejected mass while on the giant branch.

Different mechanisms operate in various regions of the MSSM parameter space to bring the relic density of the lightest neutralino, neutralino_1, assumed here to be the LSP and thus the Dark Matter (DM) particle, into the range allowed by astrophysics and cosmology. These mechanisms include coannihilation with some nearly-degenerate next-to-lightest supersymmetric particle (NLSP) such as the lighter stau (stau_1), stop (stop_1) or chargino (chargino_1), resonant annihilation via direct-channel heavy Higgs bosons H/A, the light Higgs boson h or the Z boson, and enhanced annihilation via a larger Higgsino component of the LSP in the focus-point region. These mechanisms typically select lower-dimensional subspaces in MSSM scenarios such as the CMSSM, NUHM1, NUHM2 and pMSSM10. We analyze how future LHC and direct DM searches can complement each other in the exploration of the different DM mechanisms within these scenarios. We find that the stau_1 coannihilation regions of the CMSSM, NUHM1, NUHM2 can largely be explored at the LHC via searches for missing E_T events and long-lived charged particles, whereas their H/A funnel, focus-point and chargino_1 coannihilation regions can largely be explored by the LZ and Darwin DM direct detection experiments. We find that the dominant DM mechanism in our pMSSM10 analysis is chargino_1 coannihilation: {parts of its parameter space can be explored by the LHC, and a larger portion by future direct DM searches.

We present a frequentist analysis of the parameter space of the pMSSM10, in which the following 10 soft SUSY-breaking parameters are specified independently at the mean scalar top mass scale Msusy = Sqrt[M_stop1 M_stop2]: the gaugino masses M_{1,2,3}, the 1st-and 2nd-generation squark masses M_squ1 = M_squ2, the third-generation squark mass M_squ3, a common slepton mass M_slep and a common trilinear mixing parameter A, the Higgs mixing parameter mu, the pseudoscalar Higgs mass M_A and tan beta. We use the MultiNest sampling algorithm with 1.2 x 10^9 points to sample the pMSSM10 parameter space. A dedicated study shows that the sensitivities to strongly-interacting SUSY masses of ATLAS and CMS searches for jets, leptons + MET signals depend only weakly on many of the other pMSSM10 parameters. With the aid of the Atom and Scorpion codes, we also implement the LHC searches for EW-interacting sparticles and light stops, so as to confront the pMSSM10 parameter space with all relevant SUSY searches. In addition, our analysis includes Higgs mass and rate measurements using the HiggsSignals code, SUSY Higgs exclusion bounds, the measurements B-physics observables, EW precision observables, the CDM density and searches for spin-independent DM scattering. We show that the pMSSM10 is able to provide a SUSY interpretation of (g-2)_mu, unlike the CMSSM, NUHM1 and NUHM2. As a result, we find (omitting Higgs rates) that the minimum chi^2/dof = 20.5/18 in the pMSSM10, corresponding to a chi^2 probability of 30.8 %, to be compared with chi^2/dof = 32.8/24 (31.1/23) (30.3/22) in the CMSSM (NUHM1) (NUHM2). We display 1-dimensional likelihood functions for SUSY masses, and show that they may be significantly lighter in the pMSSM10 than in the CMSSM, NUHM1 and NUHM2. We discuss the discovery potential of future LHC runs, e+e- colliders and direct detection experiments.

Previous attempts at explaining the gamma-ray excess near the Galactic Centre have focussed on dark matter annihilation directly into Standard Model particles. This results in a preferred dark matter mass of 30-40 GeV (if the annihilation is into b quarks) or 10 GeV (if it is into leptons). Here we show that the gamma-ray excess is also consistent with heavier dark matter particles; in models of secluded dark matter, dark matter with mass up to 76 GeV provides a good fit to the data. This occurs if the dark matter first annihilates to an on-shell particle that subsequently decays to Standard Model particles through a portal interaction. This is a generic process that works in models with annihilation, semi-annihilation or both. We explicitly demonstrate this in a model of hidden vector dark matter with an SU(2) gauge group in the hidden sector.

We show that it is possible for WIMP dark matter to produce a large signal in indirect dark matter searches without producing signals elsewhere. We illustrate our point by fitting the Fermi-LAT extended galactic gamma-ray excess with a simple model of Dirac dark matter that annihilates primarily into b quarks via a pseudoscalar. Current collider constraints are weak while the 14 TeV LHC run will constrain a limited portion of the parameter space. No signal is expected in additional indirect searches or at future direct detection experiments. Our results emphasise the importance of fully understanding potential indirect signals of dark matter as they may provide the only information about the dark matter particle.

Motivated by recent claims of lines in the Fermi gamma-ray spectrum, we critically examine means of enhancing neutralino annihilation into neutral gauge bosons. The signal can be boosted while remaining consistent with continuum photon constraints if a new singlet-like pseudoscalar is present. We consider singlet extensions of the MSSM, focusing on the NMSSM, where a `well-tempered' neutralino can explain the lines while remaining consistent with current constraints. We adopt a complementary numerical and analytic approach throughout in order to gain intuition for the underlying physics. The scenario requires a rich spectrum of light neutralinos and charginos leading to characteristic phenomenological signatures at the LHC whose properties we explore. Future direct detection prospects are excellent, with sizeable spin-dependent and spin-independent cross-sections.

Recent work on increasing the effective number of neutrino species (Neff) in the early universe has focussed on introducing extra relativistic species (`dark radiation'). We draw attention to another possibility: a new particle of mass less than 10 MeV that remains in thermal equilibrium with neutrinos until it becomes non-relativistic increases the neutrino temperature relative to the photons. We demonstrate that this leads to a value of Neff that is greater than three and that Neff at CMB formation is larger than at BBN. We investigate the constraints on such particles from the primordial abundance of helium and deuterium created during BBN and from the CMB power spectrum measured by ACT and SPT and find that they are presently relatively unconstrained. We forecast the sensitivity of the Planck satellite to this scenario: in addition to dramatically improving constraints on the particle mass, in some regions of parameter space it can discriminate between the new particle being a real or complex scalar.

Global frequentist fits to the CMSSM and NUHM1 using the MasterCode framework predicted m_h \simeq 119 GeV in fits incorporating the g_mu-2 constraint and \simeq 126 GeV without it. Recent results by ATLAS and CMS could be compatible with a Standard Model-like Higgs boson around m_h \simeq 125 GeV. We use the previous MasterCode analysis to calculate the likelihood for a measurement of any nominal Higgs mass within the range of 115 to 130 GeV. Assuming a Higgs mass measurement at m_h \simeq 125 GeV, we display updated global likelihood contours in the (m_0, m_{1/2}) and other parameter planes of the CMSSM and NUHM1, and present updated likelihood functions for m_gluino, m_squark, B to mu mu, and the spin-independent dark matter cross section \sigma^si. The implications of dropping g_mu-2 from the fits are also discussed. We furthermore comment on a hypothetical measurement of m_h \simeq 119 GeV.

We make frequentist analyses of the CMSSM, NUHM1, VCMSSM and mSUGRA parameter spaces taking into account all the public results of searches for supersymmetry using data from the 2010 LHC run and the Xenon100 direct search for dark matter scattering. The LHC data set includes ATLAS and CMS searches for jets + ETslash events (with or without leptons) and for the heavier MSSM Higgs bosons, and the upper limit on bs to mu mu including data from LHCb as well as CDF and D0. The absences of signals in the LHC data favour somewhat heavier mass spectra than in our previous analyses of the CMSSM, NUHM1 and VCMSSM, and somewhat smaller dark matter scattering cross sections, all close to or within the pre-LHC 68% CL ranges, but do not impact significantly the favoured regions of the mSUGRA parameter space. We also discuss the impact of the Xenon100 constraint on spin-independent dark matter scattering, stressing the importance of taking into account the uncertainty in the pi-nucleon sigma term, that affects the spin-independent scattering matrix element, and we make predictions for spin-dependent dark matter scattering. Finally, we discuss briefly the potential impact of the updated predictions for sparticle masses in the CMSSM, NUHM1, VCMSSM and mSUGRA on future e+ e- colliders.

We make a frequentist analysis of the parameter space of minimal supergravity (mSUGRA), in which, as well as the gaugino and scalar soft supersymmetry-breaking parameters being universal, there is a specific relation between the trilinear, bilinear and scalar supersymmetry-breaking parameters, A_0 = B_0 + m_0, and the gravitino mass is fixed by m_{3/2} = m_0. We also consider a more general model, in which the gravitino mass constraint is relaxed (the VCMSSM). We combine in the global likelihood function the experimental constraints from low-energy electroweak precision data, the anomalous magnetic moment of the muon, the lightest Higgs boson mass M_h, B physics and the astrophysical cold dark matter density, assuming that the lightest supersymmetric particle (LSP) is a neutralino. In the VCMSSM, we find a preference for values of m_{1/2} and m_0 similar to those found previously in frequentist analyses of the constrained MSSM (CMSSM) and a model with common non-universal Higgs masses (NUHM1). On the other hand, in mSUGRA we find two preferred regions: one with larger values of both m_{1/2} and m_0 than in the VCMSSM, and one with large m_0 but small m_{1/2}. We compare the probabilities of the frequentist fits in mSUGRA, the VCMSSM, the CMSSM and the NUHM1: the probability that mSUGRA is consistent with the present data is significantly less than in the other models. We also discuss the mSUGRA and VCMSSM predictions for sparticle masses and other observables, identifying potential signatures at the LHC and elsewhere.

We study the problem of selecting between different mechanisms of supersymmetry breaking in the Minimal Supersymmetric Standard Model using current data. We evaluate the Bayesian evidence of four supersymmetry breaking scenarios: mSUGRA, mGMSB, mAMSB and moduli mediation. The results show a strong dependence on the dark matter assumption. Using the inferred cosmological relic density as an upper bound, minimal anomaly mediation is at least moderately favoured over the CMSSM. Our fits also indicate that evidence for a positive sign of the mu parameter is moderate at best. We present constraints on the anomaly and gauge mediated parameter spaces and some previously unexplored aspects of the dark matter phenomenology of the moduli mediation scenario. We use sparticle searches, indirect observables and dark matter observables in the global fit and quantify robustness with respect to prior choice. We quantify how much information is contained within each constraint.