Modelling the Halo mass number counts in modified gravity cosmological scenarios

In this Thesis work we investigate some of different cosmological background scenarios using one of the main probes used in cosmology: the halo mass function. The observed abundance of galaxy clusters (or similarly DM haloes) can indeed be compared to its theoretical predictions to derive fundamental constrains on the cosmological scenario assumed. Given the importance of exploring and constraining models degenerate with the ΛCDM one, we test the applicability of some notable halo mass function models to these scenarios. To this purpose, we made use of the DUSTGRAIN-pathfinder N-body simulations, which assume cosmological scenarios that include modified gravity in the form of f(R) models and massive neutrinos. We carried on the analysis of 3 simulation snapshots at different redshifts, z = 0, 0.5, 1, building multiple samples of dark matter haloes by applying different overdensity thresholds during the procedure of halo identification. We started our analysis by considering the halo mass function model introduced by Despali et al. (2016), who proposed a parametrization that encapsulates the effect of the different halo mass definitions and the relative evolution with the redshift. We calibrated the main parameters of this relation by using the ΛCDM halo catalogues extracted from the DUSTGRAIN-pathfinder simulations, fitting the measured halo abundances at all redshifts and density thresholds. Afterwards we tested the same model parametrization with halo catalogues extracted from the simulations implementing both modified gravity and massive neutrinos. We repeated therefore the calibration procedure on these data to search for discrepancies with respect to the ΛCDM model. Finally we focused the analysis on the cosmological models implementing modified gravity only. We took our ΛCDM calibrated halo mass function and we modified it with the additional f (R) gravity form proposed by Gupta et al. (2022).

Back-in-time void finder: dynamical detection of cosmic voids for precision cosmology

Cosmic voids are vast and underdense regions emerging between the elements of the cosmic web and dominating the large-scale structure of the Universe. Void number counts and density profiles have been demonstrated to provide powerful cosmological probes. Indeed, thanks to their low-density nature and they very large sizes, voids represent natural laboratories to test alternative dark energy scenarios, modifications of gravity and the presence of massive neutrinos. Despite the increasing use of cosmic voids in Cosmology, a commonly accepted definition for these objects has not yet been reached. For this reason, different void finding algorithms have been proposed during the years. Voids finder algorithms based on density or geometrical criteria are affected by intrinsic uncertainties. In recent years, new solutions have been explored to face these issues. The most interesting is based on the idea of identify void positions through the dynamics of the mass tracers, without performing any direct reconstruction of the density field. The goal of this Thesis is to provide a performing void finder algorithm based on dynamical criteria. The Back-in-time void finder (BitVF) we present use tracers as test particles and their orbits are reconstructed from their actual clustered configuration to an homogeneous and isotropic distribution, expected for the Universe early epoch. Once the displacement field is reconstructed, the density field is computed as its divergence. Consequently, void centres are identified as local minima of the field. In this Thesis work we applied the developed void finding algorithm to simulations. From the resulting void samples we computed different void statistics, comparing the results to those obtained with VIDE, the most popular void finder. BitVF proved to be able to produce a more reliable void samples than the VIDE ones. The BitVF algorithm will be a fundamental tool for precision cosmology, especially with upcoming galaxy-survey.

Mock HI observations of Milky Way-like simulated galaxies

L'obiettivo di questo lavoro di tesi consiste nel descrivere sia il processo necessario per la creazione di osservazioni sintetiche di galassie simulate simili alla Via Lattea nella riga di emissione a 21 cm dell'idrogeno neutro (HI), sia il lavoro di analisi fondamentale che serve a confrontare in modo efficace l'output generato con delle osservazioni di galassie reali. Come prima cosa è descritta la teoria quantistica che sta alla base dell'emissione a 21 cm di HI, illustrando l'importanza di tale riga di emissione nell'ambito dell'astronomia e come si possano ottenere informazioni fondamentali sulle sorgenti di questa radiazione a partire dai dati osservativi. Il lavoro poi si focalizza sull'utilizzo del software MARTINI per la creazione di osservazioni sintetiche della linea a 21 cm per una galassia simulata con proprietà simili alla Via Lattea generata utilizzando il modello numerico SMUGGLE. Infine, si passa ad una breve descrizione dell'analisi dei dati sintetici creati, e al loro confronto con dei dati provenienti da osservazioni reali di galassie con proprietà simili, per ottenere una valutazione qualitativa della bontà del modello SMUGGLE impiegato nella simulazione numerica.

Dynamics of galaxies in the cluster environment: a resonant origin for the ISP?

The internal dynamics of elliptical galaxies in clusters depends on many factors, including the environment in which the galaxy is located. In addition to the strong encounters with the other galaxies, we can also consider the gravitational interaction with the ubiquitous Cluster Tidal Field (CTF). As recognized in many studies, one possible way in which CTF affects the dynamics of galaxies inside the cluster is related to the fact that they may start oscillating as “rigid bodies” around their equilibrium positions in the field, with the periods of these oscillations curiously similar to those of stellar orbits in the outer parts of galaxies. Resonances between the two motions are hence expected and this phenomenon could significantly contribute to the formation of the Intracluster Stellar Population (ISP), whose presence is abundantly confirmed by observations. In this thesis work, we propose to study the motion of an elliptical galaxy, modelled as a rigid body, in the CTF, especially when its center of mass traces a quasi-circular orbit in the cluster gravitational potential. This case extends and generalizes the previous models and findings, proceeding towards a much more realistic description of galaxy motion. In addition to this, the presence of a further oscillation, namely that of the entire galaxy along its orbit, will possibly increase the probability of having resonances and, consequently, the rate of ISP production nearly to observed values. Thus, after reviewing the dynamics of a rigid body in a generic force field, we will assess some physically relevant studies and report their main results, discussing their implications with respect to our problem. We will conclude our discussion focusing on the more realistic scenario of an elliptical galaxy whose center of mass moves on a quasi-circular orbit in a spherically symmetric potential. The derivation of the fundamental equations of motion will serve as the basis for future modelling and discussions.

Magnetic field in the outskirts of PSZ2 G096.88+24.18 from depolarization analysis of radio relics

Radio relics are one of the different types of diffuse radio sources present in a fraction of galaxy clusters. They are characterized by elongated arc-like shapes, with sizes that range between 0.5 and 2 Mpc, and highly polarized emission (up to ∼60%) at GHz frequencies The linearly polarized radiation of relics, moving through a magnetized plasma which is the ICM, is affected by the rotation of the linear polarization vector. This effect, known as “Faraday rotation”, can cause depolarization. The study of this effect allows us to constrain the magnetic field projected along the line of sight. The aim of this thesis work is to constrain the magnetic field intensity and distribution in the periphery of the cluster PSZ2 G096.88+24.18: this cluster hosts a pair of radio relics that can be used for polarization analysis. To analyse the polarization properties of the relics in PSZ2 G096.88+24.18 radio relics we used new Jansky Very Large Array (VLA) observations together with archival observations. The polarization study has been performed using the Rotation Measure Synthesis technique, which allows us to recover polarization, minimizing the bandwidth depolarization. Thanks to this technique, we recovered more polarization from the southern relic (with respect to provious works), We studied also the depolarization trend with the resolution for the southern relic, and found that the polarization fraction decreases with the beamsize. Finally, we have produced simulated magnetic fields models, varying the auto-correlation lengths of the magnetic field, in order to reproduce the observed depolarization trend in the southern relic. Comparing our observational results and model predictions, we were able to constrain the scales over which the turbulent magnetic field varies within the cluster. We conclude that the depolarization observed in the southern relic is likely due to external depolarization caused by the magnetized ICM distribution within the cluster.

Jeans modeling of massive early-type galaxies with observed stellar kinematic profiles up to 2-3 effective radii

A recent integral-field spectroscopic (IFS) survey, the MASSIVE survey (Ma et al. 2014), observed the 116 most massive (MK < −25.3 mag, stellar mass M∗ > 10^11.6 M⊙) early-type galaxies (ETGs) within 108 Mpc, out to radii as large as 40 kpc, that correspond to ∼ 2 − 3 effective radii (Re). One of the major findings of the MASSIVE survey is that the galaxy sample is split nearly equally among three groups showing three different velocity dispersion profiles σ(R) outer of a radius ∼ 5 kpc (falling, flat and rising with radius). The purpose of this thesis is to model the kinematic profiles of six ETGs included in the MASSIVE survey and representative of the three observed σ(R) shapes, with the aim of investigating their dynamical structure. Models for the chosen galaxies are built using the numerical code JASMINE (Posacki, Pellegrini, and Ciotti 2013). The code produces models of axisymmetric galaxies, based on the solution of the Jeans equations for a multicomponent gravitational potential (supermassive black hole, stars and dark matter halo). With the aim of having a good agreement between the kinematics obtained from the Jeans equations, and the observed σ and rotation velocity V of MASSIVE (Veale et al. 2016, 2018), I derived constraints on the dark matter distribution and orbital anisotropy. This work suggests a trend of the dark matter amount and distribution with the shape of the velocity dispersion profiles in the outer regions: the models of galaxies with flat or rising velocity dispersion profiles show higher dark matter fractions fDM both within 1 Re and 5 Re. Orbital anisotropy alone cannot account for the different observed trends of σ(R) and has a minor effect compared to variations of the mass profile. Galaxies with similar stellar mass M∗ that show different velocity dispersion profiles (from falling to rising) are successfully modelled with a variation of the halo mass Mh.

Constraining FIR size, dust temperature and dust mass in ALPINE galaxies at redshift 4

The study of galaxies at high redshift plays a crucial role to understand the mechanism of galaxy formation and evolution. At redshifts just after the epoch of re-ionization (4

Spectral distortions of the cosmic microwave background and small-scale gravitational waves as new windows on the Early Universe

In this thesis, we explore constraints which can be put on the primordial power spectrum of curvature perturbations beyond the scales probed by anisotropies of the cosmic microwave background and galaxy surveys. We exploit present and future measurements of CMB spectral distortions, and their synergy with CMB anisotropies, as well existing and future upper limits on the stochastic background of gravitational waves. We derive for the first time phenomenological templates that fit small-scale bumps in the primordial power spectrum generated in multi-field models of inflation. By using such templates, we study for the first time imprints of primordial peaks on anisotropies and spectral distortions of the cosmic microwave background and we investigate their contribution to the stochastic background of gravitational waves. Through a Monte Carlo Markov Chain analysis we infer for the first time the constraints on the amplitude, the width and the location of such bumps using Planck and FIRAS data. We also forecast how a future spectrometer like PIXIE could improve FIRAS boundaries. The results derived in this thesis have implications for the possibility of primordial black holes from inflation.