Abstract
We present spatially resolved morphological properties of [C II] 158 μm, [O III] 88 μm, dust, and rest-frame ultraviolet (UV) continuum emission for A1689-zD1, a strongly lensed, sub-L* galaxy at z = 7.13, by utilizing deep Atacama Large Millimeter/submillimeter Array (ALMA) and Hubble Space Telescope (HST) observations. While the [O III] line and UV continuum are compact, the [C II] line is extended up to a radius of r ∼ 12 kpc. Using multi-band rest-frame far-infrared continuum data ranging from 52 to 400 μm, we find an average dust temperature and emissivity index of K and
, respectively, across the galaxy. We find slight differences in the dust continuum profiles at different wavelengths, which may indicate that the dust temperature decreases with distance. We map the star formation rate (SFR) via IR and UV luminosities and determine a total SFR of 37 ± 1M⊙yr−1 with an obscured fraction of 87%. While the [O III] line is a good tracer of the SFR, the [C II] line shows deviation from the local L[C II]-SFR relations in the outskirts of the galaxy. Finally, we observe a clear difference in the line profile between [C II] and [O III], with significant residuals (∼5σ) in the [O III] line spectrum after subtracting a single Gaussian model. This suggests a possible origin of the extended [C II] structure from the cooling of hot ionized outflows. The extended [C II] and high-velocity [O III] emission may both contribute in part to the high L[O III]/L[C II] ratios recently reported in z > 6 galaxies.
Export citation and abstract BibTeX RIS
![](https://content.cld.iop.org/journals/0004-637X/934/1/64/revision1/apjac795blicense.gif)
Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
1. Introduction
The formation and evolution of galaxies across cosmic time is shaped by the "baryon cycle," i.e., the complex interplay between gas accretion, star formation, stellar and active galactic nucleus feedback, and gas outflows (e.g., Davé et al. 2012; Anglés-Alcázar et al. 2017; Tumlinson et al. 2017). These various processes serve to transfer gas between the interstellar medium (ISM) and the surrounding circumgalactic medium (CGM). Key to understanding the role of the baryon cycle in the formation of the first galaxies are studies investigating the morphology and structure of the ISM and CGM at high redshift, where robust observational constrains are rare.
At high redshift, far-infrared (FIR) fine-structure lines are typically used to trace the ISM and the galaxies' star formation activity. In particular, the [C II] 2 P3/2 → 2 P1/2 transition at 157.74 μm (hereafter [C II]), with its low ionization potential of 11.2 eV (compared to 13.6 eV of hydrogen) is a dominant coolant of the neutral ISM (Stacey et al. 1991; Wolfire et al. 2003), and dense photodissociation regions (PDRs; Hollenbach & Tielens 1999) associated with giant molecular clouds. The [C II] line, arising from singly ionized carbon, has been used to probe ISM properties (e.g., Capak et al. 2015; Knudsen et al. 2016; Pentericci et al. 2016; Carniani et al. 2017; Bakx et al. 2020; Matthee et al. 2020) and kinematics (e.g., Jones et al. 2017; Smit et al. 2018; Rizzo et al. 2020, 2021) in dozens of high-redshift galaxies. While the [C II] line traces multiple phases of the ISM, it has been shown to primarily trace cool, neutral gas through PDRs (Cormier et al. 2019). By contrast, the [O III] 3 P1 → 3 P0 transition at 88.4 μm (hereafter [O III]) primarily traces hot, diffuse, ionized gas, as the ionization potential of O+ (35.1 eV) is significantly higher than that of hydrogen. Moreover, both the [C II] and [O III] lines are known to trace the star formation rate (SFR; De Looze et al. 2014; Kapala et al. 2014; Herrera-Camus et al. 2015).
The advent of the Atacama Large Millimeter/submillimeter Array (ALMA) opens our observational window to detect both the [C II] and [O III] lines at z > 4, and the number of detections has risen rapidly in recent years (see Hodge & da Cunha 2020). Of particular relevance, recent ALMA stacking studies have discovered extended [C II] structures up to r ∼ 10–15 kpc around high-redshift galaxies (Fujimoto et al. 2019; Ginolfi et al. 2020a). The extended structure, called the "[C II] halo" (Fujimoto et al. 2019), shows an effective radius ∼5 times larger than that of the central galaxy as seen in the UV, as well as a potential association with the Lyα halo (e.g., Leclercq et al. 2017). Similar extended [C II] structures have also been identified in observations of individual star-forming galaxies at z > 4 (Fujimoto et al. 2020; Herrera-Camus et al. 2021) and massive quasar host galaxies at z > 6 (Maiolino et al. 2012; Cicone et al. 2015; see Meyer et al. 2022).
In several cases, observations of the extended [C II] structure have additionally detected the "broad-wing" feature in the [C II] spectrum (Gallerani et al. 2018; Ginolfi et al. 2020a; Herrera-Camus et al. 2021), characteristic of ongoing high-velocity outflows, which have been hypothesized to drive carbon-enriched gas out into the CGM. The existence of extended [C II] emission is not fully reproduced by the current galaxy formation models (Katz et al. 2019a, 2019b; Fujimoto et al. 2019; Pallottini et al. 2019; Arata et al. 2020; Katz et al. 2022), suggesting that there may be missing physics in simulations at the Epoch of Reionization (Pizzati et al. 2020).
Fujimoto et al. (2019) discuss several possible scenarios for what physical processes may power the extended [C II] emission, including outflows, satellite galaxies (in which [C II] emission traces their SFR), and photoionization via radiation from the central galaxy. Ginolfi et al. (2020a) and Fujimoto et al. (2020) find correlations between galaxy SFR and the detection of the broad-wing feature and [C II] extension, suggestive of the SF-driven outflow scenario. However, this result may be biased by the signal-to-noise ratio, as objects with higher SFR are likely brighter, making faint components in the spectra easier to detect (Ginolfi et al. 2020a; Fujimoto et al. 2020). Ginolfi et al. (2020b) introduce an additional scenario, tidal disruption, which can produce [C II] emission on CGM scales through turbulence and shocks through gravitational interactions. Moreover, while the existence of the extended [C II] structure implies the presence of dust in the CGM, high-redshift ALMA observations have generally found a compact dust continuum (Fujimoto et al. 2019; Ginolfi et al. 2020a). In contrast, observations of local galaxies suggest that dust entrained by galactic winds can reach ∼10 kpc in scale (Kaneda et al. 2009; Meléndez et al. 2015; McCormick et al. 2018; Yoon et al. 2021), and theory suggests that stellar feedback is efficient at driving dust-enriched, multiphase winds (Kannan et al. 2021). ALMA observations of quasar host galaxies at z ≳ 6 have found [C II] and dust emission extended out to r ∼ 10 kpc, and interpret both as tracing the star-forming ISM of these massive galaxies (Novak et al. 2020). A proper determination of the physical processes producing the extended [C II] structure in line with the outflow thus depends on our ability to measure, with sufficient depth, the distribution of neutral/ionized gas and dust continuum emission. Simultaneous ALMA observations of the [C II] and [O III] fine-structure lines, along with their underlying continua, are well-suited to accomplish this task.
In this work, we utilize deep ALMA observations of the galaxy A1689-zD1, a strongly lensed, sub-L* galaxy at z = 7.13 (Bradley et al. 2008; Watson et al. 2015; Knudsen et al. 2017). The high magnification (μ ∼ 9.3) of this galaxy makes it an ideal target for a spatially resolved investigation of the ISM in the abundant population of sub-L* galaxies in the Epoch of Reionization, and indeed, recent high-resolution follow-up observations have successfully detected bright [C II] 158 μm and [O III] 88 μm emission lines (K. Knudsen et al. 2022, in preparation). Moreover, rich FIR continuum data have also been taken in bands 3, 6, 7, 8, and 9, providing valuable constraints on the FIR luminosity and dust temperature (Bakx et al. 2021). The availability of deep, high-resolution, multiwavelength ALMA data makes A1689-zD1 the ideal target to address the aforementioned issues and examine the extended [C II] emission around high-redshift galaxies.
The paper is structured as follows: In Section 2, we describe the data reduction steps for both ALMA and HST data. We report the results for the spatial extent of the [C II], [O III], UV, and FIR emission in Sections 3.1 and 3.2. In Section 3.3, we spatially resolve the dust temperature, which we use in Section 3.4 to derive the spatially resolved L[C II]-SFR and L[O III]-SFR relations for A1689-zD1. We present the line spectra of [C II] and [O III] in Section 3.5. Finally, in Section 4 we discuss our results, including the interpretation of the extended emission as a "[C II] halo" versus an extended disk, and the physical origin of the emission. Throughout this work, we assume a ΛCDM cosmology with H0 = 67.4 km s−1 Mpc−1, ΩM = 0.315, ΩΛ = 0.685 (Planck Collaboration et al. 2020) and a Kroupa (2001) initial mass function (IMF).
2. Data Processing
2.1. ALMA Data
We utilize ALMA observations of A1689-zD1 (R.A. 13h 11m 299, decl.
) in bands 3, 6, 7, 8, and 9. The observations in bands 6 and 8 (ALMA program IDs 2015.1.01406.S and 2017.1.00775.S, respectively) will be presented in K. Knudsen et al. 2022 (in preparation). We use archival observations for bands 7 (Knudsen et al. 2017), and 9 (Bakx et al. 2021), and refer the reader to these papers for more details. While previous studies of A1689-zD1 have adopted the optical redshift of z = 7.5 determined from an apparent continuum break (interpreted as the Lyα; Watson et al. 2015), new observations of the [C II] and [O III] lines at high S/N place the spectroscopic redshift at z = 7.132. As such, the band 6 observations were tuned to the [C II] 158 μm line at νobs = 233.71 GHz. Similarly, bands 8 and 9 were tuned to the [O III] 88 μm, and [O III] 52 μm lines, respectively. Analysis of the [O III] 52 μm line will be presented in a separate, upcoming paper. Band 3 observations were tuned to cover the CO(7-6) and CO(6-5) lines.
The Common Astronomy Software Applications package (CASA; McMullin et al. 2007) was used for reduction, calibration, and imaging. The data were first processed via the standard ALMA pipeline reduction, which was sufficient for these observations. All images were produced with the tclean task in CASA, using natural weighting to maximize sensitivity. Cleaning was done using the "auto-multithresh" automasking algorithm (Kepley et al. 2020), initially down to the 3σ level but expanded to a lownoisethreshold of 1σ. In order to image both compact and extended emission, we apply the multiscale deconvolver (Cornwell 2008) with scales of zero (delta-function) one, and three times the beam size (i.e., the psf). Continuum maps were produced using the multifrequency synthesis (mfs; Conway et al. 1990) mode in CASA over a frequency range identified as lacking line emission. To produce moment 0 maps for the [C II] and [O III] lines, we first subtract the continuum using the CASA task uvcontsub, with order 0, and produce a map using the mfs mode. While this mode is typically used for continuum imaging, we employ it to produce an average line intensity map, as a pseudo-moment 0 map, in order to maximize the signal-to-noise ratio of the map prior to cleaning and ensure that we detect and deconvolve any faint and/or diffuse signals in the tclean process. We integrate channels between 233.4 and 234.0 GHz for the [C II] map and between 416.6 and 417.8 GHz for the [O III] map. Both ranges correspond to roughly ±400 km s−1 from the line center.
For the [C II] and [O III] line intensity maps, we obtain an rms of 18.8 μJy beam−1 and 63.0 μJy beam−1. For the band 3, 6, 7, 8, and 9 continuum maps, we obtain respective rms values of 10.4, 5.3, 44.9, 24.5, and 181.2 μJy beam−1. All maps are primary beam corrected, though the reported rms is derived prior to primary beam correction. All ALMA and HST images are aligned astrometrically using the reproject package in Python.
With natural weighting, the synthesized beam sizes for band 6 and 8 observations are 030 × 0
28 and 0
46 × 0
41, respectively. To perform a fair comparison based on the same spatial resolution, we apply a uv-taper of 0
3 for the band 6 data, which increases the beam size to 0
44 × 0
40. For the analysis of extended emission (Section 3.2), we apply a uv-taper of 0
5 (0
7) to observations in band 8 (band 6). Table 1 summarizes the beam size and depth of both the fiducial (uv-taper of 0
3 for band 6, untapered for other bands) and tapered (0
7 for band 6, 0
5 for other bands) continuum data in each band.
Table 1. Summary of ALMA Continuum Observations for A1689-zD1
Band | λrest |
![]() ![]() | θfid | θtap | σfid | σtap | References |
---|---|---|---|---|---|---|---|
(μm) | (μm) | (arcsec) | (arcsec) | (μJy beam−1) | (μJy beam−1) | ||
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) |
3 | 402.5 | 375.8, 429.2 | 0.55 × 0.46 | ⋯ | 10.4 | ⋯ | K. Knudsen et al. 2022 (in preparation) |
6 a | 163.0 | 156.2, 169.9 | 0.44 × 0.40 | 0.82 × 0.73 | 5.3 | 8.3 | K. Knudsen et al. 2022 (in preparation) b |
7 | 107.2 | 104.7, 109.7 | 0.62 × 0.58 | 0.78 × 0.74 | 44.9 | 47.1 | Knudsen et al. (2017) |
8 | 89.5 | 87.8, 91.1 | 0.46 × 0.42 | 0.72 × 0.69 | 24.5 | 28.7 | K. Knudsen et al. 2022 (in preparation) b |
9 | 53.0 | 52.6, 53.4 | 0.52 × 0.43 | 0.67 × 0.61 | 181.2 | 205.0 | Bakx et al. (2021) |
Notes. Columns: (1) ALMA band, (2 and 3) Central wavelength and range, (4 and 5) Beam size, for "fiducial" maps (with no uv-taper) and "tapered" maps (with a uv-taper of 05), (6 and 7) rms noise, (8) Reference for ALMA data.
![farcs](https://cdn.images.iop.org/Entities/farcs.gif)
![farcs](https://cdn.images.iop.org/Entities/farcs.gif)
Download table as: ASCIITypeset image
2.2. HST Data
In addition to deep ALMA observations of FIR line and continuum emission, we use archival Hubble Space Telescope (HST) data to probe the rest-frame UV continuum of A1689-zD1 (Watson et al. 2015). We obtain images taken with Hubble WFC3/IR in the J (F125W) and H (F160W) bands. At z = 7.13, these bands trace the UV continuum at 1350 Å < λ < 2050Å. These wavelengths are redder than the Lyman limit and unaffected by the Lyα forest, making this a good probe of the rest-frame UV continuum and the unobscured SFR.
Recent works (e.g., Rujopakarn et al. 2016; Dunlop et al. 2017; Fujimoto et al. 2019; Tamura et al. 2019; Herrera-Camus et al. 2021) have shown that there sometimes exist offsets of 01–0
3 between HST and ALMA astrometry, which can significantly impact results of comparisons between high-resolution observations. The HST data we use in this work have been corrected for this astrometric offset by comparison to astrometry from the Gaia survey (Gaia Collaboration Collaboration et al. 2018), which is accurate to ALMA astrometry to the milliarcsecond scale. The correction shifts the HST astrometry − 0
11 in R.A. and + 0
06 in decl.
As noted in previous analyses of A1689-zD1 (Watson et al. 2015; Knudsen et al. 2017), there is an additional source ∼ 15 westward of the target that has been identified as a low-redshift (z ∼ 2) galaxy. Other sources near A1689-zD1 can be identified as low-z interlopers galaxies by their detection in shorter-wavelength ACS/F814W observations, where A1689-zD1 disappears due to the Lyman-break nature of the galaxy. If a pixel has S/N > 3 in F125W+F160W, is more than 1
1 from A1689-zD1, and is also detected in ACS/F814W, we regard it as an interloper. We replace all such interlopers with randomly drawn values from the image rms noise. Finally, in order to present a fair comparison with the same spatial resolution between HST and ALMA images of A1689-zD1, we convolve the HST map with a Gaussian kernel constructed to match the idealized (i.e., Gaussian model) ALMA beam for the band 8 data. For the fiducial (tapered) maps, this kernel has a size of 0
44 × 0
39 (0
70 × 0
67).
2.3. Source Plane Reconstruction
To analyze the intrinsic ISM morphology of A1689-zD1, we reconstruct the ALMA and HST images into the source plane, correcting for the lens magnification and deflection.
The magnification of A1689-zD1 is derived to be a factor of μ = 9.3 ± 0.5 based on the mass model originally published in Limousin et al. (2007), and refined to include spectroscopic information of multiple systems as described in Bina et al. (2016). This parametric strong lensing model is constructed using the lenstool (Jullo et al. 2007) software, which provides a set of models sampling the posterior distribution of each parameter of the mass distribution, allowing us to derive the statistical error on the magnification. The uncertainty of 0.025 dex of the magnification at the location of the source was also confirmed independently as described in Watson et al. (2015), and the magnification is quite uniform across the source, with variations around 5% (see Figure A1 in Appendix). We show source plane reconstructed maps for each emission in Figure A2; we employ these source plane maps when measuring distance scales (i.e., Sections 3.2.2 and 3.2.3). This is because the asymmetrical deflection of the source confuses a simple relation between angular and physical scale in the image plane—the full reconstruction is necessary to assess the size/shape of the source.
When quoting the luminosity or flux from the source, we employ the notation that, for example, LIR represents the observed IR luminosity and LIR μ−1 represents the intrinsic IR luminosity, corrected for magnification.
3. Results
3.1. Two-dimensional Maps of [C II], [O III], UV, and FIR Emission
The top row of Figure 1 presents intensity maps of the [O III] and [C II] lines as well as the dust continuum at ∼90 μm and ∼163 μm and the UV continuum as probed by stacking the HST J and H bands. Blue crosses in the maps indicate the positions of the two UV continuum peaks, one in the southwest and one in the northeast. The peaks of the [O III] and [C II] maps are likely to be co-spatial with the southwest UV-bright region. This is also where the peak of the dust continuum emission lies, suggesting that this southwest component is the more actively star-forming region of the two. The consistency of the peak positions of the [O III] and [C II] with the FIR and UV continuum indicates that the majority of the line emission likely arises from the young star-forming regions visible in UV.
Figure 1. Top: image plane maps of emission from [O III] 88 μm, [C II] 158 μm, band 8 dust continuum, band 6 dust continuum, and the rest-frame UV continuum in A1689-zD1. In each panel, white dashed contours show the −3σ and −2σ levels, while colored contours show the 2σ and 3σ levels and then increase following the Fibonacci sequence. A slight uv-taper of 03 is applied to both band 6 maps to match the resolution of the band 8 maps, and the ellipses on the bottom left indicate the beam size. The rest-frame UV continuum is measured by stacking the HST H and J bands, masking out nearby sources, and convolving with a Gaussian kernel constructed to match the ALMA beam in band 8. For clarity, we show the unconvolved (high-resolution) HST map in the inset panel for the image plane. We pinpoint the two peaks of the UV continuum emission on each panel in the image plane with two blue crosses. All panels show the same 5
6 × 5
6 area on the sky. Bottom: visibility profiles of each ALMA observation. The amplitude is extracted from the real part of the complex visibility in bins of uv-distance, which corresponds inversely to the observed angular scale θ. We fit each visibility profile to a single Gaussian, and find that this is a good fit for the [O III] line and dust continuum. For [C II], however, a double Gaussian model provides a better fit, indicating the presence of significant extended emission. The bottom right panel shows histograms of uv-dist for the two ALMA bands. The detection of the extended [C II] structure is not due to an increased sensitivity to large-scale structure (i.e., small uv-dist) in band 6.
Download figure:
Standard image High-resolution imageThese images indicate that the [C II] emission is likely extended relative to the [O III], FIR, and UV. [C II] emission is detected at >2σ out to r ∼ 2– from the galaxy center. Furthermore, the dust continuum at 90 μm (band 8) appears more compact than at 163 μm (band 6), which we will revisit in Section 3.3.
3.2. Spatial Extent of the [C II] Emission
3.2.1. uv-plane Visibility Profiles
We also analyze the spatial structures of the different emission in the uv-visibility plane using the Python package uvplot (Tazzari 2017). In the bottom panel of Figure A2, we show the real part of the visibility amplitude in bins of the uv-distance, which is inversely proportional to the observed angular scale θ.
We apply a single and double Gaussian fits to the [O III], [C II], and dust continuum profiles (shown in blue and red, respectively). Based on the reduced chi-squared () statistic, the single Gaussian is a reasonably good model for the [O III] line and the dust continuum in both bands, with consistent FWHMs of ∼1
0; on the other hand, for the [C II] line, we find a significant offset from the best-fit single Gaussian at the shortest uv-distance regime far beyond the error scales. The double Gaussian shows a better fit for the [C II] line, lowering the
. For the compact and extended [C II] components, we derive a best-fit FWHM of 1
05 ± 0
03 and 4
0 ± 0
6, respectively. The FWHM of the compact [C II] component is almost the same as those of the [O III] line and the dust continuum, while the extended [C II] component is ∼4 times larger than the compact component. These results suggest that there exists an independent extended component surrounding the central galaxy, which is reminiscent of the Lyα halo universally identified around normal star-forming galaxies at z ∼ 3–6 (e.g., Momose et al. 2014; Leclercq et al. 2017). The size ratio of the extended and compact components is comparable to the previous ALMA deep [C II] stacking results of ∼5–6 estimated from galaxies and quasars at z ∼ 6 (Fujimoto et al. 2019; Novak et al. 2020).
We note that the existence of the extended [C II] emission, despite there being no extended component in the [O III] emission, is likely not just a result of the data depth or maximum recoverable scale. The S/N at the peak pixel in the [C II] intensity map is 42.9, almost comparable to [O III] (37.7), with a difference of only ∼10%. Moreover, as the histogram of uv-dist in Figure 1 indicates, the band 8 observations cover more short baselines (larger angular scales) than band 6, due to the different antenna configurations (C-4 and C-3 for bands 6 and 8, respectively).
Despite the larger error bars, we note that the profiles of the dust continuum in both bands show some excess emission at large angular scale (low uv-dist), similarly to the [C II] line. Although the lower S/N of this excess trend prevents us from giving a definitive conclusion with regard to the existence of the extended dust component, we examine the potential extended dust structure also in the image-based analysis (Section 3.2.2).
3.2.2. Image-based Radial Profiles
We derive radial surface brightness profiles, based on images for both emission lines and the UV and FIR continuum. We conduct this analysis in the source plane, as described in Section 2.3. For the radial profiles, we use the midpoint between the two UV-bright peaks as the center.
Figure 2 shows radial profiles in the source plane for [C II], [O III], and UV continuum (top) and for the dust continuum at 90 μm and 163 μm (bottom), normalized at their peaks. Points show the mean surface brightness in circular annuli whose widths increase at increasing distances, from 0.4 kpc to 1 kpc. Error bars show the standard error of the mean. For ease of readability, we show these radial profiles out to the first point at which the mean surface brightness drops below 0. At this point and beyond, we consider the profile dominated by noise fluctuation and instead show only 1σ upper limits.
Figure 2. Radial surface brightness profiles for A1689-zD1 in the source plane. In the top panel, emissions from the [C II], [O III], and UV continua are shown in red, magenta, and blue, as in Figure 1. In the bottom panel, we show emission from the dust continuum at rest-frame 163 μm (green) and 90 μm (orange) over the profiles for [C II], [O III], and the UV. The gray dotted line show the profile of the ALMA PSF. Error bars show 1σ errors on the mean surface brightness in the annulus, computed as the image , where N is the number of full beams in the annulus. When the profile first drops below 0, we consider it dominated by noise, and thus show 1σ upper limits. We detect extended [C II] emission out to 12 kpc from the galaxy center, inconsistent with the [O III] and UV profiles out to 6–8 kpc. The dust continuum follows a surface brightness profile similar to that of [C II].
Download figure:
Standard image High-resolution imageBoth [O III] and UV continuum emission are compact (though extended relative to the PSF), dropping to below 0 at radii of ∼4.5 and ∼3 kpc, respectively. In contrast, [C II] emission is detected out to r ∼ 12 kpc. Furthermore, the [C II] profile shape is inconsistent with the [O III] profile out to ∼6 kpc, and inconsistent with the UV profile out to ∼8 kpc.
Dust continuum emission appears to be extended with a radial profile similar to that of the [C II] line. This is true at both 163 μm and 90 μm, though the latter wavelength drops to 0 at ∼8 kpc while the former extends out to ∼12 kpc. 14 Both continuum profiles are inconsistent with UV and [O III] at radii of r ∼ 2–4 kpc, more in line with the [C II] profile. As with the visibility-based profiles, the large error bars at larger scales prevent robust evidence of extended emission. Nevertheless, the consistent positive surface brightness in the dust continuum may suggest some real extended emission. The presence of such extended dust continuum emission would be in contrast to previous ALMA results for star-forming galaxies at z ∼ 4–7, and we discuss this further in Section 4.
3.2.3. Gini Coefficients
While radial profiles give a quantitative sense as to the extent of [C II] emission, they may also be sensitive to the "clumpiness" of the emission. Recent analysis has shown that A1689-zD1 is a highly clumpy object, requiring multiple spatial/spectral components to fit its [C II] and [O III] data cubes (Wong et al. 2022, Knudsen et al. in prep). As such, the precise choice of center may impact the shape of the surface brightness profile. In the imaging analysis presented thus far, we have adopted the midpoint between the two UV-bright peaks as the center; however, in order to test against this possible bias, we additionally compute the Gini coefficient for each emission. The Gini coefficient (G) is a metric borrowed from econometrics, and was originally created to summarize the degree of income inequality in a population. It was first used to assess galaxy morphology by Abraham et al. (2003), who showed that it can, to first order, be interpreted similarly to other morphological indicators such as concentration. The advantage of using G in this work is that it is nonparametric: it does not depend on the galaxy having a well-defined center, but rather only indicates whether the bulk of the total intensity is contained within a few pixels or spread evenly across many. A high value of G (∼0.75) indicates that the emission from the object is concentrated in one or more bright nuclei, while a low value of G (∼0.25) indicates that the emission is more uniform.
We compute the Gini coefficient following Lotz et al. (2004) by first sorting all n pixel values Xi into increasing order and then computing
![Equation (1)](https://content.cld.iop.org/journals/0004-637X/934/1/64/revision1/apjac795beqn1.gif)
We compute the error on G in a Monte Carlo fashion by varying the underlying map by a randomized "noise map," convolved with the ALMA beam. Since the Gini coefficient is sensitive to the number of pixels included, we compute G based on all pixels within a circular aperture of radius r = 2'' for each map. Moreover, while the Gini coefficient can be sensitive to the S/N of the map, this effect is negligible above a per-pixel S/N of ∼10 (Lisker 2008), and the sensitivity to the S/N is accounted for by the Monte Carlo error computation.
Figure 3 shows the Gini coefficient computed in the source plane for each map. Our Gini coefficient results reaffirm the conclusions from the radial profiles: [O III] and UV emission are compact and concentrated in one or two bright nuclei, while [C II] and dust continuum emission are spread more evenly. Interestingly, the Gini coefficient for the dust continuum suggests a difference between the different wavelengths: the dust continuum at shorter wavelengths is more concentrated than at higher wavelengths. This may suggest that the dust temperature may vary in space, with hotter dust (which peaks at shorter wavelengths) more concentrated than colder dust (which peaks at higher wavelengths). The spatially resolved dust temperature distribution is addressed further in Section 3.3.
Figure 3. Gini coefficients for A1689-zD1. We compute the Gini coefficient for each map following Equation (1) and using a circular aperture of radius 2''. Error bars show 1σ uncertainty on the Gini coefficient measured in a Monte Carlo fashion; specifically, for n = 1000 iterations, we sum the underlying image with a randomized noise map (corresponding to the image rms and convolved with the ALMA beam), and compute the resulting Gini coefficient. Error on the Gini coefficient is inversely correlated with the S/N of the map, but all maps have sufficient S/N to constrain the Gini coefficient within ∼5%. While the [O III] line and UV continuum are very concentrated, the [C II] line FIR continuum is more uniform.
Download figure:
Standard image High-resolution image3.3. Spatially Resolved Dust Temperature
The rich FIR continuum data for A1689-zD1 facilitates accurate determinations of the dust temperature. Importantly, band 9 (rest-frame 53 μm) data are taken, which helps significantly to constrain the FIR SED at shorter wavelengths. However, the low S/N of the band 9 observations prohibits using this lower-wavelength data to constrain the dust temperature in a spatially resolved sense. We opt instead to use the ratio between bands 6 (163 μm) and 8 (90 μm) to spatially resolve the dust temperature; however, this method depends on the assumed dust emissivity index, β, which usually ranges from 1 to 2.5 (Casey 2012). Therefore, we first determine the best-fit β value via a modified blackbody fit on the galaxy as a whole, using continuum data from bands 3, 6, 7, 8, and 9. We then derive the spatially resolved dust temperature map, assuming that β remains constant in space.
Figure 4 shows the FIR SED of A1689-zD1. Measurements of the dust continuum flux from bands 6 (163 μm), 7 (107 μm), 8 (90 μm), and 9 (53 μm) are shown, and a 3σ upper limit is shown for the continuum in band 3 (∼400 μm), which is not detected. All photometry is computed in the 28 × 1
8 elliptical aperture shown over the dust continuum images on the right. This aperture is adopted consistently across all bands and is determined by eye to correspond roughly to the [O III]-emitting region. The aperture is large enough to cover the bulk of the dust continuum in bands 6–8, but not so large as to include outskirt noise fluctuation in band 9 and significantly reduce the S/N in that band. Since band 9 holds most of the constraining power in our MBB fit, it is critical to focus on this central region where the band 9 continuum is robustly detected. Error bars in Figure 4 show the combined uncertainty, including both statistical uncertainty measured from the image rms and an additional systematic uncertainty of 10% (20% in band 9) arising from flux calibration, as specified in the ALMA Cycle 8 proposer's guide.
15
We fit a modified blackbody spectrum to the data points, correcting for the decreasing temperature contrast between the target and background due to the heating of dust by the CMB following da Cunha et al. (2013). We perform this fitting using a Monte Carlo approach, allowing Tdust and β to vary freely with flat priors. We determine a best-fit dust temperature of
K and dust emissivity index
. The inset panel shows the distribution of Tdust and β values derived from the MCMC fitting, and highlights the inherent degeneracy between the two parameters. The larger uncertainty in our Tdust and β measurements compared to Bakx et al. (2021) is due to our use of a 20% uncertainty on the band 9 measurement.
Figure 4. The far-infrared SED of A1689-zD1. We fit a modified blackbody spectrum to ALMA photometry from bands 9, 8, 7, and 6, derived from 28 × 1
8 elliptical apertures placed over the continuum maps shown in the right panels. We additionally constrain the SED fitting according to the 3σ upper limit from the nondetection of continuum in band 3 (rest-frame ∼400 μm). The gray shaded region shows the 68% range of best-fit SEDs determined from varying the flux values according to their uncertainties in a Monte Carlo fashion. The inset panel plots the 1σ contour of the probability distribution of Tdust and β values derived from the Monte Carlo method, and highlights the degeneracy between the two parameters. We find a best-fit dust temperature of
K and dust emissivity index of
, where the error represents the 68% confidence interval. The images on the right show the four continuum detections, with contours denoting the 2, 3, 5, 8, 12, 16, 20, and 25 × σ levels, and dashed yellow lines denoting the aperture.
Download figure:
Standard image High-resolution imageTo extend the analysis in a spatially resolved fashion, we then fix β = 1.7 as a fiducial value. We note that, with variations in the dust chemical composition and crystallinity, it would be possible for β to vary in space with Tdust fixed. However, we have fit the full FIR SED in two spatially resolved apertures and found <1σ difference in β between the two. Therefore, we consider it reasonable to fix β and explore the resulting variation in Tdust. We do so by deriving a relationship between the observed rest-frame 90 μm to 163 μm continuum ratio (corrected for CMB effects) and the dust temperature. Figure 5 shows the spatially resolved dust temperature map for pixels where the dust continuum is detected (at 2σ) at both 163 μm and 90 μm. We show in the right panel the fractional uncertainty on the inferred dust temperature. While we see three regions of high dust temperature, the two on the far northeast and northwest sides of the galaxy are likely a product of noise fluctuation, as the uncertainty is >60%.
Figure 5. Dust temperature map for A1689-zD1. Blue contours show the UV continuum as in Figure 1. The right panel shows the fractional uncertainty on the inferred dust temperature. We see high Tdust in the southeast UV-bright region, corresponding to the bright [C II] and [O III]-emitting region, but the two other high-Tdust regions are likely a product of noise.
Download figure:
Standard image High-resolution imageWe overplot contours showing the stacked HST J + H-band image. The southeast UV-bright region corresponds to a peak in the dust temperature of ∼50 K, while the northwest region corresponds primarily to a noisy peak. However, the high dust temperature in the galaxy center stretches out toward both noisy peaks, with low errors. Similarly, the far southwest of the map shows a higher temperature than the minimum, with low errors. These regions likely correspond to distinct clumps in the galaxy ISM, which have been identified in the galaxy spectrum (Wong et al. 2022, Knudsen et al. in prep).
In addition to estimating Tdust on a pixel-by-pixel basis, we also endeavor to estimate the spatially resolved Tdust on kpc scales. Figure 6 shows the inferred dust temperature as a function of distance, in arcsec, based on image-plane radial profiles of the continuum maps at both wavelengths. We use circular annuli of variable width out to 2'' from the galaxy center, which corresponds to ∼7 kpc when accounting for the gravitational lensing. We use image-plane maps for this analysis, to minimize any additional uncertainty brought about by the source-plane reconstruction. Error bars show 1σ uncertainty on the dust temperature. The radial profile suggests that the dust temperature decreases with increasing distance from the galaxy center, starting at ∼50 K within 025 and dropping to ∼30 K beyond 1''. At further distances, it appears to increase again, though our signal-to-noise becomes too low to robustly measure the dust temperature.
Figure 6. Radial profile of the dust temperature in A1689-zD1. We derive the dust temperature from the ratio of the continuum at 90 μm to the continuum at 163 μm, correcting for CMB effects. As in Section 3.2, we use uv-tapered maps. Error bars show 1σ uncertainty on the inferred dust temperature derived in a Monte Carlo fashion from the uncertainty on the continuum ratio. We show the farthest bin as an open diamond, as the S/N is too low to provide a stringent constraint on the dust temperature. Despite this uncertainty, the dust temperature appears to decrease to larger radii, approaching the CMB temperature.
Download figure:
Standard image High-resolution imageWe explore the feasibility of a Tdust gradient by applying a constant model in which the dust temperature is fixed at the best-fit spatially integrated value (Tdust = 41 K) everywhere. This model results in a reduced chi-squared of . In contrast, a simple linear model for the decreasing trend agrees within the error in every annulus. These results indicate that a model with constant Tdust is unlikely to be consistent with the data. While the errors in our results are large, this analysis suggests that deep ALMA observations, with the aid of gravitational lensing, may be able to observe the Tdust gradient at z > 7.
3.4. The Σ[C II]–ΣSFR and Σ[C II]–ΣSFR Relations
Spatially resolving the dust temperature facilitates robust estimates of the spatially resolved total IR luminosity (8–1000 μm) and obscured SFR. For each pixel, we compute the total IR luminosity by integrating a modified blackbody with β = 1.7 and Tdust set to the inferred dust temperature shown in Figure 5, or if the temperature was not well-determined, 16 with Tdust = 30 K. We then compute the obscured SFR following the Murphy et al. (2011) calibration as
![Equation (2)](https://content.cld.iop.org/journals/0004-637X/934/1/64/revision1/apjac795beqn2.gif)
We additionally compute the unobscured (UV-based) SFR following the Murphy et al. (2011) calibration as
![Equation (3)](https://content.cld.iop.org/journals/0004-637X/934/1/64/revision1/apjac795beqn3.gif)
where LUV is the UV luminosity measured from the stacked HST J+H-band image. This probes the rest-frame wavelength range 1350 Å < λ < 2050Å, comparable to the GALEX FUV transmission curve used in the calibration of Murphy et al. (2011). The calibrations from Murphy et al. (2011) assume a Kroupa (2001) IMF and a constant star formation history with an age of ∼100 Myr. Uncertainties on the derived SFRs include the uncertainties on the ALMA/HST maps as well as the uncertainty in the dust temperature (for SFR IR). Summing over the elliptical aperture used for photometry in Section 3.3, and assuming β = 1.7, we derive an obscured SFR of SFRIR = 32 ± 1M⊙yr−1 and an unobscured SFR of SFRUV = 4.95 ± 0.03M⊙yr−1, resulting in a total SFR of ∼ 37 ± 1M⊙yr−1 and an obscured fraction of ∼87%. The high obscured SFR is in line with recent estimates for A1689-zD1 from Bakx et al. (2021). We derive a total infrared luminosity of LIR μ−1 = (2.16 ± 0.07) × 1011 L⊙.
To compare the derived SFR to [C II] and [O III] emission, we place circular annuli of varying width across the tapered maps, and compute the SFR and line surface densities in each annulus. Figure 7 shows the spatially resolved Σ[O III]–ΣSFR and Σ[C II]–ΣSFR relations with the empirical, local universe relation from De Looze et al. (2014) and (for Σ[C II]) the results for z ∼ 2–6 dusty star-forming galaxies from Spilker et al. (2016). We see that the [O III] line roughly traces the SFR following the relation for local, low-Z, dwarf galaxies, with a slight (∼0.3 dex) excess in regions with high ΣSFR (i.e., the galaxy center). By contrast, we see a ∼0.7 dex deficit of [C II] in the galaxy center, relative to local dwarfs. This deficit is consistent with the [C II]-FIR deficit in high-redshift DSFGs observed by Spilker et al. (2016), which has also been observed in local IR luminous galaxies (Díaz-Santos et al. 2014). The deficit at high ΣSFR diminishes at lower ΣSFR, with the two largest annuli showing [C II] detections despite nondetections of UV or IR continua. The [O III] line is also not detected in these annuli.
Figure 7. The spatially resolved Σ[O III]–ΣSFR and Σ[C II]–ΣSFR relations for A1689-zD1. Red points show spatially resolved data derived from circular annuli of varying widths. Gold lines show empirical spatially resolved, local universe relations from De Looze et al. (2014), while blue points show measurements of individual z ∼ 2–6 star-forming galaxies from Spilker et al. (2016), assuming ΣSFR = 1.48 × 10−10 ΣFIR. The right panels show the [O III] and [C II] contours plotted over the SFR map, with the borders of the circular annuli.
Download figure:
Standard image High-resolution image3.5. Central [C ii] and [O iii] Line Profiles
We also investigate the [C II] and [O III] line profiles in a spatially resolved manner. In these emission line spectra, high-velocity outflows in different gas phases may be identified via the "broad wings" that extend past the extent of the best-fit Gaussian profile (see, e.g., Veilleux et al. 2020). We produce line profiles from the [C II] and [O III] cubes in a small (025 radius) aperture centered at the the peak of the ΣSFR map, i.e., over the UV+FIR-bright region. This small radius is roughly the beam size of the fiducial ALMA maps. We focus on this region for two reasons: 1) the clumpy nature of A1689-zD1 confuses identification of the broad-wing feature with distinct spectral components, and 2) in a SF-driven outflow scenario, the faint broad-wing feature would be most easily detected in regions of highest ΣSFR (Davies et al. 2019). We have verified that this aperture is small enough to probe only one of the five kinematic component components identified in Knudsen et al. (in prep). Moreover, as the same aperture is applied to both the [C II] and [O III] maps, any differences in the line profiles corresponds to different kinematics between [C II] and [O III] rather than contamination from different components, which would be seen in both emission lines.
Figure 8 shows the line profiles for [C II] and [O III] in this small aperture, with Gaussian fits applied. Before performing any Gaussian fits, we determine the FWHM of each line and examine the differences in the line widths. 17 The [C II] line width of 252 ± 21 km s−1 is significantly larger than the [O III] line width of 182 ± 31 km s−1. The smaller line width for the [O III] line is consistent with hydrodynamic simulations (e.g., Katz et al. 2019b, 2022), as the [O III] tends to trace younger star-forming regions, which tend to have a lower velocity dispersion than the system as a whole. The larger line width for [C II] may also be related to the spatial extent of the [C II] emission. If the extended emission were spherically symmetric, then the larger volume of gas observed may capture more complicated kinematics along the line of sight and thus contribute to the larger line width.
Figure 8. Line profiles for [O III] (left) and [C II] (right) in a 025 aperture centered on the high-ΣSFR region. We apply Gaussian fits to both lines with the FWHM fixed at the measured value. The bottom panel shows the residuals from these fits, annotated with the S/N of the residuals integrated across the high-velocity regions at ±100–400 km s−1. While the [C II] line is broader than the [O III] line, the [O III] line shows significant high-velocity residuals, suggesting the possibility of ongoing hot ionized outflows from the galactic core.
Download figure:
Standard image High-resolution imageTo test for the outflow signature, we fit a Gaussian to each line with the FWHM fixed at the measured value. We see that the [C II] line is fit well by a single Gaussian, with a reduced chi-squared value of . By contrast, the [O III] line is not fit well by a single Gaussian, with a reduced chi-squared value of
(p < 10−7). This is owed to significant excess flux in high-velocity regions, despite the narrow line width overall. We show in the bottom panels of Figure 8 the residuals on each Gaussian fit. For [O III], we find significant (>5σ) residuals when integrating across the high-velocity tails from ±100 to 400 km s−1. As [C II] is fit well by a single Gaussian, we find no significant residuals in these high-velocity tails.
4. Discussion
We have shown in Section 3 that the [C II] emission (and possibly also the dust continuum emission) in A1689-zD1 is extended on ∼12 kpc scales relative to [O III] and UV emission. The [C II] line emission in uv space is fit well by a double Gaussian model, with an extended/compact component ratio of ∼4. This [C II] emission deviates slightly from the local Σ[C II]–ΣSFR relation, with a "deficit" of [C II] in the galaxy center and a slight excess in the furthest apertures. We additionally detect significant large-velocity residuals in the [O III] line spectrum at the galactic core, indicating the possible presence of outflow dominated by hot, ionized gas. We now discuss in depth the physical interpretation of these observations; in particular, (1) the description of the extended [C II] emission as a part of the central galaxy versus a "[C II] halo," separate from a single galaxy disk and reminiscent of the Lyα halo, and (2) the physical origin of the extended emission.
4.1. [C ii] Halo versus Extended Disk
The results presented thus far are reminiscent of the "Lyα halo" universally identified around normal star-forming galaxies at z ∼ 3–6 (e.g., Momose et al. 2014; Leclercq et al. 2017). However, another possibility is that we are observing the outskirts of the central galaxy disk even out to r ∼ 10 kpc, and that star formation in this disk is the source of the extended [C II] emission (e.g., Novak et al. 2020). Due to the detection of extended FIR emission (Figure 2), we cannot immediately rule out the latter scenario.
While the rest-UV and [O III] sizes are similarly compact, we see potential evidence of extended dust emission in both uv- and image-based analyses (e.g., Figures 1 and 2). The presence of extended dust continuum emission would be in contrast to previous ALMA results for star-forming galaxies at z ∼ 4–7 (Fujimoto et al. 2019; Ginolfi et al. 2020a; Herrera-Camus et al. 2021), and would, at first glance, suggest significant obscured star formation in the outskirts of the galaxy.
One possible explanation for this inconsistency with previous observations is that this extended dust emission is common in the early universe, but too cold and faint to be detected, being embedded in the warm CMB (e.g., da Cunha et al. 2013; Vallini et al. 2015; Zhang et al. 2016; Pallottini et al. 2017; Lagache et al. 2018). The dust continuum profiles for A1689-zD1 are similar to the UV and [O III] line up to r ∼ 3 kpc, and only diverge at larger distances. This faint extended structure may have been missed in previous stacking results, especially if the observed dust temperature gradient is common, as it is in local galaxies (Hunt et al. 2015) and as is predicted in cosmological zoom-in simulations for early galaxies (e.g., Behrens et al. 2018; Sommovigo et al. 2020). For this study, the deep integration times in bands 6 and 8, combined with the strong lensing effect (with magnification μ ∼ 10) yield an effective sensitivity ∼10 times deeper than previous stacking results (e.g., Fujimoto et al. 2019; Ginolfi et al. 2020a). Under this interpretation, many high-z star-forming galaxies may host extended dust continuum emission corresponding to their extended [C II] emission. In fact, recent stacking results for z ≳ 6 massive quasar host galaxies find similar radial profiles for [C II] and dust (Novak et al. 2020), although the physical mechanisms producing extended dust emission might be different between massive quasar host galaxies and sub-L* galaxies.
It is worth noting that the detection of extended dust continuum emission implies a larger dust mass than what has previously been estimated (e.g., from Bakx et al. 2021). Indeed, the elliptical aperture applied in Section 3.3 fails to capture ∼13% of the band 6 continuum flux observed out to r ∼ 1
5. Based on our spatially resolved Tdust estimates, fitting this excess flux to a modified blackbody with Tdust = 30 K and β = 1.7 yields a dust mass of Mdust ∼7 × 106
M⊙. This suggests that the dust mass estimate could be increased by ∼ 40% compared with the previous estimate of Bakx et al. (2021) by including the extended dust component. If the extended dust component ubiquitously exists around early galaxies, this additional, previously unaccounted for dust mass may have implications for dust evolution in the early universe (e.g., Mancini et al. 2015; Michałowski et al. 2019; Dayal et al. 2022).
The extended [C II] and dust continuum emission could be explained by localized star formation, i.e., an extended SF disk rather than an independent halo component. Given that the dust temperature at the outskirts of the galaxy approaches the CMB temperature (Figure 6), it is possible that we underestimate ΣIR in these regions despite correcting for the decreasing contrast following da Cunha et al. (2013). However, even in this scenario, the compact rest-UV size of the galaxy implies that the obscured fraction increases up to ∼100% toward the outskirts of the galaxy. We show in Figure 9 the UV and IR SFRs (top) and obscured fraction (bottom) as a function of radius in several circular annuli (as in Figure 7). While the obscured fraction decreases with increasing radius up to r ∼ 1'', the detections of IR SFR despite nondetections of UV SFR yield 3σ lower limits on the obscured fraction of ∼93%. Such a reversal of the decreasing trend is in contrast to the metallicity gradient expected from simulations, which would continue to decrease to ∼1% at the outskirts (Pallottini et al. 2017).
Figure 9. ΣSFR and obscured fraction, assuming that the extended FIR continuum emission is solely due to localized SF. Top: surface density of star formation derived from the UV (blue) and IR (red), as a function of radius in several circular annuli. Bottom: obscured fraction of star formation (SFRIR/SFRtot) as a function of radius. In both panels, arrows indicate 3σ upper/lower limits. Nondetections of the UV SFR in the outer annuli, where IR emission is detected, constrain the obscured fraction to ≳93%. Such a high obscured fraction in the outskirts is unexpected, and may suggest that the extended FIR emission is not entirely due to SF.
Download figure:
Standard image High-resolution imageTherefore, a realistic picture might include some dust heating mechanisms other than localized SF activity in the outskirts of the galaxy, such as heating from the interstellar radiation field (ISRF) or by hot outflows. In this scenario, we would need to correct the SFR estimates by removing the extended FIR emission counted as SFR. Then, the Σ[C II]–ΣSFR relation in the furthest annuli would deviate even more from the local relation, strengthening the interpretation that extended [C II] emission is powered by heating mechanisms other than the localized SF.
At the same time, the [C II] emission in the outskirts of the galaxy can also be suppressed by the CMB, depending on the gas temperature and density (Vallini et al. 2015; Kohandel et al. 2019). Since in general the gas temperature ( ≳102 K) is much higher than TCMB, the detection of the extended [C II] emission may indicate (i) the extended C+ gas maintains a relatively high density that is sufficient to prevent it from being significantly affected by the CMB suppression, or (ii) the majority of the [C II] emission in the outskirt regions is actually suppressed and the intrinsic [C II] luminosity is even greater. The latter case would yield even further deviation from the local Σ[C II]–ΣSFR relation (right panel of Figure 7.
As A1689-zD1 is a low-mass ( ∼109
M⊙), sub-L* galaxy (Watson et al. 2015), the existence of an extended disk of r ∼ 10 kpc would be striking. The effective radius for such sub-L* galaxies is generally estimated to be <0.8 kpc (e.g., Shibuya et al. 2015).
18
Following studies of the Lyα halo (e.g., Hayes et al. 2014; Momose et al. 2014; Leclercq et al. 2017), we fit the [C II] and [O III] profiles using a combined Sérsic+exponential model to fit the central and halo components, respectively. Figure 10 shows these fits. For the [O III] line, we find that including a second component is unnecessary, as the emission is fit well by a single Sérsic profile with Reff = 1.21 ±0.01 kpc and n = 0.59 ± 0.02. This is generally consistent with the rest-frame optical and UV sizes measured by Shibuya et al. (2015), though perhaps a bit more extended. For the [C II] line, the single Sérsic model fails to capture the extended emission; instead, we need a central component with and n = 0.55 ± 0.03 and an exponential halo component with
. In addition to the visibility-based analysis (Section 3.2.1), the fitting clearly confirms the existence of two separate components, which is reminiscent of the Lyα halo identified around high-z SFGs.
Figure 10. Radial profiles of the [O III] (left) and [C II] (right) emission, with two-component Sérsic+exponential profiles applied. The [O III] emission is fit well by a single Sérsic profile with Reff = 1.21 ± 0.01 kpc and n = 0.59 ± 0.02. By contrast, the [C II] line requires a central Sérsic profile with and n = 0.55 ± 0.03 as well as an exponential profile with
. The clear secondary component suggests the existence of a [C II] halo.
Download figure:
Standard image High-resolution imageIn summary, the extended dust continuum and [C II] emission might be explained by the outskirt of an extended galaxy disk. But this interpretation suggests an increasing trend of the obscured fraction as a function of radius up to ∼100% at large radii. Instead, the extended [C II] structure is more likely explained by an independent gaseous halo, which is evident from both visibility- and image-based profile fitting for the [C II] emission, and we may call it the "[C II] halo." The [C II] halo and the extended dust emission would have formed via some physical mechanism(s) other than localized SF, which we will discuss in the next section.
4.2. The Physical Origin of the [C ii] Halo
Assuming that the extended [C II] emission represents a distinct [C II] halo, we now discuss the physical origin of this feature. Fujimoto et al. (2019) present five possible explanations for [C II] emission on circumgalactic scales: A) star formation from satellite galaxies; B) circumgalactic-scale photodissociation region (CG-PDR); C) circumgalactic-scale H II region (CG-H II); D) cold streams; and E) outflow. For merging galaxies in the center of the protocluster, Ginolfi et al. (2020b) argue for a sixth possible explanation: F) turbulence and shocks associated with tidal stripping.
In all cases, some physical process is required in order to enrich the CGM with carbon and produce the halo of C+ that we observe. Assuming that cold streams are composed primarily of pristine gas, only scenarios A), E), and F) include some mechanism to enrich the CGM with carbon. Scenarios B) (CG-PDR) or C) (CG-H II) may certainly ionize the carbon and power [C II] emission, provided the CGM was pre-enriched by some other process. However, given that A1689-zD1 is a sub-L*, pre-mature galaxy at z = 7.13 whose formation history is relatively short, we assume that the metal-enrichment process is related to ongoing activity and focus our discussion on the three scenarios of (A) star formation, (E) outflow, and (F) tidal disruption.
4.2.1. Star Formation from Satellite Galaxies
As discussed in Section 4.1, we find that the spatially resolved Σ[C II]–ΣSFR relation shows a slight excess of [C II] emission relative to the local relation at low SFR (at the outskirts of the galaxy). Such deviation from the local relation suggests that the extended [C II] emission is unlikely to be explained by individual, low-mass, low-SFR satellite galaxies. Moreover, following the mass–metallicity relation (e.g., Mannucci et al. 2010), such faint, low-mass satellite galaxies generally have low metallicity. At this low-Z, the [C II] line emissivity for a given stellar continuum decreases (Vallini et al. 2015; Olsen et al. 2017; Katz et al. 2022), which would in fact have an effect opposite to what we observe.
4.2.2. Outflows
Another possible scenario for the origin of the [C II] halo is that high-velocity outflows drive carbon out to r ∼ 10 kpc. These outflows are responsible for the metal enrichment of the CGM, and this interpretation has been supported by observations of ongoing outflows, as traced by the [C II] line, in high-SFR galaxies (Gallerani et al. 2018; Ginolfi et al. 2020a; Herrera-Camus et al. 2021). Since the majority of [C II] is thought to be radiated from PDRs in the cold, neutral hydrogen gas clouds, the broad-wing feature observed in [C II] suggests that the outflowing gas is dominated by cool, neutral gas.
However, galactic outflows in general are known to be made up of multiple phases of gas, with the archetypal example of M82 showing a hot phase observed in X-ray emission (T ≳ 106 K), a warm ionized phase (T ∼ 104 K), cold and warm molecular phases, and a cold atomic phase (Heckman & Thompson 2017). The multiphase nature of galactic outflows suggests that the [C II] line may not reliably trace outflows in all galaxies. In fact, Spilker et al. (2020) find that molecular outflows are ubiquitous in z > 4 galaxies, but not always identifiable by [C II] emission. We show in Section 3.5 that the [O III] line shows significant ∼ 5σ residuals in the high-velocity region after fitting a single Gaussian model. By contrast, we find no significant residuals in [C II]. This may suggest the presence of ongoing outflows dominated by hot, ionized gas, rather than cool, neutral gas in A1689-zD1.
The theoretical connection between these hot outflows and the observed, cool, [C II] halo depends on the cooling time of the hot ionized gas (tcool). Optically thin intergalactic gas at z ∼ 3.5 (with temperature T ∼ 106 K and an overdensity δ ∼ 1–10) is predicted to have a cooling time of tcool ≳ 4 Gyr, which is longer than the Hubble time at this redshift (Madau et al. 2001). However, this depends on the specific properties of the gas, e.g., the density, metallicity, and multiphase structure of the outflows. In fact, recent "cooling outflow" models implemented into the analytical models of Pizzati et al. (2020) have successfully reproduced the [C II] halo out to r ∼ 10 kpc. In their models, hot ionized outflows experience catastrophic cooling shortly after their initial launch. Depending on the timescale on which this hot outflowing gas cools, we may or may not expect to observe the broad wings in [C II] under this scenario. This might contribute to the diversity of the presence or the absence of the [C II] broad-wing feature in recent ALMA high-z studies (e.g., Decarli et al. 2018; Gallerani et al. 2018; Bischetti et al. 2019; Fujimoto et al. 2019; Stanley et al. 2019; Ginolfi et al. 2020a; Novak et al. 2020; Spilker et al. 2020; Izumi et al. 2021a, 2021b; Herrera-Camus et al. 2021). It is worth noting that, in the case of A1689-zD1, the [C II] line profile may also be composed of narrow+broad components, as with the [O III] line, but degenerate with a single Gaussian profile. In this scenario, we are observing multiphase outflows in an early galaxy including both cold neutral and hot ionized gases.
We also note that outflows may also be capable of forming the extended dust structure potentially observed around A1689-zD1 (Section 3.2). Though the outflow may also contribute to destruction of dust, the existence of dust from ∼20 kpc to a few Mpc from galaxies has been confirmed from the reddening measurements of background objects (e.g., Ménard et al. 2010; Peek et al. 2015; Tumlinson et al. 2017). Fujimoto et al. (2020) report that 1 out of 23 isolated ALPINE galaxies shows a dust continuum profile similar to the [C II] emission, despite compact UV continuum. Interestingly, the object also shows a large Lyα EW even with the large amount of dust, and the authors argue for a potential mechanism of outflows dominated by hot ionized gas, which may serve to eject dust out to the CGM, keep the dust hot, and allow for the Lyα line to escape from the system. If the hot-mode outflows last only for a short timescale, this extended dust continuum structure could be somewhat rare in the early universe.
4.2.3. Tidal Stripping
As previously noted, A1689-zD1 is a clumpy object with multiple spatial/spectral components necessary to fit the [C II] and [O III] data cubes (Wong et al. 2022, K. Knudsen et al. 2022, in preparation). This raises the possibility that A1689-zD1 is actually comprised of multiple merging systems, and that the [C II] halo arises from the dissipation of mechanical energy via turbulence and shocks associated with tidal stripping. Such tidal effects would also be expected to eject material from the galaxy, including dust and old stars that would contribute to the ISRF in the CGM. These effects have been explored by Ginolfi et al. (2020b), who observe [C II] emission from circumgalactic gas structure around a massive ( ∼ 1010 M⊙), merging galaxy system at z = 4.57 in the center of a protocluster environment.
The close projected distances ( ∼ 05) of the different components in A1689-zD1 make it difficult to robustly distinguish between multiple, merging systems and a single complex, clumpy object. Moreover, A1689-zD1 is a sub-L*, low-mass ( ∼109
M⊙) galaxy, and the galaxy overdensity is not identified around A1689-zD1. As such, the situation may be different from a system of massive merging galaxies in the center of a high-redshift protocluster. Nevertheless, minor mergers have been predicted to drive significant kinematic disruption of gas in low-mass galaxies (see, e.g., the simulations of Kohandel et al. 2020), which could power [C II] emission. While further analysis of the kinematics of the object is needed in order to explore the merging system scenario (and will be explored in K. Knudsen et al. 2022, in preparation), it should be noted that tidal disruption in a merging system may play a role in producing the [C II] halo around A1689-zD1.
4.2.4. Summary: Origin of the [C ii] Halo
Our observational results suggest the presence of outflows, dominated by hot, ionized gas, in and around the star-forming core of A1689-zD1. In conjunction with the analytical models of Pizzati et al. (2020), our results imply that the [C II] halo likely formed via some combination of a) hot, ionized outflows and subsequent catastrophic cooling, or b) past outflow activity including the cool, neutral gas. Our results do not rule out the possibility that the CGM was pre-enriched by past outflow/merger activity and the observed [C II] halo is powered by strong ionizing photons from the galaxy core or extended H II regions. In any case, our results support the connection between the [C II] halo and the outflow history of the galaxy, and highlight the importance of considering the multiphase nature of galactic outflows.
However, it remains a possibility that the central and secondary components trace a kinematically complex bulge and a diffuse disk of the galaxy, respectively, where the latter is missed in previous HST observations of the rest-frame UV continuum. Rest-frame optical observations with JWST, which have been scheduled for Cycle 1, will map out the entire stellar distribution and help to answer definitively whether the extended [C II] emission is part of the galaxy or the CGM.
4.3. Origin of High L[O III]/L[C II] Ratios at z > 6
Recent evidence from ALMA has shown that high-redshift galaxies tend to have (∼3–10 times) higher L[O III]/L[C II] ratios than local galaxies (Hashimoto et al. 2019; Bakx et al. 2020; Carniani et al. 2020; Harikane et al. 2020). This may suggest different properties of the ISM at high-z, such as a higher ionization parameter, PDR deficit, higher burstiness parameter (Vallini et al. 2021), or lower C/O abundance ratio than solar (Becker et al. 2011; Katz et al. 2022).
Using the best sensitivity [O III] and [C II] maps so far obtained at this redshift, we have shown the existence of the [C II] halo and [O III] outflows in A1689-zD1, both of which may be additional factors to explain the high L[O III]/L[C II] ratios in z > 6 galaxies. Although our results rely on one galaxy, A1689-zD1 is just one of the abundant population of sub-L* galaxies whose properties may be a representative picture in the Epoch of Reionization. Therefore, we discuss here the contribution from the [C II] halo and [O III] outflows to the L[O III]/L[C II] ratio in A1689-zD1. A full analysis of the spatially resolved L[O III]/L[C II] ratio in A1689-zD1 will be presented in K. Knudsen et al. 2022, (in preparation).
4.3.1. Contribution from [C ii] Halo
One possible explanation for the high L[O III]/L[C II] ratios in z > 6 galaxies is the systematic existence of, and failure to account for, the extended [C II] emission on CGM scales. This has been investigated by Carniani et al. (2020), who find that correcting for surface brightness dimming in high-resolution, low-S/N observations can lower the L[O III]/L[C II] ratios by a factor of ∼2.
Figure 11 shows the L[O III]/L[C II] ratio as a function of SFR for A1689-zD1, in multiple circular apertures of different sizes. These different apertures each encompass different amounts of the [C II] halo. We additionally show the nine galaxies studied in Harikane et al. (2020). By accounting for the full extent of the [C II] halo ( <35), we find a lower L[O III]/L[C II] ratio of ∼2, compared to ∼4 for the central region of the galaxy ( <0
5). That is, accounting for extended [C II] emission can help to reconcile the high L[O III]/L[C II] ratios with lower ratios found in local galaxies, but cannot bring these values completely in line with local relations.
Figure 11. The L[O III]/L[C II] ratio vs. SFR for A1689-zD1 in multiple apertures (red diamonds), compared to the literature sample from Harikane et al. (2020) (black circles). The SFR values are computed as described in Section 3.4. While the L[O III]/L[C II] ratio is lowered by a factor of ∼2 by accounting for the full [C II] halo, this cannot account for the entirety of the high L[O III]/L[C II] ratios in high-redshift galaxies.
Download figure:
Standard image High-resolution image4.3.2. Contribution from [O iii] Outflows
Another possible explanation for the high L[O III]/L[C II] ratio, at least in the case of A1689-zD1, is the contribution of the "broad-wing" feature in the [O III] spectrum to the total [O III] luminosity. If this emission indeed corresponds to hot, ionized gas in outflows, rather than H II regions around massive stars, we would expect to see a higher L[O III]/L[C II] ratio here than in local universe galaxies not undergoing this outflow.
To examine this, we analyze the [O III] line profile in A1689-zD1. Since the single Gaussian fitting in [O III] shows significant residuals (Section 3.5), we perform two-component (narrow+broad) Gaussian fitting to the [O III] line profile. Figure 12 shows the [O III] line profile in the same central aperture as Figure 8, with a two-component Gaussian fit applied. We see that the broad component dominates the overall spectrum, with the narrow component responsible for only about 13% of the total luminosity.
Figure 12. Top: [O III] line profile in the central 025 region as in Section 3.5, but with a double Gaussian fit applied. Bottom: the L[O III]/L[C II] ratio vs. SFR for A1689-zD1, as in Figure 11, for the total [O III] luminosity as well as only the narrow component.
Download figure:
Standard image High-resolution imageThe bottom panel of Figure 12 shows how the L[O III]/L[C II] ratio changes when calculated using the total [O III] luminosity (narrow+broad) versus the luminosity of only the narrow component. We use the same SFR for both points, derived from the UV+FIR emission as described in Section 3.4. We see that, while this central region shows an overall L[O III]/L[C II] ratio of ∼3, the narrow component alone shows a ratio of ∼0.4. Assuming the low-Z dwarf galaxies in the De Looze et al. (2014) sample do not show significant outflows, this result suggests that outflowing gas in the hot, ionized phase could be a factor contributing to the high L[O III]/L[C II] ratios in high-z galaxies. While the outflow signature in the [O III] line is not universally identified around the high-z galaxies with high L[O III]/L[C II] ratios, it has been tentatively detected in the galaxy NB1006 (Inoue et al. 2016). Importantly, the S/N of the [O III] line in our data is at least 4–8 times greater than in previous studies (e.g., Harikane et al. 2020), which could explain the nondetections of the [O III] outflow signature.
5. Summary
In this paper, we present deep ALMA observations of the strongly lensed, z = 7.13 galaxy A1689-zD1, focusing on the [C II] 158 μm and [O III] 88 μm fine-structure lines as well as the FIR continuum. Utilizing archival HST data to probe the UV continuum, we compare the morphology and spatial extent of each wavelength of emission in order to examine the multiphase nature of the ISM. Our major findings are summarized below.
- 1.We detect both [C II] and [O III] emission at high S/N, as well as dust continuum emission in both bands. The [C II] emission is extended out to r ∼ 12 kpc from the galaxy center, while the [O III] and UV continuum emission extend only to r ∼ 4 kpc.
- 2.Despite low S/N, we find that the dust continuum follows a profile similar to that of the [C II] line and is more extended than the [O III] line or UV continuum.
- 3.Utilizing multiwavelength FIR continuum data, we measure a dust temperature of
with a dust emissivity index of
. We find slight (∼1–2σ) differences in the dust continuum radial profiles, indicating that dust temperature might vary in space; in fact, we find tentative indications that the temperature decreases from ∼50 K at the galaxy center to ∼30 K at r ∼ 1
5.
- 4.The spatially resolved Σ[O III]–ΣSFR relation shows an excess of [O III] emission in the galaxy core as compared to local dwarf galaxies. By contrast, the Σ[C II]–ΣSFR relations shows a slight deficit of [C II] emission in the galaxy core and a slight excess in the outskirts, consistent with z ∼ 2–4 SFGs.
- 5.The central, high-ΣSFR region shows a [C II] line width of 252 ± 21 km s−1 and an [O III] line width of 182 ± 31 km s−1, indicative of the larger spatial extent of the [C II] line along the line of sight. We find evidence for ongoing hot ionized outflows traced by [O III], with ∼ 5σ residuals from a single Gaussian model at ±100–400 km s−1. This suggests a possible origin of the [C II] halo from hot ionized outflows and subsequent cooling.
Much remains uncertain about the ISM/CGM properties, kinematics, and evolutionary stage of A1689-zD1. As A1689-zD1 will be observed in JWST Cycle 1 (Prop. ID 1840, PI: Alvarez-Marquez), future work will be able to take full advantage of imaging from the rest-frame UV to the optical. The strongly gravitationally lensed nature of this galaxy makes it a prime target with which to constrain the processes of early galaxy evolution.
The authors are grateful to Francesca Rizzo, Andrea Pallottini, Livia Vallini, Mahsa Kohandel, Andrea Ferrara, David Frayer, and Harley Katz for helpful discussions regarding the physical interpretations of the extended [C II] emission around A1689-zD1. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2015.1.01406.S, ADS/JAO.ALMA#2017.1.00775.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. This work was supported by the U.S. NSF under grant OISE-2005578 (DAWN-IRES). This project has recieved funding from the European Union's Horizon 2020 research and innovation program under the Marie Sklodowska-Curie grant agreement No. 847523 "INTERACTIONS." H.B.A. acknowledges the support and collaboration of the mentors and fellow students involved in the 2021 DAWN-IRES program. S.F. and D.W. acknowledge support from the European Research Council (ERC) Consolidator Grant funding scheme (project ConTExt, grant No. 648179) and Independent Research Fund Denmark grant DFF7014-00017. The Cosmic Dawn Center is funded by the Danish National Research Foundation under grant No. 140. K.K. acknowledges support from the Knut and Alice Wallenberg Foundation. T.H. was supported by Leading Initiative for Excellent Young Researchers, MEXT, Japan (HJH02007) and KAKENHI (20K22358). A.K.I. is supported by NAOJ ALMA Scientific Research Grant Code 2020-16B. M.J.M. acknowledges the support of the National Science Centre, Poland through the SONATA BIS grant 2018/30/E/ST9/00208. Much of this work was conducted on the ancestral land of the Meskwaki, Sauk, and Ioway Peoples.
Appendix: Supplementary Figures
Figure A1 shows a map of the magnification μ across the field. While the magnification is highest in the southwest corner of the field, and lowest in the northeast, it is quite uniform across the source.
Figure A1. Map of the lensing magnification of A1689-zD1. Contours indicate the [C II] 158 μm line intensity, and the background is colored according to the magnification at that point. The magnification is highest in the southwest corner of the field, but largely uniform (within ∼5%) over the galaxy.
Download figure:
Standard image High-resolution imageWe show in Figure A2 the 2D maps of each emission, in the image plane (as in Figure 1) as well as reconstructed into the source plane, with a uv-taper of 05 (0
7 for band 6) applied. Much of the analysis in this paper is conducted on these uv-tapered maps, to optimize the S/N of the extended, diffuse emission. Moreover, we derive radial profiles, Gini coefficients, and the dust temperature profile from the uv-tapered source plane maps, as they provide a direct measurement of the radial distance, corrected for the lens deflection.
Figure A2. Top: maps of each emission, as in Figure 1, but with a uv-taper of 05 (0
7) applied to band 8 (band 6) maps. Bottom: maps of each emission reconstructed into the source plane.
Download figure:
Standard image High-resolution imageFootnotes
- 14
We note that drop off in the dust continuum profile at 90 μm is not due to a lack of sensitivity. The peak S/N is about the same for the band 6 and band 8 dust continuum maps (∼26 versus ∼23, respectively). Moreover, our band 8 data set is more sensitive to large-scale structure, as it samples more short baselines (Figure 1).
- 15
- 16
We consider the temperature undetermined if either one or both of the band 6 or 8 continuum maps were not detected at S/N > 2, or if the resulting fractional uncertainty on the dust temperature exceeds 30%.
- 17
We determine the FWHM via linear interpolation between the data above and below 50%. We determine the error on the FWHM in a Monte Carlo fashion, by varying each data point in the line spectrum according to the uncertainty.
- 18