Massive Star Cluster Formation and Destruction in Luminous Infrared Galaxies in GOALS. II. An ACS/WFC3 Survey of Nearby LIRGs

We present the results of a Hubble Space Telescope WFC3 near-UV and Advanced Camera for Surveys Wide Field Channel optical study into the star cluster populations of a sample of 10 luminous infrared galaxies (LIRGs) in the Great Observatories All-Sky LIRG Survey. Through integrated broadband photometry we have derived ages, masses, and extinctions for a total of 1027 star clusters in galaxies with d L < 110 Mpc in order to avoid issues related to cluster bending. The measured cluster age distribution slope of dN / d τ ∝ τ − 0.5 + / − 0.12 is steeper than what has been observed in lower-luminosity star-forming galaxies. Further, differences in the slope of the observed cluster age distribution between inner- ( dN / d τ ∝ τ − 1.07 + / − 0.12 ) and outer-disk ( dN / d τ ∝ τ − 0.37 + / − 0.09 ) star clusters provide evidence of mass-dependent cluster destruction in the central regions of LIRGs driven primarily by the combined effect of strong tidal shocks and encounters with massive giant molecular clouds. Excluding the nuclear ring surrounding the Seyfert 1 nucleus in NGC 7469, the derived cluster mass function (CMF; dN / dM ∝ M α ) offers marginal evidence for a truncation in the power law at M t ∼ 2×106 M ⊙ for our three most cluster-rich sources, which are all classified as early stage mergers. Finally, we find evidence of a flattening of the CMF slope of dN / dM ∝ M − 1.42 ± 0.1 for clusters in late-stage mergers relative to early stage (α = −1.65 ± 0.02), which we attribute to an increase in the formation of massive clusters over the course of the interaction.

Anomalous Hydrogen Recombination Line Ratios in Ultraluminous Infrared Galaxies

We conducted systematic observations of the H i Brα (4.05 μm) and Brβ (2.63 μm) lines in 52 nearby (z < 0.3) ultraluminous infrared galaxies (ULIRGs) with AKARI. Among 33 ULIRGs wherein the lines are detected, 3 galaxies show anomalous Brβ/Brα line ratios (∼1.0), which are significantly higher than those for case B (0.565). Our observations also show that ULIRGs have a tendency to exhibit higher Brβ/Brα line ratios than those observed in Galactic H ii regions. The high Brβ/Brα line ratios cannot be explained by a combination of dust extinction and case B since dust extinction reduces the ratio. We explore possible causes for the high Brβ/Brα line ratios and show that the observed ratios can be explained by a combination of an optically thick Brα line and an optically thin Brβ line. We simulated the H ii regions in ULIRGs with the Cloudy code, and our results show that the high Brβ/Brα line ratios can be explained by high-density conditions, wherein the Brα line becomes optically thick. To achieve a column density large enough to make the Brα line optically thick within a single H ii region, the gas density must be as high as n ∼ 108 cm−3. We therefore propose an ensemble of H ii regions, in each of which the Brα line is optically thick, to explain the high Brβ/Brα line ratio.

Comprehensive Broadband X-Ray and Multiwavelength Study of Active Galactic Nuclei in 57 Local Luminous and Ultraluminous Infrared Galaxies Observed with NuSTAR and/or Swift/BAT

We perform a systematic X-ray spectroscopic analysis of 57 local luminous and ultraluminous infrared galaxy systems (containing 84 individual galaxies) observed with the Nuclear Spectroscopic Telescope Array and/or Swift/BAT. Combining soft X-ray data obtained with Chandra, XMM-Newton, Suzaku, and/or Swift/XRT, we identify 40 hard (>10 keV) X-ray–detected active galactic nuclei (AGNs) and constrain their torus parameters with the X-ray clumpy torus model XCLUMPY. Among the AGNs at z < 0.03, for which sample biases are minimized, the fraction of Compton-thick (N H ≥ 1024 cm−2) AGNs reaches 64 − 15 + 14 % (6/9 sources) in late mergers, while it is 24 − 10 + 12 % (3/14 sources) in early mergers, consistent with the tendency reported by Ricci et al. We find that the bolometric AGN luminosities derived from the infrared data increase but the X-ray to bolometric luminosity ratios decrease with merger stage. The X-ray-weak AGNs in late mergers ubiquitously show massive outflows at subparsec to kiloparsec scales. Among them, the most luminous AGNs (L bol,AGN ∼ 1046 erg s−1) have relatively small column densities of ≲1023 cm−2 and almost super-Eddington ratios (λ Edd ∼ 1.0). Their torus covering factors (C T (22) ∼ 0.6) are larger than those of Swift/BAT-selected AGNs with similarly high Eddington ratios. These results suggest a scenario where, in the final stage of mergers, multiphase strong outflows are produced due to chaotic quasi-spherical inflows, and the AGN becomes extremely X-ray weak and deeply buried due to obscuration by inflowing and/or outflowing material.

Gas Phase Metallicities of Local Ultra-Luminous Infrared Galaxies Follow Normal Star-Forming Galaxies

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 appears to show that the gas-phase metallicities of local ultra-luminous infrared galaxies (ULIRGs) are significantly lower than those of normal star-forming galaxies1–3. This difference has resulted in the claim that ULIRGs are fueled by metal-poor gas accretion from the outskirts4. Here we report on a new set of gas-phase metallicity measurements making use of the far-infrared spectral lines of [Oiii]52 μm, [Oiii]88 μm, and [Niii]57 μ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 reduces the scatter of the correlation significantly5. 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 ∼ 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.

An ACA Survey of [C i] 3 P 1−3 P 0, CO J = 4 − 3, and Dust Continuum in Nearby U/LIRGs

We present the results of surveying [C i] 3 P 1–3 P 0, 12CO J = 4 − 3, and 630 μm dust continuum emission for 36 nearby ultra/luminous infrared galaxies (U/LIRGs) using the Band 8 receiver mounted on the Atacama Compact Array of the Atacama Large Millimeter/submillimeter Array. We describe the survey, observations, data reduction, and results; the main results are as follows. (i) We confirmed that [C i] 3 P 1–3 P 0 has a linear relationship with both the 12CO J = 4 − 3 and 630 μm continuum. (ii) In NGC 6052 and NGC 7679, 12CO J = 4 − 3 was detected but [C i] 3 P 1–3 P 0 was not detected with a [C i] 3 P 1–3 P 0/12CO J = 4 − 3 ratio of ≲0.08. Two possible scenarios of weak [C i] 3 P 1–3 P 0 emission are C0-poor/CO-rich environments and an environment with an extremely large [C i] 3 P 1–3 P 0 missing flux. (iii) There is no clear evidence showing that galaxy mergers, AGNs, and dust temperatures control the ratios of [C i] 3 P 1–3 P 0/12CO J = 4 − 3 and L [ C I ] ( 1 − 0 ) ′ / L 630 μ m . (iv) We compare our nearby U/LIRGs with high-z galaxies, such as galaxies on the star formation main sequence (MS) at z ∼ 1 and submillimeter galaxies (SMGs) at z = 2–4. We found that the mean value for the [C i] 3 P 1–3 P 0/12CO J = 4 − 3 ratio of U/LIRGs is similar to that of SMGs but smaller than that of galaxies on the MS.