skip to main content


Title: Clump survival and migration in VDI galaxies: an analytical model versus simulations and observations
ABSTRACT

We address the nature of the giant clumps in high-z galaxies that undergo violent disc instability, distinguishing between long-lived and short-lived clumps. We study the evolution of long-lived clumps during migration through the disc via an analytical model tested by simulations and confront theory with CANDELS-HST observations. The clump ‘bathtub’ model, which considers gas and stellar gain and loss, involves four parameters: the accretion efficiency α, the star formation rate (SFR) efficiency ϵd, and the outflow mass-loading factors for gas and stars, η and ηs. The corresponding time-scales are comparable to the migration time, two-three orbital times. The accretion-rate dependence on clump mass, gas, and stars, allows an analytical solution involving exponential growing and decaying modes. For the fiducial parameter values there is a main evolution phase where the SFR and gas mass are constant and the stellar mass is rising linearly with time. This makes the inverse specific SFR an observable proxy for clump age. When η or ϵd are high, or α is low, the decaying mode induces a decline of SFR and gas mass till the migration ends. Later, the masses and SFR approach an hypothetical exponential growth with a constant specific SFR. The model matches simulations with different, moderate feedbacks, both in isolated and cosmological settings. The observed clumps agree with our predictions, indicating that the massive clumps are long-lived and migrating. A challenge is to model feedback that is non-disruptive in massive clumps but suppresses SFR to match the galactic stellar-to-halo mass ratio.

 
more » « less
NSF-PAR ID:
10362410
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
511
Issue:
1
ISSN:
0035-8711
Page Range / eLocation ID:
p. 316-340
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT We study the spatially resolved (sub-kpc) gas velocity dispersion (σ)–star formation rate (SFR) relation in the FIRE-2 (Feedback in Realistic Environments) cosmological simulations. We specifically focus on Milky Way-mass disc galaxies at late times (z ≈ 0). In agreement with observations, we find a relatively flat relationship, with σ ≈ 15–30 km s−1 in neutral gas across 3 dex in SFRs. We show that higher dense gas fractions (ratios of dense gas to neutral gas) and SFRs are correlated at constant σ. Similarly, lower gas fractions (ratios of gas to stellar mass) are correlated with higher σ at constant SFR. The limits of the σ–ΣSFR relation correspond to the onset of strong outflows. We see evidence of ‘on-off’ cycles of star formation in the simulations, corresponding to feedback injection time-scales of 10–100 Myr, where SFRs oscillate about equilibrium SFR predictions. Finally, SFRs and velocity dispersions in the simulations agree well with feedback-regulated and marginally stable gas disc (Toomre’s Q = 1) model predictions, and the simulation data effectively rule out models assuming that gas turns into stars at (low) constant efficiency (i.e. 1 per cent per free-fall time). And although the simulation data do not entirely exclude gas accretion/gravitationally powered turbulence as a driver of σ, it appears to be subdominant to stellar feedback in the simulated galaxy discs at z ≈ 0. 
    more » « less
  2. Abstract

    The Milky Way (MW) stellar disk has both a thin and a thick component. The thin disk is composed mostly of younger stars (≲8 Gyr) with a lower abundance ofα-elements, while the thick disk contains predominantly older stars (≳8–12 Gyr) with a higherαabundance, giving rise to anα-bimodality most prominent at intermediate metallicities. A proposed explanation for the bimodality is an episode of clumpy star formation, where high-αstars form in massive clumps that appear in the first few billion years of the MW’s evolution, while low-αstars form throughout the disk and over a longer time span. To better understand the evolution of clumps, we track them and their constituent stars in two clumpy MW simulations that reproduce theα-abundance bimodality, one with 10% and the other with 20% supernova feedback efficiency. We investigate the paths that these clumps take in the chemical space ([O/Fe]–[Fe/H]) as well as their mass, star formation rate (SFR), formation location, lifetime, and merger history. The clumps in the simulation with lower feedback last longer on average, with several lasting hundreds of millions of years. Some of the clumps do not reach high-α, but the ones that do on average have a higher SFR, longer lifetime, greater mass, and form closer to the Galactic center than the ones that do not. Most clumps that reach high-αmerge with others and eventually spiral into the Galactic center, but shed stars along the way to form most of the thick-disk component.

     
    more » « less
  3. ABSTRACT In this paper, we use Hubble Space Telescope/WFC3 observations of six galaxies from the DYnamics of Newly Assembled Massive Object (DYNAMO) survey, combined with stellar population modelling of the SED, to determine the stellar masses of DYNAMO clumps. The DYNAMO sample has been shown to have properties similar to z ≈ 1.5 turbulent, clumpy discs. DYNAMO sample clump masses offer a useful comparison for studies of z > 1 in that the galaxies have the same properties, yet the observational biases are significantly different. Using DYNAMO, we can more easily probe rest-frame near-IR wavelengths and also probe finer spatial scales. We find that the stellar mass of DYNAMO clumps is typically 107−108M⊙. We employ a technique that makes non-parametric corrections in removal of light from nearby clumps, and carries out a locally determined disc subtraction. The process of disc subtraction is the dominant effect, and can alter clump masses at the 0.3 dex level. Using these masses, we investigate the stellar mass function (MF) of clumps in DYNAMO galaxies. DYNAMO stellar MFs follow a declining power law with slope α ≈ −1.4, which is slightly shallower than, but similar to what is observed in z > 1 lensed galaxies. We compare DYNAMO clump masses to results of simulations. The masses and galactocentric position of clumps in DYNAMO galaxies are more similar to long-lived clumps in simulations. Similar to recent DYNAMO results on the stellar population gradients, these results are consistent with simulations that do not employ strong ‘early’ radiative feedback prescriptions. 
    more » « less
  4. ABSTRACT

    Stellar-mass binary black holes (BBHs) embedded in active galactic nucleus (AGN) discs are possible progenitors of black hole mergers detected in gravitational waves by LIGO/VIRGO. To better understand the hydrodynamical evolution of BBHs interacting with the disc gas, we perform a suite of high-resolution 2D simulations of binaries in local disc (shearing-box) models, considering various binary mass ratios, eccentricities and background disc properties. We use the γ-law equation of state and adopt a robust post-processing treatment to evaluate the mass accretion rate, torque and energy transfer rate on the binary to determine its long-term orbital evolution. We find that circular comparable-mass binaries contract, with an orbital decay rate of a few times the mass doubling rate. Eccentric binaries always experience eccentricity damping. Prograde binaries with higher eccentricities or smaller mass ratios generally have slower orbital decay rates, with some extreme cases exhibiting orbital expansion. The averaged binary mass accretion rate depends on the physical size of the accretor. The accretion flows are highly variable, and the dominant variability frequency is the apparent binary orbital frequency (in the rotating frame around the central massive BH) for circular binaries but gradually shifts to the radial epicyclic frequency as the binary eccentricity increases. Our findings demonstrate that the dynamics of BBHs embedded in AGN discs is quite different from that of isolated binaries in their own circumbinary discs. Furthermore, our results suggest that the hardening time-scales of the binaries are much shorter than their migration time-scales in the disc, for all reasonable binary and disc parameters.

     
    more » « less
  5. ABSTRACT

    Previous studies of fueling black holes in galactic nuclei have argued (on scales ${\sim}0.01{-}1000\,$pc) accretion is dynamical with inflow rates $\dot{M}\sim \eta \, M_{\rm gas}/t_{\rm dyn}$ in terms of gas mass Mgas, dynamical time tdyn, and some η. But these models generally neglected expulsion of gas by stellar feedback, or considered extremely high densities where expulsion is inefficient. Studies of star formation, however, have shown on sub-kpc scales the expulsion efficiency fwind = Mejected/Mtotal scales with the gravitational acceleration as $(1-f_{\rm wind})/f_{\rm wind}\sim \bar{a}_{\rm grav}/\langle \dot{p}/m_{\ast }\rangle \sim \Sigma _{\rm eff}/\Sigma _{\rm crit}$ where $\bar{a}_{\rm grav}\equiv G\, M_{\rm tot}(\lt r)/r^{2}$ and $\langle \dot{p}/m_{\ast }\rangle$ is the momentum injection rate from young stars. Adopting this as the simplest correction for stellar feedback, $\eta \rightarrow \eta \, (1-f_{\rm wind})$, we show this provides a more accurate description of simulations with stellar feedback at low densities. This has immediate consequences, predicting the slope and normalization of the MBH − σ and MBH − Mbulge relation, LAGN −SFR relations, and explanations for outliers in compact Es. Most strikingly, because star formation simulations show expulsion is efficient (fwind ∼ 1) below total-mass surface density $M_{\rm tot}/\pi \, r^{2}\lt \Sigma _{\rm crit}\sim 3\times 10^{9}\, \mathrm{M}_{\odot }\, {\rm kpc^{-2}}$ (where $\Sigma _{\rm crit}=\langle \dot{p}/m_{\ast }\rangle /(\pi \, G)$), BH mass is predicted to specifically trace host galaxy properties above a critical surface brightness Σcrit (B-band $\mu _{\rm B}^{\rm crit}\sim 19\, {\rm mag\, arcsec^{-2}}$). This naturally explains why BH masses preferentially reflect bulge properties or central surface densities (e.g. $\Sigma _{1\, {\rm kpc}}$), not ‘total’ galaxy properties.

     
    more » « less