We study the minimal ‘bathtub’ toy model as a tool for capturing key processes of galaxy evolution and identifying robust successes and challenges in reproducing high-z observations. The source and ...sink terms of the continuity equations for gas and stars are expressed in simple terms from first principles. The assumed dependence of star formation rate (SFR) on gas mass self-regulates the system into a unique asymptotic behaviour, which is approximated by an analytic quasi-steady-state (QSS) solution. We address the validity of the QSS at different epochs independent of earlier conditions. At high z, where the accretion is gaseous, the specific SFR (sSFR) is predicted to be sSFR ≃ (1 + z)/35/2 Gyr−1, slightly above the cosmological specific accretion rate, as observed at z = 3–8. The gas fraction is expected to decline slowly, and the observations constrain the SFR efficiency per dynamical time to ϵ ≃ 0.02. The stellar-to-virial mass ratio f
sv is predicted to be constant in time, and the observed value requires an outflow mass-loading factor η ≃ 1–3, depending on the penetration efficiency of gas into the galaxy. However, at z ∼ 2, where stars are also accreted through mergers, there is a conflict between model and observations. The model that maximizes the sSFR, with the outflows fully recycled, underestimates the sSFR by a factor of ∼3 and overestimates f
sv. With strong outflows, the model can match the observed f
sv but then it underestimates the sSFR by an order of magnitude. We discuss potential remedies including a bias due to the exclusion of quenched galaxies.
Compaction-driven black hole growth Lapiner, Sharon; Dekel, Avishai; Dubois, Yohan
Monthly notices of the Royal Astronomical Society,
07/2021, Volume:
505, Issue:
1
Journal Article
Peer reviewed
Open access
ABSTRACT
We study the interplay between galaxy evolution and central black hole (BH) growth using the NewHorizon cosmological simulation. BH growth is slow when the dark-matter halo is below a golden ...mass of $M_{\rm v}\sim 10^{12}\, \rm M_\odot$, and rapid above it. The early suppression is primarily due to gas removal by supernova (SN) feedback in the shallow potential well, predicting that BHs of ${\sim}10^5\, \rm M_\odot$ tend to lie below the linear relation with bulge mass. Rapid BH growth is allowed when the halo is massive enough to lock in the SN ejecta by its deep potential well and its heated circumgalactic medium (CGM). The onset of BH growth between these two zones is triggered by a wet-compaction event, caused, e.g. by mergers or counter-rotating streams. It brings gas that lost angular momentum into the inner-$1\, {\rm kpc}$ ‘blue nugget’ and causes major transitions in the galaxy structural, kinematic, and compositional properties, including the onset of star-formation quenching. The compaction events are confined to the golden mass by the same mechanisms of SN feedback and hot CGM. The onset of BH growth is associated with its sinkage to the centre due to the compaction-driven deepening of the potential well and the associated dynamical friction. The galaxy golden mass is thus imprinted as a threshold for rapid BH growth, allowing the AGN feedback to keep the CGM hot and maintain long-term quenching. AGN feedback is not causing the onset of quenching; they are both caused by a compaction event when the mass is between the SN and hot-CGM zones.
Abstract
We find, using cosmological simulations of galaxy clusters, that the hot X-ray emitting intracluster medium (ICM) enclosed within the outer accretion shock extends out to Rshock ∼ (2–3)Rvir, ...where Rvir is the standard virial radius of the halo. Using a simple analytic model for satellite galaxies in the cluster, we evaluate the effect of ram-pressure stripping on the gas in the inner discs and in the haloes at different distances from the cluster centre. We find that significant removal of star-forming disc gas occurs only at r ≲ 0.5Rvir, while gas removal from the satellite halo is more effective and can occur when the satellite is found between Rvir and Rshock. Removal of halo gas sets the stage for quenching of the star formation by starvation over 2–3 Gyr, prior to the satellite entry to the inner cluster halo. This scenario explains the presence of quenched galaxies, preferentially discs, at the outskirts of galaxy clusters, and the delayed quenching of satellites compared to central galaxies.
ABSTRACT
We study the effects of Kelvin–Helmholtz Instability (KHI) on the cold streams that feed massive haloes at high redshift, generalizing our earlier results to include the effects of radiative ...cooling and heating from a UV background, using analytic models and high resolution idealized simulations. We currently do not consider self-shielding, thermal conduction, or gravity. A key parameter in determining the fate of the streams is the ratio of the cooling time in the turbulent mixing layer which forms between the stream and the background following the onset of the instability, $t_{\rm cool,\, mix}$, to the time in which the mixing layer expands to the width of the stream in the non-radiative case, tshear. This can be converted into a critical stream radius, Rs, crit, such that $R_{\rm s}/R_{\rm s,crit}=t_{\rm shear}/t_{\rm cool,\, mix}$. If Rs < Rs, crit, the non-linear evolution proceeds similarly to the non-radiative case studied by Mandelker et al. If Rs > Rs,crit, which we find to almost always be the case for astrophysical cold streams, the stream is not disrupted by KHI. Rather, background mass cools and condenses on to the stream, and can increase the mass of cold gas by a factor of ∼3 within 10 stream sound crossing times. The mass entrainment induces thermal energy losses from the background and kinetic energy losses from the stream, which we model analytically. Roughly half of the dissipated energy is radiated away from gas with $T \lt 5\times 10^4\, {\rm K}$, suggesting much of it will be emitted in Ly α.
We use cosmological simulations to study a characteristic evolution pattern of high-redshift galaxies. Early, stream-fed, highly perturbed, gas-rich discs undergo phases of dissipative contraction ...into compact, star-forming systems (‘blue’ nuggets) at z ∼ 4–2. The peak of gas compaction marks the onset of central gas depletion and inside-out quenching into compact ellipticals (red nuggets) by z ∼ 2. These are sometimes surrounded by gas rings or grow extended dry stellar envelopes. The compaction occurs at a roughly constant specific star formation rate (SFR), and the quenching occurs at a constant stellar surface density within the inner kpc (Σ1). Massive galaxies quench earlier, faster, and at a higher Σ1 than lower mass galaxies, which compactify and attempt to quench more than once. This evolution pattern is consistent with the way galaxies populate the SFR-size–mass space, and with gradients and scatter across the main sequence. The compaction is triggered by an intense inflow episode, involving (mostly minor) mergers, counter-rotating streams or recycled gas, and is commonly associated with violent disc instability. The contraction is dissipative, with the inflow rate >SFR, and the maximum Σ1 anticorrelated with the initial spin parameter. The central quenching is triggered by the high SFR and stellar/supernova feedback (maybe also active galactic nucleus feedback) due to the high central gas density, while the central inflow weakens as the disc vanishes. Suppression of fresh gas supply by a hot halo allows the long-term maintenance of quenching once above a threshold halo mass, inducing the quenching downsizing.
We consider a simple gravitational heating mechanism for the long-term quenching of cooling flows and star formation in massive dark matter haloes hosting elliptical galaxies and clusters. We showed ...earlier that the virial shock heating in haloes ≥1012M⊙ triggers natural quenching in 1012–1013M⊙ haloes. Our present analytic estimates and simple simulations argue that the long-term quenching in haloes ≥Mmin∼ 7 × 1012M⊙ could be due to the gravitational energy of cosmological accretion delivered to the inner halo hot gas by cold flows via ram-pressure drag and local shocks. Mmin is obtained by comparing the gravitational power of infall into the potential well with the overall radiative cooling rate. The heating wins if the gas inner density cusp is not steeper than r−0.5 and if the masses in the cold and hot phases are comparable. The effect is stronger at higher redshifts, making the maintenance easier also at later times. Particular energy carriers into the halo core are cold gas clumps of ∼105–108M⊙. Clumps ≥105M⊙ penetrate to the inner halo with sufficient kinetic energy before they disintegrate, but they have to be ≤108M⊙ for the drag to do enough work in a Hubble time. Pressure-confined ∼104K clumps are stable against their own gravity and remain gaseous once below the Bonnor–Ebert mass ∼108M⊙. Such clumps are also immune to tidal disruption. Clumps in the desired mass range could emerge by thermal instability in the outer halo or in the filaments that feed it if the conductivity is not too high. Alternatively, such clumps may be embedded in dark matter subhaloes if the ionizing flux is ineffective, but they separate from their subhaloes by ram pressure before entering the inner halo. Heating by dynamical friction becomes dominant for massive satellites, which can contribute up to one-third of the total gravitational heating. We conclude that gravitational heating by cosmological accretion is a viable alternative to active galactic nucleus feedback as a long-term quenching mechanism.
Abstract
We address the origin of ultra-diffuse galaxies (UDGs), which have stellar masses typical of dwarf galaxies but effective radii of Milky Way-sized objects. Their formation mechanism, and ...whether they are failed L
⋆ galaxies or diffuse dwarfs, are challenging issues. Using zoom-in cosmological simulations from the Numerical Investigation of a Hundred Astrophysical Objects (NIHAO) project, we show that UDG analogues form naturally in dwarf-sized haloes due to episodes of gas outflows associated with star formation. The simulated UDGs live in isolated haloes of masses 1010–11 M⊙, have stellar masses of 107–8.5 M⊙, effective radii larger than 1 kpc and dark matter cores. They show a broad range of colours, an average Sérsic index of 0.83, a typical distribution of halo spin and concentration, and a non-negligible H i gas mass of 107 − 9 M⊙, which correlates with the extent of the galaxy. Gas availability is crucial to the internal processes which form UDGs: feedback-driven gas outflows, and subsequent dark matter and stellar expansion, are the key to reproduce faint, yet unusually extended, galaxies. This scenario implies that UDGs represent a dwarf population of low surface brightness galaxies and should exist in the field. The largest isolated UDGs should contain more H i gas than less extended dwarfs of similar M
⋆.
ABSTRACT
We study ultra-diffuse galaxies (UDGs) in zoom in cosmological simulations, seeking the origin of UDGs in the field versus galaxy groups. We find that while field UDGs arise from dwarfs in a ...characteristic mass range by multiple episodes of supernova feedback (Di Cintio et al.), group UDGs may also form by tidal puffing up and they become quiescent by ram-pressure stripping. The field and group UDGs share similar properties, independent of distance from the group centre. Their dark-matter haloes have ordinary spin parameters and centrally dominant dark-matter cores. Their stellar components tend to have a prolate shape with a Sérsic index n ∼ 1 but no significant rotation. Ram pressure removes the gas from the group UDGs when they are at pericentre, quenching star formation in them and making them redder. This generates a colour/star-formation-rate gradient with distance from the centre of the dense environment, as observed in clusters. We find that ∼20 per cent of the field UDGs that fall into a massive halo survive as satellite UDGs. In addition, normal field dwarfs on highly eccentric orbits can become UDGs near pericentre due to tidal puffing up, contributing about half of the group-UDG population. We interpret our findings using simple toy models, showing that gas stripping is mostly due to ram pressure rather than tides. We estimate that the energy deposited by tides in the bound component of a satellite over one orbit can cause significant puffing up provided that the orbit is sufficiently eccentric. We caution that while the simulations produce UDGs that match the observations, they under-produce the more compact dwarfs in the same mass range, possibly because of the high threshold for star formation or the strong feedback.
Any successful model of galaxy formation needs to explain the low rate of star formation in the small progenitors of today's galaxies. This inefficiency is necessary for reproducing the low ...stellar-to-virial mass fractions, suggested by current abundance matching models. A possible driver of this low efficiency is the radiation pressure exerted by ionizing photons from massive stars. The effect of radiation pressure in cosmological, zoom-in galaxy formation simulations is modelled as a non-thermal pressure that acts only in dense and optically thick star-forming regions. We also include photoionization and photoheating by massive stars. The full photoionization of hydrogen reduces the radiative cooling in the 104-4.5 K regime. The main effect of radiation pressure is to regulate and limit the high values of gas density and the amount of gas available for star formation. This maintains a low star formation rate of ∼1 M⊙ yr−1 in haloes with masses about 1011 M⊙ at z ≃ 3. Infrared trapping and photoionization/photoheating processes are secondary effects in this mass range. The galaxies residing in these low-mass haloes contain only ∼0.6 per cent of the total virial mass in stars, roughly consistent with abundance matching. Radiative feedback maintains an extended galaxy with a rising circular velocity profile.
We study the scaling relations between global properties of dwarf galaxies in the local group. In addition to quantifying the correlations between pairs of variables, we explore the ‘shape’ of the ...distribution of galaxies in log parameter space using standardized principal component analysis, the analysis is performed first in the 3D structural parameter space of stellar mass M*, internal velocity V and characteristic radius R* (or surface brightness μ*). It is then extended to a 4D space that includes a stellar population parameter such as metallicity Z or star formation rate . We find that the local group dwarfs basically define a one-parameter ‘fundamental line’ (FL), primarily driven by stellar mass, M*. A more detailed inspection reveals differences between the star formation properties of dwarf irregulars (dI's) and dwarf ellipticals (dE's), beyond the tendency of the latter to be more massive. In particular, the metallicities of dI's are typically lower by a factor of 3 at a given M* and they grow faster with increasing M*, showing a tighter FL in the 4D space for the dE's. The structural scaling relations of dI's resemble those of the more massive spirals, but the dI's have lower star formation rates for a given M* which also grow faster with increasing M*. On the other hand, the FL of the dE's departs from the fundamental plane of bigger ellipticals. While the one-parameter nature of the FL and the associated slopes of the scaling relations are consistent with the general predictions of supernova feedback from Dekel & Woo, the differences between the FL's of the dE's and the dI's remain a challenge and should serve as a guide for the secondary physical processes responsible for these two types.