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.

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.

Is redshift-dependent evolution of galaxies a theoretical artifact?

The physical validity of the hypothesis of (redshift-dependent) luminosity evolution in galaxies is tested by statistical analysis of an intensively studied complete high-redshift sample of normal galaxies. The necessity of the evolution hypothesis in the frame of big-bang cosmology is confirmed at a high level of statistical significance; however, this evolution is quantitatively just as predicted by chronometric cosmology, in which there is no such evolution. Since there is no direct observational means to establish the evolution postulated in big-bang studies of higher-redshift galaxies, and the chronometric predictions involve no adjustable parameters (in contrast to the two in big-bang cosmology), the hypothesized evolution appears from the standpoint of conservative scientific methodology as a possible theoretical artifact.

The differential evolution of hot- and cold-dust in galaxies in the last 7.5 GYR

Tese de mestrado em Física, apresentada à Universidade de Lisboa, através da Faculdade de Ciências, 2016

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.

Examining the growth and evolution of magnetic fields in clusters of galaxies: a numerical perspective

Magnetic fields are ubiquitous in galaxy cluster atmospheres and have a variety of astrophysical and cosmological consequences. Magnetic fields can contribute to the pressure support of clusters, affect thermal conduction, and modify the evolution of bubbles driven by active galactic nuclei. However, we currently do not fully understand the origin and evolution of these fields throughout cosmic time. Furthermore, we do not have a general understanding of the relationship between magnetic field strength and topology and other cluster properties, such as mass and X-ray luminosity. We can now begin to answer some of these questions using large-scale cosmological magnetohydrodynamic (MHD) simulations of the formation of galaxy clusters including the seeding and growth of magnetic fields. Using large-scale cosmological simulations with the FLASH code combined with a simplified model of the acceleration of cosmic rays responsible for the generation of radio halos, we find that the galaxy cluster frequency distribution and expected number counts of radio halos from upcoming low-frequency sur- veys are strongly dependent on the strength of magnetic fields. Thus, a more complete understanding of the origin and evolution of magnetic fields is necessary to understand and constrain models of diffuse synchrotron emission from clusters. One favored model for generating magnetic fields is through the amplification of weak seed fields in active galactic nuclei (AGN) accretion disks and their subsequent injection into cluster atmospheres via AGN-driven jets and bubbles. However, current large-scale cosmological simulations cannot directly include the physical processes associated with the accretion and feedback processes of AGN or the seeding and merging of the associated SMBHs. Thus, we must include these effects as subgrid models. In order to carefully study the growth of magnetic fields in clusters via AGN-driven outflows, we present a systematic study of SMBH and AGN subgrid models. Using dark-matter only cosmological simulations, we find that many important quantities, such as the relationship between SMBH mass and galactic bulge velocity dispersion and the merger rate of black holes, are highly sensitive to the subgrid model assumptions of SMBHs. In addition, using MHD calculations of an isolated cluster, we find that magnetic field strengths, extent, topology, and relationship to other gas quantities such as temperature and density are also highly dependent on the chosen model of accretion and feedback. We use these systematic studies of SMBHs and AGN inform and constrain our choice of subgrid models, and we use those results to outline a fully cosmological MHD simulation to study the injection and growth of magnetic fields in clusters of galaxies. This simulation will be the first to study the birth and evolution of magnetic fields using a fully closed accretion-feedback cycle, with as few assumptions as possible and a clearer understanding of the effects of the various parameter choices.

Searching for star formation outside galaxies: Multiwavelength analysis of the intragroup medium of Hickson Compact Group 100

We use multiwavelength data (H I, FUV, NUV, R) to search for evidence of star formation in the intragroup medium of the Hickson Compact Group 100. We find that young star-forming regions are located in the intergalactic H I clouds of the compact group which extend to over 130 kpc away from the main galaxies. A tidal dwarf galaxy (TDG) candidate is located in the densest region of the H I tail, 61 kpc from the brightest group member and its age is estimated to be only 3.3 Myr. Fifteen other intragroup H II regions and TDG candidates are detected in the Galaxy Evolution Explorer (GALEX) FUV image and within a field 10' x 10' encompassing the H I tail. They have ages <200 Myr, H I masses of 10(9.2-10.4) M(circle dot), 0.001

Two formation channels of ultra-compact dwarf galaxies in Hickson compact groups

Context. The formation of ultra-compact dwarf galaxies (UCDs) is believed to be driven by interaction, and UCDs are abundant in the cores of galaxy clusters, environments that mark the end-point of galaxy evolution. Nothing is known about the properties of UCDs in compact groups of galaxies, environments where most of galaxy evolution and interaction is believed to occur and where UCDs in an intermediate stage in their evolution may be expected. Aims. The main goal of this study is to detect and characterize, for the first time, the UCD population of compact groups of galaxies. For that, two nearby groups in different evolutionary stages, HCG22 and HCG90, were targeted. Methods. We selected about 40 UCD candidates from pre-existing photometry of both groups, and obtained spectra of these candidates using the VLT FORS2 instrument in MXU mode. Archival HST/ACS imaging was used to measure their structural parameters. Results. We detect 16 and 5 objects belonging to HCG22 and HCG90, respectively, covering the magnitude range -10.0 > M(R) > -11.5 mag. Their integrated colours are consistent with old ages covering a broad range in metallicities (metallicities confirmed by the spectroscopic measurements). Photometric mass estimates put 4 objects in HCG90 and 9 in HCG22 in the mass range of UCDs (> 2 x 10(6) M(circle dot)) for an assumed age of 12Gyr. These UCDs are on average 2-3 times larger than the typical size of Galactic GCs, covering a range of 2 less than or similar to r(h) less than or similar to 21 pc. The UCDs in HCG22 are more concentrated around the central galaxy than in HCG90, at the 99% confidence level. They cover a broad range in [alpha/Fe] abundances from sub-to super-solar. The spectra of 3 UCDs (2 in HCG22, 1 in HCG90) show tentative evidence of intermediate age stellar populations. The clearest example is the largest and most massive UCD (similar to 10(7) M(circle dot)) in our sample, which is detected in HCG22. Its properties are most consistent with a stripped dwarf galaxy nucleus. We calculate the specific frequency (S(N)) of UCDs for both groups, finding that HCG22 has about three times higher S(N) than HCG90. Conclusions. The ensemble properties of the detected UCDs supports two co-existing formation channels: a star cluster origin (low-luminosity, compact sizes, old ages, super-solar alpha/Fe), and an origin as tidally stripped dwarf nuclei (more extended and younger stellar populations). Our results imply that the UCDs detected in both groups do not, in their majority, originate from relatively recent galaxy interactions. Most of the detected UCDs have likely been brought into the group along with their host galaxies.

Chemical evolution of the Galactic bulge as traced by microlensed dwarf and subgiant stars III. Detection of lithium in the metal-poor bulge dwarf MOA-2010-BLG-285S

Context. To study the evolution of Li in the Galaxy it is necessary to observe dwarf or subgiant stars. These are the only long-lived stars whose present-day atmospheric chemical composition reflects their natal Li abundances according to standard models of stellar evolution. Although Li has been extensively studied in the Galactic disk and halo, to date there has only been one uncertain detection of Li in an unevolved bulge star. Aims. Our aim with this study is to provide the first clear detection of Li in the Galactic bulge, based on an analysis of a dwarf star that has largely retained its initial Li abundance. Methods. We performed a detailed elemental abundance analysis of the bulge dwarf star MOA-2010-BLG-285S using a high-resolution and high signal-to-noise spectrum obtained with the UVES spectrograph at the VLT when the object was optically magnified during a gravitational microlensing event (visual magnification A similar to 550 during observation). The Li abundance was determined through synthetic line profile fitting of the (7)Li resonance doublet line at 670.8 nm. The results have been corrected for departures from LTE. Results. MOA-2010-BLG-285S is, at [Fe/H] = -1.23, the most metal-poor dwarf star detected so far in the Galactic bulge. Its old age (12.5 Gyr) and enhanced [alpha/Fe] ratios agree well with stars in the thick disk at similar metallicities. This star represents the first unambiguous detection of Li in a metal-poor dwarf star in the Galactic bulge. We find an NLTE corrected Li abundance of log epsilon(Li) = 2.16, which is consistent with values derived for Galactic disk and halo dwarf stars at similar metallicities and temperatures. Conclusions. Our results show that there are no signs of Li enrichment or production in the Galactic bulge during its earliest phases. Observations of Li in other galaxies (omega Cen) and other components of the Galaxy suggest further that the Spite plateau is universal.

Cluster and cluster galaxy evolution history from IR to X-ray observations of the young cluster RX J1257.2+4738 at z=0.866

Context. The cosmic time around the z similar to 1 redshift range appears crucial in the cluster and galaxy evolution, since it is probably the epoch of the first mature galaxy clusters. Our knowledge of the properties of the galaxy populations in these clusters is limited because only a handful of z similar to 1 clusters are presently known. Aims. In this framework, we report the discovery of a z similar to 0.87 cluster and study its properties at various wavelengths. Methods. We gathered X-ray and optical data (imaging and spectroscopy), and near and far infrared data (imaging) in order to confirm the cluster nature of our candidate, to determine its dynamical state, and to give insight on its galaxy population evolution. Results. Our candidate structure appears to be a massive z similar to 0.87 dynamically young cluster with an atypically high X-ray temperature as compared to its X-ray luminosity. It exhibits a significant percentage (similar to 90%) of galaxies that are also detected in the 24 mu m band. Conclusions. The cluster RXJ1257.2+4738 appears to be still in the process of collapsing. Its relatively high temperature is probably the consequence of significant energy input into the intracluster medium besides the regular gravitational infall contribution. A significant part of its galaxies are red objects that are probably dusty with on-going star formation.

An analytical approach to the dwarf galaxies cusp problem

Aims. An analytical solution for the discrepancy between observed core-like profiles and predicted cusp profiles in dark matter halos is studied. Methods. We calculate the distribution function for Navarro-Frenk-White halos and extract energy from the distribution, taking into account the effects of baryonic physics processes. Results. We show with a simple argument that we can reproduce the evolution of a cusp to a flat density profile by a decrease of the initial potential energy.

Galaxy threshing and the formation of ultracompact dwarf galaxies

Recent spectroscopic and morphological observational studies of galaxies around NGC 1399 in the Fornax Cluster have discovered several ultracompact dwarf galaxies with intrinsic sizes of similar to 100 pc and absolute B-band magnitudes ranging from -13 to -11 mag. In order to elucidate the origin of these enigmatic objects, we perform numerical simulations on the dynamical evolution of nucleated dwarf galaxies orbiting NGC 1399 and suffering from its strong tidal gravitational field. Adopting a plausible scaling relation for dwarf galaxies, we find that the outer stellar components of a nucleated dwarf are totally removed. This is due to them being tidally stripped over the course of several passages past the central region of NGC 1399. The nucleus, however, manages to survive. We also find that the size and luminosity of the remnant are similar to those observed for ultracompact dwarf galaxies, if the simulated precursor nucleated dwarf has a mass of similar to 10(8) M.. These results suggest that ultracompact dwarf galaxies could have previously been more luminous dwarf spheroidal or elliptical galaxies with rather compact nuclei.

The Besançon Galaxy Model renewed. I. Constraints on the Galactic thin disc evolution from Tycho data

Context. The understanding of Galaxy evolution can be facilitated by the use of population synthesis models, which allow to test hypotheses on the star formation history, star evolution, as well as chemical and dynamical evolution of the Galaxy. Aims. The new version of the Besanc¸on Galaxy Model (hereafter BGM) aims to provide a more flexible and powerful tool to investigate the Initial Mass Function (IMF) and Star Formation Rate (SFR) of the Galactic disc. Methods. We present a new strategy for the generation of thin disc stars which assumes the IMF, SFR and evolutionary tracks as free parameters. We have updated most of the ingredients for the star count production and, for the first time, binary stars are generated in a consistent way. We keep in this new scheme the local dynamical self-consistency as in Bienayme et al (1987). We then compare simulations from the new model with Tycho-2 data and the local luminosity function, as a first test to verify and constrain the new ingredients. The effects of changing thirteen different ingredients of the model are systematically studied. Results. For the first time, a full sky comparison is performed between BGM and data. This strategy allows to constrain the IMF slope at high masses which is found to be close to 3.0, excluding a shallower slope such as Salpeter"s one. The SFR is found decreasing whatever IMF is assumed. The model is compatible with a local dark matter density of 0.011 M pc−3 implying that there is no compelling evidence for significant amount of dark matter in the disc. While the model is fitted to Tycho2 data, a magnitude limited sample with V<11, we check that it is still consistent with fainter stars. Conclusions. The new model constitutes a new basis for further comparisons with large scale surveys and is being prepared to become a powerful tool for the analysis of the Gaia mission data.

The Besançon Galaxy Model renewed. I. Constraints on the Galactic thin disc evolution from Tycho data

Context. The understanding of Galaxy evolution can be facilitated by the use of population synthesis models, which allow to test hypotheses on the star formation history, star evolution, as well as chemical and dynamical evolution of the Galaxy. Aims. The new version of the Besanc¸on Galaxy Model (hereafter BGM) aims to provide a more flexible and powerful tool to investigate the Initial Mass Function (IMF) and Star Formation Rate (SFR) of the Galactic disc. Methods. We present a new strategy for the generation of thin disc stars which assumes the IMF, SFR and evolutionary tracks as free parameters. We have updated most of the ingredients for the star count production and, for the first time, binary stars are generated in a consistent way. We keep in this new scheme the local dynamical self-consistency as in Bienayme et al (1987). We then compare simulations from the new model with Tycho-2 data and the local luminosity function, as a first test to verify and constrain the new ingredients. The effects of changing thirteen different ingredients of the model are systematically studied. Results. For the first time, a full sky comparison is performed between BGM and data. This strategy allows to constrain the IMF slope at high masses which is found to be close to 3.0, excluding a shallower slope such as Salpeter"s one. The SFR is found decreasing whatever IMF is assumed. The model is compatible with a local dark matter density of 0.011 M pc−3 implying that there is no compelling evidence for significant amount of dark matter in the disc. While the model is fitted to Tycho2 data, a magnitude limited sample with V<11, we check that it is still consistent with fainter stars. Conclusions. The new model constitutes a new basis for further comparisons with large scale surveys and is being prepared to become a powerful tool for the analysis of the Gaia mission data.

Environments of Acitive Galaxies

This dissertation presents studies on the environments of active galaxies. Paper I is a case study of a cluster of galaxies containing BL Lac object RGB 1745+398. We measured the velocity dispersion, mass, and richness of the cluster. This was one of the most thorough studies of the environments of a BL Lac object. Methods used in the paper could be used in the future for studying other clusters as well. In Paper II we studied the environments of nearby quasars in the Sloan Digital Sky Survey (SDSS). We found that quasars have less neighboring galaxies than luminous inactive galaxies. In the large-scale structure, quasars are usually located at the edges of superclusters or even in void regions. We concluded that these low-redshift quasars may have become active only recently because the galaxies in low-density environments evolve later to the phase where quasar activity can be triggered. In Paper III we extended the analysis of Paper II to other types of AGN besides quasars. We found that different types of AGN have different large-scale environments. Radio galaxies are more concentrated in superclusters, while quasars and Seyfert galaxies prefer low-density environments. Different environments indicate that AGN have different roles in galaxy evolution. Our results suggest that activity of galaxies may depend on their environment on the large scale. Our results in Paper III raised questions of the cause of the environment-dependency in the evolution of galaxies. Because high-density large-scale environments contain richer groups and clusters than the underdense environments, our results could reflect smaller-scale effects. In Paper IV we addressed this problem by studying the group and supercluster scale environments of galaxies together. We compared the galaxy populations in groups of different richnesses in different large-scale environments. We found that the large-scale environment affects the galaxies independently of the group richness. Galaxies in low-density environments on the large scale are more likely to be star-forming than those in superclusters even if they are in groups with the same richness. Based on these studies, the conclusion of this dissertation is that the large-scale environment affects the evolution of galaxies. This may be caused by different “speed” of galaxy evolution in low and high-density environments: galaxies in dense environments reach certain phases of evolution earlier than galaxies in underdense environments. As a result, the low-density regions at low redshifts are populated by galaxies in earlier phases of evolution than galaxies in high-density regions.