|
AbdusSalam, S. S. et al, & Eberhardt, O. (2022). Simple and statistically sound recommendations for analysing physical theories. Rep. Prog. Phys., 85(5), 052201–11pp.
Abstract: 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 Zenodo.
|
|
Adolf, P., Hirsch, M., Krieg, S., Pas, H., & Tabet, M. (2024). Fitting the DESI BAO data with dark energy driven by the Cohen-Kaplan-Nelson bound. J. Cosmol. Astropart. Phys., 08(8), 048–18pp.
Abstract: Gravity constrains the range of validity of quantum field theory. As has been pointed out by Cohen, Kaplan, and Nelson (CKN), such effects lead to interdependent ultraviolet (UV) and infrared (IR) cutoffs that may stabilize the dark energy of the universe against quantum corrections, if the IR cutoff is set by the Hubble horizon. As a consequence of the cosmic expansion, this argument implies a time-dependent dark energy density. In this paper we confront this idea with recent data from DESI BAO, Hubble and supernova measurements. We find that the CKN model provides a better fit to the data than the Lambda CDM model and can compete with other models of time-dependent dark energy that have been studied so far.
|
|
Amerio, A., Calore, F., Serpico, P. D., & Zaldivar, B. (2024). Deepening gamma-ray point-source catalogues with sub-threshold information. J. Cosmol. Astropart. Phys., 03(3), 055–18pp.
Abstract: We propose a novel statistical method to extend Fermi-LAT catalogues of highlatitude -y-ray sources below their nominal threshold. To do so, we rely on the determination of the differential source -count distribution of sub -threshold sources which only provides the statistical flux distribution of faint sources. By simulating ensembles of synthetic skies, we assess quantitatively the likelihood for pixels in the sky with relatively low -test statistics to be due to sources, therefore complementing the source -count distribution with spatial information. Besides being useful to orient efforts towards multi -messenger and multi -wavelength identification of new -y-ray sources, we expect the results to be especially advantageous for statistical applications such as cross -correlation analyses.
|
|
Jueid, A., Kip, J., Ruiz de Austri, R., & Skands, P. (2023). Impact of QCD uncertainties on antiproton spectra from dark-matter annihilation. J. Cosmol. Astropart. Phys., 04(4), 068–15pp.
Abstract: Dark-matter particles that annihilate or decay can undergo complex sequences of processes, including strong and electromagnetic radiation, hadronisation, and hadron de-cays, before particles that are stable on astrophysical time scales are produced. Antiprotons produced in this way may leave footprints in experiments such as AMS-02. Several groups have reported an excess of events in the antiproton flux in the rigidity range of 10-20 GV. However, the theoretical modeling of baryon production is not straightforward and relies in part on phenomenological models in Monte Carlo event generators. In this work, we assess the impact of QCD uncertainties on the spectra of antiprotons from dark-matter annihila-tion. As a proof-of-principle, we show that for a two-parameter model that depends only on the thermally-averaged annihilation cross section ((o -v)) and the dark-matter mass (Mx), QCD uncertainties can affect the best-fit mass by up to ti 14% (with large uncertainties for large DM masses), depending on the choice of Mx and the annihilation channel (bb over bar or W+W-), and (o -v) by up to ti 10%. For comparison, changes to the underlying diffusion parameters are found to be within 1%-5%, and the results are also quite resilient to the choice of cosmic-ray propagation model. These findings indicate that QCD uncertainties need to be included in future DM analyses. To facilitate full-fledged analyses, we provide the spectra in tabulated form including QCD uncertainties and code snippets to perform mass interpolations and quick DM fits. The code can be found in this GitHub [1] repository.
|
|
Reid, B. A. et al, & de Putter, R. (2012). The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: measurements of the growth of structure and expansion rate at z=0.57 from anisotropic clustering. Mon. Not. Roy. Astron. Soc., 426(4), 2719–2737.
Abstract: We analyse the anisotropic clustering of massive galaxies from the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) sample, which consists of 264-283 galaxies in the redshift range 0.43 < z < 0.7 spanning 3275 deg(2). Both peculiar velocities and errors in the assumed redshiftdistance relation (AlcockPaczynski effect) generate correlations between clustering amplitude and orientation with respect to the line of sight. Together with the sharp baryon acoustic oscillation (BAO) standard ruler, our measurements of the broad-band shape of the monopole and quadrupole correlation functions simultaneously constrain the comoving angular diameter distance (2190 +/- 61 Mpc) to z = 0.57, the Hubble expansion rate at z = 0.57 (92.4 +/- 4.5 km s(-1) Mpc(-1)) and the growth rate of structure at that same redshift (d(sigma 8)/d ln a = 0.43 +/- 0.069). Our analysis provides the best current direct determination of both DA and H in galaxy clustering data using this technique. If we further assume a cold dark matter expansion history, our growth constraint tightens to d(sigma 8)/d ln a = 0.415 +/- 0.034. In combination with the cosmic microwave background, our measurements of D-A,H and d(sigma 8)/d ln a all separately require dark energy at z > 0.57, and when combined imply Omega(A) = 0.74 +/- 0.016, independent of the Universe's evolution at z < 0.57. All of these constraints assume scale-independent linear growth, and assume general relativity to compute both O(10 per cent) non-linear model corrections and our errors. In our companion paper, Samushia et al., we explore further cosmological implications of these observations.
|
|
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.
|
|
Villaescusa-Navarro, F. et al, & Villanueva-Domingo, P. (2023). The CAMELS Project: Public Data Release. Astrophys. J. Suppl. Ser., 265(2), 54–14pp.
Abstract: The Cosmology and Astrophysics with Machine Learning Simulations (CAMELS) project was developed to combine cosmology with astrophysics through thousands of cosmological hydrodynamic simulations and machine learning. CAMELS contains 4233 cosmological simulations, 2049 N-body simulations, and 2184 state-of-the-art hydrodynamic simulations that sample a vast volume in parameter space. In this paper, we present the CAMELS public data release, describing the characteristics of the CAMELS simulations and a variety of data products generated from them, including halo, subhalo, galaxy, and void catalogs, power spectra, bispectra, Lya spectra, probability distribution functions, halo radial profiles, and X-rays photon lists. We also release over 1000 catalogs that contain billions of galaxies from CAMELS-SAM: a large collection of N-body simulations that have been combined with the Santa Cruz semianalytic model. We release all the data, comprising more than 350 terabytes and containing 143,922 snapshots, millions of halos, galaxies, and summary statistics. We provide further technical details on how to access, download, read, and process the data at .
|