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.

Development of a new galaxy classification technique and its application to the zCOSMOS survey

The purpose of this Thesis is to develop a robust and powerful method to classify galaxies from large surveys, in order to establish and confirm the connections between the principal observational parameters of the galaxies (spectral features, colours, morphological indices), and help unveil the evolution of these parameters from $z \sim 1$ to the local Universe. Within the framework of zCOSMOS-bright survey, and making use of its large database of objects ($\sim 10\,000$ galaxies in the redshift range $0 < z \lesssim 1.2$) and its great reliability in redshift and spectral properties determinations, first we adopt and extend the \emph{classification cube method}, as developed by Mignoli et al. (2009), to exploit the bimodal properties of galaxies (spectral, photometric and morphologic) separately, and then combining together these three subclassifications. We use this classification method as a test for a newly devised statistical classification, based on Principal Component Analysis and Unsupervised Fuzzy Partition clustering method (PCA+UFP), which is able to define the galaxy population exploiting their natural global bimodality, considering simultaneously up to 8 different properties. The PCA+UFP analysis is a very powerful and robust tool to probe the nature and the evolution of galaxies in a survey. It allows to define with less uncertainties the classification of galaxies, adding the flexibility to be adapted to different parameters: being a fuzzy classification it avoids the problems due to a hard classification, such as the classification cube presented in the first part of the article. The PCA+UFP method can be easily applied to different datasets: it does not rely on the nature of the data and for this reason it can be successfully employed with others observables (magnitudes, colours) or derived properties (masses, luminosities, SFRs, etc.). The agreement between the two classification cluster definitions is very high. ``Early'' and ``late'' type galaxies are well defined by the spectral, photometric and morphological properties, both considering them in a separate way and then combining the classifications (classification cube) and treating them as a whole (PCA+UFP cluster analysis). Differences arise in the definition of outliers: the classification cube is much more sensitive to single measurement errors or misclassifications in one property than the PCA+UFP cluster analysis, in which errors are ``averaged out'' during the process. This method allowed us to behold the \emph{downsizing} effect taking place in the PC spaces: the migration between the blue cloud towards the red clump happens at higher redshifts for galaxies of larger mass. The determination of $M_{\mathrm{cross}}$ the transition mass is in significant agreement with others values in literature.

Stellar Relics of the hierarchical assembly of the Galaxy

This PhD Thesis is part of a long-term wide research project, carried out by the "Osservatorio Astronomico di Bologna (INAF-OABO)", that has as primary goal the comprehension and reconstruction of formation mechanism of galaxies and their evolution history. There is now substantial evidence, both from theoretical and observational point of view, in favor of the hypothesis that the halo of our Galaxy has been at least partially, built up by the progressive accretion of small fragments, similar in nature to the present day dwarf galaxies of the Local Group. In this context, the photometric and spectroscopic study of systems which populate the halo of our Galaxy (i.e. dwarf spheroidal galaxy, tidal streams, massive globular cluster, etc) permits to discover, not only the origin and behaviour of these systems, but also the structure of our Galactic halo, combined with its formation history. In fact, the study of the population of these objects and also of their chemical compositions, age, metallicities and velocity dispersion, permit us not only an improvement in the understanding of the mechanisms that govern the Galactic formation, but also a valid indirect test for cosmological model itself. Specifically, in this Thesis we provided a complete characterization of the tidal Stream of the Sagittarius dwarf spheroidal galaxy, that is the most striking example of the process of tidal disruption and accretion of a dwarf satellite in to our Galaxy. Using Red Clump stars, extracted from the catalogue of the Sloan Digital Sky Survey (SDSS) we obtained an estimate of the distance, the depth along the line of sight and of the number density for each detected portion of the Stream (and more in general for each detected structure along our line of sight). Moreover comparing the relative number (i.e. the ratio) of Blue Horizontal Branch stars and Red Clump stars (the two features are tracers of different age/different metallicity populations) in the main body of the galaxy and in the Stream, in order to verify the presence of an age-metallicity gradient along the Stream. We also report the detection of a population of Red Clump stars probably associated with the recently discovered Bootes III stellar system. Finally, we also present the results of a survey of radial velocities over a wide region, extending from r ~ 10' out to r ~ 80' within the massive star cluster Omega Centauri. The survey was performed with FLAMES@VLT, to study the velocity dispersion profile in the outer regions of this stellar system. All the results presented in this Thesis, have already been published in refeered journals.

Hydrodynamical simulations of galaxy clusters in dark energy cosmologies

In this work we investigate the influence of dark energy on structure formation, within five different cosmological models, namely a concordance $\Lambda$CDM model, two models with dynamical dark energy, viewed as a quintessence scalar field (using a RP and a SUGRA potential form) and two extended quintessence models (EQp and EQn) where the quintessence scalar field interacts non-minimally with gravity (scalar-tensor theories). We adopted for all models the normalization of the matter power spectrum $\sigma_{8}$ to match the CMB data. For each model, we perform hydrodynamical simulations in a cosmological box of $(300 \ {\rm{Mpc}} \ h^{-1})^{3}$ including baryons and allowing for cooling and star formation. We find that, in models with dynamical dark energy, the evolving cosmological background leads to different star formation rates and different formation histories of galaxy clusters, but the baryon physics is not affected in a relevant way. We investigate several proxies for the cluster mass function based on X-ray observables like temperature, luminosity, $M_{gas}$, and $Y_{X}$. We confirm that the overall baryon fraction is almost independent of the dark energy models within few percentage points. The same is true for the gas fraction. This evidence reinforces the use of galaxy clusters as cosmological probe of the matter and energy content of the Universe. We also study the $c-M$ relation in the different cosmological scenarios, using both dark matter only and hydrodynamical simulations. We find that the normalization of the $c-M$ relation is directly linked to $\sigma_{8}$ and the evolution of the density perturbations for $\Lambda$CDM, RP and SUGRA, while for EQp and EQn it depends also on the evolution of the linear density contrast. These differences in the $c-M$ relation provide another way to use galaxy clusters to constrain the underlying cosmology.

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.

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.

Galactic fountain flows in high-resolution hydrodynamical simulations of galaxy disks

Feedback from the most massive components of a young stellar cluster deeply affects the surrounding ISM driving an expanding over-pressured hot gas cavity in it. In spiral galaxies these structures may have sufficient energy to break the disk and eject large amount of material into the halo. The cycling of this gas, which eventually will fall back onto the disk, is known as galactic fountains. We aim at better understanding the dynamics of such fountain flow in a Galactic context, frame the problem in a more dynamic environment possibly learning about its connection and regulation to the local driving mechanism and understand its role as a metal diffusion channel. The interaction of the fountain with a hot corona is hereby analyzed, trying to understand the properties and evolution of the extraplanar material. We perform high resolution hydrodynamical simulations with the moving-mesh code AREPO to model the multi-phase ISM of a Milky Way type galaxy. A non-equilibrium chemical network is included to self consistently follow the evolution of the main coolants of the ISM. Spiral arm perturbations in the potential are considered so that large molecular gas structures are able to dynamically form here, self shielded from the interstellar radiation field. We model the effect of SN feedback from a new-born stellar cluster inside such a giant molecular cloud, as the driving force of the fountain. Passive Lagrangian tracer particles are used in conjunction to the SN energy deposition to model and study diffusion of freshly synthesized metals. We find that both interactions with hot coronal gas and local ISM properties and motions are equally important in shaping the fountain. We notice a bimodal morphology where most of the ejected gas is in a cold $10^4$ K clumpy state while the majority of the affected volume is occupied by a hot diffuse medium. While only about 20\% of the produced metals stay local, most of them quickly diffuse through this hot regime to great scales.

Observing the epoch of galaxy formation

Significant observational progress in addressing the question of the origin and early evolution of galaxies has been made in the past few years, allowing for direct comparison of the epoch when most of the stars in the universe were forming to prevailing theoretical models. There is currently broad consistency between theoretical expectations and the observations, but rapid improvement in the data will provide much more critical tests of theory in the coming years.

Magnetars: Properties, Origin and Evolution

Magnetars are neutron stars in which a strong magnetic field is the main energy source. About two dozens of magnetars, plus several candidates, are currently known in our Galaxy and in the Magellanic Clouds. They appear as highly variable X-ray sources and, in some cases, also as radio and/or optical pulsars. Their spin periods (2–12 s) and spin-down rates (∼10−13–10−10 s s−1) indicate external dipole fields of ∼1013−15 G, and there is evidence that even stronger magnetic fields are present inside the star and in non-dipolar magnetospheric components. Here we review the observed properties of the persistent emission from magnetars, discuss the main models proposed to explain the origin of their magnetic field and present recent developments in the study of their evolution and connection with other classes of neutron stars.

Galaxy threshing and the origin of ultra-compact dwarf galaxies in the Fornax cluster

A recent all-object spectroscopic survey centred on the Fornax cluster of galaxies has discovered a population of subluminous and extremely compact members, called 'ultra-compact dwarf' (UCD) galaxies. In order to clarify the origin of these objects, we have used self-consistent numerical simulations to study the dynamical evolution a nucleated dwarf galaxy would undergo if orbiting the centre of the Fornax cluster and suffering from its strong tidal gravitational field. We find that the outer stellar components of a nucleated dwarf are removed by the strong tidal field of the cluster, whereas the nucleus manages to survive as a result of its initially compact nature. The developed naked nucleus is found to have physical properties (e. g. size and mass) similar to those observed for UCDs. We also find that although this formation process does not have a strong dependence on the initial total luminosity of the nucleated dwarf, it does depend on the radial density profile of the dark halo in the sense that UCDs are less likely to be formed from dwarfs embedded in dark matter haloes with central 'cuspy' density profiles. Our simulations also suggest that very massive and compact stellar systems can be rapidly and efficiently formed in the central regions of dwarfs through the merging of smaller star clusters. We provide some theoretical predictions on the total number and radial number density profile of UCDs in a cluster and their dependencies on cluster masses.

Diffuse optical light in galaxy clusters. I. Abell 3888

We are undertaking a program to measure the characteristics of the intracluster light ( ICL; total flux, profile, color, and substructure) in a sample of 10 galaxy clusters with a range of cluster mass, morphology, and redshift. We present here the methods and results for the first cluster in that sample, A3888. We have identified an ICL component in A3888 in V and r that contains 13% +/- 5% of the total cluster light and extends to 700 h(70)(-1) kpc (similar to 0.3r(200)) from the center of the cluster. The ICL color in our smallest radial bin is V - r 0.3 +/- 0.1, similar to the central cluster elliptical galaxies. The ICL is redder than the galaxies at 400 h(70)(-1) kpc < r < 700 h(70)(-1) kpc, although the uncertainty in any one radial bin is high. Based on a comparison of V - r color with simple stellar models, the ICL contains a component that formed more than 7 Gyr ago ( at z less than 1) with a high-metallicity ( 1.0 Z(circle dot) < Z(ICL) less than or similar to 2.5 Z(circle dot)) and a more centralized component that contains stars formed within the past 5 Gyr ( at z similar to 1). The profile of the ICL can be roughly fitted by a shallow exponential in the outer regions and a steeper exponential in the central region. We also find a concentration of diffuse light around a small group of galaxies 1.4 h(70)(-1) Mpc from the center of the cluster. In addition, we find three low surface brightness features near the cluster center that are blue ( V - r 0.0) and contain a total flux of 0.1M*. Based on these observations and X-ray and galaxy morphology, we suggest that this cluster is entering a phase of significant merging of galaxy groups in the core, whereupon we expect the ICL fraction to grow significantly with the formation of a cD galaxy, as well as the infall of groups.

The effect of local galaxy surface density on star formation for HI selected galaxies

We present the result of investigations into two theories to explain the star formation rate (SFR)-density relationship. For regions of high galaxy density, either there are fewer star-forming galaxies or galaxies capable of forming stars are present but some physical process is suppressing their star formation. We use H I Parkes All-Sky Survey's (HIPASS) HI detected galaxies and infrared and radio fluxes to investigate SFRs and efficiencies with respect to local surface density. For nearby (vel < 10 000 km s(-1)) H I galaxies, we find a strong correlation between H I mass and SFR. The number of H I galaxies decreases with increasing local surface density. For H I galaxies (1000 < vel < 6000 km s(-1)), there is no significant change in the SFR or the efficiency of star formation with respect to local surface density. We conclude that the SFR-density relation is due to a decrease in the number of H I star-forming galaxies in regions of high galaxy density and not to the suppression of star formation.

The large peculiar velocity of the cD galaxy in Abell 3653

We present a catalogue of galaxies in Abell 3653 from observations made with the 2-degree field (2dF) spectrograph at the Anglo-Australian Telescope. Of the 391 objects observed, we find 111 are bona fide members of Abell 3653. We show that the cluster has a velocity of cz= 32 214 +/- 83 km s(-1) (z= 0.10 738 +/- 0.00 027), with a velocity dispersion typical of rich, massive clusters of sigma(cz)= 880(-54)(+66). We find that the cD galaxy has a peculiar velocity of 683 +/- 96 km s(-1) in the cluster rest frame - some 7 sigma away from the mean cluster velocity, making it one of the largest and most significant peculiar velocities found for a cD galaxy to date. We investigate the cluster for signs of substructure, but do not find any significant groupings on any length scale. We consider the implications of our findings on cD formation theories.

The Las Campanas/Anglo-Australian Telescope Rich Cluster Survey - III. Spectroscopic studies of X-ray bright galaxy clusters at z similar to 0.1

We present the analysis of the spectroscopic and photometric catalogues of 11 X-ray luminous clusters at 0.07 < z < 0.16 from the Las Campanas/Anglo-Australian Telescope Rich Cluster Survey. Our spectroscopic data set consists of over 1600 galaxy cluster members, of which two-thirds are outside r(200). These spectra allow us to assign cluster membership using a detailed mass model and expand on our previous work on the cluster colour-magnitude relation ( CMR) where membership was inferred statistically. We confirm that the modal colours of galaxies on the CMR become progressively bluer with increasing radius d( B - R)/dr(p) = - 0.011 +/- 0.003 and with decreasing local galaxy density d( B - R)/dlog ( Sigma)= - 0.062 +/- 0.009. Interpreted as an age effect, we hypothesize that these trends in galaxy colour should be reflected in mean H delta equivalent width. We confirm that passive galaxies in the cluster increase in Hd line strength as dH delta/dr(p) = 0.35 +/- 0.06. Therefore, those galaxies in the cluster outskirts may have younger luminosity-weighted stellar populations; up to 3 Gyr younger than those in the cluster centre assuming d( B - R)/dt = 0.03 mag per Gyr. A variation of star formation rate, as measured by [ O II]lambda 3727 angstrom, with increasing local density of the environment is discernible and is shown to be in broad agreement with previous studies from the 2dF Galaxy Redshift Survey and the Sloan Digital Sky Survey. We divide our spectra into a variety of types based upon the MORPHs classification scheme. We find that clusters at z similar to 0.1 are less active than their higher-redshift analogues: about 60 per cent of the cluster galaxy population is non-star forming, with a further 20 per cent in the post-starburst class and 20 per cent in the currently active class, demonstrating that evolution is visible within the past 2 - 3 Gyr. We also investigate unusual populations of blue and very red non-star forming galaxies and we suggest that the former are likely to be the progenitors of galaxies which will lie on the CMR, while the colours of the latter possibly reflect dust reddening. We show that the cluster galaxies at large radii consist of both backsplash ones and those that are infalling to the cluster for the first time. We make a comparison to the field population at z similar to 0.1 and examine the broad differences between the two populations. Individually, the clusters show significant variation in their galaxy populations which we suggest reflects their recent infall histories.

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.