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 .

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).

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.

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.

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.