|
Villanueva-Domingo, P., Villaescusa-Navarro, F., Genel, S., Angles-Alcazar, D., Hernquist, L., Marinacci, F., et al. (2023). Weighing the Milky Way and Andromeda galaxies with artificial intelligence. Phys. Rev. D, 107(10), 103003–8pp.
Abstract: We present new constraints on the masses of the halos hosting the Milky Way and Andromeda galaxies derived using graph neural networks. Our models, trained on 2,000 state-of-the-art hydrodynamic simulations of the CAMELS project, only make use of the positions, velocities and stellar masses of the galaxies belonging to the halos, and are able to perform likelihood-free inference on halo masses while accounting for both cosmological and astrophysical uncertainties. Our constraints are in agreement with estimates from other traditional methods, within our derived posterior standard deviation.
|
|
|
Villanueva-Domingo, P., Villaescusa-Navarro, F., Angles-Alcazar, D., Genel, S., Marinacci, F., Spergel, D. N., et al. (2022). Inferring Halo Masses with Graph Neural Networks. Astrophys. J., 935(1), 30–15pp.
Abstract: Understanding the halo-galaxy connection is fundamental in order to improve our knowledge on the nature and properties of dark matter. In this work, we build a model that infers the mass of a halo given the positions, velocities, stellar masses, and radii of the galaxies it hosts. In order to capture information from correlations among galaxy properties and their phase space, we use Graph Neural Networks (GNNs), which are designed to work with irregular and sparse data. We train our models on galaxies from more than 2000 state-of-the-art simulations from the Cosmology and Astrophysics with MachinE Learning Simulations project. Our model, which accounts for cosmological and astrophysical uncertainties, is able to constrain the masses of the halos with a similar to 0.2 dex accuracy. Furthermore, a GNN trained on a suite of simulations is able to preserve part of its accuracy when tested on simulations run with a different code that utilizes a distinct subgrid physics model, showing the robustness of our method. The PyTorch Geometric implementation of the GNN is publicly available on GitHub (https://github.com/PabloVD/HaloGraphNet).
|
|
|
Villanueva-Domingo, P., & Villaescusa-Navarro, F. (2021). Removing Astrophysics in 21 cm Maps with Neural Networks. Astrophys. J., 907(1), 44–14pp.
Abstract: Measuring temperature fluctuations in the 21 cm signal from the epoch of reionization and the cosmic dawn is one of the most promising ways to study the universe at high redshifts. Unfortunately, the 21 cm signal is affected by both cosmology and astrophysics processes in a nontrivial manner. We run a suite of 1000 numerical simulations with different values of the main astrophysical parameters. From these simulations we produce tens of thousands of 21 cm maps at redshifts 10 <= z <= 20. We train a convolutional neural network to remove the effects of astrophysics from the 21 cm maps and output maps of the underlying matter field. We show that our model is able to generate 2D matter fields not only that resemble the true ones visually but whose statistical properties agree with the true ones within a few percent down to scales 2 Mpc(-1). We demonstrate that our neural network retains astrophysical information that can be used to constrain the value of the astrophysical parameters. Finally, we use saliency maps to try to understand which features of the 21 cm maps the network is using in order to determine the value of the astrophysical parameters.
|
|
|
Villanueva-Domingo, P., Mena, O., & Palomares-Ruiz, S. (2021). A Brief Review on Primordial Black Holes as Dark Matter. Front. Astron. Space Sci., 8, 681084–10pp.
Abstract: Primordial black holes (PBHs) represent a natural candidate for one of the components of the dark matter (DM) in the Universe. In this review, we shall discuss the basics of their formation, abundance and signatures. Some of their characteristic signals are examined, such as the emission of particles due to Hawking evaporation and the accretion of the surrounding matter, effects which could leave an impact in the evolution of the Universe and the formation of structures. The most relevant probes capable of constraining their masses and population are discussed.
|
|
|
Villanueva-Domingo, P., Mena, O., & Miralda-Escude, J. (2020). Maximum amplitude of the high-redshift 21-cm absorption feature. Phys. Rev. D, 101(8), 083502–8pp.
Abstract: We examine the maximum possible strength of the global 21-cm absorption dip on the cosmic background radiation at high-redshift caused by the atomic intergalactic medium, when the Lyman-alpha coupling is maximum, assuming no exotic cooling mechanisms from interactions with dark matter. This maximum absorption is limited by three inevitable factors that need to be accounted for: (a) heating by energy transferred from the cosmic background radiation to the hydrogen atoms via 21-cm transitions, dubbed as 21-cm heating; (b) Ly alpha heating by scatterings of Ly alpha photons from the first stars; (c) the impact of the expected density fluctuations in the intergalactic gas in standard cold dark matter theory, which reduces the mean 21-cm absorption signal. Inclusion of this third novel effect reduces the maximum global 21-cm absorption by similar to 10%. Overall, the three effects studied here reduce the 21-cm global absorption by similar to 20% at z similar or equal to 17.
|
|
|
Villanueva-Domingo, P., & Ichiki, K. (2023). 21 cm forest constraints on primordial black holes. Publ. Astron. Soc. Jpn., 75(SP1), S33–S49.
Abstract: Primordial black holes (PBHs) as part of the dark matter (DM) would modify the evolution of large-scale structures and the thermal history of the universe. Future 21 cm forest observations, sensitive to small scales and the thermal state of the intergalactic medium (IGM), could probe the existence of such PBHs. In this article, we show that the shot noise isocurvature mode on small scales induced by the presence of PBHs can enhance the amount of low-mass halos, or minihalos, and thus, the number of 21 cm absorption lines. However, if the mass of PBHs is as large as M-PBH greater than or similar to 10 M-circle dot, with an abundant enough fraction of PBHs as DM, f(PBH), the IGM heating due to accretion on to the PBHs counteracts the enhancement due to the isocurvature mode, reducing the number of absorption lines instead. The concurrence of both effects imprints distinctive signatures on the number of absorbers, allowing the abundance of PBHs to be bound. We compute the prospects for constraining PBHs with future 21 cm forest observations, finding achievable competitive upper limits on the abundance as low as f(PBH) similar to 10(-3) at M-PBH = 100 M-circle dot, or even lower at larger masses, in regions of the parameter space unexplored by current probes. The impact of astrophysical X-ray sources on the IGM temperature is also studied, which could potentially weaken the bounds.
|
|
|
Villaescusa-Navarro, F. et al, & Villanueva-Domingo, P. (2022). The CAMELS Multifield Data Set: Learning the Universe's Fundamental Parameters with Artificial Intelligence. Astrophys. J. Suppl. Ser., 259(2), 61–14pp.
Abstract: We present the Cosmology and Astrophysics with Machine Learning Simulations (CAMELS) Multifield Data set (CMD), a collection of hundreds of thousands of 2D maps and 3D grids containing many different properties of cosmic gas, dark matter, and stars from more than 2000 distinct simulated universes at several cosmic times. The 2D maps and 3D grids represent cosmic regions that span similar to 100 million light-years and have been generated from thousands of state-of-the-art hydrodynamic and gravity-only N-body simulations from the CAMELS project. Designed to train machine-learning models, CMD is the largest data set of its kind containing more than 70 TB of data. In this paper we describe CMD in detail and outline a few of its applications. We focus our attention on one such task, parameter inference, formulating the problems we face as a challenge to the community. We release all data and provide further technical details at https://camels-multifield-dataset.readthedocs.io.
|
|
|
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 .
|
|
|
Vijande, J., Tedgren, A. C., Ballester, F., Baltas, D., Papagiannis, P., Rivard, M. J., et al. (2021). Source strength determination in iridium-192 and cobalt-60 brachytherapy: A European survey on the level of agreement between clinical measurements and manufacturer certificates. Phys. Imag. Radiat. Oncol., 19, 108–111.
Abstract: Background and purpose: Brachytherapy treatment outcomes depend on the accuracy of the delivered dose distribution, which is proportional to the reference air-kerma rate (RAKR). Current societal recommendations require the medical physicist to compare the measured RAKR values to the manufacturer source calibration certificate. The purpose of this work was to report agreement observed in current clinical practice in the European Union. Materials and methods: A European survey was performed for high- and pulsed-dose-rate (HDR and PDR) highenergy sources (Ir-192 and Co-60), to quantify observed RAKR differences. Medical physicists at eighteen hospitals from eight European countries were contacted, providing 1,032 data points from 2001 to 2020. Results: Over the survey period, 77% of the Ir-192 measurements used a well chamber instead of the older Krieger phantom method. Mean differences with the manufacturer calibration certificate were 0.01% +/- 1.15% for Ir-192 and -0.1% +/- 1.3% for Co-60. Over 95% of RAKR measurements in the clinic were within 3% of the manufacturer calibration certificate. Conclusions: This study showed that the agreement level was generally better than that reflected in prior societal recommendations positing 5%. Future recommendations on high-energy HDR and PDR source calibrations in the clinic may consider tightened agreements levels.
|
|
|
Vijande, J., Carmona, V., Lliso, F., Ballester, F., & Perez-Calatayud, J. (2024). An efficient component of the redundancy calibration program to ensure equipment stability by assaying HDR Ir-192 sources at the time of replacement. J. Appl. Clin. Med. Phys, , 5pp.
Abstract: BackgroundBrachytherapy (BT) treatments involving temporary high-dose rate (HDR) sources are extensively employed in clinical practice. Ensuring the consistency of all measurement equipment at the hospital level is crucial, requiring a robust redundancy and consistency program. This enables the institution to verify the stability of the dosimetry system over time.PurposeTo describe, justify, and analyze a component of the redundancy program of the calibration protocols followed by the Radiotherapy Department of the Hospital Universitari i Polit & egrave;cnic La Fe (Val & egrave;ncia, Spain) during the last 10 years for the case of HDR BT as an additional component to ensure long term stability of the measurement equipment.MethodsAt the time the HDR BT source is replaced, its Air Kerma Strength (SK) is measured. By comparing this value with the one obtained at the time of installation (corrected by decay), a clear determination of the stability of the measurement equipment can be performed.ResultsDifference between SK,vendor and SK,hosp as a function of the measurement date is reported for a 10 years' period. All measurements are well within the +/- 3% tolerance level recommended in current international guidelines. Percentage differences of SK,hosp values at the time of replacement compared to SK,hosp ones at the time when the source was installed are within the +/- 0.5% range, reflecting oscillations around a null deviation.ConclusionsThe method proposed allows any hospital to ensure a redundancy component of the long-term stability of all equipment involved in BT measurements in a very simple and time efficient manner. Additionally, it enables the hospital to maintain a detailed log of historical differences, facilitating the identification and correction of potential systematic deviations over time.
|
|