Probing Radiation Pressure and Hot Gas Feedback through Spectral Simulation of Mid-IR to Submillimeter Fine-Structure Lines in Ultraluminous Infrared Galaxies
Abstract
Recent observational studies have shown that the deeply buried phase found in local ultraluminous infrared galaxies (ULIRGs) is often characterized by powerful feedback thought to be an important mechanism involved in the transformation of gas-rich mergers into gas-poor red and dead galaxies. Based on Herschel studies of multilevel OH transitions in ULIRGs, we have shown that many of the molecular outflows are driven by a combination of an active galactic nucleus (AGN) and a nuclear starburst, as a result of radiation pressure, winds, and supernova remnants. In some sources, however, powerful AGN feedback is required in which the pressure supplied by radiation alone appears insufficient to supply the large outflow momentum fluxes. These outflows appear to be stochastic, strong AGN feedback events that occur throughout the merging process and may involve forces due to hot gas pressure in order to supply the needed momentum fluxes. Numerous theoretical studies have shown that the ratio of the ionizing photon number density to the particle density, commonly denoted as the ionization parameter, is an important diagnostic of the relative dynamical importance of radiation pressure and hot gas pressure in quasar and starburst feedback in galaxies. In optically selected quasars, measured ionization parameters indicate that the current average hot gas pressures are insufficient to power the observed outflows, and optical depths of the AGN radiation are not expected to be high enough to provide the necessary momentum boosts. It has been suggested, however, that during the buried stage of quasar and host galaxy evolution exemplified by ULIRGs, the hot gas pressures could be higher, prior to the development of leaky paths from which the gas can escape, and that infrared fine-structure lines can be used to probe and test this mechanism. The goal of this work is to further our understanding of this phase of galactic evolution by assembling and analyzing the highest quality observations from Spitzer IRS, Herschel PACS, and Herschel SPIRE that heretofore have never been analyzed together in a comprehensive fashion. The publically available Cloudy spectral synthesis code includes photoionization physics, as well as the physics of the transition regions between HII regions and molecular clouds (photodissociation regions (PDRs) and X-ray dominated regions (XDRs) existing beyond the hot ionized centers of AGNs and the cold conditions in molecular clouds. In previous work using Cloudy, we modeled the far-infrared spectra in a sample of the 10 brightest galaxies at 60 microns with infrared luminosities greater than 10^10 Lsol. We showed that both high ionization parameters and high far-infrared optical depths are necessary to explain the far-infrared line deficits found in the only ULIRG in that sample. Here we propose to further develop these models and to apply them for the first time to the full set of available mid-IR, far-IR, and submillimeter spectroscopic observations of 20 ULIRGs and a small comparison sample of LIRGs in the Spitzer and Herschel Archives with multilevel OH observations. With the constraints provided by the numerous mid-IR through submillimeter fine-structure line diagnostics for each galaxy we will evaluate the relative dynamical importance of radiation pressure and hot gas pressure during this pivotal stage of galaxy evolution, and will determine which of these mechanisms is responsible for the momentum fluxes we derived from the multilevel OH Herschel observations of galaxies with powerful molecular outflows.
- Publication:
-
NASA ADAP Proposal
- Pub Date:
- 2017
- Bibcode:
- 2017adap.prop..207F