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.

The molecular gas content of star-forming galaxies in the ALPINE fields.

In this thesis, I aim to study the evolution with redshift of the gas mass fraction of a sample of 53 sources (from z ∼ 0.5 to z > 5) serendipitously detected in ALMA band 7 as part of the ALMA Large Program to INvestigate C II at Early Times (ALPINE). First, I used SED-fitting software CIGALE, which is able to implement energy balancing between the optical and the far infrared part, to produce a best-fit template of my sources and to have an estimate of some physical properties, such as the star formation rate (SFR), the total infrared luminosity and the total stellar mass. Then, using the tight correlation found by Scoville et al. (2014) between the ISM molecular gas mass and the rest-frame 850 μm luminosity, I used the latter, extrapolating it from the best-fit template using a code that I wrote in Python, as a tracer for the molecular gas. For my sample, I then derived the most important physical properties, such as molecular gas mass, gas mass fractions, specific star formation rate and depletion timescales, which allowed me to better categorize them and find them a place within the evolutionary history of the Universe. I also fitted our sources, via another code I wrote again in Python, with a general modified blackbody (MBB) model taken from the literature (Gilli et al. (2014), D’Amato et al. (2020)) to have a direct method of comparison with similar galaxies. What is evident at the end of the paper is that the methods used to derive the physical quantities of the sources are consistent with each other, and these in turn are in good agreement with what is found in the literature.

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 (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 (4

The structure and dynamics of Abell 1942

Aims. We present a dynamical analysis of the galaxy cluster Abell 1942 based on a set of 128 velocities obtained at the European Southern Observatory. Methods. Data on individual galaxies are presented and the accuracy of the determined velocities as some properties of the cluster are discussed. We have also made use of publicly available Chandra X-ray data. Results. We obtained an improved mean redshift value z = 0.22513 +/- 0.0008 and velocity dispersion sigma = 908(139)(+147) km s(-1). Our analysis indicates that inside a radius of similar to 1.5 h(70)(-1) Mpc (similar to 7 arcmin) the cluster is well relaxed, without any remarkable features and the X-ray emission traces the galaxy distribution fairly well. Two possible optical substructures are seen at similar to 5 arcmin from the centre in the northwest and the southwest directions, but are not confirmed by the velocity field. These clumps are, however, kinematically bound to the main structure of Abell 1942. X-ray spectroscopic analysis of Chandra data resulted in a temperature kT = 5.5+/-0.5 keV and metal abundance Z = 0.33 +/- 0.15 Z(circle dot). The velocity dispersion corresponding to this temperature using the T(X-sigma) scaling relation is in good agreement with the measured galaxy velocities. Our photometric redshift analysis suggests that the weak lensing signal observed to the south of the cluster and previously attributed to a ""dark clump"" is produced by background sources, possibly distributed as a filamentary structure.

NGC 7097: THE ACTIVE GALACTIC NUCLEUS AND ITS MIRROR, REVEALED BY PRINCIPAL COMPONENT ANALYSIS TOMOGRAPHY

Three-dimensional spectroscopy techniques are becoming more and more popular, producing an increasing number of large data cubes. The challenge of extracting information from these cubes requires the development of new techniques for data processing and analysis. We apply the recently developed technique of principal component analysis (PCA) tomography to a data cube from the center of the elliptical galaxy NGC 7097 and show that this technique is effective in decomposing the data into physically interpretable information. We find that the first five principal components of our data are associated with distinct physical characteristics. In particular, we detect a low-ionization nuclear-emitting region (LINER) with a weak broad component in the Balmer lines. Two images of the LINER are present in our data, one seen through a disk of gas and dust, and the other after scattering by free electrons and/or dust particles in the ionization cone. Furthermore, we extract the spectrum of the LINER, decontaminated from stellar and extended nebular emission, using only the technique of PCA tomography. We anticipate that the scattered image has polarized light due to its scattered nature.

Abundance variations in the globular cluster M71 (NGC 6838)

Context. Abundance variations in moderately metal-rich globular clusters can give clues about the formation and chemical enrichment of globular clusters. Aims. CN, CH, Na, Mg and Al indices in spectra of 89 stars of the template metal-rich globular cluster M71 are measured and implications on internal mixing are discussed. Methods. Stars from the turn-off up to the Red Giant Branch (0.87 < log g < 4.65) observed with the GMOS multi-object spectrograph at the Gemini-North telescope are analyzed. Radial velocities, colours, effective temperatures, gravities and spectral indices are determined for the sample. Results. Previous findings related to the CN bimodality and CN-CH anticorrelation in stars of M71 are confirmed. We also find a CN-Na correlation, and Al-Na, as well as an Mg(2)-Al anticorrelation. Conclusions. A combination of convective mixing and a primordial pollution by AGB or massive stars in the early stages of globular cluster formation is required to explain the observations.

Star formation efficiency in galaxy clusters

Context. The luminous material in clusters of galaxies exists in two forms: the visible galaxies and the X-ray emitting intra-cluster medium. The hot intra-cluster gas is the major observed baryonic component of clusters, about six times more massive than the stellar component. The mass contained within visible galaxies is approximately 3% of the dynamical mass. Aims. Our aim was to analyze both baryonic components, combining X-ray and optical data of a sample of five galaxy clusters (Abell 496, 1689, 2050, 2631 and 2667), within the redshift range 0.03 < z < 0.3. We determined the contribution of stars in galaxies and the intra-cluster medium to the total baryon budget. Methods. We used public XMM-Newton data to determine the gas mass and to obtain the X-ray substructures. Using the optical counterparts from SDSS or CFHT we determined the stellar contribution. Results. We examine the relative contribution of galaxies, intra-cluster light and intra-cluster medium to baryon budget in clusters through the stellar-to-gas mass ratio, estimated with recent data. We find that the stellar-to-gas mass ratio within r(500) (the radius within which the mean cluster density exceeds the critical density by a factor of 500), is anti-correlated with the ICM temperature, which range from 24% to 6% while the temperature ranges from 4.0 to 8.3 keV. This indicates that less massive cold clusters are more prolific star forming environments than massive hot clusters.

Rotation curve bifurcations as indicators of close recent galaxy encounters

Context. Rotation curves of interacting galaxies often show that velocities are either rising or falling in the direction of the companion galaxy. Aims. We seek to reproduce and analyse these features in the rotation curves of simulated equal-mass galaxies suffering a one-to-one encounter as possible indicators of close encounters. Methods. Using simulations of major mergers in 3D, we study the time evolution of these asymmetries in a pair of galaxies during the first passage. Results. Our main results are: (a) the rotation curve asymmetries appear right at pericentre of the first passage, (b) the significant disturbed rotation velocities occur within a small time interval, of similar to 0.5 Gyr h(-1), and, therefore, the presence of bifurcation in the velocity curve could be used as an indicator of the pericentre occurrence. These results are in qualitative agreement with previous findings for minor mergers and flybys.

Origin of primordial magnetic fields

Magnetic fields of intensities similar to those in our galaxy are also observed in high redshift galaxies, where a mean field dynamo would not have had time to produce them. Therefore, a primordial origin is indicated. It has been suggested that magnetic fields were created at various primordial eras: during inflation, the electroweak phase transition, the quark-hadron phase transition (QHPT), during the formation of the first objects, and during reionization. We suggest here that the large-scale fields similar to mu G, observed in galaxies at both high and low redshifts by Faraday rotation measurements (FRMs), have their origin in the electromagnetic fluctuations that naturally occurred in the dense hot plasma that existed just after the QHPT. We evolve the predicted fields to the present time. The size of the region containing a coherent magnetic field increased due to the fusion of smaller regions. Magnetic fields (MFs) similar to 10 mu G over a comoving similar to 1 pc region are predicted at redshift z similar to 10. These fields are orders of magnitude greater than those predicted in previous scenarios for creating primordial magnetic fields. Line-of-sight average MFs similar to 10(-2) mu G, valid for FRMs, are obtained over a 1 Mpc comoving region at the redshift z similar to 10. In the collapse to a galaxy (comoving size similar to 30 kpc) at z similar to 10, the fields are amplified to similar to 10 mu G. This indicates that the MFs created immediately after the QHPT (10(-4) s), predicted by the fluctuation-dissipation theorem, could be the origin of the similar to mu G fields observed by FRMs in galaxies at both high and low redshifts. Our predicted MFs are shown to be consistent with present observations. We discuss the possibility that the predicted MFs could cause non-negligible deflections of ultrahigh energy cosmic rays and help create the observed isotropic distribution of their incoming directions. We also discuss the importance of the volume average magnetic field predicted by our model in producing the first stars and in reionizing the Universe.

Alpha element abundances and gradients in the Milky Way bulge from FLAMES-GIRAFFE spectra of 650 K giants

Aims. We present the analysis of the [alpha/Fe] abundance ratios for a large number of stars at several locations in the Milky Way bulge with the aim of constraining its formation scenario. Methods. We obtained FLAMES-GIRAFFE spectra (R = 22 500) at the ESO Very Large Telescope for 650 bulge red giant branch (RGB) stars and performed spectral synthesis to measure Mg, Ca, Ti, and Si abundances. This sample is composed of 474 giant stars observed in 3 fields along the minor axis of the Galactic bulge and at latitudes b = -4 degrees, b = -6 degrees, b = -12 degrees. Another 176 stars belong to a field containing the globular cluster NGC 6553, located at b = -3 degrees and 5 degrees away from the other three fields along the major axis. Stellar parameters and metallicities for these stars were presented in Zoccali et al. (2008, A&A, 486, 177). We have also re-derived stellar parameters and abundances for the sample of thick and thin disk red giants analyzed in Alves-Brito et al. (2010, A&A, 513, A35). Therefore using a homogeneous abundance database for the bulge, thick and thin disk, we have performed a differential analysis minimizing systematic errors, to compare the formation scenarios of these Galactic components. Results. Our results confirm, with large number statistics, the chemical similarity between the Galactic bulge and thick disk, which are both enhanced in alpha elements when compared to the thin disk. In the same context, we analyze [alpha/Fe] vs. [Fe/H] trends across different bulge regions. The most metal rich stars, showing low [alpha/Fe] ratios at b = -4 degrees disappear at higher Galactic latitudes in agreement with the observed metallicity gradient in the bulge. Metal-poor stars ([Fe/H] < -0.2) show a remarkable homogeneity at different bulge locations. Conclusions. We have obtained further constrains for the formation scenario of the Galactic bulge. A metal-poor component chemically indistinguishable from the thick disk hints for a fast and early formation for both the bulge and the thick disk. Such a component shows no variation, neither in abundances nor kinematics, among different bulge regions. A metal-rich component showing low [alpha/Fe] similar to those of the thin disk disappears at larger latitudes. This allows us to trace a component formed through fast early mergers (classical bulge) and a disk/bar component formed on a more extended timescale.

Dynamics and constraints of the massive graviton dark matter flat cosmologies

We discuss the dynamics of the Universe within the framework of the massive graviton cold dark matter scenario (MGCDM) in which gravitons are geometrically treated as massive particles. In this modified gravity theory, the main effect of the gravitons is to alter the density evolution of the cold dark matter component in such a way that the Universe evolves to an accelerating expanding regime, as presently observed. Tight constraints on the main cosmological parameters of the MGCDM model are derived by performing a joint likelihood analysis involving the recent supernovae type Ia data, the cosmic microwave background shift parameter, and the baryonic acoustic oscillations as traced by the Sloan Digital Sky Survey red luminous galaxies. The linear evolution of small density fluctuations is also analyzed in detail. It is found that the growth factor of the MGCDM model is slightly different (similar to 1-4%) from the one provided by the conventional flat Lambda CDM cosmology. The growth rate of clustering predicted by MGCDM and Lambda CDM models are confronted to the observations and the corresponding best fit values of the growth index (gamma) are also determined. By using the expectations of realistic future x-ray and Sunyaev-Zeldovich cluster surveys we derive the dark matter halo mass function and the corresponding redshift distribution of cluster-size halos for the MGCDM model. Finally, we also show that the Hubble flow differences between the MGCDM and the Lambda CDM models provide a halo redshift distribution departing significantly from the those predicted by other dark energy models. These results suggest that the MGCDM model can observationally be distinguished from Lambda CDM and also from a large number of dark energy models recently proposed in the literature.

Cosmic distance duality relation and the shape of galaxy clusters

Context. Observations in the cosmological domain are heavily dependent on the validity of the cosmic distance-duality (DD) relation, eta = D(L)(z)(1+ z)(2)/D(A)(z) = 1, an exact result required by the Etherington reciprocity theorem where D(L)(z) and D(A)(z) are, respectively, the luminosity and angular diameter distances. In the limit of very small redshifts D(A)(z) = D(L)(z) and this ratio is trivially satisfied. Measurements of Sunyaev-Zeldovich effect (SZE) and X-rays combined with the DD relation have been used to determine D(A)(z) from galaxy clusters. This combination offers the possibility of testing the validity of the DD relation, as well as determining which physical processes occur in galaxy clusters via their shapes. Aims. We use WMAP (7 years) results by fixing the conventional Lambda CDM model to verify the consistence between the validity of DD relation and different assumptions about galaxy cluster geometries usually adopted in the literature. Methods. We assume that. is a function of the redshift parametrized by two different relations: eta(z) = 1+eta(0)z, and eta(z) = 1+eta(0)z/(1+z), where eta(0) is a constant parameter quantifying the possible departure from the strict validity of the DD relation. In order to determine the probability density function (PDF) of eta(0), we consider the angular diameter distances from galaxy clusters recently studied by two different groups by assuming elliptical (isothermal) and spherical (non-isothermal) beta models. The strict validity of the DD relation will occur only if the maximum value of eta(0) PDF is centered on eta(0) = 0. Results. It was found that the elliptical beta model is in good agreement with the data, showing no violation of the DD relation (PDF peaked close to eta(0) = 0 at 1 sigma), while the spherical (non-isothermal) one is only marginally compatible at 3 sigma. Conclusions. The present results derived by combining the SZE and X-ray surface brightness data from galaxy clusters with the latest WMAP results (7-years) favors the elliptical geometry for galaxy clusters. It is remarkable that a local property like the geometry of galaxy clusters might be constrained by a global argument provided by the cosmic DD relation.

Fossil groups in the Millennium simulation Their environment and its evolution

Context. Fossil systems are defined to be X- ray bright galaxy groups ( or clusters) with a two- magnitude difference between their two brightest galaxies within half the projected virial radius, and represent an interesting extreme of the population of galaxy agglomerations. However, the physical conditions and processes leading to their formation are still poorly constrained. Aims. We compare the outskirts of fossil systems with that of normal groups to understand whether environmental conditions play a significant role in their formation. We study the groups of galaxies in both, numerical simulations and observations. Methods. We use a variety of statistical tools including the spatial cross- correlation function and the local density parameter Delta(5) to probe differences in the density and structure of the environments of "" normal"" and "" fossil"" systems in the Millennium simulation. Results. We find that the number density of galaxies surrounding fossil systems evolves from greater than that observed around normal systems at z = 0.69, to lower than the normal systems by z = 0. Both fossil and normal systems exhibit an increment in their otherwise radially declining local density measure (Delta(5)) at distances of order 2.5 r(vir) from the system centre. We show that this increment is more noticeable for fossil systems than normal systems and demonstrate that this difference is linked to the earlier formation epoch of fossil groups. Despite the importance of the assembly time, we show that the environment is different for fossil and non- fossil systems with similar masses and formation times along their evolution. We also confirm that the physical characteristics identified in the Millennium simulation can also be detected in SDSS observations. Conclusions. Our results confirm the commonly held belief that fossil systems assembled earlier than normal systems but also show that the surroundings of fossil groups could be responsible for the formation of their large magnitude gap.