-
Galactic disc heating by density granulation in fuzzy dark matter simulations
Authors:
Hsun-Yeong Yang,
Barry T. Chiang,
Guan-Ming Su,
Hsi-Yu Schive,
Tzihong Chiueh,
Jeremiah P. Ostriker
Abstract:
Fuzzy dark matter (FDM), an attractive dark matter candidate comprising ultralight bosons (axions) with a particle mass $m_a\sim10^{-22}$ eV, is motivated by the small-scale challenges of cold dark matter and features a kpc-size de Broglie wavelength. Quantum wave interference inside an FDM halo gives rise to stochastically fluctuating density granulation; the resulting gravitational perturbations…
▽ More
Fuzzy dark matter (FDM), an attractive dark matter candidate comprising ultralight bosons (axions) with a particle mass $m_a\sim10^{-22}$ eV, is motivated by the small-scale challenges of cold dark matter and features a kpc-size de Broglie wavelength. Quantum wave interference inside an FDM halo gives rise to stochastically fluctuating density granulation; the resulting gravitational perturbations could drive significant disc thickening, providing a natural explanation for galactic thick discs. Here we present the first self-consistent simulations of FDM haloes and stellar discs, exploring $m_a=0.2-1.2\times10^{-22}$ eV and halo masses $M_\text{h} = 0.7-2.8\times10^{11}$ M$_\odot$. Disc thickening is observed in all simulated systems. The disc heating rates are approximately constant in time and increase substantially with decreasing $m_a$, reaching $dh/dt \simeq 0.04$ ($0.4$) kpc Gyr$^{-1}$ and $dσ_z^2/dt \simeq4$ ($150$) km$^2$s$^{-2}$Gyr$^{-1}$ for $m_a=1.2$ ($0.2$) $\times10^{-22}$ eV and $M_\text{h} =7\times10^{10} \text{M}_\odot$, where $h$ is the disc scale height and $σ_z$ is the vertical velocity dispersion. These simulated heating rates agree within a factor of two with the theoretical estimates of Chiang et al., confirming that the rough estimate of Church et al. overpredicts the granulation-driven disc heating rate by two orders of magnitude. However, the simulation-inferred heating rates scale less steeply than the theoretically predicted relation $dσ^2_z/dt \propto m_a^{-3}$. Finally, we examine the applicability of the Fokker-Planck approximation in FDM granulation modelling and the robustness of the $m_a$ exclusion bound derived from the Galactic disc kinematics.
△ Less
Submitted 14 March, 2024;
originally announced March 2024.
-
The Origins of Gas Accreted by Supermassive Black Holes: the Importance of Recycled Gas
Authors:
Ena Choi,
Rachel S. Somerville,
Jeremiah P. Ostriker,
Michaela Hirschmann,
Thorsten Naab
Abstract:
We investigate the fueling mechanisms of supermassive black holes (SMBHs) by analyzing ten zoom-in cosmological simulations of massive galaxies, with stellar masses $10^{11-12} M_{\odot}$ and SMBH masses $10^{8.9-9.7}$ at $z=0$ and featuring various major and minor merger events. By tracing the gas history in these simulations, we categorize the gas accreted by the central SMBHs based on its origi…
▽ More
We investigate the fueling mechanisms of supermassive black holes (SMBHs) by analyzing ten zoom-in cosmological simulations of massive galaxies, with stellar masses $10^{11-12} M_{\odot}$ and SMBH masses $10^{8.9-9.7}$ at $z=0$ and featuring various major and minor merger events. By tracing the gas history in these simulations, we categorize the gas accreted by the central SMBHs based on its origin. Gas that belonged to a different galaxy before accretion onto the BH is labeled as (1) ``external," while smoothly accreted cosmic gas is classified as (2) ``smooth." Gas produced within the primary halo through stellar evolution and subsequently accreted by the SMBH is classified as (3) ``recycled." Our analysis, which included stellar feedback, reveals that the primary fuel source for SMBHs is the recycled gas from dying stars. This recycled gas from stars in the inner region of the galaxy readily collapses toward the center, triggering starbursts, and simultaneously fueling the SMBH. Galaxy mergers also play a crucial role in fueling SMBHs in massive galaxies as SMBHs in massive halos tend to accrete a higher fraction of external gas from mergers compared to smoothly accreted gas. However, on average, it takes approximately 1.85 Gyr for external gas to enter the main galaxy and accrete onto the SMBH. Considering the presence of various other gas triggers for AGN activity alongside this time delay, the association between AGN and mergers may not always be obvious.
△ Less
Submitted 19 February, 2024; v1 submitted 13 December, 2023;
originally announced December 2023.
-
Active galactic nuclei feedback in an elliptical galaxy (III): the impacts and fate of cosmological inflow
Authors:
Bocheng Zhu,
Feng Yuan,
Suoqing Ji,
Yingjie Peng,
Luis C. Ho,
Jeremiah P. Ostriker,
Luca Ciotti
Abstract:
The cosmological inflow of a galaxy is speculated to be able to enter the galaxy and enhance the star formation rate (SFR) and black hole accretion rate (BHAR). In this paper, by performing high-resolution hydrodynamic simulations in the framework of {\it MACER}, we investigate the fate of the inflow and its impacts on the evolution of a massive elliptical galaxy. The inflow properties are adopted…
▽ More
The cosmological inflow of a galaxy is speculated to be able to enter the galaxy and enhance the star formation rate (SFR) and black hole accretion rate (BHAR). In this paper, by performing high-resolution hydrodynamic simulations in the framework of {\it MACER}, we investigate the fate of the inflow and its impacts on the evolution of a massive elliptical galaxy. The inflow properties are adopted from the cosmological simulation IllustrisTNG. We find that, the inflow gas hardly enters but is blocked beyond $\sim 20$ kpc from the central galaxy and becomes part of the circumgalactic medium (CGM). The gas pressure gradient, mainly contributed by the thermalized stellar wind and subdominantly by the energy input from the AGN, balances gravity and prevents the inflow from entering the galaxy. The SFR and BHAR are almost not affected by the normal inflow. However, if the rate of cosmological inflow were increased by a factor of 3, a small fraction of the inflow would enter the galaxy and contribute about 10\% of the gas in the galaxy. In this case, the gas density in the galaxy would increase by a factor of $\ga$ 20. This increase is not because of the additional gas supply by the inflow but the increase of gas density in the CGM caused by the inflow. Consequently, the SFR and BHAR would increase by a factor of $\sim$ 5 and $\sim 1000$ respectively. Finally, AGN feedback can perturb the motion of the inflow and heat the CGM through its intermittent outbursts.
△ Less
Submitted 5 July, 2023; v1 submitted 7 March, 2023;
originally announced March 2023.
-
The growth of intermediate mass black holes through tidal captures and tidal disruption events
Authors:
Francesco Paolo Rizzuto,
Thorsten Naab,
Antti Rantala,
Peter H. Johansson,
Jeremiah P. Ostriker,
Nicholas C. Stone,
Shihong Liao,
Dimitrios Irodotou
Abstract:
We present $N\mathrm{-body} $ simulations, including post-Newtonian dynamics, of dense clusters of low-mass stars harbouring central black holes (BHs) with initial masses of 50, 300, and 2000 $\mathrm{M_{\odot}}$. The models are evolved with the $N\mathrm{-body} $ code \textsc{bifrost} to investigate the possible formation and growth of massive BHs by the tidal capture of stars and tidal disruptio…
▽ More
We present $N\mathrm{-body} $ simulations, including post-Newtonian dynamics, of dense clusters of low-mass stars harbouring central black holes (BHs) with initial masses of 50, 300, and 2000 $\mathrm{M_{\odot}}$. The models are evolved with the $N\mathrm{-body} $ code \textsc{bifrost} to investigate the possible formation and growth of massive BHs by the tidal capture of stars and tidal disruption events (TDEs). We model star-BH tidal interactions using a velocity-dependent drag force, which causes orbital energy and angular momentum loss near the BH. About $\sim 20-30$ per cent of the stars within the spheres of influence of the black holes form Bahcall-Wolf cusps and prevent the systems from core collapse. Within the first 40 Myr of evolution, the systems experience 500 up to 1300 TDEs, depending on the initial cluster structure. Most ($> 95$ per cent) of the TDEs originate from stars in the Bahcall-Wolf cusp. We derive an analytical formula for the TDE rate as a function of the central BH mass, density and velocity dispersion of the clusters ($\dot{N}_{\mathrm{TDE}} \propto M\mathrm{_{BH}} ρσ^{-3}$). We find that TDEs can lead a 300 $\mathrm{M_{\odot}}$ BH to reach $\sim 7000 \mathrm{M_{\odot}}$ within a Gyr. This indicates that TDEs can drive the formation and growth of massive BHs in sufficiently dense environments, which might be present in the central regions of nuclear star clusters.
△ Less
Submitted 23 November, 2022;
originally announced November 2022.
-
Can ultralight dark matter explain the age-velocity dispersion relation of the Milky Way disc: A revised and improved treatment
Authors:
Barry T. Chiang,
Jeremiah P. Ostriker,
Hsi-Yu Schive
Abstract:
Ultralight axion-like particles $m_a \sim 10^{-22}$ eV, or Fuzzy Dark Matter (FDM), behave comparably to cold dark matter (CDM) on cosmological scales and exhibit a kpc-size de Broglie wavelength capable of alleviating established (sub-)galactic-scale problems of CDM. Substructures inside an FDM halo incur gravitational potential perturbations, resulting in stellar heating sufficient to account fo…
▽ More
Ultralight axion-like particles $m_a \sim 10^{-22}$ eV, or Fuzzy Dark Matter (FDM), behave comparably to cold dark matter (CDM) on cosmological scales and exhibit a kpc-size de Broglie wavelength capable of alleviating established (sub-)galactic-scale problems of CDM. Substructures inside an FDM halo incur gravitational potential perturbations, resulting in stellar heating sufficient to account for the Galactic disc thickening over a Hubble time, as first demonstrated by Church et al. We present a more sophisticated treatment that incorporates the full baryon and dark matter distributions of the Milky Way and adopts stellar disc kinematics inferred from recent Gaia, APOGEE, and LAMOST surveys. Ubiquitous density granulation and subhalo passages respectively drive inner disc thickening and flaring of the outer disc, resulting in an observationally consistent `U-shaped' disc vertical velocity dispersion profile with the global minimum located near the solar radius. The observed age-velocity dispersion relation in the solar vicinity can be explained by the FDM-substructure-induced heating and places an exclusion bound $m_a \gtrsim 0.4\times10^{-22}$ eV. We assess non-trivial uncertainties in the empirical core-halo relation, FDM subhalo mass function and tidal stripping, and stellar heating estimate. The mass range $m_a\simeq 0.5-0.7\times10^{-22}$ eV favoured by the observed thick disc kinematics is in tension with several exclusion bounds inferred from dwarf density profiles, stellar streams, and Milky Way satellite populations, which could be significantly relaxed due to the aforesaid uncertainties. Additionally, strongly anisotropic heating could help explain the formation of ultra-thin disc galaxies.
△ Less
Submitted 5 December, 2022; v1 submitted 14 November, 2022;
originally announced November 2022.
-
Loss Cone Shielding
Authors:
Odelia Teboul,
Nicholas C. Stone,
Jeremiah P. Ostriker
Abstract:
A star wandering close enough to a massive black hole (MBH) can be ripped apart by the tidal forces of the black hole. The advent of wide-field surveys at many wavelengths has quickly increased the number of tidal disruption events (TDEs) observed, and has revealed that i) observed TDE rates are lower than theoretical predictions and ii) E+A galaxies are significantly overrepresented. This overrep…
▽ More
A star wandering close enough to a massive black hole (MBH) can be ripped apart by the tidal forces of the black hole. The advent of wide-field surveys at many wavelengths has quickly increased the number of tidal disruption events (TDEs) observed, and has revealed that i) observed TDE rates are lower than theoretical predictions and ii) E+A galaxies are significantly overrepresented. This overrepresentation further worsens the tension between observed and theoretically predicted TDEs for non-E+A galaxies. Classical loss cone theory focuses on the cumulative effect of many weak scatterings. However, a strong scattering can remove a star from the distribution before it can get tidally disrupted. Most stars undergoing TDEs come from within the radius of influence, the densest environments of the universe. In such environments, close encounters rare elsewhere become non-negligible. We revise the standard loss cone theory to take into account classical two-body interactions as well as strong scattering, collisions, tidal captures, and study under which conditions close encounters can shield the loss cone. We i) analytically derive the impact of strong scattering and other close encounters, ii) compute time-dependent loss cone dynamics including both weak and strong encounters, and iii) derive analytical solutions to the Fokker-Planck equation with strong scattering. We find that i) TDE rates can be reduced to up to an order of magnitude and ii) strong shielding preferentially reduces deeply plunging stars. We also show that stellar overdensities, one possible explanation for the E+A preference, can fail to increase TDE rates when taking into account strong scattering.
△ Less
Submitted 26 October, 2023; v1 submitted 10 November, 2022;
originally announced November 2022.
-
The Metallicity Distribution Function in Outer Halo Fields of Simulated Elliptical Galaxies Compared to Observations of NGC 5128
Authors:
Ena Choi,
Jeremiah P. Ostriker,
Michaela Hirschmann,
Rachel S. Somerville,
Thorsten Naab
Abstract:
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…
▽ More
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.
△ Less
Submitted 14 March, 2022;
originally announced March 2022.
-
Massive black hole formation in dense stellar environments: Enhanced X-ray detection rates in high velocity dispersion nuclear star clusters
Authors:
Vivienne F. Baldassare,
Nicholas C. Stone,
Adi Foord,
Elena Gallo,
Jeremiah P. Ostriker
Abstract:
We analyze Chandra X-ray Observatory imaging of 108 galaxies hosting nuclear star clusters (NSCs) to search for signatures of massive black holes (BHs). NSCs are extremely dense stellar environments with conditions that can theoretically facilitate massive BH formation. Recent work by Stone et al. (2017) finds that sufficiently dense NSCs should be unstable to the runaway growth of a stellar mass…
▽ More
We analyze Chandra X-ray Observatory imaging of 108 galaxies hosting nuclear star clusters (NSCs) to search for signatures of massive black holes (BHs). NSCs are extremely dense stellar environments with conditions that can theoretically facilitate massive BH formation. Recent work by Stone et al. (2017) finds that sufficiently dense NSCs should be unstable to the runaway growth of a stellar mass BH into a massive BH via tidal captures. Furthermore, there is a velocity dispersion threshold ($40\;\rm{km\;s^{-1}}$) above which NSCs should inevitably form a massive BH. To provide an observational test of these theories, we measure X-ray emission from NSCs and compare to the measured velocity dispersion and tidal capture runaway timescale. We find that NSCs above the $40\;\rm{km\;s^{-1}}$ threshold are X-ray detected at roughly twice the rate of those below (after accounting for contamination from X-ray binaries). These results are consistent with a scenario in which dense, high-velocity NSCs can form massive BHs, providing a formation pathway that does not rely on conditions found only at high redshift.
△ Less
Submitted 4 March, 2022;
originally announced March 2022.
-
A Parameter Space Exploration of High Resolution Numerically Evolved Early Type Galaxies Including AGN Feedback and Accurate Dynamical Treatment of Stellar Orbits
Authors:
Luca Ciotti,
Jeremiah P. Ostriker,
Zhaoming Gan,
Brian Xing Jiang,
Silvia Pellegrini,
Caterina Caravita,
Antonio Mancino
Abstract:
An extensive exploration of the model parameter space of axisymmetric Early-Type Galaxies (ETGs) hosting a central supermassive Black Hole (SMBH) is conducted by means of high resolution hydrodynamical simulations performed with our code MACER. Global properties such as 1) total SMBH accreted mass, 2) final X-ray luminosity and temperature of the X-ray emitting halos, 3) total amount of new stars…
▽ More
An extensive exploration of the model parameter space of axisymmetric Early-Type Galaxies (ETGs) hosting a central supermassive Black Hole (SMBH) is conducted by means of high resolution hydrodynamical simulations performed with our code MACER. Global properties such as 1) total SMBH accreted mass, 2) final X-ray luminosity and temperature of the X-ray emitting halos, 3) total amount of new stars formed from the cooling gas, 4) total ejected mass in form of supernovae and AGN feedback induced galactic winds, are obtained as a function of galaxy structure and internal dynamics. In addition to the galactic dark matter halo, the model galaxies are also embedded in a group/cluster dark matter halo; finally cosmological accretion is also included, with amount and time dependence derived from cosmological simulations. Angular momentum conservation leads to the formation of cold HI disks; these disks further evolve under the action of star formation induced by disk instabilities, of the associated mass discharge onto the central SMBH, and of the consequent AGN feedback. At the end of the simulations, the hot (metal enriched) gas mass is roughly $10\%$ the mass in the old stars, with twice as much having been ejected into the intergalactic medium. The cold gas disks are a $\approx$ kpc in size, and the metal rich new stars are in $0.1$ kpc disks. The masses of cold gas and new stars are roughly $0.1\%$ the mass of the old stars. Overall, the final systems appear to reproduce quite successfully the main global properties of real ETGs.
△ Less
Submitted 25 May, 2022; v1 submitted 11 January, 2022;
originally announced January 2022.
-
Black hole mergers in compact star clusters and massive black hole formation beyond the mass-gap
Authors:
Francesco Paolo Rizzuto,
Thorsten Naab,
Rainer Spurzem,
Manuel Arca-Sedda,
Mirek Giersz,
Jeremiah Paul Ostriker,
Sambaran Banerjee
Abstract:
We present direct N-body simulations, carried out with Nbody6++GPU, of young and compact low metallicity star clusters with $1.1\times 10^5$ stars, a velocity dispersion of $\sim$ 10 $\mathrm{km\,s^{-1}}$, a half mass radius $R_h=0.6$ pc, and a binary fraction of $10\%$ including updated evolution models for stellar winds and pair-instability supernovae (PISNe). Within the first tens of megayears…
▽ More
We present direct N-body simulations, carried out with Nbody6++GPU, of young and compact low metallicity star clusters with $1.1\times 10^5$ stars, a velocity dispersion of $\sim$ 10 $\mathrm{km\,s^{-1}}$, a half mass radius $R_h=0.6$ pc, and a binary fraction of $10\%$ including updated evolution models for stellar winds and pair-instability supernovae (PISNe). Within the first tens of megayears of evolution, each cluster hosts several black hole (BH) merger events which nearly cover the complete mass range of primary and secondary BH masses for current LIGO/Virgo/Kagra gravitational wave detections. The importance of gravitational recoil is estimated statistically. We present several possible formation paths of massive BHs above the assumed lower PISNe mass-gap limit ($45 M_\odot$) into the intermediate-mass BH (IMBH) regime ($> 100 M_\odot$) which include collisions of stars and BHs as well as the direct collapse of stellar merger remnants with low mass cores. The stellar evolution updates result in the early formation of higher mass stellar BHs than for the previous model. The resulting higher collision rates with massive stars support the rapid formation of massive BHs. For models assuming a high accretion efficiency for star-BH mergers, we present a first-generation formation scenario for GW190521-like events, a merger of two BHs in the PISN mass-gap, which is dominated by star-BH mergers. This IMBH formation path is independent of gravitational recoil and therefore conceivable in dense stellar systems with low escape velocities. One simulated cluster even forms an IMBH binary (153$M_\odot$,173$M_\odot$) which is expected to merge within a Hubble time.
△ Less
Submitted 25 August, 2021;
originally announced August 2021.
-
An Improved and Physically-Motivated Scheme for Matching Galaxies with Dark Matter Halos
Authors:
Stephanie Tonnesen,
Jeremiah P. Ostriker
Abstract:
The simplest scheme for predicting real galaxy properties after performing a dark matter simulation is to rank order the real systems by stellar mass and the simulated systems by halo mass and then simply assume monotonicity - that the more massive halos host the more massive galaxies. This has had some success, but we study here if a better motivated and more accurate matching scheme is easily co…
▽ More
The simplest scheme for predicting real galaxy properties after performing a dark matter simulation is to rank order the real systems by stellar mass and the simulated systems by halo mass and then simply assume monotonicity - that the more massive halos host the more massive galaxies. This has had some success, but we study here if a better motivated and more accurate matching scheme is easily constructed by looking carefully at how well one could predict the simulated IllustrisTNG galaxy sample from its dark matter computations. We find that using the dark matter rotation curve peak velocity, $v_{max}$, for normal galaxies reduces the error of the prediction by 30% (18% for central galaxies and 60% for satellite systems) - following expectations from the physics of monolithic collapse. For massive systems with halo mass $>$ 10$^{12.5}$ M$_{\odot}$ hierarchical merger driven formation is the better model and dark matter halo mass remains the best single metric. Using a new single variable that combines these effects, $φ$ $=$ $v_{max}$/$v_{max,12.7}$ + M$_{peak}$/(10$^{12.7}$ M$_{\odot}$) allows further improvement and reduces the error, as compared to ranking by dark matter mass at $z=0$ by another 6% from $v_{max}$ ranking. Two parameter fits -- including environmental effects produce only minimal further impact.
△ Less
Submitted 20 October, 2021; v1 submitted 25 February, 2021;
originally announced February 2021.
-
The Connection between Mergers and AGN Activity in Simulated and Observed Massive Galaxies
Authors:
Ray S. Sharma,
Ena Choi,
Rachel S. Somerville,
Gregory F. Snyder,
Dale D. Kocevski,
Michaela Hirschmann,
Benjamin P. Moster,
Thorsten Naab,
Desika Narayanan,
Jeremiah P. Ostriker,
David J. Rosario
Abstract:
We analyze a suite of $30$ high resolution zoom-in cosmological hydrodynamic simulations of massive galaxies with stellar masses $M_{\ast} > 10^{10.9} M_\odot$, with the goal of better understanding merger activity in AGN, AGN activity in merging systems, SMBH growth during mergers, and the role of gas content. Using the radiative transfer code \textsc{Powderday}, we generate HST-WFC3 F160W synthe…
▽ More
We analyze a suite of $30$ high resolution zoom-in cosmological hydrodynamic simulations of massive galaxies with stellar masses $M_{\ast} > 10^{10.9} M_\odot$, with the goal of better understanding merger activity in AGN, AGN activity in merging systems, SMBH growth during mergers, and the role of gas content. Using the radiative transfer code \textsc{Powderday}, we generate HST-WFC3 F160W synthetic observations of redshift $0.5 < z < 3$ central galaxies, add noise properties similar to the CANDELS survey, and measure morphological properties from the synthetic images using commonly adopted non-parametric statistics. We compare the distributions of morphological properties measured from the synthetic images with a sample of inactive galaxies and X-ray selected AGN hosts from CANDELS. We study the connection between mergers and AGN activity in the simulations, the synthetic images, and the observed CANDELS sample. We find that, in both the simulations and CANDELS, even the most luminous $(L_{\rm bol} > 10^{45}$ erg s$^{-1})$ AGN in our sample are no more likely than inactive galaxies $(L_{\rm bol} < 10^{43}$ erg s$^{-1})$ to be found in merging systems. We also find that AGN activity is not overall enhanced by mergers, nor enhanced at any specific time in the $1$ Gyr preceding and following a merger. Even gas rich major mergers (stellar mass ratio $>$1:4) do not necessarily enhance AGN activity or significantly grow the central SMBH. We conclude that in the simulated massive galaxies studied here, mergers are not the primary drivers of AGN.
△ Less
Submitted 18 December, 2023; v1 submitted 5 January, 2021;
originally announced January 2021.
-
Active Galactic Nucleus Feedback in an Elliptical Galaxy with the Most Updated AGN Physics: Parameter Explorations
Authors:
Zhiyuan Yao,
Feng Yuan,
Jeremiah P. Ostriker
Abstract:
In a previous work, we have proposed a sub-grid model of active galactic nucleus (AGN) feedback by taking into account the state-of-the-art AGN physics, and used that model to study the effect of AGN feedback on the evolution of an isolated elliptical galaxy by performing two dimensional high-resolution (i.e., the Bondi radius is well resolved) simulations. In that work, typical values of model pa…
▽ More
In a previous work, we have proposed a sub-grid model of active galactic nucleus (AGN) feedback by taking into account the state-of-the-art AGN physics, and used that model to study the effect of AGN feedback on the evolution of an isolated elliptical galaxy by performing two dimensional high-resolution (i.e., the Bondi radius is well resolved) simulations. In that work, typical values of model parameters were adopted. In the present work, we extend that study by exploring the effects of uncertainties of parameter values. Such a study is also useful for us to understand the respective roles of various components of the model. These parameters include the mass flux and velocity of AGN wind and radiative efficiency in both the hot and cold feedback modes, and the initial black hole (BH) mass. We find that the velocity of AGN wind in the hot mode is the most important quantity to control the typical accretion rate and luminosity of AGN, and the mass growth of the BH. The effect of the wind on star formation is less sensitive. Within the limited parameter range explored in the current work, a stronger AGN wind suppresses star formation within ~100 pc but enhances star formation beyond this radius, while the star formation integrated over the evolution time and the whole galaxy roughly remain unchanged. AGN radiation suppresses the BH accretion in a mild way, but dust is not considered here. Finally, a smaller initial BH mass results in a more violent evolution of the BH accretion rate. The corresponding AGN spends more time in the high-luminosity state and the percentage of BH mass growth is higher. Our results indicate the robustness of AGN feedback in keeping the galaxy quenched.
△ Less
Submitted 30 November, 2020;
originally announced November 2020.
-
What is the Halo Mass Function in a Fuzzy Dark Matter Cosmology?
Authors:
Mihir Kulkarni,
Jeremiah P. Ostriker
Abstract:
Fuzzy dark matter (FDM) or wave dark matter is an alternative theory designed to solve the small-scale problems faced by the standard cold dark matter proposal for the primary material component of the universe. It is made up of ultra-light axions having mass $\sim 10^{-22}$ eV that typically have de Broglie wavelength of several kpc, alleviating some of the apparent small-scale discrepancies face…
▽ More
Fuzzy dark matter (FDM) or wave dark matter is an alternative theory designed to solve the small-scale problems faced by the standard cold dark matter proposal for the primary material component of the universe. It is made up of ultra-light axions having mass $\sim 10^{-22}$ eV that typically have de Broglie wavelength of several kpc, alleviating some of the apparent small-scale discrepancies faced by the standard $Λ$CDM paradigm. In this paper, we calculate the halo mass function for the fuzzy dark matter using a sharp-k window function and compare it with one calculated using numerical simulations, finding the peak mass at roughly $10^{10} {M_{\odot}}$ for a particle mass of $2\times 10^{-22}$ eV. We also constrain the mass of FDM particle to be $\gtrapprox 2\times10^{-22}$ eV using the observations of high-redshift ($z\sim 10$) lensed galaxies from CLASH survey.
△ Less
Submitted 6 December, 2021; v1 submitted 3 November, 2020;
originally announced November 2020.
-
The Impact of Outflows driven by Active Galactic Nuclei on Metals in and around Galaxies
Authors:
Ena Choi,
Ryan Brennan,
Rachel S. Somerville,
Jeremiah P. Ostriker,
Michaela Hirschmann,
Thorsten Naab
Abstract:
Metals in the hot gaseous halos of galaxies encode the history of star formation as well as the feedback processes that eject metals from the galaxies. X-ray observations suggest that massive galaxies have extended distributions of metals in their gas halos. We present predictions for the metal properties of massive galaxies and their gaseous halos from recent high resolution zoom-in simulations t…
▽ More
Metals in the hot gaseous halos of galaxies encode the history of star formation as well as the feedback processes that eject metals from the galaxies. X-ray observations suggest that massive galaxies have extended distributions of metals in their gas halos. We present predictions for the metal properties of massive galaxies and their gaseous halos from recent high resolution zoom-in simulations that include mechanical and radiation driven feedback from Active Galactic Nuclei (AGN). In these simulations, AGN launch high-velocity outflows, mimicking observed broad absorption line winds. By comparing two sets of simulations with and without AGN feedback, we show that our prescription for AGN feedback is capable of driving winds and enriching halo gas `inside-out' by spreading centrally enriched metals to the outskirts of galaxies, into the halo and beyond. The metal (iron) profiles of halos simulated with AGN feedback have a flatter slope than those without AGN feedback, consistent with recent X-ray observations. The predicted gas iron abundance of group scale galaxies simulated with AGN feedback is $Z_{\rm Fe} = 0.23$ $Z_{\rm Fe,\odot}$ at $0.5 r_{500}$, which is 2.5 times higher than that in simulations without AGN feedback. In these simulations, AGN winds are also important for the metal enrichment of the intergalactic medium, as the vast majority of metals ejected from the galaxy by AGN-driven winds end up beyond the halo virial radius.
△ Less
Submitted 21 September, 2020;
originally announced September 2020.
-
Intermediate Mass Black Hole Formation in compact Young Massive Star Clusters
Authors:
Francesco Paolo Rizzuto,
Thorsten Naab,
Rainer Spurzem,
Mirek Giersz,
J. P. Ostriker,
N. C. Stone,
Long Wang,
Peter Berczik,
M. Rampp
Abstract:
Young dense massive star clusters are a promising environment for the formation of intermediate mass black holes (IMBHs) through collisions. We present a set of 80 simulations carried out with Nbody6++GPU of 10 initial conditions for compact $\sim 7 \times 10^4 M_{\odot}$ star clusters with half-mass radii $R_\mathrm{h} \lesssim 1 pc$, central densities…
▽ More
Young dense massive star clusters are a promising environment for the formation of intermediate mass black holes (IMBHs) through collisions. We present a set of 80 simulations carried out with Nbody6++GPU of 10 initial conditions for compact $\sim 7 \times 10^4 M_{\odot}$ star clusters with half-mass radii $R_\mathrm{h} \lesssim 1 pc$, central densities $ρ_\mathrm{core} \gtrsim 10^5 M_\odot pc^{-3}$, and resolved stellar populations with 10\% primordial binaries. Very massive stars (VMSs) with masses up to $\sim 400 M_\odot$ grow rapidly by binary exchange and three-body scattering events with main sequences stars in hard binaries. Assuming that in VMS - stellar BH collisions all stellar material is accreted onto the BH, IMBHs with masses up to $M_\mathrm{BH} \sim 350 M_\odot$ can form on timescales of $\lesssim 15$ Myr. This process was qualitatively predicted from Monte Carlo MOCCA simulations. Despite the stochastic nature of the process - typically not more than 3/8 cluster realisations show IMBH formation - we find indications for higher formation efficiencies in more compact clusters. Assuming a lower accretion fraction of 0.5 for VMS - BH collisions, IMBHs can also form. The process might not work for accretion fractions as low as 0.1. After formation, the IMBHs can experience occasional mergers with stellar mass BHs in intermediate mass-ratio inspiral events (IMRIs) on a 100 Myr timescale. Realised with more than $10^5$ stars, 10 \% binaries, the assumed stellar evolution model with all relevant evolution processes included and 300 Myr simulation time, our large suite of simulations indicates that IMBHs of several hundred solar masses might form rapidly in massive star clusters right after their birth while they are still compact.
△ Less
Submitted 21 August, 2020;
originally announced August 2020.
-
Infra-Red Emission from Cold Gas Dusty Disks in Massive Ellipticals
Authors:
Zhaoming Gan,
Brandon S. Hensley,
Jeremiah P. Ostriker,
Luca Ciotti,
David Schiminovich,
Silvia Pellegrini
Abstract:
What is the expected infrared output of elliptical galaxies? Here we report the latest findings obtained in this high time resolution (~10 years) and high spatial resolution (2.5 parsec at center) study. We add a set of grain physics to the MACER code, including (a) dust grains made in passive stellar evolution; (b) dust grain growth due to collision and sticking; (c) grain destruction due to ther…
▽ More
What is the expected infrared output of elliptical galaxies? Here we report the latest findings obtained in this high time resolution (~10 years) and high spatial resolution (2.5 parsec at center) study. We add a set of grain physics to the MACER code, including (a) dust grains made in passive stellar evolution; (b) dust grain growth due to collision and sticking; (c) grain destruction due to thermal sputtering; (d) dust cooling of hot gas via inelastic collisions; and (e) radiation pressure on dust grains. The code improvements enable us to analyze the effects of dust on metal depletion and AGN obscuration, and also to assess the infrared output of the modeled galaxies. We simulate a representative massive elliptical galaxy of a central stellar velocity dispersion ~ 260 km/s and modest rotation. We find that: (1) the circumnuclear disk (of a size <~ 1 kpc) is dusty in its outer region where most of the metals are in dust grains, while in the inner disk most of the dust grains are destroyed by the AGN irradiation; (2) the dusty disk is optically thick to both the starlight within the disk and the radiation from the central AGN. Thus the AGN will be obscured behind the disk, and the latter is of a covering factor ~ 0.2; (3) the dust infrared emission is mainly due to the AGN irradiation. The median infrared luminosity is ~ 2e44 erg/s, and it can reach >~ 1e46 erg/s during outbursts; (4) the duty cycles of the AGN activities, star formation, and the dust infrared luminosity roughly match observations, e.g., in most of its lifetime, the simulated galaxy is a stereotypical "quiescent" elliptical galaxy with L_{IR} ~ 1e11*L_{solar}, while the star formation rate can exceed 250 M_{solar}/yr during central outbursts.
△ Less
Submitted 17 August, 2020; v1 submitted 12 December, 2019;
originally announced December 2019.
-
Metal abundances in the MACER simulations of the hot interstellar medium
Authors:
S. Pellegrini,
Z. Gan,
J. P. Ostriker,
L. Ciotti
Abstract:
A hot plasma is the dominant phase of the interstellar medium of early-type galaxies. Its origin can reside in stellar mass losses, residual gas from the formation epoch, and accretion from outside of the galaxies. Its evolution is linked to the dynamical structure of the host galaxy, to the supernova and AGN feedback, and to (late-epoch) star formation, in a way that has yet to be fully understoo…
▽ More
A hot plasma is the dominant phase of the interstellar medium of early-type galaxies. Its origin can reside in stellar mass losses, residual gas from the formation epoch, and accretion from outside of the galaxies. Its evolution is linked to the dynamical structure of the host galaxy, to the supernova and AGN feedback, and to (late-epoch) star formation, in a way that has yet to be fully understood. Important clues about the origin and evolution of the hot gas come from the abundances of heavy metals, that have been studied with increasing detail with XMM-Newton and Chandra. We present recent high resolution hydrodynamical simulations of the hot gas evolution that include the above processes, and where several chemical species, originating in AGB stars and supernovae of type Ia and II, have also been considered. The high resolution, of few parsecs in the central galactic region, allows us to track the metal enrichment, transportation and dilution throughout the galaxy. The comparison of model results with observed abundances reveals a good agreement for the region enriched by the AGN wind, but also discrepancies for the diffuse hot gas; the latter indicate the need for a revision of standard assumptions, and/or the importance of neglected effects as those due to the dust, and/or residual uncertainties in deriving abundances from the X-ray spectra.
△ Less
Submitted 12 December, 2019;
originally announced December 2019.
-
Radiation-driven outflows in AGNs: Revisiting feedback effects of scattered and reprocessed photons
Authors:
Amin Mosallanezhad,
Feng Yuan,
Jeremiah P. Ostriker,
Fatemeh Zahra Zeraatgari,
De-Fu Bu
Abstract:
We perform two-dimensional hydrodynamical simulations of slowly rotating accretion flows in the region of $ 0.01-7\, \mathrm{pc} $ around a supermassive black holes with $ M_\mathrm{BH} = 10^{8} M_{\odot} $. The accretion flow is irradiated by the photons from the central active galactic nucleus (AGN). In addition to the direct radiation from the AGN, we have also included the "re-radiation", i.e.…
▽ More
We perform two-dimensional hydrodynamical simulations of slowly rotating accretion flows in the region of $ 0.01-7\, \mathrm{pc} $ around a supermassive black holes with $ M_\mathrm{BH} = 10^{8} M_{\odot} $. The accretion flow is irradiated by the photons from the central active galactic nucleus (AGN). In addition to the direct radiation from the AGN, we have also included the "re-radiation", i.e., the locally produced radiation by Thomson scattering, line and bremsstrahlung radiation. Compare to our previous work, we have improved the calculation of radiation force due to the Thomson scattering of X-ray photons from the central AGN. We find that this improvement can significantly increase the mass flux and velocity of outflow. We have compared the properties of outflow --- including mass outflow rate, velocity, and kinetic luminosity of outflow --- in our simulation with the observed properties of outflow in AGNs and found that they are in good consistency. This implies that the combination of line and re-radiation forces is the possible origin of observed outflow in luminous AGNs.
△ Less
Submitted 1 October, 2019;
originally announced October 2019.
-
NIHAO-UHD: The properties of MW-like stellar disks in high resolution cosmological simulations
Authors:
Tobias Buck,
Aura Obreja,
Andrea V. Macciò,
Ivan Minchev,
Aaron A. Dutton,
Jeremiah P. Ostriker
Abstract:
Simulating thin and extended galactic disks has long been a challenge in computational astrophysics. We introduce the NIHAO-UHD suite of cosmological hydrodynamical simulations of Milky Way mass galaxies and study stellar disk properties such as stellar mass, size and rotation velocity which agree well with observations of the Milky Way and local galaxies. In particular, the simulations reproduce…
▽ More
Simulating thin and extended galactic disks has long been a challenge in computational astrophysics. We introduce the NIHAO-UHD suite of cosmological hydrodynamical simulations of Milky Way mass galaxies and study stellar disk properties such as stellar mass, size and rotation velocity which agree well with observations of the Milky Way and local galaxies. In particular, the simulations reproduce the age-velocity dispersion relation and a multi-component stellar disk as observed for the Milky Way. Half of our galaxies show a double exponential vertical profile, while the others are well described by a single exponential model which we link to the disk merger history. In all cases, mono-age populations follow a single exponential whose scale height varies monotonically with stellar age and radius. The scale length decreases with stellar age while the scale height increases. The general structure of the stellar disks is already set at time of birth as a result of the inside-out and upside-down formation. Subsequent evolution modifies this structure by increasing both the scale length and height of all mono-age populations. Thus, our results put tight constraints on how much dynamical memory stellar disks can retain over cosmological timescales. Our simulations demonstrate that it is possible to form thin galactic disks in cosmological simulations provided there are no significant stellar mergers at low redshifts. Most of the stellar mass is formed in-situ with only a few percent ($\lesssim5\%$) brought in by merging satellites at early times. Redshift zero snapshots and halo catalogues are publicly available.
△ Less
Submitted 6 January, 2020; v1 submitted 12 September, 2019;
originally announced September 2019.
-
AGN-driven quenching of satellite galaxies
Authors:
Gohar Dashyan,
Ena Choi,
Rachel S. Somerville,
Thorsten Naab,
Amanda C. N. Quirk,
Michaela Hirschmann,
Jeremiah P. Ostriker
Abstract:
We explore the effect of active galactic nucleus (AGN) feedback from central galaxies on their satellites by comparing two sets of cosmological zoom-in runs of 27 halos with masses ranging from $10^{12}$ to $10^{13.4}$ solar masses at z=0, with (wAGN) and without (noAGN) AGN feedback. Both simulations include stellar feedback from multiple processes, including powerful winds from supernovae, stell…
▽ More
We explore the effect of active galactic nucleus (AGN) feedback from central galaxies on their satellites by comparing two sets of cosmological zoom-in runs of 27 halos with masses ranging from $10^{12}$ to $10^{13.4}$ solar masses at z=0, with (wAGN) and without (noAGN) AGN feedback. Both simulations include stellar feedback from multiple processes, including powerful winds from supernovae, stellar winds from young massive stars, AGB stars, radiative heating within Strömgren spheres and photoelectric heating. Our wAGN model is identical to the noAGN model except that it also includes a model for black hole seeding and accretion, as well as AGN feedback via high-velocity broad absorption line winds and Compton/photoionization heating. We show that the inclusion of AGN feedback from the central galaxy significantly affects the star formation history and the gas content of the satellite galaxies. AGN feedback starts to affect the gas content and the star formation of the satellites as early as z=2. The mean gas rich fraction of satellites at z=0 decreases from 15% in the noAGN simulation to 5% in the wAGN simulation. The difference between the two sets extends as far out as five times the virial radius of the central galaxy at z=1. We investigate the quenching mechanism by studying the physical conditions in the surroundings of pairs of satellites matched across the wAGN and noAGN simulations and find an increase in the temperature and relative velocity of the intergalactic gas.
△ Less
Submitted 18 June, 2019;
originally announced June 2019.
-
MIND THE GAP: Is The Too Big To Fail Problem Resolved?
Authors:
Jeremiah P. Ostriker,
Ena Choi,
Anthony Chow,
Kundan Guha
Abstract:
The faintness of satellite systems in galaxy groups has contributed to the widely discussed "missing satellite" and "too big to fail" issues. Using techniques based on Tremaine & Richstone (1977), we show that there is no problem with the luminosity function computed from modern codes per se, but that the gap between first and second brightest systems is too big {\it given} the luminosity function…
▽ More
The faintness of satellite systems in galaxy groups has contributed to the widely discussed "missing satellite" and "too big to fail" issues. Using techniques based on Tremaine & Richstone (1977), we show that there is no problem with the luminosity function computed from modern codes per se, but that the gap between first and second brightest systems is too big {\it given} the luminosity function, that the same large gap is found in modern, large scale baryonic $Λ$CDM simulations such as EAGLE and IllustrisTNG, is even greater in dark matter only simulations, and finally, that this is most likely due to gravitationally induced merging caused by classical dynamical friction. Quantitatively the gap is larger in the computed simulations than in the randomized ones by $1.79 \pm 1.04$, $1.51 \pm 0.93$, $3.43 \pm 1.44$ and $3.33 \pm 1.35$ magnitudes in the EAGLE, IllustrisTNG, and dark matter only simulations of EAGLE and IllustrisTNG respectively. Furthermore the anomalous gaps in the simulated systems are even larger than in the real data by over half a magnitude and are still larger in the dark matter only simulations. Briefly stated, $Λ$CDM does not have a problem with an absence of "too big to fail" galaxies. Statistically significant large gaps between first and second brightest galaxies are to be expected.
△ Less
Submitted 16 July, 2019; v1 submitted 23 April, 2019;
originally announced April 2019.
-
Transition of BH feeding from the quiescent regime into star-forming cold disk regime
Authors:
Kohei Inayoshi,
Kohei Ichikawa,
Jeremiah P. Ostriker,
Rolf Kuiper
Abstract:
We study the properties of rotating accretion flows onto supermassive black holes (SMBHs) using axisymmetric two-dimensional hydrodynamical simulations with radiative cooling and BH feedback. The simulations resolve the accretion dynamics of gas outside from the BH influence radius through an inner accretion disk. For lower Bondi accretion rates in units of the Eddington rate (…
▽ More
We study the properties of rotating accretion flows onto supermassive black holes (SMBHs) using axisymmetric two-dimensional hydrodynamical simulations with radiative cooling and BH feedback. The simulations resolve the accretion dynamics of gas outside from the BH influence radius through an inner accretion disk. For lower Bondi accretion rates in units of the Eddington rate ($\dot{M}_{\rm B}\ll 10^{-3}~\dot{M}_{\rm Edd}$), the BH feeding is suppressed due to turbulent motion by several orders of magnitudes from the Bondi rate. Thus, the radiative luminosity results in as low as $\sim 10^{-10}-10^{-7}~L_{\rm Edd}$, where $L_{\rm Edd}$ is the Eddington luminosity. For higher rates of $\dot{M}_{\rm B}> 10^{-3}~\dot{M}_{\rm Edd}$, the optically-thin accreting gas cools via free-free emission and forms a geometrically-thin disk, which feeds the BH efficiently and increases the radiative luminosity to $> 10^{-3}~L_{\rm Edd}$. The transitional behavior of accreting BHs in galactic nuclei from radiatively inefficient phases to cold disk accretion naturally explains (1) the reason for the offset between the observed luminosities and theoretical predictions for nearby quiescent SMBHs, and (2) the conditions to fuel gas into the nuclear SMBH. In addition, the cold disk formed in galactic nuclei tends to be gravitationally unstable and leads to star formation when the Bondi rate is as high as $ \dot{M}_{\rm B} > 10^{-2}~M_\odot~{\rm yr}^{-1}$. This is a plausible explanation of the correlation observed between star formation rates and BH feeding rates in Seyfert galaxies.
△ Less
Submitted 15 May, 2019; v1 submitted 6 February, 2019;
originally announced February 2019.
-
On the Role of Hot Feedback Mode in Active Galactic Nuclei Feedback in an Elliptical Galaxy
Authors:
Doosoo Yoon,
Feng Yuan,
Jeremiah P. Ostriker,
Luca Ciotti,
Bocheng Zhu
Abstract:
Depending on the value of the accretion rate, black hole accretion is divided into cold and hot modes. The two modes have distinctly different physics and correspond to two feedback modes. Most previous feedback works either focus only on one mode, or the accretion physics is not always properly adopted. Here, by performing hydrodynamical numerical simulations of AGN feedback in an elliptical gala…
▽ More
Depending on the value of the accretion rate, black hole accretion is divided into cold and hot modes. The two modes have distinctly different physics and correspond to two feedback modes. Most previous feedback works either focus only on one mode, or the accretion physics is not always properly adopted. Here, by performing hydrodynamical numerical simulations of AGN feedback in an elliptical galaxy, we show that including both is important, and gives different results from including just one or the other. We specifically focus on the wind and radiation (but neglecting the jet) feedback in the hot mode, to explore whether this particular mode of feedback can play any important feedback role. For this aim, we have run two test models. In one model, we always adopt the cold mode no matter what value the accretion rate is; in another model, we turn off the AGN once it enters into the hot mode. We have calculated the AGN light curves, black hole growth, star formation, and AGN energy duty-cycle; and compared the results with the model in which the two modes are correctly included. Important differences are found. For example, if we were to adopt only the cold mode, the total mass of newly formed stars would become two orders of magnitude smaller, and the fraction of energy ejected within the high accretion regime (i.e., $L_{\rm BH} > 2\% L_{\rm Edd}$) would be too small to be consistent with observations. We have also investigated the respective roles of wind and radiation in the hot mode.
△ Less
Submitted 17 September, 2019; v1 submitted 22 January, 2019;
originally announced January 2019.
-
Adding a Suite of Chemical Abundances to the MACER Code for the Evolution of Massive Elliptical Galaxies
Authors:
Zhaoming Gan,
Ena Choi,
Jeremiah P. Ostriker,
Luca Ciotti,
Silvia Pellegrini
Abstract:
We add a suite of chemical abundances to the MACER (Massive AGN Controlled Ellipticals Resolved) 2D code, by solving 12 additional continuity equations for H, He, C, N, O, Ne, Mg, Si, S, Ca, Fe and Ni respectively with sources from AGB stars and supernovae of type Ia and II with metal yields based on standard stellar physics. New stars, formed in Toomre unstable circumnuclear disks (of a size…
▽ More
We add a suite of chemical abundances to the MACER (Massive AGN Controlled Ellipticals Resolved) 2D code, by solving 12 additional continuity equations for H, He, C, N, O, Ne, Mg, Si, S, Ca, Fe and Ni respectively with sources from AGB stars and supernovae of type Ia and II with metal yields based on standard stellar physics. New stars, formed in Toomre unstable circumnuclear disks (of a size $\lesssim150$ parsec), are assumed to have a top-heavy initial mass function with a power index of 1.65. The metal dilution effects due to cosmic accretion are also included. With the high resolution of few parsecs in central regions, resolved black hole accretion and AGN feedback, we can track the metal enrichment, transportation and dilution throughout the modeled massive elliptical galaxy of velocity dispersion $\sim280$ km/s. We retrieve the chemical composition of the BAL winds launched by the central AGN, synthesize the X-ray features of the hot ISM, and find that (1) the simulated metallicity in the BAL winds could be up to $\sim 8 Z_\odot$, while that of the hot ISM in the host galaxy is $\sim 2.3 Z_\odot$, matching well with SDSS observations of BLR gas; (2) the X-ray emitting hot gas is metal enriched with a typical value $\sim2.5 Z_\odot$; (3) the circumunuclear cold gas disk, where the metals are condensed, further enriched and recycled, plays a critical role in the metal enrichment; (4) the black hole accretion rate $\dot{M}_{\rm BH}$ linearly correlates with the star formation rate $\dot{M}_\star^+$ in the circumnuclear disk, i.e, $\dot{M}_\star^+ \sim 7.7 \dot{M}_{\rm BH}$, but lagged in time by roughly $10^6$ years.
△ Less
Submitted 8 April, 2019; v1 submitted 25 December, 2018;
originally announced December 2018.
-
Abundance matching for low mass galaxies in the CDM and FDM models
Authors:
Pierre Cristofari,
Jeremiah P. Ostriker
Abstract:
Abundance matching between galaxies and halos has proven to be an informative technique, less dependent on detailed physical approximations than some other methods. We extend the discussion to the dwarf galaxies realm and to the study of the alternative hypothesis that the dark matter is composed of ultra--light particles: Fuzzy Dark Matter. We find that, given current observations, both CDM and F…
▽ More
Abundance matching between galaxies and halos has proven to be an informative technique, less dependent on detailed physical approximations than some other methods. We extend the discussion to the dwarf galaxies realm and to the study of the alternative hypothesis that the dark matter is composed of ultra--light particles: Fuzzy Dark Matter. We find that, given current observations, both CDM and FDM have difficulties with dwarf systems in the local group, but that, if mechanisms are proposed to alleviate these difficulties, they make radically different predictions for FDM and CDM at high redshift. These differences should allow clear observational tests with currently planned experiments, with the number of galaxies per unit volume with stellar mass over $10^6$ M$_{\odot}$ lower in FDM by factors of ($10^{-2.7},10^{-2.0},10^{-1.7},10^{-0.7},10^{-0.1}$) at $z=(10,8,6,4,2)$.
△ Less
Submitted 6 November, 2018; v1 submitted 30 October, 2018;
originally announced October 2018.
-
Heating of Milky Way disc Stars by Dark Matter Fluctuations in Cold Dark Matter and Fuzzy Dark Matter Paradigms
Authors:
Benjamin V. Church,
Jeremiah P. Ostriker,
Philip Mocz
Abstract:
Although highly successful on cosmological scales, Cold Dark Matter (CDM) models predict unobserved over-dense `cusps' in dwarf galaxies and overestimate their formation rate. We consider an ultra-light axion-like scalar boson which promises to reduce these observational discrepancies at galactic scales. The model, known as Fuzzy Dark Matter (FDM), avoids cusps, suppresses small-scale power, and d…
▽ More
Although highly successful on cosmological scales, Cold Dark Matter (CDM) models predict unobserved over-dense `cusps' in dwarf galaxies and overestimate their formation rate. We consider an ultra-light axion-like scalar boson which promises to reduce these observational discrepancies at galactic scales. The model, known as Fuzzy Dark Matter (FDM), avoids cusps, suppresses small-scale power, and delays galaxy formation via macroscopic quantum pressure. We compare the substructure and density fluctuations of galactic dark matter haloes comprised of ultra-light axions to conventional CDM results. Besides self-gravitating subhaloes, FDM includes non-virialized over-dense wavelets formed by quantum interference patterns which are an efficient source of heating to galactic discs. We find that, in the solar neighborhood, wavelet heating is sufficient to give the oldest disc stars a velocity dispersion of $\sim 30 \: \mathrm{km} \: \mathrm{s}^{-1}$ within a Hubble time if energy is not lost from the disc, the velocity dispersion increasing with stellar age as $σ_D \propto t^{0.4}$ in agreement with observations. Furthermore, we calculate the radius-dependent velocity dispersion and corresponding scale height caused by the heating of this dynamical substructure in both CDM and FDM with the determination that these effects will produce a flaring that terminates the Milky Way disc at $15 - 20 \: \mathrm{kpc}$. Although the source of thickened discs is not known, the heating due to perturbations caused by dark substructure cannot exceed the total disc velocity dispersion. Therefore, this work provides a lower bound on the FDM particle mass of $m_a > 0.6 \times 10^{-22} \mathrm{eV}$. Furthermore, FDM wavelets with this particle mass should be considered a viable mechanism for producing the observed disc thickening with time.
△ Less
Submitted 20 February, 2019; v1 submitted 12 September, 2018;
originally announced September 2018.
-
MACER improved: AGN feedback computed in rotating early-type galaxies at high resolution
Authors:
Zhaoming Gan,
Luca Ciotti,
Jeremiah P. Ostriker,
Feng Yuan
Abstract:
Based on our previous modeling of AGN feedback in isolated elliptical galaxies (Gan et al. 2014) using the MACER (Massive AGN Controlled Ellipticals Resolved) code, we extend and improve the model to include rotation, to limit star formation to regions of high density and low temperature, to facilitate angular momentum transfer via the Toomre instability in gaseous disks, and to improve the treatm…
▽ More
Based on our previous modeling of AGN feedback in isolated elliptical galaxies (Gan et al. 2014) using the MACER (Massive AGN Controlled Ellipticals Resolved) code, we extend and improve the model to include rotation, to limit star formation to regions of high density and low temperature, to facilitate angular momentum transfer via the Toomre instability in gaseous disks, and to improve the treatment of hot mode (low accretion rate) AGN feedback. The model galaxy now has an extended dark matter profile that matches with standard observations, but has a resolution of parsecs in the inner region and resolves the Bondi radius. We find that the results agree reasonably well with a panoply of observations: (1) both AGN activity and star formation are primarily in central cold gaseous disks, are bursty and mainly driven by the Toomre instability; (2)AGN duty cycle agrees well with the Soltan argument, i.e., the AGN spends most of its lifetime when it is in low luminosity (half of time with $L/L_{Edd}<7\times10^{-5}$), while emitting most of its energy when it is in high luminosity (half of radiant energy emitted with $L/L_{Edd}>0.06$); (3) the total star formation is $\sim$ few percents of the initial stellar mass, occurring in the bursts that would be associated with the observed E+A phenomenon. Most of the star formation occurs in the circumnuclear disk of a size <1 kpc, which is in agreement with recent observations; (4) the ISM X-ray luminosity varies within a reasonable range (median $L_{\rm X, ISM}=9.1\times10^{39}$ erg/s) in agreement with observations.
△ Less
Submitted 25 January, 2019; v1 submitted 11 September, 2018;
originally announced September 2018.
-
The Role of Black Hole Feedback on Size and Structural Evolution in Massive Galaxies
Authors:
Ena Choi,
Rachel S. Somerville,
Jeremiah P. Ostriker,
Thorsten Naab,
Michaela Hirschmann
Abstract:
We use cosmological hydrodynamical simulations to investigate the role of feedback from accreting black holes on the evolution of sizes, compactness, stellar core density and specific star-formation of massive galaxies with stellar masses of $ M_{star} > 10^{10.9} M_{\odot}$. We perform two sets of cosmological zoom-in simulations of 30 halos to z=0: (1) without black holes and Active Galactic Nuc…
▽ More
We use cosmological hydrodynamical simulations to investigate the role of feedback from accreting black holes on the evolution of sizes, compactness, stellar core density and specific star-formation of massive galaxies with stellar masses of $ M_{star} > 10^{10.9} M_{\odot}$. We perform two sets of cosmological zoom-in simulations of 30 halos to z=0: (1) without black holes and Active Galactic Nucleus (AGN) feedback and (2) with AGN feedback arising from winds and X-ray radiation. We find that AGN feedback can alter the stellar density distribution, reduce the core density within the central 1 kpc by 0.3 dex from z=1, and enhance the size growth of massive galaxies. We also find that galaxies simulated with AGN feedback evolve along similar tracks to those characterized by observations in specific star formation versus compactness. We confirm that AGN feedback plays an important role in transforming galaxies from blue compact galaxies into red extended galaxies in two ways: (1) it effectively quenches the star formation, transforming blue compact galaxies into compact quiescent galaxies and (2) it also removes and prevents new accretion of cold gas, shutting down in-situ star formation and causing subsequent mergers to be gas-poor or mixed. Gas poor minor mergers then build up an extended stellar envelope. AGN feedback also puffs up the central region through the fast AGN driven winds as well as the slow expulsion of gas while the black hole is quiescent. Without AGN feedback, large amounts of gas accumulate in the central region, triggering star formation and leading to overly massive blue galaxies with dense stellar cores.
△ Less
Submitted 6 September, 2018;
originally announced September 2018.
-
Numerical study of active galactic nucleus feedback in an elliptical galaxy with {\it MACER}
Authors:
Feng Yuan,
Jeremiah P. Ostriker,
DooSoo Yoon,
Ya-Ping Li,
Luca Ciotti,
Zhao-Ming Gan,
Luis C. Ho,
Fulai Guo
Abstract:
This paper summarizes our recent works of studying AGN feedback in an isolated elliptical galaxy by performing high-resolution hydrodynamical numerical simulations. Bondi radius is resolved and the mass accretion rate of the black hole is calculated. The most updated AGN physics, namely the discrimination of cold and hot accretion modes and the exact descriptions of the AGN radiation and wind for…
▽ More
This paper summarizes our recent works of studying AGN feedback in an isolated elliptical galaxy by performing high-resolution hydrodynamical numerical simulations. Bondi radius is resolved and the mass accretion rate of the black hole is calculated. The most updated AGN physics, namely the discrimination of cold and hot accretion modes and the exact descriptions of the AGN radiation and wind for a given accretion rate are adopted and their interaction with the gas in the host galaxy is calculated. Physical processes such as star formation and SNe feedback are taken into account. Consistent with observation, we find the AGN spends most of the time in the low-luminosity regime. AGN feedback overall suppresses the star formation; but depending on location in the galaxy and time, it can also enhance it. The light curve of specific star formation rate is not synchronous with the AGN light curve. We find that wind usually plays a dominant role in controlling the AGN luminosity and star formation, but radiation also cannot be neglected.
△ Less
Submitted 15 July, 2018;
originally announced July 2018.
-
Momentum-driven Winds from Radiatively Efficient Black Hole Accretion and Their Impact on Galaxies
Authors:
Ryan Brennan,
Ena Choi,
Rachel S. Somerville,
Michaela Hirschmann,
Thorsten Naab,
Jeremiah P. Ostriker
Abstract:
We explore the effect of momentum-driven winds representing radiation pressure driven outflows from accretion onto supermassive black holes in a set of numerical hydrodynamical simulations. We explore two matched sets of cosmological zoom-in runs of 24 halos with masses ~$10^{12.0}-10^{13.4}$ M_sun run with two different feedback models. Our `NoAGN' model includes stellar feedback via UV heating,…
▽ More
We explore the effect of momentum-driven winds representing radiation pressure driven outflows from accretion onto supermassive black holes in a set of numerical hydrodynamical simulations. We explore two matched sets of cosmological zoom-in runs of 24 halos with masses ~$10^{12.0}-10^{13.4}$ M_sun run with two different feedback models. Our `NoAGN' model includes stellar feedback via UV heating, stellar winds and supernovae, photoelectric heating and cosmic X-ray background heating from a meta-galactic background. Our fiducial `MrAGN' model is identical except that it also includes a model for black hole seeding and accretion, as well as heating and momentum injection associated with the radiation from black hole accretion. Our MrAGN model launches galactic outflows which result in both `ejective' feedback - the outflows themselves which drive gas out of galaxies - and `preventative' feedback, which suppresses the inflow of new and recycling gas. As much as 80 % of outflowing galactic gas can be expelled, and accretion can be suppressed by as much as a factor of 30 in the MrAGN runs when compared with the NoAGN runs. The histories of NoAGN galaxies are recycling-dominated, with ~70% of material that leaves the galaxy eventually returning, and the majority of outflowing gas re-accretes on 1 Gyr timescales without AGN feedback. Outflowing gas in the MrAGN runs has higher characteristic velocity (500 - 1,000 km/s versus 100-300 km/s for outflowing NoAGN gas) and travels as far as a few Mpcs. Only ~10% of ejected material is re-accreted in the MrAGN galaxies.
△ Less
Submitted 2 May, 2018;
originally announced May 2018.
-
An Overabundance of Black Hole X-Ray Binaries in the Galactic Center from Tidal Captures
Authors:
A. Generozov,
N. C. Stone,
B. D. Metzger,
J. P. Ostriker
Abstract:
A large population of X-ray binaries (XRBs) was recently discovered within the central parsec of the Galaxy by Hailey et al. While the presence of compact objects on this scale due to radial mass segregation is, in itself, unsurprising, the fraction of binaries would naively be expected to be small because of how easily primordial binaries are dissociated in the dynamically hot environment of the…
▽ More
A large population of X-ray binaries (XRBs) was recently discovered within the central parsec of the Galaxy by Hailey et al. While the presence of compact objects on this scale due to radial mass segregation is, in itself, unsurprising, the fraction of binaries would naively be expected to be small because of how easily primordial binaries are dissociated in the dynamically hot environment of the nuclear star cluster (NSC). We propose that the formation of XRBs in the central parsec is dominated by the tidal capture of stars by black holes (BHs) and neutron stars (NSs). We model the time-dependent radial density profiles of stars and compact objects in the NSC with a Fokker-Planck approach, using the present-day stellar population and rate of in situ massive star (and thus compact object) formation as observational constraints. Of the ~10,000-40,000 BHs that accumulate in the central parsec over the age of the Galaxy, we predict that ~60 - 200 currently exist as BH-XRBs formed from tidal capture, consistent with the population seen by Hailey et al. A somewhat lower number of tidal capture NS-XRBs is also predicted. We also use our observationally calibrated models for the NSC to predict rates of other exotic dynamical processes, such as the tidal disruption of stars by the central supermassive black hole (~0.0001 per year at z=0).
△ Less
Submitted 27 June, 2018; v1 submitted 4 April, 2018;
originally announced April 2018.
-
Active Galactic Nuclei Feedback in an Elliptical Galaxy with the Most Updated AGN Physics (II): High-Angular Momentum Case
Authors:
Doosoo Yoon,
Feng Yuan,
Zhao-Ming Gan,
Jeremiah P. Ostriker,
Ya-Ping Li,
Luca Ciotti
Abstract:
This is the second paper of our series of works of studying the effects of active galactic nuclei (AGN) feedback on the cosmological evolution of an isolated elliptical galaxy by performing two-dimensional high-resolution hydrodynamical numerical simulations. In these simulations, the inner boundary is chosen so that the Bondi radius is resolved. Physical processes like star formation, SNe Ia and…
▽ More
This is the second paper of our series of works of studying the effects of active galactic nuclei (AGN) feedback on the cosmological evolution of an isolated elliptical galaxy by performing two-dimensional high-resolution hydrodynamical numerical simulations. In these simulations, the inner boundary is chosen so that the Bondi radius is resolved. Physical processes like star formation, SNe Ia and II are taken into account. Compared to previous works, the main improvements is that we adopt the most updated AGN physics, which is described in detail in the first paper of this series (Yuan et al. 2018, Paper I). These improvements include the discrimination of the two accretion modes of the central AGN and the most updated descriptions of the wind and radiation in the two modes. In Paper I, we consider the case that the specific angular momentum of the gas in the galaxy is very low. In this paper, we consider the case that the specific angular momentum of the gas is high. In the galactic scale, we adopt the gravitational torques raised due to non-axisymmetric structure in the galaxy as the mechanism of the transfer of angular momentum of gas, as proposed in some recent works. Since our simulations are axisymmetric, we make use of a parameterized prescription to mimic this mechanism. Same as Paper I, we investigate the AGN light curve, typical AGN lifetime, growth of the black hole mass, AGN duty-cycle, star formation, and the X-ray surface brightness of the galaxy. Special attention is paid to the effects of specific angular momentum of the galaxy on these properties. We find that some results are qualitatively similar to those shown in Paper I, while some results such as star formation and black hole growth do show a significant difference due to the mass concentration in the galactic disk as a consequence of galactic rotation.
△ Less
Submitted 17 July, 2018; v1 submitted 9 March, 2018;
originally announced March 2018.
-
Stellar and AGN feedback in isolated early-type galaxies: the role in regulating star formation and ISM properties
Authors:
Ya-Ping Li,
Feng Yuan,
Houjun Mo,
Doosoo Yoon,
Zhao-Ming Gan,
Luis C. Ho,
Bo Wang,
Jeremiah P. Ostriker,
Luca Ciotti
Abstract:
Understanding how galaxies maintain the inefficiency of star formation with physically self-consistent models is a central problem for galaxy evolution. Although numerous theoretical models have been proposed in recent decades, the debate still exists. By means of high-resolution two-dimensional hydrodynamical simulations, we study the three feedback effects (the stellar wind heating, SNe feedback…
▽ More
Understanding how galaxies maintain the inefficiency of star formation with physically self-consistent models is a central problem for galaxy evolution. Although numerous theoretical models have been proposed in recent decades, the debate still exists. By means of high-resolution two-dimensional hydrodynamical simulations, we study the three feedback effects (the stellar wind heating, SNe feedback, and AGN feedback) in suppressing star formation activities on the evolution of early-type galaxies with different stellar masses. AGN feedback models are updated based on \citet{Yuan2018}. The gas sources comes exclusively from the mass losses of dying low-mass stars for most of our models. We find that SNe feedback can keep star formation at a significantly low level for low mass elliptical galaxies for a cosmological evolution time. For the high mass galaxies, AGN feedback can efficiently offset the radiative cooling and thus regulate the star formation activities. Such a suppression of star formation is extremely efficient in the inner region of the galaxies. AGB heating cannot account for this suppression for low and high mass galaxies. The X-ray temperature $T_{\rm X}$ and luminosity $L_{\rm X}$ of hot plasma can be in agreement with the observed data with the inclusion of effective feedback processes. These results thus suggest that we can use $T_{\rm X}$ and $L_{\rm X}$ to probe the role of different feedback processes. The inclusion of additional gas sources can make the mass scale between SNe and AGN feedback dominating in suppressing star formation decrease to an observationally inferred value of a few $10^{10}~M_{\odot}$.
△ Less
Submitted 3 September, 2018; v1 submitted 4 March, 2018;
originally announced March 2018.
-
AGN feedback and the origin and fate of the hot gas in early-type galaxies
Authors:
S. Pellegrini,
L. Ciotti,
A. Negri,
J. P. Ostriker
Abstract:
A recent determination of the relationships between the X-ray luminosity of the ISM (Lx) and the stellar and total mass, for a sample of nearby early-type galaxies (ETGs), is used to investigate the origin of the hot gas, via a comparison with the results of hydrodynamical simulations of the ISM evolution for a large set of isolated ETGs. After the epoch of major galaxy formation (after z~2), the…
▽ More
A recent determination of the relationships between the X-ray luminosity of the ISM (Lx) and the stellar and total mass, for a sample of nearby early-type galaxies (ETGs), is used to investigate the origin of the hot gas, via a comparison with the results of hydrodynamical simulations of the ISM evolution for a large set of isolated ETGs. After the epoch of major galaxy formation (after z~2), the ISM is replenished by stellar mass losses and SN ejecta, at the rate predicted by stellar evolution, and is depleted by star formation; it is heated by the thermalization of stellar motions, SNe explosions and the mechanical (from winds) and radiative AGN feedback. The models agree well with the observed relations, even for the largely different Lx values at the same mass, thanks to the sensitivity of the gas flow to many galaxy properties; this holds for models including AGN feedback, and those without. Therefore, the mass input from the stellar population is able to account for a major part of the observed Lx; and AGN feedback, while very important to maintain massive ETGs in a time-averaged quasi-steady state, keeping low star formation and the black hole mass, does not dramatically alter the gas content originating in stellar recycled material. These conclusions are based on theoretical predictions for the stellar population contributions in mass and energy, and on a self-consistent modeling of AGN feedback.
△ Less
Submitted 6 February, 2018;
originally announced February 2018.
-
Low density, radiatively inefficient rotating-accretion flow onto a black hole
Authors:
Kohei Inayoshi,
Jeremiah P. Ostriker,
Zoltan Haiman,
Rolf Kuiper
Abstract:
We study low-density axisymmetric accretion flows onto black holes (BHs) with two-dimensional hydrodynamical simulations, adopting the $α$-viscosity prescription. When the gas angular momentum is low enough to form a rotationally supported disk within the Bondi radius ($R_{\rm B}$), we find a global steady accretion solution. The solution consists of a rotational equilibrium distribution at…
▽ More
We study low-density axisymmetric accretion flows onto black holes (BHs) with two-dimensional hydrodynamical simulations, adopting the $α$-viscosity prescription. When the gas angular momentum is low enough to form a rotationally supported disk within the Bondi radius ($R_{\rm B}$), we find a global steady accretion solution. The solution consists of a rotational equilibrium distribution at $r\sim R_{\rm B}$, where the density follows $ρ\propto (1+R_{\rm B}/r)^{3/2}$, surrounding a geometrically thick and optically thin accretion disk at the centrifugal radius, where thermal energy generated by viscosity is transported via strong convection. Physical properties of the inner solution agree with those expected in convection-dominated accretion flows (CDAF; $ρ\propto r^{-1/2}$). In the inner CDAF solution, the gas inflow rate decreases towards the center due to convection ($\dot{M}\propto r$), and the net accretion rate (including both inflows and outflows) is strongly suppressed by several orders of magnitude from the Bondi accretion rate $\dot{M}_{\rm B}$ The net accretion rate depends on the viscous strength, following $\dot{M}/\dot{M}_{\rm B}\propto (α/0.01)^{0.6}$. This solution holds for low accretion rates of $\dot{M}_{\rm B}/\dot{M}_{\rm Edd}< 10^{-3}$ having minimal radiation cooling, where $\dot{M}_{\rm Edd}$ is the Eddington rate. In a hot plasma at the bottom ($r<10^{-3}~R_{\rm B}$), thermal conduction would dominate the convective energy flux. Since suppression of the accretion by convection ceases, the final BH feeding rate is found to be $\dot{M}/\dot{M}_{\rm B} \sim 10^{-3}-10^{-2}$. This rate is as low as $\dot{M}/\dot{M}_{\rm Edd} \sim 10^{-7}-10^{-6}$ inferred for SgrA$^*$ and the nuclear BHs in M31 and M87, and can explain the low luminosities in these sources, without invoking any feedback mechanism.
△ Less
Submitted 17 September, 2018; v1 submitted 21 September, 2017;
originally announced September 2017.
-
Interactions between multiple supermassive black holes in galactic nuclei: a solution to the final parsec problem
Authors:
Taeho Ryu,
Rosalba Perna,
Zoltán Haiman,
Jeremiah P. Ostriker,
Nicholas C. Stone
Abstract:
Using few-body simulations, we investigate the evolution of supermassive black holes (SMBHs) in galaxies ($M_{\star}=10^{10}-10^{12}{\rm M}_{\odot}$ at $z=0$) at $0<z<4$. Following galaxy merger trees from the Millennium simulation, we model BH mergers with two extreme binary decay scenarios for the `hard binary' stage: a full or an empty loss cone. These two models should bracket the true evoluti…
▽ More
Using few-body simulations, we investigate the evolution of supermassive black holes (SMBHs) in galaxies ($M_{\star}=10^{10}-10^{12}{\rm M}_{\odot}$ at $z=0$) at $0<z<4$. Following galaxy merger trees from the Millennium simulation, we model BH mergers with two extreme binary decay scenarios for the `hard binary' stage: a full or an empty loss cone. These two models should bracket the true evolution, and allow us to separately explore the role of dynamical friction and that of multi-body BH interactions on BH mergers. Using the computed merger rates, we infer the stochastic gravitational wave background (GWB). Our dynamical approach is a first attempt to study the dynamical evolution of multiple SMBHs in the host galaxies undergoing mergers with various mass ratios ($10^{-4} < q_{\star} < 1$). Our main result demonstrates that SMBH binaries are able to merge in both scenarios. In the empty loss cone case, we find that BHs merge via multi-body interactions, avoiding the `final parsec' problem, and entering the PTA band with substantial orbital eccentricity. Our full loss cone treatment, albeit more approximate, suggests that the eccentricity becomes even higher when GWs become dominant, leading to rapid coalescences (binary lifetime $\lesssim1 {\rm ~Gyr}$). Despite the lower merger rates in the empty loss cone case, due to their higher mass ratios and lower redshifts, the GWB in the full/empty loss cone models are comparable ($0.70\times10^{-15}$ and $0.53\times10^{-15}$ at a frequency of $1~{\rm yr}^{-1}$, respectively). Finally, we compute the effects of high eccentricities on the GWB spectrum.
△ Less
Submitted 25 September, 2017; v1 submitted 19 September, 2017;
originally announced September 2017.
-
Synthetic nebular emission from massive galaxies I: origin of the cosmic evolution of optical emission-line ratios
Authors:
Michaela Hirschmann,
Stephane Charlot,
Anna Feltre,
Thorsten Naab,
Ena Choi,
Jeremiah P. Ostriker,
Rachel S. Somerville
Abstract:
Galaxies occupy different regions of the [OIII]$\lambda5007$/H$β$-versus-[NII]$\lambda6584$/H$α$ emission-line ratio diagram in the distant and local Universe. We investigate the origin of this intriguing result by modelling self-consistently, for the first time, nebular emission from young stars, accreting black holes (BHs) and older, post-asymptotic-giant-branch (post-AGB) stellar populations in…
▽ More
Galaxies occupy different regions of the [OIII]$\lambda5007$/H$β$-versus-[NII]$\lambda6584$/H$α$ emission-line ratio diagram in the distant and local Universe. We investigate the origin of this intriguing result by modelling self-consistently, for the first time, nebular emission from young stars, accreting black holes (BHs) and older, post-asymptotic-giant-branch (post-AGB) stellar populations in galaxy formation simulations in a full cosmological context. In post-processing, we couple new-generation nebular-emission models with high-resolution, cosmological zoom-in simulations of massive galaxies to explore which galaxy physical properties drive the cosmic evolution of the optical-line ratios [OIII]$\lambda5007$/H$β$, [NII]$\lambda6584$/H$α$, [SII]$λ\lambda6717,6731$/H$α$ and [OI]$\lambda6300$/H$α$. The line ratios of simulated galaxies agree well with observations of both star-forming and active local SDSS galaxies. Towards higher redshifts, at fixed galaxy stellar mass, the average [OIII]/H$β$ increases and [NII]/H$α$, [SII]/H$α$ and [OI]/H$α$ decrease -- widely consistent with observations. At fixed stellar mass, we identify star formation history, which controls nebular emission from young stars via the ionization parameter, as the primary driver of the cosmic evolution of [OIII]/H$β$ and [NII]/H$α$. For [SII]/H$α$ and [OI]/H$α$, this applies only to redshifts above $z=1.5$, the evolution at lower redshift being driven in roughly equal parts by nebular emission from AGN and post-AGB stars. Instead, changes in the hardness of ionizing radiation, ionized-gas density, the prevalence of BH accretion relative to star formation and the dust-to-metal mass ratio (whose impact on the gas-phase N/O ratio we model at fixed O/H) play at most a minor role in the cosmic evolution of simulated galaxy line ratios.
△ Less
Submitted 31 August, 2017; v1 submitted 31 May, 2017;
originally announced June 2017.
-
AGN feedback, quiescence and CGM metal enrichment in early-type galaxies
Authors:
Maximilian Eisenreich,
Thorsten Naab,
Ena Choi,
Jeremiah P. Ostriker,
Eric Emsellem
Abstract:
We present three-dimensional hydrodynamical simulations showing the effect of kinetic and radiative AGN feedback on a model galaxy representing a massive quiescent low-redshift early-type galaxy of $M_* = 8.41\times 10^{10} M_\odot$, harbouring a $M_\mathrm{BH} = 4\times 10^8 M_\odot $ black hole surrounded by a cooling gaseous halo. We show that, for a total baryon fraction of $\sim 20\%$ of the…
▽ More
We present three-dimensional hydrodynamical simulations showing the effect of kinetic and radiative AGN feedback on a model galaxy representing a massive quiescent low-redshift early-type galaxy of $M_* = 8.41\times 10^{10} M_\odot$, harbouring a $M_\mathrm{BH} = 4\times 10^8 M_\odot $ black hole surrounded by a cooling gaseous halo. We show that, for a total baryon fraction of $\sim 20\%$ of the cosmological value, feedback from the AGN can keep the galaxy quiescent for about 4.35 Gyr and with properties consistent with black hole mass and X-ray luminosity scaling relations. However, this can only be achieved if the AGN feedback model includes both kinetic and radiative feedback modes. The simulation with only kinetic feedback fails to keep the model galaxy fully quiescent, while one with only radiative feedback leads to excessive black-hole growth. For higher baryon fractions (e.g. 50\% of the cosmological value), the X-ray luminosities exceed observed values by at least one order of magnitude, and rapid cooling results in a star-forming galaxy. The AGN plays a major role in keeping the circumgalactic gas at observed metallicities of $Z/Z_\odot \gtrsim 0.3 $ within the central $\sim 30$ kpc by venting nuclear gas enriched with metals from residual star formation activity. As indicated by previous cosmological simulations, our results are consistent with a model for which the black hole mass and the total baryon fraction are set at higher redshifts $z > 1$ and the AGN alone can keep the model galaxy on observed scaling relations. Models without AGN feedback violate both the quiescence criterion as well as CGM metallicity constraints.
△ Less
Submitted 7 April, 2017; v1 submitted 22 February, 2017;
originally announced February 2017.
-
Modeling for Stellar Feedback in Galaxy Formation Simulations
Authors:
Alejandro Núñez,
Jeremiah P. Ostriker,
Thorsten Naab,
Ludwig Oser,
Chia-Yu Hu,
Ena Choi
Abstract:
Various heuristic approaches to model unresolved supernova (SN) feedback in galaxy formation simulations exist to reproduce the formation of spiral galaxies and the overall inefficient conversion of gas into stars. Some models, however, require resolution dependent scalings. We present a sub-resolution model representing the three major phases of supernova blast wave evolution $-$free expansion, e…
▽ More
Various heuristic approaches to model unresolved supernova (SN) feedback in galaxy formation simulations exist to reproduce the formation of spiral galaxies and the overall inefficient conversion of gas into stars. Some models, however, require resolution dependent scalings. We present a sub-resolution model representing the three major phases of supernova blast wave evolution $-$free expansion, energy conserving Sedov-Taylor, and momentum conserving snowplow$-$ with energy scalings adopted from high-resolution interstellar-medium simulations in both uniform and multiphase media. We allow for the effects of significantly enhanced SN remnant propagation in a multiphase medium with the cooling radius scaling with the hot volume fraction, $f_{\mathrm{hot}}$, as $(1 - f_{\mathrm{hot}})^{-4/5}$. We also include winds from young massive stars and AGB stars, Strömgren sphere gas heating by massive stars, and a gas cooling limiting mechanism driven by radiative recombination of dense HII regions. We present initial tests for isolated Milky-Way like systems simulated with the GADGET based code SPHgal with improved SPH prescription. Compared to pure thermal SN input, the model significantly suppresses star formation at early epochs, with star formation extended both in time and space in better accord with observations. Compared to models with pure thermal SN feedback, the age at which half the stellar mass is assembled increases by a factor of 2.4, and the mass loading parameter and gas outflow rate from the galactic disk increase by a factor of 2. Simulation results are converged for a two order of magnitude variation in particle mass in the range (1.3$-$130)$\times 10^4$ solar masses.
△ Less
Submitted 4 January, 2017;
originally announced January 2017.
-
Theoretical Challenges in Galaxy Formation
Authors:
Thorsten Naab,
Jeremiah P. Ostriker
Abstract:
Numerical simulations have become a major tool for understanding galaxy formation and evolution. Over the decades the field has made significant progress. It is now possible to simulate the formation of individual galaxies and galaxy populations from well defined initial conditions with realistic abundances and global properties. An essential component of the calculation is to correctly estimate t…
▽ More
Numerical simulations have become a major tool for understanding galaxy formation and evolution. Over the decades the field has made significant progress. It is now possible to simulate the formation of individual galaxies and galaxy populations from well defined initial conditions with realistic abundances and global properties. An essential component of the calculation is to correctly estimate the inflow to and outflow from forming galaxies since observations indicating low formation efficiency and strong circum-glactic presence of gas are persuasive. Energetic 'feedback' from massive stars and accreting super-massive black holes - generally unresolved in cosmological simulations - plays a major role for driving galactic outflows, which have been shown to regulate many aspects of galaxy evolution. A surprisingly large variety of plausible sub-resolution models succeeds in this exercise. They capture the essential characteristics of the problem, i.e. outflows regulating galactic gas flows, but their predictive power is limited. In this review we focus on one major challenge for galaxy formation theory: to understand the underlying physical processes that regulate the structure of the interstellar medium, star formation and the driving of galactic outflows. This requires accurate physical models and numerical simulations, which can precisely describe the multi-phase structure of the interstellar medium on the currently unresolved few hundred parsecs scales of large scale cosmological simulations. Such models ultimately require the full accounting for the dominant cooling and heating processes, the radiation and winds from massive stars and accreting black holes, an accurate treatment of supernova explosions as well as the non-thermal components of the interstellar medium like magnetic fields and cosmic rays.
△ Less
Submitted 20 December, 2016;
originally announced December 2016.
-
OVI Emission From the Supernovae-regulated Interstellar Medium: Simulation Vs Observation
Authors:
Miao Li,
Greg Bryan,
Jeremiah P. Ostriker
Abstract:
The OVI $λλ$1032, 1038Å doublet emission traces collisionally ionized gas with $T\approx 10^{5.5}$ K, where the cooling curve peaks for metal-enriched plasma. This warm-hot phase is usually not well-resolved in numerical simulations of the multiphase interstellar medium (ISM), but can be responsible for a significant fraction of the emitted energy. Comparing simulated OVI emission to observations…
▽ More
The OVI $λλ$1032, 1038Å doublet emission traces collisionally ionized gas with $T\approx 10^{5.5}$ K, where the cooling curve peaks for metal-enriched plasma. This warm-hot phase is usually not well-resolved in numerical simulations of the multiphase interstellar medium (ISM), but can be responsible for a significant fraction of the emitted energy. Comparing simulated OVI emission to observations is therefore a valuable test of whether simulations predict reasonable cooling rates from this phase. We calculate OVI $λ$1032Å emission, assuming collisional ionization equilibrium, for our small-box simulations of the stratified ISM regulated by supernovae. We find that the agreement is very good for our solar neighborhood model, both in terms of emission flux and mean OVI density seen in absorption. We explore runs with higher surface densities and find that, in our simulations, the OVI emission from the disk scales roughly linearly with the star formation rate. Observations of OVI emission are rare for external galaxies, but our results do not show obvious inconsistency with the existing data. Assuming the solar metallicity, OVI emission from the galaxy disk in our simulations accounts for roughly 0.5\% of supernovae heating.
△ Less
Submitted 6 January, 2017; v1 submitted 1 November, 2016;
originally announced November 2016.
-
Towards Self-Consistent Modelling of the Sgr A* Accretion Flow: Linking Theory and Observation
Authors:
Shawn R. Roberts,
Yan-Fei Jiang,
Q. Daniel Wang,
Jeremiah P. Ostriker
Abstract:
The interplay between supermassive black holes (SMBHs) and their environments is believed to command an essential role in galaxy evolution. The majority of these SMBHs are in the radiative inefficient accretion phase where this interplay remains elusive, but suggestively important, due to few observational constraints. To remedy this, we directly fit 2-D hydrodynamic simulations to Chandra observa…
▽ More
The interplay between supermassive black holes (SMBHs) and their environments is believed to command an essential role in galaxy evolution. The majority of these SMBHs are in the radiative inefficient accretion phase where this interplay remains elusive, but suggestively important, due to few observational constraints. To remedy this, we directly fit 2-D hydrodynamic simulations to Chandra observations of Sgr A* with Markov Chain Monte Carlo sampling, self-consistently modelling the 2-D inflow-outflow solution for the first time. We find the temperature and density at flow onset are consistent with the origin of the gas in the stellar winds of massive stars in the vicinity of Sgr A*. We place the first observational constraints on the angular momentum of the gas and estimate the centrifugal radius, r$_c$ $\approx$ 0.056 r$_b$ $\approx8\times10^{-3}$ pc, where r$_b$ is the Bondi radius. Less than 1\% of the inflowing gas accretes onto the SMBH, the remainder being ejected in a polar outflow. For the first time, we decouple the quiescent point-like emission from the spatially extended flow. We find this point-like emission, accounting for $\sim4$\% of the quiescent flux, is spectrally too steep to be explained by unresolved flares, nor Bremsstrahlung, but is likely a combination of a relatively steep synchrotron powerlaw and the high energy tail of inverse-Compton emission. With this self-consistent model of the accretion flow structure, we make predictions for the flow dynamics and discuss how future X-ray spectroscopic observations can further our understanding of the Sgr A* accretion flow.
△ Less
Submitted 27 December, 2016; v1 submitted 31 October, 2016;
originally announced November 2016.
-
Physics of Galactic Metals: Evolutionary Effects due to Production, Distribution, Feedback & Interaction with Black Holes
Authors:
Ena Choi,
Jeremiah P. Ostriker,
Thorsten Naab,
Rachel S. Somerville,
Michaela Hirschmann,
Alejandro Núñez,
Chia-Yu Hu,
Ludwig Oser
Abstract:
We ask how the inclusion of various physical heating processes due to the metal content of gas affect the evolution of massive galaxies and compute a suite of cosmological hydrodynamical simulations that follow these systems and their supermassive black holes. We use a smoothed particle hydrodynamics code with a pressure-entropy formulation and a more accurate treatment of the metal production, tu…
▽ More
We ask how the inclusion of various physical heating processes due to the metal content of gas affect the evolution of massive galaxies and compute a suite of cosmological hydrodynamical simulations that follow these systems and their supermassive black holes. We use a smoothed particle hydrodynamics code with a pressure-entropy formulation and a more accurate treatment of the metal production, turbulent diffusion and cooling rate based on individual element abundances. The feedback models include (1) AGN feedback via high velocity BAL winds and Compton/photoionization heating, (2) explicit stellar feedback from multiple processes including powerful winds from supernova events, stellar winds from young massive stars and AGB stars as well as radiative heating within Stromgren spheres around massive stars, and (3) additional heating effects due to the presence of metals including grain photoelectric heating, metallicity dependent X-ray heating by nearby accreting black holes and from the cosmic X-ray background, which are the major improvement in our feedback model. With a suite of zoom-in simulations of 30 halos with $M_{vir} \sim 10^{12-13.4}$, we show that energy and momentum budget from all feedback effects generate realistic galaxy properties. We explore the detailed role of each feedback model with three additional sets of simulations with varying input physics. We show that the metal induced heating mechanisms reduce the fraction of accreted stellar material by mainly suppressing the growth of diffuse small stellar systems at high redshift but overall have a relatively minor effect on the final stellar and gas properties of massive galaxies. The inclusion of AGN feedback significantly improves the ability of our cosmological simulations to yield realistic gas and stellar properties of massive galaxies with reasonable fraction of the final stellar mass accreted from other galaxies.
△ Less
Submitted 14 June, 2017; v1 submitted 28 October, 2016;
originally announced October 2016.
-
Quantifying Supernovae-Driven Multiphase Galactic Outflows
Authors:
Miao Li,
Greg L. Bryan,
Jeremiah P. Ostriker
Abstract:
Galactic outflows are ubiquitously observed in star-forming disk galaxies and are critical for galaxy formation. Supernovae (SNe) play the key role in driving the outflows, but there is no consensus as to how much energy, mass and metal they can launch out of the disk. We perform 3D, high-resolution hydrodynamic simulations to study SNe-driven outflows from stratified media. Assuming SN rate scale…
▽ More
Galactic outflows are ubiquitously observed in star-forming disk galaxies and are critical for galaxy formation. Supernovae (SNe) play the key role in driving the outflows, but there is no consensus as to how much energy, mass and metal they can launch out of the disk. We perform 3D, high-resolution hydrodynamic simulations to study SNe-driven outflows from stratified media. Assuming SN rate scales with gas surface density $Σ_{\rm{gas}}$ as in the Kennicutt-Schmidt (KS) relation, we find the mass loading factor, $η_m$, defined as the mass outflow flux divided by the star formation surface density, decreases with increasing $Σ_{\rm{gas}}$ as $η_m \propto Σ^{-0.61}_{\rm{gas}}$. Approximately $Σ_{\rm{gas}}$ $\lesssim$ 50 $M_\odot/\rm{pc}^2$ marks when $η_m \gtrsim$1. About 10-50\% of the energy and 40-80\% of the metals produced by SNe end up in the outflows. The tenuous hot phase ($T>3\times 10^5$ K), which fills 60-80\% of the volume at mid-plane, carries the majority of the energy and metals in outflows. We discuss how various physical processes, including vertical distribution of SNe, photoelectric heating, external gravitational field and SN rate, affect the loading efficiencies. The relative scale height of gas and SNe is a very important factor in determining the loading efficiencies.
△ Less
Submitted 9 May, 2017; v1 submitted 27 October, 2016;
originally announced October 2016.
-
Ultralight scalars as cosmological dark matter
Authors:
Lam Hui,
Jeremiah P. Ostriker,
Scott Tremaine,
Edward Witten
Abstract:
An intriguing alternative to cold dark matter (CDM) is that the dark matter is a light ( $m \sim 10^{-22}$ eV) boson having a de Broglie wavelength $λ\sim 1$ kpc, often called fuzzy dark matter (FDM). We describe the arguments from particle physics that motivate FDM, review previous work on its astrophysical signatures, and analyze several unexplored aspects of its behavior. In particular, (i) FDM…
▽ More
An intriguing alternative to cold dark matter (CDM) is that the dark matter is a light ( $m \sim 10^{-22}$ eV) boson having a de Broglie wavelength $λ\sim 1$ kpc, often called fuzzy dark matter (FDM). We describe the arguments from particle physics that motivate FDM, review previous work on its astrophysical signatures, and analyze several unexplored aspects of its behavior. In particular, (i) FDM halos smaller than about $10^7 (m/10^{-22} {\rm eV})^{-3/2} M_\odot$ do not form. (ii) FDM halos are comprised of a core that is a stationary, minimum-energy configuration called a "soliton", surrounded by an envelope that resembles a CDM halo. (iii) The transition between soliton and envelope is determined by a relaxation process analogous to two-body relaxation in gravitating systems, which proceeds as if the halo were composed of particles with mass $\sim ρλ^3$ where $ρ$ is the halo density. (iv) Relaxation may have substantial effects on the stellar disk and bulge in the inner parts of disk galaxies. (v) Relaxation can produce FDM disks but an FDM disk in the solar neighborhood must have a half-thickness of at least $300 (m/10^{-22} {\rm eV})^{-2/3}$ pc. (vi) Solitonic FDM sub-halos evaporate by tunneling through the tidal radius and this limits the minimum sub-halo mass inside 30 kpc of the Milky Way to roughly $10^8 (m/10^{-22} {\rm eV})^{-3/2} M_\odot$. (vii) If the dark matter in the Fornax dwarf galaxy is composed of CDM, most of the globular clusters observed in that galaxy should have long ago spiraled to its center, and this problem is resolved if the dark matter is FDM.
△ Less
Submitted 24 February, 2017; v1 submitted 26 October, 2016;
originally announced October 2016.
-
The Effect of the AGN Feedback on the Interstellar Medium of Early-Type Galaxies: 2D Hydrodynamical Simulations of the Low-Rotation Case
Authors:
L. Ciotti,
S. Pellegrini,
A. Negri,
J. P. Ostriker
Abstract:
We present 2D hydrodynamical simulations for the evolution of early-type galaxies containing central massive black holes (MBHs), starting at age 2 Gyr. The code contains accurate and physically consistent radiative and mechanical AGN wind feedback, with parsec-scale central resolution. Mass input comes from stellar evolution; energy input includes Type Ia and II supernova and stellar heating; star…
▽ More
We present 2D hydrodynamical simulations for the evolution of early-type galaxies containing central massive black holes (MBHs), starting at age 2 Gyr. The code contains accurate and physically consistent radiative and mechanical AGN wind feedback, with parsec-scale central resolution. Mass input comes from stellar evolution; energy input includes Type Ia and II supernova and stellar heating; star-formation is included. Realistic, axisymmetric dynamical models for the galaxies are built solving the Jeans' equations. The lowest mass models (Mstar = 8 10^{10}Msun) develop global outflows sustained by SNIa's heating, ending with a significantly lower amount of hot gas and new stars. In more massive models, nuclear outbursts last to the present epoch, with large and frequent fluctuations in nuclear emission and from the gas (Lx). Each burst last ~ 10^{7.5} yr, during which (for r < 2-3 kpc) cold, inflowing, and hot, outflowing gas phases coexist. The Lx-T relation for the gas matches that of local galaxies. AGN activity causes positive feedback for star formation. Roughly half of the total mass loss is recycled into new stars (DeltaMstar), just ~ 3% of it is accreted on the MBH, the remainder being ejected from the galaxy. The ratio between the mass of gas expelled to that in to new stars, the load factor, is ~0.6. Rounder galaxies shapes lead to larger final MBH masses, DeltaMstar, and Lx. Almost all the time is spent at very low nuclear luminosities, yet one quarter of the total energy is emitted at an Eddington ratio > 0.1. The duty-cycle of AGN activity approximates 4% (Abridged).
△ Less
Submitted 11 August, 2016;
originally announced August 2016.
-
Formation of Massive Black Holes in Galactic Nuclei: Runaway Tidal Encounters
Authors:
Nicholas C. Stone,
Andreas H. W. Kuepper,
Jeremiah P. Ostriker
Abstract:
Nuclear star clusters (NSCs) and supermassive black holes (SMBHs) both inhabit galactic nuclei, coexisting in a range of bulge masses, but excluding each other in the largest or smallest galaxies. We propose that the transformation of NSCs into SMBHs occurs via runaway tidal captures, once NSCs exceed a certain critical central density and velocity dispersion. The bottleneck in this process, as wi…
▽ More
Nuclear star clusters (NSCs) and supermassive black holes (SMBHs) both inhabit galactic nuclei, coexisting in a range of bulge masses, but excluding each other in the largest or smallest galaxies. We propose that the transformation of NSCs into SMBHs occurs via runaway tidal captures, once NSCs exceed a certain critical central density and velocity dispersion. The bottleneck in this process, as with all collisional runaways, is growing the first e-fold in black hole mass. The growth of a stellar mass black hole past this bottleneck occurs as tidally captured stars are consumed in repeated episodes of mass transfer at pericenter. Tidal captures may turn off as a growth channel once the black hole reaches a mass ~100-1000 solar masses, but tidal disruption events will continue and appear capable of growing the seed SMBH to larger sizes. The runaway slows (becomes sub-exponential) once the seed SMBH consumes the core of its host NSC. While the bulk of the cosmic mass density in SMBHs is ultimately produced (via the Soltan-Paczynski argument) by episodic gaseous accretion in very massive galaxies, the smallest SMBHs have probably grown from strong tidal encounters with NSC stars. SMBH seeds that grow for a time $t$ entirely through this channel will follow simple power law relations with the velocity dispersion, $σ$, of their host galaxy. In the simplest regime it is $M_\bullet \sim σ^{3/2}\sqrt{M_\star t / G} \sim 10^{6}M_\odot (σ/ 50~{\rm km~s}^{-1})^{3/2}(t/10^{10}~{\rm yr})^{1/2}$, but the exponents and prefactor can differ slightly depending on the details of loss cone refilling. Current tidal disruption event rates predicted from this mechanism are consistent with observations.
△ Less
Submitted 6 June, 2016;
originally announced June 2016.
-
Hyper-Eddington accretion flows onto massive black holes
Authors:
Kohei Inayoshi,
Zoltan Haiman,
Jeremiah P. Ostriker
Abstract:
We study very-high rate spherically symmetric accretion flows onto a massive black hole (BH; 10^2 < M_BH < 10^6 Msun) embedded in a dense gas cloud with a low abundance of metals, performing one-dimensional hydrodynamical simulations which include multi-frequency radiation transfer and non-equilibrium primordial chemistry. We find that rapid gas supply from the Bondi radius at a hyper-Eddington ra…
▽ More
We study very-high rate spherically symmetric accretion flows onto a massive black hole (BH; 10^2 < M_BH < 10^6 Msun) embedded in a dense gas cloud with a low abundance of metals, performing one-dimensional hydrodynamical simulations which include multi-frequency radiation transfer and non-equilibrium primordial chemistry. We find that rapid gas supply from the Bondi radius at a hyper-Eddington rate can occur without being impeded by radiation feedback when (n/10^5 cm^-3) > (M_BH/10^4Msun)^{-1}(T/10^4 K)^{3/2}, where n and T are the density and temperature of ambient gas outside of the Bondi radius. The resulting accretion rate in this regime is steady, and larger than 3000 times the Eddington rate. At lower Bondi rates, the accretion is episodic due to radiative feedback and the average rate is limited below the Eddington rate. For the hyper-Eddington case, the steady solution consists of two parts: a radiation-dominated central core, where photon trapping due to electron scattering is important, and an accreting envelope which follows a Bondi profile with T~8000 K. When the emergent luminosity is limited below the Eddington luminosity because of photon trapping, radiation from the central region does not affect the gas dynamics at larger scales. We apply our result to the rapid formation of massive BHs in protogalaxies with a virial temperature of T_vir> 10^4 K. Once a seed BH forms at the center of the galaxy, it can grow up to a maximum ~10^5 (T_vir/10^4 K) Msun via gas accretion independent of the initial BH mass. Finally, we discuss possible observational signatures of rapidly accreting BHs with/without allowance for dust. We suggest that these systems could explain Lya emitters without X-rays and luminous infrared sources with hot dust emission, respectively.
△ Less
Submitted 7 April, 2016; v1 submitted 3 November, 2015;
originally announced November 2015.
-
AGN feedback and star formation in ETGs: negative and positive feedback
Authors:
L. Ciotti,
J. P. Ostriker,
A. Negri,
S. Pellegrini,
S. Posacki,
G. Novak
Abstract:
AGN feedback from supermassive black holes (SMBHs) at the center of early type galaxies is commonly invoked as the explanation for the quenching of star formation in these systems. The situation is complicated by the significant amount of mass injected in the galaxy by the evolving stellar population over cosmological times. In absence of feedback, this mass would lead to unobserved galactic cooli…
▽ More
AGN feedback from supermassive black holes (SMBHs) at the center of early type galaxies is commonly invoked as the explanation for the quenching of star formation in these systems. The situation is complicated by the significant amount of mass injected in the galaxy by the evolving stellar population over cosmological times. In absence of feedback, this mass would lead to unobserved galactic cooling flows, and to SMBHs two orders of magnitude more massive than observed. By using high-resolution 2D hydrodynamical simulations with radiative transport and star formation in state-of-the-art galaxy models, we show how the intermittent AGN feedback is highly structured on spatial and temporal scales, and how its effects are not only negative (shutting down the recurrent cooling episodes of the ISM), but also positive, inducing star formation in the inner regions of the host galaxy.
△ Less
Submitted 3 November, 2015; v1 submitted 2 November, 2015;
originally announced November 2015.