Search

The connection between galaxies and dark matter halos encompasses a range of processes and play a pivotal role in our understanding of galaxy formation and evolution. Traditionally, this link has been established through physical or empirical models. Machine learning techniques are adaptable tools that handle high-dimensional data and grasp associations between numerous attributes. In particular, probabilistic models capture the stochasticity inherent to these complex relations. We compare different probabilistic machine learning methods to model the uncertainty in the halo-galaxy connection and efficiently generate galaxy catalogs that faithfully resemble the reference sample by predicting joint distributions of central galaxy properties conditioned to their host halo features. The analysis is based on the IllustrisTNG300 simulation. The methods model the distributions in different ways. We compare a multilayer perceptron that predicts the parameters of a multivariate Gaussian distribution, a multilayer perceptron classifier, and the method of normalizing flows. The classifier predicts the parameters of a Categorical distribution, which are defined in a high-dimensional parameter space through a Voronoi cell-based hierarchical scheme. We evaluate the model's performances under various sample selections based on halo properties. The three methods exhibit comparable results, with normalizing flows showing the best performance in most scenarios. The models reproduce the main features of galaxy properties distributions with high-fidelity and reproduce the results obtained with traditional, deterministic, estimators. Our results also indicate that different halos and galaxy populations are subject to varying degrees of stochasticity, which has relevant implications for studies of large-scale structure., Comment: 14 pages, 8 figures, 6 tables

It has been recently shown that a powerful way to constrain cosmological parameters from galaxy redshift surveys is to train graph neural networks to perform field-level likelihood-free inference without imposing cuts on scale. In particular, de Santi et al. (2023) developed models that could accurately infer the value of $\Omega_{\rm m}$ from catalogs that only contain the positions and radial velocities of galaxies that are robust to uncertainties in astrophysics and subgrid models. However, observations are affected by many effects, including 1) masking, 2) uncertainties in peculiar velocities and radial distances, and 3) different galaxy selections. Moreover, observations only allow us to measure redshift, intertwining galaxies' radial positions and velocities. In this paper we train and test our models on galaxy catalogs, created from thousands of state-of-the-art hydrodynamic simulations run with different codes from the CAMELS project, that incorporate these observational effects. We find that, although the presence of these effects degrades the precision and accuracy of the models, and increases the fraction of catalogs where the model breaks down, the fraction of galaxy catalogs where the model performs well is over 90 %, demonstrating the potential of these models to constrain cosmological parameters even when applied to real data., Comment: 39 pages, 25 figures. For the reference in the abstract (de Santi et al. 2023) see arXiv:2302.14101

We present CAMELS-ASTRID, the third suite of hydrodynamical simulations in the Cosmology and Astrophysics with MachinE Learning (CAMELS) project, along with new simulation sets that extend the model parameter space based on the previous frameworks of CAMELS-TNG and CAMELS-SIMBA, to provide broader training sets and testing grounds for machine-learning algorithms designed for cosmological studies. CAMELS-ASTRID employs the galaxy formation model following the ASTRID simulation and contains 2,124 hydrodynamic simulation runs that vary 3 cosmological parameters ($\Omega_m$, $\sigma_8$, $\Omega_b$) and 4 parameters controlling stellar and AGN feedback. Compared to the existing TNG and SIMBA simulation suites in CAMELS, the fiducial model of ASTRID features the mildest AGN feedback and predicts the least baryonic effect on the matter power spectrum. The training set of ASTRID covers a broader variation in the galaxy populations and the baryonic impact on the matter power spectrum compared to its TNG and SIMBA counterparts, which can make machine-learning models trained on the ASTRID suite exhibit better extrapolation performance when tested on other hydrodynamic simulation sets. We also introduce extension simulation sets in CAMELS that widely explore 28 parameters in the TNG and SIMBA models, demonstrating the enormity of the overall galaxy formation model parameter space and the complex non-linear interplay between cosmology and astrophysical processes. With the new simulation suites, we show that building robust machine-learning models favors training and testing on the largest possible diversity of galaxy formation models. We also demonstrate that it is possible to train accurate neural networks to infer cosmological parameters using the high-dimensional TNG-SB28 simulation set.

We discover analytic equations that can infer the value of $\Omega_{\rm m}$ from the positions and velocity moduli of halo and galaxy catalogues. The equations are derived by combining a tailored graph neural network (GNN) architecture with symbolic regression. We first train the GNN on dark matter halos from Gadget N-body simulations to perform field-level likelihood-free inference, and show that our model can infer $\Omega_{\rm m}$ with $\sim6\%$ accuracy from halo catalogues of thousands of N-body simulations run with six different codes: Abacus, CUBEP$^3$M, Gadget, Enzo, PKDGrav3, and Ramses. By applying symbolic regression to the different parts comprising the GNN, we derive equations that can predict $\Omega_{\rm m}$ from halo catalogues of simulations run with all of the above codes with accuracies similar to those of the GNN. We show that by tuning a single free parameter, our equations can also infer the value of $\Omega_{\rm m}$ from galaxy catalogues of thousands of state-of-the-art hydrodynamic simulations of the CAMELS project, each with a different astrophysics model, run with five distinct codes that employ different subgrid physics: IllustrisTNG, SIMBA, Astrid, Magneticum, SWIFT-EAGLE. Furthermore, the equations also perform well when tested on galaxy catalogues from simulations covering a vast region in parameter space that samples variations in 5 cosmological and 23 astrophysical parameters. We speculate that the equations may reflect the existence of a fundamental physics relation between the phase-space distribution of generic tracers and $\Omega_{\rm m}$, one that is not affected by galaxy formation physics down to scales as small as $10~h^{-1}{\rm kpc}$., Comment: 32 pages, 13 figures, summary video: https://youtu.be/STZHvDHkVgo

We train graph neural networks to perform field-level likelihood-free inference using galaxy catalogs from state-of-the-art hydrodynamic simulations of the CAMELS project. Our models are rotational, translational, and permutation invariant and do not impose any cut on scale. From galaxy catalogs that only contain $3$D positions and radial velocities of $\sim 1, 000$ galaxies in tiny $(25~h^{-1}{\rm Mpc})^3$ volumes our models can infer the value of $\Omega_{\rm m}$ with approximately $12$ % precision. More importantly, by testing the models on galaxy catalogs from thousands of hydrodynamic simulations, each having a different efficiency of supernova and AGN feedback, run with five different codes and subgrid models - IllustrisTNG, SIMBA, Astrid, Magneticum, SWIFT-EAGLE -, we find that our models are robust to changes in astrophysics, subgrid physics, and subhalo/galaxy finder. Furthermore, we test our models on $1,024$ simulations that cover a vast region in parameter space - variations in $5$ cosmological and $23$ astrophysical parameters - finding that the model extrapolates really well. Our results indicate that the key to building a robust model is the use of both galaxy positions and velocities, suggesting that the network have likely learned an underlying physical relation that does not depend on galaxy formation and is valid on scales larger than $\sim10~h^{-1}{\rm kpc}$., Comment: 34 pages, 12 figures. For a video summarizing the results, see https://youtu.be/b59ep7cyPOs

The relationship between galaxies and haloes is central to the description of galaxy formation, and a fundamental step towards extracting precise cosmological information from galaxy maps. However, this connection involves several complex processes that are interconnected. Machine Learning methods are flexible tools that can learn complex correlations between a large number of features, but are traditionally designed as deterministic estimators. In this work, we use the IllustrisTNG300-1 simulation and apply neural networks in a binning classification scheme to predict probability distributions of central galaxy properties, namely stellar mass, colour, specific star formation rate, and radius, using as input features the halo mass, concentration, spin, age, and the overdensity on a scale of 3 $h^{-1}$ Mpc. The model captures the intrinsic scatter in the relation between halo and galaxy properties, and can thus be used to quantify the uncertainties related to the stochasticity of the galaxy properties with respect to the halo properties. In particular, with our proposed method, one can define and accurately reproduce the properties of the different galaxy populations in great detail. We demonstrate the power of this tool by directly comparing traditional single-point estimators and the predicted joint probability distributions, and also by computing the power spectrum of a large number of tracers defined on the basis of the predicted colour-stellar mass diagram. We show that the neural networks reproduce clustering statistics of the individual galaxy populations with excellent precision and accuracy., Comment: 12 pages, 7 figures

Cosmological covariance matrices are fundamental for parameter inference, since they are responsible for propagating uncertainties from the data down to the model parameters. However, when data vectors are large, in order to estimate accurate and precise matrices we need huge numbers of observations, or rather costly simulations - neither of which may be viable. In this work we propose a machine learning approach to alleviate this problem in the context of the matrices used in the study of large-scale structure. With only a small amount of data (matrices built with samples of 50-200 halo power spectra) we are able to provide significantly improved matrices, which are almost indistinguishable from the ones built from much larger samples (thousands of spectra). In order to perform this task we trained convolutional neural networks to denoise the matrices, using in the training process a data set made up entirely of spectra extracted from simple, inexpensive halo simulations (mocks). We then show that the method not only removes the noise in the matrices of the cheap simulation, but it is also able to successfully denoise the matrices of halo power spectra from N-body simulations. We compare the denoised to the other matrices using several metrics, and in all of them they score better, without any signs of spurious artifacts. With the help of the Wishart distribution we derive an analytical extrapolation for the effective sample augmentation allowed by the denoiser. Finally, we show that, by using the denoised matrices, the cosmological parameters can be recovered with nearly the same accuracy as when using matrices built with a sample of 30,000 spectra in the case of the cheap simulations, and with 15,000 spectra in the case of the N-body simulations. Of particular interest is the bias in the Hubble parameter $H_0$, which was significantly reduced after applying the denoiser., Comment: Matches published version; very minor changes wrt V1

Elucidating the connection between the properties of galaxies and the properties of their hosting haloes is a key element in galaxy formation. When the spatial distribution of objects is also taken under consideration, it becomes very relevant for cosmological measurements. In this paper, we use machine learning techniques to analyse these intricate relations in the IllustrisTNG300 magnetohydrodynamical simulation, predicting baryonic properties from halo properties. We employ four different algorithms: extremely randomized trees, K-nearest neighbours, light gradient boosting machine, and neural networks, along with a unique and powerful combination of the results from all four approaches. Overall, the different algorithms produce consistent results in terms of predicting galaxy properties from a set of input halo properties that include halo mass, concentration, spin, and halo overdensity. For stellar mass, the Pearson correlation coefficient is 0.98, dropping down to 0.7-0.8 for specific star formation rate (sSFR), colour, and size. In addition, we apply, for the first time in this context, an existing data augmentation method, synthetic minority over-sampling technique for regression with Gaussian noise (SMOGN), designed to alleviate the problem of imbalanced data sets, showing that it improves the overall shape of the predicted distributions and the scatter in the halo-galaxy relations. We also demonstrate that our predictions are good enough to reproduce the power spectra of multiple galaxy populations, defined in terms of stellar mass, sSFR, colour, and size with high accuracy. Our results align with previous reports suggesting that certain galaxy properties cannot be reproduced using halo features alone., Comment: Matches published version; very minor changes wrt V1

Results from the gas gain monitoring (GGM) system for the muon detector using RPC in the CMS experiment at the LHC is presented. The system is designed to provide fast and accurate determination of any shift in the working point of the chambers due to gas mixture changes., Comment: 8 pages, 7 figures, uses lnfprepCMS.sty, presented by S. Colafranceschi at IEEE 08 Dresden (Germany) November 2008

The status of the CMS RPC Gas Gain Monitoring (GGM) system developed at the Frascati Laboratory of INFN (Istituto Nazionale di Fisica Nucleare) is reported on. The GGM system is a cosmic ray telescope based on small RPC detectors operated with the same gas mixture used by the CMS RPC system. The GGM gain and efficiency are continuously monitored on-line, thus providing a fast and accurate determination of any shift in working point conditions. The construction details and the first result of GGM commissioning are described., Comment: 8 pages, 9 figures, uses lnfprepCMS.sty, presented by L. Benussi at RPC07, Mumbai, INDIA 2008

A practical scheme for entanglement creation between distant atoms located inside a single-mode optical cavity is discussed. We show that the degree of entanglement and the time it takes for the entanglement to reach its optimum value is a sensitive function the initial conditions and the position of the atoms inside the cavity mode. It is found that the entangled properties of the two atoms can readily be extracted from dynamics of a simple two-level system. Effectively, we engineer two coupled qubits whose the dynamics are analogous to that of a driven single two-level system. It is found that spatial variations of the coupling constants actually help to create transient entanglement which may appear on the time scale much longer than that predicted for the case of equal coupling constants. When the atoms are initially prepared in an entangled state, they may remain entangled for all times. We also find that the entanglement exhibits an interesting phenomenon of diffraction when the the atoms are located between the nodes and antinodes of the cavity mode. The diffraction pattern of the entanglement varies with time and we explain this effect in terms of the quantum property of complementarity, which is manifested as a tradeoff between the knowledge of energy of the exchanged photon versus the evolution time of the system., Comment: Phys. Rev. A75, 042307 (2007)

The Resistive Plate Chambers (RPC) detector of the CMS experiment at the LHC proton collider (CERN, Switzerland) will employ an online gas analysis and monitoring system of the freon-based gas mixture used. We give an overview of the CMS RPC gas system, describe the project parameters and first results on gas-chromatograph analysis. Finally, we report on preliminary results for a set of monitor RPC., Comment: 9 pages, 8 figures. Presented by Stefano Bianco (Laboratori Nazionali di Frascati dell'INFN) at the IEEE NSS, San Diego (USA), October 2006

In this research, we report on the compositional, microstructural and crystallographic properties of a lead coin which has been regarded for many years as a genuine silver coin minted in the Southern Italy in the course of the 4th century BC. The material characterisation of this object allowed detecting an ancient forging technology, not previously reported, which was meant for the silvering of lead substrates The data collected have disclosed a contemporary counterfeiting procedure based on a metal coating process onto a Pb substrate. This coating has been identified as a bi-layer with a Cu innermost and an Ag outermost visible layer. As far as the coating application technique is concerned, the gathered evidence has clearly indicated that the original appearance of this artifact cannot be explained in terms of any of the established methods for the growth of an artificially silvered coating in classical antiquity. This technology is now being explained in terms of modern, fully non destructive scientific methods., Comment: 9 pages including 8 figures. To be presented at the Cavallino Archaeometry Workshop 2006, Cavallino (Lecce - Italy) - May 22-25, 2006