Source author record

Pat Scott

Pat Scott appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

hep-ph astro-ph.CO astro-ph.SR astro-ph.HE astro-ph.GA hep-ex physics.data-an gr-qc astro-ph.EP astro-ph.IM astro-ph hep-th Biological Physics Computation cond-mat.dis-nn Neurons and Cognition nucl-ex

Catalog footprint

What is connected

63works

17topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Global fit of 2HDM with future collider results

In this work, we summarize a global fit study of Type-II two Higgs doublet models (2HDM), and explore the impact of future SM-like Higgs and Z-pole precision measurements on the allowed parameter space. The work is based on the study results of a global fit of 2HDMs with the tool GAMBIT, utilising various current constraints including theoretical constraints (unitarity, perturbativity and vacuum stability), Higgs searches at colliders, electroweak physics and flavour constraints. We further investigate the ability of future facilities, such as the HL-LHC, CEPC, ILC and FCC-ee to explore the 2HDM parameter space.

preprint2022arXiv

Simple and statistically sound recommendations for analysing physical theories

Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Our examples can be reproduced with the code publicly available at https://doi.org/10.5281/zenodo.4322283.

preprint2022arXiv

Simulation of energy transport by dark matter scattering in stars

Asymmetric dark matter (ADM) that is captured in stars can act as an efficient conductor of heat. Small ADM-induced changes in a star's temperature gradient are known to alter neutrino fluxes and asteroseismological signatures, erase convective cores and modify a star's main sequence lifetime. The Sun's proximity to us makes it an ideal laboratory for studying these effects. However, the two formalisms commonly used to parametrize such heat transport were developed over 30 years ago, and calibrated with a single set of simulations. What's more, both are based on assumptions that break down at the Knudsen transition, where heat transport is maximized. We construct a Monte Carlo simulation to exactly solve the Boltzmann collision equation, determining the steady-state distribution and luminosity carried in stars by ADM with cross sections that depend on velocity and momentum. We find that, although the established (Gould and Raffelt) formalism based on local thermal equilibrium does well for constant cross sections, the isothermal (Spergel and Press) method actually performs better across all models with a simple, universal rescaling function. Based on simulation results, we provide recommendations on the parametrization of DM heat transport in stellar evolution models.

preprint2021arXiv

CosmoBit: A GAMBIT module for computing cosmological observables and likelihoods

We introduce $\sf{CosmoBit}$, a module within the open-source $\sf{GAMBIT}$ software framework for exploring connections between cosmology and particle physics with joint global fits. $\sf{CosmoBit}$ provides a flexible framework for studying various scenarios beyond $Λ$CDM, such as models of inflation, modifications of the effective number of relativistic degrees of freedom, exotic energy injection from annihilating or decaying dark matter, and variations of the properties of elementary particles such as neutrino masses and the lifetime of the neutron. Many observables and likelihoods in $\sf{CosmoBit}$ are computed via interfaces to $\sf{AlterBBN}$, $\sf{CLASS}$, $\sf{DarkAges}$, $\sf{MontePython}$, $\sf{MultiModeCode}$, and $\sf{plc}$. This makes it possible to apply a wide range of constraints from large-scale structure, Type Ia supernovae, Big Bang Nucleosynthesis and the cosmic microwave background. Parameter scans can be performed using the many different statistical sampling algorithms available within the $\sf{GAMBIT}$ framework, and results can be combined with calculations from other $\sf{GAMBIT}$ modules focused on particle physics and dark matter. We include extensive validation plots and a first application to scenarios with non-standard relativistic degrees of freedom and neutrino temperature, showing that the corresponding constraint on the sum of neutrino masses is much weaker than in the standard scenario.

preprint2021arXiv

The GAMBIT Universal Model Machine: from Lagrangians to Likelihoods

We introduce the GAMBIT Universal Model Machine (GUM), a tool for automatically generating code for the global fitting software framework GAMBIT, based on Lagrangian-level inputs. GUM accepts models written symbolically in FeynRules and SARAH formats, and can use either tool along with MadGraph and CalcHEP to generate GAMBIT model, collider, dark matter, decay and spectrum code, as well as GAMBIT interfaces to corresponding versions of SPheno, micrOMEGAs, Pythia and Vevacious (C++). In this paper we describe the features, methods, usage, pathways, assumptions and current limitations of GUM. We also give a fully worked example, consisting of the addition of a Majorana fermion simplified dark matter model with a scalar mediator to GAMBIT via GUM, and carry out a corresponding fit.

preprint2021arXiv

Thermal WIMPs and the Scale of New Physics: Global Fits of Dirac Dark Matter Effective Field Theories

We assess the status of a wide class of WIMP dark matter (DM) models in light of the latest experimental results using the global fitting framework $\textsf{GAMBIT}$. We perform a global analysis of effective field theory (EFT) operators describing the interactions between a gauge-singlet Dirac fermion and the Standard Model quarks, the gluons and the photon. In this bottom-up approach, we simultaneously vary the coefficients of 14 such operators up to dimension 7, along with the DM mass, the scale of new physics and several nuisance parameters. Our likelihood functions include the latest data from $\mathit{Planck}$, direct and indirect detection experiments, and the LHC. For DM masses below 100 GeV, we find that it is impossible to satisfy all constraints simultaneously while maintaining EFT validity at LHC energies. For new physics scales around 1 TeV, our results are influenced by several small excesses in the LHC data and depend on the prescription that we adopt to ensure EFT validity. Furthermore, we find large regions of viable parameter space where the EFT is valid and the relic density can be reproduced, implying that WIMPs can still account for the DM of the universe while being consistent with the latest data.

preprint2020arXiv

GAMBIT and its Application in the Search for Physics Beyond the Standard Model

The Global and Modular Beyond-Standard Model Inference Tool (GAMBIT) is an open source software framework for performing global statistical fits of particle physics models, using a wide range of particle and astroparticle data. In this review, we describe the design principles of the package, the statistical and sampling frameworks, the experimental data included, and the first two years of physics results generated with it. This includes supersymmetric models, axion theories, Higgs portal dark matter scenarios and an extension of the Standard Model to include right-handed neutrinos. Owing to the broad spectrum of physics scenarios tackled by the GAMBIT community, this also serves as a convenient, self-contained review of the current experimental and theoretical status of the most popular models of dark matter.

preprint2020arXiv

Global analyses of Higgs portal singlet dark matter models using GAMBIT

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.

preprint2020arXiv

Reinterpretation of LHC Results for New Physics: Status and Recommendations after Run 2

We report on the status of efforts to improve the reinterpretation of searches and measurements at the LHC in terms of models for new physics, in the context of the LHC Reinterpretation Forum. We detail current experimental offerings in direct searches for new particles, measurements, technical implementations and Open Data, and provide a set of recommendations for further improving the presentation of LHC results in order to better enable reinterpretation in the future. We also provide a brief description of existing software reinterpretation frameworks and recent global analyses of new physics that make use of the current data.

preprint2020arXiv

Self-consistent approach for measuring the energy spectra and composition of cosmic rays and determining the properties of hadronic interactions at high energy

Air showers, produced by the interaction of energetic cosmic rays with the atmosphere, are an excellent alternative to study particle physics at energies beyond any human-made particle accelerator. For that, it is necessary to identify first the mass composition of the primary cosmic ray (and its energy). None of the existing high energy interaction models have been able to reproduce coherently all air shower observables over the entire energy and zenith angle phase space. This is despite having tried all possible combinations for the cosmic ray mass composition. This proposal outlines a self-consistent strategy to study high energy particle interactions and identify the energy spectra and mass composition of cosmic rays. This strategy involves the participation of different particle accelerators and astrophysics experiments. This is important to cover the entire cosmic ray energy range and a larger phase-space of shower observables to probe the high energy interaction models.

preprint2019arXiv

Axion global fits with Peccei-Quinn symmetry breaking before inflation using GAMBIT

We present global fits of cosmologically stable axion-like particle and QCD axion models in the mass range 0.1 neV to 10 eV. We focus on the case where the Peccei-Quinn symmetry is broken before the end of inflation, such that the initial value of the axion field can be considered to be homogeneous throughout the visible Universe. We include detailed likelihood functions from light-shining-through-wall experiments, haloscopes, helioscopes, the axion relic density, horizontal branch stars, supernova 1987A, white dwarf cooling, and gamma-ray observations. We carry out both frequentist and Bayesian analyses, with and without the inclusion of white dwarf cooling. We explore the degree of fine-tuning present in different models and identify parameter regions where it is possible for QCD axion models to account for both the dark matter in the Universe and the cooling hints, comparing them to specific DFSZ- and KSVZ-type models. We find the most credible parameter regions, allowing us to set (prior-dependent) upper and lower bounds on the axion mass. Our analysis also suggests that QCD axions in this scenario most probably make up a non-negligible but sub-dominant component of the dark matter in the Universe.

preprint2019arXiv

ColliderBit: a GAMBIT module for the calculation of high-energy collider observables and likelihoods

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.

preprint2019arXiv

Combined collider constraints on neutralinos and charginos

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χ_1^0}, m_{\tildeχ_2^0}, m_{\tildeχ_3^0}, m_{\tildeχ_4^0})$ = (8-155, 103-260, 130-473, 219-502) GeV and chargino masses of $(m_{\tildeχ_1^{\pm}}, m_{\tildeχ_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σ$, 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$σ$. 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.

preprint2019arXiv

SpecBit, DecayBit and PrecisionBit: GAMBIT modules for computing mass spectra, particle decay rates and precision observables

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.

preprint2018arXiv

A global fit of the MSSM with GAMBIT

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.

preprint2018arXiv

DarkBit: A GAMBIT module for computing dark matter observables and likelihoods

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.

preprint2018arXiv

Global analyses of supersymmetry with GAMBIT

I provide a brief summary of the status of supersymmetric models with parameters defined at either the unification or weak scale, based on global fits using the GAMBIT framework.

preprint2018arXiv

Global fits of GUT-scale SUSY models with GAMBIT

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.

preprint2018arXiv

Impact of vacuum stability, perturbativity and XENON1T on global fits of $\mathbb{Z}_2$ and $\mathbb{Z}_3$ scalar singlet dark matter

Scalar singlet dark matter is one of the simplest and most predictive realisations of the WIMP (weakly-interacting massive particle) idea. Although the model is constrained from all directions by the latest experimental data, it still has viable regions of parameter space. Another compelling aspect of scalar singlets is their ability to stabilise the electroweak vacuum. Indeed, models of scalar dark matter are not low-energy effective theories, but can be valid all the way to the Planck scale. Using the GAMBIT framework, we present the first global fit to include both the low-energy experimental constraints and the theoretical constraints from UV physics, considering models with a scalar singlet charged under either a $\mathbb{Z}_2$ or a $\mathbb{Z}_3$ symmetry. We show that if the model is to satisfy all experimental constraints, completely stabilise the electroweak vacuum up to high scales, and also remain perturbative to those scales, one is driven to a relatively small region of parameter space. This region has a Higgs-portal coupling slightly less than 1, a dark matter mass of 1 to 2 TeV and a spin-independent nuclear scattering cross-section around 10$^{-45}$ cm$^2$.

preprint2018arXiv

Status of the scalar singlet dark matter model

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.

preprint2017arXiv

Comparison of statistical sampling methods with ScannerBit, the GAMBIT scanning module

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.

preprint2017arXiv

FlavBit: A GAMBIT module for computing flavour observables and likelihoods

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.

preprint2017arXiv

GAMBIT: The Global and Modular Beyond-the-Standard-Model Inference Tool

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.

preprint2016arXiv

A realistic assessment of the CTA sensitivity to dark matter annihilation

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.

preprint2016arXiv

Combined analysis of effective Higgs portal dark matter models

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ψ ψ$ and pseudoscalar $\overlineψ iγ_5 ψ$ 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.

preprint2016arXiv

Generalised form factor dark matter in the Sun

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 $σ_{SI} \propto q^2$ is particularly appealing, yielding more than a $6σ$ 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.

preprint2016arXiv

Implications of solar wind measurements for solar models and composition

We critically examine recent claims of a high solar metallicity by von Steiger \& Zurbuchen (2016) based on \textit{in situ} measurements of the solar wind, rather than the standard spectroscopically-inferred abundances (Asplund et al. 2009). We test the claim by Vagnozzi et al. (2016) that a composition based on the solar wind enables one to construct a standard solar model in agreement with helioseismological observations and thus solve the decades-old solar modelling problem. We show that, although some helioseismological observables are improved compared to models computed with spectroscopic abundances, most are in fact worse. The high abundance of refractory elements leads to an overproduction of neutrinos, with a predicted $^8$B flux that is nearly twice its observed value, and $^7$Be and CNO fluxes that are experimentally ruled out at high confidence. A combined likelihood analysis shows that models using the vSZ16 abundances fare worse than AGSS09 despite a higher metallicity. We also present astrophysical and spectroscopic arguments showing the vSZ16 composition to be an implausible representation of the solar interior, identifying the first ionisation potential effect in the outer solar atmosphere and wind as the likely culprit.

preprint2016arXiv

Pippi - painless parsing, post-processing and plotting of posterior and likelihood samples

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 .

preprint2016arXiv

Ultracompact minihalos as probes of inflationary cosmology

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.

preprint2016arXiv

Updated constraints on velocity and momentum-dependent asymmetric dark matter

We present updated constraints on dark matter models with momentum-dependent or velocity-dependent interactions with nuclei, based on direct detection and solar physics. We improve our previous treatment of energy transport in the solar interior by dark matter scattering, leading to significant changes in fits to many observables. Based on solar physics alone, DM with a spin-independent $q^{4}$ coupling provides the best fit to data, and a statistically satisfactory solution to the solar abundance problem. Once direct detection limits are accounted for however, the best solution is spin-dependent $v^2$ scattering with a reference cross-section of 10$^{-35}$ cm$^2$ (at a reference velocity of $v_0=220$ km s$^{-1}$), and a dark matter mass of about 5 GeV.

preprint2015arXiv

A possible indication of momentum-dependent asymmetric dark matter in the Sun

Broad disagreement persists between helioseismological observables and predictions of solar models computed with the latest surface abundances. Here we show that most of these problems can be solved by the presence of asymmetric dark matter coupling to nucleons as the square of the momentum $q$ exchanged in the collision. We compute neutrino fluxes, small frequency separations, surface helium abundances, sound speed profiles and convective zone depths for a number of models, showing more than a $6σ$ preference for $q^2$ models over others, and over the Standard Solar Model. The preferred mass (3\,GeV) and reference dark matter-nucleon cross-section ($10^{-37}$\,cm$^2$ at $q_0 = 40$\,MeV) are within the region of parameter space allowed by both direct detection and collider searches.

preprint2015arXiv

Constraints on cosmic strings from ultracompact minihalos

Cosmic strings are expected to form loops. These can act as seeds for accretion of dark matter, leading to the formation of ultracompact minihalos (UCMHs). We perform a detailed study of the accretion of dark matter onto cosmic string loops and compute the resulting mass distribution of UCMHs. We then apply observational limits on the present-day abundance of UCMHs to derive corresponding limits on the cosmic string tension $Gμ$. The bounds are strongly dependent upon the assumed distribution of loop velocities and their impacts on UCMH formation. Under the assumption that a loop can move up to a thousand times its own radius and still form a UCMH, we find a limit of $Gμ\le 5\times10^{-8}$. We show, in opposition to previous results, that strong limits on the cosmic string tension are not obtainable from UCMHs when more stringent (and realistic) requirements are placed on loop velocities.

preprint2015arXiv

Update on scalar singlet dark matter

One of the simplest models of dark matter is that where a scalar singlet field S comprises some or all of the dark matter, and interacts with the standard model through an HHSS coupling to the Higgs boson. We update the present limits on the model from LHC searches for invisible Higgs decays, the thermal relic density of S, and dark matter searches via indirect and direct detection. We point out that the currently allowed parameter space is on the verge of being significantly reduced with the next generation of experiments. We discuss the impact of such constraints on possible applications of scalar singlet dark matter, including a strong electroweak phase transition, and the question of vacuum stability of the Higgs potential at high scales.

preprint2014arXiv

Sensitivity of CTA to dark matter signals from the Galactic Center

(abridged) The Galactic Center is one of the most promising targets for indirect detection of dark matter with gamma rays. We investigate the sensitivity of the upcoming Cherenkov Telescope Array (CTA) to dark matter annihilation and decay in the Galactic Center. As the inner density profile of the Milky Way's dark matter halo is uncertain, we study the impact of the slope of the Galactic density profile, inwards of the Sun, on the prospects for detecting a dark matter signal with CTA. We find that the sensitivity achieved by CTA to annihilation signals is strongly dependent on the inner profile slope, whereas the dependence is more mild in the case of dark matter decay. Surprisingly, we find that the optimal choice of signal and background regions is virtually independent of the assumed density profile. For the fiducial case of a Navarro-Frenk-White profile, we find that CTA will be able to probe annihilation cross sections well below the canonical thermal relic value for dark matter masses from a few tens of GeV up to $\sim 5$ TeV for annihilation to $τ^{+}τ^{-}$, and will achieve only a slightly weaker sensitivity for annihilation to $b\bar{b}$ or $μ^{+}μ^{-}$. CTA will improve significantly on current sensitivity to annihilation signals for dark matter masses above $\sim 100$ GeV, covering parameter space that is complementary to that probed by searches with the Fermi Large Area Telescope. The interpretation of apparent excesses in the measured cosmic-ray electron and positron spectra as signals of dark matter decay will also be testable with CTA observations of the Galactic Center. We demonstrate that both for annihilation and for decay, including spectral information for hard channels (such as $μ^{+}μ^{-}$ and $τ^{+}τ^{-}$) leads to enhanced sensitivity for dark matter masses above $m_{\rm DM}\sim 200$ GeV.

preprint2014arXiv

The elemental composition of the Sun I. The intermediate mass elements Na to Ca

The composition of the Sun is an essential piece of reference data for astronomy, cosmology, astroparticle, space and geo-physics. This article, dealing with the intermediate-mass elements Na to Ca, is the first in a series describing the comprehensive re-determination of the solar composition. In this series we severely scrutinise all ingredients of the analysis across all elements, to obtain the most accurate, homogeneous and reliable results possible. We employ a highly realistic 3D hydrodynamic solar photospheric model, which has successfully passed an arsenal of observational diagnostics. To quantify systematic errors, we repeat the analysis with three 1D hydrostatic model atmospheres (MARCS, MISS and Holweger & Müller 1974) and a horizontally and temporally-averaged version of the 3D model ($\langle$3D$\rangle$). We account for departures from LTE wherever possible. We have scoured the literature for the best transition probabilities, partition functions, hyperfine and other data, and stringently checked all observed profiles for blends. Our final 3D+NLTE abundances are: $\logε_{\mathrm{Na}}=6.21\pm0.04$, $\logε_{\mathrm{Mg}}=7.59\pm0.04$, $\logε_{\mathrm{Al}}=6.43\pm0.04$, $\logε_{\mathrm{Si}}=7.51\pm0.03$, $\logε_{\mathrm{P}}=5.41\pm0.03$, $\log ε_{\mathrm{S}}=7.13\pm0.03$, $\logε_{\mathrm{K}}=5.04\pm0.05$ and $\logε_{\mathrm{Ca}}=6.32\pm0.03$. The uncertainties include both statistical and systematic errors. Our results are systematically smaller than most previous ones with the 1D semi-empirical Holweger & Müller model. The $\langle$3D$\rangle$ model returns abundances very similar to the full 3D calculations. This analysis provides a complete description and a slight update of the Na to Ca results presented in Asplund, Grevesse, Sauval & Scott (arXiv:0909.0948), with full details of all lines and input data.

preprint2014arXiv

The elemental composition of the Sun II. The iron group elements Sc to Ni

We redetermine the abundances of all iron group nuclei in the Sun, based on neutral and singly-ionised lines of Sc, Ti, V, Mn, Fe, Co and Ni in the solar spectrum. We employ a realistic 3D hydrodynamic model solar atmosphere, corrections for departures from local thermodynamic equilibrium (NLTE), stringent line selection procedures and high quality observational data. We have scoured the literature for the best quality oscillator strengths, hyperfine constants and isotopic separations available for our chosen lines. We find $\log ε_\mathrm{Sc}=3.16\pm0.04$, $\log ε_\mathrm{Ti}=4.93\pm0.04$, $\log ε_\mathrm{V}=3.89\pm0.08$, $\log ε_\mathrm{Cr}=5.62\pm0.04$, $\log ε_\mathrm{Mn}=5.42\pm0.04$, $\log ε_\mathrm{Fe}=7.47\pm0.04$, $\log ε_\mathrm{Co}=4.93\pm0.05$ and $\log ε_\mathrm{Ni}=6.20\pm0.04$. Our uncertainties factor in both statistical and systematic errors (the latter estimated for possible errors in the model atmospheres and NLTE line formation). The new abundances are generally in good agreement with the CI meteoritic abundances but with some notable exceptions. This analysis constitutes both a full exposition and a slight update of the preliminary results we presented in Asplund, Grevesse, Sauval & Scott (arXiv:0909.0948), including full line lists and details of all input data we employed.

preprint2014arXiv

The elemental composition of the Sun III. The heavy elements Cu to Th

We re-evaluate the abundances of the elements in the Sun from copper ($Z=29$) to thorium ($Z=90$). Our results are mostly based on neutral and singly-ionised lines in the solar spectrum. We use the latest 3D hydrodynamic solar model atmosphere, and in a few cases also correct for departures from local thermodynamic equilibrium (LTE) using non-LTE (NLTE) calculations performed in 1D. In order to minimise statistical and systematic uncertainties, we make stringent line selections, employ the highest-quality observational data and carefully assess oscillator strengths, hyperfine constants and isotopic separations available in the literature, for every line included in our analysis. Our results are typically in good agreement with the abundances in the most pristine meteorites, but there are some interesting exceptions. This analysis constitutes both a full exposition and a slight update of the relevant parts of the preliminary results we presented in Asplund, Grevesse, Sauval & Scott (arXiv:0909.0948), including full line lists and details of all input data that we have employed.

preprint2014arXiv

Thermal conduction by dark matter with velocity and momentum-dependent cross-sections

We use the formalism of Gould and Raffelt to compute the dimensionless thermal conduction coefficients for scattering of dark matter particles with standard model nucleons via cross-sections that depend on the relative velocity or momentum exchanged between particles. Motivated by models invoked to reconcile various recent results in direct detection, we explicitly compute the conduction coefficients $α$ and $κ$ for cross-sections that go as $v_{\rm rel}^2$, $v_{\rm rel}^4$, $v_{\rm rel}^{-2}$, $q^2$, $q^4$ and $q^{-2}$, where $v_{\rm rel}$ is the relative DM-nucleus velocity and $q$ is the momentum transferred in the collision. We find that a $v_{\rm rel}^{-2}$ dependence can significantly enhance energy transport from the inner solar core to the outer core. The same can true for any $q$-dependent coupling, if the dark matter mass lies within some specific range for each coupling. This effect can complement direct searches for dark matter; combining these results with state-of-the-art Solar simulations should greatly increase sensitivity to certain DM models. It also seems possible that the so-called Solar Abundance Problem could be resolved by enhanced energy transport in the solar core due to such velocity- or momentum-dependent scatterings.

preprint2013arXiv

Dark Matter CMB Constraints and Likelihoods for Poor Particle Physicists

The cosmic microwave background provides constraints on the annihilation and decay of light dark matter at redshifts between 100 and 1000, the strength of which depends upon the fraction of energy ending up in the form of electrons and photons. The resulting constraints are usually presented for a limited selection of annihilation and decay channels. Here we provide constraints on the annihilation cross section and decay rate, at discrete values of the dark matter mass, for all the annihilation and decay channels whose secondary spectra have been computed using PYTHIA in arXiv:1012.4515 ("PPPC4DMID: A Poor Particle Physicist Cookbook for Dark Matter Indirect Detection"), namely e, mu, tau, V -> e, V -> mu, V -> tau, u, d, s, c, b, t, gamma, g, W, Z and h. By interpolating in mass, these can be used to find the CMB constraints and likelihood functions from WMAP7 and Planck for a wide range of dark matter models, including those with annihilation or decay into a linear combination of different channels.

preprint2013arXiv

Detection of isolated population III stars with the James Webb Space Telescope

The first population III stars are predicted to form in minihalos at a redshift of approximately 10-30. The James Webb Space Telescope (JWST), tentatively scheduled for launch in 2018, will probably be able to detect some of the first galaxies, but whether it will also be able to detect the first stars remains more doubtful. Here, we explore the prospects of detecting an isolated population III star or a small cluster of population III stars down to redshift 2 in either lensed or unlensed fields. Our calculations are based on realistic stellar atmospheres and take into account the potential flux contribution from the surrounding HII region. We find that unlensed population III stars are beyond the reach of JWST, and that even lensed population III stars will be extremely difficult to detect. However, the main problem with the latter approach is not necessarily that the lensed stars are too faint, but that their surface number densities are too low. To detect even one population III star of 60 solar masses when pointing JWST through the galaxy cluster MACS J0717.5+3745, the lensing cluster with the largest Einstein radius detected so far, the cosmic star formation rate of population III stars would need to be approximately an order of magnitude higher than predicted by the most optimistic current models.

preprint2013arXiv

Hunting for dark halo substructure using submilliarcsecond-scale observations of macrolensed radio jets

Dark halo substructure may reveal itself through secondary, small-scale gravitational lensing effects on light sources that are macrolensed by a foreground galaxy. Here, we explore the prospects of using Very Long Baseline Interferometry (VLBI) observations of multiply-imaged quasar jets to search for submilliarcsecond-scale image distortions produced by various forms of dark substructures in the 1e3-1e8 Msolar mass range. We present lensing simulations relevant for the angular resolutions attainable with the existing European VLBI Network (EVN), the global VLBI array, and an upcoming observing mode in which the Atacama Large Millimeter Array (ALMA) is connected to the global VLBI array. While observations of this type would not be sensitive to standard cold dark matter subhalos, they can be used to detect more compact forms of halo substructure predicted in alternative structure formation scenarios. By mapping ~5 strongly lensed systems, it should be possible to detect or robustly rule out primordial black holes in the 1e3-1e6 Msolar mass range if they constitute >1% percent of the dark matter in these lenses. Ultracompact minihalos are harder to detect using this technique, but 1e6-1e8 Msolar ultracompact minihalos could in principle be detected if they constitute >10% of the dark matter.

preprint2013arXiv

Improved constraints on the primordial power spectrum at small scales from ultracompact minihalos

For a Gaussian spectrum of primordial density fluctuations, ultracompact minihalos (UCMHs) of dark matter are expected to be produced in much greater abundance than, e.g., primordial black holes. Forming shortly after matter-radiation equality, these objects would develop very dense and spiky dark matter profiles. In the standard scenario where dark matter consists of thermally-produced, weakly-interacting massive particles, UCMHs could thus appear as highly luminous gamma-ray sources, or leave an imprint in the cosmic microwave background by changing the reionisation history of the Universe. We derive corresponding limits on the cosmic abundance of UCMHs at different epochs, and translate them into constraints on the primordial power spectrum. We find the resulting constraints to be quite severe, especially at length scales much smaller than what can be directly probed by the cosmic microwave background or large-scale structure observations. We use our results to provide an updated compilation of the best available constraints on the power of density fluctuations on all scales, ranging from the present-day horizon to scales more than 20 orders of magnitude smaller.

preprint2013arXiv

Light bosons and photospheric solutions to the solar abundance problem

It is well known that current spectroscopic determinations of the chemical composition of the Sun are starkly at odds with the metallicity implied by helioseismology. We investigate whether the discrepancy may be due to conversion of photons to a new light boson in the solar photosphere. We examine the impact of particles with axion-like interactions with the photon on the inferred photospheric abundances, showing that resonant axion-photon conversion is not possible in the region of the solar atmosphere in which line-formation occurs. Although non-resonant conversion in the line-forming regions can in principle impact derived abundances, constraints from axion-photon conversion experiments rule out the couplings necessary for these effects to be detectable. We show that this extends to hidden photons and chameleons (which would exhibit similar phenomenological behaviour), ruling out known theories of new light bosons as photospheric solutions to the solar abundance problem.

preprint2013arXiv

Number Counts and Non-Gaussianity

We describe a general procedure for using number counts of any object to constrain the probability distribution of the primordial fluctuations, allowing for generic weak non-Gaussianity. We apply this procedure to use limits on the abundance of primordial black holes and dark matter ultracompact minihalos (UCMHs) to characterize the allowed statistics of primordial fluctuations on very small scales. We present constraints on the power spectrum and the amplitude of the skewness for two different families of non-Gaussian distributions, distinguished by the relative importance of higher moments. Although primordial black holes probe the smallest scales, ultracompact minihalos provide significantly stronger constraints on the power spectrum and so are more likely to eventually provide small-scale constraints on non-Gaussianity.

preprint2013arXiv

Sensitivity of IceCube-DeepCore to neutralino dark matter in the MSSM-25

We analyse the sensitivity of IceCube-DeepCore to annihilation of neutralino dark matter in the solar core, generated within a 25 parameter version of the minimally supersymmetric standard model (MSSM-25). We explore the 25-dimensional parameter space using scanning methods based on importance sampling and using DarkSUSY 5.0.6 to calculate observables. Our scans produced a database of 6.02 million parameter space points with neutralino dark matter consistent with the relic density implied by WMAP 7-year data, as well as with accelerator searches. We performed a model exclusion analysis upon these points using the expected capabilities of the IceCube-DeepCore Neutrino Telescope. We show that IceCube-DeepCore will be sensitive to a number of models that are not accessible to direct detection experiments such as SIMPLE, COUPP and XENON100, indirect detection using Fermi-LAT observations of dwarf spheroidal galaxies, nor to current LHC searches.

preprint2012arXiv

Constraints on small-scale cosmological perturbations from gamma-ray searches for dark matter

Events like inflation or phase transitions can produce large density perturbations on very small scales in the early Universe. Probes of small scales are therefore useful for e.g. discriminating between inflationary models. Until recently, the only such constraint came from non-observation of primordial black holes (PBHs), associated with the largest perturbations. Moderate-amplitude perturbations can collapse shortly after matter-radiation equality to form ultracompact minihalos (UCMHs) of dark matter, in far greater abundance than PBHs. If dark matter self-annihilates, UCMHs become excellent targets for indirect detection. Here we discuss the gamma-ray fluxes expected from UCMHs, the prospects of observing them with gamma-ray telescopes, and limits upon the primordial power spectrum derived from their non-observation by the Fermi Large Area Space Telescope.

preprint2012arXiv

Fundamental statistical limitations of future dark matter direct detection experiments

We discuss irreducible statistical limitations of future ton-scale dark matter direct detection experiments. We focus in particular on the coverage of confidence intervals, which quantifies the reliability of the statistical method used to reconstruct the dark matter parameters, and the bias of the reconstructed parameters. We study 36 benchmark dark matter models within the reach of upcoming ton-scale experiments. We find that approximate confidence intervals from a profile-likelihood analysis exactly cover or over-cover the true values of the WIMP parameters, and are hence conservative. We evaluate the probability that unavoidable statistical fluctuations in the data might lead to a biased reconstruction of the dark matter parameters, or large uncertainties on the reconstructed parameter values. We show that this probability can be surprisingly large, even for benchmark models leading to a large event rate of order a hundred counts. We find that combining data sets from two different targets leads to improved coverage properties, as well as a substantial reduction of statistical bias and uncertainty on the dark matter parameters.

preprint2012arXiv

Quantifying impacts of short-term plasticity on neuronal information transfer

Short-term changes in efficacy have been postulated to enhance the ability of synapses to transmit information between neurons, and within neuronal networks. Even at the level of connections between single neurons, direct confirmation of this simple conjecture has proven elusive. By combining paired-cell recordings, realistic synaptic modelling and information theory, we provide evidence that short-term plasticity can not only improve, but also reduce information transfer between neurons. We focus on a concrete example in rat neocortex, but our results may generalise to other systems. When information is contained in the timings of individual spikes, we find that facilitation, depression and recovery affect information transmission in proportion to their impacts upon the probability of neurotransmitter release. When information is instead conveyed by mean spike rate only, the influences of short-term plasticity critically depend on the range of spike frequencies that the target network can distinguish (its effective dynamic range). Our results suggest that to efficiently transmit information, the brain must match synaptic type, coding strategy and network connectivity during development and behaviour.

preprint2011arXiv

Dark stars: structure, evolution and impacts upon the high-redshift Universe

The most compelling and popular models for dark matter predict that it should congregate and annihilate in stellar cores. Stars where annihilation contributes substantially to the total energy budget look very different to those with which we are familiar. Here I explain the general features of stars modified by dark matter annihilation with the help of a series of grids of 'dark' stellar evolutionary models, and describe the public code with which they were computed. I go on to discuss possible impacts of dark stars on the high-redshift Universe, including the history of reionisation. The preliminary reionisation calculations reproduced here are based on dedicated models for dark star atmospheres, and for the stellar populations to which dark stars would belong.

preprint2011arXiv

How well will ton-scale dark matter direct detection experiments constrain minimal supersymmetry?

Weakly interacting massive particles (WIMPs) are amongst the most interesting dark matter (DM) candidates. Many DM candidates naturally arise in theories beyond the standard model (SM) of particle physics, like weak-scale supersymmetry (SUSY). Experiments aim to detect WIMPs by scattering, annihilation or direct production, and thereby determine the underlying theory to which they belong, along with its parameters. Here we examine the prospects for further constraining the Constrained Minimal Supersymmetric Standard Model (CMSSM) with future ton-scale direct detection experiments. We consider ton-scale extrapolations of three current experiments: CDMS, XENON and COUPP, with 1000 kg-years of raw exposure each. We assume energy resolutions, energy ranges and efficiencies similar to the current versions of the experiments, and include backgrounds at target levels. Our analysis is based on full likelihood constructions for the experiments. We also take into account present uncertainties on hadronic matrix elements for neutralino-quark couplings, and on halo model parameters. We generate synthetic data based on four benchmark points and scan over the CMSSM parameter space using nested sampling. We construct both Bayesian posterior PDFs and frequentist profile likelihoods for the model parameters, as well as the mass and various cross-sections of the lightest neutralino. Future ton-scale experiments will help substantially in constraining supersymmetry, especially when results of experiments primarily targeting spin-dependent nuclear scattering are combined with those directed more toward spin-independent interactions.

preprint2011arXiv

Impacts of Dark Stars on Reionization and Signatures in the Cosmic Microwave Background

We perform a detailed and systematic investigation of the possible impacts of dark stars upon the reionization history of the Universe, and its signatures in the cosmic microwave background (CMB). We compute hydrogen reionization histories, CMB optical depths and anisotropy power spectra for a range of stellar populations including dark stars. If dark stars capture large amounts of dark matter via nuclear scattering, reionization can be substantially delayed, leading to decreases in the integrated optical depth to last scattering and large-scale power in the EE polarization power spectrum. Using the integrated optical depth observed by WMAP7, in our canonical reionization model we rule out the section of parameter space where dark stars with high scattering-induced capture rates tie up more than ~90% of all the first star-forming baryons, and live for over ~250 Myr. When nuclear scattering delivers only moderate amounts of dark matter, reionization can instead be sped up slightly, modestly increasing the CMB optical depth. If dark stars do not obtain any dark matter via nuclear scattering, effects upon reionization and the CMB are negligible. The effects of dark stars upon reionization and its CMB markers can be largely mimicked or compensated for by changes in the existing parameters of reionization models, making dark stars difficult to disentangle from astrophysical uncertainties, but also widening the range of standard parameters in reionization models that can be made consistent with observations.

preprint2011arXiv

Implications for constrained supersymmetry of combined H.E.S.S. observations of dwarf galaxies, the Galactic halo and the Galactic Centre

In order to place limits on dark matter (DM) properties using $γ$-ray observations, previous analyses have often assumed a very simple parametrisation of the $γ$-ray annihilation yield; typically, it has been assumed that annihilation proceeds through a single channel only. In realistic DM models, annihilation may occur into many different final states, making this quite a rough ansatz. With additional processes like virtual internal bremsstrahlung and final state radiation, this ansatz becomes even more incorrect, and the need for scans of explicit model parameter spaces becomes clear. Here we present scans of the parameter space of the Constrained Minimal Supersymmetric Standard Model (CMSSM), considering $γ$-ray spectra from three dwarf galaxies, the Galactic Centre region and the broader Galactic halo, as obtained with the High-Energy Stereoscopic System (H.E.S.S.). We present a series of likelihood scans combining the H.E.S.S. data with other experimental results. We show that observations of the Sagittarius, Carina and Sculptor dwarf galaxies disfavour the coannihilation region of the CMSSM and models with large annihilation cross-sections. This is true even under reasonable assumptions about the DM density profiles, and constitutes the strongest constraint to date on coannihilation models within the CMSSM. The Galactic halo has a similar, but weaker, effect. The Galactic Centre search is complicated by a strong (unknown) $γ$-ray source, and we see no effect on the CMSSM parameter space when assuming a realistic Galactic Centre DM density profile.

preprint2011arXiv

Observational constraints on supermassive dark stars

Some of the first stars could be cooler and more massive than standard stellar models would suggest, due to the effects of dark matter annihilation in their cores. It has recently been argued that such objects may attain masses in the 10^4--10^7 solar mass range, and that such supermassive dark stars should be within reach of the upcoming James Webb Space Telescope. Notwithstanding theoretical difficulties with this proposal, we argue here that some of these objects should also be readily detectable with both the Hubble Space Telescope and ground-based 8--10 m class telescopes. Existing survey data already place strong constraints on 10^7 solar mass dark stars at z~10. We show that such objects must be exceedingly rare or short-lived to have avoided detection.

preprint2011arXiv

Searches for Particle Dark Matter: An Introduction

The identity of dark matter is one of the key outstanding problems in both particle and astrophysics. In this thesis, I describe a number of complementary searches for particle dark matter. I discuss how the impact of dark matter on stars can constrain its interaction with nuclei, focussing on main sequence stars close to the Galactic Centre, and on the first stars as seen through the upcoming James Webb Space Telescope. The mass and annihilation cross-section of dark matter particles can be probed with searches for gamma rays produced in astronomical targets. Dwarf galaxies and ultracompact, primordially-produced dark matter minihalos turn out to be especially promising in this respect. I illustrate how the results of these searches can be combined with constraints from accelerators and cosmology to produce a single global fit to all available data. Global fits in supersymmetry turn out to be quite technically demanding, even with the simplest predictive models and the addition of complementary data from a bevy of astronomical and terrestrial experiments; I show how genetic algorithms can help in overcoming these challenges.

preprint2011arXiv

Statistical coverage for supersymmetric parameter estimation: a case study with direct detection of dark matter

Models of weak-scale supersymmetry offer viable dark matter (DM) candidates. Their parameter spaces are however rather large and complex, such that pinning down the actual parameter values from experimental data can depend strongly on the employed statistical framework and scanning algorithm. In frequentist parameter estimation, a central requirement for properly constructed confidence intervals is that they cover true parameter values, preferably at exactly the stated confidence level when experiments are repeated infinitely many times. Since most widely-used scanning techniques are optimised for Bayesian statistics, one needs to assess their abilities in providing correct confidence intervals in terms of the statistical coverage. Here we investigate this for the Constrained Minimal Supersymmetric Standard Model (CMSSM) when only constrained by data from direct searches for dark matter. We construct confidence intervals from one-dimensional profile likelihoods and study the coverage by generating several pseudo-experiments for a few benchmark sets of pseudo-true parameters. We use nested sampling to scan the parameter space and evaluate the coverage for the benchmarks when either flat or logarithmic priors are imposed on gaugino and scalar mass parameters. The sampling algorithm has been used in the configuration usually adopted for exploration of the Bayesian posterior. We observe both under- and over-coverage, which in some cases vary quite dramatically when benchmarks or priors are modified. We show how most of the variation can be explained as the impact of explicit priors as well as sampling effects, where the latter are indirectly imposed by physicality conditions. For comparison, we also evaluate the coverage for Bayesian credible intervals, and observe significant under-coverage in those cases.

preprint2010arXiv

A Profile Likelihood Analysis of the Constrained MSSM with Genetic Algorithms

The Constrained Minimal Supersymmetric Standard Model (CMSSM) is one of the simplest and most widely-studied supersymmetric extensions to the standard model of particle physics. Nevertheless, current data do not sufficiently constrain the model parameters in a way completely independent of priors, statistical measures and scanning techniques. We present a new technique for scanning supersymmetric parameter spaces, optimised for frequentist profile likelihood analyses and based on Genetic Algorithms. We apply this technique to the CMSSM, taking into account existing collider and cosmological data in our global fit. We compare our method to the MultiNest algorithm, an efficient Bayesian technique, paying particular attention to the best-fit points and implications for particle masses at the LHC and dark matter searches. Our global best-fit point lies in the focus point region. We find many high-likelihood points in both the stau co-annihilation and focus point regions, including a previously neglected section of the co-annihilation region at large m_0. We show that there are many high-likelihood points in the CMSSM parameter space commonly missed by existing scanning techniques, especially at high masses. This has a significant influence on the derived confidence regions for parameters and observables, and can dramatically change the entire statistical inference of such scans.

preprint2010arXiv

Finding high-redshift dark stars with the James Webb Space Telescope

The first stars in the history of the Universe are likely to form in the dense central regions of 10^5-10^6 Msolar cold dark matter halos at z=10-50. The annihilation of dark matter particles in these environments may lead to the formation of so-called dark stars, which are predicted to be cooler, larger, more massive and potentially more long-lived than conventional population III stars. Here, we investigate the prospects of detecting high-redshift dark stars with the upcoming James Webb Space Telescope (JWST). We find that dark stars at z>6 are intrinsically too faint to be detected by JWST. However, by exploiting foreground galaxy clusters as gravitational telescopes, certain varieties of cool (Teff < 30000 K) dark stars should be within reach at redshifts up to z=10. If the lifetimes of dark stars are sufficiently long, many such objects may also congregate inside the first galaxies. We demonstrate that this could give rise to peculiar features in the integrated spectra of galaxies at high redshifts, provided that dark stars make up at least 1 percent of the total stellar mass in such objects.

preprint2010arXiv

Gamma rays from ultracompact primordial dark matter minihalos

Ultracompact minihalos have recently been proposed as a new class of dark matter structure. These minihalos would be produced by phase transitions in the early Universe or features in the inflaton potential, and constitute non-baryonic massive compact halo objects (MACHOs) today. We examine the prospect of detecting ultracompact minihalos in gamma-rays if dark matter consists of self-annihilating particles. We compute present-day fluxes from minihalos produced in the electron-positron annihilation epoch, and the QCD and electroweak phase transitions in the early Universe. Even at a distance of 100 pc, minihalos produced during the electron-positron annihilation epoch should be eminently detectable today, either by the Fermi satellite, current Air Cherenkov telescopes, or even in archival EGRET data. Within ~1 pc, minihalos formed in the QCD phase transition would have similar predicted fluxes to the dwarf spheroidal galaxies targeted by current indirect dark matter searches, so might also be detectable by present or upcoming experiments.

preprint2009arXiv

Dark stars at the Galactic centre - the main sequence

In regions of very high dark matter density such as the Galactic centre, the capture and annihilation of WIMP dark matter by stars has the potential to significantly alter their evolution. We describe the dark stellar evolution code DarkStars, and present a series of grids of WIMP-influenced stellar models for main sequence stars. We describe changes in which occur as a function of the rate of energy injection by WIMPs, for stars of 0.3-2.0 solar masses and metallicities Z = 0.0003-0.02. We show what rates of energy injection can be obtained using realistic orbital parameters for stars at the Galactic centre, including detailed consideration of the velocity and density profiles of dark matter. Capture and annihilation rates are strongly boosted when stars follow elliptical rather than circular orbits. If there is a spike of dark matter induced by the supermassive black hole at the Galactic centre, single solar-mass stars following orbits with periods as long as 50 years and eccentricities as low as 0.9 could be significantly affected. Binary systems with similar periods about the Galactic centre could be affected on even less eccentric orbits. The most striking evidence of this scenario would be the existence of a binary consisting of a low-mass protostar and a higher-mass evolved star. The observation of low-mass stars and/or binaries on such orbits would either provide a detection of WIMP dark matter, or place stringent limits on the combination of the WIMP mass, spin-dependent nuclear-scattering cross-section, halo density and velocity distribution near the Galactic centre. In some cases, the limits on the WIMP mass and spin-dependent nuclear-scattering cross-section would be of comparable sensitivity to current direct-detection experiments.

preprint2009arXiv

Direct Constraints on Minimal Supersymmetry from Fermi-LAT Observations of the Dwarf Galaxy Segue 1

The dwarf galaxy Segue 1 is one of the most promising targets for the indirect detection of dark matter. Here we examine what constraints 9 months of Fermi-LAT gamma-ray observations of Segue 1 place upon the Constrained Minimal Supersymmetric Standard Model (CMSSM), with the lightest neutralino as the dark matter particle. We use nested sampling to explore the CMSSM parameter space, simultaneously fitting other relevant constraints from accelerator bounds, the relic density, electroweak precision observables, the anomalous magnetic moment of the muon and B-physics. We include spectral and spatial fits to the Fermi observations, a full treatment of the instrumental response and its related uncertainty, and detailed background models. We also perform an extrapolation to 5 years of observations, assuming no signal is observed from Segue 1 in that time. Results marginally disfavour models with low neutralino masses and high annihilation cross-sections. Virtually all of these models are however already disfavoured by existing experimental or relic density constraints.

preprint2009arXiv

The chemical composition of the Sun

The solar chemical composition is an important ingredient in our understanding of the formation, structure and evolution of both the Sun and our solar system. Furthermore, it is an essential reference standard against which the elemental contents of other astronomical objects are compared. In this review we evaluate the current understanding of the solar photospheric composition. In particular, we present a re-determination of the abundances of nearly all available elements, using a realistic new 3-dimensional (3D), time-dependent hydrodynamical model of the solar atmosphere. We have carefully considered the atomic input data and selection of spectral lines, and accounted for departures from LTE whenever possible. The end result is a comprehensive and homogeneous compilation of the solar elemental abundances. Particularly noteworthy findings are significantly lower abundances of carbon, nitrogen, oxygen and neon compared with the widely-used values of a decade ago. The new solar chemical composition is supported by a high degree of internal consistency between available abundance indicators, and by agreement with values obtained in the solar neighborhood and from the most pristine meteorites. There is, however, a stark conflict with standard models of the solar interior according to helioseismology, a discrepancy that has yet to find a satisfactory resolution.

preprint2009arXiv

The DarkStars code: a publicly available dark stellar evolution package

We announce the public release of the 'dark' stellar evolution code DarkStars. The code simultaneously solves the equations of WIMP capture and annihilation in a star with those of stellar evolution assuming approximate hydrostatic equilibrium. DarkStars includes the most extensive WIMP microphysics of any dark evolution code to date. The code employs detailed treatments of the capture process from a range of WIMP velocity distributions, as well as composite WIMP distribution and conductive energy transport schemes based on the WIMP mean-free path in the star. We give a brief description of the input physics and practical usage of the code, as well as examples of its application to dark stars at the Galactic centre.

preprint2008arXiv

Impacts of WIMP dark matter upon stellar evolution: main-sequence stars

The presence of large amounts of WIMP dark matter in stellar cores has been shown to have significant effects upon models of stellar evolution. We present a series of detailed grids of WIMP-influenced stellar models for main sequence stars, computed using the DarkStars code. We describe the changes in stellar structure and main sequence evolution which occur for masses ranging from 0.3 to 2.0 solar masses and metallicities from Z = 0.0003-0.02, as a function of the rate of energy injection by WIMPs. We then go on to show what rates of energy injection can be obtained using realistic orbital parameters for stars near supermassive black holes, including detailed considerations of dark matter halo velocity and density profiles. Capture and annihilation rates are strongly boosted when stars follow elliptical rather than circular orbits, causing WIMP annihilation to provide up to 100 times the energy of hydrogen fusion in stars at the Galactic centre.

Pat Scott

What is connected

Connect this record

See the researcher in context

Building this map preview

63 published item(s)

Global fit of 2HDM with future collider results

Simple and statistically sound recommendations for analysing physical theories

Simulation of energy transport by dark matter scattering in stars

CosmoBit: A GAMBIT module for computing cosmological observables and likelihoods

The GAMBIT Universal Model Machine: from Lagrangians to Likelihoods

Thermal WIMPs and the Scale of New Physics: Global Fits of Dirac Dark Matter Effective Field Theories

GAMBIT and its Application in the Search for Physics Beyond the Standard Model

Global analyses of Higgs portal singlet dark matter models using GAMBIT

Reinterpretation of LHC Results for New Physics: Status and Recommendations after Run 2

Self-consistent approach for measuring the energy spectra and composition of cosmic rays and determining the properties of hadronic interactions at high energy

Axion global fits with Peccei-Quinn symmetry breaking before inflation using GAMBIT

ColliderBit: a GAMBIT module for the calculation of high-energy collider observables and likelihoods

Combined collider constraints on neutralinos and charginos

SpecBit, DecayBit and PrecisionBit: GAMBIT modules for computing mass spectra, particle decay rates and precision observables

A global fit of the MSSM with GAMBIT

DarkBit: A GAMBIT module for computing dark matter observables and likelihoods

Global analyses of supersymmetry with GAMBIT

Global fits of GUT-scale SUSY models with GAMBIT

Impact of vacuum stability, perturbativity and XENON1T on global fits of $\mathbb{Z}_2$ and $\mathbb{Z}_3$ scalar singlet dark matter

Status of the scalar singlet dark matter model

Comparison of statistical sampling methods with ScannerBit, the GAMBIT scanning module

FlavBit: A GAMBIT module for computing flavour observables and likelihoods

GAMBIT: The Global and Modular Beyond-the-Standard-Model Inference Tool

A realistic assessment of the CTA sensitivity to dark matter annihilation

Combined analysis of effective Higgs portal dark matter models

Generalised form factor dark matter in the Sun

Implications of solar wind measurements for solar models and composition

Pippi - painless parsing, post-processing and plotting of posterior and likelihood samples

Ultracompact minihalos as probes of inflationary cosmology

Updated constraints on velocity and momentum-dependent asymmetric dark matter

A possible indication of momentum-dependent asymmetric dark matter in the Sun

Constraints on cosmic strings from ultracompact minihalos

Update on scalar singlet dark matter

Sensitivity of CTA to dark matter signals from the Galactic Center

The elemental composition of the Sun I. The intermediate mass elements Na to Ca

The elemental composition of the Sun II. The iron group elements Sc to Ni

The elemental composition of the Sun III. The heavy elements Cu to Th

Thermal conduction by dark matter with velocity and momentum-dependent cross-sections

Dark Matter CMB Constraints and Likelihoods for Poor Particle Physicists

Detection of isolated population III stars with the James Webb Space Telescope

Hunting for dark halo substructure using submilliarcsecond-scale observations of macrolensed radio jets

Improved constraints on the primordial power spectrum at small scales from ultracompact minihalos

Light bosons and photospheric solutions to the solar abundance problem

Number Counts and Non-Gaussianity

Sensitivity of IceCube-DeepCore to neutralino dark matter in the MSSM-25

Constraints on small-scale cosmological perturbations from gamma-ray searches for dark matter

Fundamental statistical limitations of future dark matter direct detection experiments

Quantifying impacts of short-term plasticity on neuronal information transfer

Dark stars: structure, evolution and impacts upon the high-redshift Universe

How well will ton-scale dark matter direct detection experiments constrain minimal supersymmetry?

Impacts of Dark Stars on Reionization and Signatures in the Cosmic Microwave Background

Implications for constrained supersymmetry of combined H.E.S.S. observations of dwarf galaxies, the Galactic halo and the Galactic Centre

Observational constraints on supermassive dark stars

Searches for Particle Dark Matter: An Introduction

Statistical coverage for supersymmetric parameter estimation: a case study with direct detection of dark matter

A Profile Likelihood Analysis of the Constrained MSSM with Genetic Algorithms

Finding high-redshift dark stars with the James Webb Space Telescope

Gamma rays from ultracompact primordial dark matter minihalos

Dark stars at the Galactic centre - the main sequence

Direct Constraints on Minimal Supersymmetry from Fermi-LAT Observations of the Dwarf Galaxy Segue 1

The chemical composition of the Sun

The DarkStars code: a publicly available dark stellar evolution package

Impacts of WIMP dark matter upon stellar evolution: main-sequence stars