ACTIVE DISK BUILDING IN A LOCAL H i-MASSIVE LIRG: THE SYNERGY BETWEEN GAS, DUST, AND STAR FORMATION

, , , , , , , and

Published 2010 November 30 © 2010. The American Astronomical Society. All rights reserved.
, , Citation M. E. Cluver et al 2010 ApJ 725 1550 DOI 10.1088/0004-637X/725/2/1550

0004-637X/725/2/1550

ABSTRACT

HIZOA J0836-43 is the most H i-massive (MH i = 7.5 × 1010M) galaxy detected in the HIPASS volume (δ = −90° to + 25°, v<12,700 km s−1) and lies optically hidden behind the Milky Way. Markedly different from other extreme H i disks in the local universe, it is a luminous infrared galaxy (LIRG) with an actively star-forming disk (>50 kpc), central to its ∼130 kpc gas disk, with a total star formation rate (SFR) of ∼20.5 M yr-1. Spitzer spectroscopy reveals an unusual combination of powerful polycyclic aromatic hydrocarbon (PAH) emission coupled to a relatively weak warm dust continuum, suggesting photodissociation-region-dominated emission. Compared to a typical LIRG with similar total infrared luminosity (LTIR = 1011L), the PAHs in HIZOA J0836-43 are more than twice as strong, whereas the warm dust continuum (λ>20 μm) is best fit by a star-forming galaxy with LTIR = 1010L. Mopra CO observations suggest an extended molecular gas component (H2 + He>3.7 × 109M) and a lower limit of ∼64% for the gas-mass fraction; this is above average compared to local disk systems, but similar to that of z ∼ 1.5 BzK galaxies (∼57%). However, the star formation efficiency (SFE = LIR/L'CO) for HIZOA J0836-43 of 140 L (K km s−1 pc2)−1 is similar to that of local spirals and other disk galaxies at high redshift, in strong contrast to the increased SFE seen in merging and strongly interacting systems. HIZOA J0836-43 is actively forming stars and building a massive stellar disk. Its evolutionary phase of star formation (Mstellar, SFR, and gas fraction) compared to more distant systems suggests that it would be considered typical at redshift z ∼ 1. This galaxy provides a rare opportunity in the nearby universe for studying (at z ∼ 0.036) how disks were building and galaxies evolving at z ∼ 1, when similarly large gas fractions were likely more common.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

The mechanism by which galaxies acquire their gas and form stars is an essential component of the theoretical models of galaxy formation and evolution (Dekel et al. 2009). However, large-scale cosmological simulations currently lack the dynamic range to resolve the internal dynamics of galaxies essential for modeling gas inflows associated with elevated star formation (Robaina et al. 2009). Since the physics responsible for converting gas into stars is complex, a major aspect of understanding this evolution is the relation between gas content and star formation.

The cold H i gas in galaxies represents the reservoir of material that can fuel star formation. Therefore, H i-massive disk galaxies provide an ideal laboratory to test models of disk formation, disk instability, and star formation. However, giant gas disks (MH i > 1010M) are relatively rare in the local (z < 0.1) universe and the H i mass function drops off steeply, e.g., at a mass of MH i = 6 × 1010M, the volume density is only 1.6 × 10−5  Mpc−3 (Zwaan et al. 2005). Ideally, we would like to study galaxies that are gas-rich and at an early stage of stellar building, as determined by their gas fraction, stellar mass, and star formation rate (SFR), in order to probe a relatively "pristine" stage of evolution and limit the complexity of subsequent feedback systems.

The rare examples of local massive, gas-rich disks (MH i > 5 × 1010M) include galaxies like UGC 4288 and Malin 1 which have prominent bulges, but appear underevolved, harboring extended low surface brightness disks (McGaugh 1994). Their low SFRs (van den Hoek et al. 2000) suggest that they are in a quiescent phase, for example, Malin 1 has an SFR of only 0.38 M yr−1 (Rahman et al. 2007), providing few clues as to its past and future development.

Other examples of local H i-rich disks show a range of properties. The barred Sc galaxy NGC 6744 has an H i disk diameter of 54 kpc containing MH i = 2.2 × 1010M (Koribalski et al. 2004), but an SFR of only 0.11 M yr−1 (Bohrer-Adornes & Lüdke 2004). Circinus is a Seyfert 2 galaxy with 8.1 × 109M of H i contained in a ∼120 kpc disk (Jones et al. 1999) and experiencing a modest starburst with an SFR of 2.52 M yr−1 (Elmouttie et al. 1998). M83 is similarly H i-rich with 7.7 × 109M in a barred, ∼130 kpc disk (Huchtmeier & Bohnenstengel 1981). This late-type galaxy is experiencing significantly more star formation at a rate of 6 M yr−1 (Walker et al. 1993).

In strong contrast to the extreme H i disks, these galaxies have formed large stellar disks and, in the case of M83, a large-scale bar, and appear to be in an advanced stage of evolution. The most massive source in the HIPASS BGC (H i Parkes All-Sky Survey Brightest Galaxy Catalog), ESO 390-G004, contains MH i = 3.7 × 1010M and shows clear signs of interaction, but has a very low SFR, <0.14 M yr−1 (E. Jütte & B. S. Koribalski 2011, in preparation). It seems apparent that converting H i fuel into stars can be challenging and the process by which galaxies accrete and retain their gas, and subsequently form stars, remains poorly understood.

Large surveys probing H i content, stellar mass, and SFR in the z < 0.05 universe (e.g., HIPASS, Koribalski et al. 2004; GASS, Galex Arecibo SDSS Survey, Catinella et al. 2010) provide insight as to how cold gas responds to physical conditions in galaxies and probes mechanisms responsible for regulating gas accretion and suppressing further galaxy growth from gas converting into stars.

However, studies of the early stages of stellar disk building, expected to be ubiquitous at higher redshift, are severely hampered due to observational sensitivity limits. There is growing evidence for the existence of a significant population of systems at z > 1 with star-forming disks, where gas consumption, and not major mergers, is driving stellar mass growth (Daddi et al. 2008; Genzel et al. 2010). Ideally, we should look for counterparts in the local universe where observations are easier. The recent discovery of a nearby H i-massive galaxy undergoing a powerful starburst as it builds its stellar disk (Cluver et al. 2008, hereafter Paper I) provides a rare opportunity to study this fundamental process.

HIZOA J0836-43 is the most H i-massive galaxy (MH i = 7.5 × 1010M) in HIPASS, discovered as part of a deep survey of the zone of avoidance (Kraan-Korteweg et al. 2005). Lying at a redshift of z = 0.036, the H i disk is ∼130 kpc in diameter, comparable to Malin 1, and has a dynamical mass of 1.4 × 1012M (Donley et al. 2006). However, instead of finding a relatively quiescent or even low surface brightness galaxy, typical of giant H i disks in the local universe, Spitzer Space Telescope observations revealed a luminous starburst (LTIR = 1.2 × 1011L), with spatially extended star formation. The SFR derived from the far-infrared (FIR) luminosities is ∼20.5 M yr−1 (Paper I), but the 20 cm radio continuum derived rate is ∼35 M yr−1 (Donley et al. 2006). Intriguingly, the galaxy lacks a prominent warm dust component continuum seen in other locally observed starburst systems (Paper I).

Figure 1 is a composite view of HIZOA J0836-43 showing infrared, polycyclic aromatic hydrocarbon (PAH), and optical emission (Digitized Sky Survey, DSS, R band), as well as H i, 20 cm continuum, and [Ne ii] 12.81 μm and 70 μm contours (detailed photometric properties are presented in Paper I). Due to heavy foreground extinction (AV = 7.3 mag), the galaxy is almost invisible in the optical, but the H i contours demonstrate the notable extent of the disk (∼130 kpc). Near-infrared (NIR) imaging of the galaxy reveals an early-type spiral (S0/Sa) with a prominent bulge, first noted by Donley et al. (2006). However, the combined Infrared Array Camera (IRAC) bands, with a full width at half-maximum (FWHM) of ∼2'' (see Fazio et al. 2004 for instrument details), show that in the mid-infrared (MIR) the galaxy has a star-forming disk approximately 50 kpc in diameter. The extended star formation is best shown by the [Ne ii] 12.81 μm contours overlaid on a PAH three-color image (FWHM ∼ 2'') consisting of spectral maps of the 6.2, 7.7, and 11.3 μm emission bands. The PAH emission, excited in regions of star formation, extends into the disk where there is little emission from the old stellar population that dominates the NIR image; this suggests that the galaxy is building a disk from the inside out. The PAH map and [Ne ii] 12.81 μm contours reveal a clear asymmetry on the eastern side of the disk, a feature also seen in the 20 cm continuum contours and suggested by the H i contours. The Multiband Imaging Photometer on Spitzer (MIPS; Rieke et al. 2004) 24 μm image (FWHM ∼6'') resolves the galaxy, which appears as a point source at 70 μm (shown by the contours), and indicates vigorous star formation traced by dust reprocessing starlight. As discussed in Paper I, the combined photometric spectral energy distribution (SED) is fitted well by an Sc-type galaxy template with strong PAH and dust emission (GRASIL code; Silva et al. 1998).

Figure 1.

Figure 1. Composite image showing HIZOA J0836-43 at optical (DSS R band) and infrared wavelengths; in addition, H i (FWHM ∼ 20''), 20 cm continuum (FWHM ∼ 8''), [Ne ii] 12.81 μm (FWHM ∼ 4''), and 70 μm (FWHM ∼ 18'') contours are shown. All images have a field of view of ∼3' (or ∼130 kpc). H i and 20 cm radio continuum contours are from Donley et al. (2006). JHKs imaging (FWHM ∼ 1farcs2) is from Cluver et al. (2008).

Standard image High-resolution image

In this paper, we explore in detail the spectroscopic properties of HIZOA J0836-43 using observations from the Spitzer Space Telescope and the Mopra millimeter-wave antenna. We build on the results of Paper I and provide new measurements of cold molecular gas (through CO observations on Mopra) to illustrate the seemingly unique (locally) nature of this system and discuss the importance of such a galaxy in the context of the star formation mechanism observed in disk galaxies at z ∼ 1. HIZOA J0836-43 is located behind the southern Milky Way (08h36m51fs6, −43°37'41farcs0) and we adopt a distance of DL = 148 Mpc, as found by Donley et al. (2006), derived from its recessional velocity, vhel = 10,689 km s−1.

2. OBSERVATIONS AND DATA REDUCTION

The primary data set is from the Spitzer Space Telescope GO-3 program 30914 and consists of imaging and spectroscopy centered on the galaxy. See Paper I for details of imaging observations.

2.1. Spitzer Spectroscopy

Spectroscopy of HIZOA J0836-43 was obtained on 2007 May 2 and 3 using the Infrared Spectrograph (IRS) instrument (Houck et al. 2004) on board Spitzer. Observations were done in low-resolution, R(=λ/Δλ) ∼64–128, "mapping mode" (short-low, SL) covering 5–14 μm and consisting of 3 × 60 s integrations. This set of observations consists of three separate mappings: center, east, and west of the galaxy nucleus covering 0farcm4 × 0farcm7 (${\sim} \frac{1}{3}$ of the disk) with ∼10% overlap. Outlier regions for background subtraction were observed with the same mapping scheme using a 26'' step parallel to the slit.

High-resolution (R ∼ 650) "staring mode" observations were obtained of the nucleus, covering 10–20 μm (short-high, SH) and 19–38 μm (long-high, LH) with integrations of 4 × 30 s and 4 × 14 s, respectively. A region ∼1' south of the galaxy without confusion from foreground Galactic emission was observed for background subtraction.

The Spitzer Science Center (SSC) IRS pipeline version S16.1.0 was used for spectral reductions such as wavelength and flux calibration, ramp fitting, dark subtraction, and droop and linearity corrections. Basic Calibrated Data (BCDs) of the sky background (or "outrigger") observations for the SL center and west maps were used for background removal by means of the SSC tool, CUBISM6 (Smith et al. 2007). Contaminating emission in the east outrigger observation resulted in a combination of east, west, and center outriggers being used for the east map background removal. This was done by smoothing the east background observations and subtracting them from the original frames to create an image of the east array's behavior (with the background removed). The center and west outriggers were smoothed and combined to produce a background map. This was added to the east array map to make a suitable background for CUBISM. Outlier rejection was done using the CUBISM algorithm, as well as by visual inspection of the spectral cube. CUBISM was used to extract a spectrum of a 37'' aperture, matched in SL1 (first order) and SL2 (second order), to capture the infrared disk and nuclear region of the galaxy. SL2 was scaled to match SL1 in flux (a scale factor of 1.48).

SH and LH data were reduced by subtracting median-combined background images (from the "off" position observations) and carefully cleaning individual BCDs of bad pixels and cosmic ray contamination. CUBISM was used to combine the BCDs into a spectral cube and further outliers were found by visual inspection. An SH extraction of the nuclear region was performed using CUBISM (∼9farcs25 × 4farcs5) and a matched SL extraction performed. Due to differing slit sizes, the extracted LH aperture was slightly larger and the spectrum scaled down by 55% to match the SH continuum flux at 19.3 μm.

Spectral features were measured using ISO Spectral Analysis Package (ISAP).7

2.2. Mopra CO(1–0) Observations

Molecular line observations of HIZOA J0836-43 were carried out with the 22 m Mopra8 telescope in the CO J(1–0) line (rest frequency 115.2 GHz) during 2009 May 25–27 and 2010 August 12–18. The Mopra antenna is located near Coonabarabran (NSW, Australia) at an altitude of 866 m and is part of the Australia Telescope National Facility (ATNF). Our pointing position is 08h36m51fs54, −43°37'41farcs5, which corresponds to the 20 cm radio continuum core of HIZOA J0836-43, as measured by Donley et al. (2006). Using a 3  mm receiver, the Mopra Spectrometer (MOPS) in wide-band mode covers a frequency range between 108.2 and 115.6 GHz. At the redshift of the galaxy, the CO line is centered at 111.3 GHz; each sub-band of MOPS is 2.2 GHz wide with 8096 channels (× 2 polarizations) giving a channel width of 256.25 kHz (or 0.69 km s−1). At this frequency, the antenna efficiency is 0.42 and the beam size is 32''. With this beam size, approximately half of the star-forming disk is covered; however, the molecular gas is expected to lie mostly in the nuclear and interior regions of the disk, largely covered by the Mopra beam. Hence, we expect our observations to sample the bulk of the molecular gas in this system, but should be considered a lower limit. The galaxy was observed in an "off1-on-on-off2" position-switching mode, using two different "off" positions in order to minimize the deleterious effects of 30 MHz standing wave ripples that form between the dish and receiver. "Off1" is approximately 1' north of the galaxy (in Galactic coordinates) and "off2" approximately 1' south of the galaxy. A consistent calibration was maintained by measuring the system temperature (Tsys) every 10–15 minutes using a known-temperature blackbody "paddle" and adjusting Tsys accordingly for the target observations. Typically, Tsys during our observations ranged between 250 and 400 K.

Combining data from the two observing runs, the galaxy was observed for a total of 16.6 on-source hours. The basic observations were reduced using the ATNF-developed Python-based ATNF Spectral line Analysis package (ASAP) and scripts that were customized to our observation mode. This produced calibrated spectra: ν (GHz) versus antenna temperature (K). To remove the slowly undulating baselines, a Chebyshev polynomial was fit to the line-free channels and the baseline reassessed by measuring the residual root mean square (rms) and assigning a variance to the spectrum. The typical baseline noise for an individual spectrum was ∼50 mK. Spectra with particularly poor baseline noise, >120 mK, due to the standing wave ripple were rejected. Once the antenna temperature was converted to a main-beam temperature, the individual observations and sub-bands were combined into one deep spectrum using an inverse-variance weighted-average of the spectra and performing outlier rejection. The final spectrum frequency was converted into a velocity assuming a redshift of z = 0.0356. The achieved 1σ rms sensitivity of the combined spectrum is 3.3 mK at 111 GHz. Applying Gaussian smoothing at 20 km s−1 resolution, comparable to the resolution of the Parkes H i observations of the galaxy, the rms sensitivity increases to 1.3 mK. The flux density is related to the antenna temperature according to $S_{\nu } (\rm Jy) = 2.65 \times (T_{a} / \eta) \times \theta [\mbox{$^\prime $}]^2 / \lambda [\rm cm]^2$, where η = 0.42, Ta is the antenna temperature, θ is the solid angle of the source, and λ is the wavelength.

Additional observations centered on the southeastern disk warp feature (08h36m52fs6, −43°37'47farcs6), consisting of 2.2 on-source hours, achieved a smoothed rms of 0.70 mK (smoothing at 40 km s−1 resolution).

3. RESULTS

3.1. Spectroscopy of HIZOA J0836-43

Spectroscopy of HIZOA J08360-43 allows us to probe the conditions of the interstellar medium (ISM) and analyze the star formation mechanism. The larger SL mapping coverage allows us to study the emission in the nuclear region and the inner part of the star-forming disk, but with only short wavelength coverage (5–14 μm) we can chiefly only probe the PAH emission. The measurements for the 37'' SL extraction are included and discussed in Appendix A.

We focus here on the combined spectrum (low and high resolution) of the nuclear region, presented in Figure 2. At wavelengths <15 μm, the galaxy's emission is dominated by PAH features. The excitation of these molecules is linked to large-scale star formation (Genzel et al. 1998); this is consistent with the luminous infrared galaxy (LIRG) photometric properties and SFR determined in Paper I. Measurements of the PAH emission bands are listed in Table 1 and indicated in Figure 2; they are discussed further in Section 3.2.

Figure 2.

Figure 2. Combined low- (SL; blue) and high-resolution (SH, LH; black) spectrum of the nuclear region of HIZOA J0836-43. Vertical lines indicate the locations of known emission lines (dashed purple lines), many of which are clearly detected, and prominent PAH lines (dashed red lines).

Standard image High-resolution image

Table 1. PAH Emission and Ionic Line Measurements

Line λrest FWHM EQW νFνa νLνa
  (μm) (10−3 μm) (μm) ($10^{-17}\,\rm {W}\, \rm {m}^{-2}$) (107L)
PAH 6.23   0.722 97.58 ± 2.14 66.9 ± 1.5
PAH 7.67   0.819 219.21 ± 4.76 150.0 ± 3.3
PAH 8.59   0.227 44.00 ± 1.52 30.2 ± 1.0
PAH 11.27   0.851 59.33 ± 1.03 40.7 ± 0.7
PAH 12.60   0.419 20.43 ± 1.37 14.0 ± 0.9
PAH 16.44   0.052 4.17 ± 0.47 2.9 ± 0.3
PAH 17.04   1.573 21.79 ± 2.25 14.9 ± 1.5
PAH 17.37   0.057 1.20 ± 0.41 0.8 ± 0.3
[Ar ii] 6.99b 88.40 0.039 7.55 ± 0.62 5.12 ± 0.11
[Ne ii] 12.81 32.80 0.208 13.02 ± 0.12 8.92 ± 0.08
[Ne iii] 15.56 41.09 0.051 1.40 ± 0.08 0.96 ± 0.06
[S iii] 18.71 52.40 0.165 3.18 ± 0.11 2.12 ± 0.07
[S iii] 33.47 98.30 0.156 4.52 ± 0.28 3.09 ± 0.20
[Si ii] 34.81 92.34 0.247 8.84 ± 0.48 6.05 ± 0.33

Notes. aError computed using the flux relative to the baseline rms. b[Ar ii] is blended with H2 S(2).

Download table as:  ASCIITypeset image

In the high-resolution spectra of HIZOA J0836-43, we detect the fine-structure lines of [Ar ii] 6.99 μm, [Ne ii] 12.81 μm, [Ne iii] 15.56 μm, [S iii] 18.71 μm, [S iii] 33.38 μm, and [Si iii] 34.82 μm. With the exception of the [Si ii] 34.82 μm line, these are all nebular lines from ionized hydrogen gas regions. The [Si ii] line, however, can be produced in ionized regions as well as warm atomic gas such as found in photodissociation regions (PDRs). The [O iv] 25.89 μm and [Fe ii] 25.99 μm lines are blended, but the absence of high-excitation emission in the spectrum makes [O iv] 25.89 μm unlikely. Similarly, we do not detect the high-excitation [Ne v] 24.32 μm line, which in combination with the lack of [O iv] 25.89 μm emission implies a weak, buried, or absent active galactic nucleus (AGN; Sturm et al. 2002). This agrees with the findings of Donley et al. (2006) based on their 20 cm continuum observations of HIZOA J0836-43. The measured equivalent widths (EQWs) and fluxes are shown in Table 1 and diagnostics associated with the detected lines are presented in Section 3.3.

Conspicuously absent in the spectrum shown in Figure 2 is a steeply rising warm dust continuum typically associated with an LIRG, where dust is heated due to copious star formation. Instead, it appears that the galaxy is dominated by cold dust, with only a weak contribution from MIR continuum emission. Indeed the 6 μm and 15 μm continuum measurements are more indicative of a quiescent galaxy (Paper I). Very small grains (VSGs) are believed to radiate in the 15–60 μm wavelength region, i.e., between emission from PAHs and radiation from classical interstellar grains (Puget & Léger 1989). Ongoing studies explore the possibility that the VSGs are PAH "clusters" and that isolated PAH molecules are actually produced by the destruction of these small carbonaceous grains (Boulanger et al. 1998; Rapacioli et al. 2006). The spectrum of HIZOA J0836-43 shows a distinctive lack of emission from these VSGs. This is seen quantitatively in continuum measurements; the flux densities at 15 μm and 30 μm are 0.022 ± 0.002 Jy and 0.068 ± 0.004 Jy, respectively. In their sample of starburst systems, Brandl et al. (2006) determined a relation between the total infrared luminosity (LTIR) of a galaxy and the power in its continuum. The values obtained for HIZOA J0836-43 (above) correspond to LestIR(continuum) ≈ 1.8 × 1010L. This is ∼6 times smaller than the value determined from the broadband photometry in Paper I (LTIR ∼ 1.1 × 1011L). This illustrates the relative weakness of the continuum in the nuclear region in comparison to the total infrared emission from the galaxy.

3.2. PAH Emission Line Diagnostics

Measurements of the PAH features of HIZOA J0836-43, listed in Table 1 and indicated in Figure 2, have been made using a spline continuum as described by Peeters et al. (2002) and Spoon et al. (2007). For completeness (and further comparison), we include measurements made using the PAH decomposition method (PAHFIT; as described by Smith et al. 2007) in Appendix B.

The PAH emission in HIZOA J0836-43 appears uncommonly strong, relative to the continuum, with a 7.7 μm EQW of 0.82 μm compared to the average of 0.53 μm for the starburst sample of Brandl et al. (2006). The galaxy has a 7.7 μm/11.3 μm PAH flux ratio of 3.69 ± 0.11, which is a factor of ∼2 greater than that found for typical starbursts. For the 8.6 μm and 11.3 μm PAH feature, Brandl et al. (2006) find a range of EQWs from 0.003 to 0.192 μm and 0.180to0.931 μm, respectively. HIZOA J0836-43 has corresponding EQWs of 0.227 μm and 0.851 μm by comparison. However, it is the 17 μm PAH feature which clearly shows the difference between HIZOA J0836-43 and other starbursts. The largest value is an EQW of 0.792 μm for the galaxy NGC 7252, while in strong contrast, HIZOA J0836-43 has an EQW of 1.573 μm, almost twice as large as that found for NGC 7252. We attribute these large EQWs to a combination of power in the PAH emission lines and the absence of a sharply rising dust continuum typical of starburst systems.

PAH emission is closely related to the total infrared (TIR) luminosity of a starbursting system since they both measure star formation activity. The 6.2 μm PAH flux should scale with the TIR luminosity of the system within a factor of ∼2 (Brandl et al. 2006). Using the relation of Brandl et al. (2006) derived from their starburst sample, and applying it HIZOA J0836-43, the estimated infrared luminosity is LIRest(PAH) ≈ 4.7 × 1011L. This metric overestimates by a factor of ∼4 the total luminosity determined from photometry. This further illustrates the strong (nuclear) PAH emission in relation to the total infrared emission of the galaxy, most notably the weak emission from a warm dust continuum. Strong PAH emission coupled with similarly strong continuum emission would result in a large LTIR. HIZOA J0836-43 appears atypical, as it has powerful PAH emission, but lacks similarly strong continuum emission, thus giving rise to the observed large discrepancies in the star formation metrics as compared to other starburst systems.

Brandl et al. (2006) find a relation between L(6.2 μm) PAH and νLν(5.5 μm) for low-luminosity starbursts, log [L(6.2 μm)PAH] = −0.49 + 0.96 log [νLν(5.5 μm)]. Substituting the νLν(5.5 μm) = 6.35 × 109L yields a value of L(6.2 μm) = 8.3 × 108L, close to the measured value of L(6.2 μm) = 6.7 × 108L. HIZOA J0836-43 therefore lies on the correlation, but with a luminosity greater than that seen in local starbursts, but less than that seen in ultraluminous infrared galaxies (ULIRGs) before they "turn off" the relation; see Desai et al. (2007).

The total power in the PAH emission compared to the LTIR of the galaxy can be illustrated by summing over the luminosities of the PAH emission features (see Table 1), yielding a total PAH luminosity of LPAH = 3.2 × 109L. Given that the galaxy has a TIR luminosity of LTIR = 1.2 × 1011L, the nuclear PAH emission constitutes ∼2.7% of the total infrared emission. The 10 ULIRGs studied by Armus et al. (2007) have percentages ranging from 0.4%to2.1% by comparison. In fact, we show in Section 5 (see also Figure 4) that only hyper-ULIRGs have PAHs of similar strength (relatively) compared to HIZOA J0836-43.

3.3. Ionic Emission Line Diagnostics

In the following analysis, we use the relative strengths of the emission lines (presented in Table 1) to probe the heating mechanisms, metallicity, electron density, and star formation activity in the central region of HIZOA J0836-43.

The ratios of the fine-structure line fluxes of different ionic species of the same element (e.g., $\rm {Ne^{2+}}/\rm {Ne^{1+}}$) provide a measure of the excitation and thus an indication of the hardness of the radiation field. The ratio of [Ne iii] 15.56 μm/[Ne ii] 12.81 μm (ionization potential of 41 eV/21.6 eV) is ∼ 0.11, comparable to typical values for starbursting systems (Brandl et al. 2006; Verma et al. 2003). Although the [Ne iii]/[Ne ii] ratio is not uncommon, it is below the median value (0.28) found for star-forming regions in the SINGS (Spitzer Infrared Nearby Galaxy Survey) sample (Dale et al. 2009). This implies greater weight in the [Ne ii] 12.81 μm flux compared to [Ne iii] 15.56 μm. The [Ne ii] 12.81 μm emission acts as an indication of the surface density of H ii regions as it has a much lower ionization potential compared to [Ne iii] 15.56 μm. The low [Ne iii]/[Ne ii] ratio indicates a relatively soft radiation field resulting in H ii emission and potentially large amounts of emission from PDRs.

The ratio of [S iii] 18.71 μm to [S iii] 33.48 μm is an electron density diagnostic. The [S iii] 18.71 μm/[S iii] 33.48 μm ratio of ∼ 0.70 corresponds to an electron density of $n_{e} \sim 300\, \rm {cm}^{-3}$ (Smith et al. 2009), typical of starbursting and normal disk systems (Verma et al. 2003; Dale et al. 2009).

Dense PDRs and X-ray dominated regions, powered by AGN, show strong [Si ii] 34.82 μm emission due to the significant cooling effect of this line (Hollenbach & Tielens 1999), while the [S iii] 33.48 μm line acts as a strong tracer of H ii regions (Dale et al. 2006). Compared to the AGN and star-forming regions in Dale et al. (2009), the [Si ii] 34.82 μm/[S iii] 33.48 μm ratio of 1.96 ± 0.16 is high compared to star-forming regions in SINGS, but toward the low end for the AGN distribution. Given the lack of an AGN tracer in its spectrum, it appears that HIZOA J0836-43 has large amounts of dense PDRs producing powerful [Si ii] emission which dominates the [S iii] coming from H ii regions, causing the inflated [Si ii] 34.82 μm/[S iii] 33.48 μm ratio. Dust grain destruction, such as caused by an AGN or shocks, will increase the amount of [Si ii] 34.82 μm and [Fe ii] 25.99 μm in the ISM (Dale et al. 2009). These species are therefore enhanced in active galaxies. In HIZOA J0836-43, we do not detect significant [Fe ii] 25.99 μm  which indicates that [Si ii] 34.82 μm enhancement is coming from elsewhere. We also know that the radiation field is not hard enough to encourage such a mechanism. Since [Si ii] 34.82 μm traces strong PDR emission, we conclude that the elevated [Si ii] 34.82 μm/[S iii] 33.48 μm ratio is the result of strong PDR emission compared to emission from H ii regions. We summarize in Table 2 the main properties from Sections 3.2 and 3.3.

Table 2. Summary of Spectroscopic and Derived Properties of HIZOA J0836-43

Property Value
6.2 μm PAH EQW 0.72 μm
7.7 μm PAH EQW 0.82 μm
F(7.7 μm PAH)/F(11.3 μm PAH) 3.69 ± 0.10
LPAH 3.2 × 109L
[Ne iii] 15.56 μm/[Ne ii] 12.81 μm 0.11 ± 0.01
[S iii] 18.71 μm/[S iii] 33.48 μm 0.70 ± 0.05
[Si ii] 34.82 μm/[S iii] 33.48 μm 1.96 ± 0.16
nea ${\sim} 300 \, \rm {cm^{-3}}$
LIRest(PAH)b ≈4.7 × 1011L
LIRest(continuum)c ≈1.8 × 1010L
LTIRd $\phantom{\approx}1.2 \times 10^{11} \,L_{\odot }$

Notes. aEstimated from [S iii]18.71 μm/[S iii]33.48 μm ratio. bEstimated from 6.2 μm PAH flux. cEstimated from 15 μm and 30 μm flux densities. dCluver et al. (2008).

Download table as:  ASCIITypeset image

3.4. Star Formation Rates

In this section, we use several metrics to estimate the star formation in HIZOA J0836-43. Given the galaxy's unusual properties, we expect a range of values depending on the metric used. We have seen that the PAH features appear relatively strong compared to the warm dust continuum and we expect high SFRs from PAH metrics and low SFRs from warm dust continuum metrics. Since the SFR relations have not been calibrated using galaxies such as HIZOA J0836-43, most should be interpreted qualitatively. Our relatively conservative number of 20.5 M yr−1, determined from the FIR emission, provides a credible measure of the star formation through dust reprocessed starlight; we retain this as our principal SFR for comparisons.

The strongest PAH features are useful measures of SFR, since the mechanism of their production and excitation is closely linked to star formation. Farrah et al. (2007) used the scaling between [Ne ii] + [Ne iii] and the combined 6.2 μm and 11.3 μm luminosities for their sample of 57 ULIRGs to obtain the relation ${\rm SFR}(M_{\odot }\, {\rm yr^{-1}}) = 1.18 \times 10^{-41} L_{P} (\rm {erg\, s^{-1}})$. Substituting HIZOA J0836-43's value of $L_{P} = 4.12 \times 10^{42}\, \rm {erg\, s^{-1}}$ predicts an SFR of ∼49 M yr−1, approximately twice the SFR calculated using the FIR photometry. This suggests a relative excess of PAH emission in comparison to the total IR emission of the galaxy, already seen in Section 3.2. It is, however, closer to the 20 cm radio continuum derived quantity (∼35 M yr−1) from Donley et al. (2006), assuming that star formation dominates this emission.

Houck et al. (2007) use the starburst sample of Brandl et al. (2006) to derive an SFR relation based on the scaling between νLν(7.7 μm) and LIR, yielding log[SFR] (M yr−1) = log[νLν(7.7 μm) erg s−1] −42.57. Applying the 7.7 μm luminosity from Table 1 gives an SFR (in the nuclear region) of 24.7 M yr−1.

The sum of the [Ne ii] and [Ne iii] luminosities is tightly correlated with TIR luminosity, which in turn is correlated to the SFR (Ho & Keto 2007). This relation was modified by Farrah et al. (2007), obtaining

where LN is the combined luminosity of the two neon lines, fion is the fraction of photons contributing to the ionizing the gas, and $f_{{\rm Ne}^{+}}$ and $f_{{\rm Ne}^{++}}$ are the fractional abundances of [Ne ii] and [Ne iii], respectively. Using fion = 0.6, $f_{{\rm Ne}^{+}} = 0.75$, $f_{{\rm Ne}^{++}} = 0.1$ (Ho & Keto 2007), and $L_{N} = 3.95 \times 10^{41}\, \rm {erg \, s^{-1}}$ yields an SFR of ∼20.8 M yr-1.

This line-emission-derived SFR agrees well with the SFR obtained using the FIR photometry, 20.5 M yr−1, indicating that the star formation is consistent with the strength of the ultraviolet radiation field.

The correspondence between SFR and TIR luminosity for infrared galaxies has been investigated by Rieke et al. (2009) with SFR equations derived using the MIPS L24 measurements. They find SFR = 7.8 × 10−10L24(7.76 × 10−11L24)0.048, where L24 is in L and for galaxies with LTIR > 1010L. HIZOA J0836-43 would have an SFR of 8.5 M yr−1 using this metric. Since the 24 μm dust continuum appears weak in the spectrum of HIZOA J0836-43, we would expect this to underestimate the total star formation in the galaxy.

We have SL coverage (5–14 μm) for the nuclear region and the inner MIR disk, so we can use the larger (37'') extraction area to probe the PAH emission in the disk of HIZOA J0836-43 (Appendix A) and use these values to predict the SFR. Strong PAH emission from the star-forming regions in the disk results in high luminosities (see Table 5). Using the 7.7 μm PAH value in the relation of Houck et al. (2007; see above) yields an SFR of 47.4 M yr−1. The 6.2 μm and 11.3 μm PAH luminosities in the relation of Farrah et al. (2007), see above, predict an SFR of 134.5 M yr−1. The extended star formation in the disk of HIZOA J0836-43 and the resulting PAH emission suggest a remarkably high SFR. We summarize the SFRs mentioned here in Table 3.

Table 3. Derived Star Formation Rates for HIZOA J0836-43

SFR Input Reference
49 M yr−1 6.2 μm, 11.3 μm PAHa Farrah et al. (2007)
24.7 M yr−1 7.7 μm PAHa Houck et al. (2007)
21 M yr-1 [Ne iii]/[Ne ii] Farrah et al. (2007)
8.5 M yr-1 24 μm Rieke et al. (2009)
47.4 M yr-1 7.7 μm PAHb Houck et al. (2007)
134.5 M yr-1 6.2 μm, 11.3 μm PAHb Farrah et al. (2007)
35 M yr-1c 20 cm continuum Bell (2003)
20.5 M yr-1d FIR photometry Kennicutt (1998)

Notes. aNuclear extraction. b37'' extraction. cEstimated from 20 cm continuum; Donley et al. (2006). dEstimated from FIR photometry; Cluver et al. (2008).

Download table as:  ASCIITypeset image

4. MOLECULAR GAS

4.1. Warm Molecular Hydrogen

Emission lines from pure rotational excitation of molecular hydrogen (S(0), S(1), S(2), and S(3) transitions) are also observed in the spectrum of the galaxy (Figure 2) and likely arise in PDRs neighboring H ii regions, but other possible heating mechanisms are shock heating, X-ray heating, and cosmic ray heating (Rigopoulou et al. 2002). These emission lines act as a powerful probe of the ISM as they constrain the energy injection that excites the H2. The measured properties of these lines for HIZOA J0836-43 are summarized in Table 4 and their excitation can be used to estimate the temperature and mass of (warm) molecular hydrogen. We obtain a mass of 1.3 × 107 M and a temperature of 333 K.

Table 4. H2 Emission Line Measurements (Nucleus)

Line λrest FWHM EQW Flux
  (μm) (10−3μm) (μm) ($10^{-17} \ \rm {W} \ \rm {m}^{-2}$)
H2 S(3) 9.66 32.66 0.037 1.31
H2 S(2)a 12.28 24.63 0.009 0.920
H2 S(1) 17.04 41.0 0.186 1.57
H2 S(0) 28.22 98.0 0.042 0.61

Note. aH2 S(2) is blended with [Ar ii].

Download table as:  ASCIITypeset image

Although this appears to be a relatively "typical" amount of warm H2 (albeit on the low side) for a starburst galaxy (Rigopoulou et al. 2002), it is only 0.02% of the total H i gas in the galaxy. There exists an apparent relationship between the luminosities of the 7.7 μm PAH and H2 S(1) emission lines (Rigopoulou et al. 2002), since both types of emission are thought to originate in molecular clouds with star formation. The measured values for HIZOA J0836-43 are $L_{7.7\, \mu {\rm m}} = 5.8 \times 10^{35}\, \rm {W}$ and $L_{{H_{2} {\it S}(1)}} = 3.3 \times 10^{32}\, \rm {W}$. The starbursts in the sample of Rigopoulou et al. (2002) appear to follow the relation: log L7.7 μm = 11.01 + 0.69 log $L_{{\rm H}_2 S(1)}$. Substituting L7.7 μm implies $L_{\rm {H_{2} {\it S}(1)}} \approx 7.5 \times 10^{35}\, \rm {W}$, three orders of magnitude larger than what is observed. This suggests a relative paucity of warm H2 given the LIRG starburst in HIZOA J0836-43. This combined with the exceptionally strong PAH and distinctively weak MIR continuum suggests atypical ISM conditions or heating mechanisms compared to the local starburst galaxy population.

4.2. Molecular CO

The interplay between neutral and molecular hydrogen in HIZOA J0836-43 is of primary importance to understanding the starburst mechanism and ISM conditions we observe. We therefore need to gauge the amount of cold H2 as traced by CO gas. This key observable allows us to probe the efficiency of star formation, as well as compare to other star-forming systems.

The relatively large distance of HIZOA J0836-42, 148 Mpc, combined with the exceptionally broad velocity width of the CO (1–0) line presents a serious challenge in detecting the molecular gas with the single 22 m dish Mopra telescope—the line is mostly lost in the noise at ∼5 mK sensitivity with 1 km s−1 resolution. However, if we smooth to a scale of ∼20 km s−1, which is close to the optimal bin size for the peaks in the H i rotation curve (and the resolution of the Parkes H i observations), then we clearly detect the CO emission (see Figure 3). The molecular gas has a rotation curve that is nearly identical to that of the atomic gas, at least for the velocity peaks where we have the best signal-to-noise ratio (S/N), spanning a velocity range of ± 315 km s−1.

Figure 3.

Figure 3. 12CO J(1–0) Mopra detection corresponding to 16.6 on-source hours after applying Gaussian smoothing of 20 km s−1 (solid line); the achieved 1σ sensitivity of the combined spectrum is 0.8 mK. For comparison, the H i profile from Parkes (Donley et al. 2006), scaled to match the CO peak at 10,400 km s−1, is shown as the blue dashed line.

Standard image High-resolution image

At low rotational velocities (±100 km s−1), the CO gas is not detected, lying just below the 1σ–2σ noise limit. Integrating across the entire 630 km s−1 line width, the velocity-integrated intensity is 1.59 K km s−1 (16.5 Jy km s−1), achieved with an S/N of 3.3. Using these values, the total molecular gas can be estimated. Assuming that the line is optically thin and a standard Galactic conversion factor of N(H2)/ ICO = 2 × 1020 cm−2, Lisenfeld et al. (2002) showed that the total molecular gas mass is M(H2)[M] = 75 ICO D2 Ω, where D is in Mpc and Ω is the area in arcsec2. Using this formula, we find that the H2 column density is 3.18 × 1020 cm−2, and the total molecular gas mass (H2 + He, where the He fraction is 1.38) is 3.7 × 109 M. By comparison the relations of Solomon et al. (1992) and Evans et al. (1999), using the Milky Way conversion factor (α = 4.6), both yield 3.9 × 109 M.

4.3. Gas Diagnostics

Comparing the total molecular gas mass to the total atomic hydrogen mass, the H2 to H i gas-mass fraction is ∼5%. We do not detect the low velocity gas seen in the H i distribution (see Figure 3) and given the achieved sensitivity of our observations we find it unlikely that it is present. If this galaxy behaves like most gas-rich spiral galaxies observed in the local universe, we expect that most of the gas will be located within our beam. However, to determine a complete census of the molecular gas in the galaxy, we need to fully map the entire star-forming disk (∼1' versus 30'' Mopra beam), and in particular the gas content in the "warp" region. Additional observations centered on the southeastern disk warp feature (08h36m52fs6, −43°37'47farcs6), consisting of only 2.2 on-source hours, resulted in a weak detection of molecular gas with an S/N of ∼2 (for the −300 km s−1 velocity peak). This hints at the molecular gas reservoir being extended in the disk, but additional observations are required to verify the distribution, and obtain a complete inventory of the molecular gas content. We, therefore, treat the value obtained above as corresponding to a lower limit for the total molecular gas in the disk and adopt an H2 column density of 3.18 × 1020 cm−2, corresponding to an H2 + He mass of 3.7 × 109 M (using α = 4.6), for all further comparisons in this paper.

As shown by Gao & Solomon (2004), the FIR luminosity appears correlated with the amount of molecular gas in LIRGs. This formula applied to HIZOA J0836-43 predicts 1.3 × 1010 M of H2, approximately a factor of three higher than what is implied from CO observations, supporting the possibility that our detection is missing a significant fraction of molecular gas.

The cold gas fraction, $f_{\rm {cold}} = M{\rm (H_{2})}/[M{\rm (H_{2})}+M_{*}]$, can be used to estimate the progress of stellar building. In HIZOA J0836-43, this fraction corresponds to ∼8% of cold gas, which appears low compared to the local LIRG sample of Wang et al. (2006) given the galaxy's stellar mass. This could be due to our CO measurement not capturing the extended component of the gas, however, HIZOA J0836-43 does not appear to be a typical local LIRG and we compare to these systems with caution.

5. INTERPRETATION OF OBSERVATIONS

HIZOA J0836-43 appears to have strong PAH emission, but weak dust continuum emission. As shown in Section 3, the PAH luminosity overpredicts the TIR luminosity (LTIR = 1.2 × 1011L) of the galaxy by a factor of ∼4. Conversely, the continuum fluxes underestimate LTIR by a factor of ∼6. To illustrate this discrepancy, we calculate the (synthetic) integrated 24 μm emission using the spectrum of the nuclear region of HIZOA J0836-43 (Figure 2) and compare that to the luminosity seen in the PAH bands. The spectral flux of 0.126 Jy corresponds to a luminosity of L24 = 2.4 × 109L. Using the total PAH power from Table 1 (LPAH = 3.2 × 109L), we find a ratio of LPAH/L24 in the nuclear region of ∼1.3. If we use the PAH luminosities measured using PAHFIT (Smith et al. 2007) given in Appendix B, LPAH = 9.7 × 109L and LPAH/L24 ∼ 4.0.

The unusual combination of strong PAH emission relative to a weak warm dust continuum is shown graphically in Figure 4. This plots the SED templates for luminous and ultraluminous, purely star-forming infrared galaxies, from Rieke et al. (2009), scaled to match the spectrum of HIZOA J0836-43 (using the short wavelength continuum). The warm dust continuum (λ>20 μm) of the galaxy is best fit by a star-forming galaxy with LTIR = 1010L, however, the PAH strength is more similar to a hyper-ULIRG with LTIR = 1013L. Compared to LIRGs with similar infrared luminosity (LTIR = 1011L), the 6.2 μm and 7.7 μm PAHs in HIZOA J0836-43 are ∼2 times stronger and the 11.3 μm PAH is ∼2.5 times more luminous.

Figure 4.

Figure 4. Spitzer spectrum of HIZOA J0836-43 (thick black line) with scaled SED templates from Rieke et al. (2009) for star-forming galaxies with LTIR = 1010L (blue solid), LTIR = 1011L (red dotted), LTIR = 1012L (green solid), and LTIR = 1013L (magenta solid). The relative PAH strength in HIZOA J0836-43 significantly exceeds that seen in LIRGs of similar infrared luminosity.

Standard image High-resolution image

PAH strengths and continuum fluxes can be used to determine the dominant mechanisms producing MIR emission. Armus et al. (2007) modified the diagnostic of Laurent et al. (2000) to determine the relative contributions of MIR emission sources in ULIRGs, AGNs, and starburst galaxies. Using the 15–5.5 μm continuum flux ratio (to separate starbursts and AGNs) versus the 6.2 μm EQW to the 5.5 μm continuum (to determine the quiescent star formation contribution), one can discern between AGN, H ii, and PDR-dominated emission. Using this diagnostic (shown in Figure 5), HIZOA J0836-43's MIR emission appears strongly dominated by emission arising from PDRs to the extent that emission from AGN and H ii regions is virtually insignificant. The relative strength of the PAH bands compared to the continuum emission in HIZOA J0836-43 is amongst the largest observed in any star-forming galaxy, implying unusually strong PDR-dominated emission. It appears similar to the relative emission determined for the reflection nebulae, NGC 7023 and NGC 2023 (Peeters et al. 2004), but with stronger 6.2 μm PAH emission. Both are considered prototypical of PDR emission, arising from B stars producing a small region of ionized gas and illuminating a shell of molecular gas (Wyrowski et al. 2000; Fuente et al. 2000).

Figure 5.

Figure 5. MIR diagnostic diagram from Armus et al. (2007); the three vertices represent prototypical AGN, H ii, and PDR emission, respectively, and the dashed lines show the fractional AGN contribution to the MIR emission. HIZOA J0836-43 is shown (star), as well as the ULIRGs Markarian 231 (triangle) and Arp 220 (square), and the starburst galaxy NGC 3256 (circle) for comparison (data points from Armus et al. 2007).

Standard image High-resolution image

Since HIZOA J0836-43 appears to have dust emission characteristics similar to that of NGC 7023, a reflection nebula, we explore this similarity further and plot in Figure 6 the low-resolution (SL/LL) spectrum of NGC 7023 scaled to match HIZOA J0836-43 at 10 μm. As we would expect, the composition of the two regions is quite different with HIZOA J0836-43 showing strong emission from ionic species. However, the PAH emission of NGC 7023 is remarkably similar to that of HIZOA J0836-43, further confirming our hypothesis of PDR-dominated emission. The most noticeable difference between these two systems is that for λ>20 μm, NGC 7023 has a steeply rising continuum indicating far more warm dust emission compared to HIZOA J0836-43. This arises from the H ii region producing the reflection nebula. It is apparent that the PAH emission from HIZOA J0836-43 closely resembles that of a reflection nebula such as NGC 7023, yet the dust grain composition and environment in HIZOA J0836-43 are clearly different, seen most notably in the weak dust continuum of the galaxy.

Figure 6.

Figure 6. Spectrum of HIZOA J0836-43 in comparison to the Galactic reflection nebula NGC 7023 (green) with the low-resolution Spitzer spectrum of NGC 7023 scaled to match HIZOA J0836-43 at 10 μm.

Standard image High-resolution image

The broadband photometry of HIZOA J0836-43, with an excess of 8 μm and 70 μm emission relative to 24 μm and 160 μm (Paper I), is consistent with the strong PAH emission bands and weakly rising continuum seen in the galaxy's spectrum. The spectroscopy and photometry independently indicate this distinguishing trait of HIZOA J0836-43; the origin of this dissimilarity compared to other local extragalactic objects poses a challenge to explain.

The extended star formation in HIZOA J0836-43 was discussed in Paper I by comparing the stellar light distribution with the PAH emission surface brightness (see Figure 1); the galaxy appears to be building its disk from the inside out. This gives rise to the increasing PAH luminosities when we compare the nuclear region spectrum with the larger 37'' extraction (Appendix A). For example, the total PAH luminosity increases from LPAH ∼ 3.2 × 109L in the nuclear region to ∼6.5 × 109L in the larger extraction (measured using a spline continuum). The LPAH/LTIR ratio increases from ∼3% to ∼6% as we include more of the disk (or from ∼9% to ∼20% if we use the PAHFIT values; see Appendix B). However, we note that the 37'' extraction only includes the nuclear and inner disk. We measure the extent of the PAH emission (using the 11.3 μm PAH map; see Figure 1) and determine an angular size of ∼65'', asymmetrically distributed to the east of the nuclear region. This corresponds to a star-forming disk >50 kpc in diameter. The galaxy has a substantial stellar mass in place (Mstellar = 4.4 × 1010M) with a prominent NIR bulge. Given the fuel available (MH i = 7.5 × 1010M) and the active star formation (∼21 M yr−1), we are observing a galaxy in an uncommon evolutionary state, i.e., as it vigorously builds a massive stellar disk.

6. DISCUSSION—A DISK FROM THE PAST?

In the local universe (z < 0.1), LIRGs are relatively rare comprising <5% of the total infrared energy density of local galaxies (Soifer & Neugebauer 1991). These systems are dominated by galaxy mergers and strong interactions, as evidenced by their disturbed morphologies (see, e.g., Farrah et al. 2001; Ishida 2004).

Compared to the systems in the GOALS (Great All-Sky LIRG Survey) study (Armus et al. 2009), HIZOA J0836-43 has a low L24 μm/L70 μm ratio and a weaker MIR continuum than the "coldest" LIRG in the sample, NGC 5734 (Paper I). This is an S0 galaxy (LTIR = 1.20 × 1011L) that it is interacting with the galaxy, NGC 5743, triggering the observed star formation.

LIRGs are more common at higher redshifts with ∼70% of the infrared background at z = 1 being produced by LIRGs (Le Floc'h et al. 2005), but ∼ 50% of intermediate redshift (0.1 < z < 1) LIRGs have disks that show little evidence of a recent merger (Bell et al. 2005; Melbourne et al. 2008). These LIRGs, which are predominantly gas-rich spiral systems, could be triggered by weak interactions with neighbors, minor mergers, and internal bar instabilities, achieving heightened star formation as a result of higher gas fractions compared to local disks (Melbourne et al. 2005).

We speculate that HIZOA J0836-43's current level of active star formation comes after a relatively long period of inactivity in which the galaxy would have resembled an S0/Sa type (much as Malin 1 appears today). HIZOA J0836-43 appears to have neighboring dwarf galaxies (Cluver 2009; M. E. Cluver et al. (2011), in preparation, hereafter Paper III), but no major merger interactions. As noted earlier, the infrared emission exhibits an asymmetry, or warp, along the eastern side in the PAH spectral map and this irregularity appears discernible in the 20 cm radio continuum and H i channel maps. This feature is possibly the result of recent accretion (such as a minor merger) or tidal interaction that has exacerbated a disk instability, which may have triggered the starburst as gas flows into the dense center (Paper I).

The study of Robaina et al. (2009) has shown that major interactions/mergers were an insignificant factor in stellar mass growth for z < 1, but are, however, crucial in triggering the most intense starbursts. This emphasizes the relevance of systems such as HIZOA J0836-43, which are undergoing increased star formation without the presence of a major merger.

Gas-rich galaxies are expected to be more common in the distant universe and evidence suggests that gas consumption, and not major mergers, was driving stellar mass growth (Daddi et al. 2008; Melbourne et al. 2008). Compared to local "normal" galaxies, HIZOA J0836-43 appears to be a "scaled-up" disk galaxy in terms of its gas mass and star formation properties. This is illustrated in Figure 7 plotting SFR versus H i mass for galaxies in SINGS, including here Malin 1, M82, and HIZOA J0836-43. A linear least-squares fit through the SINGS points shows that HIZOA J0836-43 lies on the high end of this relation, as opposed to Malin 1 which is quiescent by comparison, and M82, the prototypical starburst. As discussed in Paper I, the combined photometric SED is best fit by an Sc-type galaxy template dominated by cold dust emission longward of 60 μm, although with stronger PAH and dust emission, and does not match a starburst template.

Figure 7.

Figure 7. SFR vs. H i mass for HIZOA J0836-43 and SINGS galaxies. The SFRs for the SINGS galaxies are calculated from their Hα emission (Kennicutt et al. 2003). The SFR of Malin 1 is from Rahman et al. (2007) and its H i mass from Pickering et al. (1997). The M82 data are from Yun et al. (1993) and Sanders et al. (2003) for the H i and SFR, respectively.

Standard image High-resolution image

The Spitzer IRS spectrum of HIZOA J0836-43 shows weak silicate absorption indicating that the total (internal) dust extinction is low. The galaxy shows increased PAH emission compared to local spirals, yet preserves a relatively weak warm dust continuum emission. This suggests that the starburst geometry is extended, reducing the total dust column compared to compact or embedded star-forming regions. This agrees with the extended star formation observed in the disk (see Section 5). Alternatively, if regions of star formation are gas-rich, the total dust extinction could be less as the same level of star formation is sustained for a smaller dust mass. These scenarios were similarly suggested by Farrah et al. (2008) to explain distant (z ∼ 1.7) ULIRG spectra that appear more similar to local starbursts than local ULIRGs.

The recent work by Daddi et al. (2010a) has shown that normal, NIR-selected (BzK) galaxies at z ∼ 1.5 have very high gas fractions and can be viewed as scaled-up local disk galaxies hosting low efficiency star formation in extended, low-excitation gas disks. Using their LIR versus L'CO relation, including local spirals and BzK galaxies, the L'CO lower limit for HIZOA J0836-43 (8.4 × 108 K km s−1 pc2) lies on the relation exactly between these populations. Gas-mass fractions (Mgas/(Mgas + Mstellar) with Mgas including H i, H2, and He are crucial for comparisons with local spirals and more distant systems, such as the BzKs. Determining the gas-mass fraction using our current CO data yields ∼64%, comparable to the BzK galaxies (∼57%) and far above the average of 20% for local spirals (Leroy et al. 2008), despite being a lower limit.

We calculate a star formation efficiency (SFE = LIR/L'CO) for HIZOA J0836-43 of 140 L (K km s−1 pc2)−1, comparable to the average SFE for BzK galaxies of ≈100 L (K km s−1 pc2)−1 found by Aravena et al. (2010). This is below the average determined for high-redshift submillimeter galaxies (SMGs; 560 ± 210 L (K km s−1 pc2)−1; Tacconi et al. 2006) and local ULIRGs (∼ 225 L (K km s−1 pc2)−1; Yao et al. 2003). The SFE is, however, similar to values found for local spirals (Leroy et al. 2008) and for other disk galaxies at high redshift (Tacconi et al. 2010).

Similarly, Dannerbauer et al. (2009) find evidence of low-excitation molecular gas reservoirs in NIR-selected massive galaxies at z ∼ 1.5. The galaxies have SFEs and CO excitation properties similar to local spirals and the Milky Way; using a Milky Way molecular conversion factor yields molecular gas reservoirs of ∼1011M and gas-mass fractions of ⩾0.6. They speculate that the low-excitation properties observed may be typical in distant galaxies.

Tacconi et al. (2010) have shown that normal, massive star forming z ∼ 1.2 and z ∼ 2.3 are gas-rich with high molecular gas fractions. In terms of SFR and gas mass, they are scaled-up versions of normal z ∼ 0 gas-rich disk galaxies. Compared to this sample, HIZOA J0836-43 has a molecular gas fraction at the low end of the z ∼ 1.2 galaxies, but does not appear similar to the z ∼ 2.3 population. In addition, HIZOA J0836-43 lies on the SFR–Mstellar correlation determined using the z ∼ 1.2 sample, but shows no agreement with the z ∼ 2.3 fit.

This SFR–Mstellar relation was observed by Noeske et al. (2007) who found that for 2905 galaxies in the AEGIS (All-wavelength Extended Groth Strip International Survey) field there exists a "main sequence" where a given redshift corresponds to a limited range of SFRs and Mstellar. Combining samples at different redshifts from a number of studies yields the following "main-sequence line" correlation between SFR and stellar mass as a function of redshift:

(Bouché et al. 2010). Substituting HIZOA J0836-43's SFR and stellar mass corresponds to a redshift of z ∼ 0.95.

There is an emerging picture of two evolutionary tracks in observed galaxy star formation, one merger driven and the other purely contained in star-forming disks (Genzel et al. 2010). Daddi et al. (2010b) propose two modes of star formation, a long-lasting mode for disks and a more rapid mode for starbursts. They attribute this to the fraction of molecular gas in dense clouds. We find that HIZOA J0836-43 lies on the sequence of disks (when comparing LIR versus $M_{{\rm H_2}}$ and LIR/$M{_{\rm H_2}}$ versus LIR), scaled-up from spirals and in the locus occupied by BzK galaxies.

A study of 181 LIRGs at intermediate redshift (z ∼ 0.7) by Giovannoli et al. (2010) finds the median stellar mass of the sample to be Mstellar ∼ 6.3 × 1010M and a median SFR of 23 M yr−1. HIZOA J0836-43's stellar mass (Mstellar ∼ 4.4 × 1010M) and SFR (∼21 M yr−1 from the FIR photometry) are comparable to these distant systems. This combined with the gas fractions and SFEs as previously discussed, suggest a tantalizing link to a distant epoch of star formation.

HIZOA J0836-43 appears unusual in the local universe as an H i-massive LIRG. However, many galaxies have large reservoirs of H i that are not currently being converted into stars and will indeed lie dormant unless the density threshold of neutral gas is exceeded, enabling molecular clouds to form and subsequently collapse. A recent event has likely triggered the current starburst, whereas similarly extreme H i-massive disks observed in the local universe remain quiescent. This could be due to a combination of their rarity and the "catch 22" situation that the low density environment that has preserved them (Donley et al. 2006) has also left them bereft of suitable triggers. We could speculate as to how HIZOA J0836-43 formed its huge gas disk—the possibility that it is the product of a gas-rich merger of two galaxies, or due to accretion of cold gas from the intergalactic medium—and the effect of nurture due to its environment (explored further in Paper III). But, irrespective of its origin, we are observing inside-out stellar disk building in a gas-dominated galaxy, with minimal contamination from an AGN. If this galaxy mimics the pristine, early evolution of a gas disk, so critical to our understanding of the complex systems that subsequently form, HIZOA J0836-43 could provide a singular opportunity for studying a star formation mechanism (at z ∼ 0.036), which is pervasive in the z ∼ 1 universe.

7. SUMMARY

In this paper, we have presented the results of our Spitzer infrared and molecular gas observations of HIZOA J08360-43. We summarize here the main findings.

  • 1.  
    The galaxy spectrum appears dominated by strong, extended PAH emission with LPAH ∼ 3.2 × 109L in the nuclear region increasing to ∼6.5 × 109L for a region including the inner disk (measured using a spline continuum). This corresponds to LPAH/LTIR ratio of ∼3% (nuclear) to ∼6% (nuclear+inner disk).
  • 2.  
    In strong contrast, however, the dust 24 μm continuum appears intrinsically weak compared to other starburst systems and (U)LIRGs, resembling that of an "average" disk galaxy. Using the synthetic 24 μm flux, integrated from the nuclear region spectrum, corresponds to LPAH/L24 of ∼1.3. This paucity of emission from warm dust (indicating a lack of emission from VSGs) is seen independently in the galaxy's photometry.
  • 3.  
    Compared to LIRGs with similar infrared luminosity (LTIR = 1011L), the PAHs in HIZOA J0836-43 are on average ∼2 times stronger, whereas the warm dust continuum (λ>20 μm) is best fit by a star-forming galaxy with LTIR = 1010L.
  • 4.  
    The combination of strong PAH emission and relatively weak continuum places HIZOA J0836-43 in the extreme PDR-dominated region of the MIR diagnostic diagram of Peeters et al. (2004), similar to the reflection nebula NGC 7023. We speculate that the absence of a substantial warm dust continuum may be due to the dust geometry in a gas-rich disk, as the star formation appears strongly extended across a region of >50  kpc and not severely embedded.
  • 5.  
    The [Ne iii] 15.56 μm/[Ne ii] 12.81 μm ratio indicates a typical radiation field strength compared to other starbursting systems and the electron density of $n_{e} \sim 300\, \rm {cm}^{-3}$ (from the [S iii] 18.71 μm/[S iii] 33.48 μm ratio) is also normal for disk systems. The [Si ii] 34.82 μm/[S iii] 33.48 μm ratio suggests strong PDR emission compared to emission from H ii regions.
  • 6.  
    We determine a lower limit for the molecular gas content of the galaxy (∼4 × 109M); the corresponding gas-mass fraction is ∼64%, comparable to the BzK galaxies (∼57%) and far above the average of 20% for local spirals. However, the SFR (SFE = LIR/L'CO) for HIZOA J0836-43 of 140 L (K km s−1 pc2)−1 is similar to that of local spirals, BzK galaxies, and other disk galaxies at high redshift.
  • 7.  
    HIZOA J0836-43 appears scaled-up compared to local spirals in terms of its SFR, stellar mass, and gas content. The combination of SFR and stellar mass suggests that it would be considered typical at a redshift of z ∼ 1 and its stage of active building, given its gas supply, appears analogous to these distant systems.

We thank Danny Dale and SINGS, Lee Armus and GOALS for data access. We are grateful to Steve Lord, Ray Norris, Joe Mazzarella, and Justin Howell for valuable discussions and inputs. We thank the ATNF for hosting us (in Epping and Narrabri), and in particular Balt Indermuehle for his input regarding our Mopra observations. Support for this work was provided by NASA through an award issued by JPL/Caltech. M.E.C., R.K.-K., and P.A.W. thank the NRF (South Africa) for financial support.

APPENDIX A: Emission Lines in 37'' SL Extraction

The Spitzer SL mapping coverage enables us to make a large (37'' aperture) extraction of HIZOA J0836-43 which includes the nuclear region and the inner part of the MIR disk. Given the short wavelength coverage (5–14 μm), we measure predominantly PAH emission, as shown in Figure 8. Table 5 contains the measurements made using a spline continuum.

Figure 8.

Figure 8. IRS SL spectrum of HIZOA J0836-43; 37'' aperture. The blue line represents the SL2 module and the black line the SL1 module.

Standard image High-resolution image

Table 5. Emission Line Measurements—37'' Aperture (SL)

Line λrest FWHM EQW Fluxa νLνa
  (μm) (10−3μm) (μm) ($10^{-17} \ \rm {W} \ \rm {m}^{-2}$) (107L)
PAH 6.23   0.996 295.1 ± 6.9 201.3 ± 4.7
[Ar ii]b 6.98 88.2 0.043 21.17 ± 1.57 14.4 ± 1.1
PAH 7.67   0.744 423.2 ± 9.5 288.6 ± 6.4
PAH 8.61   0.106 59.87 ± 1.27 40.8 ± 0.9
PAH 11.26   0.735 141.2 ± 2.2 96.3 ± 1.5
H2 S(2) 12.28 200.0 0.157 4.54 ± 0.91 3.1 ± 0.6
PAH 12.60   0.361 33.10 ± 1.14 22.6 ± 0.8
[Ne ii] 12.81 155.0 0.281 40.47 ± 0.71 27.6 ± 0.5

Notes. aError corresponds to line flux relative to baseline rms. b[Ar ii] is blended with H2 S(2).

Download table as:  ASCIITypeset image

Summing the individual luminosities for each feature corresponds to a total PAH luminosity of LPAH ∼ 6.5 × 109L, compared to that obtained for the 9farcs25 nuclear extraction in Section 3.2 (LPAH = 3.2 × 109L). The LPAH/LTIR ratio therefore increases (compared to the nuclear extraction) to 0.06 or 6% from ∼2.7%.

APPENDIX B: PAHFIT Measurements

An alternative method for measuring PAH emission is to robustly disentangle the emission lines, dust features (e.g., PAH complexes), silicate absorption, and dust continuum components, as described in Smith et al. (2007), using PAHFIT an IDL package developed to decompose IRS low-resolution spectra. We summarize in Table 6 the PAH EQWs and luminosities obtained using PAHFIT on the nuclear region extraction of HIZOA J0836-43; these can be used to compare to PAH studies that utilize this method (for example, Smith et al. 2007).

Table 6. PAHFIT Emission Line Equivalent Widths and Luminosities (Nuclear Region)

PAH λrest EQW Flux Luminosity
(μm) (μm) ($10^{-17} \ \rm {W} \ \rm {m}^{-2}$) (107L)
 5.27 0.186 12.7 8.7
 5.70 0.211 16.9 11.5
 6.22 2.030 182.0 124.2
 6.69 0.489 45.9 31.3
 7.42 1.420 132.0 90.1
 7.60 2.760 255.0 174.0
 7.85 2.630 236.0 161.0
 8.33 0.656 53.5 36.5
 8.61 1.610 119.0 81.2
10.68 0.179 8.2 5.59
11.23 1.020 48.0 32.8
11.33 2.230 106.0 72.3
11.99 0.975 48.0 32.8
12.62 1.700 81.7 55.7
12.69 0.097 4.66 3.2
13.48 0.527 23.7 16.2
14.19 0.124 5.31 3.6
16.45 0.224 8.03 5.5
17.04 0.761 26.1 17.8
17.38 0.151 5.02 3.4

Download table as:  ASCIITypeset image

Summing the individual luminosities for each feature corresponds to a total PAH luminosity of LPAH ∼ 9.7 × 109L, approximately three times the luminosity obtained in Section 3.2 (LPAH = 3.2 × 109L) using the spline continuum measurements. This corresponds to an LPAH/LTIR ratio of 0.09 or 9%. We compare this result to the study of Smith et al. (2007) who used the PAH decomposition method to investigate PAH emission in 59 star-forming galaxies in SINGS, determining that the majority of galaxies have LPAH/LTIR close to 0.1, with lower and upper limits of ∼0.02 and 0.11, respectively. However, given the extended star formation in HIZOA J0836-43's ratio of 0.06 is on the low side compared to the majority of their galaxies, especially since it is an LIRG with a relatively high SFR.

Table 7 compares the luminosities of the strongest PAH complexes (6.2, 7.7, and 11.3 μm), measured for the nuclear region, to the luminosities obtained using the spline continuum (from Section 3.2). The L(decomp)/L(spline) ratios are similar to what has been found for star-forming galaxies in SINGS (Smith et al. 2007) as shown in Table 7.

Table 7. Comparison between PAH Feature Luminosities—PAHFIT versus Spline Continuum

PAH Lpahfit $\frac{L_{\rm pahfit}}{L_{\rm spline}}$ $\frac{L_{\rm pahfit}}{L_{\rm spline}}$a
  (108L)    
 6.2 μm 12.42 1.86 1.70
 7.7 μmb 42.5 2.83 3.53
11.3 μmb 10.51 2.58 1.86

Notes. aSmith et al. (2007). bPAH complex.

Download table as:  ASCIITypeset image

Using PAHFIT on the large extraction (37''), shown in Figure 8, allows us to measure the PAH emission contained in the nuclear region and the inner MIR disk. The results are summarized in Table 8.

Table 8. PAHFIT Emission Line Equivalent Widths and Luminosities (37'' Extraction)

PAH λrest EQW Flux Luminosity
(μm) (μm) ($10^{-17} \ \rm {W} \ \rm {m}^{-2}$) (107L)
 5.27 4.060 58.9 40.2
 5.70 3.070 73.6 50.2
 6.22 13.50 517.0 352.6
 6.69 4.020 217.0 148.0
 7.42 5.970 470.0 320.5
 7.60 5.330 442.0 301.4
 7.85 4.910 448.0 305.5
 8.33 1.290 132.0 90.0
 8.61 2.050 219.0 149.4
10.68 0.085 11.1 7.57
11.23 0.660 93.9 64.0
11.33 1.390 201.0 137.1
11.99 0.583 89.5 61.0
12.62 0.936 147.0 100.3
12.69 0.104 16.3 11.1
13.48 0.205 31.5 21.5

Download table as:  ASCIITypeset image

The total PAH luminosity in the 37'' extraction is LPAH ∼ 2.2 × 1010L, corresponding to LPAH/LTIR ratio of 0.2 or 20%. The LPAHLν(24 μm) ratio is ∼2; this corresponds to twice as much power coming from PAH features compared to warm dust emission. Using the PAH luminosities for the two extractions listed in Tables 6 and 8, we can calculate LPAHLν(24 μm) ratios of the 6.2, 7.7, and 11.3 μm features. We list these in Table 9, as well as values from the study of Smith et al. (2007).

Table 9. Ratios of PAH Luminosity to 24 μm Luminosity

Ratio LPAHLν(24 μm) LPAHLν(24 μm)a
L6.2 μmLν(24 μm) ∼0.11b ∼0.1
L7.7 μmLν(24 μm) ∼0.46b ∼0.4
L11.3 μmLν(24 μm) ∼0.10b ∼0.1
L6.2 μmLν(24 μm) ∼0.32c ∼0.1
L7.7 μmLν(24 μm) ∼0.84c ∼0.4
L11.3 μmLν(24 μm) ∼0.18c ∼0.1

Notes. aSmith et al. (2007). bNuclear region. c37'' extraction.

Download table as:  ASCIITypeset image

Despite the luminosity ratio's in Smith et al. (2007) showing a large scatter, the large extraction of HIZOA J0836-43 exceeds the values obtained for the SINGS sample. For example, the L6.2 μmLν(24 μm) and L11.3 μmLν(24 μm) exceed the highest ratios seen in star-forming galaxies in SINGS (Smith et al. 2007). Given that we are not extracting over the entire star-forming disk, we expect this ratio to be even larger.

Footnotes

  • The ISO Spectral Analysis Package (ISAP) is a joint development by the LWS and SWS Instrument Teams and Data Centers. Contributing institutes are CESR, IAS, IPAC, MPE, RAL, and SRON.

  • The Mopra radio telescope is part of the Australia Telescope which is funded by the Commonwealth of Australia for operations as a National Facility managed by CSIRO.

Please wait… references are loading.
10.1088/0004-637X/725/2/1550