Source author record

Sebastian Wild

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

23 published item(s)

preprint2023arXiv

Multiway Powersort

We present a stable mergesort variant, Multiway Powersort, that exploits existing runs and finds nearly-optimal merging orders for k-way merges with negligible overhead. This builds on Powersort (Munro & Wild, ESA2018), which has recently replaced Timsort's suboptimal merge policy in the CPython reference implementation of Python, as well as in PyPy and further libraries. Multiway Powersort reduces the number of memory transfers, which increasingly determine the cost of internal sorting (as observed with Multiway Quicksort (Kushagra et al., ALENEX 2014; Aumüller & Dietzfelbinger, TALG 2016; Wild, PhD thesis 2016) and the inclusion of Dual-Pivot Quicksort in the Java runtime library). We demonstrate that our 4-way Powersort implementation can achieve substantial speedups over standard (2-way) Powersort and other stable sorting methods without compromising the optimally run-adaptive performance of Powersort.

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

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.

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

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$.

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

DAMA confronts null searches in the effective theory of dark matter-nucleon interactions

We examine the dark matter interpretation of the modulation signal reported by the DAMA experiment from the perspective of effective field theories displaying Galilean invariance. We consider the most general effective coupling leading to the elastic scattering of a dark matter particle with spin 0 or 1/2 off a nucleon, and we analyze the compatibility of the DAMA signal with the null results from other direct detection experiments, as well as with the non-observation of a high energy neutrino flux in the direction of the Sun from dark matter annihilation. To this end, we develop a novel semi-analytical approach for comparing experimental results in the high-dimensional parameter space of the non-relativistic effective theory. Assuming the standard halo model, we find a strong tension between the dark matter interpretation of the DAMA modulation signal and the null result experiments. We also list possible ways-out of this conclusion.

preprint2016arXiv

Signatures from Scalar Dark Matter with a Vector-like Quark Mediator

We present a comprehensive study of a model where the dark matter is composed of a singlet real scalar that couples to the Standard Model predominantly via a Yukawa interaction with a light quark and a colored vector-like fermion. A distinctive feature of this scenario is that thermal freeze-out in the early universe may be driven by annihilation both into gluon pairs at one-loop ($gg$) and by virtual internal Bremsstrahlung of a gluon ($q \bar{q} g$). Such a dark matter candidate may also be tested through direct and indirect detection and at the LHC; viable candidates have either a mass nearly degenerate with that of the fermionic mediator or a mass above about 2 TeV.

preprint2015arXiv

A novel approach to derive halo-independent limits on dark matter properties

We propose a method that allows to place an upper limit on the dark matter elastic scattering cross section with nucleons which is independent of the velocity distribution. Our approach combines null results from direct detection experiments with indirect searches at neutrino telescopes, and goes beyond previous attempts to remove astrophysical uncertainties in that it directly constrains the particle physics properties of the dark matter. The resulting halo-independent upper limits on the scattering cross section of dark matter are remarkably strong and reach $σ_{\text{SI}}^{p} \lesssim 10^{-43} \, (10^{-42})~\text{cm}^2$ and $σ_{\text{SD}}^{p} \lesssim 10^{-37} \, (3\times 10^{-37 })~\text{cm}^2$, for dark matter particles of $m_{\text{DM}}\sim 1~\text{TeV}$ annihilating into $W^+W^-$ ($b\bar b$), assuming $ρ_\text{loc}=0.3\text{ GeV}/\text{cm}^3$.

preprint2015arXiv

Analysis of Pivot Sampling in Dual-Pivot Quicksort

The new dual-pivot Quicksort by Vladimir Yaroslavskiy - used in Oracle's Java runtime library since version 7 - features intriguing asymmetries. They make a basic variant of this algorithm use less comparisons than classic single-pivot Quicksort. In this paper, we extend the analysis to the case where the two pivots are chosen as fixed order statistics of a random sample. Surprisingly, dual-pivot Quicksort then needs more comparisons than a corresponding version of classic Quicksort, so it is clear that counting comparisons is not sufficient to explain the running time advantages observed for Yaroslavskiy's algorithm in practice. Consequently, we take a more holistic approach and give also the precise leading term of the average number of swaps, the number of executed Java Bytecode instructions and the number of scanned elements, a new simple cost measure that approximates I/O costs in the memory hierarchy. We determine optimal order statistics for each of the cost measures. It turns out that the asymmetries in Yaroslavskiy's algorithm render pivots with a systematic skew more efficient than the symmetric choice. Moreover, we finally have a convincing explanation for the success of Yaroslavskiy's algorithm in practice: Compared with corresponding versions of classic single-pivot Quicksort, dual-pivot Quicksort needs significantly less I/Os, both with and without pivot sampling.

preprint2015arXiv

Average Case and Distributional Analysis of Dual-Pivot Quicksort

In 2009, Oracle replaced the long-serving sorting algorithm in its Java 7 runtime library by a new dual-pivot Quicksort variant due to Vladimir Yaroslavskiy. The decision was based on the strikingly good performance of Yaroslavskiy's implementation in running time experiments. At that time, no precise investigations of the algorithm were available to explain its superior performance - on the contrary: Previous theoretical studies of other dual-pivot Quicksort variants even discouraged the use of two pivots. Only in 2012, two of the authors gave an average case analysis of a simplified version of Yaroslavskiy's algorithm, proving that savings in the number of comparisons are possible. However, Yaroslavskiy's algorithm needs more swaps, which renders the analysis inconclusive. To force the issue, we herein extend our analysis to the fully detailed style of Knuth: We determine the exact number of executed Java Bytecode instructions. Surprisingly, Yaroslavskiy's algorithm needs sightly more Bytecode instructions than a simple implementation of classic Quicksort - contradicting observed running times. Like in Oracle's library implementation we incorporate the use of Insertionsort on small subproblems and show that it indeed speeds up Yaroslavskiy's Quicksort in terms of Bytecodes; but even with optimal Insertionsort thresholds the new Quicksort variant needs slightly more Bytecode instructions on average. Finally, we show that the (suitably normalized) costs of Yaroslavskiy's algorithm converge to a random variable whose distribution is characterized by a fixed-point equation. From that, we compute variances of costs and show that for large n, costs are concentrated around their mean.

preprint2015arXiv

Dirac dark matter with a charged mediator: a comprehensive one-loop analysis of the direct detection phenomenology

We analyze the direct detection signals of a toy model consisting of a Dirac dark matter particle which couples to one Standard Model fermion via a scalar mediator. For all scenarios, the dark matter particle scatters off nucleons via one loop-induced electromagnetic and electroweak moments, as well as via the one-loop exchange of a Higgs boson. Besides, and depending on the details of the model, the scattering can also be mediated at tree level via the exchange of the scalar mediator or at one loop via gluon-gluon interactions. We show that, for thermally produced dark matter particles, the current limits from the LUX experiment on these scenarios are remarkably strong, even for dark matter coupling only to leptons. We also discuss future prospects for XENON1T and DARWIN and we argue that multi-ton xenon detectors will be able to probe practically the whole parameter space of the model consistent with thermal production and perturbativity. We also discuss briefly the implications of our results for the dark matter interpretation of the Galactic GeV excess.

preprint2014arXiv

Analysis of Branch Misses in Quicksort

The analysis of algorithms mostly relies on counting classic elementary operations like additions, multiplications, comparisons, swaps etc. This approach is often sufficient to quantify an algorithm's efficiency. In some cases, however, features of modern processor architectures like pipelined execution and memory hierarchies have significant impact on running time and need to be taken into account to get a reliable picture. One such example is Quicksort: It has been demonstrated experimentally that under certain conditions on the hardware the classically optimal balanced choice of the pivot as median of a sample gets harmful. The reason lies in mispredicted branches whose rollback costs become dominating. In this paper, we give the first precise analytical investigation of the influence of pipelining and the resulting branch mispredictions on the efficiency of (classic) Quicksort and Yaroslavskiy's dual-pivot Quicksort as implemented in Oracle's Java 7 library. For the latter it is still not fully understood why experiments prove it 10% faster than a highly engineered implementation of a classic single-pivot version. For different branch prediction strategies, we give precise asymptotics for the expected number of branch misses caused by the aforementioned Quicksort variants when their pivots are chosen from a sample of the input. We conclude that the difference in branch misses is too small to explain the superiority of the dual-pivot algorithm.

preprint2014arXiv

Analysis of Quickselect under Yaroslavskiy's Dual-Pivoting Algorithm

There is excitement within the algorithms community about a new partitioning method introduced by Yaroslavskiy. This algorithm renders Quicksort slightly faster than the case when it runs under classic partitioning methods. We show that this improved performance in Quicksort is not sustained in Quickselect; a variant of Quicksort for finding order statistics. We investigate the number of comparisons made by Quickselect to find a key with a randomly selected rank under Yaroslavskiy's algorithm. This grand averaging is a smoothing operator over all individual distributions for specific fixed order statistics. We give the exact grand average. The grand distribution of the number of comparison (when suitably scaled) is given as the fixed-point solution of a distributional equation of a contraction in the Zolotarev metric space. Our investigation shows that Quickselect under older partitioning methods slightly outperforms Quickselect under Yaroslavskiy's algorithm, for an order statistic of a random rank. Similar results are obtained for extremal order statistics, where again we find the exact average, and the distribution for the number of comparisons (when suitably scaled). Both limiting distributions are of perpetuities (a sum of products of independent mixed continuous random variables).

preprint2014arXiv

Antihelium from Dark Matter

Cosmic-ray anti-nuclei provide a promising discovery channel for the indirect detection of particle dark matter. Hadron showers produced by the pair-annihilation or decay of Galactic dark matter generate anti-nucleons which can in turn form light anti-nuclei. Previous studies have only focused on the spectrum and flux of low energy antideuterons which, although very rarely, are occasionally also produced by cosmic-ray spallation. Heavier elements ($A\geq3$) have instead entirely negligible astrophysical background and a primary yield from dark matter which could be detectable by future experiments. Using a Monte Carlo event generator and an event-by-event phase space analysis, we compute, for the first time, the production spectrum of \antihe and \antiT for dark matter annihilating or decaying to $b\bar{b}$ and ${W^+}{W^-}$ final states. We then employ a semi-analytic model of interstellar and heliospheric propagation to calculate the \antihe flux as well as to provide tools to relate the anti-helium spectrum corresponding to an arbitrary antideuteron spectrum. Finally, we discuss prospects for current and future experiments, including GAPS and AMS-02.

preprint2014arXiv

Higher order dark matter annihilations in the Sun and implications for IceCube

Dark matter particles captured in the Sun would annihilate producing a neutrino flux that could be detected at the Earth. In some channels, however, the neutrino flux lies in the MeV range and is thus undetectable at IceCube, namely when the dark matter particles annihilate into $e^+e^-$, $μ^+μ^-$ or light quarks. On the other hand, the same interaction that mediates the annihilations into light fermions also leads, via higher order effects, to the production of weak gauge bosons (and in the case of quarks also gluons) that generate a high energy neutrino flux potentially observable at IceCube. We consider in this paper tree level annihilations into a fermion-antifermion pair with the associated emission of one gauge boson and one loop annihilations into two gauge bosons, and we calculate the limits on the scattering cross section of dark matter particles with protons in scenarios where the dark matter particle couples to electrons, muons or light quarks from the non-observation of an excess of neutrino events in the direction of the Sun. We find that the limits on the spin-dependent scattering cross section are, for some scenarios, stronger than the limits from direct detection experiments.

preprint2014arXiv

Pivot Sampling in Dual-Pivot Quicksort

The new dual-pivot Quicksort by Vladimir Yaroslavskiy - used in Oracle's Java runtime library since version 7 - features intriguing asymmetries in its behavior. They were shown to cause a basic variant of this algorithm to use less comparisons than classic single-pivot Quicksort implementations. In this paper, we extend the analysis to the case where the two pivots are chosen as fixed order statistics of a random sample and give the precise leading term of the average number of comparisons, swaps and executed Java Bytecode instructions. It turns out that - unlike for classic Quicksort, where it is optimal to choose the pivot as median of the sample - the asymmetries in Yaroslavskiy's algorithm render pivots with a systematic skew more efficient than the symmetric choice. Moreover, the optimal skew heavily depends on the employed cost measure; most strikingly, abstract costs like the number of swaps and comparisons yield a very different result than counting Java Bytecode instructions, which can be assumed most closely related to actual running time.

preprint2014arXiv

Sharp Gamma-ray Spectral Features from Scalar Dark Matter Annihilations

The search for sharp features in the gamma-ray spectrum is a promising approach to identify a signal from dark matter annihilations over the astrophysical backgrounds. In this paper we investigate the generation of gamma-ray lines and internal bremsstrahlung signals in a toy model where the dark matter particle is a real scalar that couples to a lepton and an exotic fermion via a Yukawa coupling. We show that the Fermi-LAT and H.E.S.S. searches for line-like spectral features severely constrain regions of the parameter space where the scalar dark matter is thermally produced. Finally, we also discuss the complementarity of the searches for sharp spectral features with other indirect dark matter searches, as well as with direct and collider searches.

preprint2013arXiv

Average Case Analysis of Java 7's Dual Pivot Quicksort

Recently, a new Quicksort variant due to Yaroslavskiy was chosen as standard sorting method for Oracle's Java 7 runtime library. The decision for the change was based on empirical studies showing that on average, the new algorithm is faster than the formerly used classic Quicksort. Surprisingly, the improvement was achieved by using a dual pivot approach, an idea that was considered not promising by several theoretical studies in the past. In this paper, we identify the reason for this unexpected success. Moreover, we present the first precise average case analysis of the new algorithm showing e.g. that a random permutation of length $n$ is sorted using $1.9n\ln n-2.46n+\mathcal{O}(\ln n)$ key comparisons and $0.6n\ln n+0.08n+\mathcal{O}(\ln n)$ swaps.

preprint2013arXiv

Determination of the Cosmic Antideuteron Flux in a Monte Carlo approach

We investigate in this paper the antideuteron flux produced in high energy collisions of cosmic rays with the interstellar matter. We employ the Monte Carlo generator DPMJET-III together with the coalescence model to simulate, in an event-by-event basis, the antideuteron production in cosmic ray collisions in our Galaxy. Then, we use a diffusion model to calculate the expected flux. We find a secondary antideuteron flux at the Earth which is approximately a factor of two smaller than in previous calculations.

preprint2013arXiv

High-energy neutrino signals from the Sun in dark matter scenarios with internal bremsstrahlung

We investigate the prospects to observe a high energy neutrino signal from dark matter annihilations in the Sun in scenarios where the dark matter is a Majorana fermion that couples to a quark and a colored scalar via a Yukawa coupling. In this minimal scenario, the dark matter capture and annihilation in the Sun can be studied in a single framework. We find that, for small and moderate mass splitting between the dark matter and the colored scalar, the two-to-three annihilation $q \bar q g$ plays a central role in the calculation of the number of captured dark matter particles. On the other hand, the two-to-three annihilation into $q \bar q Z$ gives, despite its small branching fraction, the largest contribution to the neutrino flux at the Earth at the highest energies. We calculate the limits on the model parameters using IceCube observations of the Sun and we discuss their interplay with the requirement of equilibrium of captures and annihilations in the Sun and with the requirement of thermal dark matter production. We also compare the limits from IceCube to the limits from direct detection, antiproton measurements and collider searches.

preprint2013arXiv

Prospects of antideuteron detection from dark matter annihilations or decays at AMS-02 and GAPS

The search for cosmic antideuterons has been proposed as a promising method to indirectly detect dark matter, due to the very small background flux from spallations expected at the energies relevant to experiments. The antideuteron flux from dark matter annihilations or decays is, however, severely constrained by the non-observation of an excess in the antiproton-to-proton fraction measured by PAMELA. In this paper we calculate, for representative dark matter annihilation and decay channels, upper limits on the number of antideuteron events at AMS-02 and GAPS from requiring that the associated antiproton flux is in agreement with the PAMELA data. To this end, we first analyze in detail the formation of antideuterons in the coalescence model using an event-by-event Monte Carlo simulation and using data from various high energy experiments. We find that the resulting coalescence momentum shows a dependence on the underlying process and on the center of mass energy involved. Then, we calculate, using a diffusion model, the flux of antideuterons at the Earth from dark matter annihilations or decays. Our results indicate that, despite the various sources of uncertainty, the observation of an antideuteron flux at AMS-02 or GAPS from dark matter annihilations or decays will be challenging.