Local ULIRGs Research Articles

Cosmic Sands: The Origin of Dusty, Star-forming Galaxies in the Epoch of Reionization

We present the Cosmic Sands suite of cosmological zoom-in simulations based on the simba galaxy formation model in order to study the buildup of the first massive and dusty galaxies in the early universe. Residing in the most massive halos, we find that the compact proto-massive galaxies undergo nearly continuous mergers with smaller subhalos, boosting star formation rates (SFRs) and the buildup of stellar mass. The galaxies are already appreciably chemically evolved by z = 7, with modeled dust masses comparable to those inferred from observations in the same epoch, except for the most extreme systems. We track gas accretion onto the galaxies to understand how extreme SFRs can be sustained by these early systems. We find that smooth gas accretion can maintain SFRs above 250 M ⊙ yr−1, but to achieve SFRs that boost galaxies well above the main sequence, a larger perturbation like a gas-rich major merger is necessary to trigger a starburst episode. Post-processing the Cosmic Sands simulations with dust RT, we find that, while the infrared luminosities of the most-dust-rich galaxies are comparable to local ULIRGs, they are substantially dimmer than classical z = 2 submillimeter galaxies. We end with a discussion on the possible reasons for this discrepancy at the highest masses and the future work we intend to carry out to study the chemical enrichment of the earliest dusty galaxies.

A technique to select the most obscured galaxy nuclei

Compact obscured nuclei (CONs) are mainly found in local luminous and ultraluminous infrared galaxies (U/LIRGs). In the local Universe, these sources are generally selected through the detection of the HCN–vib (3-2) emission line at submillimetre wavelengths. In this work, we present a diagnostic method to select deeply buried nuclei based on mid-infrared (mid-IR) polycyclic aromatic hydrocarbons (PAHs) and mid-IR continuum ratios. Using Spitzer InfraRed Spectrograph (IRS) spectra of a representative sample of local ULIRGs (z < 0.27), we examine their PAH and underlying continuum emission ratios. For deeply embedded sources, we find that the 9.7 μm silicate absorption band has a particularly pronounced effect on the 11.3 μm PAH feature. The low flux level in the nuclear silicate absorption band enhances the 11.3 μm PAH feature contrast (high PAH equivalent width) compared to that of the other PAH features. The technique has been extended to include the use of the underlying 11.3/12.7 and 11.3/6.2 μm continuum ratios. However, the latter are affected by the extinction coming from both the host galaxy and the nuclear region, whereas the foreground (host-galaxy) extinction is cancelled out when using the PAH equivalent width ratios. We apply our method to local U/LIRGs from the HERUS and GOALS samples and classify 14 ULIRGs and 10 LIRGs as CON candidates, which corresponds to 30% of the ULIRGs and 7% of the LIRGs from these samples. We find that the observed continuum ratios of CON-dominated sources can be explained by assuming torus models with a tapered disc geometry and a smooth dust distribution. This suggests that the nuclear dusty structure of deeply obscured galaxy nuclei has an extremely high dust coverage. Finally, we demonstrate that the use of mid-IR colour–colour diagrams is an effective way to select CON-dominated sources at different redshifts. In particular, the combination of filters of the James Webb Space Telescope/Mid-Infrared Instrument will enable the selection of CONs out to z ∼ 1.5. This will allow the selection of CONs to be extended to high redshifts where U/LIRGs are more numerous.

Physics of ULIRGs with MUSE and ALMA: The PUMA project

We analyze new high-resolution (400 pc) ∼220 GHz continuum and CO(2–1) Atacama Large Millimeter Array (ALMA) observations of a representative sample of 23 local (z < 0.165) ultra-luminous infrared systems (ULIRGs; 34 individual nuclei) as part of the “Physics of ULIRGs with MUSE and ALMA” (PUMA) project. The deconvolved half-light radii of the ∼220 GHz continuum sources, rcont, are between < 60 pc and 350 pc (median 80–100 pc). We associate these regions with the regions emitting the bulk of the infrared luminosity (LIR). The good agreement, within a factor of 2, between the observed ∼220 GHz fluxes and the extrapolation of the infrared gray-body as well as the small contributions from synchrotron and free–free emission support this assumption. The cold molecular gas emission sizes, rCO, are between 60 and 700 pc and are similar in advanced mergers and early interacting systems. On average, rCO are ∼2.5 times larger than rcont. Using these measurements, we derived the nuclear LIR and cold molecular gas surface densities (ΣLIR = 1011.5 − 1014.3 L⊙ kpc−2 and ΣH2 = 102.9 − 104.2 M⊙ pc−2, respectively). Assuming that the LIR is produced by star formation, the median ΣLIR corresponds to ΣSFR = 2500 M⊙ yr−1 kpc−2. This ΣSFR implies extremely short depletion times, ΣH2/ΣSFR < 1–15 Myr, and unphysical star formation efficiencies > 1 for 70% of the sample. Therefore, this favors the presence of an obscured active galactic nucleus (AGN) in these objects that could dominate the LIR. We also classify the ULIRG nuclei in two groups: (a) compact nuclei (rcont < 120 pc) with high mid-infrared excess emission (ΔL6−20 μm/LIR) found in optically classified AGN; and (b) nuclei following a relation with decreasing ΔL6−20 μm/LIR for decreasing rcont. The majority, 60%, of the nuclei in interacting systems lie in the low-rcont end (<120 pc) of this relation, while this is the case for only 30% of the mergers. This suggests that in the early stages of the interaction, the activity occurs in a very compact and dust-obscured region while, in more advanced merger stages, the activity is more extended, unless an optically detected AGN is present. Approximately two-thirds of the nuclei have nuclear radiation pressures above the Eddington limit. This is consistent with the ubiquitous detection of massive outflows in local ULIRGs and supports the importance of the radiation pressure in the outflow launching process.

AGN and star formation across cosmic time

ABSTRACT We investigate the balance of power between stars and AGN across cosmic history, based on the comparison between the infrared (IR) galaxy luminosity function (LF) and the IR AGN LF. The former corresponds to emission from dust heated by stars and AGN, whereas the latter includes emission from AGN-heated dust only. We find that at all redshifts (at least up to z ∼ 2.5), the high-luminosity tails of the two LFs converge, indicating that the most IR-luminous galaxies are AGN-powered. Our results shed light to the decades-old conundrum regarding the flatter high-luminosity slope seen in the IR galaxy LF compared to that in the UV and optical. We attribute this difference to the increasing fraction of AGN-dominated galaxies with increasing total IR luminosity (LIR). We partition the LIR−z parameter space into a star formation-dominated and an AGN-dominated region, finding that the most luminous galaxies at all epochs lie in the AGN-dominated region. This sets a potential ‘limit’ to attainable star formation rates, casting doubt on the abundance of ‘extreme starbursts’: if AGN did not exist, LIR > 1013 L⊙ galaxies would be significantly rarer than they currently are in our observable Universe. We also find that AGN affect the average dust temperatures (Tdust) of galaxies and hence the shape of the well-known LIR−Tdust relation. We propose that the reason why local ULIRGs are hotter than their high-redshift counterparts is because of a higher fraction of AGN-dominated galaxies amongst the former group.

Extreme CO Isotopologue Line Ratios in ULIRGS: Evidence for a Top-heavy IMF

Abstract We present high-resolution ALMA observations of the C18O, 13CO, and 12CO J = 1–0 isotopologues in three nearby ultraluminous infrared galaxies (ULIRGS; Arp 220, IRAS 13120-5453, and IRAS 17208-0014) and one nearby post-merger galaxy (NGC 2623). In all four systems, we measure high 12CO/C18O and 12CO/13CO integrated line ratios while the 13CO/C18O ratio is observed to be extremely low in comparison to typical star-forming disks, supporting previous work. We investigate whether these unusual line ratios are due to dynamical effects, astrochemistry within the gas, or nucleosynthesis in stars. Assuming both lines are optically thin, low 13CO/C18O values suggest that C18O is more abundant than 13CO in the interstellar medium of these systems. A plausible explanation is that local ULIRGs and their progeny have an excess in massive star formation; in other words, they are producing a top-heavy stellar initial mass function.

Resolving Star Formation on Subkiloparsec Scales in the High-redshift Galaxy SDP.11 Using Gravitational Lensing*

Abstract We investigate the properties of the interstellar medium, star formation, and the current-day stellar population in the strongly lensed star-forming galaxy H-ATLAS J091043.1-000321 (SDP.11), at z = 1.7830, using new Herschel and Atacama Large Millimeter/submillimeter Array (ALMA) observations of far-infrared fine-structure lines of carbon, oxygen, and nitrogen. We report detections of the [O iii] 52 μm, [N iii] 57 μm, and [O i] 63 μm lines from Herschel/PACS, and present high-resolution imaging of the [C ii] 158 μm line, and underlying continuum, using ALMA. We resolve the [C ii] line emission into two spatially offset Einstein rings, tracing the red and blue velocity components of the line, in the ALMA/Band 9 observations at 0.″2 resolution. The values seen in the [C ii]/far-infrared (FIR) ratio map, as low as ∼0.02% at the peak of the dust continuum, are similar to those of local ULIRGs, suggesting an intense starburst in this source. This is consistent with the high intrinsic FIR luminosity (∼3 × 1012 L ⊙), ∼16 Myr gas depletion timescale, and ≲8 Myr timescale since the last starburst episode, estimated from the hardness of the UV radiation field. By applying gravitational lensing models to the visibilities in the uv-plane, we find that the lensing magnification factor varies by a factor of two across SDP.11, affecting the observed line profiles. After correcting for the effects of differential lensing, a symmetric line profile is recovered, suggesting that the starburst present here may not be the result of a major merger, as is the case for local ULIRGs, but instead could be powered by star formation activity spread across a 3–5 kpc rotating disk.

PAHs as tracers of the molecular gas in star-forming galaxies

[Abridged] We combine new CO(1-0) line observations of 24 intermediate redshift galaxies (0.03 < z < 0.28) along with literature data of galaxies at 0<z<4 to explore scaling relations between the dust and gas content using PAH 6.2 $\mu$m ($L_{6.2}$), CO ($L'_{\rm CO}$), and infrared ($L_{\rm IR}$) luminosities for a wide range of redshifts and physical environments. Our analysis confirms the existence of a universal $L_{6.2}-L_{\rm CO}$ correlation followed by normal star-forming galaxies (SFGs) and starbursts (SBs) at all redshifts. This relation is also followed by local ULIRGs that appear as outliers in the $L_{6.2}-L_{\rm IR}$ and $L_{\rm IR}-L'_{\rm CO}$ relations from the sequence defined by normal SFGs. The emerging tight ($\sigma \approx 0.26$ dex) and linear ($\alpha = 1.03$) relation between $L_{6.2}$ and $L_{\rm CO}$ indicates a $L_{6.2}$ to molecular gas ($M_{\rm H_2}$) conversion factor of $\alpha_{6.2} = M_{\rm H2}/L_{6.2} = (2.7\pm1.3) \times \alpha_{\rm CO}$, where $\alpha_{\rm CO}$ is the $L'_{\rm CO}$ to $M_{\rm H_2}$ conversion factor. We also find that on galaxy integrated scales, PAH emission is better correlated with cold rather than with warm dust emission, suggesting that PAHs are associated with the diffuse cold dust, which is another proxy for $M_{\rm H_2}$. Focusing on normal SFGs among our sample, we employ the dust continuum emission to derive $M_{\rm H_2}$ estimates and find a constant $M_{\rm H_2}/L_{6.2}$ ratio of $\alpha_{6.2} = 12.3 \ M_{\rm H_2}/{\rm L}_{\odot}$ ($\sigma\approx 0.3$ dex). We propose that the presented $L_{6.2}-L'_{\rm CO}$ and $L_{6.2}-M_{\rm H_2}$ relations will serve as useful tools for the determination of the physical properties of high-$z$ SFGs, for which PAH emission will be routinely detected by the James Webb Space Telescope.

Far-infrared Herschel SPIRE spectroscopy of lensed starbursts reveals physical conditions of ionized gas

The most intensively star-forming galaxies are extremely luminous at far-infrared (FIR) wavelengths, highly obscured at optical and ultraviolet wavelengths, and lie at $z\ge 1-3$. We present a programme of ${\it Herschel}$ FIR spectroscopic observations with the SPIRE FTS and photometric observations with PACS, both on board ${\it Herschel}$, towards a sample of 45 gravitationally lensed, dusty starbursts across $z\sim 1-3.6$. In total, we detected 27 individual lines down to 3-$\sigma$, including nine $[\rm C{\small II}]$ 158-$\mu$m lines with confirmed spectroscopic redshifts, five possible $[\rm C{\small II}]$ lines consistent with their far-infrared photometric redshifts, and in some individual sources a few $[\rm O{\small III}]$ 88-$\mu$m, $[\rm O{\small III}]$ 52-$\mu$m, $[\rm O{\small I}]$ 145-$\mu$m, $[\rm O{\small I}]$ 63-$\mu$m, $[\rm N{\small II}]$ 122-$\mu$m, and OH 119-$\mu$m (in absorption) lines. To derive the typical physical properties of the gas in the sample, we stack all spectra weighted by their intrinsic luminosity and by their 500-$\mu$m flux densities, with the spectra scaled to a common redshift. In the stacked spectra, we detect emission lines of $[\rm C{\small II}]$ 158-$\mu$m, $[\rm N{\small II}]$ 122-$\mu$m, $[\rm O{\small III}]$ 88-$\mu$m, $[\rm O{\small III}]$ 52-$\mu$m, $[\rm O{\small I}]$ 63-$\mu$m, and the absorption doublet of OH at 119-$\mu$m, at high fidelity. We find that the average electron densities traced by the $[\rm N{\small II}]$ and $[\rm O{\small III}]$ lines are higher than the average values in local star-forming galaxies and ULIRGs, using the same tracers. From the $[\rm N{\small II}]/[\rm C{\small II}]$ and $[\rm O{\small I}]/[\rm C{\small II}]$ ratios, we find that the $[\rm C{\small II}]$ emission is likely dominated by the photo-dominated regions (PDR), instead of by ionised gas or large-scale shocks.

Quantifying the AGN-driven outflows in ULIRGs (QUADROS) – I: VLT/Xshooter observations of nine nearby objects

Although now routinely incorporated into hydrodynamic simulations of galaxy evolution, the true importance of the feedback effect of the outflows driven by active galactic nuclei (AGN) remains uncertain from an observational perspective. This is due to a lack of accurate information on the densities, radial scales and level of dust extinction of the outflow regions. Here we use the unique capabilities of VLT/Xshooter to investigate the warm outflows in a representative sample of 9 local (0.06 $<$ z $<$ 0.15) ULIRGs with AGN nuclei and, for the first time, accurately quantify the key outflow properties. We find that the outflows are compact (0.05 $<$ R$_{[OIII]}$ $<$ 1.2 kpc), significantly reddened (median E(B-V)$\sim$0.5 magnitudes), and have relatively high electron densities (3.4 $<$ log$_{10}$ n$_e$ (cm$^{-3}$) $<$ 4.8). It is notable that the latter densities -- obtained using trans-auroral [SII] and [OII] emission-line ratios -- exceed those typically assumed for the warm, emission-line outflows in active galaxies, but are similar to those estimated for broad and narrow absorption line outflow systems detected in some type 1 AGN. Even if we make the most optimistic assumptions about the true (deprojected) outflow velocities, we find relatively modest mass outflow rates ($0.07 < \dot{M} < 11$ M$_{sol}$ yr$^{-1}$) and kinetic powers measured as a fraction of the AGN bolometric luminosities ($4\times10^{-4} < \dot{E}/L_{BOL} <1$%). Therefore, although warm, AGN-driven outflows have the potential to strongly affect the star formation histories in the inner bulge regions ($r \sim$ 1kpc) of nearby ULIRGs, we lack evidence that they have a significant impact on the evolution of these rapidly evolving systems on larger scales.

Predicting HCN, HCO+, multi-transition CO, and dust emission of star-forming galaxies

High-zstar-forming galaxies have significantly higher gas fractions and star-formation efficiencies per molecular gas mass than local star-forming galaxies. In this work, we take a closer look at the gas content or fraction and the associated star-formation rate in main sequence and starburst galaxies atz= 0 andz~ 1–2 by applying an analytical model of galactic clumpy gas disks to samples of local spiral galaxies, ULIRGs, submillimeter (smm), and high-zstar-forming galaxies. The model simultaneously calculates the total gas mass, Hi/H2mass, the gas velocity dispersion, IR luminosity, IR spectral energy distribution, CO spectral line energy distribution (SLED), HCN(1–0) and HCO+(1–0) emission of a galaxy given its size, integrated star formation rate, stellar mass radial profile, rotation curve, and ToomreQparameter. The model reproduces the observed CO luminosities and SLEDs of all sample galaxies within the model uncertainties (~0.3 dex). Whereas the CO emission is robust against the variation of model parameters, the HCN and HCO+emissions are sensitive to the chemistry of the interstellar medium. The CO and HCN mass-to-light conversion factors, including CO-dark H2, are given and compared to the values found in the literature. All model conversion factors have uncertainties of a factor of two. Both the HCN and HCO+emissions trace the dense molecular gas to a factor of approximately two for the local spiral galaxies, ULIRGs and smm-galaxies. Approximately 80% of the molecular line emission of compact starburst galaxies originates in non-self-gravitating gas clouds. The effect of HCN infrared pumping is small but measurable (10–20%). The gas velocity dispersion varies significantly with the ToomreQparameter. TheQ= 1.5 model yields high-velocity dispersions (vdisp≫ 10 km s-1) consistent with available observations of high-zstar-forming galaxies and ULIRGs. However, we note that these high-velocity dispersions are not mandatory for starburst galaxies. The integrated Kennicutt-Schmidt law has a slope of approximately 1 for the local spirals, ULIRGs, and smm-galaxies, whereas the slope is 1.7 for high-zstar-forming galaxies. The model shows Kennicutt-Schmidt laws with respect to the molecular gas surface density with slopes of approximately 1.5 for local spiral galaxies, high-zstar-forming galaxies. The relation steepens for compact starburst galaxies. The model star-formation rate per unit area is, as observed, proportional to the molecular gas surface density divided by the dynamical timescale. Our relatively simple analytic model together with the recipes for the molecular line emission appears to capture the essential physics of galactic clumpy gas disks.

The Bright and Dark Sides of High-redshift Starburst Galaxies from Herschel and Subaru Observations

Abstract We present rest-frame optical spectra from the FMOS-COSMOS survey of 12 z ∼ 1.6 Herschel starburst galaxies, with star formation rate (SFR) elevated by ×8, on average, above the star-forming main sequence (MS). Comparing the Hα to IR luminosity ratio and the Balmer decrement, we find that the optically thin regions of the sources contain on average only ∼10% of the total SFR, whereas ∼90% come from an extremely obscured component that is revealed only by far-IR observations and is optically thick even in Hα. We measure the [N ii]6583/Hα ratio, suggesting that the less obscured regions have a metal content similar to that of the MS population at the same stellar masses and redshifts. However, our objects appear to be metal-rich outliers from the metallicity–SFR anticorrelation observed at fixed stellar mass for the MS population. The [S ii]6732/[S ii]6717 ratio from the average spectrum indicates an electron density n e ∼ 1100 cm−3 , larger than what was estimated for MS galaxies but only at the 1.5σ level. Our results provide supporting evidence that high-z MS outliers are analogous of local ULIRGs and are consistent with a major-merger origin for the starburst event.

The Interstellar Medium in High-redshift Submillimeter Galaxies as Probed by Infrared Spectroscopy∗

Abstract Submillimeter galaxies (SMGs) at are luminous in the far-infrared, and have star formation rates, SFR, of hundreds to thousands of solar masses per year. However, it is unclear whether they are true analogs of local ULIRGs or whether the mode of their star formation is more similar to that in local disk galaxies. We target these questions by using Herschel-PACS to examine the conditions in the interstellar medium (ISM) in far-infrared luminous SMGs at –4. We present 70–160 μm photometry and spectroscopy of the [O iv]26 μm, [Fe ii]26 μm, [S iii]33 μm, [Si ii]34 μm, [O iii]52 μm, [N iii]57 μm, and [O i]63 μm fine-structure lines and the S(0) and S(1) hydrogen rotational lines in 13 lensed SMGs identified by their brightness in early Herschel data. Most of the 13 targets are not individually spectroscopically detected; we instead focus on stacking these spectra with observations of an additional 32 SMGs from the Herschel archive—representing a complete compilation of PACS spectroscopy of SMGs. We detect [O i]63 μm, [Si ii]34 μm, and [N iii]57 μm at in the stacked spectra, determining that the average strengths of these lines relative to the far-IR continuum are , , and , respectively. Using the [O iii]52 μm/[N iii]57 μm emission line ratio, we show that SMGs have average gas-phase metallicities . By using PDR modeling and combining the new spectral measurements with integrated far-infrared fluxes and existing [C ii]158 μm data, we show that SMGs have average gas densities, n, of and FUV field strengths, (in Habing units: ), consistent with both local ULIRGs and lower luminosity star-forming galaxies.

PHIBSS: exploring the dependence of the CO–H2 conversion factor on total mass surface density at z&lt;1.5

We present an analysis of the relationship between the CO-H$_{2}$ conversion factor ($\alpha_{\rm CO}$) and total mass surface density ($\Sigma_{\rm tot}$) in star-forming galaxies at $z < 1.5$. Our sample, which is drawn from the IRAM Plateau de Bure HIgh-$z$ Blue Sequence Survey (PHIBSS) and the CO Legacy Database for GASS (COLD GASS), includes 'normal,' massive star-forming galaxies that dominate the evolution of the cosmic star formation rate (SFR) at this epoch and probe the $\Sigma_{\rm tot}$ regime where the strongest variation in $\alpha_{\rm CO}$ is observed. We constrain $\alpha_{\rm CO}$ via existing CO observations, measurements of the star formation rate, and an assumed molecular gas depletion time ($t_{\rm dep}$=$M_{\rm gas}$/SFR) --- the latter two of which establish the total molecular gas mass independent of the observed CO luminosity. For a broad range of adopted depletion times, we find that $\alpha_{\rm CO}$ is independent of total mass surface density, with little deviation from the canonical Milky Way value. This runs contrary to a scenario in which $\alpha_{\rm CO}$ decreases as surface density increases within the extended clouds of molecular gas that potentially fuel clumps of star formation in $z\sim1$ galaxies, similar to those observed in local ULIRGs. Instead, our results suggest that molecular gas, both at $z\sim0$ and $z\sim1$, is primarily in the form of self-gravitating molecular clouds. While CO observations suggest a factor of $\sim3$ reduction in the average molecular gas depletion time between $z \sim 0$ and $z\sim1$, we find that, for typical galaxies, the structure of molecular gas and the process of star formation at $z \sim 1$ is otherwise remarkably similar to that observed in local star-forming systems.

The AT-LESS CO(1-0) survey of submillimetre galaxies in the Extended Chandra Deep Field South: First results on cold molecular gas in galaxies at z ∼ 2

We present the first results from our on-going Australia Telescope Compact Array survey of CO(1-0) in ALMA-identified submillimetre galaxies in the Extended Chandra Deep Field South. Strong detections of CO(1-0) emission from two submillimetre galaxies, ALESS 122.1 (z = 2.0232) and ALESS 67.1 (z = 2.1230), were obtained. We estimate gas masses of M_gas ~ 1.3 \times 10^{11} M_odot and M_gas ~ 1.0 \times 10^{11} M_\odot for ALESS 122.1 and ALESS 67.1, respectively, adopting alpha_CO = 1.0. Dynamical mass estimates from the kinematics of the CO(1-0) line yields M_dyn (sin i)^2 = 2.1 +- 1.1 \times 10^{11} M_odot and (3.2 +- 0.9) \times 10^{11} M_\odot for ALESS 122.1 and ALESS 67.1, respectively. This is consistent with the total baryonic mass estimates of these two systems. We examine star formation efficiency using the L_FIR versus L'_CO(1-0) relation for samples of local ULIRGs and LIRGs, and more distant star-forming galaxies, with CO(1-0) detections. We find some evidence of a shallower slope for ULIRGs and SMGs compared to less luminous systems, but a larger sample is required for definite conclusions. We determine gas-to-dust ratios of 170 +- 30 and 140 +- 30 for ALESS 122.1 and ALESS 67.1, respectively, showing ALESS 122.1 has an unusually large gas reservoir. By combining the 38.1 GHz continuum detection of ALESS 122.1 with 1.4 and 5.5 GHz data, we estimate that the free-free contribution to radio emission at 38.1 GHz is 34 +- 17 microJy, yielding a star formation rate (1400 +- 700 M_\odot yr^{-1}) consistent with that from the infrared luminosity.

CLUMPY AND EXTENDED STARBURSTS IN THE BRIGHTEST UNLENSED SUBMILLIMETER GALAXIES

ABSTRACT The central structure in three of the brightest unlensed z = 3–4 submillimeter galaxies is investigated through 0.″015–0.″05 (120–360 pc) 860 μm continuum images obtained using the Atacama Large Millimeter/submillimeter Array (ALMA). The distribution in the central kiloparsec in AzTEC1 and AzTEC8 is extremely complex, and they are composed of multiple ∼200 pc clumps. AzTEC4 consists of two sources that are separated by ∼1.5 kpc, indicating a mid-stage merger. The peak star formation rate densities in the central clumps are ∼300–3000 M ⊙ yr−1 kpc−2, suggesting regions with extreme star formation near the Eddington limit. By comparing the flux obtained by ALMA and Submillimeter Array, we find that 68%–90% of the emission is extended (≳1 kpc) in AzTEC4 and 8. For AzTEC1, we identify at least 11 additional compact (∼200 pc) clumps in the extended 3–4 kpc region. Overall, the data presented here suggest that the luminosity surface densities observed at ≲150 pc scales are roughly similar to that observed in local ULIRGs, as in the eastern nucleus of Arp 220. Between 10% and 30% of the 860 μm continuum is concentrated in clumpy structures in the central kiloparsec, while the remaining flux is distributed over ≳1 kpc regions, some of which could also be clumpy. These sources can be explained by a rapid inflow of gas such as a merger of gas-rich galaxies, surrounded by extended and clumpy starbursts. However, the cold mode accretion model is not ruled out.

MULTI-WAVELENGTH LENS RECONSTRUCTION OF A PLANCK AND HERSCHEL-DETECTED STAR-BURSTING GALAXY

ABSTRACT We present a source-plane reconstruction of a Herschel and Planck-detected gravitationally lensed dusty star-forming galaxy (DSFG) at z = 1.68 using Hubble, Submillimeter Array (SMA), and Keck observations. The background submillimeter galaxy (SMG) is strongly lensed by a foreground galaxy cluster at z = 0.997 and appears as an arc with a length of ∼15″ in the optical images. The continuum dust emission, as seen by SMA, is limited to a single knot within this arc. We present a lens model with source-plane reconstructions at several wavelengths to show the difference in magnification between the stars and dust, and highlight the importance of multi-wavelength lens models for studies involving lensed DSFGs. We estimate the physical properties of the galaxy by fitting the flux densities to model spectral energy distributions leading to a magnification-corrected star-formation rate (SFR) of 390 ± 60 M yr−1 and a stellar mass of . These values are consistent with high-redshift massive galaxies that have formed most of their stars already. The estimated gas-to-baryon fraction, molecular gas surface density, and SFR surface density have values of 0.43 ± 0.13, 350 ± 200 pc−2, and M yr−1 kpc−2, respectively. The ratio of SFR surface density to molecular gas surface density puts this among the most star-forming systems, similar to other measured SMGs and local ULIRGs.

Water, methanol and dense gas tracers in the local ULIRG Arp 220: results from the new SEPIA Band 5 Science Verification campaign

Abstract We present a line survey of the ultraluminous infrared galaxy Arp 220, taken with the newly installed SEPIA (Swedish-European Southern Observatory PI receiver for APEX) Band 5 instrument on APEX (Atacama Pathfinder Experiment). We illustrate the capacity of SEPIA to detect the 183.3 GHz H2O 31,3–22,0 line against the atmospheric H2O absorption feature. We confirm the previous detection of the HCN(2–1) line, and detect new transitions of standard dense gas tracers such as HNC(2–1), HCO+(2–1), CS(4–3), C34S(4–3) and HC3N(20–19). We also detect HCN(2–1) v2 = 1 and the 193.5 GHz methanol (4–3) group for the first time. The absence of time variations in the megamaser water line compared to previous observations seems to rule out an AGN nuclear origin for the line. It could, on the contrary, favour a thermal origin instead, but also possibly be a sign that the megamaser emission is associated with star-forming cores washed out in the beam. We finally discuss how the new transitions of HCN, HNC and HCO+ refine our knowledge of the interstellar medium physical conditions in Arp 220.

PROPERTIES OF INTERSTELLAR MEDIUM IN INFRARED-BRIGHT QSOs PROBED BY [O i] 63 μm AND [C ii] 158 μm EMISSION LINES*

ABSTRACT We present a study of the interstellar medium (ISM) in the host galaxies of nine QSOs at 0.1 &lt; z &lt; 0.2 with black hole masses of to based on the far-IR spectroscopy taken with Herschel Space Observatory. We detect the [O i] 63 μm ([C ii] 158 μm) emission in 6 (8) out of 8 (9) sources. Our QSO sample has far-infrared luminosities ( ) ∼ several times . The observed line-to- ratios ( / and / ) are in the ranges of 2.6 × 10−4 to 10−2 and 2.8 × 10−4 to 2 × 10−3, respectively (including upper limits). These ratios are comparable to the values found in local ULIRGs, but higher than the average value published so far for IR-bright QSOs. One target, W0752+19, shows an additional broad velocity component (∼720 km s−1) and exceptionally strong [O i] 63 μm emission with / of 10−2, an order of magnitude higher than the average value found among local (U)LIRGs. Combining with the analyses of the Sloan Digital Sky Survey optical spectra, we conclude that the [O i] 63 μm emission in these QSOs is unlikely excited by shocks. We infer that the broad [O i] 63 μm emission in W0752+19 could arise from the warm and dense ISM in the narrow-line region of the central active galactic nucleus. Another possible explanation is the existence of a dense gas outflow with cm−3, where the corresponding broad [C ii] emission is suppressed. Based on the far-IR [O i] and [C ii] line ratios, we estimate constraints on the ISM density and UV radiation field intensity of cm−3 and , respectively. These values are consistent with those found in local Seyfert 1 ULIRGs. In contrast, the gas with broad velocity width in W0752+19 has cm−3 and .

AN ALMA SURVEY OF SUB-MILLIMETER GALAXIES IN THE EXTENDEDCHANDRADEEP FIELD SOUTH: PHYSICAL PROPERTIES DERIVED FROM ULTRAVIOLET-TO-RADIO MODELING

[abridged] The ALESS survey has followed-up a sample of 122 sub-millimeter sources in the Extended Chandra Deep Field South at 870um with ALMA, allowing to pinpoint the positions of sub-millimeter galaxies (SMGs) to 0.3'' and to find their precise counterparts at different wavelengths. This enabled the first compilation of the multi-wavelength spectral energy distributions (SEDs) of a statistically reliable survey of SMGs. In this paper, we present a new calibration of the MAGPHYS modelling code that is optimized to fit these UV-to-radio SEDs of z>1 star-forming galaxies using an energy balance technique to connect the emission from stellar populations, dust attenuation and dust emission in a physically consistent way. We derive statistically and physically robust estimates of the photometric redshifts and physical parameters for the ALESS SMGs. We find that they have a median stellar mass $M_\ast=(8.9\pm0.1)\times10^{10} M_\odot$, SFR$=280\pm70 M_\odot$/yr, overall V-band dust attenuation $A_V=1.9\pm0.2$ mag, dust mass $M_\rm{dust}=(5.6\pm1.0)\times10^8 M_\odot$, and average dust temperature Tdust~40 K. The average intrinsic SED of the ALESS SMGs resembles that of local ULIRGs in the IR range, but the stellar emission of our average SMG is brighter and bluer, indicating lower dust attenuation, possibly because they are more extended. We explore how the average SEDs vary with different parameters, and we provide a new set of SMG templates. To put the ALESS SMGs into context, we compare their stellar masses and SFRs with those of less actively star-forming galaxies at the same redshifts. At z~2, about half of the SMGs lie above the star-forming main sequence, while half are at the high-mass end of the sequence. At higher redshifts (z~3.5), the SMGs tend to have higher SFR and Mstar, but the fraction of SMGs that lie significantly above the main sequence decreases to less than a third.

Physical properties ofz&gt; 4 submillimeter galaxies in the COSMOS field

We study the physical properties of a sample of 6 SMGs in the COSMOS field, spectroscopically confirmed to lie at z>4. We use new GMRT 325 MHz and 3 GHz JVLA data to probe the rest-frame 1.4 GHz emission at z=4, and to estimate the sizes of the star-forming (SF) regions of these sources, resp. Combining our size estimates with those available in the literature for AzTEC1 and AzTEC3 we infer a median radio-emitting size for our z>4 SMGs of (0.63"+/-0.12")x(0.35"+/-0.05") or 4.1x2.3 kpc^2 (major times minor axis; assuming z=4.5) or lower if we take the two marginally resolved SMGs as unresolved. This is consistent with the sizes of SF regions in lower-redshift SMGs, and local normal galaxies, yet higher than the sizes of SF regions of local ULIRGs. Our SMG sample consists of a fair mix of compact and more clumpy systems with multiple, perhaps merging, components. With an average formation time of ~280 Myr, derived through modeling of the UV-IR SEDs, the studied SMGs are young systems. The average stellar mass, dust temperature, and IR luminosity we derive are M*~1.4x10^11 M_sun, T_dust~43 K, and L_IR~1.3x10^13L_sun, resp. The average L_IR is up to an order of magnitude higher than for SMGs at lower redshifts. Our SMGs follow the correlation between dust temperature and IR luminosity as derived for Herschel-selected 0.1<z<2 galaxies. We study the IR-radio correlation for our sources and find a deviation from that derived for z<3 ULIRGs (<q_IR>=1.95+/-0.26 for our sample, compared to q~2.6 for IR luminous galaxies at z<2). In summary, we find that the physical properties derived for our z>4 SMGs put them at the high end of the L_IR-T_dust distribution of SMGs, and that our SMGs form a morphologically heterogeneous sample. Thus, further in-depth analyses of large, statistical samples of high-redshift SMGs are needed to fully understand their role in galaxy formation and evolution.