Stars form from cold molecular interstellar gas. As this is relatively rare in the local Universe, galaxies like the Milky Way form only a few new stars per year. Typical massive galaxies in the ...distant Universe formed stars an order of magnitude more rapidly. Unless star formation was significantly more efficient, this difference suggests that young galaxies were much more molecular-gas rich. Molecular gas observations in the distant Universe have so far largely been restricted to very luminous, rare objects, including mergers and quasars, and accordingly we do not yet have a clear idea about the gas content of more normal (albeit massive) galaxies. Here we report the results of a survey of molecular gas in samples of typical massive-star-forming galaxies at mean redshifts of about 1.2 and 2.3, when the Universe was respectively 40% and 24% of its current age. Our measurements reveal that distant star forming galaxies were indeed gas rich, and that the star formation efficiency is not strongly dependent on cosmic epoch. The average fraction of cold gas relative to total galaxy baryonic mass at z = 2.3 and z = 1.2 is respectively about 44% and 34%, three to ten times higher than in today’s massive spiral galaxies. The slow decrease between z 2 and z 1 probably requires a mechanism of semi-continuous replenishment of fresh gas to the young galaxies.
Aims. We investigate the fueling and the feedback of nuclear activity in the nearby (D = 14 Mpc) Seyfert 2 barred galaxy NGC 1068 by studying the distribution and kinematics of molecular gas in the ...torus and its connections to the host galaxy disk. Methods. We used the Atacama Large Millimeter Array (ALMA ) to image the emission of a set of molecular gas tracers in the circumnuclear disk (CND) and the torus of the galaxy using the CO(2–1), CO(3–2), and HCO+(4–3) lines and their underlying continuum emission with high spatial resolutions (0.03″ − 0.09″ ≃ 2 − 6 pc). These transitions, which span a wide range of physical conditions of molecular gas (n(H2)⊂103 − 107 cm−3), are instrumental in revealing the density radial stratification and the complex kinematics of the gas in the torus and its surroundings. Results. The ALMA images resolve the CND as an asymmetric ringed disk of D ≃ 400 pc in size and ≃1.4 × 108 M⊙ in mass. The CND shows a marked deficit of molecular gas in its central ≃130 pc region. The inner edge of the ring is associated with the presence of edge-brightened arcs of NIR polarized emission, which are identified with the current working surface of the ionized wind of the active galactic nucleus (AGN). ALMA proves the existence of an elongated molecular disk/torus in NGC 1068 of Mtorusgas ≃ 3 × 105 M⊙ M torus gas ≃ 3 × 10 5 M ⊙ $ M_{\mathrm{torus}}^{\mathrm{gas}}\simeq3\times10^{5}\,M_{{\odot}} $ , which extends over a large range of spatial scales D ≃ 10 − 30 pc around the central engine. The new observations evidence the density radial stratification of the torus: the HCO+(4–3) torus, with a full size DHCO+(4 − 3) = 11 ± 0.6 pc, is a factor of between two and three smaller than its CO(2–1) and CO(3–2) counterparts, which have full sizes of DCO(3 − 2) = 26 ± 0.6 pc and DCO(2 − 1) = 28 ± 0.6 pc, respectively. This result brings into light the many faces of the molecular torus. The torus is connected to the CND through a network of molecular gas streamers detected inside the CND ring. The kinematics of molecular gas show strong departures from circular motions in the torus, the gas streamers, and the CND ring. These velocity field distortions are interconnected and are part of a 3D outflow that reflects the effects of AGN feedback on the kinematics of molecular gas across a wide range of spatial scales around the central engine. In particular, we estimate through modeling that a significant fraction of the gas inside the torus ( ≃ 0.4 − 0.6 × Mtorusgas ≃ 0.4 − 0.6 × M torus gas $ {\simeq}0.4{-}0.6 \times M_{\mathrm{torus}}^{\mathrm{gas}} $ ) and a comparable amount of mass along the gas streamers are outflowing. However, the bulk of the mass, momentum, and energy of the molecular outflow of NGC 1068 is contained at larger radii in the CND region, where the AGN wind and the radio jet are currently pushing the gas assembled at the Inner Lindblad Resonance (ILR) ring of the nuclear stellar bar. Conclusions. In our favored scenario a wide-angle AGN wind launched from the accretion disk of NGC1068 is currently impacting a sizable fraction of the gas inside the torus. However, a large gas reservoir (≃1.2 − 1.8 × 105 M⊙), which lies close to the equatorial plane of the torus, remains unaffected by the feedback of the AGN wind and can therefore continue fueling the AGN for at least ≃1 − 4 Myr. Nevertheless, AGN fueling currently seems thwarted on intermediate scales (15 pc ≤r ≤ 50 pc).
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.
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.
We use the first systematic data sets of CO molecular line emission in z∼ 1–3 normal star-forming galaxies (SFGs) for a comparison of the dependence of galaxy-averaged star formation rates on ...molecular gas masses at low and high redshifts, and in different galactic environments. Although the current high-z samples are still small and biased towards the luminous and massive tail of the actively star-forming ‘main-sequence’, a fairly clear picture is emerging. Independent of whether galaxy-integrated quantities or surface densities are considered, low- and high-z SFG populations appear to follow similar molecular gas–star formation relations with slopes 1.1 to 1.2, over three orders of magnitude in gas mass or surface density. The gas-depletion time-scale in these SFGs grows from 0.5 Gyr at z∼ 2 to 1.5 Gyr at z∼ 0. The average corresponds to a fairly low star formation efficiency of 2 per cent per dynamical time. Because star formation depletion times are significantly smaller than the Hubble time at all redshifts sampled, star formation rates and gas fractions are set by the balance between gas accretion from the halo and stellar feedback. In contrast, very luminous and ultraluminous, gas-rich major mergers at both low and high z produce on average four to 10 times more far-infrared luminosity per unit gas mass. We show that only some fraction of this difference can be explained by uncertainties in gas mass or luminosity estimators; much of it must be intrinsic. A possible explanation is a top-heavy stellar mass function in the merging systems but the most likely interpretation is that the star formation relation is driven by global dynamical effects. For a given mass, the more compact merger systems produce stars more rapidly because their gas clouds are more compressed with shorter dynamical times, so that they churn more quickly through the available gas reservoir than the typical normal disc galaxies. When the dependence on galactic dynamical time-scale is explicitly included, disc galaxies and mergers appear to follow similar gas-to-star formation relations. The mergers may be forming stars at slightly higher efficiencies than the discs.
We report ALMA observations of CO(3–2) emission in the Seyfert 2 galaxy NGC 1433 at the unprecedented spatial resolution of 0.′′5 = 24 pc. Our aim is to probe active galactic nucleus (AGN) feeding ...and feedback phenomena through the morphology and dynamics of the gas inside the central kpc. The galaxy NGC 1433 is a strongly barred spiral with three resonant rings: one at the ultra-harmonic resonance near corotation, and the others at the outer and inner Lindblad resonances (OLR and ILR). A nuclear bar of 400 pc radius is embedded in the large-scale primary bar. The CO map, which covers the whole nuclear region (nuclear bar and ring), reveals a nuclear gaseous spiral structure, inside the nuclear ring encircling the nuclear stellar bar. This gaseous spiral is well correlated with the dusty spiral seen in Hubble Space Telescope (HST) images. The nuclear spiral winds up in a pseudo-ring at ~200 pc radius, which might correspond to the inner ILR. Continuum emission is detected at 0.87 mm only at the very centre, and its origin is more likely thermal dust emission than non-thermal emission from the AGN. It might correspond to the molecular torus expected to exist in this Seyfert 2 galaxy. The HCN(4–3) and HCO+(4–3) lines were observed simultaneously, but only upper limits are derived, with a ratio to the CO(3–2) line lower than 1/60 at 3σ, indicating a relatively low abundance of very dense gas. The kinematics of the gas over the nuclear disk reveal rather regular rotation only slightly perturbed by streaming motions due to the spiral; the primary and secondary bars are too closely aligned with the galaxy major or minor axis to leave a signature in the projected velocities. Near the nucleus, there is an intense high-velocity CO emission feature redshifted to 200 km s-1 (if located in the plane), with a blue-shifted counterpart, at 2′′ (100 pc) from the centre. While the CO spectra are quite narrow in the centre, this wide component is interpreted as an outflow involving a molecular mass of 3.6 × 106 M⊙ and a flow rate ~7 M⊙/yr. The flow could be in part driven by the central star formation, but is mainly boosted by the AGN through its radio jets.
The broadening of atomic emission lines by high-velocity motion of gas near accreting supermassive black holes is an observational hallmark of quasars
. Observations of broad emission lines could ...potentially constrain the mechanism for transporting gas inwards through accretion disks or outwards through winds
. The size of regions for which broad emission lines are observed (broad-line regions) has been estimated by measuring the delay in light travel time between the variable brightness of the accretion disk continuum and the emission lines
-a method known as reverberation mapping. In some models the emission lines arise from a continuous outflow
, whereas in others they arise from orbiting gas clouds
. Directly imaging such regions has not hitherto been possible because of their small angular size (less than 10
arcseconds
). Here we report a spatial offset (with a spatial resolution of 10
arcseconds, or about 0.03 parsecs for a distance of 550 million parsecs) between the red and blue photo-centres of the broad Paschen-α line of the quasar 3C 273 perpendicular to the direction of its radio jet. This spatial offset corresponds to a gradient in the velocity of the gas and thus implies that the gas is orbiting the central supermassive black hole. The data are well fitted by a broad-line-region model of a thick disk of gravitationally bound material orbiting a black hole of 3 × 10
solar masses. We infer a disk radius of 150 light days; a radius of 100-400 light days was found previously using reverberation mapping
. The rotation axis of the disk aligns in inclination and position angle with the radio jet. Our results support the methods that are often used to estimate the masses of accreting supermassive black holes and to study their evolution over cosmic time.
Abstract
We derive 2D dust attenuation maps at ∼1 kpc resolution from the UV continuum for 10 galaxies on the
z
∼ 2 star-forming main sequence (SFMS). Comparison with IR data shows that 9 out of 10 ...galaxies do not require further obscuration in addition to the UV-based correction, though our sample does not include the most heavily obscured, massive galaxies. The individual rest-frame
V
-band dust attenuation (
A
V
) radial profiles scatter around an average profile that gently decreases from ∼1.8 mag in the center down to ∼0.6 mag at ∼3–4 half-mass radii. We use these maps to correct UV- and H
α
-based star formation rates (SFRs), which agree with each other. At masses
, the dust-corrected specific SFR (sSFR) profiles are on average radially constant at a mass-doubling timescale of ∼300 Myr, pointing at a synchronous growth of bulge and disk components. At masses
, the sSFR profiles are typically centrally suppressed by a factor of ∼10 relative to the galaxy outskirts. With total central obscuration disfavored, this indicates that at least a fraction of massive
z
∼ 2 SFMS galaxies have started their inside-out star formation quenching that will move them to the quenched sequence. In combination with other observations, galaxies above and below the ridge of the SFMS relation have, respectively, centrally enhanced and centrally suppressed sSFRs relative to their outskirts, supporting a picture where bulges are built owing to gas “compaction” that leads to a high central SFR as galaxies move toward the upper envelope of the SFMS.
In the cold dark matter cosmology, the baryonic components of galaxies-stars and gas-are thought to be mixed with and embedded in non-baryonic and non-relativistic dark matter, which dominates the ...total mass of the galaxy and its dark-matter halo. In the local (low-redshift) Universe, the mass of dark matter within a galactic disk increases with disk radius, becoming appreciable and then dominant in the outer, baryonic regions of the disks of star-forming galaxies. This results in rotation velocities of the visible matter within the disk that are constant or increasing with disk radius-a hallmark of the dark-matter model. Comparisons between the dynamical mass, inferred from these velocities in rotational equilibrium, and the sum of the stellar and cold-gas mass at the peak epoch of galaxy formation ten billion years ago, inferred from ancillary data, suggest high baryon fractions in the inner, star-forming regions of the disks. Although this implied baryon fraction may be larger than in the local Universe, the systematic uncertainties (owing to the chosen stellar initial-mass function and the calibration of gas masses) render such comparisons inconclusive in terms of the mass of dark matter. Here we report rotation curves (showing rotation velocity as a function of disk radius) for the outer disks of six massive star-forming galaxies, and find that the rotation velocities are not constant, but decrease with radius. We propose that this trend arises because of a combination of two main factors: first, a large fraction of the massive high-redshift galaxy population was strongly baryon-dominated, with dark matter playing a smaller part than in the local Universe; and second, the large velocity dispersion in high-redshift disks introduces a substantial pressure term that leads to a decrease in rotation velocity with increasing radius. The effect of both factors appears to increase with redshift. Qualitatively, the observations suggest that baryons in the early (high-redshift) Universe efficiently condensed at the centres of dark-matter haloes when gas fractions were high and dark matter was less concentrated.
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
.