Source author record

Michaela Hirschmann

Michaela Hirschmann 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

29works
2topics
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

29 published item(s)

preprint2026arXiv

From "The Cliff" to "Virgil": Mapping the Spectral Diversity of Little Red Dots with JWST/NIRSpec

One of JWST's most unexpected discoveries is the emergence of "Little Red Dots'' (LRDs): compact sources at $z \gtrsim 3$ with blue rest-frame UV continua, red optical slopes, and broad Balmer emission lines that challenge standard models and suggest a population of early, unusual active galactic nuclei (AGNs). Using a comprehensive photometric selection and public NIRSpec/PRISM spectroscopy across six JWST deep fields, we identify a large sample of 118 LRDs with high-S/N spectra, enabling a population-wide analysis of their UV-optical continuum and emission lines. We find clear correlations between rest-frame color ([0.3-0.9\,$μ$m]) and slopes: bluer LRDs have blue UV slopes ($β_{ν,\mathrm{UV}} \sim 0.3$) and red optical slopes, while redder LRDs exhibit redder UV slopes ($β_{ν,\mathrm{UV}} \sim 1.1$). The continuum shape shows a similar trend: redder LRDs display prominent Balmer breaks and curvature, while bluer LRDs follow power-law-like optical SEDs. From literature compilations, $\sim$60% of known broad-line AGNs satisfy our LRD criteria, and up to 90% of LRDs show broad Balmer lines. Emission-line diagnostics reveal a shift from high H$_α$/H$_β$ and low [OIII]$\lambda5007$/H$_β$ in redder LRDs to the opposite in bluer ones, along with stronger narrow-line equivalent widths, suggesting a transition from AGN- to host-dominated emission. We fit the spectra with a two-component model combining a gas-enshrouded black hole (BH) and a galaxy host. Redder LRDs require higher-luminosity, unreddened BHs and modestly reddened hosts; bluer LRDs require lower-luminosity, reddened BHs and dust-free galaxies. This framework reproduces the diversity in colors and spectral shape by varying BH luminosity, obscuration, and host-to-BH luminosity ratio.

preprint2026arXiv

The X-Ray Dot: Exotic Dust or a Late-Stage Little Red Dot?

JWST's "Little Red Dots" (LRDs) are increasingly interpreted as active galactic nuclei (AGN) obscured by dense thermalized gas rather than dust as evidenced by their X-ray weakness, blackbody-like continua, and Balmer line profiles. A key question is how LRDs connect to standard UV-luminous AGN and whether transitional phases exist and if they are observable. We present the "X-Ray Dot" (XRD), a compact source at $z=3.28$ observed by the NIRSpec WIDE GTO survey. The XRD exhibits LRD hallmarks: a blackbody-like ($T_{\rm eff} \simeq 6400\,$K) red continuum, a faint but blue rest-UV excess, falling mid-IR emission, and broad Balmer lines ($\rm FWHM \sim 2700-3200\,km\,s^{-1}$). Unlike LRDs, however, it is remarkably X-ray luminous ($L_\textrm{2$-$10$\,$keV} = 10^{44.18}\,$erg$\,$s$^{-1}$) and has a continuum inflection that is bluewards of the Balmer limit. We find that the red rest-optical and blue mid-IR continuum cannot be reproduced by standard dust-attenuated AGN models without invoking extremely steep extinction curves, nor can the weak mid-IR emission be reconciled with well-established X-ray--torus scaling relations. We therefore consider an alternative scenario: the XRD may be an LRD in transition, where the gas envelope dominates the optical continuum but optically thin sightlines allow X-rays to escape. The XRD may thus provide a physical link between LRDs and standard AGN, offering direct evidence that LRDs are powered by supermassive black holes and providing insight into their accretion properties.

preprint2025arXiv

The Noctua Suite of Simulations -- The Difficulty of Growing Massive Black Holes in Low-Mass Dwarf Galaxies

Aims. We study the individual and cumulative impact of stellar feedback processes on massive black hole (MBH) growth in a simulated low-mass dwarf galaxy. Methods. A suite of high-resolution radiation-hydrodynamic simulations called Noctua is performed, using the ArepoNoctua numerical framework for BHs in galaxy simulations. The chemical evolution of the gas is explicitly modelled in a time-dependent non-equilibrium way. Two types of stellar feedback are considered: individually-traced type II supernova (SNII) explosions, and radiatively transferred (on-the-fly) ionising stellar radiation (ISR) from OB stars. As part of the numerical framework, we develop and apply a novel physically-motivated model for MBH gas accretion, taking into account the angular momentum of the gas in the radiatively efficient regime, to estimate the gas accretion rate from the sub-grid accretion disc. Results. Without any stellar feedback, an initial $10^4~\mathrm{M}_\odot$ MBH is able to steadily grow over time, roughly doubling its mass after 800 Myr. Surprisingly, the growth of the MBH is more than doubled when only ISR feedback is considered, compared to the no stellar feedback run. This is due to the star formation rate (SFR) being highly suppressed (to a similar level or slightly above that when SNII feedback is considered), enabling a higher cumulative net gas inflow onto the MBH from not only the cold neutral- and molecular medium phases, but also the unstable- and warm neutral medium phases. With SNII feedback included, the gas accretion onto the MBH is episodic over time, and is suppressed by more than an order of magnitude already during the first 150 Myr. When combining SNII with ISR feedback, the growth of the MBH remains suppressed due to SNII feedback, but to a lesser extent compared to the SNII-only feedback run, due to a slightly lower SFR, and hence a reduced number of SNII events.

preprint2022arXiv

The BPT Diagram in Cosmological Galaxy Formation Simulations: Understanding the Physics Driving Offsets at High-Redshift

The Baldwin, Philips, & Terlevich diagram of [O III]/H$β$ vs. [N II]/H$α$ (hereafter N2-BPT) has long been used as a tool for classifying galaxies based on the dominant source of ionizing radiation. Recent observations have demonstrated that galaxies at $z\sim2$ reside offset from local galaxies in the N2-BPT space. In this paper, we conduct a series of controlled numerical experiments to understand the potential physical processes driving this offset. We model nebular line emission in a large sample of galaxies, taken from the SIMBA cosmological hydrodynamic galaxy formation simulation, using the CLOUDY photoionization code to compute the nebular line luminosities from H II regions. We find that the observed shift toward higher [O III]/H$β$ and [N II]/H$α$ values at high redshift arises from sample selection: when we consider only the most massive galaxies $M_* \sim 10^{10-11} M_\odot$, the offset naturally appears, due to their high metallicities. We predict that deeper observations that probe lower-mass galaxies will reveal galaxies that lie on a locus comparable to $z\sim 0$ observations. Even when accounting for sample selection effects, we find that there is a subtle mismatch between simulations and observations. To resolve this discrepancy, we investigate the impact of varying ionization parameters, H II region densities, gas-phase abundance patterns, and increasing radiation field hardness on N2-BPT diagrams. We find that either decreasing the ionization parameter or increasing the N/O ratio of galaxies at fixed O/H can move galaxies along a self-similar arc in N2-BPT space that is occupied by high-redshift galaxies.

preprint2022arXiv

The Metallicity Distribution Function in Outer Halo Fields of Simulated Elliptical Galaxies Compared to Observations of NGC 5128

Stellar metallicity distribution functions (MDF) have been measured for resolved stellar populations in the outer halos of many galaxies in nearby groups. Among them, the MDF of NGC 5128, the central giant elliptical in the Centaurus group, provides essential constraints for theories of massive galaxy formation and hierarchical assembly. To investigate the formation and chemical evolution history of the outer halo of giant elliptical galaxies, we examine the chemical properties of three zoom-in high resolution cosmological hydrodynamical simulations of an NGC 5128-like giant elliptical galaxy and compare their outer halo MDFs to the observed one of NGC 5128. Even though the simulated galaxies have different merging histories and age distributions, all predicted MDFs are in good qualitative agreement with the observed one. The median metallicity of the simulated galaxies is on average $\rm [M/H]=-0.41 \pm 0.06$ compared to the observed value of $\rm [M/H]=-0.38 \pm 0.02$ for NGC 5128, and the dispersion in metallicity is $\sim 0.77$ dex for both observed and simulated galaxies. We investigate the origin of the stars ending up in the outer halo field of simulated galaxies and show that most have an `accreted' origin, formed in other small galaxies and later accreted in mergers. Only $\sim 15$ percent of the stars are formed `in situ' within the main progenitor of galaxy and radially migrate outwards. We show that the contribution of metal-rich in situ stars is sub-dominant in the outer halos of our simulated galaxies, but can be prominent in the inner regions.

preprint2020arXiv

Finding Direct-Collapse Black Holes at Birth

Direct-collapse black holes (DCBHs) are currently one of the leading contenders for the origins of the first quasars in the universe, over 300 of which have now been found at $z >$ 6. But the birth of a DCBH in an atomically-cooling halo does not by itself guarantee it will become a quasar by $z \sim$ 7, the halo must also be located in cold accretion flows or later merge with a series of other gas-rich halos capable of fueling the BH's rapid growth. Here, we present near infrared luminosities for DCBHs born in cold accretion flows in which they are destined to grow to 10$^9$ M$_{\odot}$ by $z \sim$ 7. Our observables, which are derived from cosmological simulations with radiation hydrodynamics with Enzo, reveal that DCBHs could be found by the James Webb Space Telescope at $z \lesssim$ 20 and strongly-lensed DCBHs might be found in future wide-field surveys by Euclid and the Wide-Field Infrared Space Telescope at $z \lesssim$ 15.

preprint2020arXiv

Gas accretion regulates the scatter of the mass-metallicity relation

In this paper, we take advantage of the GAlaxy Evolution and Assembly (GAEA) semi-analytic model to analyse the origin of secondary dependencies in the local galaxy mass - gas metallicity relation. Our model reproduces quite well the trends observed in the local Universe as a function of galaxy star formation rate and different gas-mass phases. We show that the cold gas content (whose largest fraction is represented by the atomic gas phase) can be considered as the third parameter governing the scatter of the predicted mass-metallicity relation, in agreement with the most recent observational measurements. The trends can be explained with fluctuations of the gas accretion rates: a decrease of the gas supply leads to an increase of the gas metallicity due to star formation, while an increase of the available cold gas leads to a metallicity depletion. We demonstrate that the former process is responsible for offsets above the mass-metallicity relation, while the latter is responsible for deviations below the mass-metallicity relation. In low and intermediate mass galaxies, these negative offsets are primarily determined by late gas cooling dominated by material that has been previously ejected due to stellar feedback.

preprint2020arXiv

Reionization with galaxies and active galactic nuclei

In this work we investigate the properties of the sources that reionized the intergalactic medium (IGM) in the high-redshift Universe. Using a semi-analytical model aimed at reproducing galaxies and black holes in the first 1.5 Gyr of the Universe, we revisit the relative role of star formation and black hole accretion in producing ionizing photons that can escape into the IGM. Both star formation and black hole accretion are regulated by supernova feedback, resulting in black hole accretion being stunted in low-mass halos. We explore a wide range of combinations for the escape fraction of ionizing photons (redshift-dependent, constant and scaling with stellar mass) from both star formation ($\langle f_{\rm esc}^{\rm sf} \rangle$) and AGN ($f_{\rm esc}^{\rm bh}$) to find: (i) the ionizing budget is dominated by stellar radiation from low stellar mass ($M_*<10^9 {\rm M_\odot}$ ) galaxies at $z>6$ with the AGN contribution (driven by $M_{bh}>10^6 {\rm M_\odot}$ black holes in $M_* > 10^9 {\rm M_\odot}$ galaxies) dominating at lower redshifts; (ii) AGN only contribute $10-25\%$ to the cumulative ionizing emissivity by $z=4$ for the models that match the observed reionization constraints; (iii) if the stellar mass dependence of $\langle f_{\rm esc}^{\rm sf} \rangle$ is shallower than $f_{\rm esc}^{\rm bh}$, at $z<7$ a transition stellar mass exists above which AGN dominate the escaping ionizing photon production rate; (iv) the transition stellar mass decreases with decreasing redshift. While AGN dominate the escaping emissivity above the knee of the stellar mass function at $z \sim 6.8$, they take-over at stellar masses that are a tenth of the knee mass by $z=4$.

preprint2020arXiv

Stellar age gradients and inside-out star formation quenching in galaxy bulges

Radial age gradients hold the cumulative record of the multitude of physical processes driving the build-up of stellar populations and the ensuing star formation (SF) quenching process in galaxy bulges, therefore potentially sensitive discriminators between competing theoretical concepts on bulge formation and evolution. Based on spectral modeling of integral field spectroscopy data from the CALIFA survey, we derive mass- and light-weighted stellar age gradients ($\nabla$(t,B)L,M) within the photometrically determined bulge radius (RB) of a representative sample of local face-on late-type galaxies that span 2.6 dex in stellar mass. Our analysis documents a trend for decreasing $\nabla$(t,B)L,M with increasing M,T, with high-mass bulges predominantly showing negative age gradients and vice versa. The inversion from positive to negative $\nabla$(t,B)L,M occurs at logM,T ~ 10, which roughly coincides with the transition from lower-mass bulges whose gas excitation is powered by SF to bulges classified as Composite, LINER or Seyfert. We discuss two limiting cases for the origin of radial age gradients in massive LTG bulges. The first assumes that the stellar age in the bulge is initially spatially uniform, thus the observed age gradients arise from an inside-out SF quenching (ioSFQ) front that is radially expanding with a mean velocity vq. In this case, the age gradients translate into a slow ioSFQ that lasts until z~2, suggesting mild negative feedback by SF or an AGN. If negative age gradients in massive bulges are not due to ioSFQ but primarily due to their inside-out formation process, then the standard hypothesis of quasi-monolithic bulge formation has to be discarded in favor of a scenario that involves gradual buildup of stellar mass over 2-3 Gyr through, e.g., inside-out SF and inward migration of SF clumps from the disk. In this case, rapid AGN-driven ioSFQ cannot be ruled out.

preprint2020arXiv

The Fundamental Plane of Massive Quiescent Galaxies at z~2

We examine the Fundamental Plane (FP) and mass-to-light ratio ($M/L$) scaling relations using the largest sample of massive quiescent galaxies at $1.5<z<2.5$ to date. The FP ($r_{e}, σ_{e}, I_{e}$) is established using $19$ $UVJ$ quiescent galaxies from COSMOS with $Hubble$ $Space$ $Telescope$ $(HST)$ $H_{F160W}$ rest-frame optical sizes and X-shooter absorption line measured stellar velocity dispersions. For a very massive, ${\rm{log}}(M_{\ast}/M_{\odot})>11.26$, subset of 8 quiescent galaxies at $z>2$, from Stockmann et al. (2020), we show that they cannot passively evolve to the local Coma cluster relation alone and must undergo significant structural evolution to mimic the sizes of local massive galaxies. The evolution of the FP and $M/L$ scaling relations, from $z=2$ to present-day, for this subset are consistent with passive aging of the stellar population and minor merger structural evolution into the most massive galaxies in the Coma cluster and other massive elliptical galaxies from the MASSIVE Survey. Modeling the luminosity evolution from minor merger added stellar populations favors a history of merging with "dry" quiescent galaxies.

preprint2020arXiv

The influence of environment on satellite galaxies in the GAEA semi-analytic model

Reproducing the observed quenched fraction of satellite galaxies has been a long-standing issue for galaxy formation models. We modify the treatment of environmental effects in our state-of-the-art GAlaxy Evolution and Assembly (GAEA) semi-analytic model to improve our modelling of satellite galaxies. Specifically, we implement gradual stripping of hot gas, ram-pressure stripping of cold gas, and an updated algorithm to account for angular momentum exchanges between the gaseous and stellar disc components of model galaxies. Our updated model is in good agreement with the observed quenched fractions for central and satellite galaxies, and their dependencies on stellar mass and halo mass. We also find consistency between model predictions and observational estimates of quenching time for satellite galaxies; HI, H$_2$ fractions of central galaxies; and deficiencies of HI, H$_2$, SFR of galaxies in cluster haloes. Analysing results from our updated model, we find that the dominant quenching mechanisms are hot-gas stripping for low-mass galaxies, and AGN feedback for massive galaxies. The Ram-pressure stripping of cold gas has a limited impact on the predicted quenched fractions, nonetheless, is needed to reproduce the observed HI deficiencies.

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

Ionized gas kinematics of massive elliptical galaxies in CALIFA and in cosmological zoom-in simulations

(Abridged) We present an investigation of kinematical imprints of AGN feedback on the Warm Ionized gas Medium (WIM) of massive early-type galaxies (ETGs). To this end, we take a two-fold approach that involves a comparative analysis of Halpha velocity fields in 123 local ETGs from the CALIFA integral field spectroscopy survey with 20 simulated galaxies from high-resolution hydrodynamic cosmological SPHgal simulations. The latter were re-simulated for two modeling setups, one with and another without AGN feedback. In order to quantify the effects of AGN feedback on gas kinematics we measure three parameters that probe deviations from simple regular rotation using the kinemetry package. These indicators trace the possible presence of distinct kinematic components in Fourier space (k3,5/k1), variations in the radial profile of the kinematic major axis (sigma_PA), and offsets between the stellar and gas velocity fields (Delta Phi). These quantities are monitored in the simulations from a redshift 3 to 0.2 to assess the connection between black hole accretion history, stellar mass growth and kinematical perturbation of the WIM. Observed local massive galaxies show a broad range of irregularities, indicating disturbed warm gas motions, irrespective of being classified via diagnostic lines as AGN or not. Simulations of massive galaxies with AGN feedback generally exhibit higher irregularity parameters than without AGN feedback, more consistent with observations. Besides AGN feedback, other processes like major merger events or infalling gas clouds can lead to elevated irregularity parameters, but they are typically of shorter duration. More specifically, k3,5/k1 is most sensitive to AGN feedback, whereas Delta Phi is most strongly affected by gas infall.

preprint2019arXiv

X-shooter Spectroscopy and HST Imaging of 15 Ultra Massive Quiescent Galaxies at $z\gtrsim2$

We present a detailed analysis of a large sample of spectroscopically confirmed ultra-massive quiescent galaxies (${\rm{log}}(M_{\ast}/M_{\odot})\sim11.5$) at $z\gtrsim2$. This sample comprises 15 galaxies selected in the COSMOS and UDS fields by their bright K-band magnitudes and followed up with VLT/X-shooter spectroscopy and HST/WFC3 $H_{F160W}$ imaging. These observations allow us to unambiguously confirm their redshifts ascertain their quiescent nature and stellar ages, and to reliably assess their internal kinematics and effective radii. We find that these galaxies are compact, consistent with the high mass end of the mass-size relation for quiescent galaxies at $z=2$. Moreover, the distribution of the measured stellar velocity dispersions of the sample is consistent with the most massive local early-type galaxies from the MASSIVE Survey showing that evolution in these galaxies, is dominated by changes in size. The HST images reveal, as surprisingly high, that $40\ \%$ of the sample have tidal features suggestive of mergers and companions in close proximity, including three galaxies experiencing ongoing major mergers. The absence of velocity dispersion evolution from $z=2$ to $0$, coupled with a doubling of the stellar mass, with a factor of four size increase and the observed disturbed stellar morphologies support dry minor mergers as the primary drivers of the evolution of the massive quiescent galaxies over the last 10 billion years.

preprint2016arXiv

Galaxy assembly, stellar feedback and metal enrichment: the view from the GAEA model

One major problem of current theoretical models of galaxy formation is given by their inability to reproduce the apparently `anti-hierarchical' evolution of galaxy assembly: massive galaxies appear to be in place since $z\sim 3$, while a significant increase of the number densities of low mass galaxies is measured with decreasing redshift. In this work, we perform a systematic analysis of the influence of different stellar feedback schemes, carried out in the framework of GAEA, a new semi-analytic model of galaxy formation. It includes a self-consistent treatment for the timings of gas, metal and energy recycling, and for the chemical yields. We show this to be crucial to use observational measurements of the metallicity as independent and powerful constraints for the adopted feedback schemes. The observed trends can be reproduced in the framework of either a strong ejective or preventive feedback model. In the former case, the gas ejection rate must decrease significantly with cosmic time (as suggested by parametrizations of the cosmological `FIRE' simulations). Irrespective of the feedback scheme used, our successful models always imply that up to 60-70 per cent of the baryons reside in an `ejected' reservoir and are unavailable for cooling at high redshift. The same schemes predict physical properties of model galaxies (e.g. gas content, colour, age, and metallicity) that are in much better agreement with observational data than our fiducial model. The overall fraction of passive galaxies is found to be primarily determined by internal physical processes, with environment playing a secondary role.

preprint2016arXiv

HI-selected Galaxies in Hierarchical Models of Galaxy Formation and Evolution

In this work, we study the basic statistical properties of HI-selected galaxies extracted from six different semi-analytic models, all run on the same cosmological N-body simulation. One model includes an explicit treatment for the partition of cold gas into atomic and molecular hydrogen. All models considered agree nicely with the measured HI mass function in the local Universe, with the measured scaling relations between HI and galaxy stellar mass, and with the predicted 2-point correlation function for HI rich galaxies. One exception is given by one model that predicts very little HI associated with galaxies in haloes above 10^12 Msun: we argue this is due to a too efficient radio-mode feedback for central galaxies, and to a combination of efficient stellar feedback and instantaneous stripping of hot gas for satellites. We demonstrate that treatment of satellite galaxies introduces large uncertainties at low HI masses. While models assuming non instantaneous stripping of hot gas tend to form satellite galaxies with HI masses slightly smaller than those of centrals with the same stellar mass, instantaneous gas stripping does not translate necessarily in lower HI masses. In fact, the adopted stellar feedback and star formation affect the satellites too. We analyze the relation between HI content and spin of simulated haloes: low spin haloes tend to host HI poor galaxies while high spin haloes are populated by galaxies in a wide range of HI mass. In our simulations, this is due to a correlation between the initial gas disk size and the halo spin.

preprint2016arXiv

Origin and properties of dual and offset active galactic nuclei in a cosmological simulation at z=2

In the last few years, it became possible to observationally resolve galaxies with two distinct nuclei in their centre. For separations smaller than 10kpc, dual and offset active galactic nuclei (AGN) are distinguished: in dual AGN, both nuclei are active, whereas in offset AGN only one nucleus is active. To study the origin of such AGN pairs, we employ a cosmological, hydrodynamic simulation with a large volume of (182Mpc)^3 from the set of Magneticum Pathfinder Simulations. The simulation self-consistently produces 35 resolved black hole (BH) pairs at redshift z=2, with a comoving distance smaller than 10kpc. 14 of them are offset AGN and nine are dual AGN, resulting in a fraction of (1.2 \pm 0.3)% AGN pairs with respect to the total number of AGN. In this paper, we discuss fundamental differences between the BH and galaxy properties of dual AGN, offset AGN and inactive BH pairs and investigate their different triggering mechanisms. We find that in dual AGN the BHs have similar masses and the corresponding BH from the less massive progenitor galaxy always accretes with a higher Eddington ratio. In contrast, in offset AGN the active BH is typically more massive than its non-active counterpart. Furthermore, dual AGN in general accrete more gas from the intergalactic medium than offset AGN and non-active BH pairs. This highlights that merger events, particularly minor mergers, do not necessarily lead to strong gas inflows and thus, do not always drive strong nuclear activity.

preprint2016arXiv

Variations of the initial mass function in semi-analytical models: implications for the mass assembly and the chemical enrichment of galaxies in the GAEA model

In this work, we investigate the implications of the Integrated Galaxy-wide stellar Initial Mass Function (IGIMF) approach in the framework of the semi-analytic model GAEA (GAlaxy Evolution and Assembly), which features a detailed treatment of chemical enrichment and stellar feedback. The IGIMF provides an analytic description of the dependence of the stellar IMF shape on the rate of star formation in galaxies. We find that our model with a universal IMF predicts a rather flat [$α$/Fe]-stellar mass relation. The model assuming the IGIMF, instead, is able to reproduce the observed increase of $α$-enhancement with stellar mass, in agreement with previous studies. This is mainly due to the fact that massive galaxies are characterized by larger star formation rates at high-redshift, leading to stronger $α$-enhancement with respect to low-mass galaxies. At the same time, the IGIMF hypothesis does not affect significantly the trend for shorter star formation timescales for more massive galaxies. We argue that in the IGIMF scenario the [$α$/Fe] ratios are good tracers of the highest star formation events. The final stellar masses and mass-to-light-ratio of our model massive galaxies are larger than those estimated from the synthetic photometry assuming a universal IMF, providing a self-consistent interpretation of similar recent results, based on dynamical analysis of local early type galaxies.

preprint2015arXiv

A refined sub-grid model for black hole accretion and AGN feedback in large cosmological simulations

In large scale cosmological hydrodynamic simulations simplified sub-grid models for gas accretion onto black holes and AGN feedback are commonly used. Such models typically depend on various free parameters, which are not well constrained. We present a new advanced model containing a more detailed description of AGN feedback, where those parameters reflect the results of recent observations. The model takes the dependency of these parameters on the black hole properties into account and describes a continuous transition between the feedback processes acting in the so-called radio-mode and quasar-mode. In addition, we implement a more detailed description of the accretion of gas onto black holes by distinguishing between hot and cold gas accretion. Our new implementations prevent black holes from gaining too much mass, particularly at low redshifts so that our simulations are now very successful in reproducing the observed present-day black hole mass function. Our new model also suppresses star formation in massive galaxies slightly more efficiently than many state-of-the-art models. Therefore, the simulations that include our new implementations produce a more realistic population of quiescent and star-forming galaxies compared to recent observations, even if some discrepancies remain. In addition, the baryon conversion efficiencies in our simulation are - except for the high mass end - consistent with observations presented in literature over the mass range resolved by our simulations. Finally, we discuss the significant impact of the feedback model on the low-luminous end of the AGN luminosity function.

preprint2015arXiv

nIFTy Cosmology: Comparison of Galaxy Formation Models

We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The participating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data sets, stellar models, and merger trees. In this paper we apply them without recalibration and this leads to a wide variety of predictions for the stellar mass function, specific star formation rates, stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter is much larger than seen in previous comparison studies primarily because the codes have been used outside of their native environment within which they are well tested and calibrated. The purpose of the `nIFTy comparison of galaxy formation models' is to bring together as many different galaxy formation modellers as possible and to investigate a common approach to model calibration. This paper provides a unified description for all participating models and presents the initial, uncalibrated comparison as a baseline for our future studies where we will develop a common calibration framework and address the extent to which that reduces the scatter in the model predictions seen here.

preprint2015arXiv

On the dependence of galaxy morphologies on galaxy mergers

The distribution of galaxy morphological types is a key test for models of galaxy formation and evolution, providing strong constraints on the relative contribution of different physical processes responsible for the growth of the spheroidal components. In this paper, we make use of a suite of semi-analytic models to study the efficiency of galaxy mergers in disrupting galaxy discs and building galaxy bulges. In particular, we compare standard prescriptions usually adopted in semi-analytic models, with new prescriptions proposed by Kannan et al., based on results from high-resolution hydrodynamical simulations, and we show that these new implementations reduce the efficiency of bulge formation through mergers. In addition, we compare our model results with a variety of observational measurements of the fraction of spheroid-dominated galaxies as a function of stellar and halo mass, showing that the present uncertainties in the data represent an important limitation to our understanding of spheroid formation. Our results indicate that the main tension between theoretical models and observations does not stem from the survival of purely disc structures (i.e. bulgeless galaxies), rather from the distribution of galaxies of different morphological types, as a function of their stellar mass.

preprint2015arXiv

The stellar accretion origin of stellar population gradients in massive galaxies at large radii

We investigate the evolution of stellar population gradients from $z=2$ to $z=0$ in massive galaxies at large radii ($r > 2R_{\mathrm{eff}}$) using ten cosmological zoom simulations of halos with $6 \times 10^{12} M_{\odot} < M_{\mathrm{halo}} < 2 \times 10^{13}M_{\odot}$. The simulations follow metal cooling and enrichment from SNII, SNIa and AGB winds. We explore the differential impact of an empirical model for galactic winds that reproduces the mass-metallicity relation and its evolution with redshift. At larger radii the galaxies, for both models, become more dominated by stars accreted from satellite galaxies in major and minor mergers. In the wind model, fewer stars are accreted, but they are significantly more metal poor resulting in steep global metallicity ($\langle \nabla Z_{\mathrm{stars}} \rangle= -0.35$ dex/dex) and color (e.g. $\langle \nabla g-r \rangle = -0.13$ dex/dex) gradients in agreement with observations. In contrast, colour and metallicity gradients of the models without winds are inconsistent with observations. Age gradients are in general mildly positive at $z=0$ ($\langle \nabla Age_{\mathrm{stars}} \rangle= 0.04$ dex/dex) with significant differences between the models at higher redshift. We demonstrate that for the wind model, stellar accretion is steepening existing in-situ metallicity gradients by about 0.2 dex by the present day and helps to match observed gradients of massive early-type galaxies at large radii. Colour and metallicity gradients are significantly steeper for systems which have accreted stars in minor mergers, while galaxies with major mergers have relatively flat gradients, confirming previous results. This study highlights the importance of stellar accretion for stellar population properties of massive galaxies at large radii, which can provide important constraints for formation models.

preprint2014arXiv

Cosmological simulations of black hole growth: AGN luminosities and downsizing

In this study, we present a detailed, statistical analysis of black hole growth and the evolution of active galactic nuclei (AGN) using cosmological hydrodynamic simulations run down to $z=0$. The simulations self-consistently follow radiative cooling, star formation, metal enrichment, black hole growth and associated feedback processes from both supernovae typeII/Ia and AGN. We consider two simulation runs, one with a large co-moving volume of $(500\ \mathrm{Mpc})^3$ and one with a smaller volume of $(68\ \mathrm{Mpc})^3$ but with a by a factor of almost 20 higher mass resolution. Consistently with previous results, our simulations can widely match observed black hole properties of the local Universe. Furthermore, our simulations can successfully reproduce the evolution of the bolometric AGN luminosity function for both the low-luminosity and the high-luminosity end up to $z=3.0$. In addition, the smaller but higher resolution run is able to match the observational data of the low bolometric luminosity end at higher redshifts $z=3-4$. We also perform a direct comparison with the observed soft and hard X-ray luminosity functions of AGN, including an empirical correction for a torus-level obscuration, and find a similarly good agreement. These results nicely demonstrate that the observed "anti-hierarchical" trend in the AGN number density evolution (i.e. the number densities of luminous AGN peak at higher redshifts than those of faint AGN) is self-consistently predicted by our simulations. Implications of this downsizing behaviour on active black holes, their masses and Eddington-ratios are discussed. Overall, the downsizing behaviour in the AGN number density as a function of redshift can be mainly attributed to the evolution of the gas density in the resolved vicinity of a (massive) black hole. (shortened)

preprint2014arXiv

The Atlas3D project - XXVI. HI discs in real and simulated fast and slow rotators

One quarter of all nearby early-type galaxies (ETGs) outside Virgo host a disc/ring of HI with size from a few to tens of kpc and mass up to ~1e+9 solar masses. Here we investigate whether this HI is related to the presence of a stellar disc within the host making use of the classification of ETGs in fast and slow rotators (FR/SR). We find a large diversity of HI masses and morphologies within both families. Surprisingly, SRs are detected as often, host as much HI and have a similar rate of HI discs/rings as FRs. Accretion of HI is therefore not always linked to the growth of an inner stellar disc. The weak relation between HI and stellar disc is confirmed by their frequent kinematical misalignment in FRs, including cases of polar and counterrotating gas. In SRs the HI is usually polar. This complex picture highlights a diversity of ETG formation histories which may be lost in the relative simplicity of their inner structure and emerges when studying their outer regions. We find that LCDM hydrodynamical simulations have difficulties reproducing the HI properties of ETGs. The gas discs formed in simulations are either too massive or too small depending on the star formation feedback implementation. Kinematical misalignments match the observations only qualitatively. The main point of conflict is that nearly all simulated FRs and a large fraction of all simulated SRs host corotating HI. This establishes the HI properties of ETGs as a novel challenge to simulations.

preprint2014arXiv

The definition of environment and its relation to the quenching of galaxies at z=1-2 in a hierarchical Universe

A well calibrated method to describe the environment of galaxies at all redshifts is essential for the study of structure formation. Such a calibration should include well understood correlations with halo mass, and the possibility to identify galaxies which dominate their potential well (centrals), and their satellites. Focusing on z = 1 and 2 we propose a method of environmental calibration which can be applied to the next generation of low to medium resolution spectroscopic surveys. Using an up-to-date semi-analytic model of galaxy formation, we measure the local density of galaxies in fixed apertures on different scales. There is a clear correlation of density with halo mass for satellite galaxies, while a significant population of low mass centrals is found at high densities in the neighbourhood of massive haloes. In this case the density simply traces the mass of the most massive halo within the aperture. To identify central and satellite galaxies, we apply an observationally motivated stellar mass rank method which is both highly pure and complete, especially in the more massive haloes where such a division is most meaningful. Finally we examine a test case for the recovery of environmental trends: the passive fraction of galaxies and its dependence on stellar and halo mass for centrals and satellites. With careful calibration, observationally defined quantities do a good job of recovering known trends in the model. This result stands even with reduced redshift accuracy, provided the sample is deep enough to preserve a wide dynamic range of density.

preprint2014arXiv

The influence of the environmental history on quenching star formation in a $Λ$CDM universe

We present a detailed analysis of the influence of the environment and of the environmental history on quenching star formation in central and satellite galaxies in the local Universe. We take advantage of publicly available galaxy catalogues obtained from applying a galaxy formation model to the Millennium simulation. In addition to halo mass, we consider the local density of galaxies within various fixed scales. Comparing our model predictions to observational data (SDSS), we demonstrate that the models are failing to reproduce the observed density dependence of the quiescent galaxy fraction in several aspects: for most of the stellar mass ranges and densities explored, models cannot reproduce the observed similar behaviour of centrals and satellites, they slightly under-estimate the quiescent fraction of centrals and significantly over-estimate that of satellites. We show that in the models, the density dependence of the quiescent central galaxies is caused by a fraction of "backsplash" centrals which have been satellites in the past (and were thus suffering from environmental processes). Turning to satellite galaxies, the density dependence of their quiescent fractions reflects a dependence on the time spent orbiting within a parent halo of a particular mass, correlating strongly with halo mass and distance from the halo centre. Comparisons with observational estimates suggest relatively long gas consumption time scales of roughly 5 Gyr in low mass satellite galaxies. The quenching time scales decrease with increasing satellite stellar mass. Overall, a change in modelling both internal processes (star formation and feedback) and environmental processes (e.g. making them dependent on dynamical friction time-scales and preventing the re-accretion of gas onto backsplash galaxies) is required for improving currently used galaxy formation models.

preprint2013arXiv

Isolated galaxies in hierarchical galaxy formation models - present-day properties and environmental histories

In this study, we have carried out a detailed, statistical analysis of isolated model galaxies, taking advantage of publicly available hierarchical galaxy formation models. To select isolated galaxies, we employ 2D methods widely used in the observational literature, as well as a more stringent 3D isolation criterion that uses the full 3D-real space information. In qualitative agreement with observational results, isolated model galaxies have larger fractions of late-type, star forming galaxies with respect to randomly selected samples of galaxies with the same mass distribution. We also find that the samples of isolated model galaxies typically contain a fraction of less than 15 per cent of satellite galaxies, that reside at the outskirts of their parent haloes where the galaxy number density is low. Projection effects cause a contamination of 2D samples of about 18 per cent, while we estimate a typical completeness of 65 per cent. Our model isolated samples also include a very small (few per cent) fraction of bulge dominated galaxies (B/T > 0.8) whose bulges have been built mainly by minor mergers. Our study demonstrates that about 65-70 per cent of 2D isolated galaxies that are classified as isolated at z = 0 have indeed been completely isolated since z = 1 and only 7 per cent have had more than 3 neighbours within a comoving radius of 1 Mpc. Irrespectively of the isolation criteria, roughly 45 per cent of isolated galaxies have experienced at least one merger event in the past (most of the mergers are minor, with mass ratios between 1:4 and 1:10). The latter point validates the approximation that isolated galaxies have been mainly influenced by internal processes.

preprint2013arXiv

The effect of metal enrichment and galactic winds on galaxy formation in cosmological zoom simulations

We investigate the differential effects of metal cooling and galactic stellar winds on the cosmological formation of individual galaxies with three sets of cosmological, hydrodynamical zoom simulations of 45 halos in the mass range 10^11<M_halo<10^13M_sun. Models including both galactic winds and metal cooling (i) suppress early star formation at z>1 and predict reasonable star formation histories, (ii) produce galaxies with high cold gas fractions (30-60 per cent) at high redshift, (iii) significantly reduce the galaxy formation efficiencies for halos (M_halo<10^12M_sun) at all redshifts in agreement with observational and abundance matching constraints, (iv) result in high-redshift galaxies with reduced circular velocities matching the observed Tully-Fisher relation at z~2, and (v) significantly increase the sizes of low-mass galaxies (M_stellar<3x10^10M_sun) at high redshift resulting in a weak size evolution - a trend in agreement with observations. However, the low redshift (z<0.5) star formation rates of massive galaxies are higher than observed (up to ten times). No tested model predicts the observed size evolution for low-mass and high-mass galaxies simultaneously. Due to the delayed onset of star formation in the wind models, the metal enrichment of gas and stars is delayed and agrees well with observational constraints. Metal cooling and stellar winds are both found to increase the ratio of in situ formed to accreted stars - the relative importance of dissipative vs. dissipationless assembly. For halo masses below ~10^12M_sun, this is mainly caused by less stellar accretion and compares well to predictions from semi-analytical models but still differs from abundance matching models. For higher masses, the fraction of in situ stars is over-predicted due to the unrealistically high star formation rates at low redshifts.

preprint2012arXiv

Origin of the anti-hierarchical growth of black holes

Observational studies have revealed a "downsizing" trend in black hole (BH) growth: the number densities of luminous AGN peak at higher redshifts than those of faint AGN. This would seem to imply that massive black holes formed before low mass black holes, in apparent contradiction to hierarchical clustering scenarios. We investigate whether this observed "downsizing" in BH growth is reproduced in a semi-analytic model for the formation and evolution of galaxies and black holes, set within the hierarchical paradigm for structure formation (Somerville et al. 2008; S08). In this model, black holes evolve from light seeds (\sim100M\odot) and their growth is merger-driven. The original S08 model (baseline model) reproduces the number density of AGN at intermediate redshifts and luminosities, but underproduces luminous AGN at very high redshift (z > 3) and overproduces them at low redshift (z < 1). In addition, the baseline model underproduces low-luminosity AGN at low redshift (z < 1). To solve these problems we consider several modifications to the physical processes in the model: (1) a 'heavy' black hole seeding scenario (2) a sub-Eddington accretion rate ceiling that depends on the cold gas fraction, and (3) an additional black hole accretion mode due to disk instabilities. With these three modifications, the models can explain the observed downsizing, successfully reproduce the bolometric AGN luminosity function and simultaneously reproduce galaxy and black hole properties in the local Universe. We also perform a comparison with the observed soft and hard X-ray luminosity functions of AGN, including an empirical correction for torus-level obscuration, and reach similar conclusions. Our best-fit model suggests a scenario in which disk instabilities are the main driver for moderately luminous Seyfert galaxies at low redshift, while major mergers are the main trigger for luminous AGN.