The Universe 10 Billion years ago: Morphologies, Star-formation rates and Galactic-scale winds in z ∼ 2 Galaxies

This thesis is devoted to the study of the properties of high-redsfhit galaxies in the epoch 1 < z < 3, when a substantial fraction of galaxy mass was assembled, and when the evolution of the star-formation rate density peaked. Following a multi-perspective approach and using the most recent and high-quality data available (spectra, photometry and imaging), the morphologies and the star-formation properties of high-redsfhit galaxies were investigated. Through an accurate morphological analyses, the built up of the Hubble sequence was placed around z ~ 2.5. High-redshift galaxies appear, in general, much more irregular and asymmetric than local ones. Moreover, the occurrence of morphological k-­correction is less pronounced than in the local Universe. Different star-formation rate indicators were also studied. The comparison of ultra-violet and optical based estimates, with the values derived from infra-red luminosity showed that the traditional way of addressing the dust obscuration is problematic, at high-redshifts, and new models of dust geometry and composition are required. Finally, by means of stacking techniques applied to rest-frame ultra-violet spectra of star-forming galaxies at z~2, the warm phase of galactic-scale outflows was studied. Evidence was found of escaping gas at velocities of ~ 100 km/s. Studying the correlation of inter-­stellar absorption lines equivalent widths with galaxy physical properties, the intensity of the outflow-related spectral features was proven to depend strongly on a combination of the velocity dispersion of the gas and its geometry.

Tracing AGN accretion and star formation in Herschel galaxies

In this Thesis, we study the physical properties and the cosmic evolution of AGN and their host galaxies since z∼3. Our analysis exploits samples of star forming galaxies detected with Herschel at far-IR wavelengths (from 70 up to 500 micron) in different extragalactic surveys, such as COSMOS and the deep GOODS (South and North) fields. The broad-band ancillary data available in COSMOS and the GOODS fields, allows us to implement Herschel and Spitzer photometry with multi-wavelength ancillary data. We perform a multicomponent SED-fitting decomposition to decouple the emission due to star formation from that due to AGN accretion, and to estimate both host-galaxy parameters (such as stellar mass, M* and star formation rate, SFR), and nuclear intrinsic bolometric luminosities. We use the individual estimates of AGN bolometric luminosity obtained through SED-fitting decomposition to reconstruct the redshit evolution of the AGN bolometric luminosity function since z∼3. The resulting trends are used to estimate the overall AGN accretion rate density at different cosmic epochs and to trace the first ever estimate of the AGN accretion history from an IR survey. Later on, we focus our study on the connection between AGN accretion and integrated galaxy properties. We analyse the relationships of AGN accretion with galaxy properties in the SFR-M* plane and at different cosmic epochs. Finally, we infer what is the parameter that best correlates with AGN accretion, comparing our results with previous studies and discussing their physical implications in the context of current scenarios of AGN/galaxy evolution.

Dynamics of the halo gas in disc galaxies

This Thesis studies the dynamics of hot and cold gas outside the plane in galaxies like the Milky-Way (extra-planar gas) and focuses on the interaction between disc and halo material. Stationary models for the cold phase of the extra-planar gas are presented. They show that the kinematics of this phase must be influenced by the interaction with an ambient medium that we identify as the hot cosmological corona that surrounds disc galaxies. To study this interaction a novel hydrodynamical code has been implemented and a series of hydrodynamical simulations has been run to investigate the mass and momentum exchange between the cold extra-planar gas clouds and the hot corona. These simulations show that the coronal gas can condense efficiently in the turbulent wakes that form behind the cold clouds and it can be accreted by the disc to sustain star formation. They also predict that the corona cannot be a static structure but it must rotate and lag by approximately 80-120 km/s with respect to the disc. Implications of the results of this Thesis for the evolution of star-forming galaxies and for the large-scale dynamics of galactic coronae are also briefly discussed.

The evolution of massive clumps in star forming regions

In this thesis two related arguments are investigated: - The first stages of the process of massive star formation, investigating the physical conditions and -properties of massive clumps in different evolutionary stages, and their CO depletion; - The influence that high-mass stars have on the nearby material and on the activity of star formation. I characterise the gas and dust temperature, mass and density of a sample of massive clumps, and analyse the variation of these properties from quiescent clumps, without any sign of active star formation, to clumps likely hosting a zero-age main sequence star. I briefly discuss CO depletion and recent observations of several molecular species, tracers of Hot Cores and/or shocked gas, of a subsample of these clumps. The issue of CO depletion is addressed in more detail in a larger sample consisting of the brightest sources in the ATLASGAL survey: using a radiative tranfer code I investigate how the depletion changes from dark clouds to more evolved objects, and compare its evolution to what happens in the low-mass regime. Finally, I derive the physical properties of the molecular gas in the photon-dominated region adjacent to the HII region G353.2+0.9 in the vicinity of Pismis 24, a young, massive cluster, containing some of the most massive and hottest stars known in our Galaxy. I derive the IMF of the cluster and study the star formation activity in its surroundings. Much of the data analysis is done with a Bayesian approach. Therefore, a separate chapter is dedicated to the concepts of Bayesian statistics.

X-ray properties and evolution of high-redshift AGN, and the gas content of host galaxies

In this thesis, I have investigated the evolution of the high-redshift (z > 3) AGN population by collecting data from some of the major Chandra and XMM-Newton surveys. The final sample (141 sources) is one of the largest selected at z> 3 in the X- rays and it is characterised by a very high redshift completeness (98%). I derived the spectral slopes and obscurations through a spectral anaysis and I assessed the high-z evolution by deriving the luminosity function and the number counts of the sample. The best representation of the AGN evolution is a pure density evolution (PDE) model: the AGN space density is found to decrease by a factor of 10 from z=3 to z=5. I also found that about 50% of AGN are obscured by large column densities (logNH > 23). By comparing these data with those in the Local Universe, I found a positive evolution of the obscured AGN fraction with redshift, especially for luminous (logLx > 44) AGN. I also studied the gas content of z < 1 AGN-hosting galaxies and compared it with that of inactive galaxies. For the first time, I applied to AGN a method to derive the gas mass previously used for inactive galaxies only. AGN are found to live preferentially in gas-rich galaxies. This result on the one hand can help us in understanding the AGN triggering mechanisms, on the other hand explains why AGN are preferentially hosted by star-forming galaxies.

Millimeter characteristics of Galactic and extragalactic starbursts

The goal of the thesis was to compare Galactic and extragalactic star-forming environments, to understand whether we can derive scaling relations for the extragalactic star formation from our understanding of the Galactic star formation; and to analyze the effect of the angular resolution of the observations and the molecular tracer used in extragalactic studies. It is therefore necessary to perform observations at the same linear resolution in both cases. With this in mind we have started a program aimed at comparing properties of mini-starburst regions in our Galaxy and starbursts in nearby galaxies at similar linear resolutions.

Improving the cosmic distance ladder. Distance and structure of the Large Magellanic Cloud

The Large Magellanic Cloud (LMC) is widely considered as the first step of the cosmological distance ladder, since it contains many different distance indicators. An accurate determination of the distance to the LMC allows one to calibrate these distance indicators that are then used to measure the distance to far objects. The main goal of this thesis is to study the distance and structure of the LMC, as traced by different distance indicators. For these purposes three types of distance indicators were chosen: Classical Cepheids,``hot'' eclipsing binaries and RR Lyrae stars. These objects belong to different stellar populations tracing, in turn, different sub-structures of the LMC. The RR Lyrae stars (age >10 Gyr) are distributed smoothly and likely trace the halo of the LMC. Classical Cepheids are young objects (age 50-200 Myr), mainly located in the bar and spiral arm of the galaxy, while ``hot'' eclipsing binaries mainly trace the star forming regions of the LMC. Furthermore, we have chosen these distance indicators for our study, since the calibration of their zero-points is based on fundamental geometric methods. The ESA cornerstone mission Gaia, launched on 19 December 2013, will measure trigonometric parallaxes for one billion stars with an accuracy of 20 micro-arcsec at V=15 mag, and 200 micro-arcsec at V=20 mag, thus will allow us to calibrate the zero-points of Classical Cepheids, eclipsing binaries and RR Lyrae stars with an unprecedented precision.

Modeling the cosmological co-evolution of supermassive black holes and galaxies

In this Thesis, we investigate the cosmological co-evolution of supermassive black holes (BHs), Active Galactic Nuclei (AGN) and their hosting dark matter (DM) halos and galaxies, within the standard CDM scenario. We analyze both analytic, semi-analytic and hybrid techniques and use the most recent observational data available to constrain the assumptions underlying our models. First, we focus on very simple analytic models where the assembly of BHs is directly related to the merger history of DM haloes. For this purpose, we implement the two original analytic models of Wyithe & Loeb 2002 and Wyithe & Loeb 2003, compare their predictions to the AGN luminosity function and clustering data, and discuss possible modifications to the models that improve the match to the observation. Then we study more sophisticated semi-analytic models in which however the baryonic physics is neglected as well. Finally we improve the hybrid simulation of De Lucia & Blaizot 2007, adding new semi-analytical prescriptions to describe the BH mass accretion rate during each merger event and its conversion into radiation, and compare the derived BH scaling relations, fundamental plane and mass function, and the AGN luminosity function with observations. All our results support the following scenario: • The cosmological co-evolution of BHs, AGN and galaxies can be well described within the CDM model. • At redshifts z & 1, the evolution history of DM halo fully determines the overall properties of the BH and AGN populations. The AGN emission is triggered mainly by DM halo major mergers and, on average, AGN shine at their Eddington luminosity. • At redshifts z . 1, BH growth decouples from halo growth. Galaxy major mergers cannot constitute the only trigger to accretion episodes in this phase. • When a static hot halo has formed around a galaxy, a fraction of the hot gas continuously accretes onto the central BH, causing a low-energy “radio” activity at the galactic centre, which prevents significant gas cooling and thus limiting the mass of the central galaxies and quenching the star formation at late time. • The cold gas fraction accreted by BHs at high redshifts seems to be larger than at low redshifts.

The gaseous environment of radio galaxies: a new perspective from high-resolution x-ray spectroscopy

It is known that massive black holes have a profound effect on the evolution of galaxies, and possibly on their formation by regulating the amount of gas available for the star formation. However, how black hole and galaxies communicate is still an open problem, depending on how much of the energy released interacts with the circumnuclear matter. In the last years, most studies of feedback have primarily focused on AGN jet/cavity systems in the most massive galaxy clusters. This thesis intends to investigate the feedback phenomena in radio--loud AGNs from a different perspective studying isolated radio galaxies, through high-resolution spectroscopy. In particular one NLRG and three BLRG are studied, searching for warm gas, both in emission and absorption, in the soft X-ray band. I show that the soft spectrum of 3C33 originates from gas photoionized by the central engine. I found for the first time WA in 3C382 and 3C390.3. I show that the observed warm emitter/absorbers is not uniform and probably located in the NLR. The detected WA is slow implying a mass outflow rate and kinetic luminosity always well below 1% the L(acc) as well as the P(jet). Finally the radio--loud properties are compared with those of type 1 RQ AGNs. A positive correlation is found between the mass outflow rate/kinetic luminosity, and the radio loudness. This seems to suggest that the presence of a radio source (the jet?) affects the distribution of the absorbing gas. Alternatively, if the gas distribution is similar in Seyferts and radio galaxies, the M(out) vs rl relation could simply indicate a major ejection of matter in the form of wind in powerful radio AGNs.

Hydrodynamical simulations of early-type galaxies: effects of galaxy shape and stellar dynamics on hot coronae

Early-Type galaxies (ETGs) are embedded in hot (10^6-10^7 K), X-ray emitting gaseous haloes, produced mainly by stellar winds and heated by Type Ia supernovae explosions, by the thermalization of stellar motions and occasionally by the central super-massive black hole (SMBH). In particular, the thermalization of the stellar motions is due to the interaction between the stellar and the SNIa ejecta and the hot interstellar medium (ISM) already residing in the ETG. A number of different astrophysical phenomena determine the X-ray properties of the hot ISM, such as stellar population formation and evolution, galaxy structure and internal kinematics, Active Galactic Nuclei (AGN) presence, and environmental effects. With the aid of high-resolution hydrodynamical simulations performed on state-of-the-art galaxy models, in this Thesis we focus on the effects of galaxy shape, stellar kinematics and star formation on the evolution of the X-ray coronae of ETGs. Numerical simulations show that the relative importance of flattening and rotation are functions of the galaxy mass: at low galaxy masses, adding flattening and rotation induces a galactic wind, thus lowering the X-ray luminosity; at high galaxy masses the angular momentum conservation keeps the central regions of rotating galaxies at low density, whereas in non-rotating models a denser and brighter atmosphere is formed. The same dependence from the galaxy mass is present in the effects of star formation (SF): in light galaxies SF contributes to increase the spread in Lx, while at high galaxy masses the halo X-ray properties are marginally sensitive to SF effects. In every case, the star formation rate at the present epoch quite agrees with observations, and the massive, cold gaseous discs are partially or completely consumed by SF on a time-scale of few Gyr, excluding the presence of young stellar discs at the present epoch.

Accretion of Mass and Angular Momentum onto The Discs of Spiral Galaxies

In this Thesis, we study the accretion of mass and angular momentum onto the disc of spiral galaxies from a global and a local perspective and comparing theory predictions with several observational data. First, we propose a method to measure the specific mass and radial growth rates of stellar discs, based on their star formation rate density profiles and we apply it to a sample of nearby spiral galaxies. We find a positive radial growth rate for almost all galaxies in our sample. Our galaxies grow in size, on average, at one third of the rate at which they grow in mass. Our results are in agreement with theoretical expectations if known scaling relations of disc galaxies are not evolving with time. We also propose a novel method to reconstruct accretion profiles and the local angular momentum of the accreting material from the observed structural and chemical properties of spiral galaxies. Applied to the Milky Way and to one external galaxy, our analysis indicates that accretion occurs at relatively large radii and has a local deficit of angular momentum with respect to the disc. Finally, we show how structure and kinematics of hot gaseous coronae, which are believed to be the source of mass and angular momentum of massive spiral galaxies, can be reconstructed from their angular momentum and entropy distributions. We find that isothermal models with cosmologically motivated angular momentum distributions are compatible with several independent observational constraints. We also consider more complex baroclinic equilibria: we describe a new parametrization for these states, a new self-similar family of solution and a method for reconstructing structure and kinematics from the joint angular momentum/entropy distribution.

On the role of AGN feedback in shaping galaxies across cosmic time

Understanding how Active Galactic Nuclei (AGN) shape galaxy evolution is a key challenge of modern astronomy. In the framework where black hole (BH) and galaxy growth are linked, AGN feedback must be tackled both at its “causes” (e.g. AGN-driven winds) and its “effects” (alteration of the gas reservoir in AGN hosts). The most informative cosmic time is z~1-3, at the peak of AGN activity and galaxy buildup, the so-called cosmic noon. The aim of this thesis is to provide new insights regarding some key questions that still remain open in this research field: i) What are the properties of AGN-driven sub-pc scale winds at z>1? ii) Are AGN-driven winds effective in influencing the life of galaxies? iii) Do AGN impact directly on star formation (SF) and gas content of their hosts? I first address AGN feedback as “caught in the act” by studying ultra-fast outflows (UFOs), X-ray AGN-driven winds, in gravitationally lensed quasars. I build the first statistically robust sample of high-z AGN, not preselected based on AGN-driven winds. I derive a first estimate of the high-z UFO detection fraction and measure the UFO duty cycle of a single high-z quasar for the first time. I also address the “effects” of AGN feedback on the life of host galaxies. If AGN influence galaxy growth, then they will reasonably impact the molecular gas reservoir first, and SF as a consequence. Through a comparative study of the molecular gas content in cosmic-noon AGN hosts and matched non-active galaxies (i.e., galaxies not hosting an AGN), we find that the host galaxies of more regular AGN (not selected to be the most luminous) are generally similar to non-active galaxies. However, we report on the possibility of a luminosity effect regulating the efficiency by which AGN might impact on galaxy growth.

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.

A high energy perspective of the co-evolution of black holes and their host galaxies

Thanks to the Chandra and XMM–Newton surveys, the hard X-ray sky is now probed down to a flux limit where the bulk of the X-ray background is almost completely resolved into discrete sources, at least in the 2–8 keV band. Extensive programs of multiwavelength follow-up observations showed that the large majority of hard X–ray selected sources are identified with Active Galactic Nuclei (AGN) spanning a broad range of redshifts, luminosities and optical properties. A sizable fraction of relatively luminous X-ray sources hosting an active, presumably obscured, nucleus would not have been easily recognized as such on the basis of optical observations because characterized by “peculiar” optical properties. In my PhD thesis, I will focus the attention on the nature of two classes of hard X-ray selected “elusive” sources: those characterized by high X-ray-to-optical flux ratios and red optical-to-near-infrared colors, a fraction of which associated with Type 2 quasars, and the X-ray bright optically normal galaxies, also known as XBONGs. In order to characterize the properties of these classes of elusive AGN, the datasets of several deep and large-area surveys have been fully exploited. The first class of “elusive” sources is characterized by X-ray-to-optical flux ratios (X/O) significantly higher than what is generally observed from unobscured quasars and Seyfert galaxies. The properties of well defined samples of high X/O sources detected at bright X–ray fluxes suggest that X/O selection is highly efficient in sampling high–redshift obscured quasars. At the limits of deep Chandra surveys (∼10−16 erg cm−2 s−1), high X/O sources are generally characterized by extremely faint optical magnitudes, hence their spectroscopic identification is hardly feasible even with the largest telescopes. In this framework, a detailed investigation of their X-ray properties may provide useful information on the nature of this important component of the X-ray source population. The X-ray data of the deepest X-ray observations ever performed, the Chandra deep fields, allows us to characterize the average X-ray properties of the high X/O population. The results of spectral analysis clearly indicate that the high X/O sources represent the most obscured component of the X–ray background. Their spectra are harder (G ∼ 1) than any other class of sources in the deep fields and also of the XRB spectrum (G ≈ 1.4). In order to better understand the AGN physics and evolution, a much better knowledge of the redshift, luminosity and spectral energy distributions (SEDs) of elusive AGN is of paramount importance. The recent COSMOS survey provides the necessary multiwavelength database to characterize the SEDs of a statistically robust sample of obscured sources. The combination of high X/O and red-colors offers a powerful tool to select obscured luminous objects at high redshift. A large sample of X-ray emitting extremely red objects (R−K >5) has been collected and their optical-infrared properties have been studied. In particular, using an appropriate SED fitting procedure, the nuclear and the host galaxy components have been deconvolved over a large range of wavelengths and ptical nuclear extinctions, black hole masses and Eddington ratios have been estimated. It is important to remark that the combination of hard X-ray selection and extreme red colors is highly efficient in picking up highly obscured, luminous sources at high redshift. Although the XBONGs do not present a new source population, the interest on the nature of these sources has gained a renewed attention after the discovery of several examples from recent Chandra and XMM–Newton surveys. Even though several possibilities were proposed in recent literature to explain why a relatively luminous (LX = 1042 − 1043erg s−1) hard X-ray source does not leave any significant signature of its presence in terms of optical emission lines, the very nature of XBONGs is still subject of debate. Good-quality photometric near-infrared data (ISAAC/VLT) of 4 low-redshift XBONGs from the HELLAS2XMMsurvey have been used to search for the presence of the putative nucleus, applying the surface-brightness decomposition technique. In two out of the four sources, the presence of a nuclear weak component hosted by a bright galaxy has been revealed. The results indicate that moderate amounts of gas and dust, covering a large solid angle (possibly 4p) at the nuclear source, may explain the lack of optical emission lines. A weak nucleus not able to produce suffcient UV photons may provide an alternative or additional explanation. On the basis of an admittedly small sample, we conclude that XBONGs constitute a mixed bag rather than a new source population. When the presence of a nucleus is revealed, it turns out to be mildly absorbed and hosted by a bright galaxy.

A spectroscopic and photometric study of stellar populations in a sample of clusters in the Magellanic Clouds

This Ph.D. Thesis has been carried out in the framework of a long-term and large project devoted to describe the main photometric, chemical, evolutionary and integrated properties of a representative sample of Large and Small Magellanic Cloud (LMC and SMC respectively) clusters. The globular clusters system of these two Irregular galaxies provides a rich resource for investigating stellar and chemical evolution and to obtain a detailed view of the star formation history and chemical enrichment of the Clouds. The results discussed here are based on the analysis of high-resolution photometric and spectroscopic datasets obtained by using the last generation of imagers and spectrographs. The principal aims of this project are summarized as follows: • The study of the AGB and RGB sequences in a sample of MC clusters, through the analysis of a wide near-infrared photometric database, including 33 Magellanic globulars obtained in three observing runs with the near-infrared camera SOFI@NTT (ESO, La Silla). • The study of the chemical properties of a sample of MCs clusters, by using optical and near-infrared high-resolution spectra. 3 observing runs have been secured to our group to observe 9 LMC clusters (with ages between 100 Myr and 13 Gyr) with the optical high-resolution spectrograph FLAMES@VLT (ESO, Paranal) and 4 very young (<30 Myr) clusters (3 in the LMC and 1 in the SMC) with the near-infrared high-resolution spectrograph CRIRES@VLT. • The study of the photometric properties of the main evolutive sequences in optical Color- Magnitude Diagrams (CMD) obtained by using HST archive data, with the final aim of dating several clusters via the comparison between the observed CMDs and theoretical isochrones. The determination of the age of a stellar population requires an accurate measure of the Main Sequence (MS) Turn-Off (TO) luminosity and the knowledge of the distance modulus, reddening and overall metallicity. For this purpose, we limited the study of the age just to the clusters already observed with high-resolution spectroscopy, in order to date only clusters with accurate estimates of the overall metallicity.