Source author record

Francesco Shankar

Francesco Shankar appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

53works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

53 published item(s)

preprint2022arXiv

An Eddington ratio-driven origin for the ${\rm L}_{\rm X}-{\rm M}_{*}$ relation in quiescent and star forming active galaxies

A mild correlation exists in active galaxies between the mean black hole accretion, as traced by the mean X-ray luminosity $\left<{\rm L}_{\rm X}\right>$, and the host galaxy stellar mass M$_*$, characterised by a normalisation steadily decreasing with cosmic time and lower in more quiescent galaxies. We create comprehensive semi-empirical mock catalogues of active black holes to pin down which parameters control the shape and evolution of the $\left<{\rm L}_{\rm X}\right>-{\rm M}_*$ relation of X-ray detected active galaxies. We find that the normalisation of the $\left<{\rm L}_{\rm X}\right>-{\rm M}_*$ relation is largely independent of the fraction of active galaxies (the duty cycle), but strongly dependent on the mean Eddington ratio, when adopting a constant underlying M$_{\rm BH}-{\rm M}_*$ relation as suggested by observational studies. The data point to a decreasing mean Eddington ratio with cosmic time and with galaxy stellar mass at fixed redshift. Our data can be reproduced by black holes and galaxies evolving on similar M$_{\rm BH}-{\rm M}_*$ relations but progressively decreasing their average Eddington ratios, mean X-ray luminosities, and specific star formation rates, when moving from the starburst to the quiescent phase. Models consistent with the observed $\left<{\rm L}_{\rm X}\right>-{\rm M}_*$ relation and independent measurements of the mean Eddington ratios, are characterised by M$_{\rm BH}-{\rm M}_*$ relations lower than those derived from dynamically measured local black holes. Our results point to the $\left<{\rm L}_{\rm X}\right>-{\rm M}_*$ relation as a powerful diagnostic to: 1) probe black hole-galaxy scaling relations and the level of accretion onto black holes; 2) efficiently break the degeneracies between duty cycles and accretion rates in cosmological models of black holes.

preprint2022arXiv

Coincidence between morphology and star-formation activity through cosmic time: the impact of the bulge growth

The origin of the quenching in galaxies is still highly debated. Different scenarios and processes are proposed. We use multi-band (400-1600 nm) bulge-disc decompositions of massive galaxies in the redshift range 0<z<2 to explore the distribution and the evolution of galaxies in the log SFR-logM* plane as a function of the stellar mass-weighted bulge-to-total ratio (BTM) and also for internal galaxy components (bulge/disc) separately. We find evidence of a clear link between the presence of a bulge and the flattening of the Main Sequence in the high-mass end. All bulgeless galaxies (BTM<0.2) lie on the main-sequence, and there is little evidence of a quenching channel without bulge growth. Galaxies with a significant bulge component (BTM>0.2) are equally distributed in number between star forming and passive regions. The vast majority of bulges in the Main Sequence galaxies are quiescent, while star formation is localized in the disc component. Our current findings underline a strong correlation between the presence of the bulge and the star formation state of the galaxy. A bulge, if present, is often quiescent, independently of the morphology or the star formation activity of the host galaxy. Additionally, if a galaxy is quiescent, with a large probability, is hosting a bulge. Conversely, if the galaxy has a disky shape is highly probable to be star forming.

preprint2022arXiv

The low-end of the black hole mass function at cosmic dawn

Understanding the formation and growth of supermassive black holes (SMBHs) at high redshift represents a major challenge for theoretical models. In this work we investigate the early evolution of the first SMBHs by constraining their distribution in mass and luminosity at $z > 4$. In particular, we focus on the poorly explored low-mass end of the nuclear black hole (BH) distribution down to $z \simeq 4$, and explore its connection with the nature of the first BH seeds and the processes governing their mass growth. To this aim, we have developed CAT (Cosmic Archaeology Tool), a new semi-analytic model that describes the formation of the first stars and black holes in a self-consistent way and follows the co-evolution of nuclear BHs and their host galaxies for a representative population at $z > 4$. We find that current observational constraints favour models where the growth of BH seeds is Eddington limited and occurs at the Bondi-Hoyle-Lyttleton rate or where super-Eddington accretion occurs via a slim disk during gas rich galaxy mergers. The main difference between these two model variants lies at the low-end of the predicted mass and luminosity functions at $4 \le z \le 6$, where a clear gap appears in the first model, reflecting the stunted growth of light BH seeds formed as remnants of the first stars. Detecting this signature will be extremely challenging even for the future generation of space observatories, such as JWST, Athena and Lynx.

preprint2021arXiv

The Black Hole Mass Function Across Cosmic Times I. Stellar Black Holes and Light Seed Distribution

This is the first paper in a series aimed at modeling the black hole (BH) mass function, from the stellar to the intermediate to the (super)massive regime. In the present work we focus on stellar BHs and provide an ab-initio computation of their mass function across cosmic times. Specifically, we exploit the state-of-the-art stellar and binary evolutionary code \texttt{SEVN}, and couple its outputs with redshift-dependent galaxy statistics and empirical scaling relations involving galaxy metallicity, star-formation rate and stellar mass. The resulting relic mass function ${\rm d}N/{\rm d}V{\rm d}\log m_\bullet$ as a function of the BH mass $m_\bullet$ features a rather flat shape up to $m_\bullet\approx 50\, M_\odot$ and then a log-normal decline for larger masses, while its overall normalization at a given mass increases with decreasing redshift. We highlight the contribution to the local mass function from isolated stars evolving into BHs and from binary stellar systems ending up in single or binary BHs. We also include the distortion on the mass function induced by binary BH mergers, finding that it has a minor effect at the high-mass end. We estimate a local stellar BH relic mass density of $ρ_\bullet\approx 5\times 10^7\, M_\odot$ Mpc$^{-3}$, which exceeds by more than two orders of magnitude that in supermassive BHs; this translates into an energy density parameter $Ω_\bullet\approx 4\times 10^{-4}$, implying that the total mass in stellar BHs amounts to $\lesssim 1\%$ of the local baryonic matter. We show how our mass function for merging BH binaries compares with the recent estimates from gravitational wave observations by LIGO/Virgo, and discuss the possible implications for dynamical formation of BH binaries in dense environments like star clusters. [abridged]

preprint2021arXiv

The weak dependence of velocity dispersion on disk fractions, mass-to-light ratio and redshift: Implications for galaxy and black hole evolution

Velocity dispersion ($σ$) is a key driver for galaxy structure and evolution. We here present a comprehensive semi-empirical approach to compute $σ$ via detailed Jeans modelling assuming both a constant and scale-dependent mass-to-light ratio $M^*/L$. We compare with a large sample of local galaxies from MaNGA and find that both models can reproduce the Faber-Jackson (FJ) relation and the weak dependence of $σ$ on bulge-to-total ratio $B/T$ (for $B/T\gtrsim 0.25$). The dynamical-to-stellar mass ratio within $R\lesssim R_e$ can be fully accounted for by a gradient in $M^*/L$. We then build velocity dispersion evolutionary tracks $σ_{ap}[M^*,z]$ (within an aperture) along the main progenitor dark matter haloes assigning stellar masses, effective radii and Sersic indices via a variety of abundance matching and empirically motivated relations. We find: 1) clear evidence for downsizing in $σ_{ap}[M^*,z]$ along the progenitor tracks; 2) at fixed stellar mass $σ\propto(1+z)^{0.2-0.3}$ depending on the presence or not of a gradient in $M^*/L$. We extract $σ_{ap}[M^*,z]$ from the TNG50 hydrodynamic simulation and find very similar results to our models with constant $M^*/L$. The increasing dark matter fraction within $R_e$ tends to flatten the $σ_{ap}[M^*,z]$ along the progenitors at $z \gtrsim 1$ in constant $M^*/L$ models, while $σ_{ap}[M^*,z]$ have a steeper evolution in the presence of a stellar gradient. We then show that a combination of mergers and gas accretion are likely responsible for the constant or increasing $σ_{ap}[M^*,z]$ with time. Finally, our $σ_{ap}[M^*,z]$ are consistent with a nearly constant and steep $M_{bh}-σ$ relation at $z\lesssim 2$, with black hole masses derived from the $L_X-M^*$ relation.

preprint2020arXiv

Correlation of structure and stellar properties of galaxies in Stripe 82

Establishing a correlation (or lack thereof) between the bimodal colour distribution of galaxies and their structural parameters is crucial to understand the origin of bimodality. To achieve that, we have performed 2D mass-based structural decomposition (bulge+disc) of all disc galaxies (total$=$1263) in the Herschel imaging area of the Stripe 82 region using $K_s$ band images from the VICS82 survey. The scaling relations thus derived are found to reflect the internal kinematics and are employed in combination to select an indubitable set of classical and pseudo bulge hosting disc galaxies. The rest of the galaxies ($<20\%$) are marked as discs with "ambiguous" bulges. Pseudo and classical bulge disc galaxies exhibit clear bimodality in terms of all stellar parameters ($M_*$, sSFR, $r-K_s$). All pseudo bulge disc galaxies are blue and star-forming and all classical bulge disc galaxies are red and quiescent with less than $5\%$ digressions. Ambiguous bulge disc galaxies are intermittent to pseudo and classical bulge disc galaxies in the distribution of all structural and stellar parameters. $Δ$$\langleμ_{eb}\rangle$ - based on the placement of bulges on the Kormendy relation - is found to be the most efficient single structural indicator of both bulge type and stellar activity. The placement of ambiguous bulge disc galaxies on scaling relations and fundamental plane, in addition to their peculiar stellar properties suggest that they are dominantly a part of the green valley.

preprint2020arXiv

The case for the fundamental $M_{\rm BH}$-$σ$ relation

Strong scaling relations between host galaxy properties (such as stellar mass, bulge mass, luminosity, effective radius etc) and their nuclear supermassive black hole's mass point towards a close co-evolution. In this work, we first review previous efforts supporting the fundamental importance of the relation between supermassive black hole mass and stellar velocity dispersion ($M_{\rm BH}$-$σ_{\rm e}$). We then present further original work supporting this claim via analysis of residuals and principal component analysis applied to some among the latest compilations of local galaxy samples with dynamically measured supermassive black hole masses. We conclude with a review of the main physical scenarios in favour of the existence of a $M_{\rm BH}$-$σ_{\rm e}$ relation, with a focus on momentum-driven outflows.

preprint2020arXiv

The relationship between fine galaxy stellar morphology and star formation activity in cosmological simulations: a deep learning view

Hydrodynamical simulations of galaxy formation and evolution attempt to fully model the physics that shapes galaxies. The agreement between the morphology of simulated and real galaxies, and the way the morphological types are distributed across galaxy scaling relations are important probes of our knowledge of galaxy formation physics. Here we propose an unsupervised deep learning approach to perform a stringent test of the fine morphological structure of galaxies coming from the Illustris and IllustrisTNG (TNG100 and TNG50) simulations against observations from a subsample of the Sloan Digital Sky Survey. Our framework is based on PixelCNN, an autoregressive model for image generation with an explicit likelihood. We adopt a strategy that combines the output of two PixelCNN networks in a metric that isolates the fine morphological details of galaxies from the sky background. We are able to \emph{quantitatively} identify the improvements of IllustrisTNG, particularly in the high-resolution TNG50 run, over the original Illustris. However, we find that the fine details of galaxy structure are still different between observed and simulated galaxies. This difference is driven by small, more spheroidal, and quenched galaxies which are globally less accurate regardless of resolution and which have experienced little improvement between the three simulations explored. We speculate that this disagreement, that is less severe for quenched disky galaxies, may stem from a still too coarse numerical resolution, which struggles to properly capture the inner, dense regions of quenched spheroidal galaxies.

preprint2020arXiv

The Rise of Active Galactic Nuclei in the GAlaxy Evolution and Assembly semi-analytic model

We present a new implementation of the GAlaxy Evolution and Assembly (GAEA) semi-analytic model, that features an improved modelling of the process of cold gas accretion onto supermassive black hole (SMBHs), derived from both analytic arguments and high-resolution simulations. We consider different scenarios for the loss of angular momentum required for the available cold gas to be accreted onto the central SMBHs, and we compare different combinations of triggering mechanisms, including galaxy mergers and disc instabilities in star forming discs. We compare our predictions with the luminosity function (LF) observed for Active Galactic Nuclei (AGN) and we confirm that a non-instantaneous accretion timescale (either in the form of a low-angular momentum reservoir or as an assumed light curve evolution) is needed in order to reproduce the measured evolution of the AGN-LF and the so-called AGN-downsizing trend. Moreover, we also study the impact of AGN feedback, in the form of AGN-driven outflows, on the SF properties of model galaxies, using prescriptions derived both from empirical studies or from numerical experiments. We show that AGN-driven outflows are effective in suppressing the residual star formation rate in massive galaxies ($> 10^{11} M_\odot$) without changing their overall assembly history. These winds also affect the SFR of lower mass galaxies, resulting in a too large fraction of passive galaxies at $< 10^{10} M_\odot$. Finally, we study the Eddington ratio distribution as a function of SMBH mass, showing that only objects more massive than $10^8 M_\odot$ are already in a self-regulated state as inferred from observations.

preprint2019arXiv

Galaxy sizes and the galaxy-halo connection -- I: the remarkable tightness of the size distributions

The mass and structural assembly of galaxies is a matter of intense debate. Current theoretical models predict the existence of a linear relationship between galaxy size ($R_e$) and the host dark matter halo virial radius ($R_h$).\\ By making use of semi-empirical models compared to the size distributions of central galaxies from the Sloan Digital Sky Survey, we provide robust constraints on the normalization and scatter of the $R_e-R_h$ relation. We explore the parameter space of models in which the $R_e-R_h$ relation is mediated by either the spin parameter or the concentration of the host halo, or a simple constant the nature of which is in principle unknown. We find that the data require extremely tight relations for both early-type and late-type galaxies (ETGs,LTGs), especially for more massive galaxies. These constraints challenge models based solely on angular momentum conservation, which predict significantly wider distributions of galaxy sizes and no trend with stellar mass, if taken at face value. We discuss physically-motivated alterations to the original models that bring the predictions into better agreement with the data. We argue that the measured tight size distributions of SDSS disk galaxies can be reproduced by semi-empirical models in which the $R_e-R_h$ connection is mediated by the \emph{stellar} specific angular momenta $j_{star}.$ We find that current cosmological models of galaxy formation broadly agree with our constraints for LTGs, and justify the strong link between $R_e$ and $j_{star}$ that we propose, however the tightness of the $R_e-R_h$ relation found in such ab-initio theoretical models for ETGs is in tension with our semi-empirical findings.

preprint2019arXiv

No significant evolution of relations between Black hole mass and Galaxy total stellar mass up to z~2.5

We investigate the cosmic evolution of the ratio between black hole mass (MBH) and host galaxy total stellar mass (Mstellar) out to z~2.5 for a sample of 100 X-ray-selected moderate-luminosity, broad-line active galactic nuclei (AGNs) in the Chandra-COSMOS Legacy Survey. By taking advantage of the deep multi-wavelength photometry and spectroscopy in the COSMOS field, we measure in a uniform way the galaxy total stellar mass using a SED decomposition technique and the black hole mass based on broad emission line measurements and single-epoch virial estimates. Our sample of AGN host galaxies has total stellar masses of 10^10-12Msun, and black hole masses of 10^7.0-9.5Msun. Combining our sample with the relatively bright AGN samples from the literature, we find no significant evolution of the MBH-Mstellar relation with black hole-to-host total stellar mass ratio of MBH/Mstellar~0.3% at all redshifts probed. We conclude that the average black hole-to-host stellar mass ratio appears to be consistent with the local value within the uncertainties, suggesting a lack of evolution of the MBH-Mstellar relation up to z~2.5.

preprint2019arXiv

Probing black hole accretion tracks, scaling relations and radiative efficiencies from stacked X-ray active galactic nuclei

The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh-Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh-Mstar relation requires a mean radiative efficiency ~0.15, in line with theoretical expectations for accretion onto spinning black holes. However, matching the "raw" observed relation for inactive black holes requires a mean radiative efficiency around 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve a mean radiative efficiency ~0.12-0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized.

preprint2016arXiv

An assessment of the "too big to fail" problem for field dwarf galaxies in view of baryonic feedback effects

Recent studies have established that extreme dwarf galaxies --whether satellites or field objects-- suffer from the so called "too big to fail" (TBTF) problem. Put simply, the TBTF problem consists of the fact that it is difficult to explain both the measured kinematics of dwarfs and their observed number density within the LCDM framework. The most popular proposed solutions to the problem involve baryonic feedback processes. For example, reionization and baryon depletion can decrease the abundance of halos that are expected to host dwarf galaxies. Moreover, feedback related to star formation can alter the dark matter density profile in the central regions of low-mass halos. In this article we assess the TBTF problem for field dwarfs, taking explicitly into account the baryonic effects mentioned above. We find that 1) reionization feedback cannot resolve the TBTF problem on its own, because the halos in question are too massive to be affected by it, and that 2) the degree to which profile modification can be invoked as a solution to the TBTF problem depends on the radius at which galactic kinematics are measured. Based on a literature sample of about 90 dwarfs with interferometric observations in the 21cm line of atomic hydrogen (HI), we conclude that the TBTF problem persists despite baryonic effects. However, the preceding statement assumes that the sample under consideration is representative of the general population of field dwarfs. In addition, the unexplained excess of dwarf galaxies in LCDM could be as small as a factor of ~ 1.8, given the current uncertainties in the measurement of the galactic velocity function. Both of these caveats highlight the importance of upcoming uniform surveys with HI interferometers for advancing our understanding of the issue.

preprint2016arXiv

Setting firmer constraints on the evolution of the most massive, central galaxies from their local abundances and ages

There is still much debate surrounding how the most massive, central galaxies in the local universe have assembled their stellar mass, especially the relative roles of in-situ growth versus later accretion via mergers. In this paper, we set firmer constraints on the evolutionary pathways of the most massive central galaxies by making use of empirical estimates on their abundances and stellar ages. The most recent abundance matching and direct measurements strongly favour that a substantial fraction of massive galaxies with Mstar>3x10^11 Msun reside at the centre of clusters with mass Mhalo>3x10^13 Msun. Spectral analysis supports ages >10 Gyrs, corresponding to a formation redshift z_form >2. We combine these two pieces of observationally-based evidence with the mass accretion history of their host dark matter haloes. We find that in these massive haloes, the stellar mass locked up in the central galaxy is comparable to, if not greater than, the total baryonic mass at z_form. These findings indicate that either only a relatively minor fraction of their present-day stellar mass was formed in-situ at z_form, or that these massive, central galaxies form in the extreme scenario where almost all of the baryons in the progenitor halo are converted into stars. Interestingly, the latter scenario would not allow for any substantial size growth since the galaxy's formation epoch either via mergers or expansion. We show our results hold irrespective of systematic uncertainties in stellar mass, abundances, galaxy merger rates, stellar initial mass function, star formation rate and dark matter accretion histories.

preprint2016arXiv

The optical-UV emissivity of quasars: dependence on black hole mass and radio loudness

We analyzed a large sample of radio-loud and radio-quiet quasar spectra at redshift 1.0 < z < 1.2 to compare the inferred underlying quasar continuum slopes (after removal of the host galaxy contribution) with accretion disk models. The latter predict redder (decreasing) alpha_3000 continuum slopes (L_ν~ν^alpha at 3000Ang) with increasing black hole mass, bluer alpha_3000 with increasing luminosity at 3000Ang, and bluer alpha_3000 with increasing spin of the black hole, when all other parameters are held fixed. We find no clear evidence for any of these predictions in the data. In particular we find that: (i) alpha_3000 shows no significant dependence on black hole mass or luminosity. Dedicated Monte Carlo tests suggest that the substantial observational uncertainties in the black hole virial masses can effectively erase any intrinsic dependence of alpha_3000 on black hole mass, in line with some previous studies. (ii) The mean slope alpha_3000 of radio-loud sources, thought to be produced by rapidly spinning black holes, is comparable to, or even redder than, that of radio-quiet quasars. Indeed, although quasars appear to become more radio loud with decreasing luminosity, we still do not detect any significant dependence of alpha_3000 on radio loudness. The predicted mean alpha_3000 slopes tend to be bluer than in the data. Disk models with high inclinations and dust extinction tend to produce redder slopes closer to empirical estimates. Our mean alpha_3000 values are close to the ones independently inferred at z<0.5 suggesting weak evolution with redshift, at least for moderately luminous quasars.

preprint2015arXiv

A matter of measurement: rotation velocities and the velocity function of dwarf galaxies

The velocity function derived from large scale surveys can be compared with the predictions of LCDM cosmology, by matching the measured rotation velocities Vrot of galaxies to the maximum circular velocity of dark matter (DM) halos Vmax. For Vrot<50km/s, a major discrepancy arises between the observed and LCDM velocity functions. However, the manner in which different observational measures of Vrot are associated with Vmax is not straight forward in dwarf galaxies. We instead relate galaxies to DM halos using the empirical baryon- mass to halo-mass relation, and show that different observational measures of Vrot result in very different velocity functions. We show how the W50 velocity function, i.e. using the HI profile line width at 50% of peak HI flux to measure Vrot, can be reconciled with a LCDM cosmology. Our semi-empirical methodology allows us to determine the region of rotation curves that are probed by HI measurements (RHI), and shows that the Vrot of dwarfs are generally measured at a fraction of Rmax, explaining their tendency to have rising rotation curves. We provide fitting formulae for relating RHI and Reff (the effective radius) to the virial radius of DM halos. To continue to use velocity functions as a probe of LCDM cosmology, it is necessary to be precise about how the different measures of rotation velocity are probing the mass of the DM halos, dropping the assumption that any measure of rotational velocity can be equally used as a proxy for Vmax.

preprint2015arXiv

Avoiding progenitor bias: The structural and mass evolution of Brightest Group and Cluster Galaxies in Hierarchical models since z~1

The mass and structural evolution of massive galaxies is one of the hottest topics in galaxy formation. This is because it may reveal invaluable insights into the still debated evolutionary processes governing the growth and assembly of spheroids. However, direct comparison between models and observations is usually prevented by the so-called "progenitor bias", i.e., new galaxies entering the observational selection at later epochs, thus eluding a precise study of how pre-existing galaxies actually evolve in size. To limit this effect, we here gather data on high-redshift brightest group and cluster galaxies, evolve their (mean) host halo masses down to z=0 along their main progenitors, and assign as their "descendants" local SDSS central galaxies matched in host halo mass. At face value, the comparison between high redshift and local data suggests a noticeable increase in stellar mass of a factor of >2 since z~1, and of >2.5 in mean effective radius. We then compare the inferred stellar mass and size growth with those predicted by hierarchical models for central galaxies, selected at high redshifts to closely match the halo and stellar mass bins as in the data. Only hierarchical models characterized by very limited satellite stellar stripping and parabolic orbits are capable of broadly reproducing the stellar mass and size increase of a factor ~2-4 observed in cluster galaxies since z ~1. The predicted, average (major) merger rate since z~1 is in good agreement with the latest observational estimates.

preprint2015arXiv

Interpreting the possible break in the Black Hole - Bulge mass relation

Recent inspections of local available data suggest that the almost linear relation between the stellar mass of spheroids ($M_{\rm sph}$) and the mass of the super massive Black Holes (BHs) residing at their centres, shows a break below $M_{\rm sph} \sim 10^{10}\ {\rm M}_\odot$, with a steeper, about quadratic relation at smaller masses. We investigate the physical mechanisms responsible for the change in slope of this relation, by comparing data with the results of the semi-analytic model of galaxy formation MORGANA, which already predicted such a break in its original formulation. We find that the change of slope is mostly induced by effective stellar feedback in star-forming bulges. The shape of the relation is instead quite insensitive to other physical mechanisms connected to BH accretion such as disc instabilities, galaxy mergers, Active Galactic Nucleus (AGN) feedback, or even the exact modelling of accretion onto the BH, direct or through a reservoir of low angular momentum gas. Our results support a scenario where most stars form in the disc component of galaxies and are carried to bulges through mergers and disc instabilities, while accretion onto BHs is connected to star formation in the spheroidal component. Therefore, a model of stellar feedback that produces stronger outflows in star-forming bulges than in discs will naturally produce a break in the scaling relation. Our results point to a form of co-evolution especially at lower masses, below the putative break, mainly driven by stellar feedback rather than AGN feedback.

preprint2015arXiv

Is there a "too big to fail" problem in the field?

We use the Arecibo legacy fast ALFA (ALFALFA) 21cm survey to measure the number density of galaxies as a function of their rotational velocity, $V_\mathrm{rot,HI}$ (as inferred from the width of their 21cm emission line). Based on the measured velocity function we statistically connect galaxies with their host halo, via abundance matching. In a lambda cold dark matter ($Λ$CDM) cosmology, dwarf galaxies are expected to be hosted by halos that are significantly more massive than indicated by the measured galactic velocity; if smaller halos were allowed to host galaxies, then ALFALFA would measure a much higher galactic number density. We then seek observational verification of this predicted trend by analyzing the kinematics of a literature sample of gas-rich dwarf galaxies. We find that galaxies with $V_\mathrm{rot,HI} \lesssim 25$ $\mathrm{km} \, \mathrm{s}^{-1}$ are kinematically incompatible with their predicted $Λ$CDM host halos, in the sense that hosts are too massive to be accommodated within the measured galactic rotation curves. This issue is analogous to the "too big to fail" problem faced by the bright satellites of the Milky Way, but here it concerns extreme dwarf galaxies in the field. Consequently, solutions based on satellite-specific processes are not applicable in this context. Our result confirms the findings of previous studies based on optical survey data and addresses a number of observational systematics present in these works. Furthermore, we point out the assumptions and uncertainties that could strongly affect our conclusions. We show that the two most important among them -namely baryonic effects on the abundances of halos and on the rotation curves of halos- do not seem capable of resolving the reported discrepancy.

preprint2015arXiv

The morphologies of massive galaxies from z~3 - Witnessing the 2 channels of bulge growth

[abridged] We quantify the morphological evolution of z~0 massive galaxies ($M*/M_\odot\sim10^{11}$) from z~3 in the 5 CANDELS fields. The progenitors are selected using abundance matching techniques to account for the mass growth. The morphologies strongly evolve from z~3. At z<1, the population matches the massive end of the Hubble sequence, with 30% of spheroids, 50% of galaxies with equally dominant disk and bulge components and 20% of disks. At z~2-3 there is a majority of irregular systems (~60-70%) with still 30% of spheroids. We then analyze the SFRs, gas fractions and structural properties for the different morphologies independently. Our results suggest two distinct channels for the growth of bulges in massive galaxies. Around 30-40% were already bulges at z~2.5, with low average SFRs and gas-fractions (10-15%), high Sersic indices (n>3-4) and small effective radii ($R_e$~1 kpc) pointing towards an early formation through gas-rich mergers or VDI. Between z~ 2.5 and z~0, they rapidly increase their size by a factor of ~4-5, become all passive but their global morphology remains unaltered. The structural evolution is independent of the gas fractions, suggesting that it is driven by ex-situ events. The remaining 60% experience a gradual morphological transformation, from clumpy disks to more regular bulge+disks systems, essentially happening at z>1. It results in the growth of a significant bulge component (n~3) for 2/3 of the systems possibly through the migration of clumps while the remaining 1/3 keeps a rather small bulge (n~1.5-2). The transition phase between disturbed and relaxed systems and the emergence of the bulge is correlated with a decrease of the star formation activity and the gas fractions. The growth of the effective radii scales roughly with $H(z)^{-1}$ and it is therefore consistent with the expected growth of disks in galaxy haloes.

preprint2014arXiv

Environmental dependence of bulge-dominated galaxy sizes in hierarchical models of galaxy formation. Comparison with the local Universe

We compare state-of-the-art semi-analytic models of galaxy formation as well as advanced sub-halo abundance matching models with a large sample of early-type galaxies from SDSS at z < 0.3. We focus our attention on the dependence of median sizes of central galaxies on host halo mass. The data do not show any difference in the structural properties of early-type galaxies with environment, at fixed stellar mass. All hierarchical models considered in this work instead tend to predict a moderate to strong environmental dependence, with the median size increasing by a factor of about 1.5-3 when moving from low to high mass host haloes. At face value the discrepancy with the data is highly significant, especially at the cluster scale, for haloes above log Mhalo > 14. The convolution with (correlated) observational errors reduces some of the tension. Despite the observational uncertainties, the data tend to disfavour hierarchical models characterized by a relevant contribution of disc instabilities to the formation of spheroids, strong gas dissipation in (major) mergers, short dynamical friction timescales, and very short quenching timescales in infalling satellites. We also discuss a variety of additional related issues, such as the slope and scatter in the local size-stellar mass relation, the fraction of gas in local early-type galaxies, and the general predictions on satellite galaxies.

preprint2014arXiv

On the intermediate-redshift central stellar mass-halo mass relation, and implications for the evolution of the most massive galaxies since z~1

The stellar mass-halo mass relation is a key constraint in all semi-analytic, numerical, and semi-empirical models of galaxy formation and evolution. However, its exact shape and redshift dependence remain debated. Several recent works support a relation in the local Universe steeper than previously thought. Based on the comparisons with a variety of data on massive central galaxies, we show that this steepening holds up to z~1, for stellar masses Mstar>2e11 Msun. Specifically, we find significant evidence for a high-mass end slope of β>0.35-0.70, instead of the usual β~0.20-0.30 reported by a number of previous results. When including the independent constraints from the recent BOSS clustering measurements, the data, independent of any systematic errors in stellar masses, tend to favor a model with a very small scatter (< 0.15 dex) in stellar mass at fixed halo mass, in the redshift range z < 0.8 and for Mstar>3e11 Msun, suggesting a close connection between massive galaxies and host halos even at relatively recent epochs. We discuss the implications of our results with respect to the evolution of the most massive galaxies since z~1.

preprint2014arXiv

Unveiling the Intrinsic X-ray Properties of Broad Absorption Line Quasars with a Relatively Unbiased Sample

There is growing evidence of a higher intrinsic fraction of broad absorption line quasars (BALQSOs) than that obtained in optical surveys, on which most previous X-ray studies of BALQSOs have focused. Here we present Chandra observations of 18 BALQSOs at $z\sim2$, selected from a near-infrared (2MASS) sample, where the BALQSO fraction is likely to be close to the intrinsic fraction. We measure photon indices using the stacked spectra of the optically-faint ($i-K_s\geq 2.3$ mag) and optically-bright ($i-K_s < 2.3$ mag) samples to be $Γ\simeq 1.5$--$2.1$. We constrain their intrinsic column density by modelling the X-ray fractional hardness ratio, finding a mean column density of $3.5\times10^{22}$ \cmsq\ assuming neutral absorption. We incorporate SDSS optical measurements (rest frame UV) to study the broadband spectral index between the X-ray and UV bands, and compare this to a large sample of normal quasars. We estimate that the optically-faint BALQSOs are X-ray weaker than the optically-bright ones, and the entire sample of BALQSOs are intrinsically X-ray weak when compared to normal AGN. Correcting for magnification of X-ray emission via gravitational lensing by the central black hole viewed at large inclination angles makes these BALQSOs even more intrinsically X-ray weak. Finally, we estimate AGN kinetic feedback efficiencies of a few percent for an X-ray wind of $0.3c$ in high-ionization BALQSOs. Combined with energy carried by low-ionization BALQSOs and UV winds, the total kinetic energy in BALQSOs can be sufficient to provide AGN kinetic feedback required to explain the co-evolution between black holes and host galaxies.

preprint2013arXiv

Black Hole Demography: From scaling relations to models

In this contributed paper I review our current knowledge of the local Black Hole (BH) scaling relations, and their impact on the determination of the local BH mass function. I particularly emphasize the remaining systematic uncertainties impinging upon a secure determination of the BH mass function and how progress can be made. I then review and discuss the evidence for a different time evolution for separate BH-galaxy scaling relations, and how these independent empirical evidences can be reconciled with the overall evolution of the structural properties of the host galaxies. I conclude discussing BH demography in the context of semi-empirical continuity accretion models, as well as more complex evolutionary models, emphasizing the general constraints we can set on them.

preprint2013arXiv

No evidence for a dependence of the mass size relation of early-type galaxies on environment in the local Universe

The early--type galaxy (ETG) mass--size relation has been largely studied to understand how these galaxies have assembled their mass. One key observational result of the last years is that massive galaxies increased their size by a factor of a few at fixed stellar mass from $z\sim2$. Hierarchical models favor minor mergers as a plausible driver of this size growth. Some of these models, predict a significant environmental dependence in the sense that galaxies residing in more massive halos tend to be larger than galaxies in lower mass halos, at fixed stellar mass and redshift. At present, observational results of this environmental dependence have been contradictory. In this paper we revisit this issue in the local Universe, by investigating how the sizes of massive ETGs depend on large-scale environment using an updated and accurate sample of ETGs in different environments - field, group, clusters - from the Sloan Digital Sky Survey DR7. Our analysis does not show any significant environmental dependence of the sizes of central and satellites ETGs at fixed stellar mass at $z\sim0$. The size-mass relation of early-type galaxies at $z\sim0$ seems to be universal, i.e., independent of the mass of the host halo and of the position of the galaxy in that halo (central or satellite). The result is robust to different galaxy selections based on star formation, morphology or central density. In fact, considering our observational errors and the size of the sample, any size ratio larger than $30-40%$ between massive galaxies ($log(M_*/M_\odot)>11$) living in clusters and in the field can be ruled out at $3σ$ level.

preprint2013arXiv

Supermassive Black Holes and Their Host Galaxies - II. The correlation with near-infrared luminosity revisited

We present an investigation of the scaling relations between Supermassive Black Hole (SMBH) masses (Mbh), and their host galaxies' K-band bulge (Lbul) and total (Ltot) luminosities. The wide-field WIRCam imager at the Canada-France-Hawaii-Telescope (CFHT) was used to obtain the deepest and highest resolution near infrared images available for a sample of 35 galaxies with securely measured Mbh, selected irrespective of Hubble type. For each galaxy, we derive bulge and total magnitudes using a two-dimensional image decomposition code that allows us to account, if necessary, for large- and small-scale disks, cores, bars, nuclei, rings, envelopes and spiral arms. We find that the present-day Mbh-Lbul and Mbh-Ltot relations have consistent intrinsic scatter, suggesting that Mbh correlates equally well with bulge and total luminosity of the host. Our analysis provides only mild evidence of a decreased scatter if the fit is restricted to elliptical galaxies. The log-slopes of the Mbh-Lbul and Mbh-Ltot relations are 0.75+/-0.10 and 0.92+/-0.14, respectively. However, while the slope of the Mbh-Lbul relation depends on the detail of the image decomposition, the characterization of Mbh-Ltot does not. Given the difficulties and ambiguities of decomposing galaxy images into separate components, our results indicate that Ltot is more suitable as a tracer of SMBH mass than Lbul, and that the Mbh-Ltot relation should be used when studying the co-evolution of SMBHs and galaxies.

preprint2013arXiv

The Hot and Energetic Universe: A White Paper presenting the science theme motivating the Athena+ mission

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

preprint2013arXiv

The Hot and Energetic Universe: The formation and growth of the earliest supermassive black holes

A crucial challenge in astrophysics over the coming decades will be to understand the origins of supermassive black holes (SMBHs) that lie at the centres of most, if not all, galaxies. The processes responsible for the initial formation of these SMBHs and their early growth via accretion - when they are seen as Active Galactic Nuclei (AGN) - remain unknown. To address this challenge, we must identify low luminosity and obscured z>6 AGNs, which represent the bulk of early SMBH growth. Sensitive X-ray observations are a unique signpost of accretion activity, uncontaminated by star formation processes, which prevent reliable AGN identification at other wavelengths (e.g. optical, infrared). The Athena+ Wide Field Imager will enable X-ray surveys to be carried out two orders of magnitude faster than with Chandra or XMM-Newton, opening a new discovery space and identifying over 400 z>6 AGN, including obscured sources. Athena+ will also play a fundamental role to enhance the scientific return of future multiwavelength facilities that will probe the physical conditions within the host galaxies of early SMBHs, which is vital for understanding how SMBHs form, what fuels their subsequent growth, and to assess their impact on the early Universe. Follow-up of samples of z>6 galaxies with the Athena+ X-ray Integral Field Unit could also reveal the presence of highly obscured AGNs, thanks to the detection of strong iron lines. Thus, Athena+ will enable the first quantitative measurements of the extent and distribution of SMBH accretion in the early Universe.

preprint2012arXiv

Accretion-Driven Evolution of Black Holes: Eddington Ratios, Duty Cycles, and Active Galaxy Fractions

We develop semi-empirical models of the supermassive black hole and active galactic nucleus (AGN) populations, which incorporate the black hole growth implied by the observed AGN luminosity function assuming a radiative efficiency ε, and a distribution of Eddington ratios λ. By generalizing these continuity-equation models to allow a distribution P(λ|mbh,z) we are able to draw on constraints from observationally estimated P(λ) distributions and active galaxy fractions while accounting for the luminosity thresholds of observational samples. We consider models with a Gaussian distribution of log λ, and Gaussians augmented with a power-law tail to low λ. Within our framework, reproducing the high observed AGN fractions at low redshift requires a characteristic Eddington ratio λ_c that declines at late times, and matching observed Eddington ratio distributions requires P(λ) that broadens at low redshift. To reproduce the observed increase of AGN fraction with black hole or galaxy mass, we also require a λ_c that decreases with increasing black hole mass, reducing the AGN luminosity associated with the most massive black holes. Finally, achieving a good match to the high mass end of the local black hole mass function requires an increased radiative efficiency at high black hole mass. We discuss the potential impact of black hole mergers or a λ-dependent bolometric correction, and we compute evolutionary predictions for black hole and galaxy specific accretion rates. Despite the flexibility of our framework, no one model provides a good fit to all the data we consider.

preprint2012arXiv

Early-type galaxies at z = 1.3. I. The Lynx supercluster: cluster and groups at z=1.3. Morphology and color-magnitude relation

We confirm the detection of 3 groups in the Lynx supercluster, at z~1.3, and give their redshifts and masses. We study the properties of the group galaxies as compared to the central clusters, RXJ0849+4452 and RXJ0848+4453, selecting 89 galaxies in the clusters and 74 galaxies in the groups. We morphologically classify galaxies by visual inspection, noting that our early-type galaxy (ETG) sample would have been contaminated at the 30% -40% level by simple automated classification methods (e.g. based on Sersic index). In luminosity selected samples, both clusters and groups show high fractions of Sa galaxies. The ETG fractions never rise above ~50% in the clusters, which is low compared to the fractions observed in clusters at z~1. However, ETG plus Sa fractions are similar to those observed for ETGs in clusters at z~1. Bulge-dominated galaxies visually classified as Sas might also be ETGs with tidal features or merger remnants. They are mainly red and passive, and span a large range in luminosity. Their star formation seems to have been quenched before experiencing a morphological transformation. Because their fraction is smaller at lower redshifts, they might be the spiral population that evolves into ETGs. For mass-selected samples, the ETG fraction show no significant evolution with respect to local clusters, suggesting that morphological transformations occur at lower masses and densities. The ETG mass-size relation shows evolution towards smaller sizes at higher redshift in both clusters and groups, while the late-type mass-size relation matches that observed locally. The group ETG red sequence shows lower zero points and larger scatters than in clusters, both expected to be an indication of a younger galaxy population. The estimated age difference is small when compared to the difference in age at different galaxy masses.

preprint2012arXiv

Seeking for the leading actor on the cosmic stage: Galaxies versus Supermassive Black Holes

We present a Special Issue on the interplay of galaxies and Supermassive Black Holes (SMBHs) recently published in Advances in Astronomy. This is the introductory paper containing the motivation for this Special Issue together with a brief description of the articles which are part of the manuscript and the link to the entire book (http://www.hindawi.com/journals/aa/si/610485/). We hope this Special Issue will be useful for many astronomers who want to get an update on the current status of the AGN-Galaxy coevolution topic.

preprint2012arXiv

Size Evolution of Spheroids in a Hierarchical Universe

We present basic predictions of an updated version of the Munich semi-analytic hierarchical galaxy formation model that grows bulges via mergers and disk instabilities. Overall, we find that while spheroids below Ms ~ 10^11 Msun grow their sizes via a mixture of disk instability and mergers, galaxies above it mainly evolve via mergers. Including gas dissipation in major mergers, efficiently shrinks galaxies, especially those with final mass Ms < 10^11 Msun that are the most gas-rich, improving the match with different observables. We find that the predicted scatter in sizes at fixed stellar mass is still larger than the observed one by up to <40%. Spheroids are, on average, more compact at higher redshifts at fixed stellar mass, and at fixed redshift and stellar mass larger galaxies tend to be more starforming. More specifically, while for bulge-dominated galaxies the model envisages a nearly mass-independent decrease in sizes, the predicted size evolution for intermediate-mass galaxies is more complex. The z=2 progenitors of massive galaxies with mass around Ms and B/T>0.7 at z=0, are found to be mostly disc-dominated galaxies with a median B/T ~ 0.3, with only ~20% remaining bulge-dominated. The model also predicts that central spheroids living in more massive haloes tend to have larger sizes at fixed stellar mass. Including host halo mass dependence in computing velocity dispersions, allows the model to properly reproduce the correlations with stellar mass. We also discuss the fundamental plane, the correlations with galaxy age, the structural properties of pseudobulges, and the correlations with central black holes.

preprint2012arXiv

The evolution of the mass-size relation for early type galaxies from z~1 to the present: dependence on environment, mass-range and detailed morphology

[abridged] We study the dependence of the galaxy size evolution on morphology, stellar mass and large scale environment for a sample of 298 group and 384 field quiescent early-type galaxies from the COSMOS survey, selected from z~1 to the present, and with masses $log(M/M_\odot)>10.5$. The galaxy size growth depends on galaxy mass and early-type galaxy morphology, e.g., elliptical galaxies evolve differently than lenticular galaxies. At the low mass end -$10.5<Log(M/M_\odot)<11$, ellipticals do not show strong size growth from $z\sim1$ to the present (10% to 30% depending on the morphological classification). On the other end, massive ellipticals -log(M/M_\odot)>11.2$- approximately doubled their size. Interestingly, lenticular galaxies display different behavior: they appear more compact on average and they do show a size growth of \sim60% since z=1 independent of stellar mass. We compare our results with state-of-the art semi-analytic models. While major and minor mergers can account for most of the galaxy size growth, we find that with present data and the theoretical uncertainties in the modeling we cannot state clear evidence favoring either merger or mass loss via quasar and/or stellar winds as the primary mechanism driving the evolution. The galaxy mass--size relation and the size growth do not depend on environment in the halo mass range explored in this work (field to group mass $log(M_h/M_\odot)<14$), i.e., group and field galaxies follow the same trends, which is at variance with predictions from current hierarchical models that show a clear dependence of size growth on halo mass for massive ellipticals -$log(M_*/M_\odot)>11.2$.

preprint2011arXiv

A tale of two populations: the stellar mass of central and satellite galaxies

We develop a new empirical methodology to study the relation between the stellar mass of galaxies and the mass of their host subhaloes. Our approach is similar to abundance matching, and is based on assigning a stellar mass to each subhalo within a large cosmological N-body simulation. However, we significantly extend the method to include a different treatment for central and satellite galaxies, allowing the stellar mass of satellite galaxies to depend on both the host subhalo mass, and on the halo mass. Unlike in previous studies, our models are constrained by two observations: the stellar mass function of galaxies, and their auto-correlation functions (CFs). We present results for ~10^6 different successful models, showing that the uncertainty in the host subhalo mass reaches a factor of ~10 for a given stellar mass. Our results thus indicate that the host subhalo mass of central and satellite galaxies is poorly constrained by using information coming solely from the abundance and CFs of galaxies. In addition, we demonstrate that the fraction of stellar mass relative to the universal baryon fraction can vary between different models, and can reach ~0.6 for a specific population of low mass galaxies. We conclude that additional observational constraints are needed, in order to better constrain the mass relation between haloes and galaxies. These might include weak lensing, satellite kinematics, or measures of clustering other than the CFs. Moreover, CFs based on larger surveys with a better estimate of the systematic uncertainties are needed.

preprint2011arXiv

Black Hole growth and AGN obscuration by instability-driven inflows in high-redshift disk galaxies fed by cold streams

Disk galaxies at high redshift have been predicted to maintain high gas surface densities due to continuous feeding by intense cold streams leading to violent gravitational instability, transient features and giant clumps. Gravitational torques between the perturbations drive angular momentum out and mass in, and the inflow provides the energy for keeping strong turbulence. We use analytic estimates of the inflow for a self-regulated unstable disk at a Toomre stability parameter Q~1, and isolated galaxy simulations capable of resolving the nuclear inflow down to the central parsec. We predict an average inflow rate ~10 Msun/yr through the disk of a 10^11 Msun galaxy, with conditions representative of z~2 stream-fed disks. The inflow rate scales with disk mass and (1+z)^{3/2}. It includes clump migration and inflow of the smoother component, valid even if clumps disrupt. This inflow grows the bulge, while only a fraction ~ 10^-3 of it needs to accrete onto a central black hole (BH), in order to obey the observed BH-bulge relation. A galaxy of 10^11 Msun at z~2 is expected to host a BH of ~10^8 Msun, accreting on average with moderate sub-Eddington luminosity L_X ~ 10^42-43 erg/s, accompanied by brighter episodes when dense clumps coalesce. We note that in rare massive galaxies at z~6, the same process may feed 10^9 Msun BH at the Eddington rate. High central gas column densities can severely obscure AGN in high-redshift disks, possibly hindering their detection in deep X-ray surveys.

preprint2011arXiv

Constraints on star-formation driven galaxy winds from the mass-metallicity relation at z=0

We extend a chemical evolution model relating galaxy stellar mass and gas-phase oxygen abundance (the mass-metallicity relation) to explicitly consider the mass-dependence of galaxy gas fractions and outflows. Using empirically derived scalings of galaxy mass with halo virial velocity in conjunction with the most recent observations of z~0 total galaxy cold gas fractions and the mass-metallicity relation, we place stringent global constraints on the magnitude and scaling of the efficiency with which star forming galaxies expel metals. We demonstrate that under the assumptions that metal accretion is negligible and the stellar initial mass function does not vary, efficient outflows are required to reproduce the mass-metallicity relation; without winds, gas-to-stellar mass ratios >~ 0.3 dex higher than observed are needed. Moreover, z=0 gas fractions are low enough that while they have some effect on the magnitude of outflows required, the slope of the gas fraction--stellar mass relation does not strongly affect our conclusions on how the wind efficiencies must scale with galaxy mass. Despite systematic uncertainties in the normalization and slope of the mass-metallicity relation, we show that the metal expulsion efficiency zetaw=(Zw/Zg)etaw (where Zw is the wind metallicitiy and Zg is the interstellar medium metallicity) must be both high and scale steeply with mass. Specifically, we show that zetaw >> 1 and zetaw proportional to vvir^-3 or steeper. In contrast, momentum- or energy-driven outflow models suggest that etaw should scale as vvir^-1 or vvir^-2, respectively, implying that the Zw-Mstar relation should be shallower than the Zg-Mstar relation. [abridged]

preprint2011arXiv

Suzaku Observations of Three FeLoBAL QSOs, SDSS J0943+5417, J1352+4239, and J1723+5553

We present Suzaku observations of three iron low-ionization broad absorption line quasars (FeLoBALs). We detect J1723+5553 (3σ) in the observed 2--10 keV band, and constrain its intrinsic nh column density to nh > 6\times10^{23} \cmsq by modeling its X-ray hardness ratio. We study the broadband spectral index, aox, between the X-ray and UV bands by combining the X-ray measurements and the UV flux extrapolated from 2MASS magnitudes, assuming a range of intrinsic column densities, and then comparing the aox values for the three FeLoBALs with those from a large sample of normal quasars. We find that the FeLoBALs are consistent with the spectral energy distribution (SED) of normal quasars if the intrinsic nh column densities are nh > 7\times10^{23} \cmsq for J0943+5417, nh > 2\times10^{24} \cmsq for J1352+4293, and 6\times10^{23} < \nh < 3\times10^{24} \cmsq for J1723+5553. At these large intrinsic column densities, the optical depth from Thomson scattering can reach \sim6, which will significantly modulate the UV flux. Our results suggest that the X-ray absorbing material is located at a different place from the UV absorbing wind, likely between the X-ray and UV emitting regions. We find a significant kinetic feedback efficiency for FeLoBALs, indicating that the outflows are an important feedback mechanism in quasars.

preprint2011arXiv

The Intrinsic Fractions and Radio Properties of Low Ionization Broad Absorption Line Quasars

Low-ionization (MgII, FeII, FeIII) broad absorption line quasars (LoBALs) probe a relatively obscured quasar population, and could be at an early evolutionary stage for quasars. We study the intrinsic fractions of LoBALs using the SDSS, 2MASS, and FIRST surveys. We find that the LoBAL fractions of the near infra-red (NIR) and radio samples are approximately 5--7 times higher than those measured in the optical sample. This suggests that the fractions measured in the NIR and radio bands are closer to the intrinsic fractions of the populations, and that the optical fractions are significantly biased due to obscuration effects, similar to high-ionization broad absorption line quasars (HiBALs). We also find that the LoBAL fractions decrease with increasing radio luminosities, again, similar to HiBALs. In addition, we find tentative evidence for high fractions of LoBALs at high NIR luminosities, especially for FeLoBALs with a fraction of ~18 per cent at M_K_s < -31 mag. This population of NIR luminous LoBALs may be at an early evolutionary stage of quasar evolution. We use a two-component model of LoBALs including a pure geometric component and a luminosity dependent component at high NIR luminosities, and obtain better fits than those from a pure geometric model. Therefore, the LoBAL population can be modelled as a hybrid of both the geometric and evolutionary models, where the geometric component constitutes 3.4\pm0.3, 5.8\pm0.4, and 1.5\pm0.3 per cent of the quasar population for BI-LoBALs, AI-LoBALs, and FeLoBALs, respectively. Considering a population of obscured quasars that do not enter the SDSS survey, which could have a much higher LoBAL fraction, we expect that intrinsic fraction of LoBALs could be even higher.

preprint2010arXiv

Constraints on black hole duty cycles and the black hole-halo relation from SDSS quasar clustering

We use Shen et al.'s (2009) measurements of luminosity-dependent clustering in the SDSS Data Release 5 Quasar Catalog, at redshifts 0.4 < z < 2.5, to constrain the relation between quasar luminosity and host halo mass and to infer the duty cycle f_opt, the fraction of black holes that shine as optically luminous quasars at a given time. We assume a monotonic mean relation between quasar luminosity and host halo mass, with log-normal scatter Σ. For specified f_opt and Σ, matching the observed quasar space density determines the normalization of the luminosity-halo mass relation, from which we predict the clustering bias. The data show no change of bias between the faint and bright halves of the quasar sample but a modest increase in bias for the brightest 10%. At the mean redshift z=1.45 of the sample, the data can be well described either by models with small intrinsic scatter (Σ=0.1 dex) and a duty cycle f_opt=6*10^(-4) or by models with much larger duty cycles and larger values of the scatter. "Continuity equation" models of the black hole mass population imply f_opt > 2*10^(-3) in this range of masses and redshifts, and the combination of this constraint with the clustering measurements implies scatter Σ> 0.4 dex. These findings contrast with those inferred from the much stronger clustering of high-luminosity quasars at z~4, which require minimal scatter between luminosity and halo mass and duty cycles close to one.

preprint2010arXiv

Further constraining galaxy evolution models through the Size Function of SDSS Early-type galaxies

We discuss how the effective radius Phi(Re) function (ERF) recently worked out by Bernardi et al. (2009) represents a new testbed to improve the current understanding of Semi-analytic Models of Galaxy formation. In particular, we here show that a detailed hierarchical model of structure formation can broadly reproduce the correct peak in the size distribution of local early-type galaxies, although it significantly overpredicts the number of very compact and very large galaxies. This in turn is reflected in the predicted size-mass relation, much flatter than the observed one, due to too large (~3 kpc) low-mass galaxies (<10^11 \msun), and to a non-negligible fraction of compact (< 0.5-1 kpc) and massive galaxies (> 10^11 \msun). We also find that the latter discrepancy is smaller than previously claimed, and limited to only ultracompact (Re < 0.5 kpc) galaxies when considering elliptical-dominated samples. We explore several causes behind these effects. We conclude that the former problem might be linked to the initial conditions, given that large and low-mass galaxies are present at all epochs in the model. The survival of compact and massive galaxies might instead be linked to their very old ages and peculiar merger histories. Overall, knowledge of the galactic stellar mass {\em and} size distributions allows a better understanding of where and how to improve models.

preprint2010arXiv

On the radiative efficiencies, Eddington ratios, and duty cycles of luminous high-redshift quasars

We investigate the characteristic radiative efficiency ε, Eddington ratio λ, and duty cycle P_0 of high-redshift active galactic nuclei (AGN), drawing on measurements of the AGN luminosity function at z=3-6 and, especially, on recent measurements of quasar clustering at z=3-4.5 from the Sloan Digital Sky Survey. The free parameters of our models are ε, λ, and the normalization, scatter, and redshift evolution of the relation between black hole mass \mbh and halo virial velocity V_vir. We compute the luminosity function from the implied growth of the black hole mass function and the quasar correlation length from the bias of the host halos. We test our adopted formulae for the halo mass function and halo bias against measurements from the large N-body simulation developed by the MICE collaboration. The strong clustering of AGNs observed at z=3 and, especially, at z=4 implies that massive black holes reside in rare, massive dark matter halos. Reproducing the observed luminosity function then requires high efficiency εand/or low Eddington ratio λ, with a lower limit (based on 2σagreement with the measured z=4 correlation length) ε> 0.7λ/(1+0.7λ), implying ε> 0.17 for λ> 0.25. Successful models predict high duty cycles, P_0~0.2, 0.5, and 0.9 at z=3.1, 4.5 and 6, respectively, and they require that the fraction of halo baryons locked in the central black hole is much larger than the locally observed value. The rapid drop in the abundance of the massive and rare host halos at z>7 implies a proportionally rapid decline in the number density of luminous quasars, much stronger than simple extrapolations of the z=3-6 luminosity function would predict. (abridged)

preprint2009arXiv

Binary Quasars at High Redshift I: 24 New Quasar Pairs at z ~ 3-4

The clustering of quasars on small scales yields fundamental constraints on models of quasar evolution and the buildup of supermassive black holes. This paper describes the first systematic survey to discover high redshift binary quasars. Using color-selection and photometric redshift techniques, we searched 8142 deg^2 of SDSS imaging data for binary quasar candidates, and confirmed them with follow-up spectroscopy. Our sample of 27 high redshift binaries (24 of them new discoveries) at redshifts 2.9 < z < 4.3 with proper transverse separations 10 kpc < R_{\perp} < 650 kpc increases the number of such objects known by an order of magnitude. Eight members of this sample are very close pairs with R_{\perp} < 100 kpc, and of these close systems four are at z > 3.5. The completeness and efficiency of our well-defined selection algorithm are quantified using simulated photometry and we find that our sample is ~ 50% complete. Our companion paper uses this knowledge to make the first measurement of the small scale clustering (R < 1 Mpc/h comoving) of high-redshift quasars. High redshift binaries constitute exponentially rare coincidences of two extreme (M >~ 10^9 Msun) supermassive black holes. At z ~ 4 there is about one close binary per 10 Gpc^3, thus these could be the highest sigma peaks, the analogs of superclusters, in the early Universe.

preprint2009arXiv

Binary Quasars at High Redshift II: Sub-Mpc Clustering at z ~ 3-4

We present measurements of the small-scale (0.1<~ r <~ 1 Mpc/h) quasar two-point correlation function at z>2.9, for a flux-limited (i<21) sample of 15 binary quasars compiled by Hennawi et al. (2009). The amplitude of the small-scale clustering increases from z ~ 3 to z ~ 4. The small-scale clustering amplitude is comparable to or lower than power-law extrapolations (with slope gamma=2) from the large-scale correlation function of the i<20.2 quasar sample from the Sloan Digital Sky Survey. Using simple prescriptions relating quasars to dark matter halos, we model the observed small-scale clustering with halo occupation models. Reproducing the large-scale clustering amplitude requires that the active fraction of the black holes in the central galaxies of halos is near unity, but the level of small-scale clustering favors an active fraction of black holes in satellite galaxies 0.1 <~ f_s <~ 0.5 at z >~ 3.

preprint2009arXiv

On merger bias and the clustering of quasars

We use the large catalogues of haloes available for the Millennium Simulation to test whether recently merged haloes exhibit stronger large-scale clustering than other haloes of the same mass. This effect could help to understand the very strong clustering of quasars at high redshift. However, we find no statistically significant excess bias for recently merged haloes over the redshift range 2 < z < 5, with the most massive haloes showing an excess of at most ~5%. We also consider galaxies extracted from a semianalytic model built on the Millennium Simulation. At fixed stellar mass, we find an excess bias of ~ 20-30% for recently merged objects, decreasing with increasing stellar mass. The fact that recently-merged galaxies are found in systematically more massive haloes than other galaxies of the same stellar mass accounts for about half of this signal, and perhaps more for high-mass galaxies. The weak merger bias of massive systems suggests that objects of merger-driven nature, such as quasars, do not cluster significantly differently than other objects of the same characteristic mass. We discuss the implications of these results for the interpretation of clustering data with respect to quasar duty cycles, visibility times, and evolution in the black hole-host mass relation.

preprint2009arXiv

Sizes and ages of SDSS ellipticals: Comparison with hierarchical galaxy formation models

In a sample of about 45,700 early-type galaxies extracted from SDSS, we find that the shape, normalization, and dispersion around the mean size-stellar mass relation is the same for young and old systems, provided the stellar mass is greater than 3*10^10 Msun. This is difficult to reproduce in pure passive evolution models, which generically predict older galaxies to be much more compact than younger ones of the same stellar mass. However, this aspect of our measurements is well reproduced by hierarchical models of galaxy formation. Whereas the models predict more compact galaxies at high redshifts, subsequent minor, dry mergers increase the sizes of the more massive objects, resulting in a flat size-age relation at the present time. At lower masses, the models predict that mergers are less frequent, so that the expected anti-correlation between age and size is not completely erased. This is in good agreement with our data: below 3*10^10 Msun, the effective radius R_e is a factor of ~2 lower for older galaxies. These successes of the models are offset by the fact that the predicted sizes have other serious problems, which we discuss.

preprint2009arXiv

The age dependence of the size-stellar mass relation and some implications

We use a sample of about 48,000 SDSS early-type galaxies to show that older galaxies have smaller half-light radii re and larger velocity dispersions sigma than younger ones of the same stellar mass Mstar. We use the age-corrected luminosity Lcorr as a proxy for Mstar to minimize biases: below Lcorr~1e11 Lsun, galaxies with age ~11 Gyrs have re smaller by 40% and sigma larger by 25%, compared to galaxies that are ~4 Gyrs younger. The sizes and velocity dispersions of more luminous galaxies vary by less than 15%, whatever their age, a challenge for current galaxy formation models. A closer check reveals that the lowering in the dispersion is caused by older galaxies that show a significant departure from the re--Lcorr and sigma--Lcorr relations at high Lcorr. Such features might find an explanation in models where more massive galaxies undergo more minor mergers than less massive galaxies at late times, thus causing a break in the homology. In terms of the Fundamental Plane of early-type galaxies, the data indicate that all galaxies show a significant and similar increase in the dynamical-to-stellar mass ratio with increasing mass, independent of their age. However, older galaxies have smaller Mdyn/Mstar ratios than objects which formed more recently. These findings may suggest that lower mass galaxies and, at fixed stellar mass, higher redshift galaxies, formed from gas-richer progenitors, thus underwent more dissipation and contraction.

preprint2009arXiv

The Demography of Super-Massive Black Holes: Growing Monsters at the Heart of Galaxies

Supermassive black holes (BHs) appear to be ubiquitous at the center of all galaxies which have been observed at high enough sensitivities and resolution with the Hubble Space Telescope. Their masses are found to be tightly linked with the masses and velocity dispersions of their host galaxies. On the other hand, BHs are widely held to constitute the central engines of quasars and active galactic nuclei (AGN) in general. It is however still unclear how BHs have grown, and whether they have co-evolved with their hosts. In this Review I discuss how, in ways independent of specific models, constraints on the growth history of BHs and their host galaxies have been set by matching the statistics of local BHs to the emissivity, number density, and clustering properties of AGNs at different cosmological epochs. I also present some new results obtained through a novel numerical code which evolves the BH mass function and clustering adopting broad distributions of Eddington ratios. I finally review BH evolution in a wider cosmological context, connecting BH growth to galaxy evolution.

preprint2009arXiv

The relative growth of Optical and Radio Quasars in SDSS

We cross-correlate the SDSS DR3 quasar sample with FIRST and the Vestergaard et al. black hole (BH) mass sample to compare the mean accretion histories of optical and radio quasars. We find significant statistical evidence that radio quasars have a higher mean Eddington ratio Lambda at z > 2 with respect to optical quasars, while the situation is clearly reverse at z < 1. At z > 2 radio quasars happen to be less massive than optical quasars; however, as redshift decreases radio quasars appear in increasingly more massive BHs with respect to optical quasars. These two trends imply that radio sources are not a mere random subsample of optical quasars. No clear correlation between radio activity and BH mass and/or accretion rate is evident from our data, pointing to other BH properties, possibly the spin, as the driver of radio activity. We have checked that our main results do not depend on any evident bias. We perform detailed modelling of reasonable accretion histories for optical and radio quasars, finding that radio quasars grow by a factor of a few, at the most, since z ~ 4. The comparison between the predicted mass function of active radio quasars and the observed optical luminosity function of radio quasars, implies a significantly lower probability for lower mass BHs to be radio loud at all epochs, in agreement with what is observed in the local universe.

preprint2008arXiv

Quasar Clustering from SDSS DR5: Dependences on Physical Properties

Using a homogenous sample of 38,208 quasars with a sky coverage of $4000 {\rm deg^2}$ drawn from the SDSS Data Release Five quasar catalog, we study the dependence of quasar clustering on luminosity, virial black hole mass, quasar color, and radio loudness. At $z<2.5$, quasar clustering depends weakly on luminosity and virial black hole mass, with typical uncertainty levels $\sim 10%$ for the measured correlation lengths. These weak dependences are consistent with models in which substantial scatter between quasar luminosity, virial black hole mass and the host dark matter halo mass has diluted any clustering difference, where halo mass is assumed to be the relevant quantity that best correlates with clustering strength. However, the most luminous and most massive quasars are more strongly clustered (at the $\sim 2σ$ level) than the remainder of the sample, which we attribute to the rapid increase of the bias factor at the high-mass end of host halos. We do not observe a strong dependence of clustering strength on quasar colors within our sample. On the other hand, radio-loud quasars are more strongly clustered than are radio-quiet quasars matched in redshift and optical luminosity (or virial black hole mass), consistent with local observations of radio galaxies and radio-loud type 2 AGN. Thus radio-loud quasars reside in more massive and denser environments in the biased halo clustering picture. Using the Sheth et al.(2001) formula for the linear halo bias, the estimated host halo mass for radio-loud quasars is $\sim 10^{13} h^{-1}M_\odot$, compared to $\sim 2\times 10^{12} h^{-1}M_\odot$ for radio-quiet quasar hosts at $z\sim 1.5$.

preprint2008arXiv

Self-Consistent Models of the AGN and Black Hole Populations: Duty Cycles, Accretion Rates, and the Mean Radiative Efficiency

We construct evolutionary models of the populations of AGN and supermassive black holes, in which the black hole mass function grows at the rate implied by the observed luminosity function, given assumptions about the radiative efficiency and the Eddington ratio. We draw on a variety of recent X-ray and optical measurements to estimate the bolometric AGN luminosity function and compare to X-ray background data and the independent estimate of Hopkins et al. (2007) to assess remaining systematic uncertainties. The integrated AGN emissivity closely tracks the cosmic star formation history, suggesting that star formation and black hole growth are closely linked at all redshifts. Observational uncertainties in the local black hole mass function remain substantial, with estimates of the integrated black hole mass density ρ_BH spanning the range 3-5.5x10^5 Msun/Mpc^3. We find good agreement with estimates of the local mass function for a reference model where all active black holes have efficiency \eps=0.065 and L_bol/L_Edd~0.4. In this model, the duty cycle of 10^9 Msun black holes declines from 0.07 at z=3 to 0.004 at z=1 and 0.0001 at z=0. The decline is shallower for less massive black holes, a signature of "downsizing" evolution in which more massive black holes build their mass earlier. The predicted duty cycles and AGN clustering bias in this model are in reasonable accord with observational estimates. If the typical Eddington ratio declines at z<2, then the "downsizing" of black hole growth is less pronounced. Matching the integrated AGN emissivity to the local black hole mass density implies \eps=0.075 (ρ_BH/4.5x10^5 Msun/Mpc^3)^{-1} for our standard luminosity function estimate (25% higher for Hopkins et al.'s), lower than the values \eps=0.16-0.20 predicted by MHD simulations of disk accretion.

preprint2008arXiv

The Evolution of the M_BH-sigma relation Inferred from the Age Distribution of Local Early-Type Galaxies and AGN Evolution

We utilize the local velocity dispersion function (VDF) of spheroids, together with their inferred age--distributions, to predict the VDF at higher redshifts (0<z<6), under the assumption that (i) most of the stars in each nearby spheroid formed in a single episode, and (ii) the velocity dispersion sigma remained nearly constant afterward. We assume further that a supermassive black hole (BH) forms concurrently with the stars, and within ~1 Gyr of the formation of the potential well of the spheroid, and that the relation between the mass of the BH and host velocity dispersion maintains the form M_BH ~ sigma^{beta} with beta~4, but with the normalization allowed to evolve with redshift as ~(1+z)^{alpha}. We compute the BH mass function associated with the VDF at each redshift, and compare the accumulated total BH mass density with that inferred from the integrated quasar luminosity function (LF; the so--called Soltan argument). This comparison is insensitive to the assumed duty cycle or Eddington ratio of quasar activity, and we find that the match between the two BH mass densities favors a relatively mild redshift evolution, with alpha ~ 0.26, with a positive evolution as strong as alpha>1.3 excluded at the 99% confidence level. A direct match between the characteristic BH mass in the VDF--based and quasar LF--based BH mass functions also yields a mean Eddington ratio of lambda ~ 0.5-1 that is roughly constant within 0<z<3. A strong positive evolution in the M_BH-sigma relation is still allowed by the data if galaxies increase, on average, their velocity dispersions since the moment of formation, due to dissipative processes. If we assume that the mean velocity dispersion of the host galaxies evolves as sigma(z)=sigma(0)*(1+z)^{-gamma}, we find a lower limit of gamma>0.23 for alpha>1.5. abridged

preprint2007arXiv

2MASS Reveals a Large Intrinsic Fraction of BALQSOs

The intrinsic fraction of broad absorption line quasars (BALQSOs) is important in constraining geometric and evolutionary models of quasars. We present the fraction of BALQSOs in 2MASS detected quasars within the SDSS DR3 sample in the redshift range of 1.7 < z < 4.38. The fraction of BALQSOs is 40.4^{+3.4}_{-3.3}% in the 2MASS 99% database K_s band completeness sample, and 38.5^{+1.7}_{-1.7}% in the larger 2MASS sample extending below the completeness limit. These fractions are significantly higher than the 26% reported in the optical bands for the same parent sample. We also present the fraction of BALQSOs as functions of apparent magnitudes, absolute magnitudes, and redshift in the 2MASS and SDSS bands. The 2MASS fractions are consistently higher than the SDSS fractions in every comparison, and the BALQSO fractions steadily increase with wavelength from the SDSS u to the 2MASS K_s bands. Furthermore, the i - K_s color distributions of BALQSOs and non-BALQSOs indicate that BALQSOs are redder than non-BALQSOs, with a K-S test probability of 2e-12. These results are consistent with the spectral difference between BALQSOs and non-BALQSOs including both the absorption troughs and dust extinction in BALQSOs, which leads to significant selection biases against BALQSOs in the optical bands. Using a simple simulation incorporating the luminosity function of quasars and the amount of obscuration for BALQSOs, we simultaneously fit the BALQSO fractions in the SDSS and 2MASS bands. We obtain a true BALQSO fraction of 43\pm2% for luminous quasars (M_{K_s} \lesssim -30.1 mag).

preprint2007arXiv

On the faint end of the high redshift AGN luminosity function

Using the results of recent optical surveys we conclude that the {\it non}-detection of quasars down to faint magnitudes implies a significant flattening of the high redshift (z~6) optical active galactic nuclei (AGN) luminosity function for M_{1450}>-24.7. We find that all the data are consistent with a faint-end slope for the optical AGN luminosity function of β=-2.2 and β=-2.8, at the 90% and 99% confidence level respectively, flatter than the bright-end slope of β'~ -3.2. We also show that X-ray deep surveys have probed even fainter magnitudes than the optical ones yielding more significant constraints on the shallow faint-end slope of the optical luminosity function. The inclusion of Type II AGN candidates, detected in the Chandra deep fields, hints towards an higher normalization for the total AGN luminosity function, if these sources lie at 5<z <6.5. We then discuss simple theoretical models of AGN formation and evolution in the context of cold dark matter cosmology. The comparison with the total AGN luminosity function favors a redshift-dependent relation between black hole and dark matter halo masses of the type M_BH~M_ halo^α, with 1.3< α< 1.7, compatible with independent studies from statistical analysis and rotation curve measurements. Finally we compute the quasar contribution to reionization to be <9% at z~6, up to ~30% when integrated within 5.5< z<6.5, significantly smaller than that from galaxies.