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
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.
ABSTRACT
We report on our combined analysis of HST, VLT/MUSE, VLT/SINFONI, and ALMA observations of the local Seyfert 2 galaxy, NGC 5728 to investigate in detail the feeding and feedback of the ...active galactic nucleus (AGN). The data sets simultaneously probe the morphology, excitation, and kinematics of the stars, ionized gas, and molecular gas over a large range of spatial scales (10 pc to 10 kpc). NGC 5728 contains a large stellar bar that is driving gas along prominent dust lanes to the inner 1 kpc where the gas settles into a circumnuclear ring. The ring is strongly star forming and contains a substantial population of young stars as indicated by the lowered stellar velocity dispersion and gas excitation consistent with H ii regions. We model the kinematics of the ring using the velocity field of the CO (2–1) emission and stars and find it is consistent with a rotating disc. The outer regions of the disc, where the dust lanes meet the ring, show signatures of inflow at a rate of 1 M$\odot$ yr−1. Inside the ring, we observe three molecular gas components corresponding to the circular rotation of the outer ring, a warped disc, and the nuclear stellar bar. The AGN is driving an ionized gas outflow that reaches a radius of 250 pc with a mass outflow rate of 0.08 M$\odot$ yr−1 consistent with its luminosity and scaling relations from previous studies. While we observe distinct holes in CO emission which could be signs of molecular gas removal, we find that largely the AGN is not disrupting the structure of the circumnuclear region.
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.
We report new detections and limits from a NOEMA and ALMA CO(1-0) search for molecular outflows in 13 local galaxies with high far-infrared surface brightness, and combine these with local universe ...CO outflow results from the literature. The CO line ratios and spatial outflow structure of our targets provide some constraints on the conversion steps from observables to physical quantities such as molecular mass outflow rates. Where available, ratios between outflow emission in higher J CO transitions and in CO(1-0) are typically consistent with excitation
R
i
1
≲ 1. However, for IRAS 13120−5453,
R
31
= 2.10 ± 0.29 indicates optically thin CO in the outflow. Like much of the outflow literature, we use
α
CO(1 − 0)
= 0.8, and we present arguments for using
C
= 1 in deriving molecular mass outflow rates
Ṁ
out
=
C
M
out
v
out
/
R
out
. We compare the two main methods for molecular outflow detection: CO millimeter interferometry and
Herschel
OH-based spectroscopic outflow searches. For 26 sources studied with both methods, we find an 80% agreement in detecting
v
out
≳ 150 km s
−1
outflows, and non-matches can be plausibly ascribed to outflow geometry and signal-to-noise ratio. For a published sample of 12 bright ultraluminous infrared galaxies with detailed OH-based outflow modeling, CO outflows are detected in all but one. Outflow masses, velocities, and sizes for these 11 sources agree well between the two methods, and modest remaining differences may relate to the different but overlapping regions sampled by CO emission and OH absorption. Outflow properties correlate better with active galactic nucleus (AGN) luminosity and with bolometric luminosity than with far-infrared surface brightness. The most massive outflows are found for systems with current AGN activity, but significant outflows in nonAGN systems must relate to star formation or to AGN activity in the recent past. We report scaling relations for the increase of outflow mass, rate, momentum rate, and kinetic power with bolometric luminosity. Short flow times of ∼10
6
yr and some sources with resolved multiple outflow episodes support a role of intermittent driving, likely by AGNs.
I present a technique to remove the residual OH airglow emission from near-infrared spectra. Historically, the need to subtract out the strong and variable OH airglow emission lines from 1 to 2.5-μm ...spectra has imposed severe restrictions on observational strategy. For integral-field spectroscopy, where the field of view is limited, the standard technique is to observe blank sky frames at regular intervals. However, even this does not usually provide sufficient compensation if individual exposure times are longer than 2–3 min due to (1) changes in the absolute flux of the OH lines, (2) variations in flux among the individual OH lines and (3) effects of instrumental flexure which can lead to ‘P-Cygni’ type residuals. The data processing method presented here takes all of these effects into account and serendipitously also improves background subtraction between the OH lines. It allows one, in principle, to use sky frames taken hours or days previously so that observations can be performed in a quasi-stare mode. As a result, the observing efficiency (i.e. fraction of time spent on a source) at the telescope can be dramatically increased.
We present a census of ionized gas outflows in 599 normal galaxies at redshift 0.6 < z < 2.7, mostly based on integral field spectroscopy of Hα, N ii, and S ii line emission. The sample fairly ...homogeneously covers the main sequence of star-forming galaxies with masses 9.0 < log(M */M ⊙) < 11.7, and probes into the regimes of quiescent galaxies and starburst outliers. About one-third exhibits the high-velocity component indicative of outflows, roughly equally split into winds driven by star formation (SF) and active galactic nuclei (AGNs). The incidence of SF-driven winds correlates mainly with SF properties. These outflows have typical velocities of ∼450 km s−1, local electron densities of n e ∼ 380 cm−3, modest mass loading factors of ∼0.1–0.2 at all galaxy masses, and energetics compatible with momentum driving by young stellar populations. The SF-driven winds may escape from log(M */M ⊙) ≲ 10.3 galaxies, but substantial mass, momentum, and energy in hotter and colder outflow phases seem required to account for low galaxy formation efficiencies in the low-mass regime. Faster AGN-driven outflows (∼1000–2000 km s−1) are commonly detected above log(M */M ⊙) ∼ 10.7, in up to ∼75% of log(M */M ⊙) ≳ 11.2 galaxies. The incidence, strength, and velocity of AGN-driven winds strongly correlates with stellar mass and central concentration. Their outflowing ionized gas appears denser (n e ∼ 1000 cm−3), and possibly compressed and shock-excited. These winds have comparable mass loading factors as the SF-driven winds but carry ∼10 (∼50) times more momentum (energy). The results confirm our previous findings of high-duty-cycle, energy-driven outflows powered by AGN above the Schechter mass, which may contribute to SF quenching.
We present the completed KMOS3D survey, an integral field spectroscopic survey of 739 \(\mathrm{log}({M}_{\star }/{M}_{\odot })\gt 9\) galaxies at 0.6 < z < 2.7 using the K-band Multi Object ...Spectrograph (KMOS) at the Very Large Telescope. The KMOS3D survey provides a population-wide census of kinematics, star formation, outflows, and nebular gas conditions both on and off the star-forming galaxy main sequence through the spatially resolved and integrated properties of Hα, N ii, and S ii emission lines. We detect Hα emission for 91% of galaxies on the main sequence of star formation and 79% overall. The depth of the survey has allowed us to detect galaxies with star formation rates below 1 M ⊙ yr−1, as well as to resolve 81% of detected galaxies with ≥3 resolution elements along the kinematic major axis. The detection fraction of Hα is a strong function of both color and offset from the main sequence, with the detected and nondetected samples exhibiting different spectral energy distribution shapes. Comparison of Hα and UV+IR star formation rates reveal that dust attenuation corrections may be underestimated by 0.5 dex at the highest masses (\(\mathrm{log}({M}_{\star }/{M}_{\odot })\gt 10.5\)). We confirm our first year results of a high rotation-dominated fraction (monotonic velocity gradient and v rot/\({\sigma }_{0}\gt \sqrt{3.36}\)) of 77% for the full KMOS3D sample. The rotation-dominated fraction is a function of both stellar mass and redshift, with the strongest evolution measured over the redshift range of the survey for galaxies with \(\mathrm{log}({M}_{\star }/{M}_{\odot })\lt 10.5\). With this paper, we include a final data release of all 739 observed objects (http://www.mpe.mpg.de/ir/KMOS3D).
Using the cosmological baryonic accretion rate and normal star formation (SF) efficiencies, we present a very simple model for star-forming galaxies that accounts for the mass and redshift ...dependences of the star formation rate (SFR)-mass and Tully-Fisher (TF) relations from z ∼ 2 to the present. The time evolution follows from the fact that each modeled galaxy approaches a steady state where the SFR follows the (net) cold gas accretion rate. The key feature of the model is a halo mass floor M min 10 11 M below which accretion is quenched in order to simultaneously account for the observed slopes of the SFR-mass and TF relations. The same successes cannot be achieved via an SF threshold (or delay) nor by varying the SF efficiency or the feedback efficiency. Combined with the mass ceiling for cold accretion due to virial shock heating, the mass floor M min explains galaxy "downsizing," where more massive galaxies formed earlier and over a shorter period of time. It turns out that the model also accounts for the observed galactic baryon and gas fractions as a function of mass and time, and the cosmic SFR density, which are all resulting from the mass floor M min. The model helps us to understand that it is the cosmological decline of accretion rate that drives the decrease of cosmic SFR density between z ∼ 2 and z = 0 and the rise of the cosmic SFR density from z ∼ 6 to z ∼ 2 that allows us to put a constraint on our main parameter M min 10 11 M. Among the physical mechanisms that could be responsible for the mass floor, our view is that photoionization feedback (from first in situ hot stars) lowering the cooling efficiency is likely to play a large role.
Abstract
We present a high-resolution kinematic study of the massive main-sequence star-forming galaxy (SFG) SDSS J090122.37+181432.3 (J0901) at
z
= 2.259, using ∼0.″36 Atacama Large ...Millimeter/submillimeter Array CO(3–2) and ∼0.″1–0.″5 SINFONI/VLT H
α
observations. J0901 is a rare, strongly lensed but otherwise normal massive (
log
(
M
⋆
/
M
⊙
)
∼
11
) main-sequence SFG, offering a unique opportunity to study a typical massive SFG under the microscope of lensing. Through forward dynamical modeling incorporating lensing deflection, we fit the CO and H
α
kinematics in the image plane out to about one disk effective radius (
R
e
∼ 4 kpc) at an ∼600 pc delensed physical resolution along the kinematic major axis. Our results show high intrinsic dispersions of the cold molecular and warm ionized gas (
σ
0,mol.
∼ 40 km s
−1
and
σ
0,ion.
∼ 66 km s
−1
) that remain constant out to
R
e
; a moderately low dark matter fraction (
f
DM
∼ 0.3–0.4) within
R
e
; and a centrally peaked Toomre
Q
parameter—agreeing well with the previously established
σ
0
versus
z
,
f
DM
versus Σ
baryon
, and
Q
's radial trends using large-sample non-lensed main-sequence SFGs. Our data further reveal a high stellar mass concentration within ∼1–2 kpc with little molecular gas, and a clumpy molecular gas ring-like structure at
R
∼ 2–4 kpc, in line with the inside-out quenching scenario. Our further analysis indicates that J0901 had assembled half of its stellar mass only ∼400 Myr before its observed cosmic time, and the cold gas ring and dense central stellar component are consistent with signposts of a recent wet compaction event of a highly turbulent disk found in recent simulations.