Abstract
We present Atacama Large Millimeter/submillimeter Array observations of CO lines and dust continuum emission of the source RCSGA 032727–132609, a young z = 1.7 low-metallicity starburst galaxy. The CO(3–2) and CO(6–5) lines and continuum at rest-frame 450 μm are detected and show a resolved structure in the image plane. We use the corresponding lensing model to obtain a source plane reconstruction of the detected emissions revealing an intrinsic flux density of μJy and intrinsic CO luminosities and . We used the resolved properties in the source plane to obtain molecular gas and star formation rate surface densities of and , respectively. The intrinsic properties of RCSGA 032727–132609 show an enhanced star formation activity compared to local spiral galaxies with similar molecular gas densities, supporting the ongoing merger–starburst phase scenario. RCSGA 032727–132609 also appears to be a low-density starburst galaxy similar to local blue compact dwarf galaxies, which have been suggested as local analogs to high-redshift low-metallicity starburst systems. Finally, the CO excitation level in the galaxy is consistent with having the peak at , with a higher excitation concentrated in the star-forming clumps.
Export citation and abstract BibTeX RIS
1. Introduction
The study of the distribution of molecular hydrogen (H2) and star formation rate (SFR) in high-redshift galaxies has been a prolific field for the past 20 years. The CO emission lines, good tracers of the molecular gas (Omont 2007), and the cold dust continuum emission have been resolved in a number of galaxies at (Carilli & Walter 2013; Hodge et al. 2015). These observations have shown a picture where the SFR is strongly linked to the distribution of molecular gas, indicating that two mechanisms for star formation could be in place, the quiescent, normal phase, and the starburst phase (Daddi et al. 2010; Genzel et al. 2010, 2015). Most galaxies at all redshifts appear to be in the main-sequence (MS) phase (e.g. Brinchmann et al. 2004), where the SFR is proportional to the stellar and molecular gas masses. These galaxies would be producing stars following the standard mechanism (still not fully understood) of star formation. The starburst phase corresponds to those galaxies that, for a given stellar and/or molecular gas mass, have SFRs of at least 10 times higher than that of MS galaxies. The starburst galaxies show a more efficient process of star formation, expected to be triggered by other factors such as galaxy interactions or particular environmental conditions (Dekel et al. 2009; Engel et al. 2010).
It has been found that starburst galaxies at high redshift () show similar star formation conditions and efficiencies as those observed in particular regions of local (ultra)luminous infrared galaxies (Carilli & Walter 2013). The unresolved nature of most of the high-redshift observations does not allow us to see whether the high-redshift starburst galaxies have higher efficiencies across the whole galaxy or the star formation is dominated by small, highly efficient regions as in the local galaxies. The need for higher-resolution observations is clear if we want to understand the starburst process and its conditions.
The best way to resolve the emission in galaxies without increasing substantially the observation time is by targeting gravitationally lensed galaxies, particularly the spatially resolved bright gravitational arcs. Observations of several lensed sources already show that different star-forming regions show different star formation conditions (Hezaveh et al. 2013; ALMA Partnership et al. 2015; Hatsukade et al. 2015; Spilker et al. 2016). In this Letter, we present the CO and continuum emission of the second-brightest optical giant arc discovered to date, the arc RCSGA 032727–132609 (hereafter RCS0327), strongly lensed by the foreground galaxy cluster RCS2 032727–132623 at z = 0.564 (Wuyts et al. 2010), to reveal the star formation process in a starburst at z = 1.7. The kinematic analysis of Hα on RCS0327 strongly suggests an interaction consistent with a merger of galaxies (Wuyts et al. 2014). The SFRs measured in individual clumps and across the galaxy fall well above the MS relation for galaxies at z = 1.7 and are consistent with a young low-metallicity () starburst enhanced by the ongoing merger (Whitaker et al. 2014; Wuyts et al. 2014). Furthermore, Bordoloi et al. (2016) measured resolved galactic winds in RCS0327, showing that the outflows are comparable to those observed in local starbursts.
Throughout this Letter, we adopt a cosmology with km s−1 Mpc−1, , and (Planck Collaboration et al. 2016).
2. Observations
The target RCS0327 was observed as part of the ALMA project 2015.1.00920.S, which aimed to detect the emission lines CO(3–2), CO(6–5), and CO(8–7) and the underlying continuum. At the redshift for RCS0327 of (Wuyts et al. 2014) the lines are centered at 127.895 GHz for CO(3–2), 255.746 GHz for CO(6–5), and 340.934 GHz for CO(8–7). The observations of the lines CO(6–5) and CO(3–2) were carried out during 2016 January 1 and 10. The observations of CO(8–7), unfortunately, were not completed during cycle 3.
In each observation, two overlapping spectral windows (SPWs) were placed to detect the line and two were placed to detect the continuum. The passband and amplitude calibrator for both observations was J0423–0120, with J0336–1302 being used as phase calibrator. The reduction of the data was performed using the scripts provided by ALMA and used the pipeline version for cycle 3 and casa (McMullin et al. 2007) version 4.5.1.
Both emission lines were observed using the array configuration C36-1. The imaging of the calibrated data was done using the casa task clean. The data for CO(3–2) were imaged using natural weighting returning a synthesized beam size of 254 × 185 and position angle of −8568. The natural weighting synthesized beam for the CO(6–5) observations was 149 × 098 at a position angle of 7975. To achieve a beam size similar to CO(3–2), a taper is applied to the CO(6–5) data to obtain a synthesized beam of 216 × 208 at a position angle of 5964.
Data cubes using a spectral resolution of 20 km s−1 were created for each of the emission lines. The continuum emission is estimated by using the line free channels and then is subtracted from the uv data using uvcontsub, and new continuum subtracted cubes are created using the same procedure as before. The continuum and line images were interactively cleaned by manually masking out the emission.
3. Results
The CO(3–2), CO(6–5), and rest-frame 450 μm continuum emission are detected in the brightest regions of the arc (Figure 1). The continuum emission at rest-frame 830 μm is also detected but with lower significance. To find the CO(3–2) total flux we extract the spectra on the positions e1, e2, and u1+u2, clumps identified by Sharon et al. (2012) that are near the peaks of the observed emission, which are plotted in the bottom panel of Figure 1. A small astrometric offset was applied to the Hubble Space Telescope (HST) coordinates to match the 830 μm continuum detection of the brightest cluster galaxy to its near-infrared counterpart. We find that the region going from 127.874 to 127.976 GHz ( km s−1) provides a good range for the total emission of the line. The CO(3–2) line flux measured in the image plane for the section of the arc going from e1 to e2 (corresponding to lensed images 1 and 2) is Jy km s−1, while for CO(6–5) is Jy km s−1.
The spatially integrated continuum emission measured in the same region as CO(3–2) is and . A second emitting region is detected over the arc, outside the primary beam (PB), marked with a red circle in the middle panel of Figure 1. The continuum emission for that component is .
4. Discussion
4.1. Source Plane Reconstruction
Working with a galaxy that is magnified by strong gravitational lensing requires additional analysis to recover the galaxy's intrinsic properties. Knowledge of the lensing deflection field is required to reconstruct the galaxy in the source plane, and to account for the magnification—which can be highly spatially variable. Different galaxy regions will be stretched and magnified by a different factor, resulting in different physical scales over the arc. To fully understand our detected emission we need to take them to the source plane, where the physical scale is unique.
To do the source plane reconstruction for RCS0327, we used uvmcmcfit,10 which is a Python implementation to fit emission models to interferometric data in the uv plane (Bussmann et al. 2013, 2015; González-López et al. 2017). Exploring the emission in the uv plane should return the maximum amount of information from the observations and it has proven to be useful in revealing the source plane emission in bright SMGs discovered by the South Pole Telescope and observed with ALMA (Hezaveh et al. 2013; Spilker et al. 2016).
Uvmcmcfit can fit the source plane emission of a galaxy together with the lensing potential. The galaxy emission is fitted assuming a 2D elliptical Gaussian (special case of a Sérsic profile with n = 0.5), while the lensing potential is fitted by a singular isothermal ellipsoid (SIE). Recent high-resolution imaging of high-redshift SMGs have shown no strong preference between fitting the dust continuum emission with Gaussian or Sérsic profiles (Hodge et al. 2016; Spilker et al. 2016), supporting the usage of a simple 2d elliptical Gaussian function for the source model. We modified uvmcmcfit so that it uses a user provided lensing model deflection field allowing us to incorporate the detailed lensing model that is already available, based on HST imaging (Sharon et al. 2012) and fit the ALMA data only to constrain the properties of the source plane emission. During the fitting of the emission, the lensing model is held fixed and only the source plane emission model is allowed to vary. The model found by using high resolution HST observations of multiple images of strongly lensed galaxies at different redshifts together with the cluster members information outperforms in complexity and quality to the model we could find by fitting the lensing potential and source structure together using only the ALMA observations (Sharon et al. 2012).
For the case of CO(3–2), a single 2D elliptical Gaussian was needed to fit the observed emission. We measure an intrinsic flux of μJy (for the frequency range of km s−1) and an effective radius arcsec ( kpc). The flux weighted magnification value for CO(3–2) is and an intrinsic CO(3–2) luminosity of .
In the case of CO(6–5), a single component gives μJy (for the frequency range of km s−1) and an effective radius arcsec ( kpc). The flux weighted magnification value for CO(6–5) is and an intrinsic CO(6–5) luminosity of .
To fit the continuum emission at 450 μm, we needed two Gaussians in the source plane, as a single component was not enough to account for the whole observed flux outside the PB in Figure 1. The main component is well described by a Gaussian with μJy and arcsec ( kpc). The second component returns μJy and arcsec. The flux weighted magnification value for the main component of continuum emission is , while for the second component is .
In Figure 2, we present the observed emission for CO(3–2), CO(6–5), and continuum at 450 μm (left panels); the image plane representation of the best-fit model found for each case (middle panels); and the residual images obtained after subtracting the best model simulated visibilities from the observed ones (right panels). In all cases, the best-fit models appear to account for most of the observed emission.
Download figure:
Standard image High-resolution image4.2. Gas and Dust Distribution
We now use the CO(3–2) to estimate the amount and distribution of molecular gas in the galaxy. We first need to estimate the amount of CO(1–0) luminosity based on the observed CO(3–2) luminosity and use the CO conversion factor () to convert it to molecular gas mass (Bolatto et al. 2013; Carilli & Walter 2013). The CO excitation level depends mainly on the density and temperature of the gas, and it has been found that the excitation level differs for different source populations. Similar results have been found for the CO conversion factor, where a typical value of has been used for nuclear starbursts such as submillimeter galaxies (SMGs) and quasi-stellar objects (QSOs) and the Milky Way value of in MS high-redshift CSGs. Because of RCS0327 being cataloged as a starburst, to estimate the molecular mass we will use the values estimated for SMGs and starbursts at high redshift.
Assuming a CO excitation level valid for SMGs of (similar to ∼0.6 for MS CSGs; Carilli & Walter 2013), we obtain an intrinsic CO line luminosity of and . Using the effective radius measured for the CO(3-2) emission, we estimate the total area for the molecular gas surface of kpc2 and obtained an H2 surface density of .
In Figure 3, we present the source plane emission from the optical, CO(3–2), and continuum at 450 μm. The regions where the CO(3–2) and continuum emission are produced corresponds to clumps a–f (see Figure 4 in Sharon et al. 2012), which have a combined spectral energy distribution of (Wuyts et al. 2014). We use the estimated total in the same region derived by Sharon et al. (2012) as an upper limit to the total SFR produced by clumps and the interstellar medium (ISM) combined. The SFR in combination with the molecular gas mass give a depletion time of ∼40 Myr and an SFR surface density of .
Download figure:
Standard image High-resolution imageWe also derive resolved properties for individual regions plotted in Figure 3(b). These regions have a size of 025 (2.2 kpc) in the source plane, which given the variable magnification some will be larger than the beam projected at the source plane. We use a set of 50 source plane reconstructions taken from the fitting iterations to estimate the significance of the detection in each region. We also take the image plane noise map into the source plane to give a proper upper limit for the emission in the regions with . The results for all the regions are presented in Table 1.
Table 1. Flux Density Values Measured in the Source Plane
Region | a | a | b | |||
---|---|---|---|---|---|---|
(1) | (2) | (3) | (4) | (5) | (6) | (7) |
Main Componentc | 1.1 ± 0.2 | |||||
01 | 4.0 ± 0.9 | 6.7 ± 3.3 | 3.6 ± 3.4 | 0.5 ± 0.6 | 7.1 ± 3.5 | |
02 (e, f, u) | 6.9 ± 1.0 | 12.3 ± 1.8 | 16.3 ± 6.4 | 1.3 ± 0.6 | 13.0 ± 1.9 | |
03 | 1.4 ± 0.6 | 6.6 ± 2.1 | 0.0 ± 1.6 | <0.5 | 7.0 ± 2.2 | |
04 | 0.0 ± 0.3 | 1.2 ± 1.1 | 0.0 ± 1.6 | ⋯ | <2.3 | <0.10 |
05 | 0.0 ± 0.2 | 0.1 ± 0.7 | 0.0 ± 1.5 | ⋯ | <1.5 | <0.07 |
06 (t) | 0.0 ± 0.2 | 0.0 ± 0.7 | 0.0 ± 1.6 | ⋯ | <1.5 | <0.07 |
07 (t) | 0.0 ± 0.3 | 0.0 ± 0.7 | 0.0 ± 1.8 | ⋯ | <1.5 | <0.10 |
08 | 0.5 ± 0.3 | 5.0 ± 3.1 | 10.2 ± 9.7 | ⋯ | <6.6 | <0.10 |
09 (d) | 3.3 ± 1.0 | 12.8 ± 1.9 | 38.7 ± 9.9 | 3.0 ± 0.9 | 13.5 ± 2.0 | |
10 (b, c) | 3.0 ± 1.9 | 12.7 ± 2.7 | 0.0 ± 1.6 | <0.3 | 13.4 ± 2.9 | <0.64 |
11 (r) | 0.4 ± 0.3 | 4.3 ± 2.5 | 0.0 ± 1.6 | ⋯ | <5.3 | <0.10 |
12 (r) | 0.0 ± 0.3 | 0.5 ± 0.7 | 0.0 ± 1.7 | ⋯ | <1.5 | <0.10 |
13 (t) | 0.0 ± 0.3 | 0.0 ± 0.7 | 0.0 ± 1.8 | ⋯ | <1.5 | <0.10 |
14 | 0.0 ± 0.4 | 0.0 ± 0.7 | 0.0 ± 2.1 | ⋯ | <1.5 | <0.13 |
15 | 0.0 ± 0.3 | 0.0 ± 0.8 | 0.0 ± 1.9 | ⋯ | <1.7 | <0.10 |
16 | 0.3 ± 0.3 | 0.1 ± 0.7 | 0.0 ± 1.8 | ⋯ | <1.5 | <0.10 |
17 (a) | 0.9 ± 0.6 | 0.2 ± 0.7 | 0.0 ± 1.8 | ⋯ | <1.5 | <0.20 |
18 (r) | 0.7 ± 0.7 | 0.1 ± 0.7 | 0.0 ± 1.9 | ⋯ | <1.5 | <0.24 |
19 (r) | 0.1 ± 0.3 | 0.0 ± 0.7 | 0.0 ± 1.9 | ⋯ | <1.5 | <0.10 |
20 | 0.0 ± 0.3 | 0.0 ± 0.7 | 0.0 ± 2.1 | ⋯ | <1.5 | <0.10 |
21 | 0.2 ± 0.5 | 0.0 ± 0.7 | 0.0 ± 2.8 | ⋯ | <1.5 | <0.17 |
22 | 0.0 ± 0.3 | 0.0 ± 0.8 | 0.0 ± 2.1 | ⋯ | <1.7 | <0.10 |
23 | 0.0 ± 0.3 | 0.0 ± 0.8 | 0.0 ± 2.0 | ⋯ | <1.7 | <0.10 |
24 | 0.0 ± 0.3 | 0.0 ± 0.7 | 0.0 ± 2.0 | ⋯ | <1.5 | <0.10 |
25 | 0.1 ± 0.3 | 0.0 ± 0.7 | 0.0 ± 2.1 | ⋯ | <1.5 | <0.10 |
26 (s) | 0.1 ± 0.3 | 0.0 ± 0.7 | 0.0 ± 2.2 | ⋯ | <1.5 | <0.10 |
27 | 0.1 ± 0.4 | 0.0 ± 0.7 | 0.0 ± 2.5 | ⋯ | <1.5 | <0.13 |
28 | 0.4 ± 0.5 | 0.0 ± 0.8 | 0.0 ± 3.1 | ⋯ | <1.7 | <0.17 |
Notes. Upper limits correspond to .
aFor the frequency range of km s−1. bUsing . cAnd source plane emission knots (Sharon et al. 2012).Download table as: ASCIITypeset image
In Figure 4, we present the resolved properties for the total emission of RCS0327 and for each of the regions described above. We compare these results with those obtained for Milky Way molecular clouds (Heiderman et al. 2010; Evans et al. 2014), local spirals (Kennicutt 1998; Bigiel et al. 2010), local starbursts (Kennicutt 1998), local blue compact dwarfs galaxies (BCDs; Amorín et al. 2016), low-redshift dusty normal star-forming galaxies (VALES; Villanueva et al. 2017), z = 1–3 star-forming galaxies (SFGs; Genzel et al. 2010; Freundlich et al. 2013; Tacconi et al. 2013), SMGs (Bothwell et al. 2010), and the lensed SMGs observed in high resolution by ALMA SDP.81 (Hatsukade et al. 2015). We notice that RCS0327 falls above the relation found by Daddi et al. (2010) for MS and starburst galaxies (blue dashed lines), supporting the starburst nature of RCS0327. We point out that the position of RCS0327 on the diagram depends strongly on the assumed value for , but even when using the Milky value of , the galaxy would at most move to be on top of the starburst relation.
Download figure:
Standard image High-resolution imageFor our estimate of the molecular mass we see that RCS0327 properties are similar to the local BCDs, which are low-metallicity starbursts galaxies showing higher star-forming efficiencies when compared to normal disk galaxies (Hunt et al. 2015; Amorín et al. 2016). BCDs have already been found to work well as local analogs to similar redshift low-metallicity starbursts based on properties derived using optical spectroscopy (Brammer et al. 2012). Based on the latter, we can use the –metallicity relation found for BCDs to estimate an value for RCS0327. The relation presented by Amorín et al. (2016) returns a value of for RCS0327, consistent with the value given by other relations (Hunt et al. 2015). The new value would increase the estimate molecular gas mass for RCS0327 putting it near the MS relation with a depletion time of ∼1 Gyr, in the same region as the z = 1–3 SFGs.
4.3. CO Excitation Level
We can use the detected CO(3–2) and CO(6–5) emission lines to constrain the CO excitation level in RCS0327. The total intrinsic flux densities return a fraction of . This value is consistent with having the peak at , similar to some CO excitation levels measured on SMGs and lower than the excitation level measured for QSOs (Carilli & Walter 2013).
We can use the same method presented in the previous section to obtain resolved measurements of in the regions presented in Figure 3. We have five regions (see Table 1) with detections of CO(3–2) and constraints on CO(6–5). The ratios go from in region 10 to in region 9. Our results are consistent with those found for the high-z disk galaxy presented by Bournaud et al. (2015), where the CO excitation level is for the main star-forming clumps and for the inter-clump gas. The warmer and denser gas associated with the main clumps allows for a higher CO excitation level when compared to the more extended gas. In the case of RCS0327, the main star-forming clumps identified by Sharon et al. (2012) and Wuyts et al. (2014) correspond to the regions 2 and 9, which have ratios of 1.3 ± 0.6 and 3.0 ± 0.9. Regions 1, 3, and 10 show values, consistent with the inter-clump gas of the simulations (Figure 3).
5. Conclusion
We have presented ALMA observations of the emission lines CO(3–2), CO(6–5), and 450 μm rest-frame continuum emission in the z = 1.7 young low-metallicity starburst strongly lensed galaxy RCS0327. The source plane reconstruction of the detected emission reveals that the molecular gas, traced by CO(3–2), is located on top of the star-forming clumps showing the gas reservoir fueling the ongoing star-forming process. The continuum dust emission follows a similar angular extension and distribution as the molecular gas but also extending toward parts of the galaxy not as bright in CO(3–2), showing a clear spatial offset. The molecular gas and SFR surface density show that RCS0327 is a low-density starburst, similar to local BCDs and probably triggered by an ongoing merger. These results support the scenario where BCDs are identified as a local counterpart to the low-metallicity starburst galaxies at high redshift.
The detected CO(3–2) and CO(6–5) return a CO excitation level consistent with having the peak at at large scales. The total excitation appears to be the result of the combination of higher-excitation regions near the star-forming clumps and lower-excitation regions over the more extended gas phase. This is one of the first times where the CO excitation level is resolved in detail on a galaxy at high redshift.
We have shown that giant gravitational arcs offer an excellent opportunity to resolve in detail the different phases of the ISM in the cases where a good lensing model is in hand. Our coarse observations already show that RCS0327 is not well described by a single-mode star-forming galaxy, showing different CO excitation levels, molecular gas reservoirs, and dust obscuration across . Future high angular resolution observations and the extension to other bright arcs will take us one step closer to understanding the star formation process in galaxies at high redshift.
This Letter makes use of the following ALMA data: ADS/JAO.ALMA#2015.1.00920.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada) and NSC and ASIAA (Taiwan), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. This research has been supported by CONICYT-Chile grant Basal-CATA PFB-06/2007, FONDECYT Regular 1141218, and ALMA-CONICYT project 31160033. M.A. acknowledges partial support from FONDECYT through grant 114009.