|
|
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.
|
|
|
|
Casas, J. A., Moreno, J. M., Rius, N., Ruiz de Austri, R., & Zaldivar, B. (2011). Fair scans of the seesaw. Consequences for predictions on LFV processes. J. High Energy Phys., 03(3), 034–22pp.
Abstract: We give a straightforward procedure to scan the seesaw parameter-space, using the common “R-parametrization”, in a complete way. This includes a very simple rule to incorporate the perturbativity requirement as a condition for the entries of the R-matrix. As a relevant application, we show that the somewhat propagated belief that BR(mu -> e, gamma) in supersymmetric seesaw models depends strongly on the value of theta(13) is an “optical effect” produced by incomplete scans, and does not hold after a careful analytical and numerical study. When the complete scan is done, BR(mu -> e, gamma) gets very insensitive to theta(13). This holds even if the right-handed neutrino masses are kept constant or under control (as is required for succesful leptogenesis). In most cases the values of BR(mu -> e, gamma) are larger than the experimental upper bound. Including (unflavoured) leptogenesis does not introduce any further dependence on theta(13), although decreases the typical value of BR(mu -> e, gamma).
|
|
|
|
Catumba, G., Ramos, A., & Zaldivar, B. (2025). Stochastic automatic differentiation for Monte Carlo processes. Comput. Phys. Commun., 307, 109396–13pp.
Abstract: Monte Carlo methods represent a cornerstone of computer science. They allow sampling high dimensional distribution functions in an efficient way. In this paper we consider the extension of Automatic Differentiation (AD) techniques to Monte Carlo processes, addressing the problem of obtaining derivatives (and in general, the Taylor series) of expectation values. Borrowing ideas from the lattice field theory community, we examine two approaches. One is based on reweighting while the other represents an extension of the Hamiltonian approach typically used by the Hybrid Monte Carlo (HMC) and similar algorithms. We show that the Hamiltonian approach can be understood as a change of variables of the reweighting approach, resulting in much reduced variances of the coefficients of the Taylor series. This work opens the door to finding other variance reduction techniques for derivatives of expectation values.
|
|
|
|
de los Rios, M., Petac, M., Zaldivar, B., Bonaventura, N. R., Calore, F., & Iocco, F. (2023). Determining the dark matter distribution in simulated galaxies with deep learning. Mon. Not. Roy. Astron. Soc., 525(4), 6015–6035.
Abstract: We present a novel method of inferring the dark matter (DM) content and spatial distribution within galaxies, using convolutional neural networks (CNNs) trained within state-of-the-art hydrodynamical simulations (Illustris-TNG100). Within the controlled environment of the simulation, the framework we have developed is capable of inferring the DM mass distribution within galaxies of mass similar to 10(11)-10(13)M(circle dot) from the gravitationally baryon-dominated internal regions to the DM-rich, baryon-depleted outskirts of the galaxies, with a mean absolute error always below approximate to 0.25 when using photometrical and spectroscopic information. With respect to traditional methods, the one presented here also possesses the advantages of not relying on a pre-assigned shape for the DM distribution, to be applicable to galaxies not necessarily in isolation, and to perform very well even in the absence of spectroscopic observations.
|
|
|
|
Dimitriou, A., Figueroa, D. G., & Zaldivar, B. (2024). Fast likelihood-free reconstruction of gravitational wave backgrounds. J. Cosmol. Astropart. Phys., 09(9), 032–51pp.
Abstract: based) techniques for reconstructing the spectral shape of a gravitational wave background (GWB). We focus on the reconstruction of an arbitrarily shaped signal (approximated by a piecewise power-law in many frequency bins) by the LISA detector, but the method can be easily extended to either template-dependent signals, or to other detectors, as long as a characterisation of the instrumental noise is available. As proof of the technique, we quantify the ability of LISA to reconstruct signals of arbitrary spectral shape (blind reconstruction), considering a diversity of frequency profiles, and including astrophysical backgrounds in some cases. As a teaser of how the method can reconstruct signals characterised by a parameter-dependent template (template reconstruction), we present a dedicated study for power-law signals. While our technique has several advantages with respect to traditional MCMC methods, we validate it with the latter for concrete cases. This work opens the door for both fast and accurate Bayesian parameter estimation of GWBs, with essentially no computational overhead during the inference step. Our set of tools are integrated into the package GWBackFinder, which is publicly available in GitHub.
|
|
|
|
Gammaldi, V., Zaldivar, B., Sanchez-Conde, M. A., & Coronado-Blazquez, J. (2023). A search for dark matter among Fermi-LAT unidentified sources with systematic features in machine learning. Mon. Not. Roy. Astron. Soc., 520(1), 1348–1361.
Abstract: Around one-third of the point-like sources in the Fermi-LAT catalogues remain as unidentified sources (unIDs) today. Indeed, these unIDs lack a clear, univocal association with a known astrophysical source. If dark matter (DM) is composed of weakly interacting massive particles (WIMPs), there is the exciting possibility that some of these unIDs may actually be DM sources, emitting gamma-rays from WIMPs annihilation. We propose a new approach to solve the standard, machine learning (ML) binary classification problem of disentangling prospective DM sources (simulated data) from astrophysical sources (observed data) among the unIDs of the 4FGL Fermi-LAT catalogue. We artificially build two systematic features for the DM data which are originally inherent to observed data: the detection significance and the uncertainty on the spectral curvature. We do it by sampling from the observed population of unIDs, assuming that the DM distributions would, if any, follow the latter. We consider different ML models: Logistic Regression, Neural Network (NN), Naive Bayes, and Gaussian Process, out of which the best, in terms of classification accuracy, is the NN, achieving around 93 . 3 per cent +/- 0 . 7 per cent performance. Other ML evaluation parameters, such as the True Ne gativ e and True Positive rates, are discussed in our work. Applying the NN to the unIDs sample, we find that the de generac y between some astrophysical and DM sources can be partially solved within this methodology. None the less, we conclude that there are no DM source candidates among the pool of 4FGL Fermi-LAT unIDs.
|
|
