Probing the innermost regions of AGN via time-resolved spectral analysis

The dynamics and geometry of the material inflowing and outflowing close to the supermassive black hole in active galactic nuclei are still uncertain. X-rays are the most suitable way to study the AGN innermost regions because of the Fe Kα emission line, a proxy of accretion, and Fe absorption lines produced by outflows. Winds are typically classified as Warm Absorbers (slow and mildly ionized) and Ultra Fast Outflows (fast and highly ionized). Transient Obscurers -optically thick winds that produce strong spectral hardening in X-rays, lasting from days to months- have been observed recently. Emission and absorption features vary on time-scales from hours to years, probing phenomena at different distances from the SMBH. In this work, we use time-resolved spectral analysis to investigate the accretion and ejection flows, to characterize them individually and search for correlations. We analyzed XMM-Newtomn data of a set of the brightest Seyfert 1 galaxies that went through an obscuration event: NGC 3783, NGC 3227, NGC 5548, and NGC 985. Our aim is to search for emission/absorption lines in short-duration spectra (∼ 10ks), to explore regions as close as the SMBH as the statistics allows for, and possibly catch transient phenomena. First we run a blind search to detect emission/absorption features, then we analyze their evolution with Residual Maps: we visualize simultaneously positive and negative residuals from the continuum in the time-energy plane, looking for patterns and relative time-scales. In NGC 3783 we were able to ascribe variations of the Fe Kα emission line to absorptions at the same energy due to clumps in the obscurer, whose presence is detected at >3σ, and to determine the size of the clumps. In NGC 3227 we detected a wind at ∼ 0.2c at ∼ 2σ, briefly appearing during an obscuration event.

Towards a full cosmological exploitation of cosmic voids

A stately fraction of the Universe volume is dominated by almost empty space. Alongside the luminous filamentary structures that make it up, there are vast and smooth regions that have remained outside the Cosmology spotlight during the past decades: cosmic voids. Although essentially devoid of matter, voids enclose fundamental information about the cosmological framework and have gradually become an effective and competitive cosmological probe. In this Thesis work we present fundamental results about the cosmological exploitation of voids. We focused on the number density of voids as a function of their radius, known as void size function, developing an effective pipeline for its cosmological usage. We proposed a new parametrisation of the most used theoretical void size function to model voids identified in the distribution of biased tracers (i.e. dark matter haloes, galaxies and galaxy clusters), a step of fundamental importance to extend the analysis to real data surveys. We then applied our built methodology to study voids in alternative cosmological scenarios. Firstly we exploited voids with the aim of breaking the degeneracies between cosmological scenarios characterised by modified gravity and the inclusion of massive neutrinos. Secondly we analysed voids in the perspective of the Euclid survey, focusing on the void abundance constraining power on dynamical dark energy models with massive neutrinos. Moreover we explored other void statistics like void profiles and clustering (i.e. the void-galaxy and the void-void correlation), providing cosmological forecasts for the Euclid mission. We finally focused on the probe combination, highlighting the incredible potential of the joint analysis of multiple void statistics and of the combination of the void size function with different cosmological probes. Our results show the fundamental role of the void analysis in constraining the fundamental parameters of the cosmological model and pave the way for future studies on this topic.

An efficient Jeans modelling of axisymmetric galaxies with multiple stellar components

Dynamical models of stellar systems represent a powerful tool to study their internal structure and dynamics, to interpret the observed morphological and kinematical fields, and also to support numerical simulations of their evolution. We present a method especially designed to build axisymmetric Jeans models of galaxies, assumed as stationary and collisionless stellar systems. The aim is the development of a rigorous and flexible modelling procedure of multicomponent galaxies, composed of different stellar and dark matter distributions, and a central supermassive black hole. The stellar components, in particular, are intended to represent different galaxy structures, such as discs, bulges, halos, and can then have different structural (density profile, flattening, mass, scale-length), dynamical (rotation, velocity dispersion anisotropy), and population (age, metallicity, initial mass function, mass-to-light ratio) properties. The theoretical framework supporting the modelling procedure is presented, with the introduction of a suitable nomenclature, and its numerical implementation is discussed, with particular reference to the numerical code JASMINE2, developed for this purpose. We propose an approach for efficiently scaling the contributions in mass, luminosity, and rotational support, of the different matter components, allowing for fast and flexible explorations of the model parameter space. We also offer different methods of the computation of the gravitational potentials associated of the density components, especially convenient for their easier numerical tractability. A few galaxy models are studied, showing internal, and projected, structural and dynamical properties of multicomponent galaxies, with a focus on axisymmetric early-type galaxies with complex kinematical morphologies. The application of galaxy models to the study of initial conditions for hydro-dynamical and $N$-body simulations of galaxy evolution is also addressed, allowing in particular to investigate the large number of interesting combinations of the parameters which determine the structure and dynamics of complex multicomponent stellar systems.

Unveiling the unknown of cool stars with high-resolution spectroscopy

Cool giant and supergiant stars are among the brightest populations in any stellar system and they are easily observable out to large distances, especially at infrared wavelengths. These stars also dominate the integrated light of star clusters in a wide range of ages, making them powerful tracers of stellar populations in more distant galaxies. High-resolution near-IR spectroscopy is a key tool for quantitatively investigating their kinematic, evolutionary and chemical properties. However, the systematic exploration and calibration of the NIR spectral diagnostics to study these cool stellar populations based on high-resolution spectroscopy is still in its pioneering stage. Any effort to make progress in the field is innovative and of impact on stellar archaeology and stellar evolution. This PhD project takes the challenge of exploring that new parameter space and characterizing the physical properties, the chemical content and the kinematics of cool giants and supergiants in selected disc fields and clusters of our Galaxy, with the ultimate goal of tracing their past and recent star formation and chemical enrichment history. By using optical HARPS-N and near-infrared GIANO-B high-resolution stellar spectra in the context of the large program SPA-Stellar Population Astrophysics: the detailed, age-resolved chemistry of the Milky Way disk” (PI L. Origlia), an extensive study of Arcturus, a standard calibrator for red giant stars, has been performed. New diagnostics of stellar parameters as well as optimal linelists for chemical analysis have been provided. Then, such diagnostics have been used to determine evolutionary properties, detailed chemical abundances of almost 30 different elements and mixing processes of a homogeneous sample of red supergiant stars in the Perseus complex.

A hydrodynamical view on the mass accretion onto massive halos

This Thesis presents the results of my work on how galaxy clusters form by the accretion of sub-clumps and diffuse materials, and how the accreted energy is distributed in the X-ray emitting plasma. Indeed, on scales larger than tens of millions of light years, the Universe is self-organised by gravity into a spiderweb, the Cosmic Web. Galaxy clusters are the knots of this Cosmic Web, but a strong definition of filaments (which link different knots) and their physical proprieties, is still uncertain. Even if this pattern was determined by studying the spatial distribution of galaxies in the optical band, recently, also in the X-rays probes of filamentary structures around galaxy clusters were obtained. Therefore, given these observational facilities, the galaxy clusters’ outskirts are the best candidate regions to detect filaments and study their physical characteristics. However, from X-rays observations, we have only a few detections of cosmic filaments to date. On the other hand, it is crucial to understand how the accreted energy is dissipated in the baryon content of galaxy clusters and groups. Indeed, it is well known that in the central region of galaxy clusters and groups, the baryon fraction increases with the halo mass. On the outer region, the lack of X-rays constraints influences our understanding of the evolution of baryons in the halos volume. The standard assumption of “closed-box” system, for which the baryon fraction should approach the cosmological ratio Omega_bar/Omega_m, for galaxy clusters and groups seems to be too strong, especially for less massive objects. Moreover, a complete redshift evolution of baryons in galaxy clusters and groups is still missing.

Exploring the internal kinematics of Galactic Globular Cluster cores

Globular clusters (GCs) are traditionally described as simple quasi-relaxed non-rotating stellar systems, characterized by spherical symmetry and isotropy in velocity space. However, recent studies have shown deviations from isotropic velocity distributions and significant internal rotation in many GCs, suggesting that their internal structure and kinematics are more complex than previously thought. The aim of this thesis is to investigate the internal kinematics of Galactic Globular Clusters (GGCs) as part of the Multi-Instrument Kinematic Survey (MIKiS), which exploits the capabilities of different ESO-VLT spectrographs to obtain comprehensive velocity dispersion (VD) and rotation profiles of GGCs. Moreover, this thesis has the particular goal of unraveling the kinematics of GC cores, which are still largely unexplored, by taking advantage of the exceptional spatial resolution of the adaptive-optics assisted integral-field spectrograph MUSE/NFM. The thesis presents a thorough kinematic study of three GGCs NGC 1904, NGC 6440, and NGC 6569. By combining the data sets acquired with four different spectrographs, we obtained the radial velocity (RV) of more than 1000 individual stars in each cluster, sampling from the innermost to the outermost regions. This allowed us to obtain the entire VD profile of each cluster and exclude the presence of an intermediate-mass black hole in the core of NGC 1904, at odds with previous findings obtained from integrated-light spectra. The studies also revealed signatures of internal rotation in each of the GCs studied. These results, supported by those of N-body simulations, prove that GCs were born with a significant initial rotation that they gradually lost through internal two-body relaxation and angular momentum loss carried away by escaping stars. Furthermore, we derived the structural parameters of NGC 6440 and NGC 6569, obtaining a comprehensive overview of the internal kinematics and structure of these GCs, which is necessary to properly reconstruct the evolutionary history of these systems.

Cross-population radio galaxies as a test of the jet-accretion relationship

Radio galaxies (RGs) are extremely relevant in addressing important unknowns concerning the interaction among black hole accretion, radio jets, and the environment. In the classical scheme, their accretion rate and ejection of relativistic jets are directly linked: efficient accretion (HERG) is associated with powerful edge-brightened jets (FRIIs); inefficient accretion (LERG) is associated with weak edge-darkened jets (FRIs). The observation of RGs with an inefficient engine associated with edge-brightened radio emission (FRII-LERGs) broke this scheme. FRII-LERGs constitute a suitable population to explore how accretion and ejection are linked and evaluate the environment's role in shaping jets. To this aim, we performed a multiwavelength study of different RGs catalogs spanning from Jy to mJy flux densities. At first, we investigated the X-ray properties of a sample of 51 FRIIs belonging to the 3CR catalog at z<0.3. Two hypotheses were invoked to explain FRII-LERGs behavior: evolution from classical FRIIs; the role of the environment. Next, we explored the mJy sky by studying the optical-radio properties of hundreds of RGs at z<0.15 (Best & Heckman 2012 sample). FRII-LERGs appear more similar to the old FRI-LERGs than to the young FRII-HERGs. These results point towards an evolutive scenario, however, nuclear time scale changes, star population aging, and kpc-Mpc radio structure modification do not agree. The role of the Mpc environment was then investigated. The Wen et al. 2015 galaxy clusters sample, built exploiting the SDSS survey, allowed us to explore the habitat of 7219 RGs at z<0.3. Most RGs are found to live in outside clusters. For these sources, differences among RG classes are still present. Thus, the environment is not the key parameter, and the possibility of intrinsic differences was reconsidered: we speculated that different black hole properties (spin and magnetic field at its horizon) could determine the observed spread in jet luminosity.

New indicators of star cluster dynamical evolution

This thesis presents a study of globular clusters (GCs), based on analysis of Monte Carlo simulations of globular clusters (GCs) with the aim to define new empirical parameters measurable from observations and able to trace the different phases of their dynamical evolution history. During their long term dynamical evolution, due to mass segregation and and dynamical friction, massive stars transfer kinetic energy to lower-mass objects, causing them to sink toward the cluster center. This continuous transfer of kinetic energy from the core to the outskirts triggers the runaway contraction of the core, known as "core collapse" (CC), followed by episodes of expansion and contraction called gravothermal oscillations. Clearly, such an internal dynamical evolution corresponds to significant variations also of the structure of the system. Determining the dynamical age of a cluster can be challenging as it depends on various internal and external properties. The traditional classification of GCs as CC or post-CC systems relies on detecting a steep power-law cusp in the central density profile, which may not always be reliable due to post-CC oscillations or other processes. In this thesis, based on the normalized cumulative radial distribution (nCRD) within a fraction of the half-mass radius is analyzed, and three diagnostics (A5, P5, and S2.5) are defined. These diagnostics show sensitivity to dynamical evolution and can distinguish pre-CC clusters from post-CC clusters.The analysis performed using multiple simulations with different initial conditions, including varying binary fractions and the presence of dark remnants showed the time variations of the diagnostics follow distinct patterns depending on the binary fraction and the retention or ejection of black holes. This analysis is extended to a larger set of simulations matching the observed properties of Galactic GCs, and the parameters show a potential to distinguish the dynamical stages of the observed clusters as well.

Unveiling the expansion history of the universe with cosmic chronometers and gravitational waves

This Thesis explores two novel and independent cosmological probes, Cosmic Chronometers (CCs) and Gravitational Waves (GWs), to measure the expansion history of the Universe. CCs provide direct and cosmology-independent measurements of the Hubble parameter H(z) up to z∼2. In parallel, GWs provide a direct measurement of the luminosity distance without requiring additional calibration, thus yielding a direct measurement of the Hubble constant H0=H(z=0). This Thesis extends the methodologies of both of these probes to maximize their scientific yield. This is achieved by accounting for the interplay of cosmological and astrophysical parameters to derive them jointly, study possible degeneracies, and eventually minimize potential systematic effects. As a legacy value, this work also provides interesting insights into galaxy evolution and compact binary population properties. The first part presents a detailed study of intermediate-redshift passive galaxies as CCs, with a focus on the selection process and the study of their stellar population properties using specific spectral features. From their differential aging, we derive a new measurement of the Hubble parameter H(z) and thoroughly assess potential systematics. In the second part, we develop a novel methodology and pipeline to obtain joint cosmological and astrophysical population constraints using GWs in combination with galaxy catalogs. This is applied to GW170817 to obtain a measurement of H0. We then perform realistic forecasts to predict joint cosmological and astrophysical constraints from black hole binary mergers for upcoming gravitational wave observatories and galaxy surveys. Using these two probes we provide an independent reconstruction of H(z) with direct measurements of H0 from GWs and H(z) up to z∼2 from CCs and demonstrate that they can be powerful independent probes to unveil the expansion history of the Universe.

Cosmology from probe combination: synergies between the number counts of galaxy clusters and cosmic voids

This Thesis work concerns the complementary study of the abundance of galaxy clusters and cosmic voids identified in cosmological simulations, at different redshifts. In particular, we focus our analyses on the combination of the cosmological constraints derived from these probes, which can be considered statistically independent, given the different aspects of Universe density field they map. Indeed, we aim at showing the orthogonality of the derived cosmological constraints and the resulting impressive power of the combination of these probes. To perform this combination we apply three newly implemented algorithms that allow us to combine independent probes. These algorithms represent a flexible and user-friendly tool to perform different techniques for probe combination and are implemented within the environment provided by the large set of free software C++/Python CosmoBolognaLib. All the new implemented codes provide simple and flexible tools that will be soon applied to the data coming from currently available and next-generation wide-field surveys to perform powerful combined cosmological analyses.

Investigating the origin of radio emission from the intra-cluster bridge A399-A401

Radio bridges are a newly observed class of phenomena, they are large scale filamentary structures that span between pairs of clusters in a pre-merger phase. In my Thesis, starting from the observed synchrotron emission in the bridge A399-A401, I apply a Fermi I shock re-acceleration model under different conditions to the bridge. The purpose of my work is to check the likelihood of the above mentioned model. In particular the Inverse Compton emission that the model predicts is above the current observational constraints on Inverse Compton for the bridge A399-A401, this means that a Fermi I re-acceleration model is unlikely to describe the emission from the bridge, while a Fermi II turbulent re-acceleration model could still explain the origin of the emission in the bridge.

Georges Lemaître: a visionary contribute to the quest of the world model

The aim of this dissertation is to present the sequence of events which brought the scientific community of the early 20th century to conceive an expanding Universe born from a single origin. Among the facts here reported, some are well-known, some others instead are little-known backstories, not so easy neither to obtain nor to trust. Indeed, several matters shown in this thesis, now as then, create a battleground among scientists. Amid the numerous personalities whose contributions are discussed in this work, the main protagonist is surely Georges Lemaître, who managed to combine – without overlapping – his being both a priest and a scientist. The first chapter is dedicated to his biography, from his childhood in Belgium, to his early adulthood between England and the USA, to his success in the scientific community. The second and the third chapter explain how the race to the understanding of a Universe which not only expands, but also originated from a singularity, developed. The Belgian priest’s discoveries, as shown, were challenged by other important scientists, who, in several cases, Lemaître had a friendly relationship with. As a consequence, the fourth and final chapter deals with the multiple relations that the priest managed to build, thanks to his politeness and kindness. Moreover, it is also covered Lemaître’s personal connection with the Church and religion, without forgetting the personalities that influenced him – above all, Saint Thomas Aquinas. As a conclusion to this thesis, two appendices gather not only a summary of Lemaître’s works which are not already described in the chapters, but also the biographies of all the characters presented in this dissertation.

Exploiting the clustering of cosmic voids as a novel cosmological probe

The investigations of the large-scale structure of our Universe provide us with extremely powerful tools to shed light on some of the open issues of the currently accepted Standard Cosmological Model. Until recently, constraining the cosmological parameters from cosmic voids was almost infeasible, because the amount of data in void catalogues was not enough to ensure statistically relevant samples. The increasingly wide and deep fields in present and upcoming surveys have made the cosmic voids become promising probes, despite the fact that we are not yet provided with a unique and generally accepted definition for them. In this Thesis we address the two-point statistics of cosmic voids, in the very first attempt to model its features with cosmological purposes. To this end, we implement an improved version of the void power spectrum presented by Chan et al. (2014). We have been able to build up an exceptionally robust method to tackle with the void clustering statistics, by proposing a functional form that is entirely based on first principles. We extract our data from a suite of high-resolution N-body simulations both in the LCDM and alternative modified gravity scenarios. To accurately compare the data to the theory, we calibrate the model by accounting for a free parameter in the void radius that enters the theory of void exclusion. We then constrain the cosmological parameters by means of a Bayesian analysis. As far as the modified gravity effects are limited, our model is a reliable method to constrain the main LCDM parameters. By contrast, it cannot be used to model the void clustering in the presence of stronger modification of gravity. In future works, we will further develop our analysis on the void clustering statistics, by testing our model on large and high-resolution simulations and on real data, also addressing the void clustering in the halo distribution. Finally, we also plan to combine these constraints with those of other cosmological probes.

Analytic models of self-gravitating rotating gaseous tori with central black hole

In this thesis project, I present stationary models of rotating fluids with toroidal distributions that can be used to represent the active galactic nuclei (AGN) central obscurers, i.e. molecular tori (Combes et al., 2019), as well as geometrically thick accretion discs, like ADAF discs (Narayan and Yi, 1995) or Polish doughnuts (Abramowicz, 2005). In particular, I study stationary rotating systems with a more general baroclinic distribution (with a vertical gradient of the angular velocity), which are often more realistic and less studied, due to their complexity, than the barotropic ones (with cylindrical rotation), which are easier to construct. In the thesis, I compute analytically the main intrinsic and projected properties of the power-law tori based on the potential-density pairs of Ciotti and Bertin (2005). I study the density distribution and the resulting gravitational potential for different values of α, in the range 2 < α < 5. For the same models, I compute the surface density of the systems when seen face-on and edge-on. I then apply the stationary Euler equations to obtain rotational velocity and temperature distributions of the self-gravitating models in the absence of an external gravitational potential. In the thesis I also consider the power-law tori with the presence of a central black hole in addition to the gas self-gravity, and solving analytically the stationary Euler equations, I compute how the properties of the system are modified by the black hole and how they vary as a function of the black hole mass. Finally, applying the Solberg-Høiland criterion, I show that these baroclinic stationary models are linearly stable in the absence of the black hole. In the presence of the black hole I derive the analytical condition for stability, which depends on α and on the black hole mass. I also study the stability of the tori in the hypothesis that they are weakly magnetized, finding that they are always unstable to this instability.

Modellizzazione dell'evento di microlensing OB171226

Un'anomalia della curva di luce di un evento di microlensing può indicare la presenza di una lente composta da un sistema stella-pianeta. Quindi questo effetto si può utilizzare come tecnica per la ricerca di esopianeti. Partendo dallo studio teorico dei fondamenti del lensing gravitazionale e del microlensing, la tesi di laurea presentata include il lavoro di analisi svolto durante il tirocinio curriculare. Lo scopo di questo lavoro è stato sviluppare un modello per la curva di luce dell'evento OB171226 rilevato nel 2017 dai dati di quattro telescopi terrestri (OGLE e KMTNet) ed uno spaziale (Spitzer), grazie al quale si è visto come questo probabilmente non sia un sistema di tipo planetario ma composto da due stelle compagne di diversa massa. Questa evento era già stato studiato ma lo sviluppo ulteriore fatto è stato modellizzarlo considerando anche il moto orbitale dei due corpi e l'effetto di parallasse da satellite per i dati spaziali.