Local Ultraluminous Infrared Galaxies Research Articles

A possible relation between global CO excitation and massive molecular outflows in local ULIRGs

Local ultra-luminous infrared galaxies (ULIRGs) have been observed to host ubiquitous molecular outflows, including the most massive and powerful ever detected. These sources have also exceptionally excited global, galaxy-integrated CO ladders. A connection between outflows and molecular gas excitation has however never been established, since previous multi-J CO surveys were limited in spectral resolution and sensitivity and so they could only probe the global molecular gas conditions. In this work, we address this question using new, ground-based, sensitive heterodyne spectroscopy of multiple CO rotational lines (up to CO(7−6)) in a sample of 17 local ULIRGs. We used the Atacama Pathfinder Experiment (APEX) telescope to survey the CO(Jup ≥ 4) lines at a high signal-to-noise ratio, and complemented these data with CO(Jup ≤ 3) APEX and Atacama Large Millimeter Array (ALMA and ACA) observations presented in Montoya Arroyave et al. (2023, A&A, 673, A13). We detected a total of 74 (out of 75) CO lines, with up to six transitions per source. The resulting CO spectral line energy distributions (SLEDs) show a wide range in gas excitation, in agreement with previous studies on ULIRGs. Some CO SLEDs peak at Jup ∼ 3, 4, which we classify as “lower excitation”, while others plateau or keep increasing up to the highest-J CO transition probed, and we classify these as “higher excitation”. Our analysis includes for completeness the results of CO SLED fits performed with a single large velocity gradient component, but our main focus is the investigation of possible links between global CO excitation and the presence of broad and/or high-velocity CO spectral components that can contain outflowing gas. We discovered an increasing trend of line width as a function of Jup of the CO transition, which is significant at the 4σ level and appears to be driven by the eight sources that we classified as higher excitation. We further analyzed such higher-excitation ULIRGs, by performing a decomposition of their CO spectral profiles into multiple components, and we derived CO ladders that are clearly more excited for the spectral components characterized by higher velocities and/or velocity dispersion. Because these sources are known to host widespread molecular outflows, we favor an interpretation whereby the highly excited CO-emitting gas in ULIRGs resides in galactic-scale massive molecular outflows whose emission fills a large fraction of the beam of our APEX high-J CO observations. On the other hand, our results challenge alternative scenarios for which the high CO excitation in ULIRGs can be explained by classical component of the interstellar medium, such as photon- or X-ray dominated regions around the nuclear sources.

SMART: Spectral energy distributions Markov chain Analysis with Radiative Transfer models

Abstract In this paper we present the publicly available open-source spectral energy distribution (SED) fitting code SMART (Spectral energy distributions Markov chain Analysis with Radiative Transfer models). Implementing a Bayesian Markov chain Monte Carlo (MCMC) method, SMART fits the ultraviolet to millimetre SEDs of galaxies exclusively with radiative transfer models that currently constitute four types of pre-computed libraries, which describe the starburst, active galactic nucleus (AGN) torus, host galaxy and polar dust components. An important novelty of SMART is that, although it fits SEDs exclusively with radiative transfer models, it takes comparable time to popular energy balance methods to run. Here we describe the key features of SMART and test it by fitting the multi-wavelength SEDs of the 42 local ultraluminous infrared galaxies (ULIRGs) that constitute the HERschel Ultraluminous Infrared Galaxy Survey (HERUS) sample. The Spitzer spectroscopy data of the HERUS ULIRGs are included in the fitting at a spectral resolution, which is matched to that of the radiative transfer models. We also present other results that highlight the performance and versatility of SMART. SMART promises to be a useful tool for studying galaxy evolution in the James Webb Space Telescope (JWST) era. SMART is developed in Python and is available at https://github.com/ch-var/SMART.git.

Radio-only and Radio-to-far-ultraviolet Spectral Energy Distribution Modeling of 14 ULIRGs: Insights into the Global Properties of Infrared Bright Galaxies

We present detailed spectral energy distribution (SED) modeling of 14 local ultraluminous infrared galaxies (ULIRGs) with outstanding photometric data from the literature covering the ultraviolet–infrared (FIR) and radio bands (∼50 MHz to ∼30 GHz). We employ the CIGALE SED fitting code to model the ultraviolet–FIR–radio SED. For the radio-only SED modeling, we use the UltraNest package, leveraging its nested sampling algorithm. Combining the results from our previous study on 11 luminous infrared galaxies (LIRGs), we discuss the global astrophysical properties of a sample of 25 starburst galaxies (z < 0.5). Their radio spectra are frequently characterized by bends and turnovers, with no indication of ULIRGs exhibiting more complicated SEDs than LIRGs despite showing more signs of interactions. Including radio measurements in the CIGALE modeling constrained the dust luminosity and star formation rate (SFR) estimates by more than 1 order of magnitude better than previously reported for starburst galaxies. We show that total and nonthermal radio luminosity at 1.4 and 4.8 GHz frequencies can be good estimators of recent SFRs for all LIRGs and those ULIRGS with an insignificant influence of active galactic nuclei. A weaker but still significant correlation is observed between radio SFRs at 1.4 GHz and old (averaged over 100 Myr) SFRs based on SED modeling, indicative of multiple episodes of starburst activity during their lifetime. The thermal radio luminosity at 4.8 GHz is a better tracer of recent star formation than the thermal luminosity at 1.4 GHz. Statistically, our modeled nonthermal radio spectral indices do not significantly correlate with redshift, stellar mass, SFR, specific SFR, and dust mass.

The first ground-based detection of the 752 GHz water line in local ultra-luminous infrared galaxies using APEX-SEPIA

Abstract We report the first ground-based detection of the water line p-H2O (211–202) at 752.033 GHz in three z &amp;lt; 0.08 ultra-luminous infrared galaxies (ULIRGs): IRAS 06035-7102, IRAS 17207-0014 and IRAS 09022-3615. Using the Atacama Pathfinder EXperiment (APEX), with its Swedish-ESO PI Instrument for APEX (SEPIA) band-9 receiver, we detect this H2O line with overall signal-to-noise ratios of 8–10 in all three galaxies. Notably, this is the first detection of this line in IRAS 06035-7102. Our new APEX-measured fluxes, between 145 to 705 Jy km s−1, are compared with previous values taken from Herschel SPIRE FTS. We highlight the great capabilities of APEX for resolving the H2O line profiles with high spectral resolutions while also improving by a factor of two the significance of the detection within moderate integration times. While exploring the correlation between the p-H2O(211–202) and the total infrared luminosity, our galaxies are found to follow the trend at the bright end of the local ULIRG’s distribution. The p-H2O(211–202) line spectra are compared to the mid-J CO and HCN spectra, and dust continuum previously observed with ALMA. In the complex interacting system IRAS 09022-3615, the profile of the water emission line is offset in velocity with respect to the ALMA CO(J = 4 − 3) emission. For IRAS 17207-0014 and IRAS 06035-7102, the profiles between the water line and the CO lines are spectroscopically aligned. This pilot study demonstrates the feasibility of directly conducting ground-based high-frequency observations of this key water line, opening the possibility of detailed follow-up campaigns to tackle its nature.

Graph Theoretical Analysis of local ultraluminous infrared galaxies and quasars

We present a methodological framework for studying galaxy evolution by utilizing Graph Theory and network analysis tools. We study the evolutionary processes of local ultraluminous infrared galaxies (ULIRGs) and quasars and the underlying physical processes, such as star formation and active galactic nucleus (AGN) activity, through the application of Graph Theoretical analysis tools. We extract, process and analyze mid-infrared spectra of local (z ¡ 0.4) ULIRGs and quasars between 5-38μm through internally developed Python routines, in order to generate similarity graphs, with the nodes representing ULIRGs being grouped together based on the similarity of their spectra. Additionally, we extract and compare physical features from the mid-IR spectra, such as the polycyclic aromatic hydrocarbons (PAHs) emission and silicate depth absorption features, as indicators of the presence of star-forming regions and obscuring dust, in order to understand the underlying physical mechanisms of each evolutionary stage of ULIRGs. Our analysis identifies five groups of local ULIRGs based on their mid-IR spectra, which is quite consistent with the well established fork classification diagram by providing a higher level classification. We demonstrate how graph clustering algorithms and network analysis tools can be utilized as unsupervised learning techniques for revealing direct or indirect relations between various galaxy properties and evolutionary stages, which provides an alternative methodology to previous works for classification in galaxy evolution. Additionally, our methodology compares the output of several graph clustering algorithms in order to demonstrate the best-performing Graph Theoretical tools for studying galaxy evolution.

Classification of local ultraluminous infrared galaxies and quasars with kernel principal component analysis

ABSTRACT We present a new diagnostic diagram for local ultraluminous infrared galaxies (ULIRGs) and quasars, analysing particularly the Spitzer Space Telescope’s infrared spectrograph spectra of 102 local ULIRGs and 37 Palomar Green quasars. Our diagram is based on a special non-linear mapping of these data, employing the kernel principal component analysis method. The novelty of this map lies in the fact that it distributes the galaxies under study on the surface of a well-defined ellipsoid, which, in turn, links basic concepts from geometry to physical properties of the galaxies. Particularly, we have found that the equatorial direction of the ellipsoid corresponds to the evolution of the power source of ULIRGs, starting from the pre-merger phase, moving through the starburst-dominated coalescing stage towards the active galactic nucleus-dominated phase, and finally terminating with the post-merger quasar phase. On the other hand, the meridian directions distinguish deeply obscured power sources of the galaxies from unobscured ones. These observations have also been verified by comparison with simulated ULIRGs and quasars using radiative transfer models. The diagram correctly identifies unique galaxies with extreme features that lie distinctly away from the main distribution of the galaxies. Furthermore, special two-dimensional projections of the ellipsoid recover almost monotonic variations of the two main physical properties of the galaxies, the silicate and polycyclic aromatic hydrocarbon features. This suggests that our diagram naturally extends the well-known Spoon diagram and it can serve as a diagnostic tool for existing and future infrared spectroscopic data, such as those provided by the James Webb Space Telescope.

Low gas-phase metallicities of ultraluminous infrared galaxies are a result of dust obscuration

Despite advances in observational data, theoretical models, and computational techniques to simulate key physical processes in the formation and evolution of galaxies, the stellar mass assembly of galaxies still remains an unsolved problem today. Optical spectroscopic measurements appear to show that the gas-phase metallicities of local ultra-luminous infrared galaxies (ULIRGs) are significantly lower than those of normal star-forming galaxies. This difference has resulted in the claim that ULIRGs are fueled by metal-poor gas accretion from the outskirts\cite{Mannucci10}. Here we report on a new set of gas-phase metallicity measurements making use of the far-infrared spectral lines of [O{\sc iii}]52 $\mu$m, [O{\sc iii}]88 $\mu$m, and [N{\sc iii}]57 $\mu$m instead of the usual optical lines. Photoionization models have resulted in a metallicity diagnostic based on these three lines that break the electron density degeneracy and reduce the scatter of the correlation significantly. Using new data from SOFIA and archival data from Herschel Space Observatory, we find that local ULIRGs lie on the mass-metallicity relation of star-forming galaxies and have metallicities comparable to other galaxies with similar stellar masses and star formation rates. The lack of a departure suggests that ULIRGs follow the same mass assembly mechanism as luminous star-forming galaxies and $\sim 0.3$ dex under-abundance in metallicities derived from optical lines is a result of heavily obscured metal-rich gas which has a negligible effect when using the FIR line diagnostics.

Quantifying the cool ISM in radio AGNs: evidence for late-time retriggering by galaxy mergers and interactions

ABSTRACTWe use deep Herschel observations of the complete 2Jy sample of powerful radio active galactic nuclei (AGNs) in the local Universe (0.05 &amp;lt; z &amp;lt; 0.7) to probe their cool interstellar medium (ISM) contents and star-forming properties, comparing them against other samples of nearby luminous AGNs and quiescent galaxies. This allows us to investigate triggering and feedback mechanisms. We find that the dust masses of the strong-line radio galaxies (SLRGs) in our sample are similar to those of radio-quiet quasars, and that their median dust mass (Mdust = 2 × 107 M⊙) is enhanced by a factor of ∼200 compared to that of non-AGN ellipticals, but lower by a factor of ∼16 relative to that of local ultraluminous infrared galaxies (ULIRGs). Along with compelling evidence for merger signatures in optical images, the SLRGs in our sample also show relatively high star formation efficiencies, despite the fact that many of them fall below the main sequence for star-forming galaxies. Together, these results suggest that most of our SLRGs have been retriggered by late-time mergers that are relatively minor in terms of their gas contents. In comparison with the SLRGs, the radio AGNs with weak optical emission lines (weak-line radio galaxies – WLRGs) and edge-darkened radio jets (Fanaroff–Riley Class I radio sources – FRIs) have both lower cool ISM masses and star formation rates (by a factor of &amp;gt;30), consistent with being fuelled by a different mechanism (e.g. the direct accretion of hot gas).

A new look at local ultraluminous infrared galaxies: the atlas and radiative transfer models of their complex physics

ABSTRACT We present the ultraviolet to submillimetre spectral energy distributions (SEDs) of the HERschel Ultra Luminous Infrared Galaxy Survey (HERUS) sample of 42 local ultraluminous infrared galaxies (ULIRGs) and fit them with a Markov chain Monte Carlo code using the CYprus models for Galaxies and their NUclear Spectra (CYGNUS) radiative transfer models for starbursts, active galactic nucleus (AGN) tori, and host galaxy. The Spitzer IRS spectroscopy data are included in the fitting. Our Bayesian SED fitting method takes comparable time to popular energy balance methods but it is more physically motivated and versatile. All HERUS galaxies harbour high rates of star formation but we also find bolometrically significant AGN in all of the galaxies of the sample. We estimate the correction of the luminosities of the AGN in the ULIRGs due to the anisotropic emission of the torus and find that it could be up to a factor of ∼10 for nearly edge-on tori. We present a comparison of our results with the smooth torus model of Fritz et al. and the two-phase models of Siebenmorgen et al. and SKIRTOR. We find that the CYGNUS AGN torus models fit significantly better the SEDs of our sample compared to all other models. We find no evidence that strong AGN appear either at the beginning or end of a starburst episode or that starbursts and AGN affect each other. IRAS 01003−2238 and Mrk 1014 show evidence of dual AGNs in their SED fits, suggesting a minimum dual AGN fraction in the sample of 5 per cent.

Discovery of a luminous starburst galaxy with hundreds of thousands of Wolf–Rayet stars

ABSTRACT This paper reports the detection of a strong Wolf–Rayet (W-R) emission bump feature at 4400–4700 Å in the optical spectrum of a distant galaxy SDSS J150009.81+452844.4 at a redshift of 0.453. The estimated number of W-R stars is 1.5 × 105 without applying any internal extinction correction, and at least 4.5 × 105 after correcting for dust extinction. Such a number, though with uncertainties inherent from the extinction correction and others, appears to outnumber those in previous W-R galaxies after correcting for intrinsic dust extinction with the numbers available in the literature. These massive stars must have formed in an instantaneous star-forming episode lasting less than about 5 Myr. We estimate a star formation rate of at least ∼80 $\mathrm{M}_{\odot }\, \mathrm{yr^{-1}}$, indicative of a violent starburst. Its mid-infrared spectrum resembles closely those of typical local ultraluminous infrared galaxies (ULIRGs). From optical to ultraviolet, it possesses a luminous blue continuum, indicating that the starburst is not heavily obscured, in contrast to that found in most ULIRGs. There is evidence for strong outflows based on the detection of systematically blue-shifted broad wings in the nebular emission lines (a bulk velocity −190 km s−1), as well as in the blue-shifted Mg ii absorption lines. The emission-line gases show a wide velocity range, from −1200 to ∼−2000 km s−1 in blueshift to 570 km s−1 in redshift. We interpret this as a galactic fountain, of which part of the outflowing gas may be falling back to the galactic plane. Our method could be used to find a sample of similar objects, which would help understand the star formation history and stellar feedback in starburst galaxies.

The physics of gas phase metallicity gradients in galaxies

ABSTRACT We present a new model for the evolution of gas phase metallicity gradients in galaxies from first principles. We show that metallicity gradients depend on four ratios that collectively describe the metal equilibration time-scale, production, transport, consumption, and loss. Our model finds that most galaxy metallicity gradients are in equilibrium at all redshifts. When normalized by metal diffusion, metallicity gradients are governed by the competition between radial advection, metal production, and accretion of metal-poor gas from the cosmic web. The model naturally explains the varying gradients measured in local spirals, local dwarfs, and high-redshift star-forming galaxies. We use the model to study the cosmic evolution of gradients across redshift, showing that the gradient in Milky Way-like galaxies has steepened over time, in good agreement with both observations and simulations. We also predict the evolution of metallicity gradients with redshift in galaxy samples constructed using both matched stellar masses and matched abundances. Our model shows that massive galaxies transition from the advection-dominated to the accretion-dominated regime from high to low redshifts, which mirrors the transition from gravity-driven to star formation feedback-driven turbulence. Lastly, we show that gradients in local ultraluminous infrared galaxies (major mergers) and inverted gradients seen both in the local and high-redshift galaxies may not be in equilibrium. In subsequent papers in this series, we show that the model also explains the observed relationship between galaxy mass and metallicity gradients, and between metallicity gradients and galaxy kinematics.

Planck’s Dusty GEMS

We present ALMA, NOEMA, and IRAM-30 m/EMIR observations of the high-density tracer molecules HCN, HCO+, and HNC in three of the brightest lensed dusty star-forming galaxies at z ≃ 3–3.5, part of the Planck’s Dusty Gravitationally Enhanced subMillimetre Sources (GEMS), with the aim of probing the gas reservoirs closely associated with their exceptional levels of star formation. We obtained robust detections of ten emission lines between Jup = 4 and 6, as well as several additional upper flux limits. In PLCK_G244.8+54.9, the brightest source at z = 3.0, the HNC(5–4) line emission at 0.1″ resolution, together with other spatially-integrated line profiles, suggests comparable distributions of dense and more diffuse gas reservoirs, at least over the most strongly magnified regions. This rules out any major effect from differential lensing. This line is blended with CN(4–3) and in this source, we measure a HNC(5–4)/CN(4–3) flux ratio of 1.76 ±0. 86. Dense-gas line profiles generally match those of mid-J CO lines, except in PLCK_G145.2+50.8, which also has dense-gas line fluxes that are relatively lower, perhaps due to fewer dense cores and more segregated dense and diffuse gas phases in this source. The HCO+/HCN ≳ 1 and HNC/HCN ∼ 1 line ratios in our sample are similar to those of nearby ultraluminous infrared galaxies (ULIRGs) and consistent with photon-dominated regions without any indication of important mechanical heating or active galactic nuclei feedback. We characterize the dense-gas excitation in PLCK_G244.8+54.9 using radiative transfer models assuming pure collisional excitation and find that mid-J HCN, HCO+, and HNC lines arise from a high-density phase with an H2 density of n ∼ 105–106 cm−3, although important degeneracies hinder a determination of the exact conditions. The three GEMS are consistent with extrapolations of dense-gas star-formation laws derived in the nearby Universe, adding further evidence that the extreme star-formation rates observed in the most active galaxies at z ∼ 3 are a consequence of their important dense-gas contents. The dense-gas-mass fractions traced by HCN/[CI] and HCO+/[CI] line ratios are elevated, but not exceptional as compared to other lensed dusty star-forming galaxies at z &gt; 2, and they fall near the upper envelope of local ULIRGs. Despite the higher overall gas fractions and local gas-mass surface densities observed at high redshift, the dense-gas budget of rapidly star-forming galaxies seems to have evolved little between z ∼ 3 and z ∼ 0. Our results favor constant dense-gas depletion times in these populations, which is in agreement with theoretical models of star formation.

Kiloparsec-scale ALMA Imaging of [C ii] and Dust Continuum Emission of 27 Quasar Host Galaxies at z ∼ 6

Abstract We present a study of the [C ii] 158 μm line and underlying far-infrared (FIR) continuum emission of 27 quasar host galaxies at z ∼ 6, traced by the Atacama Large Millimeter/submillimeter Array at a spatial resolution of ∼1 physical kpc. The [C ii] emission in the bright, central regions of the quasars have sizes of 1.0–4.8 kpc. The dust continuum emission is typically more compact than [C ii]. We find that 13/27 quasars (approximately one-half) have companion galaxies in the field, at projected separations of 3–90 kpc. The position of dust emission and the Gaia-corrected positions of the central accreting black holes are cospatial (typical offsets ≲0.″1). This suggests that the central black holes are located at the bottom of the gravitational wells of the dark matter halos in which the quasar hosts reside. Some outliers with offsets of ∼500 pc can be linked to disturbed morphologies, most likely due to ongoing or recent mergers. We find no correlation between the central brightness of the FIR emission and the bolometric luminosity of the accreting black hole. The FIR-derived star formation rate densities (SFRDs) in the host galaxies peak at the galaxies’ centers, at typical values between 100 and 1000 M ⊙ yr−1 kpc−2. These values are below the Eddington limit for star formation, but similar to those found in local ultraluminous infrared galaxies. The SFRDs drop toward larger radii by an order of magnitude. Likewise, the [C ii]/FIR luminosity ratios of the quasar hosts are lowest in their centers (few ×10−4) and increase by a factor of a few toward the galaxies’ outskirts, consistent with resolved studies of lower-redshift sources.

An ALMA Survey of the SCUBA-2 Cosmology Legacy Survey UKIDSS/UDS Field: The Far-infrared/Radio Correlation for High-redshift Dusty Star-forming Galaxies

Abstract We study the radio properties of 706 submillimeter galaxies (SMGs) selected at 870 μm with the Atacama Large Millimeter Array from the SCUBA-2 Cosmology Legacy Survey map of the Ultra Deep Survey field. We detect 273 SMGs at &gt;4σ in deep Karl G. Jansky Very Large Array 1.4 GHz observations, of which a subset of 45 SMGs are additionally detected in 610 MHz Giant Metre-Wave Radio Telescope imaging. We quantify the far-infrared/radio correlation (FIRRC) through parameter q IR, defined as the logarithmic ratio of the far-infrared and radio luminosity, and include the radio-undetected SMGs through a stacking analysis. We determine a median q IR = 2.20 ± 0.03 for the full sample, independent of redshift, which places these z ∼ 2.5 dusty star-forming galaxies 0.44 ± 0.04 dex below the local correlation for both normal star-forming galaxies and local ultra-luminous infrared galaxies (ULIRGs). Both the lack of redshift evolution and the offset from the local correlation are likely the result of the different physical conditions in high-redshift starburst galaxies, compared to local star-forming sources. We explain the offset through a combination of strong magnetic fields (B ≳ 0.2 mG), high interstellar medium (ISM) densities and additional radio emission generated by secondary cosmic rays. While local ULIRGs are likely to have similar magnetic field strengths, we find that their compactness, in combination with a higher ISM density compared to SMGs, naturally explains why local and high-redshift dusty star-forming galaxies follow a different FIRRC. Overall, our findings paint SMGs as a homogeneous population of galaxies, as illustrated by their tight and nonevolving far-infrared/radio correlation.

A systematic study of ULIRGs using near-infrared absorption bands reveals a strong UV environment in their star-forming regions

Abstract We present a systematic study of the 3.0 μm H2O ice and the 3.4 μm aliphatic carbon absorption features toward 48 local ultraluminous infrared galaxies (ULIRGs) using spectra obtained by the AKARI Infrared Camera to investigate the UV environment in their star-forming regions. All the ULIRGs in our sample exhibit a ratio of optical depth of H2O ice to silicate dust (τ3.0/τ9.7) that is lower than that in the Taurus dark cloud. This implies that ULIRGs cannot be described as an ensemble of low-mass star-forming regions and that a significant amount of high-mass star-forming regions contribute to star-forming clouds in local ULIRGs. The results also show that the ratios of optical depth of aliphatic carbon to silicate dust, τ3.4/τ9.7, exhibit diverse values. We investigate two effects that can affect this ratio: the geometric temperature gradient (which increases the ratio) and the intense UV environment (which decreases it). The geometric temperature gradient is typically considered as a sign of active galactic nuclei (AGN). ULIRGs with AGN signs (optical classification, near-infrared color, and a polycyclic aromatic hydrocarbon emission strength of $3.3\, \mu$m) indeed tend to exhibit a large τ3.4/τ9.7 ratio. However, we find that the presence of buried AGN is not the only cause of the geometric temperature gradient, because the enhancement of the ratio is also evident in pure starburst-like ULIRGs without these AGN signs. Regarding the intense UV environment in star-forming regions, the correlation between the aliphatic carbon ratio and the ratio of the [C ii] 158 μm line luminosity to the far-infrared luminosity ($L_{[\mathrm{C}\,\,\small {ii}]}/L_{\mathrm{FIR}}$), which represents the UV environment in photodissociation regions, implies that the intense UV environment causes the decrease of the aliphatic carbon ratio. We find that an intense UV environment (G/nH &amp;gt; 3) in star-forming regions is needed for the aliphatic carbon ratio to be suppressed.

Deciphering an evolutionary sequence of merger stages in infrared-luminous starburst galaxies atz∼ 0.7

Based on optical and near-IRMagellanFIRE spectra of 25 starburst galaxies at 0.5 &lt; z &lt; 0.9, a recent publication showed that their attenuation properties can be explained by a single-parameter sequence of total obscurations ranging fromAV = 2 toAV = 30 toward the starburst core centers in a mixed stars and dust configuration. We investigate here the origin of this sequence for the same sample. We show that total attenuations anticorrelate with the starburst sizes in radio (3 GHz) with a significance larger than 5σand a scatter of 0.26 dex. More obscured and compact starbursts also show enhanced N2 (=[NII]/Hα) ratios and larger line velocity widths that we attribute to an increasing shock contribution toward later merger phases, driven by deeper gravitational potential wells at the coalescence. Additionally, the attenuation is also linked to the equivalent width (EW) of hydrogen recombination lines, which is sensitive to the luminosity weighted age of the relatively unobscured stellar populations. Overall, the correlations amongAV, tot, radio size, line width, N2 and EW of Balmer and Paschen lines converge toward suggesting an evolutionary sequence of merger stages: all of these quantities are likely to be good time-tracers of the merger phenomenon, and their large spanned range appears to be characteristic of the different merger phases. Half of our sample at higher obscurations have radio sizes approximately 3 times smaller than early type galaxies at the same redshift, suggesting that, in analogy with local ultraluminous infrared galaxies (ULIRGs), these cores cannot be directly forming elliptical galaxies. Finally, we detect mid-IR AGN torus for half of our sample and additional X-ray emission for 6 starbursts; intriguingly, the latter have systematically more compact sizes, suggestive of emerging AGNs toward later merger stages, possibly precursors of a later QSO phase.

What determines the shape of the local (z &amp;lt; 0.1) infrared galaxy luminosity function?

ABSTRACT We investigate what shapes the infrared luminosity function of local galaxies by comparing it to the local infrared active galactic nucleus (AGN) luminosity function. The former corresponds to emission from dust heated by stars and AGN, whereas the latter includes emission from AGN-heated dust only. Our results show that infrared emission from AGN starts mixing into the galaxy luminosity function in the luminous infrared galaxy (LIRG) regime and becomes significant in the ultraluminous infrared galaxy (ULIRG) regime, with the luminosity above which local ULIRGs become AGN-dominated being in the log$L_{\rm IR}/\rm L_{\odot }\sim 12.2$–12.7 range. We propose that as a result of the AGN contribution, the infrared galaxy luminosity function has a flatter high-luminosity slope than UV/optical galaxy luminosity functions. Furthermore, we note that the increased AGN contribution as a function of LIR is reflected in the average dust temperature (Tdust) of local galaxies, and may be responsible for the local LIR–Tdust relation. However, although our results show that AGN play a central role in defining the properties of local ULIRGs, we find that the dominant power source in the local ULIRG population is star formation.

The Dual Role of Starbursts and Active Galactic Nuclei in Driving Extreme Molecular Outflows

Abstract We report molecular gas observations of IRAS 20100−4156 and IRAS 03158+4227, two local ultraluminous infrared galaxies (ULIRGs) hosting some of the fastest and most massive molecular outflows known. Using Atacama Large Millimeter Array and Plateau de Bure Interferometer observations, we spatially resolve the CO (1−0) emission from the outflowing molecular gas in both and find maximum outflow velocities of v max ∼ 1600 and ∼1700 km s−1 for IRAS 20100−4156 and IRAS 03158+4227, respectively. We find total gas mass outflow rates of and ∼350 M ⊙ yr−1, respectively, corresponding to molecular gas depletion timescales and ∼16 Myr. This is nearly 3 times shorter than the depletion timescales implied by star formation, and ∼46 Myr, respectively. To determine the outflow driving mechanism, we compare the starburst luminosity (L *) and active galactic nucleus (AGN) luminosity (L AGN) to the outflowing energy and momentum fluxes, using mid-infrared spectral decomposition to discern L AGN. Comparison to other molecular outflows in ULIRGs reveals that outflow properties correlate similarly with L * and L IR as with L AGN, indicating that AGN luminosity alone may not be a good tracer of feedback strength and that a combination of AGN and starburst activity may be driving the most powerful molecular outflows. We also detect the OH 1.667 GHz maser line from both sources and demonstrate its utility in detecting molecular outflows.

HERUS: the far-IR/submm spectral energy distributions of local ULIRGs and photometric atlas

We present the Herschel-SPIRE photometric atlas for a complete flux limited sample of 43 local ultraluminous infrared galaxies (ULIRGs), selected at 60 μm by IRAS, as part of the HERschel ULIRG Survey (HERUS). Photometry observations were obtained using the SPIRE instrument at 250, 350, and 500 μm. We describe these observations, present the results, and combine the new observations with data from IRAS to examine the far-infrared spectral energy distributions (SEDs) of these sources. We fit the observed SEDs of HERUS objects with a simple parametrized modified blackbody model, where temperature and emissivity β are free parameters. We compare the fitted values to those of non-ULIRG local galaxies, and find, in agreement with earlier results, that HERUS ULIRGs have warmer dust (median temperature T = 37.9 ± 4.7 K compared to 21.3 ± 3.4 K) but a similar β distribution (median β = 1.7 compared to 1.8) to the Herschel reference sample (HRS, Cortese et al. 2014) galaxies. Dust masses are found to be in the range of 107.5–109 M⊙, significantly higher than that of HRS sources. We compare our results for local ULIRGs with higher redshift samples selected at 250 and 850 μm. These latter sources generally have cooler dust and/or redder 100-to-250 μm colours than our 60 μm-selected ULIRGs. We show that this difference may in part be the result of the sources being selected at different wavelengths rather than being a simple indication of rapid evolution in the properties of the population.

A dusty star-forming galaxy at z = 6 revealed by strong gravitational lensing

Since their discovery, submillimetre-selected galaxies (SMGs) have revolutionized the field of galaxy formation and evolution. From the hundreds of square degrees mapped at submillimetre wavelengths, only a handful of sources have been confirmed to lie at z>5 and only two at z>6. All of these SMGs are rare examples of extreme starburst galaxies with star formation rates (SFRs) of >1000 M_sun/yr and therefore are not representative of the general population of dusty star-forming galaxies. Consequently, our understanding of the nature of these sources, at the earliest epochs, is still incomplete. Here we report the spectroscopic identification of a gravitationally amplified (mu = 9.3 +/- 1.0) dusty star-forming galaxy at z=6.027. After correcting for gravitational lensing, we derive an intrinsic less-extreme SFR of 380 +/- 50 M_sun/yr for this source and find that its gas and dust properties are similar to those measured for local Ultra Luminous Infrared Galaxies (ULIRGs), extending the local trends to a poorly explored territory in the early Universe. The star-formation efficiency of this galaxy is similar to those measured in its local analogues, despite a ~12 Gyr difference in cosmic time.