• ### SpecBit, DecayBit and PrecisionBit: GAMBIT modules for computing mass spectra, particle decay rates and precision observables(1705.07936)

Jan. 17, 2019 hep-ph
We present the GAMBIT modules SpecBit, DecayBit and PrecisionBit. Together they provide a new framework for linking publicly available spectrum generators, decay codes and other precision observable calculations in a physically and statistically consistent manner. This allows users to automatically run various combinations of existing codes as if they are a single package. The modular design allows software packages fulfilling the same role to be exchanged freely at runtime, with the results presented in a common format that can be easily passed to downstream dark matter, collider and flavour codes. These modules constitute an essential part of the broader GAMBIT framework, a major new software package for performing global fits. In this paper we present the observable calculations, data, and likelihood functions implemented in the three modules, as well as the conventions and assumptions used in interfacing them with external codes. We also present 3-BIT-HIT, a command-line utility for computing mass spectra, couplings, decays and precision observables in the MSSM, which shows how the three modules can be easily used independently of GAMBIT.
• ### Global fits of GUT-scale SUSY models with GAMBIT(1705.07935)

Oct. 23, 2018 hep-ph, astro-ph.CO
We present the most comprehensive global fits to date of three supersymmetric models motivated by grand unification: the Constrained Minimal Supersymmetric Standard Model (CMSSM), and its Non-Universal Higgs Mass generalisations NUHM1 and NUHM2. We include likelihoods from a number of direct and indirect dark matter searches, a large collection of electroweak precision and flavour observables, direct searches for supersymmetry at LEP and Runs I and II of the LHC, and constraints from Higgs observables. Our analysis improves on existing results not only in terms of the number of included observables, but also in the level of detail with which we treat them, our sampling techniques for scanning the parameter space, and our treatment of nuisance parameters. We show that stau co-annihilation is now ruled out in the CMSSM at more than 95\% confidence. Stop co-annihilation turns out to be one of the most promising mechanisms for achieving an appropriate relic density of dark matter in all three models, whilst avoiding all other constraints. We find high-likelihood regions of parameter space featuring light stops and charginos, making them potentially detectable in the near future at the LHC. We also show that tonne-scale direct detection will play a largely complementary role, probing large parts of the remaining viable parameter space, including essentially all models with multi-TeV neutralinos.
• ### A global fit of the MSSM with GAMBIT(1705.07917)

Oct. 22, 2018 hep-ph, astro-ph.CO
We study the seven-dimensional Minimal Supersymmetric Standard Model (MSSM7) with the new GAMBIT software framework, with all parameters defined at the weak scale. Our analysis significantly extends previous weak-scale, phenomenological MSSM fits, by adding more and newer experimental analyses, improving the accuracy and detail of theoretical predictions, including dominant uncertainties from the Standard Model, the Galactic dark matter halo and the quark content of the nucleon, and employing novel and highly-efficient statistical sampling methods to scan the parameter space. We find regions of the MSSM7 that exhibit co-annihilation of neutralinos with charginos, stops and sbottoms, as well as models that undergo resonant annihilation via both light and heavy Higgs funnels. We find high-likelihood models with light charginos, stops and sbottoms that have the potential to be within the future reach of the LHC. Large parts of our preferred parameter regions will also be accessible to the next generation of direct and indirect dark matter searches, making prospects for discovery in the near future rather good.
• ### Status of the scalar singlet dark matter model(1705.07931)

Oct. 19, 2018 hep-ph, astro-ph.CO
One of the simplest viable models for dark matter is an additional neutral scalar, stabilised by a $\mathbb{Z}_2$ symmetry. Using the GAMBIT package and combining results from four independent samplers, we present Bayesian and frequentist global fits of this model. We vary the singlet mass and coupling along with 13 nuisance parameters, including nuclear uncertainties relevant for direct detection, the local dark matter density, and selected quark masses and couplings. We include the dark matter relic density measured by Planck, direct searches with LUX, PandaX, SuperCDMS and XENON100, limits on invisible Higgs decays from the Large Hadron Collider, searches for high-energy neutrinos from dark matter annihilation in the Sun with IceCube, and searches for gamma rays from annihilation in dwarf galaxies with the Fermi-LAT. Viable solutions remain at couplings of order unity, for singlet masses between the Higgs mass and about 300 GeV, and at masses above $\sim$1 TeV. Only in the latter case can the scalar singlet constitute all of dark matter. Frequentist analysis shows that the low-mass resonance region, where the singlet is about half the mass of the Higgs, can also account for all of dark matter, and remains viable. However, Bayesian considerations show this region to be rather fine-tuned.
• ### Dark matter substructure cannot explain properties of the Fermi Galactic Centre excess(1612.01539)

Sept. 2, 2018 hep-ph, astro-ph.CO, astro-ph.HE
An excess of gamma rays has been identified at the centre of the Milky Way, and annihilation of dark matter has been posited as a potential source. This hypothesis faces significant challenges: difficulty characterizing astrophysical backgrounds, the need for a non-trivial adiabatic contraction of the inner part of the Milky Way's dark matter halo, and recent observations of photon clustering, which suggest that the majority of the excess is due to unresolved point sources. Here we point out that the apparent point-like nature of the emission rules out the dark matter interpretation of the excess entirely. Attempting to model the emission with dark matter point sources either worsens the problem with the inner slope, requires an unrealistically large minihalo fraction toward the Galactic Centre, or overproduces the observed emission at higher latitudes.
• ### Pitfalls of iterative pole mass calculation in electroweak multiplets(1710.01511)

April 15, 2018 hep-ph
The radiatively-induced mass splitting between components of an electroweak multiplet is typically of order 100 MeV. This is sufficient to endow the charged components with macroscopically-observable lifetimes, and ensure an electrically-neutral dark matter particle. We show that a commonly used iterative procedure to compute radiatively-corrected pole masses can lead to very different mass splittings than a non-iterative calculation at the same loop order. By estimating the uncertainties of the two one-loop results, we show that the iterative procedure is significantly more sensitive to the choice of renormalisation scale and gauge parameter than the non-iterative method. This can cause the lifetime of the charged component to vary by up to 12 orders of magnitude if iteration is employed. We show that individual pole masses exhibit similar scale-dependence regardless of the procedure, but that the leading scale-dependent terms cancel when computing the mass splitting if and only if the non-iterative procedure is employed. We show that this behaviour persists at two-loop order: the precision of the mass splitting improves in the non-iterative approach, but our results suggest that higher-order corrections do not reduce the uncertainty in the iterative calculation enough to resolve the problem at two-loop order. We conclude that the iterative procedure should not be used for computing pole masses in situations where electroweak mass splittings are phenomenologically relevant.
• ### Two-loop mass splittings in electroweak multiplets: winos and minimal dark matter(1712.00968)

April 1, 2018 hep-ph
The radiatively-induced splitting of masses in electroweak multiplets is relevant for both collider phenomenology and dark matter. Precision two-loop corrections of $\mathcal{O}$(MeV) to the triplet mass splitting in the wino limit of the minimal supersymmetric standard model can affect particle lifetimes by up to $40\%$. We improve on previous two-loop self-energy calculations for the wino model by obtaining consistent input parameters to the calculation via two-loop renormalisation-group running, and including the effect of finite light quark masses. We also present the first two-loop calculation of the mass splitting in an electroweak fermionic quintuplet, corresponding to the viable form of minimal dark matter (MDM). We place significant constraints on the lifetimes of the charged and doubly-charged fermions in this model. We find that the two-loop mass splittings in the MDM quintuplet are not constant in the large-mass limit, as might naively be expected from the triplet calculation. This is due to the influence of the additional heavy fermions in loop corrections to the gauge boson propagators.
• ### DarkBit: A GAMBIT module for computing dark matter observables and likelihoods(1705.07920)

Jan. 9, 2018 hep-ph, astro-ph.CO, astro-ph.IM
We introduce DarkBit, an advanced software code for computing dark matter constraints on various extensions to the Standard Model of particle physics, comprising both new native code and interfaces to external packages. This release includes a dedicated signal yield calculator for gamma-ray observations, which significantly extends current tools by implementing a cascade decay Monte Carlo, as well as a dedicated likelihood calculator for current and future experiments (gamlike). This provides a general solution for studying complex particle physics models that predict dark matter annihilation to a multitude of final states. We also supply a direct detection package that models a large range of direct detection experiments (DDcalc), and provides the corresponding likelihoods for arbitrary combinations of spin-independent and spin-dependent scattering processes. Finally, we provide custom relic density routines along with interfaces to DarkSUSY, micrOMEGAs, and the neutrino telescope likelihood package nuLike. DarkBit is written in the framework of the Global And Modular Beyond the Standard Model Inference Tool (GAMBIT), providing seamless integration into a comprehensive statistical fitting framework that allows users to explore new models with both particle and astrophysics constraints, and a consistent treatment of systematic uncertainties. In this paper we describe its main functionality, provide a guide to getting started quickly, and show illustrative examples for results obtained with DarkBit (both as a standalone tool and as a GAMBIT module). This includes a quantitative comparison between two of the main dark matter codes (DarkSUSY and micrOMEGAs), and application of DarkBit's advanced direct and indirect detection routines to a simple effective dark matter model.
• ### GAMBIT: The Global and Modular Beyond-the-Standard-Model Inference Tool(1705.07908)

Nov. 3, 2017 hep-ph, astro-ph.CO
We describe the open-source global fitting package GAMBIT: the Global And Modular Beyond-the-Standard-Model Inference Tool. GAMBIT combines extensive calculations of observables and likelihoods in particle and astroparticle physics with a hierarchical model database, advanced tools for automatically building analyses of essentially any model, a flexible and powerful system for interfacing to external codes, a suite of different statistical methods and parameter scanning algorithms, and a host of other utilities designed to make scans faster, safer and more easily-extendible than in the past. Here we give a detailed description of the framework, its design and motivation, and the current models and other specific components presently implemented in GAMBIT. Accompanying papers deal with individual modules and present first GAMBIT results. GAMBIT can be downloaded from gambit.hepforge.org.
• ### Comparison of statistical sampling methods with ScannerBit, the GAMBIT scanning module(1705.07959)

We introduce ScannerBit, the statistics and sampling module of the public, open-source global fitting framework GAMBIT. ScannerBit provides a standardised interface to different sampling algorithms, enabling the use and comparison of multiple computational methods for inferring profile likelihoods, Bayesian posteriors, and other statistical quantities. The current version offers random, grid, raster, nested sampling, differential evolution, Markov Chain Monte Carlo (MCMC) and ensemble Monte Carlo samplers. We also announce the release of a new standalone differential evolution sampler, Diver, and describe its design, usage and interface to ScannerBit. We subject Diver and three other samplers (the nested sampler MultiNest, the MCMC GreAT, and the native ScannerBit implementation of the ensemble Monte Carlo algorithm T-Walk) to a battery of statistical tests. For this we use a realistic physical likelihood function, based on the scalar singlet model of dark matter. We examine the performance of each sampler as a function of its adjustable settings, and the dimensionality of the sampling problem. We evaluate performance on four metrics: optimality of the best fit found, completeness in exploring the best-fit region, number of likelihood evaluations, and total runtime. For Bayesian posterior estimation at high resolution, T-Walk provides the most accurate and timely mapping of the full parameter space. For profile likelihood analysis in less than about ten dimensions, we find that Diver and MultiNest score similarly in terms of best fit and speed, outperforming GreAT and T-Walk; in ten or more dimensions, Diver substantially outperforms the other three samplers on all metrics.
• ### Neutrino telescope searches for dark matter in the Sun(1710.05190)

Oct. 14, 2017 hep-ph, astro-ph.CO
I give a brief review of a few recent developments and future directions in the search for dark matter using high-energy neutrinos from the Sun. This includes the ability to recast neutrino telescope limits on nuclear scattering of dark matter to arbitrary new theories, and new calculations of the solar atmospheric background relevant to such searches. I also touch on applications to global searches for new physics, and prospects for improving searches for asymmetric dark matter in the Sun.
• In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.
• ### Heating of galactic gas by dark matter annihilation in ultracompact minihalos(1611.08619)

May 26, 2017 hep-ph, astro-ph.CO, astro-ph.HE
The existence of substructure in halos of annihilating dark matter would be expected to substantially boost the rate at which annihilation occurs. Ultracompact minihalos of dark matter (UCMHs) are one of the more extreme examples of this. The boosted annihilation can inject significant amounts of energy into the gas of a galaxy over its lifetime. Here we determine the impact of the boost factor from UCMH substructure on the heating of galactic gas in a Milky Way-type galaxy, by means of N-body simulation. If $1\%$ of the dark matter exists as UCMHs, the corresponding boost factor can be of order $10^5$. For reasonable values of the relevant parameters (annihilation cross section $3\times10^{-26} ~\textrm{cm}^3~ \textrm{s}^{-1}$, dark matter mass 100 GeV, 10% heating efficiency), we show that the presence of UCMHs at the 0.1% level would inject enough energy to eject significant amounts of gas from the halo, potentially preventing star formation within $\sim$1 kpc of the halo centre.
• ### ColliderBit: a GAMBIT module for the calculation of high-energy collider observables and likelihoods(1705.07919)

April 2, 2019 hep-ph, hep-ex
We describe ColliderBit, a new code for the calculation of high energy collider observables in theories of physics beyond the Standard Model (BSM). ColliderBit features a generic interface to BSM models, a unique parallelised Monte Carlo event generation scheme suitable for large-scale supercomputer applications, and a number of LHC analyses, covering a reasonable range of the BSM signatures currently sought by ATLAS and CMS. ColliderBit also calculates likelihoods for Higgs sector observables, and LEP searches for BSM particles. These features are provided by a combination of new code unique to ColliderBit, and interfaces to existing state-of-the-art public codes. ColliderBit is both an important part of the GAMBIT framework for BSM inference, and a standalone tool for efficiently applying collider constraints to theories of new physics.
• ### FlavBit: A GAMBIT module for computing flavour observables and likelihoods(1705.07933)

May 22, 2017 hep-ph, hep-ex
Flavour physics observables are excellent probes of new physics up to very high energy scales. Here we present FlavBit, the dedicated flavour physics module of the global-fitting package GAMBIT. FlavBit includes custom implementations of various likelihood routines for a wide range of flavour observables, including detailed uncertainties and correlations associated with LHCb measurements of rare, leptonic and semileptonic decays of B and D mesons, kaons and pions. It provides a generalised interface to external theory codes such as SuperIso, allowing users to calculate flavour observables in and beyond the Standard Model, and then test them in detail against all relevant experimental data. We describe FlavBit and its constituent physics in some detail, then give examples from supersymmetry and effective field theory illustrating how it can be used both as a standalone library for flavour physics, and within GAMBIT.
• ### Evaporation and scattering of momentum- and velocity-dependent dark matter in the Sun(1703.07784)

March 22, 2017 hep-ph, astro-ph.CO
Dark matter with momentum- or velocity-dependent interactions with nuclei has shown significant promise for explaining the so-called Solar Abundance Problem, a longstanding discrepancy between solar spectroscopy and helioseismology. The best-fit models are all rather light, typically with masses in the range of 3-5 GeV. This is exactly the mass range where dark matter evaporation from the Sun can be important, but to date no detailed calculation of the evaporation of such models has been performed. Here we carry out this calculation, for the first time including arbitrary velocity- and momentum-dependent interactions, thermal effects, and a completely general treatment valid from the optically thin limit all the way through to the optically thick regime. We find that depending on the dark matter mass, interaction strength and type, the mass below which evaporation is relevant can vary from 1 to 4 GeV. This has the effect of weakening some of the better-fitting solutions to the Solar Abundance Problem, but also improving a number of others. As a by-product, we also provide an improved derivation of the capture rate that takes into account thermal and optical depth effects, allowing the standard result to be smoothly matched to the well-known saturation limit.
• ### Effect of electromagnetic dipole dark matter on energy transport in the solar interior(1610.06737)

Feb. 28, 2017 hep-ph, astro-ph.SR
In recent years, a revised set of solar abundances has led to a discrepancy in the sound-speed profile between helioseismology and theoretical solar models. Conventional solutions require additional mechanisms for energy transport within the Sun. Vincent et al. have recently suggested that dark matter with a momentum or velocity dependent cross section could provide a solution. In this work, we consider three models of dark matter with such cross sections and their effect on the stellar structure. In particular, the three models incorporate dark matter particles interacting through an electromagnetic dipole moment: an electric dipole, a magnetic dipole or an anapole. Each model is implemented in the \texttt{DarkStec} stellar evolution program, which incorporates the effects of dark matter capture and heat transport within the solar interior. We show that dark matter with an anapole moment of $\sim1\mathrm{GeV}^{-2}$ or magnetic dipole moment of $\sim10^{-3}\mu_p$ can improve the sound-speed profile, small frequency separations and convective zone radius with respect to the Standard Solar Model. However, the required dipole moments are strongly excluded by direct detection experiments.
• ### Towards the next generation of simplified Dark Matter models(1607.06680)

Feb. 25, 2017 hep-ph, hep-ex
This White Paper is an input to the ongoing discussion about the extension and refinement of simplified Dark Matter (DM) models. Based on two concrete examples, we show how existing simplified DM models (SDMM) can be extended to provide a more accurate and comprehensive framework to interpret and characterise collider searches. In the first example we extend the canonical SDMM with a scalar mediator to include mixing with the Higgs boson. We show that this approach not only provides a better description of the underlying kinematic properties that a complete model would possess, but also offers the option of using this more realistic class of scalar mixing models to compare and combine consistently searches based on different experimental signatures. The second example outlines how a new physics signal observed in a visible channel can be connected to DM by extending a simplified model including effective couplings. This discovery scenario uses the recently observed excess in the high-mass diphoton searches of ATLAS and CMS for a case study to show that such a pragmatic approach can aid the experimental search programme to verify/falsify a potential signal and to study its underlying nature. In the next part of the White Paper we outline other interesting options for SDMM that could be studied in more detail in the future. Finally, we discuss important aspects of supersymmetric models for DM and how these could help to develop of more complete SDMM.
• ### Pippi - painless parsing, post-processing and plotting of posterior and likelihood samples(1206.2245)

Interpreting samples from likelihood or posterior probability density functions is rarely as straightforward as it seems it should be. Producing publication-quality graphics of these distributions is often similarly painful. In this short note I describe pippi, a simple, publicly-available package for parsing and post-processing such samples, as well as generating high-quality PDF graphics of the results. Pippi is easily and extensively configurable and customisable, both in its options for parsing and post-processing samples, and in the visual aspects of the figures it produces. I illustrate some of these using an existing supersymmetric global fit, performed in the context of a gamma-ray search for dark matter. Pippi can be downloaded and followed at http://github.com/patscott/pippi .
• ### Ultracompact minihalos as probes of inflationary cosmology(1512.04597)

Oct. 20, 2016 hep-ph, gr-qc, astro-ph.CO
Cosmological inflation generates primordial density perturbations on all scales, including those far too small to contribute to the cosmic microwave background. At these scales, isolated ultracompact minihalos of dark matter can form well before standard structure formation, if the perturbations have sufficient amplitude. Minihalos affect pulsar timing data and are potentially bright sources of gamma rays. The resulting constraints significantly extend the observable window of inflation in the presence of cold dark matter, coupling two of the key problems in modern cosmology.
• ### A realistic assessment of the CTA sensitivity to dark matter annihilation(1408.4131)

Sept. 22, 2016 hep-ph, astro-ph.CO, astro-ph.HE
We estimate the sensitivity of the upcoming CTA gamma-ray telescope to DM annihilation at the Galactic centre, improving on previous analyses in a number of significant ways. First, we perform a detailed analyses of all backgrounds, including diffuse astrophysical emission for the first time in a study of this type. Second, we present a statistical framework for including systematic errors and estimate the consequent degradation in sensitivity. These errors may come from e.g. event reconstruction, Monte Carlo determination of the effective area or uncertainty in atmospheric conditions. Third, we show that performing the analysis on a set of suitably optimised regions of interest makes it possible to partially compensate for the degradation in sensitivity caused by systematics and diffuse emission. To probe dark matter with the canonical thermal annihilation cross-section, CTA systematics like non-uniform variations in acceptance over a single field of view must be kept below the 0.3% level, unless the dark matter density rises more steeply in the centre of the Galaxy than predicted by a typical Navarro-Frenk-White or Einasto profile. For a contracted $r^{-1.3}$ profile, and systematics at the 1% level, CTA can probe annihilation to $b\bar{b}$ at the canonical thermal level for dark matter masses between 100 GeV and 10 TeV.
• ### Combined analysis of effective Higgs portal dark matter models(1512.06458)

June 14, 2016 hep-ph, astro-ph.CO
We combine and extend the analyses of effective scalar, vector, Majorana and Dirac fermion Higgs portal models of dark matter (DM), in which DM couples to the Standard Model (SM) Higgs boson via an operator of the form $\mathcal{O}_{\textrm{DM}}\, H^\dagger H$. For the fermion models, we take an admixture of scalar $\overline{\psi} \psi$ and pseudoscalar $\overline{\psi} i\gamma_5 \psi$ interaction terms. For each model, we apply constraints on the parameter space based on the Planck measured DM relic density and the LHC limits on the Higgs invisible branching ratio. For the first time, we perform a consistent study of the indirect detection prospects for these models based on the WMAP7/Planck observations of the cosmic microwave background, a combined analysis of 15 dwarf spheroidal galaxies by Fermi-LAT and the upcoming Cherenkov Telescope Array (CTA). We also perform a correct treatment of the momentum-dependent direct search cross section that arises from the pseudoscalar interaction term in the fermionic DM theories. We find, in line with previous studies, that current and future direct search experiments such as LUX and XENON1T can exclude much of the parameter space, and we demonstrate that a joint observation in both indirect and direct searches is possible for high mass weakly interacting massive particles. In the case of a pure pseudoscalar interaction of a fermionic DM candidate, future gamma-ray searches are the only class of experiment capable of probing the high mass range of the theory.
• ### Generalised form factor dark matter in the Sun(1504.04378)

May 24, 2016 hep-ph, astro-ph.CO, astro-ph.SR
We study the effects of energy transport in the Sun by asymmetric dark matter with momentum and velocity-dependent interactions, with an eye to solving the decade-old Solar Abundance Problem. We study effective theories where the dark matter-nucleon scattering cross-section goes as $v_{\rm rel}^{2n}$ and $q^{2n}$ with $n = -1, 0, 1$ or $2$, where $v_{\rm rel}$ is the dark matter-nucleon relative velocity and $q$ is the momentum exchanged in the collision. Such cross-sections can arise generically as leading terms from the most basic nonstandard DM-quark operators. We employ a high-precision solar simulation code to study the impact on solar neutrino rates, the sound speed profile, convective zone depth, surface helium abundance and small frequency separations. We find that the majority of models that improve agreement with the observed sound speed profile and depth of the convection zone also reduce neutrino fluxes beyond the level that can be reasonably accommodated by measurement and theory errors. However, a few specific points in parameter space yield a significant overall improvement. A 3-5 GeV DM particle with $\sigma_{SI} \propto q^2$ is particularly appealing, yielding more than a $6\sigma$ improvement with respect to standard solar models, while being allowed by direct detection and collider limits. We provide full analytical capture expressions for $q$- and $v_{\rm rel}$-dependent scattering, as well as complete likelihood tables for all models.
• ### Investigating dark matter substructure with pulsar timing: I. Constraints on ultracompact minihalos(1509.02938)

Nov. 19, 2015 gr-qc, astro-ph.CO, astro-ph.GA
Small-scale dark matter structure within the Milky Way is expected to affect pulsar timing. The change in gravitational potential induced by a dark matter halo passing near the line of sight to a pulsar would produce a varying delay in the light travel time of photons from the pulsar. Individual transits produce an effect that would either be too rare or too weak to be detected in 30-year pulsar observations. However, a population of dark matter subhalos would be expected to produce a detectable effect on the measured properties of pulsars if the subhalos constitute a significant fraction of the total halo mass. The effect is to increase the dispersion of measured period derivatives across the pulsar population. By statistical analysis of the ATNF pulsar catalogue, we place an upper limit on this dispersion of $\log \sigma_{\dot{P}} \leq -17.05$. We use this to place strong upper limits on the number density of ultracompact minihalos within the Milky Way. These limits are completely independent of the particle nature of dark matter.
• ### Investigating dark matter substructure with pulsar timing: II. Improved limits on small-scale cosmology(1509.02941)

Nov. 19, 2015 hep-ph, gr-qc, astro-ph.CO
Ultracompact Minihalos (UCMHs) have been proposed as a type of dark matter sub-structure seeded by large-amplitude primordial perturbations and topological defects. UCMHs are expected to survive to the present era, allowing constraints to be placed on their cosmic abundance using observations within our own Galaxy. Constraints on their number density can be linked to conditions in the early universe that impact structure formation, such as increased primordial power on small scales, generic weak non-Gaussianity, and the presence of cosmic strings. We use new constraints on the abundance of UCMHs from pulsar timing to place generalised limits on the parameters of each of these cosmological scenarios. At some scales, the limits are the strongest to date, exceeding those from dark matter annihilation. Our new limits have the added advantage of being independent of the particle nature of dark matter, as they are based only on gravitational effects.