Study of cosmic nuclei fluxes with AMS-02: implication for dark matter search

L’Alpha Magnetic Spectrometer (AMS-02) é un rivelatore per raggi cosmici (CR) progettato e costruito da una collaborazione internazionale di 56 istituti e 16 paesi ed installato il 19 Maggio del 2011 sulla Stazione Spaziale Internazionale (ISS). Orbitando intorno alla Terra, AMS-02 sará in grado di studiare con un livello di accuratezza mai raggiunto prima la composizione dei raggi cosmici, esplorando nuove frontiere nella fisica delle particelle, ricercando antimateria primordiale ed evidenze indirette di materia oscura. Durante il mio lavoro di tesi, ho utilizzato il software GALPROP per studiare la propagazione dei CR nella nostra Galassia attraverso il mezzo interstellare (ISM), cercando di individuare un set di parametri in grado di fornire un buon accordo con i dati preliminari di AMS-02. In particolare, mi sono dedicata all’analisi del processo di propagazione di nuclei, studiando i loro flussi e i relativi rapporti. Il set di propagazione ottenuto dall’analisi é stato poi utilizzato per studiare ipotetici flussi da materia oscura e le possibili implicazioni per la ricerca indiretta attraverso AMS-02.

Zoomed simulations of Halo segregation in cosmological models with two species of coupled dark matter

La materia ordinaria copre soli pochi punti percentuali della massa-energia totale dell'Universo, che è invece largamente dominata da componenti “oscure”. Il modello standard usato per descriverle è il modello LambdaCDM. Nonostante esso sembri consistente con la maggior parte dei dati attualmente disponibili, presenta alcuni problemi fondamentali che ad oggi restano irrisolti, lasciando spazio per lo studio di modelli cosmologici alternativi. Questa Tesi mira a studiare un modello proposto recentemente, chiamato “Multi-coupled Dark Energy” (McDE), che presenta interazioni modificate rispetto al modello LambdaCDM. In particolare, la Materia Oscura è composta da due diversi tipi di particelle con accoppiamento opposto rispetto ad un campo scalare responsabile dell'Energia Oscura. L'evoluzione del background e delle perturbazioni lineari risultano essere indistinguibili da quelle del modello LambdaCDM. In questa Tesi viene presentata per la prima volta una serie di simulazioni numeriche “zoomed”. Esse presentano diverse regioni con risoluzione differente, centrate su un singolo ammasso di interesse, che permettono di studiare in dettaglio una singola struttura senza aumentare eccessivamente il tempo di calcolo necessario. Un codice chiamato ZInCo, da me appositamente sviluppato per questa Tesi, viene anch'esso presentato per la prima volta. Il codice produce condizioni iniziali adatte a simulazioni cosmologiche, con differenti regioni di risoluzione, indipendenti dal modello cosmologico scelto e che preservano tutte le caratteristiche dello spettro di potenza imposto su di esse. Il codice ZInCo è stato usato per produrre condizioni iniziali per una serie di simulazioni numeriche del modello McDE, le quali per la prima volta mostrano, grazie all'alta risoluzione raggiunta, che l'effetto di segregazione degli ammassi avviene significativamente prima di quanto stimato in precedenza. Inoltre, i profili radiale di densità ottenuti mostrano un appiattimento centrale nelle fasi iniziali della segregazione. Quest'ultimo effetto potrebbe aiutare a risolvere il problema “cusp-core” del modello LambdaCDM e porre limiti ai valori dell'accoppiamento possibili.

Properties of cosmic voids in dark energy simulations

The last decade has witnessed the establishment of a Standard Cosmological Model, which is based on two fundamental assumptions: the first one is the existence of a new non relativistic kind of particles, i. e. the Dark Matter (DM) that provides the potential wells in which structures create, while the second one is presence of the Dark Energy (DE), the simplest form of which is represented by the Cosmological Constant Λ, that sources the acceleration in the expansion of our Universe. These two features are summarized by the acronym ΛCDM, which is an abbreviation used to refer to the present Standard Cosmological Model. Although the Standard Cosmological Model shows a remarkably successful agreement with most of the available observations, it presents some longstanding unsolved problems. A possible way to solve these problems is represented by the introduction of a dynamical Dark Energy, in the form of the scalar field ϕ. In the coupled DE models, the scalar field ϕ features a direct interaction with matter in different regimes. Cosmic voids are large under-dense regions in the Universe devoided of matter. Being nearby empty of matter their dynamics is supposed to be dominated by DE, to the nature of which the properties of cosmic voids should be very sensitive. This thesis work is devoted to the statistical and geometrical analysis of cosmic voids in large N-body simulations of structure formation in the context of alternative competing cosmological models. In particular we used the ZOBOV code (see ref. Neyrinck 2008), a publicly available void finder algorithm, to identify voids in the Halos catalogues extraxted from CoDECS simulations (see ref. Baldi 2012 ). The CoDECS are the largest N-body simulations to date of interacting Dark Energy (DE) models. We identify suitable criteria to produce voids catalogues with the aim of comparing the properties of these objects in interacting DE scenarios to the standard ΛCDM model, at different redshifts. This thesis work is organized as follows: in chapter 1, the Standard Cosmological Model as well as the main properties of cosmic voids are intro- duced. In chapter 2, we will present the scalar field scenario. In chapter 3 the tools, the methods and the criteria by which a voids catalogue is created are described while in chapter 4 we discuss the statistical properties of cosmic voids included in our catalogues. In chapter 5 the geometrical properties of the catalogued cosmic voids are presented by means of their stacked profiles. In chapter 6 we summarized our results and we propose further developments of this work.

Asymmetric drift and halo flattening in disk galaxies

The main result in this work is the solution of the Jeans equations for an axisymmetric galaxy model containing a baryonic component (distributed according to a Miyamoto-Nagai profile) and a dark matter halo (described by the Binney logarithmic potential). The velocity dispersion, azimuthal velocity and some other interesting quantities such as the asymmetric drift are studied, along with the influence of the model parameters on these (observable) quantities. We also give an estimate for the velocity of the radial flow, caused by the asymmetric drift. Other than the mathematical beauty that lies in solving a model analytically, the interest of this kind of results can be mainly found in numerical simulations that study the evolution of gas flows. For example, it is important to know how certain parameters such as the shape (oblate, prolate, spherical) of a dark matter halo, or the flattening of the baryonic matter, or the mass ratio between dark and baryonic matter, have an influence on observable quantities such as the velocity dispersion. In the introductory chapter, we discuss the Jeans equations, which provide information about the velocity dispersion of a system. Next we will consider some dynamical quantities that will be useful in the rest of the work, e.g. the asymmetric drift. In Chapter 2 we discuss in some more detail the family of galaxy models we studied. In Chapter 3 we give the solution of the Jeans equations. Chapter 4 describes and illustrates the behaviour of the velocity dispersion, as a function of the several parameters, along with asymptotic expansions. In Chapter 5 we will investigate the behaviour of certain dynamical quantities for this model. We conclude with a discussion in Chapter 6.

Study of the sensitivity of the XENON1T experiment with the profile likelihood method

Oggi sappiamo che la materia ordinaria rappresenta solo una piccola parte dell'intero contenuto in massa dell'Universo. L'ipotesi dell'esistenza della Materia Oscura, un nuovo tipo di materia che interagisce solo gravitazionalmente e, forse, tramite la forza debole, è stata avvalorata da numerose evidenze su scala sia galattica che cosmologica. Gli sforzi rivolti alla ricerca delle cosiddette WIMPs (Weakly Interacting Massive Particles), il generico nome dato alle particelle di Materia Oscura, si sono moltiplicati nel corso degli ultimi anni. L'esperimento XENON1T, attualmente in costruzione presso i Laboratori Nazionali del Gran Sasso (LNGS) e che sarà in presa dati entro la fine del 2015, segnerà un significativo passo in avanti nella ricerca diretta di Materia Oscura, che si basa sulla rivelazione di collisioni elastiche su nuclei bersaglio. XENON1T rappresenta la fase attuale del progetto XENON, che ha già realizzato gli esperimenti XENON10 (2005) e XENON100 (2008 e tuttora in funzione) e che prevede anche un ulteriore sviluppo, chiamato XENONnT. Il rivelatore XENON1T sfrutta circa 3 tonnellate di xeno liquido (LXe) e si basa su una Time Projection Chamber (TPC) a doppia fase. Dettagliate simulazioni Monte Carlo della geometria del rivelatore, assieme a specifiche misure della radioattività dei materiali e stime della purezza dello xeno utilizzato, hanno permesso di predire con accuratezza il fondo atteso. In questo lavoro di tesi, presentiamo lo studio della sensibilità attesa per XENON1T effettuato tramite il metodo statistico chiamato Profile Likelihood (PL) Ratio, il quale nell'ambito di un approccio frequentista permette un'appropriata trattazione delle incertezze sistematiche. In un primo momento è stata stimata la sensibilità usando il metodo semplificato Likelihood Ratio che non tiene conto di alcuna sistematica. In questo modo si è potuto valutare l'impatto della principale incertezza sistematica per XENON1T, ovvero quella sulla emissione di luce di scintillazione dello xeno per rinculi nucleari di bassa energia. I risultati conclusivi ottenuti con il metodo PL indicano che XENON1T sarà in grado di migliorare significativamente gli attuali limiti di esclusione di WIMPs; la massima sensibilità raggiunge una sezione d'urto σ=1.2∙10-47 cm2 per una massa di WIMP di 50 GeV/c2 e per una esposizione nominale di 2 tonnellate∙anno. I risultati ottenuti sono in linea con l'ambizioso obiettivo di XENON1T di abbassare gli attuali limiti sulla sezione d'urto, σ, delle WIMPs di due ordini di grandezza. Con tali prestazioni, e considerando 1 tonnellata di LXe come massa fiduciale, XENON1T sarà in grado di superare gli attuali limiti (esperimento LUX, 2013) dopo soli 5 giorni di acquisizione dati.

Axion-like particles and the 3.5 keV line in 4D string models

This work is focused on axions and axion like particles (ALPs) and their possible relation with the 3.55 keV photon line detected, in recent years, from galaxy clusters and other astrophysical objects. We focus on axions that come from string compactification and we study the vacuum structure of the resulting low energy 4D N=1 supergravity effective field theory. We then provide a model which might explain the 3.55 keV line through the following processes. A 7.1 keV dark matter axion decays in two light axions, which, in turn, are transformed into photons thanks to the Primakoff effect and the existence of a kinetic mixing between two U(1)s gauge symmetries belonging respectively to the hidden and the visible sector. We present two models, the first one gives an outcome inconsistent with experimental data, while the second can yield the desired result.

Modelling the size distribution of cosmic voids

The Standard Cosmological Model is generally accepted by the scientific community, there are still an amount of unresolved issues. From the observable characteristics of the structures in the Universe,it should be possible to impose constraints on the cosmological parameters. Cosmic Voids (CV) are a major component of the LSS and have been shown to possess great potential for constraining DE and testing theories of gravity. But a gap between CV observations and theory still persists. A theoretical model for void statistical distribution as a function of size exists (SvdW) However, the SvdW model has been unsuccesful in reproducing the results obtained from cosmological simulations. This undermines the possibility of using voids as cosmological probes. The goal of our thesis work is to cover the gap between theoretical predictions and measured distributions of cosmic voids. We develop an algorithm to identify voids in simulations,consistently with theory. We inspecting the possibilities offered by a recently proposed refinement of the SvdW (the Vdn model, Jennings et al., 2013). Comparing void catalogues to theory, we validate the Vdn model, finding that it is reliable over a large range of radii, at all the redshifts considered and for all the cosmological models inspected. We have then searched for a size function model for voids identified in a distribution of biased tracers. We find that, naively applying the same procedure used for the unbiased tracers to a halo mock distribution does not provide success- full results, suggesting that the Vdn model requires to be reconsidered when dealing with biased samples. Thus, we test two alternative exten- sions of the model and find that two scaling relations exist: both the Dark Matter void radii and the underlying Dark Matter density contrast scale with the halo-defined void radii. We use these findings to develop a semi-analytical model which gives promising results.

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

Primordial black holes from multi-field inflation

In this master's thesis, the formation of Primordial Black Holes (PBHs) in the context of multi-field inflation is studied. In these models, the interaction of isocurvature and curvature perturbations can lead to a significant enhancement of the latter, and to the subsequent production of PBHs. Depending on their mass, these can account for a significant fraction (or, in some cases, the entirety) of the universe's Dark Matter content. After studying the theoretical framework of generic N-field inflationary models, the focus is restricted to the two-field case, for which a few concrete realisations are analysed. A numerical code (written in Wolfram Mathematica) is developed to make quantitative predictions for the main inflationary observables, notably the scalar power spectra. Parallelly, the production of PBHs due to the dynamics of 2-field inflation is examined: their mass, as well as the fraction of Dark Matter they represent, is calculated for the models considered previously.

Amplification mechanisms of inflationary spectra in induced gravity

Inflation is the primordial stage of accelerated expansion of the Universe which solves the issues of the initial conditions of a decelerating Universe (horizon, flatness and entropy problems). Moreover, it is supposed that quantum fluctuations originated during the first moments after the Big Bang gave rise to the formation of galaxies and other structures of the Universe when inflation ends. Among these structures also primordial black holes (PBHs) may have been generated. The interest in PBHs relies on their possible connection with dark matter: they could constitute a portion or even the whole dark matter content of our Universe.\\ In this work we consider inflation in the Induced Gravity (IR) context and study possible mechanisms of amplification of the curvature perturbations generated during the cosmic acceleration. In particular we consider the possibility of a period of Constant Roll (CR). Starting from the previous work of Starobinsky et al. Our final purpose is to analyse the power spectrum of the scalar perturbations and to find in which conditions there is an enhancement of the power spectrum possibly leading to PBHs formation.

Studio degli eventi di muone nell'esperimento XENONnT

La Materia Oscura (Dark Matter, DM) deve il suo nome al fatto che non interagisce elettromagneticamente, ma solo gravitazionalmente e debolmente (proprietà che ne complica particolarmente la rivelazione). Molti sforzi, sia sperimentali che teorici, sono stati dedicati alla sua ricerca a tal punto da essere considerata uno dei più grandi misteri della fisica moderna. I candidati più promettenti a costituire la materia oscura sono da ricercarsi in teorie oltre il Modello Standard e fra essi figurano le WIMPs (Weakly Interacting Massive Particles). Poiché le particelle di DM sono caratterizzate da sezioni d’urto molto piccole, per poterle osservare si necessita di grandi rivelatori, ultrapuri e situati in ambienti a bassa radioattività. XENONnT è attualmente uno degli esperimenti più sensibili al mondo per la ricerca diretta di WIMPs, grazie all’utilizzo di una camera di proiezione temporale (TPC) a doppia fase (liquido-gas) che presenta una massa bersaglio di 5.9 t di xenon liquido (LXe). XENONnT mira a rilevare lo scattering di WIMPs contro i nuclei bersaglio, sfruttando i segnali osservabili di luce e carica prodotti dai rinculi nel LXe. Per raggiungere tale risultato sono fondamentali il sistema di veto di neutroni (Neutron Veto, NV) ed il sistema di veto di muoni (Muon Veto, MV): il NV ha lo scopo di rivelare i neutroni radiogenici provenienti dai materiali dell’esperimento mentre il MV (già impiegato in XENON1T) ha la finalità di eliminare il rumore di fondo dovuto ai raggi cosmici. L'obiettivo di questa tesi è lo studio degli eventi muone osservati dal MV ed il NV di XENONnT; misurarne il rate e verificare la congruenza dei risultati ottenuti con quanto atteso da simulazioni Monte Carlo. Questa prova ha costituito un'ulteriore verifica del corretto funzionamento dei sistemi di veto, ed ha permesso di approfondire la conoscenza sulla risposta del Neutron Veto al passaggio di un muone all’interno del suo volume attivo.

Gravitational microscope in MACSJ1206: unveiling galaxy cluster member mass distribution

Gravitational lensing is a powerful tool to investigate the properties of the distribution of matter, be it barionic or dark. In this work we take advantage of Strong Gravitational Lensing to infer the properties of one of the galaxy-scale substructures that makes up the cluster MACSJ1206. It is relatively easy to model the morphology of the visible components of a galaxy, while the morphology of the dark matter distribution cannot be so easily constrained. Being sensitive to the whole mass, strong lensing provides a way to probe DM distribution, and this is the reason why it is the best tool to study the substructure. The goal of this work consists of performing an analysis of the substructure previously mentioned, an early type galaxy (ETG), by analyzing the highly magnified Einstein ring around it, in order to put stringent constraints on its matter distribution, that, for an ETG, is commonly well described by an isothermal profilele. This turns out to be interesting for three main different reasons. It is well known that galaxies in clusters are subject to interaction processes, both dynamic and hydrodynamic, that can significantly modify the distribution of matter within them. Therefore, finding a different profile from the one usually expected could be a sign that the galaxy has undergone processes that have changed its structure. Studying the mass distribution also means studying the dark matter component, which not only still presents great questions today, but which is also not obviously distributed in the same way as in an isolated galaxy. What emerges from the analysis is that the total mass distribution of the galaxy under examination turns out to have a slope much steeper than the isothermal usually expected.

Constraints on the duty cycle of the remnant radio galaxy NGC 6086 in the galaxy group Abell 2162

Galaxy clusters and groups are the most massive bounded structures and the knots of the large-scale structure of the Universe. These structures reside in dark matter haloes, hosting tens to hundreds of galaxies and they are filled with hot and rarefied gas. Radio Galaxies are a peculiar class of galaxies with a luminosity in the radio band up to 10^46 erg/s between 10 MHz and 100 GHz. These galaxies are a subclass of AGN in which there is accretion on the Super Massive Black Hole. The accretion generates jets of relativistic particles and magnetic fields which lose energy through synchrotron radiation, best observable at radio frequencies. The study of the spectral ageing of the AGN plasma is fundamental to understand its evolution, interaction with the environment and to constrain the AGN duty cycle. n this thesis, we have investigated the duty cycle of the nearby remnant radio galaxy NGC 6086, located in the centre of the galaxy group Abell 2162. We have made major steps forward thanks to the new high-sensitivity interferometers in the low-frequency radio band. We have detected for the first time three filaments of emission and a second couple of lobes. We have performed an integrated and resolved analysis on the previously known inner lobes, the new filaments and the older outer lobes. We have performed an age estimate of the two pairs of lobes to give constraints on the duty cycle of the source and an estimate of its active time.

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.

Stellar mass loss from the Fornax dwarf spheroidal galaxy

Dwarf galaxies often experience gravitational interactions from more massive companions. These interactions can deform galaxies, turn star formation on or off, or give rise to mass loss phenomena. In this thesis work we propose to study, through N-body simulations, the stellar mass loss suffered by the dwarf spheroid galaxy (dSph) Fornax orbiting in the Milky Way gravitational potential. Which is a key phenomenon to explain the mass budget problem: the Fornax globular clusters together have a stellar mass comparable to that of Fornax itself. If we look at the stellar populations which they are made of and we apply the scenarios of stellar population formation we find that, originally, they must have been >= 5 times more massive. For this reason, they must have lost or ejected stars through dynamic interactions. However, as presented in Larsen et al (2012), field stars alone are not sufficient to explain this scenario. We may assume that some of those stars fell into Fornax, and later were stripped by Milky Way. In order to study this solution we built several illustrative single component simulations, with a tabulated density model using the P07ecc orbit studied from Battaglia et al (2015). To divide the single component into stellar and dark matter components we have defined a posterior the probability function P(E), where E is the initial energy distribution of the particles. By associating each particle with a fraction of stellar mass and dark matter. In this way we built stellar density profiles without repeating simulations. We applied the method to Fornax using the profile density tables obtained in Pascale et al (2018) as observational constraints and to build the model. The results confirm the results previously obtained with less flexible models by Battaglia et al (2015). They show a stellar mass loss < 4% within 1.6 kpc and negligible within 3 kpc, too small to solve the mass budget problem.