The SHINING survey offers a great opportunity to study the properties of the ionized and neutral media of galaxies from prototypical starbursts and active galactic nuclei (AGNs) to heavily obscured ...objects. Based on Herschel/PACS observations of the main far-infrared (FIR) fine-structure lines, in this paper, we analyze the physical mechanisms behind the observed line deficits in galaxies, the apparent offset of luminous infrared galaxies (LIRGs) from the mass-metallicity relation, and the scaling relations between C ii 158 m line emission and star formation rate (SFR). Based on a toy model and the Cloudy code, we conclude that the increase in the ionization parameter with FIR surface brightness can explain the observed decrease in the line-to-FIR continuum ratio of galaxies. In the case of the C ii line, the increase in the ionization parameter is accompanied by a reduction in the photoelectric heating efficiency and the inability of the line to track the increase in the FUV radiation field as galaxies become more compact and luminous. In the central approximately kiloparsec regions of AGN galaxies, we observe a significant increase in the O i 63 m/C ii line ratio; the AGN impact on the line-to-FIR ratios fades on global scales. Based on extinction-insensitive metallicity measurements of LIRGs, we confirm that they lie below the mass-metallicity relation, but the offset is smaller than those reported in studies that use optical-based metal abundances. Finally, we present scaling relations between C ii emission and SFR in the context of the main sequence of star-forming galaxies.
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.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
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 Herschel/PACS spectrometer to study the global and spatially resolved far-infrared (FIR) fine-structure line emission in a sample of 52 galaxies that constitute the SHINING survey. These ...galaxies include star-forming, active-galactic nuclei (AGNs), and luminous infrared galaxies (LIRGs). We find an increasing number of galaxies (and kiloparsec-size regions within galaxies) with low line-to-FIR continuum ratios as a function of increasing FIR luminosity (LFIR), dust infrared color, LFIR to molecular gas mass ratio (LFIR/Mmol), and FIR surface brightness ( FIR). The correlations between the C ii/FIR or O i/FIR ratios with FIR are remarkably tight (∼0.3 dex scatter over almost four orders of magnitude in FIR). We observe that galaxies with and FIR 1011 L kpc−2 tend to have weak fine-structure line-to-FIR continuum ratios, and that LIRGs with infrared sizes 1 kpc have line-to-FIR ratios comparable to those observed in typical star-forming galaxies. We analyze the physical mechanisms driving these trends in Paper II. The combined analysis of the C ii, N ii 122 m, and O iii 88 m lines reveals that the fraction of the C ii line emission that arises from neutral gas increases from 60% to 90% in the most active star-forming regions and that the emission originating in the ionized gas is associated with low-ionization, diffuse gas rather than with dense gas in H ii regions. Finally, we report the global and spatially resolved line fluxes of the SHINING galaxies to enable the comparison and planning of future local and high-z studies.
We use Herschel 70 to 160 μm images to study the size of the far-infrared emitting region in about 400 local galaxies and quasar (QSO) hosts. The sample includes normal “main-sequence” star-forming ...galaxies, as well as infrared luminous galaxies and Palomar-Green QSOs, with different levels and structures of star formation. Assuming Gaussian spatial distribution of the far-infrared (FIR) emission, the excellent stability of the Herschel point spread function (PSF) enables us to measure sizes well below the PSF width, by subtracting widths in quadrature. We derive scalings of FIR size and surface brightness of local galaxies with FIR luminosity, with distance from the star-forming main-sequence, and with FIR color. Luminosities LFIR~ 1011 L⊙ can be reached with a variety of structures spanning 2 dex in size. Ultraluminous LFIR≳ 1012 L⊙ galaxies far above the main-sequence inevitably have small Re,70~ 0.5 kpc FIR emitting regions with large surface brightness, and can be close to optically thick in the FIR on average over these regions. Compared to these local relations, first ALMA sizes for the dust emission regions in high redshift galaxies, measured at somewhat longer rest wavelengths, suggest larger sizes at the same IR luminosity. We report a remarkably tight relation with 0.15 dex scatter between FIR surface brightness and the ratio of Cii 158 μm emission and FIR emission – the so-called Cii-deficit is more tightly linked to surface brightness than to FIR luminosity or FIR color. Among 33 z ≤ 0.1 PG QSOs with typical LFIR/LBol,AGN ≈ 0.1, 19 have a measured 70 μm half light radius, with median Re,70 = 1.1 kpc. This is consistent with the FIR size for galaxies with similar LFIR but lacking a QSO, in accordance with a scenario where the rest FIR emission of these types of QSOs is, in most cases, due to host star formation.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
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.
In this first paper on the results of our Herschel PACS survey of local ultra luminous infrared galaxies (ULIRGs), as part of our SHINING survey of local galaxies, we present far-infrared ...spectroscopy of Mrk 231, the most luminous of the local ULIRGs, and a type 1 broad absorption line AGN. For the first time in a ULIRG, all observed far-infrared fine-structure lines in the PACS range were detected and all were found to be deficient relative to the far infrared luminosity by 1–2 orders of magnitude compared with lower luminosity galaxies. The deficits are similar to those for the mid-infrared lines, with the most deficient lines showing high ionization potentials. Aged starbursts may account for part of the deficits, but partial covering of the highest excitation AGN powered regions may explain the remaining line deficits. A massive molecular outflow, discovered in OH and 18OH, showing outflow velocities out to at least 1400 km s-1, is a unique signature of the clearing out of the molecular disk that formed by dissipative collapse during the merger. The outflow is characterized by extremely high ratios of 18O/16O suggestive of interstellar medium processing by advanced starbursts.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
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.
We report the results from a systematic search for molecular (OH 119 micron) outflows with Herschel/PACS in a sample of 43 nearby (z < 0.3) galaxy mergers, mostly ultraluminous infrared galaxies ...(ULIRGs) and QSOs. We find that the character of the OH feature (strength of the absorption relative to the emission) correlates with that of the 9.7 micron silicate feature, a measure of obscuration in ULIRGs. Unambiguous evidence for molecular outflows, based on the detection of OH absorption profiles with median velocities more blueshifted than−50 km/s, is seen in 26 (70%) of the 37 OH-detected targets, suggesting a wide-angle (approx. 145 deg.) outflow geometry. Conversely, unambiguous evidence for molecular inflows, based on the detection of OH absorption profiles with median velocities more redshifted than +50 km/s is seen in only four objects, suggesting a planar or filamentary geometry for the inflowing gas. Terminal outflow velocities of approx. −1000 km/s are measured in several objects, but median outflow velocities are typically approx.−200 km/s−1. While the outflow velocities show no statistically significant dependence on the star formation rate, they are distinctly more blueshifted among systems with large active galactic nucleus (AGN) fractions and luminosities log (L(sub AGN)/L(sub solar)) => 11.8 +/- 0.3. The quasars in these systems play a dominant role in driving the molecular outflows. However, the most AGN dominated systems, where OH is seen purely in emission, show relatively modest OH line widths, despite their large AGN luminosities, perhaps indicating that molecular outflows subside once the quasar has cleared a path through the obscuring material.
We combine two approaches to isolate the AGN luminosity at near-IR wavelengths and relate the near-IR pure AGN luminosity to other tracers of the AGN. Using integral-field spectroscopic data of an ...archival sample of 51 local AGNs, we estimate the fraction of non-stellar light by comparing the nuclear equivalent width of the stellar 2.3 μm CO absorption feature with the intrinsic value for each galaxy. We compare this fraction to that derived from a spectral decomposition of the integrated light in the central arcsecond and find them to be consistent with each other. Using our estimates of the near-IR AGN light, we find a strong correlation with presumably isotropic AGN tracers. We show that a significant offset exists between type 1 and type 2 sources in the sense that type 1 sources are 7 (10) times brighter in the near-IR at \hbox{$\log \lmir{}$}logLAGNMIR = 42.5 (\hbox{$\log \lx{}$}logLAGNX = 42.5). These offsets only become clear when treating infrared type 1 sources as type 1 AGNs. All AGNs have very red near- to mid-IR dust colors. This, as well as the range of observed near-IR temperatures, can be explained with a simple model with only two free parameters: the obscuration to the hot dust and the ratio between the warm and hot dust areas. We find obscurations of \hbox{$A_{V}^{\rm hot} = 5 \ldots 15$}AVhot=5...15 mag for infrared type 1 sources and \hbox{$A_{V}^{\rm hot} = 15 \ldots 35$}AVhot=15...35 mag for type 2 sources. The ratio of hot dust to warm dust areas of about 1000 is nicely consistent with the ratio of radii of the respective regions as found by infrared interferometry.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK