Nuclear outflows driven by accreting massive black holes are one of the main feedback mechanisms invoked at high-z to reproduce the distinct separation between star-forming disk galaxies and ...quiescent spheroidal systems. Yet our knowledge of feedback at high-z remains limited by the lack of observations of the multiple gas phases in galaxy outflows. In this work, we use new deep, high spatial resolution ALMA CO(3-2) and archival Very Large Telescope/SINFONI H observations to study the molecular and ionized components of the active galactic nucleus (AGN)-driven outflow in zC400528, a massive main-sequence galaxy at z = 2.3 in the process of quenching. We detect a powerful molecular outflow that shows a positive velocity gradient before a turnover and extends for at least ∼10 kpc from the nuclear region, about three times the projected size of the ionized wind. The molecular gas in the outflow does not reach velocities high enough to escape the galaxy and is therefore expected to be reaccreted. Keeping in mind the various assumptions involved in the analysis, we find that the mass and energetics of the outflow are dominated by the molecular phase. The AGN-driven outflow in zC400528 is powerful enough to deplete the molecular gas reservoir on a timescale comparable to that needed to exhaust it by star formation. This suggests that the nuclear outflow is one of the main quenching engines at work in the observed suppression of the central star formation activity in zC400528.
Abstract
We investigate the relationship between star formation activity and outflow properties on kiloparsec scales in a sample of 28 star-forming galaxies at
z
∼ 2–2.6, using adaptive optics ...assisted integral field observations from SINFONI on the Very Large Telescope. The narrow and broad components of the H
α
emission are used to simultaneously determine the local star formation rate surface density (
), and the outflow velocity
and mass outflow rate
, respectively. We find clear evidence for faster outflows with larger mass loading factors at higher
. The outflow velocities scale as
∝
0.34±0.10
, which suggests that the outflows may be driven by a combination of mechanical energy released by supernova explosions and stellar winds, as well as radiation pressure acting on dust grains. The majority of the outflowing material does not have sufficient velocity to escape from the galaxy halos, but will likely be re-accreted and contribute to the chemical enrichment of the galaxies. In the highest
regions the outflow component contains an average of ∼45% of the H
α
flux, while in the lower
regions only ∼10% of the H
α
flux is associated with outflows. The mass loading factor,
η
=
/SFR, is positively correlated with
but is relatively low even at the highest
:
η
≲ 0.5 × (380 cm
−3
/
n
e
). This may be in tension with the
η
≳ 1 required by cosmological simulations, unless a significant fraction of the outflowing mass is in other gas phases and has sufficient velocity to escape the galaxy halos.
We report high-quality, H or CO rotation curves (RCs) to several Re for 41 large, massive, star-forming disk galaxies (SFGs) across the peak of cosmic galaxy evolution (z ∼ 0.67-2.45), taken with the ...ESO-VLT, the LBT and IRAM-NOEMA. Most RC41 SFGs have reflection-symmetric RCs plausibly described by equilibrium dynamics. We fit the major axis position-velocity cuts using beam-convolved forward modeling generated in three dimensions, with models that include a bulge and turbulent disk component embedded in a dark matter (DM) halo. We include priors for stellar and molecular gas masses, optical light effective radii and inclinations, and DM masses from abundance-matching scaling relations. Two-thirds or more of the z ≥ 1.2 SFGs are baryon dominated within a few Re of typically 5.5 kpc and have DM fractions less than maximal disks (median 〈 f DM ( R e ) 〉 = 0.12 ). At lower redshift (z < 1.2), that fraction is less than one-third. DM fractions correlate inversely with the baryonic angular momentum parameter, baryonic surface density, and bulge mass. Inferred low DM fractions cannot apply to the entire disk and halo but more plausibly reflect a flattened, or cored, inner DM density distribution. The typical central "DM deficit" in these cores relative to Navarro-Frenk-White (NFW) distributions is ∼30% of the bulge mass. The observations are consistent with rapid radial transport of baryons in the first-generation massive gas-rich halos forming globally gravitationally unstable disks and leading to efficient build-up of massive bulges and central black holes. A combination of heating due to dynamical friction and AGN feedback may drive DM out of the initial cusps.
Abstract
This paper provides an update of our previous scaling relations between galaxy-integrated molecular gas masses, stellar masses, and star formation rates (SFRs), in the framework of the star ...formation main sequence (MS), with the main goal of testing for possible systematic effects. For this purpose our new study combines three independent methods of determining molecular gas masses from CO line fluxes, far-infrared dust spectral energy distributions, and ∼1 mm dust photometry, in a large sample of 1444 star-forming galaxies between
z
= 0 and 4. The sample covers the stellar mass range log(
M
*
/
M
⊙
) = 9.0–11.8, and SFRs relative to that on the MS,
δ
MS = SFR/SFR(MS), from 10
−1.3
to 10
2.2
. Our most important finding is that all data sets, despite the different techniques and analysis methods used, follow the same scaling trends, once method-to-method zero-point offsets are minimized and uncertainties are properly taken into account. The molecular gas depletion time
t
depl
, defined as the ratio of molecular gas mass to SFR, scales as (1 +
z
)
−0.6
× (
δ
MS)
−0.44
and is only weakly dependent on stellar mass. The ratio of molecular to stellar mass
μ
gas
depends on (
1
+
z
)
2.5
×
(
δ
MS
)
0.52
×
(
M
*
)
−
0.36
, which tracks the evolution of the specific SFR. The redshift dependence of
μ
gas
requires a curvature term, as may the mass dependences of
t
depl
and
μ
gas
. We find no or only weak correlations of
t
depl
and
μ
gas
with optical size
R
or surface density once one removes the above scalings, but we caution that optical sizes may not be appropriate for the high gas and dust columns at high
z
.
We make publicly available a catalog of calibrated environmental measures for galaxies in the five 3D-Hubble Space Telescope (HST)/CANDELS deep fields. Leveraging the spectroscopic and grism ...redshifts from the 3D-HST survey, multiwavelength photometry from CANDELS, and wider field public data for edge corrections, we derive densities in fixed apertures to characterize the environment of galaxies brighter than mag in the redshift range . By linking observed galaxies to a mock sample, selected to reproduce the 3D-HST sample selection and redshift accuracy, each 3D-HST galaxy is assigned a probability density function of the host halo mass, and a probability that it is a central or a satellite galaxy. The same procedure is applied to a z = 0 sample selected from Sloan Digital Sky Survey. We compute the fraction of passive central and satellite galaxies as a function of stellar and halo mass, and redshift, and then derive the fraction of galaxies that were quenched by environment specific processes. Using the mock sample, we estimate that the timescale for satellite quenching is it is longer at lower stellar mass or lower redshift, but remarkably independent of halo mass. This indicates that, in the range of environments commonly found within the 3D-HST sample ( ), satellites are quenched by exhaustion of their gas reservoir in the absence of cosmological accretion. We find that the quenching times can be separated into a delay phase, during which satellite galaxies behave similarly to centrals at fixed stellar mass, and a phase where the star formation rate drops rapidly ( Gyr), as shown previously at z = 0. We conclude that this scenario requires satellite galaxies to retain a large reservoir of multi-phase gas upon accretion, even at high redshift, and that this gas sustains star formation for the long quenching times observed.
Abstract
We analyze H
α
or CO rotation curves extending out to several galaxy effective radii for 100 massive, large, star-forming disk galaxies (SFGs) across the peak of cosmic galaxy star formation ...(
z
∼ 0.6–2.5), more than doubling the previous sample presented by Genzel et al. and Price et al. The observations were taken with SINFONI and KMOS integral-field spectrographs at the ESO-Very Large Telescope, LUCI-LBT, NOEMA-IRAM, and Atacama Large Millimeter/submillimeter Array. We fit the major-axis kinematics with beam-convolved, forward models of turbulent rotating disks with bulges embedded in dark matter (DM) halos, including the effects of pressure support. The fraction of dark to total matter within the disk effective radius (
R
e
∼ 5 kpc),
f
DM
(
R
e
) =
V
2
DM
(
R
e
)/
V
2
circ
(
R
e
) decreases with redshift: at
z
∼ 1 (
z
∼ 2) the median DM fraction is 0.38 ± 0.23 (0.27 ± 0.18), and a third (half) of all galaxies are
maximal
disks with
f
DM
(
R
e
) < 0.28. DM fractions correlate inversely with the baryonic surface density, and the low DM fractions can be explained with a flattened, or cored, inner DM density distribution. At
z
∼ 2, there is ≈40% less DM mass on average within
R
e
compared to expected values based on cosmological stellar-mass–halo-mass relations. The DM deficit is more evident at high star formation rate surface densities (≳2.5
M
⊙
yr
−1
kpc
2
) and galaxies with massive bulges (≥10
10
M
⊙
). A combination of stellar or active galactic nucleus feedback, and/or heating due to dynamical friction, may drive the DM from cuspy into cored mass distributions, pointing to an efficient buildup of massive bulges and central black holes at
z
∼ 2 SFGs.
Abstract
We present a follow-up analysis examining the dynamics and structures of 41 massive, large star-forming galaxies at
z
∼ 0.67 − 2.45 using both ionized and molecular gas kinematics. We fit ...the galaxy dynamics with models consisting of a bulge, a thick, turbulent disk, and an NFW dark matter halo, using code that fully forward-models the kinematics, including all observational and instrumental effects. We explore the parameter space using Markov Chain Monte Carlo (MCMC) sampling, including priors based on stellar and gas masses and disk sizes. We fit the full sample using extracted 1D kinematic profiles. For a subset of 14 well-resolved galaxies, we also fit the 2D kinematics. The MCMC approach robustly confirms the results from least-squares fitting presented in Paper I: the sample galaxies tend to be baryon-rich on galactic scales (within one effective radius). The 1D and 2D MCMC results are also in good agreement for the subset, demonstrating that much of the galaxy dynamical information is captured along the major axis. The 2D kinematics are more affected by the presence of noncircular motions, which we illustrate by constructing a toy model with constant inflow for one galaxy that exhibits residual signatures consistent with radial motions. This analysis, together with results from Paper I and other studies, strengthens the finding that massive, star-forming galaxies at
z
∼ 1 − 2 are baryon-dominated on galactic scales, with lower dark matter fractions toward higher baryonic surface densities. Finally, we present details of the kinematic fitting code used in this analysis.
Abstract
We investigate the resolved kinematics of the molecular gas, as traced by the Atacama Large Millimeter/submillimeter Array in CO (2−1), of 25 cluster member galaxies across three different ...clusters at a redshift of
z
∼ 1.6. This is the first large-scale analysis of the molecular gas kinematics of cluster galaxies at this redshift. By separately estimating the rotation curve of the approaching and receding sides of each galaxy via kinematic modeling, we quantify the difference in total circular velocity to characterize the overall kinematic asymmetry of each galaxy. 3/14 of the galaxies in our sample that we are able to model have similar degrees of asymmetry as that observed in galaxies in the field at similar redshift based on observations of mainly ionized gas. However, this leaves 11/14 galaxies in our sample with significantly higher asymmetry, and some of these galaxies have degrees of asymmetry of up to ∼50 times higher than field galaxies observed at similar redshift. Some of these extreme cases also have one-sided tail-like morphology seen in the molecular gas, supporting a scenario of tidal and/or ram pressure interaction. Such stark differences in the kinematic asymmetry in clusters versus the field suggest the evolutionary influence of dense environments, established as being a major driver of galaxy evolution at low redshift, is also active in the high-redshift universe.
Abstract
We present deep observations of a
z
= 1.4 massive, star-forming galaxy (SFG) in molecular and ionized gas at comparable spatial resolution (CO 3–2, NOrthern Extended Millimeter Array ...(NOEMA); H
α
, Large Binocular Telescope (LBT)). The kinematic tracers agree well, indicating that both gas phases are subject to the same gravitational potential and physical processes affecting the gas dynamics. We combine the one-dimensional velocity and velocity dispersion profiles in CO and H
α
to forward-model the galaxy in a Bayesian framework, combining a thick exponential disk, a bulge, and a dark matter halo. We determine the dynamical support due to baryons and dark matter, and find a dark matter fraction within one effective radius of
f
DM
(
≤
R
e
)
=
0.18
−
0.04
+
0.06
. Our result strengthens the evidence for strong baryon-dominance on galactic scales of massive
z
∼ 1–3 SFGs recently found based on ionized gas kinematics alone.