|
|
Strege, C., Bertone, G., Besjes, G. J., Caron, S., Ruiz de Austri, R., Strubig, A., et al. (2014). Profile likelihood maps of a 15-dimensional MSSM. J. High Energy Phys., 09(9), 081–59pp.
Abstract: We present statistically convergent profile likelihood maps obtained via global fits of a phenomenological Minimal Supersymmetric Standard Model with 15 free parameters (the MSSM-15), based on over 250M points. We derive constraints on the model parameters from direct detection limits on dark matter, the Planck relic density measurement and data from accelerator searches. We provide a detailed analysis of the rich phenomenology of this model, and determine the SUSY mass spectrum and dark matter properties that are preferred by current experimental constraints. We evaluate the impact of the measurement of the anomalous magnetic moment of the muon (g – 2) on our results, and provide an analysis of scenarios in which the lightest neutralino is a subdominant component of the dark matter. The MSSM-15 parameters are relatively weakly constrained by current data sets, with the exception of the parameters related to dark matter phenomenology (M-1, M-2, mu), which are restricted to the sub-TeV regime, mainly due to the relic density constraint. The mass of the lightest neutralino is found to be < 1.5TeV at 99% C.L., but can extend up to 3 TeV when excluding the g – 2 constraint from the analysis. Low-mass bino-like neutralinos are strongly favoured, with spin-independent scattering cross-sections extending to very small values, similar to 10(-20) pb. ATLAS SUSY null searches strongly impact on this mass range, and thus rule out a region of parameter space that is outside the reach of any current or future direct detection experiment. The best-fit point obtained after inclusion of all data corresponds to a squark mass of 2.3 TeV, a gluino mass of 2.1 TeV and a 130 GeV neutralino with a spin-independent cross-section of 2.4 x 10(-10) pb, which is within the reach of future multi-ton scale direct detection experiments and of the upcoming LHC run at increased centre-of-mass energy.
|
|
|
|
Ghosh, P., Lopez-Fogliani, D. E., Mitsou, V. A., Muñoz, C., & Ruiz de Austri, R. (2014). Probing the μnu SSM with light scalars, pseudoscalars and neutralinos from the decay of a SM-like Higgs boson at the LHC. J. High Energy Phys., 11(11), 102–57pp.
Abstract: The “mu from nu” supersymmetric standard model (mu nu SSM) can accommodate the newly discovered Higgs-like scalar boson with a mass around 125GeV. This model provides a solution to the mu-problem and simultaneously reproduces correct neutrino physics by the simple use of right-handed neutrino superfields. These new superfields together with the introduced R-parity violation can produce novel and characteristic signatures of the μnu SSM at the LHC. We explore the signatures produced through two-body Higgs decays into the new states, provided that these states lie below in the mass spectrum. For example, a pair produced light neutralinos depending on the associated decay length can give rise to displaced multi-leptons/taus/jets/photons with small/moderate missing transverse energy. In the same spirit, a Higgs-like scalar decaying to a pair of scalars/pseudoscalars can produce final states with prompt multi-leptons/taus/jets/photons.
|
|
|
|
Trotta, R., Johannesson, G., Moskalenko, I. V., Porter, T. A., Ruiz de Austri, R., & Strong, A. W. (2011). Constraints on Cosmic-Ray Propagation Models from a Global Bayesian Analysis. Astrophys. J., 729(2), 106–16pp.
Abstract: Research in many areas of modern physics such as, e. g., indirect searches for dark matter and particle acceleration in supernova remnant shocks rely heavily on studies of cosmic rays (CRs) and associated diffuse emissions (radio, microwave, X-rays, gamma-rays). While very detailed numerical models of CR propagation exist, a quantitative statistical analysis of such models has been so far hampered by the large computational effort that those models require. Although statistical analyses have been carried out before using semi-analytical models (where the computation is much faster), the evaluation of the results obtained from such models is difficult, as they necessarily suffer from many simplifying assumptions. The main objective of this paper is to present a working method for a full Bayesian parameter estimation for a numerical CR propagation model. For this study, we use the GALPROP code, the most advanced of its kind, which uses astrophysical information, and nuclear and particle data as inputs to self-consistently predict CRs, gamma-rays, synchrotron, and other observables. We demonstrate that a full Bayesian analysis is possible using nested sampling and Markov Chain Monte Carlo methods (implemented in the SuperBayeS code) despite the heavy computational demands of a numerical propagation code. The best-fit values of parameters found in this analysis are in agreement with previous, significantly simpler, studies also based on GALPROP.
|
|
|
|
Balazs, C. et al, Mamuzic, J., & Ruiz de Austri, R. (2021). A comparison of optimisation algorithms for high-dimensional particle and astrophysics applications. J. High Energy Phys., 05(5), 108–46pp.
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.
|
|
|
|
Stoppa, F., Ruiz de Austri, R., Vreeswijk, P., Bhattacharyya, S., Caron, S., Bloemen, S., et al. (2023). AutoSourceID-FeatureExtractor Optical image analysis using a two-step mean variance estimation network for feature estimation and uncertainty characterisation. Astron. Astrophys., 680, A108–14pp.
Abstract: Aims. In astronomy, machine learning has been successful in various tasks such as source localisation, classification, anomaly detection, and segmentation. However, feature regression remains an area with room for improvement. We aim to design a network that can accurately estimate sources' features and their uncertainties from single-band image cutouts, given the approximated locations of the sources provided by the previously developed code AutoSourceID-Light (ASID-L) or other external catalogues. This work serves as a proof of concept, showing the potential of machine learning in estimating astronomical features when trained on meticulously crafted synthetic images and subsequently applied to real astronomical data.Methods. The algorithm presented here, AutoSourceID-FeatureExtractor (ASID-FE), uses single-band cutouts of 32x32 pixels around the localised sources to estimate flux, sub-pixel centre coordinates, and their uncertainties. ASID-FE employs a two-step mean variance estimation (TS-MVE) approach to first estimate the features and then their uncertainties without the need for additional information, for example the point spread function (PSF). For this proof of concept, we generated a synthetic dataset comprising only point sources directly derived from real images, ensuring a controlled yet authentic testing environment.Results. We show that ASID-FE, trained on synthetic images derived from the MeerLICHT telescope, can predict more accurate features with respect to similar codes such as SourceExtractor and that the two-step method can estimate well-calibrated uncertainties that are better behaved compared to similar methods that use deep ensembles of simple MVE networks. Finally, we evaluate the model on real images from the MeerLICHT telescope and the Zwicky Transient Facility (ZTF) to test its transfer learning abilities.
|
|
|
|
Stoppa, F., Vreeswijk, P., Bloemen, S., Bhattacharyya, S., Caron, S., Johannesson, G., et al. (2022). AutoSourceID-Light Fast optical source localization via U-Net and Laplacian of Gaussian. Astron. Astrophys., 662, A109–8pp.
Abstract: Aims. With the ever-increasing survey speed of optical wide-field telescopes and the importance of discovering transients when they are still young, rapid and reliable source localization is paramount. We present AutoSourceID-Light (ASID-L), an innovative framework that uses computer vision techniques that can naturally deal with large amounts of data and rapidly localize sources in optical images. Methods. We show that the ASID-L algorithm based on U-shaped networks and enhanced with a Laplacian of Gaussian filter provides outstanding performance in the localization of sources. A U-Net network discerns the sources in the images from many different artifacts and passes the result to a Laplacian of Gaussian filter that then estimates the exact location. Results. Using ASID-L on the optical images of the MeerLICHT telescope demonstrates the great speed and localization power of the method. We compare the results with SExtractor and show that our method outperforms this more widely used method. ASID-L rapidly detects more sources not only in low- and mid-density fields, but particularly in areas with more than 150 sources per square arcminute. The training set and code used in this paper are publicly available.
|
|
|
|
Stoppa, F., Bhattacharyya, S., Ruiz de Austri, R., Vreeswijk, P., Caron, S., Zaharijas, G., et al. (2023). AutoSourceID-Classifier Star-galaxy classification using a convolutional neural network with spatial information. Astron. Astrophys., 680, A109–16pp.
Abstract: Aims. Traditional star-galaxy classification techniques often rely on feature estimation from catalogs, a process susceptible to introducing inaccuracies, thereby potentially jeopardizing the classification's reliability. Certain galaxies, especially those not manifesting as extended sources, can be misclassified when their shape parameters and flux solely drive the inference. We aim to create a robust and accurate classification network for identifying stars and galaxies directly from astronomical images.Methods. The AutoSourceID-Classifier (ASID-C) algorithm developed for this work uses 32x32 pixel single filter band source cutouts generated by the previously developed AutoSourceID-Light (ASID-L) code. By leveraging convolutional neural networks (CNN) and additional information about the source position within the full-field image, ASID-C aims to accurately classify all stars and galaxies within a survey. Subsequently, we employed a modified Platt scaling calibration for the output of the CNN, ensuring that the derived probabilities were effectively calibrated, delivering precise and reliable results.Results. We show that ASID-C, trained on MeerLICHT telescope images and using the Dark Energy Camera Legacy Survey (DECaLS) morphological classification, is a robust classifier and outperforms similar codes such as SourceExtractor. To facilitate a rigorous comparison, we also trained an eXtreme Gradient Boosting (XGBoost) model on tabular features extracted by SourceExtractor. While this XGBoost model approaches ASID-C in performance metrics, it does not offer the computational efficiency and reduced error propagation inherent in ASID-C's direct image-based classification approach. ASID-C excels in low signal-to-noise ratio and crowded scenarios, potentially aiding in transient host identification and advancing deep-sky astronomy.
|
|
|
|
Cabrera, M. E., Casas, A., Ruiz de Austri, R., & Bertone, G. (2014). LHC and dark matter phenomenology of the NUGHM. J. High Energy Phys., 12(12), 114–39pp.
Abstract: We present a Bayesian analysis of the NUGHM, a supersymmetric scenario with non-universal gaugino masses and Higgs masses, including all the relevant experimental observables and dark matter constraints. The main merit of the NUGHM is that it essentially includes all the possibilities for dark matter (DM) candidates within the MSSM, since the neutralino and chargino spectrum -and composition- are as free as they can be in the general MSSM. We identify the most probable regions in the NUHGM parameter space, and study the associated phenomenology at the LHC and the prospects for DM direct detection. Requiring that the neutralino makes all of the DM in the Universe, we identify two preferred regions around m(chi 10) = 1 TeV, 3 TeV, which correspond to the (almost) pure Higgsino and wino case. There exist other marginal regions (e.g. Higgs-funnel), but with much less statistical weight. The prospects for detection at the LHC in this case are quite pessimistic, but future direct detection experiments like LUX and XENON1T, will be able to probe this scenario. In contrast, when allowing other DM components, the prospects for detection at the LHC become more encouraging – the most promising signals being, beside the production of gluinos and squarks, the production of the heavier chargino and neutralino states, which lead to WZ and same-sign WW final states – and direct detection remains a complementary, and even more powerful, way to probe the scenario.
|
|
|
|
Jueid, A., Kip, J., Ruiz de Austri, R., & Skands, P. (2024). The Strong Force meets the Dark Sector: a robust estimate of QCD uncertainties for anti-matter dark matter searches. J. High Energy Phys., 02(2), 119–48pp.
Abstract: In dark-matter annihilation channels to hadronic final states, stable particles – such as positrons, photons, antiprotons, and antineutrinos – are produced via complex sequences of phenomena including QED/QCD radiation, hadronisation, and hadron decays. These processes are normally modelled by Monte Carlo (MC) event generators whose limited accuracy imply intrinsic QCD uncertainties on the predictions for indirect-detection experiments like Fermi-LAT, Pamela, IceCube or Ams-02. In this article, we perform a comprehensive analysis of QCD uncertainties, meaning both perturbative and nonperturbative sources of uncertainty are included – estimated via variations of MC renormalization-scale and fragmentation-function parameters, respectively – in antimatter spectra from dark-matter annihilation, based on parametric variations of the Pythia 8 event generator. After performing several retunings of light-quark fragmentation functions, we define a set of variations that span a conservative estimate of the QCD uncertainties. We estimate the effects on antimatter spectra for various annihilation channels and final-state particle species, and discuss their impact on fitted values for the dark-matter mass and thermally-averaged annihilation cross section. We find dramatic impacts which can go up to O(10%) for the annihilation cross section. We provide the spectra in tabulated form including QCD uncertainties and code snippets to perform fast dark-matter fits, in this github repository.
|
|
|
|
Caron, S., Casas, J. A., Quilis, J., & Ruiz de Austri, R. (2018). Anomaly-free dark matter with harmless direct detection constraints. J. High Energy Phys., 12(12), 126–24pp.
Abstract: Dark matter (DM) interacting with the SM fields via a Z-boson (Z-portal') remains one of the most attractive WIMP scenarios, both from the theoretical and the phenomenological points of view. In order to avoid the strong constraints from direct detection and dilepton production, it is highly convenient that the Z has axial coupling to DM and leptophobic couplings to the SM particles, respectively. The latter implies that the associated U(1) coincides with baryon number in the SM sector. In this paper we completely classify the possible anomaly-free leptophobic Z with minimal dark sector, including the cases where the coupling to DM is axial. The resulting scenario is very predictive and perfectly viable from the present constraints from DM detection, EW observables and LHC data (di-lepton, di-jet and mono-jet production). We analyze all these constraints, obtaining the allowed areas in the parameter space, which generically prefer mZ less than or similar to 500 GeV, apart from resonant regions. The best chances to test these viable areas come from future LHC measurements.
|
|
