A panchromatic view of the evolution of supermassive black holes

This PhD Thesis is devoted to the accurate analysis of the physical properties of Active Galactic Nuclei (AGN) and the AGN/host-galaxy interplay. Due to the broad-band AGN emission (from radio to hard X-rays), a multi-wavelength approach is mandatory. Our research is carried out over the COSMOS field, within the context of the XMM-Newton wide-field survey. To date, the COSMOS field is a unique area for comprehensive multi-wavelength studies, allowing us to define a large and homogeneous sample of QSOs with a well-sampled spectral coverage and to keep selection effects under control. Moreover, the broad-band information contained in the COSMOS database is well-suited for a detailed analysis of AGN SEDs, bolometric luminosities and bolometric corrections. In order to investigate the nature of both obscured (Type-2) and unobscured (Type-1) AGN, the observational approach is complemented with a theoretical modelling of the AGN/galaxy co-evolution. The X-ray to optical properties of an X-ray selected Type-1 AGN sample are discussed in the first part. The relationship between X-ray and optical/UV luminosities, parametrized by the spectral index αox, provides a first indication about the nature of the central engine powering the AGN. Since a Type-1 AGN outshines the surrounding environment, it is extremely difficult to constrain the properties of its host-galaxy. Conversely, in Type-2 AGN the host-galaxy light is the dominant component of the optical/near-IR SEDs, severely affecting the recovery of the intrinsic AGN emission. Hence a multi-component SED-fitting code is developed to disentangle the emission of the stellar populationof the galaxy from that associated with mass accretion. Bolometric corrections, luminosities, stellar masses and star-formation rates, correlated with the morphology of Type-2 AGN hosts, are presented in the second part, while the final part concerns a physically-motivated model for the evolution of spheroidal galaxies with a central SMBH. The model is able to reproduce two important stages of galaxy evolution, namely the obscured cold-phase and the subsequent quiescent hot-phase.

Solving the cooling flow problem through mechanical AGN feedback

A fundamental gap in the current understanding of collapsed structures in the universe concerns the thermodynamical evolution of the ordinary, baryonic component. Unopposed radiative cooling of plasma would lead to the cooling catastrophe, a massive inflow of condensing gas toward the centre of galaxies, groups and clusters. The last generation of multiwavelength observations has radically changed our view on baryons, suggesting that the heating linked to the active galactic nucleus (AGN) may be the balancing counterpart of cooling. In this Thesis, I investigate the engine of the heating regulated by the central black hole. I argue that the mechanical feedback, based on massive subrelativistic outflows, is the key to solving the cooling flow problem, i.e. dramatically quenching the cooling rates for several billion years without destroying the cool-core structure. Using an upgraded version of the parallel 3D hydrodynamic code FLASH, I show that anisotropic AGN outflows can further reproduce fundamental observed features, such as buoyant bubbles, cocoon shocks, sonic ripples, metals dredge-up, and subsonic turbulence. The latter is an essential ingredient to drive nonlinear thermal instabilities, which cause cold gas condensation, a residual of the quenched cooling flow and, later, fuel for the AGN feedback engine. The self-regulated outflows are systematically tested on the scales of massive clusters, groups and isolated elliptical galaxies: in lighter less bound objects the feedback needs to be gentler and less efficient, in order to avoid drastic overheating. In this Thesis, I describe in depth the complex hydrodynamics, involving the coupling of the feedback energy to that of the surrounding hot medium. Finally, I present the merits and flaws of all the proposed models, with a critical eye toward observational concordance.

Testing future weak lensing surveys through simulations of observations

Weak lensing experiments such as the future ESA-accepted mission Euclid aim to measure cosmological parameters with unprecedented accuracy. It is important to assess the precision that can be obtained in these measurements by applying analysis software on mock images that contain many sources of noise present in the real data. In this Thesis, we show a method to perform simulations of observations, that produce realistic images of the sky according to characteristics of the instrument and of the survey. We then use these images to test the performances of the Euclid mission. In particular, we concentrate on the precision of the photometric redshift measurements, which are key data to perform cosmic shear tomography. We calculate the fraction of the total observed sample that must be discarded to reach the required level of precision, that is equal to 0.05(1+z) for a galaxy with measured redshift z, with different ancillary ground-based observations. The results highlight the importance of u-band observations, especially to discriminate between low (z < 0.5) and high (z ~ 3) redshifts, and the need for good observing sites, with seeing FWHM < 1. arcsec. We then construct an optimal filter to detect galaxy clusters through photometric catalogues of galaxies, and we test it on the COSMOS field, obtaining 27 lensing-confirmed detections. Applying this algorithm on mock Euclid data, we verify the possibility to detect clusters with mass above 10^14.2 solar masses with a low rate of false detections.

The gaseous halo of the Milky Way

In the last decade, sensitive observations have revealed that disc galaxies are surrounded by multiphase gaseous halos produced by the circulation of gas from the discs to the environment and vice-versa. This Thesis is a study of the gaseous halo of the Milky Way carried out via the modelling of the HI emission and the available absorption-line data. We fitted simple kinematical models to the HI LAB Survey and found that the Galaxy has a massive (~3x10^8 Mo) HI halo extending a few kiloparsecs above the plane. This layer rotates more slowly than the disc and shows a global inflow motion, a kinematics similar to that observed in the HI halos of nearby galaxies. We built a dynamical model of the galactic fountain to reproduce the properties of this layer. In this model, fountain clouds are ejected from the disc by SN feedback and - as suggested by hydrodynamical simulations - triggers the cooling of coronal gas, which is entrained by the cloud wakes and accretes onto the disc when the clouds fall back. For a proper choice of the parameters, the model reproduces well the HI data and predicts an accretion of coronal gas onto the disc at a rate of 2 Mo/yr. We extended this model to the warm-hot component of the halo, showing that most of the ion absorption features observed towards background sources are consistent with being produced in the turbulent wakes that lag behind the fountain clouds. Specifically, the column densities, positions, and velocities of the absorbers are well reproduced by our model. Finally, we studied the gas content of galaxies extracted from a cosmological N-body+SPH simulation, and found that an HI halo with the forementioned properties is not observed, probably due ti the relatively low resolution of the simulations.

Resolved stellar populations in milky way and andromeda's new satellites

This thesis concerns the study of the variable stars and resolved stellar populations in four recently discovered dSphs, namely, Hercules and Ursa Major I (UMa I), which are UFD satellites of the MW; Andromeda XIX (And XIX) and Andromeda XXI (And XXI), which are satellites of M31. The main aim is to obtain detailed informations on the properties (age, metallicity, distance, and Oosterhoff type) of the stellar populations in these galaxies, to compare them with those of other satellites around the MW and M31, both ''classical'' dSphs and UFDs. The observables used to achieve these goals are the pulsating variables, especially the RR Lyrae stars, and the color magnitude diagram (CMD) of the resolved stellar populations. In particular, for UMa I, we combined B, V time-series observations from four different ground-based telescopes (Cassini, TLS, TT1 and Subaru) and for Hercules, we used archival data acquired with the Advanced Camera for Surveys (ACS) on board the HST. We used, instead B and V times-series photometry obtained with the Large Binocular Telescope (LBT) for And XIX and And XXI .

Constraints on neutrino mass fraction using Redshift Space Distortions

Redshift Space Distortions (RSD) are an apparent anisotropy in the distribution of galaxies due to their peculiar motion. These features are imprinted in the correlation function of galaxies, which describes how these structures distribute around each other. RSD can be represented by a distortions parameter $\beta$, which is strictly related to the growth of cosmic structures. For this reason, measurements of RSD can be exploited to give constraints on the cosmological parameters, such us for example the neutrino mass. Neutrinos are neutral subatomic particles that come with three flavours, the electron, the muon and the tau neutrino. Their mass differences can be measured in the oscillation experiments. Information on the absolute scale of neutrino mass can come from cosmology, since neutrinos leave a characteristic imprint on the large scale structure of the universe. The aim of this thesis is to provide constraints on the accuracy with which neutrino mass can be estimated when expoiting measurements of RSD. In particular we want to describe how the error on the neutrino mass estimate depends on three fundamental parameters of a galaxy redshift survey: the density of the catalogue, the bias of the sample considered and the volume observed. In doing this we make use of the BASICC Simulation from which we extract a series of dark matter halo catalogues, characterized by different value of bias, density and volume. This mock data are analysed via a Markov Chain Monte Carlo procedure, in order to estimate the neutrino mass fraction, using the software package CosmoMC, which has been conveniently modified. In this way we are able to extract a fitting formula describing our measurements, which can be used to forecast the precision reachable in future surveys like Euclid, using this kind of observations.

Galaxy evolution in the VIMOS public extragalactic redshift survey (VIPERS)

In this work I present the first measurements of the galaxy stellar mass function (GSMF) from the first public release of the VIPERS catalogue, containing ∼55,000 objects. First, I present the survey design, its scientific goal, the redshift measurements and validation. Then, I provide details about the estimate of galaxy stellar masses, star formation rates, and other physical quantities. I derive the GSMF of different galaxy types (e.g. active and passive galaxies) and as a function of the environment (defined through the local galaxy density contrast). These estimates represent new observational evidence useful to characterise the mechanism of galaxy evolution.

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.

Black Hole and Galaxy Growth over Cosmic Time: the Chandra COSMOS Legacy Survey

The study of supermassive black hole (SMBH) accretion during their phase of activity (hence becoming active galactic nuclei, AGN), and its relation to the host-galaxy growth, requires large datasets of AGN, including a significant fraction of obscured sources. X-ray data are strategic in AGN selection, because at X-ray energies the contamination from non-active galaxies is far less significant than in optical/infrared surveys, and the selection of obscured AGN, including also a fraction of heavily obscured AGN, is much more effective. In this thesis, I present the results of the Chandra COSMOS Legacy survey, a 4.6 Ms X-ray survey covering the equatorial COSMOS area. The COSMOS Legacy depth (flux limit f=2x10^(-16) erg/s/cm^(-2) in the 0.5-2 keV band) is significantly better than that of other X-ray surveys on similar area, and represents the path for surveys with future facilities, like Athena and X-ray Surveyor. The final Chandra COSMOS Legacy catalog contains 4016 point-like sources, 97% of which with redshift. 65% of the sources are optically obscured and potentially caught in the phase of main BH growth. We used the sample of 174 Chandra COSMOS Legacy at z>3 to place constraints on the BH formation scenario. We found a significant disagreement between our space density and the predictions of a physical model of AGN activation through major-merger. This suggests that in our luminosity range the BH triggering through secular accretion is likely preferred to a major-merger triggering scenario. Thanks to its large statistics, the Chandra COSMOS Legacy dataset, combined with the other multiwavelength COSMOS catalogs, will be used to answer questions related to a large number of astrophysical topics, with particular focus on the SMBH accretion in different luminosity and redshift regimes.

Suche nach Dunkler Materie in Galaxien und Galaxienhaufen mit dem Neutrinoteleskop IceCube

Trotz zahlreicher Hinweise auf die Existenz von Dunkler Materie, konnten entsprechende Teilchen bisher nicht nachgewiesen werden. Eine große Anzahl an Experimenten wird durchgeführt, um die Eigenschaften möglicher Kandidatenteilchen zu untersuchen. Eine Strategie ist die Suche nach einem Neutrinosignal aus den Annihilationen von schwach wechselwirkenden massiven Teilchen (WIMPs) in Regionen mit hoher Dichte von Dunkler Materie. Mögliche Zielobjekte dieser Suchen sind die Erde, die Sonne, das Zentrum und der Halo der Milchstraße sowie entfernte Objekte, die einen hohen Anteil Dunkler Materie aufweisen.In der vorliegenden Arbeit wird die erste Suche nach einem Neutrinosignal von Zwerggalaxien, größeren Galaxien und Galaxienhaufen beschrieben. Da kein Signal nachgewiesen wurde, konnten obere Grenzen auf den Annihilationsquerschnitt von WIMPs gesetzt werden. Die stärksten Grenzen wurden aus der Beobachtung des Virgo-Haufens unter der Annahme einer großen Signalverstärkung durch Unterstrukturen in der Dichteverteilung abgeleitet. Für WIMP-Massen oberhalb von einigen TeV ist das Ergebnis vergleichbar mit Grenzen, die aus der Suche mit Gammateleskopen abgeleitet wurden. Für den direkten Annihilationskanal in zwei Neutrinos konnte der Wirkungsquerschnitt stärker eingeschränkt werden, als in bisherigen Analysen.

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.

Multifrequency VLBA observations of two compact symmetric objects

In this thesis the results of the multifrequency VLBA observations of the GPS 1944+5448 and the HFP J0111+3906 are presented. They are compact objects smaller than about 100 pc, completely embedded in the host galaxy. The availability of multi-epoch VLBI observations spanning more than 10 years, allowed us to compute the hot spot advance speed in order to obtain the kinematic age of both sources. Both radio sources are young, in agreement with the idea that they are in an early evolutionary stage. The spectral analysis of each source component, such as the lobes, the hot spots, the core and the jets, making a comparison with the theoretical ones is described. In addition the physical parameters derived from VLBA images as the magnetic field, the luminosity, the energy and the ambient medium density of both sources are discussed.

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.