We have mapped the superwind/halo region of the nearby starburst galaxy M82 in the mid-infrared with $Spitzer-IRS$. The spectral regions covered include the H$_2 S(1)-S(3)$, [NeII], [NeIII] emission lines and PAH features. We estimate the total warm
H$_2$ mass and the kinetic energy of the outflowing warm molecular gas to be between $M_{warm}sim5-17times10^6$ M$_{odot}$ and $E_{K}sim6-20times10^{53}$ erg. Using the ratios of the 6.2, 7.7 and 11.3 micron PAH features in the IRS spectra, we are able to estimate the average size and ionization state of the small grains in the superwind. There are large variations in the PAH flux ratios throughout the outflow. The 11.3/7.7 and the 6.2/7.7 PAH ratios both vary by more than a factor of five across the wind region. The Northern part of the wind has a significant population of PAHs with smaller 6.2/7.7 ratios than either the starburst disk or the Southern wind, indicating that on average, PAH emitters are larger and more ionized. The warm molecular gas to PAH flux ratios (H$_2/PAH$) are enhanced in the outflow by factors of 10-100 as compared to the starburst disk. This enhancement in the H$_2/PAH$ ratio does not seem to follow the ionization of the atomic gas (as measured with the [NeIII]/[NeII] line flux ratio) in the outflow. This suggests that much of the warm H$_2$ in the outflow is excited by shocks. The observed H$_2$ line intensities can be reproduced with low velocity shocks ($v < 40$ km s$^{-1}$) driven into moderately dense molecular gas ($10^2 <n_H < 10^4$ cm$^{-3}$) entrained in the outflow.
The manner in which gas accretes and orbits within circumnuclear rings has direct implications for the star formation process. In particular, gas may be compressed and shocked at the inflow points, resulting in bursts of star formation at these locat
ions. Afterwards the gas and young stars move together through the ring. In addition, star formation may occur throughout the ring, if and when the gas reaches sufficient density to collapse under gravity. These two scenarios for star formation in rings are often referred to as the `pearls on a string and `popcorn paradigms. In this paper, we use new Herschel PACS observations, obtained as part of the KINGFISH Open Time Key Program, along with archival Spitzer and ground-based observations from the SINGS Legacy project, to investigate the heating and cooling of the interstellar medium in the nearby star-forming ring galaxy, NGC4736. By comparing spatially resolved estimates of the stellar FUV flux available for heating, with the gas and dust cooling derived from the FIR continuum and line emission, we show that while star formation is indeed dominant at the inflow points in NGC 4736, additional star formation is needed to balance the gas heating and cooling throughout the ring. This additional component most likely arises from the general increase in gas density in the ring over its lifetime. Our data provide strong evidence, therefore, for a combination of the two paradigms for star formation in the ring in NGC4736.
We present Herschel/PACS observations of extended [CII]157.7{mu}m line emission detected on ~ 1 - 10 kpc scales in 60 local luminous infrared galaxies (LIRGs) from the Great Observatories All-sky LIRG Survey (GOALS). We find that most of the extra-nu
clear emission show [CII]/FIR ratios >~ 4 x 10^-3, larger than the mean ratio seen in the nuclei, and similar to those found in the extended disks of normal star-forming galaxies and the diffuse inter-stellar medium (ISM) of our Galaxy. The [CII] deficits found in the most luminous local LIRGs are therefore restricted to their nuclei. There is a trend for LIRGs with warmer nuclei to show larger differences between their nuclear and extra-nuclear [CII]/FIR ratios. We find an anti-correlation between [CII]/FIR and the luminosity surface density, {Sigma}_IR, for the extended emission in the spatially-resolved galaxies. However, there is an offset between this trend and that found for the LIRG nuclei. We use this offset to derive a beam filling-factor for the star-forming regions within the LIRG disks of ~ 6 % relative to their nuclei. We confront the observed trend to photo-dissociation region (PDR) models and find that the slope of the correlation is much shallower than the model predictions. Finally, we compare the correlation found between [CII]/FIR and {Sigma}_IR with measurements of high-redshift starbursting IR-luminous galaxies.
NGC 2146, a nearby luminous infrared galaxy (LIRG), presents evidence for outflows along the disk minor axis in all gas phases (ionized, neutral atomic and molecular). We present an analysis of the multi-phase starburst driven superwind in the centra
l 5 kpc as traced in spatially resolved spectral line observations, using far-IR Herschel PACS spectroscopy, to probe the effects on the atomic and ionized gas, and optical integral field spectroscopy to examine the ionized gas through diagnostic line ratios. We observe an increased ~250 km/s velocity dispersion in the [OI] 63 micron, [OIII] 88 micron, [NII] 122 micron and [CII] 158 micron fine-structure lines that is spatially coincident with high excitation gas above and below the disk. We model this with a slow ~200 km/s shock and trace the superwind to the edge of our field of view 2.5 kpc above the disk. We present new SOFIA 37 micron observations to explore the warm dust distribution, and detect no clear dust entrainment in the outflow. The stellar kinematics appear decoupled from the regular disk rotation seen in all gas phases, consistent with a recent merger event disrupting the system. We consider the role of the superwind in the evolution of NGC 2146 and speculate on the evolutionary future of the system. Our observations of NGC 2146 in the far-IR allow an unobscured view of the wind, crucial for tracing the superwind to the launching region at the disk center, and provide a local analog for future ALMA observations of outflows in high redshift systems.
We present the data and our analysis of MIR fine-structure emission lines detected in Spitzer/IRS high-res spectra of 202 local LIRGs observed as part of the GOALS project. We detect emission lines of [SIV], [NeII], [NeV], [NeIII], [SIII]18.7, [OIV],
[FeII], [SIII]33.5, and [SiII]. Over 75% of our galaxies are classified as starburst (SB) sources in the MIR. We compare ratios of the emission line fluxes to stellar photo- and shock-ionization models to constrain the gas properties in the SB nuclei. Comparing the [SIV]/[NeII] and [NeIII]/[NeII] ratios to the Starburst99-Mappings III models with an instantaneous burst history, the line ratios suggest that the SB in our LIRGs have ages of 1-4.5Myr, metallicities of 1-2Z_sun, and ionization parameters of 2-8e7cm/s. Based on the [SIII]/[SIII] ratios, the electron density in LIRG nuclei has a median electron density of ~300cm-3 for sources above the low density limit. We also find that strong shocks are likely present in 10 SB sources. A significant fraction of the GOALS sources have resolved neon lines and 5 show velocity differences of >200km/s in [NeIII] or [NeV] relative to [NeII]. Furthermore, 6 SB and 5 AGN LIRGs show a trend of increasing line width with ionization potential, suggesting the possibility of a compact energy source and stratified ISM in their nuclei. We confirm a strong correlation between the [NeII]+[NeIII] emission, as well as [SIII]33.5, with both the IR luminosity and the 24um warm dust emission measured from the spectra. Finally, we find no correlation between the hardness of the radiation field or the line width and the ratio of the total IR to 8um emission (IR8). This may be because the IR luminosity and the MIR fine-structure lines are sensitive to different timescales over the SB, or that IR8 is more sensitive to the geometry of the warm dust region than the radiation field producing the HII region emission.
We present the first results of a survey of the [CII]158um emission line in 241 luminous infrared galaxies (LIRGs) comprising the Great Observatories All-sky Survey (GOALS) sample, obtained with the PACS instrument on board Herschel. The [CII] lumino
sities of the LIRGs in GOALS range from ~10^7 to 2x10^9 Lsun. We find that LIRGs show a tight correlation of [CII]/FIR with far-IR flux density ratios, with a strong negative trend spanning from ~10^-2 to 10^-4, as the average temperature of dust increases. We find correlations between the [CII]/FIR ratio and the strength of the 9.7um silicate absorption feature as well as with the luminosity surface density of the mid-IR emitting region (Sigma_MIR), suggesting that warmer, more compact starbursts have substantially smaller [CII]/FIR ratios. Pure star-forming (SF) LIRGs have a mean [CII]/FIR ~ 4x10^-3, while galaxies with low 6.2um PAH equivalent widths (EWs), indicative of the presence of active galactic nuclei (AGN), span the full range in [CII]/FIR. However, we show that even when only pure SF galaxies are considered, the [CII]/FIR ratio drops by an order of magnitude, from 10^-2 to 10^-3, with Sigma_MIR and Sigma_IR, implying that the [CII] luminosity is not a good indicator of the star formation rate (SFR) for most LIRGs, for it does not scale linearly with the warm dust emission. Moreover, even in LIRGs in which we detect an AGN in the mid-IR, the majority (2/3) of galaxies show [CII]/FIR >= 10^-3 typical of high 6.2um PAH EW sources, suggesting that most AGNs do not contribute significantly to the far-IR emission. We provide an empirical relation between the [CII]/FIR and the specific SFR (SSFR) for SF LIRGs. Finally, we present predictions for the starburst size based on the observed [CII] and far-IR luminosities which should be useful for comparing with results from future surveys of high-redshift galaxies with ALMA and CCAT.
Nuclear stellar cusps are defined as central excess light component in the stellar light profiles of galaxies and are suggested to be stellar relics of intense compact starbursts in the central ~100-500 pc region of gas-rich major mergers. Here we pr
obe the build-up of nuclear cusps during the actual starburst phase for a complete sample of Luminous Infrared Galaxy systems (85 LIRGs, with 11.4<log[LIR/L_sun]<12.5) in the GOALS sample. Cusp properties are derived via 2-dimensional fitting of the nuclear stellar light imaged in the near-infrared by the Hubble Space Telescope and have been combined with mid-IR diagnostics for AGN/starburst characterization. We find that nuclear stellar cusps are resolved in 76% of LIRGs (merger and non-interacting galaxies). The cusp strength and luminosity increases with far-IR luminosity (excluding AGN) and merger stage, confirming theoretical models that starburst activity is associated with the build-up of nuclear stellar cusps. Evidence for ultra compact nuclear starbursts is found in ~13% of LIRGs, which have a strong unresolved central NIR light component but no significant contribution of an AGN. The nuclear near-IR surface density (measured within 1 kpc radius) increases by a factor of ~5 towards late merger stages. A careful comparison to local early-type galaxies with comparable masses reveals (a) that local (U)LIRGs have a significantly larger cusp fraction and (b) that the majority of the cusp LIRGs have host galaxy luminosities (H-band) similar to core ellipticals which is roughly one order in magnitude larger than for cusp ellipticals.
NGC 1097 is a nearby Seyfert 1 galaxy with a bright circumnuclear starburst ring, a strong large-scale bar and an active nucleus. We present a detailed study of the spatial variation of the far infrared (FIR) [CII]158um and [OI]63um lines and mid-inf
rared H2 emission lines as tracers of gas cooling, and of the polycyclic aromatic hydrocarbon (PAH) bands as tracers of the photoelectric heating, using Herschel-PACS, and Spitzer-IRS infrared spectral maps. We focus on the nucleus and the ring, and two star forming regions (Enuc N and Enuc S). We estimated a photoelectric gas heating efficiency ([CII]158um+[OI]63um)/PAH in the ring about 50% lower than in Enuc N and S. The average 11.3/7.7um PAH ratio is also lower in the ring, which may suggest a larger fraction of ionized PAHs, but no clear correlation with [CII]158{mu}m/PAH(5.5 - 14um) is found. PAHs in the ring are responsible for a factor of two more [CII]158um and [OI]63um emission per unit mass than PAHs in the Enuc S. SED modeling indicates that at most 25% of the FIR power in the ring and Enuc S can come from high intensity photodissociation regions (PDRs), in which case G0 ~ 10^2.3 and nH ~ 10^3.5 cm^-3 in the ring. For these values of G0 and nH PDR models cannot reproduce the observed H2 emission. Much of the the H2 emission in the starburst ring could come from warm regions in the diffuse ISM that are heated by turbulent dissipation or shocks.
The Spitzer Space Telescope has identified a population of ultra-luminous infrared galaxies (ULIRGs) at z ~ 2 that may play an important role in the evolution of massive galaxies. We measure the stellar masses of two populations of Spitzer-selected U
LIRGs, both of which have extremely red R-[24] colors (dust-obscured galaxies, or DOGs) and compare our results with sub-millimeter selected galaxies (SMGs). One set of 39 DOGs has a local maximum in their mid-IR spectral energy distribution (SED) at rest-frame 1.6um associated with stellar emission (bump DOGs), while the other set of 51 DOGs has a power-law dominated mid-IR SED with spectral features typical of obscured AGN (power-law DOGs). We use stellar population synthesis models applied self-consistently to broad-band photometry in the rest-frame ultra-violet, optical, and near-infrared of each of these populations and test a variety of stellar population synthesis codes, star-formation histories (SFHs), and initial mass functions (IMFs). Assuming a simple stellar population SFH and a Chabrier IMF, we find that the median and inner quartile stellar masses of SMGs, bump DOGs and power-law DOGs are given by log(M_*/M_sun) = 10.42_-0.36^+0.42, 10.62_-0.32^+0.36, and 10.71_-0.34^+0.40, respectively. Implementing more complicated SFHs with multiple age components increases these mass estimates by up to 0.5 dex. Our stellar mass estimates are consistent with physical mechanisms for the origin of z~2 ULIRGs that result in high star-formation rates for a given stellar mass. Such mechanisms are usually driven by a major merger of two gas-rich systems, rather than smooth accretion of gas and small satellites.
We have mapped the key mid-IR diagnostics in eight major merger systems of the Toomre Sequence (NGC4676, NGC7592, NGC6621, NGC2623, NGC6240, NGC520, NGC3921, and NGC7252) using the Spitzer Infrared Spectrograph (IRS). With these maps, we explore the
variation of the ionized-gas, PAH, and warm-gas (H_2) properties across the sequence and within the galaxies. While the global PAH interband strength and ionized gas flux ratios ([Ne III]/[Ne II]) are similar to those of normal star forming galaxies, the distribution of the spatially resolved PAH and fine structure line flux ratios is significant different from one system to the other. Rather than a constant H_2/PAH flux ratio, we find that the relation between the H_2 and PAH fluxes is characterized by a power law with a roughly constant exponent (0.61+/-0.05) over all merger components and spatial scales. While following the same power law on local scales, three galaxies have a factor of ten larger integrated (i.e. global) H_2/PAH flux ratio than the rest of the sample, even larger than what it is in most nearby AGNs. These findings suggest a common dominant excitation mechanism for H_2 emission over a large range of global H_2/PAH flux ratios in major mergers. Early merger systems show a different distribution between the cold (CO J=1-0) and warm (H_2) molecular gas component, which is likely due to the merger interaction. Strong evidence for buried star formation in the overlap region of the merging galaxies is found in two merger systems (NGC6621 and NGC7592) as seen in the PAH, [Ne II], [Ne III], and warm gas line emission, but with no apparent corresponding CO (J=1-0) emission. Our findings also demonstrate that the variations of the physical conditions within a merger are much larger than any systematic trends along the Toomre Sequence.
