Abstract: Higgs-pair production via gluon fusion is the dominant production mechanism of Higgs-boson pairs at hadron colliders. In this work, we present details of our numerical determination of the full next-to-leading-order (NLO) QCD corrections to the leading top-quark loops. Since gluon fusion is a loop-induced process at leading order, the NLO calculation requires the calculation of massive two-loop diagrams with up to four different mass/energy scales involved. With the current methods, this can only be done numerically, if no approximations are used. We discuss the setup and details of our numerical integration. This will be followed by a phenomenological analysis of the NLO corrections and their impact on the total cross section and the invariant Higgs-pair mass distribution. The last part of our work will be devoted to the determination of the residual theoretical uncertainties with special emphasis on the uncertainties originating from the scheme and scale dependence of the (virtual) top mass. The impact of the trilinear Higgs-coupling variation on the total cross section will be discussed.

Abstract: Additional Higgs bosons appear in many extensions of the Standard Model (SM). While most existing searches for additional Higgs bosons concentrate on final states consisting of SM particles, final states containing beyond the SM (BSM) particles play an important role in many BSM models. In order to facilitate future searches for such final states, we develop a simplified model framework for heavy Higgs boson decays to a massive SM boson as well as one or more invisible particles. Allowing one kind of BSM mediator in each decay chain, we classify the possible decay topologies for each final state, taking into account all different possibilities for the spin of the mediator and the invisible particles. Our comparison of the kinematic distributions for each possible model realization reveals that the distributions corresponding to the different simplified model topologies are only mildly affected by the different spin hypotheses, while there is significant sensitivity for distinguishing between the different decay topologies. As a consequence, we point out that expressing the results of experimental searches in terms of the proposed simplified model topologies will allow one to constrain wide classes of different BSM models. The application of the proposed simplified model framework is explicitly demonstrated for the example of a mono-Higgs search. For each of the simplified models that are proposed in this paper we provide all necessary ingredients for performing Monte-Carlo simulations such that they can readily be applied in experimental analyses.

Abstract: Optimisation problems are ubiquitous in particle and astrophysics, and involve locating the optimum of a complicated function of many parameters that may be computationally expensive to evaluate. We describe a number of global optimisation algorithms that are not yet widely used in particle astrophysics, benchmark them against random sampling and existing techniques, and perform a detailed comparison of their performance on a range of test functions. These include four analytic test functions of varying dimensionality, and a realistic example derived from a recent global fit of weak-scale supersymmetry. Although the best algorithm to use depends on the function being investigated, we are able to present general conclusions about the relative merits of random sampling, Differential Evolution, Particle Swarm Optimisation, the Covariance Matrix Adaptation Evolution Strategy, Bayesian Optimisation, Grey Wolf Optimisation, and the PyGMO Artificial Bee Colony, Gaussian Particle Filter and Adaptive Memory Programming for Global Optimisation algorithms.