The gas mass fraction and the dynamical state in x-ray luminous clusters of galaxies at low redshift

Gli ammassi di galassie sono le strutture gravitazionalmente legate con le più profonde buche di potenziale, pertanto è previsto che questi contengano una frazione di barioni non molto diversa da quella cosmologica. Con l’introduzione di modelli sempre più accurati di fisica barionica all’interno di simulazioni idrodinamiche è stato possibile predire la percentuale cosmica di barioni presente negli ammassi di galassie. Unendo questi modelli previsionali con misure della frazione di gas in ammassi e informazioni sulla densità di barioni dell’Universo si può ottenere una stima della densità di materia cosmica Ωm. L'obiettivo di questo lavoro di Tesi è la stima di Ωm a partire dalla frazione di gas osservata in questi sistemi. Questo lavoro era stato già fatto in precedenza, ma tenendo in considerazione solo gli ammassi più massivi e dinamicamente rilassati. Usando parametri che caratterizzano la morfologia della distribuzione di brillanza superficiale nei raggi X, abbiamo classificato i nostri oggetti come rilassati o disturbati, laddove presentassero evidenze di recenti attività di interazione. Abbiamo dunque valutato l’impatto degli oggetti disturbati sulla stima del parametro cosmologico Ωm, computando il Chi2 tra la frazione di massa barionica nell’Universo e quella da noi ricavata. Infine abbiamo investigato una relazione tra il valore della frazione di gas degli ammassi rilassati e quello dei disturbati, in modo da correggere quindi questi ultimi, riportandoli nei dintorni del valore medio per i rilassati e usarli per ampliare il campione e porre un vincolo più stringente su Ωm. Anche con il limitato campione a nostra disposizione, è stato possibile porre un vincolo più stretto su Ωm, utilizzando un maggior numero di oggetti e riducendo così l’errore statistico.

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.

High resolution study of narrow tailed radio galaxies in clusters

Clusters of galaxies are the most massive and large gravitationally bounded systems in the whole Universe. Their study is of fundamental importance to constrain cosmological parameters and to obtain informations regarding various kind of emission in different wavebands. In particular, in the radio domain, beside the diffuse emission, the study is focused on the radio galaxies emission. Radio galaxies in clusters can have peculiar morphology, since they interact with the intracluster medium (ICM) in which they are embedded. Particularly, in this thesis we focused our attention on the so-called Narrow-Angle Tailed radio galaxies (NAT), which present radio jets that are bent at extreme angle, up to 90 degrees, from their original orientation. Some NAT show a narrow extended structure and the two radio tails are not resolved even with high resolution radio observations. An example is provided by the source IC310, in the Perseus Cluster, whose structure has been recently interpreted as due to Doppler boosting effects of a relativistic jet oriented at a small angle with respect to the line of sight. If the structure is due to relativistic effects, this implies that the jets are relativistic at about 400 kpc from the core, but this is in contrast with unified models, which predict that for low-power radio source (NAT are classified as FRI radio galaxies) the jets decelerate to sub-relativistic speed within a few kpc from the core. To investigate this scientific topic, in this thesis we have analyzed the innermost structure of a sample of eleven radio galaxies showing a very narrow NAT structure. We can conclude that the structure of these radio galaxies is different from that of IC310. These radio galaxies are indeed strongly influenced by environmental effects and are similar to classical NAT sources.

Non-thermal properties of two galaxy filament candidates from JVLA observations

Large-scale structures can be considered an interesting and useful "laboratory" to better investigate the Universe; in particular the filaments connecting clusters and superclusters of galaxies can be a powerful tool for this intent, since they are not virialised systems yet. The large structures in the Universe have been studied in different bands, in particular the present work takes into consideration the emission in the radio band. In the last years both compact and diffuse radio emission have been detected, revealing to be associated to single objects and clusters of galaxies respectively. The detection of these sources is important, because the radiation process is the synchrotron emission, which in turn is linked to the presence of a magnetic field: therefore studying these radio sources can help in investigating the magnetic field which permeates different portions of space. Furthermore, radio emission in optical filaments have been detected recently, opening new chances to further improve the understanding of structure formation. Filaments can be seen as the net which links clusters and superclusters. This work was made with the aim of investigating non-thermal properties in low-density regions, looking for possible filaments associated to the diffuse emission. The analysed sources are 0917+75, which is located at a redshift z = 0.125, and the double cluster system A399-A401, positioned at z = 0.071806 and z = 0.073664 respectively. Data were taken from VLA/JVLA observations, and reduced and calibrated with the package AIPS, following the standard procedure. Isocountour and polarisation maps were yielded, allowing to derive the main physical properties. Unfortunately, because of a low quality data for A399-A401, it was not possible to see any radio halo or bridge.

Heavily obscured AGN: an ideal laboratory to study the early co-evolution of galaxies and black holes

Obscured AGN are a crucial ingredient to understand the full growth history of super massive black holes and the coevolution with their host galaxies, since they constitute the bulk of the BH accretion. In the distant Universe, many of them are hosted by submillimeter galaxies (SMGs), characterized by a high production of stars and a very fast consumption of gas. Therefore, the analysis of this class of objects is fundamental to investigate the role of the ISM in the early coevolution of galaxies and black holes. We present a multiwavelength study of a sample of six obscured X-ray selected AGN at z>2.5 in the CDF-S, detected in the far-IR/submm bands. We performed the X-ray spectral analysis based on the 7Ms Chandra dataset, which provides the best X-ray spectral information currently available for distant AGN. We were able to place constraints on the obscuring column densities and the intrinsic luminosities of our targets. Moreover, we built up the UV to FIR spectral energy distributions (SEDs) by combining the broad-band photometry from CANDELS and the Herschel catalogs, and analyzed them by means of an SED decomposition technique. Therefore, we derived important physical parameters of both the host galaxy and the AGN. In addition, we obtained, through an empirical calibration, the gas mass in the host galaxy and assessed the galaxy sizes in order to estimate the column density associated with the host ISM. The comparison of the ISM column densities with the values measured from the X-ray spectral analysis pointed out that the contribution of the host ISM to the obscuration of the AGN emission can be substantial, ranging from ~10% up to ~100% of the value derived from the X-ray spectra. The absorption may occur at different physical scales in these sources and, in particular, the medium in the host galaxy is an ingredient that should be taken into account, since it may have a relevant role in driving the early co-evolution of galaxies with their black holes.

Constraining mass and shape of galaxy clusters through large scale structures

The mass estimation of galaxy clusters is a crucial point for modern cosmology, and can be obtained by several different techniques. In this work we discuss a new method to measure the mass of galaxy clusters connecting the gravitational potential of the cluster with the kinematical properties of its surroundings. We explore the dynamics of the structures located in the region outside virialized cluster, We identify groups of galaxies, as sheets or filaments, in the cluster outer region, and model how the cluster gravitational potential perturbs the motion of these structures from the Hubble fow. This identification is done in the redshift space where we look for overdensities with a filamentary shape. Then we use a radial mean velocity profile that has been found as a quite universal trend in simulations, and we fit the radial infall velocity profile of the overdensities found. The method has been tested on several cluster-size haloes from cosmological N-body simulations giving results in very good agreement with the true values of virial masses of the haloes and orientation of the sheets. We then applied the method to the Coma cluster and even in this case we found a good correspondence with previous. It is possible to notice a mass discrepancy between sheets with different alignments respect to the center of the cluster. This difference can be used to reproduce the shape of the cluster, and to demonstrate that the spherical symmetry is not always a valid assumption. In fact, if the cluster is not spherical, sheets oriented along different axes should feel a slightly different gravitational potential, and so give different masses as result of the analysis described before. Even this estimation has been tested on cosmological simulations and then applied to Coma, showing the actual non-sphericity of this cluster.

Angular momentum of S0 galaxies

The aim of this thesis is to study the angular momentum of a sample of S0 galaxies. In the quest to understand whether the formation of S0 galaxies is more closely linked to that of ellipticals or that of spirals, our goal is to compare the amount of their specific angular momentum as a function of stellar mass with respect to spirals. Through kinematic comparison between these different classes of galaxies we aim to understand if a scenario of passive evolution, in which the galaxy’s gas is consumed and the star formation is quenched, can be considered as plausible mechanism to explain the transformation from spirals to S0s. In order to derive the structural and photometric parameters of galaxy sub-components we performed a bulge-disc decomposition of optical images using GALFIT. The stellar kinematic of the galaxies was measured using integral field spectroscopic data from CALIFA survey. The development of new original software, based on a Monte Carlo Markov Chain algorithm, allowed us to obtain the values of the line of sight velocity and velocity dispersion of disc and bulge components. The result that we obtained is that S0 discs have a distribution of stellar specific angular momentum that is in full agreement with that of spiral discs, so the mechanism of simple fading can be considered as one of the most important for transformation from spirals to S0s.

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.

Parallelization of the algorithm WHAM with NVIDIA CUDA.

The aim of my thesis is to parallelize the Weighting Histogram Analysis Method (WHAM), which is a popular algorithm used to calculate the Free Energy of a molucular system in Molecular Dynamics simulations. WHAM works in post processing in cooperation with another algorithm called Umbrella Sampling. Umbrella Sampling has the purpose to add a biasing in the potential energy of the system in order to force the system to sample a specific region in the configurational space. Several N independent simulations are performed in order to sample all the region of interest. Subsequently, the WHAM algorithm is used to estimate the original system energy starting from the N atomic trajectories. The parallelization of WHAM has been performed through CUDA, a language that allows to work in GPUs of NVIDIA graphic cards, which have a parallel achitecture. The parallel implementation may sensibly speed up the WHAM execution compared to previous serial CPU imlementations. However, the WHAM CPU code presents some temporal criticalities to very high numbers of interactions. The algorithm has been written in C++ and executed in UNIX systems provided with NVIDIA graphic cards. The results were satisfying obtaining an increase of performances when the model was executed on graphics cards with compute capability greater. Nonetheless, the GPUs used to test the algorithm is quite old and not designated for scientific calculations. It is likely that a further performance increase will be obtained if the algorithm would be executed in clusters of GPU at high level of computational efficiency. The thesis is organized in the following way: I will first describe the mathematical formulation of Umbrella Sampling and WHAM algorithm with their apllications in the study of ionic channels and in Molecular Docking (Chapter 1); then, I will present the CUDA architectures used to implement the model (Chapter 2); and finally, the results obtained on model systems will be presented (Chapter 3).

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.

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 (4galaxies transform from newly-formed objects to more complex structures. In a few hundred Myrs, galaxies increase their stellar mass, change their morphological type and progressively become more obscured due to increased dust attenuation of the UV light. Therefore, determining physical parameters regarding dust is essential to trace the history of the star formation rate (SFR). The main purpose of this thesis is to determine the spatial extent of the dust emission in high-redshift galaxies and to provide a lower limit on dust temperature, to constrain the dust mass. This is achieved by studying 23 FIR continuum detected main-sequence galaxies of the ALMA Large Program to INvestigate (ALPINE) survey, performed at high redshift (4of the dust emission, i.e. the dust size, with the stellar and gas distribution, traced by the UV and [CII] emission, respectively. Finally, we put these results in a broader context, by studying the dust size evolution as a function of cosmic time. We derive dust size measurements via a Gaussian fit in the image and uv plane. Out of the 23 FIR-continuum-detected targets, 20 have been considered in this work since they are isolated systems. Of these 20, 7 are spatially resolved; for each of the remaining 13, we provide an upper limit to the dust size. We find that the gas emission is more extended than the dust spatial scale, by a factor of 1.40±0.29, while the latter appears to be larger than the stellar emission size. Moreover, we do not find any significant trend for dust size as a function of the stellar mass and the redshift. In addition, we provide a minimum dust temperature estimate for the 7 resolved sources, for which we find Tmin∼16−19K. We also derive dust masses for the resolved sources, logMdust∼7−8M⊙.

Stationary states in random walks on networks

In this thesis we dealt with the problem of describing a transportation network in which the objects in movement were subject to both finite transportation capacity and finite accomodation capacity. The movements across such a system are realistically of a simultaneous nature which poses some challenges when formulating a mathematical description. We tried to derive such a general modellization from one posed on a simplified problem based on asyncronicity in particle transitions. We did so considering one-step processes based on the assumption that the system could be describable through discrete time Markov processes with finite state space. After describing the pre-established dynamics in terms of master equations we determined stationary states for the considered processes. Numerical simulations then led to the conclusion that a general system naturally evolves toward a congestion state when its particle transition simultaneously and we consider one single constraint in the form of network node capacity. Moreover the congested nodes of a system tend to be located in adjacent spots in the network, thus forming local clusters of congested nodes.