EDP Sciences logo
Subscriber Authentication Point
Search
Menu
Home All issues Volume 631 (November 2019) A&A, 631 (2019) A167 Full HTML
Free Access
Issue
A&A
Volume 631, November 2019
Article Number A167
Number of page(s) 8
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/201936169
Published online 15 November 2019

© ESO 2019

1. Introduction

The interaction between the interstellar medium (ISM) and stars is key to understanding galaxy formation and evolution. Stars form from dense, cool gas, and have a profound influence on their surrounding ISM by heating and ionizing the surrounding gas, which, in turn, inhibits further star formation. Stellar winds and supernovae enrich the ISM with newly formed metals, influencing the gas cooling rates, and hence the star-formation efficiency. A detailed knowledge of the physical parameters of the ISM, such as the metallicity, ionization state, density, and temperature is therefore essential to study galaxy evolution. The main diagnostics available over a wide range of redshifts are dust and stellar continuum emission as well as spectral emission and absorption lines.

The most frequently used methods to determine the star-formation rate (SFR) include direct far-ultraviolet (far-UV) emission from young stars, far-infrared (far-IR) dust continuum emission from dust heated by young stars, and recombination lines (e.g. Hα) ionized by the emission field of the young stellar population (see Kennicutt & Evans 2012, for a review). To determine the gas metallicity, rest-frame optical emission lines are traditionally the most frequently used tracers (e.g. Maiolino et al. 2008). The density, ionization parameter, and temperature of the gas can be derived from line ratios of fine structure lines (FSL; e.g. Villar-Martin et al. 1997; De Breuck et al. 2000). One common disadvantage of these mostly optical tracers is that they can be severely affected by dust obscuration. Such problems can be avoided, or at least minimized by observing rest-frame far-IR FSL. While some galaxies such as Arp 220 are optically thick out to λrest = 240 μm (Rangwala et al. 2011), the extinction corrections to far+IR FSL are at most ∼50% for the shortest wavelength lines (Uzgil et al. 2016). Moreover, these far-IR FSL are often brighter than their optical counterparts (e.g. Palay et al. 2012; Moriwaki et al. 2018); the brightest FSL of carbon, oxygen, and nitrogen are the major coolants of major phases in the ISM, and can each carry up to ∼1% of the total far-IR luminosity.

The far-IR FSL originate both from photo-dissociation regions (PDR) and from HII regions, depending if their ionization potentials are lower or higher than that of hydrogen (13.6 eV), respectively. All atomic FSL ([OI] 63, 145 μm and [CI] 370, 609 μm) mainly originate from the neutral PDR. But also the ionic [CII] 158 μm line has an ionization potential of 11.26 eV, and therefore originates from both the PDR and HII regions. Any diagnostic line ratios including such multiple-origin lines should therefore invoke a self-consistent model that treats both PDR and HII regions in pressure equilibrium (e.g. Abel et al. 2005; Nagao et al. 2012; Cormier et al. 2019). The relative contributions to composite lines such as [CII] can also be determined observationally by comparing the line flux with that from another FSL with a similar critical density that originates fully in the HII region. One such line is [NII] 205 μm, which has n crit e − $ n_{{\rm{crit}}}^{{e^ - }} $ = 48 cm−3, compared to n crit e − $ n_{{\rm{crit}}}^{{e^ - }} $ = 50 cm−3 for [CII]1 at the same temperature (assuming electrons are the main collision partners, see Fernández-Ontiveros et al. 2016). This technique has been used on individual galaxies at high redshift (e.g. Decarli et al. 2014), as well as on samples of nearby galaxies (Kamenetzky et al. 2014; Díaz-Santos et al. 2017; Croxall et al. 2017), yielding [CII] PDR fractions of ∼75%. The latter authors did find a clear dependence on metallicity in this PDR/HII fraction. If the fraction can be found by other means, this dependence can also be turned around to determine the gas metallicity independently of classical optical diagnostics (e.g. Nagao et al. 2011, 2012; De Breuck et al. 2014; Béthermin et al. 2016; Pereira-Santaella et al. 2017). An easier solution is to consider only HII-dominated FSL to determine the physical parameters. For example the [NII] 122 μm/[NII] 205 μm and [OIII] 52 μm/[OIII] 88 μm ratios are powerful density tracers in the range 1–103 and 10–104.5 cm−3, respectively (Fernández-Ontiveros et al. 2016). The most promising metallicity tracers are the [OIII] 88 μm/[NII] 122 μm or [OIII] 88 μm/[NII] 205 μm ratios (Pereira-Santaella et al. 2017). Using Herschel Spectral and Photometric Imaging Receiver (SPIRE) spectroscopy of three gravitationally lensed dusty star-forming galaxies (DSFG), Rigopoulou et al. (2018) use the [OIII] 88 μm/[NII] 122 μm ratio to argue that these DSFGs must already have close to solar metallicities. Tadaki et al. (2019) confirm a similar result on an unlensed DSFG observed with the Atacama Large Millimeter/submillimeter Array (ALMA).

The applicability of mid-IR and far-IR FSL as diagnostic tools of the ISM has received a boost thanks to the publication of samples of nearby galaxies observed with the infrared space observatory (ISO) and Herschel (e.g. Brauher et al. 2008; Farrah et al. 2013; Sargsyan et al. 2014; Kamenetzky et al. 2014; Cormier et al. 2015; Cigan et al. 2016; Herrera-Camus et al. 2016; Fernández-Ontiveros et al. 2016; Zhao et al. 2016; Díaz-Santos et al. 2017; Zhang et al. 2018). At high redshift (z ≳ 1), the far-IR FSL conveniently shift into the (sub)millimetre atmospheric windows. The most popular line is clearly [CII] 158 μm, followed by the [CI] 370, 609 μm lines (e.g. Walter et al. 2011; Bothwell et al. 2017). Both single dish submillimetre telescopes and interferometers have also detected the [NII] 122 μm and 205 μm lines at high redshift (Ferkinhoff et al. 2011, 2015; Combes et al. 2012; Nagao et al. 2012; Decarli et al. 2012, 2014; Béthermin et al. 2016; Pavesi et al. 2016; Lu et al. 2017; Tadaki et al. 2019; Novak et al. 2019). After a slow start, the [OIII] 88 μm line is quickly becoming a popular line to confirm redshifts of galaxies in the epoch of reionization (z ≳ 6), where it shifts into the submillimetre atmospheric windows below 500 GHz (Ferkinhoff et al. 2010; Inoue et al. 2016; Carniani et al. 2017; Marrone et al. 2018; Vishwas et al. 2018; Hashimoto et al. 2018, 2019; Walter et al. 2018; Tamura et al. 2019; Tadaki et al. 2019; Novak et al. 2019). Deep Herschel/SPIRE spectroscopy has also revealed a number of FSL detections, either in individual objects (Valtchanov et al. 2011; Coppin et al. 2012; George et al. 2013; Uzgil et al. 2016; Rigopoulou et al. 2018; Zhang et al. 2018), or in stacked spectra (Wardlow et al. 2017; Wilson et al. 2017; Zhang et al. 2018). These include the only detections of the [OI] 63 μm line at high redshift reported thus far. This [OI]3P2−3P1 line is arguably the best tracer for the star-forming gas, as it traces the very dense ( n crit H $ n^{\mathrm{H}}_{\mathrm{crit}} $ = 5×105 cm−3) neutral gas (one caveat being the frequeny presence of self-absorption observed in local ultra-luminous IR galaxies, Rosenberg et al. 2015). Like the [CII] line, the [OI] 63 μm line also shows a “deficit” in the most luminous far-IR sources, though with a higher scatter (Graciá-Carpio et al. 2011; Cormier et al. 2015; Díaz-Santos et al. 2017). Surprisingly, this bright FSL has not been frequently observed with ALMA, probably because it is only observable in the highest frequency bands. At least as surprising is that the fainter, but more accessible [OI]3P1−3P0 transition at λrest = 145 μm ( n crit H $ n^{H}_{\mathrm{crit}} $ = 9.5×104 cm−3) has only recently been detected at high redshifts (Novak et al. 2019 report a tentative detection in a z = 7.5 quasar). Also in nearby galaxies, this [OI] 145 μm line has not been observed very frequently as in most cases, it is fainter than the nearby [CII] 158 μm line. After initial detections with ISO (Malhotra et al. 2001; Brauher et al. 2008), Herschel has now detected [OI] 145 μm in significant samples of nearby galaxies (Spinoglio et al. 2015; Cormier et al. 2015; Fernández-Ontiveros et al. 2016; Herrera-Camus et al. 2018), and recently in a z = 6.5 lensed quasar (Yang et al. 2019).

In this paper, we report the second such [OI] 145 μm detection in SPT-S J041839−4751.8 (hereinafter SPT 0418−47) at z = 4.2248 (Weiß et al. 2013). With a lensing magnification μ = 32.7 ± 2.7 (Spilker et al. 2016), SPT 0418−47 is one of the most strongly lensed DSFGs from our sample selected from the South Pole Telescope, making its FSL accessible to the Atacama Pathfinder EXperiment (APEX). We also report the detection of the [OIII] 88 μm line, and use the five detected FSL in SPT 0418−47 to constrain the metallicity in this DSFG. Throughout this paper, we assume a ΛCDM cosmology with H0 = 67.8 s−1 Mpc−1, Ωm = 0.308, and ΩΛ = 0.692 (Planck Collaboration XIII 2016). At z = 4.2248, this corresponds to a luminosity distance DL = 39.150 Gpc and a scale of 7 kpc/″.

2. Observations

2.1. New data from APEX

We observed the [OI] 145 μm (νrest = 2060.07 GHz) using the First-Light APEX Submillimeter Heterodyne instrument (FLASH; Klein et al. 2014) on the APEX telescope (Güsten et al. 2006). The observations were taken under the Max Planck project 093.F-9512, from 2014 May 14 to 20. The total on-source integration time was 7.3 h, with precipitable water vapour (PWV) in the range 0.4–0.8 mm. The receiver was tuned to 394.286 GHz in the lower sideband, and used the wobbler with an amplitude of 30″ and frequency of 1.5 Hz. We reduced the data using the standard procedures in the Continuum and Line Analysis Single-dish Software (CLASS; Pety 2005), and used an antenna gain of Jy K−1 = 43. We estimate the overall calibration uncertainty at 15%.

To observe the [OIII] 88 μm line (νrest = 3393.01 GHz), we used the Swedish-ESO PI Instrument for APEX (SEPIA; Belitsky et al. 2018). The observations were taken under the ESO project E-098.A-0513, on 2016 August 10 and 2017 April 26 with the first double sideband (DSB) incarnation of the receiver and the 4 GHz wide XFFTS backend, and under the Max Planck account M-0101.F-9522 on 2019 June 8 with the final sideband-separed (2SB) version of the receiver and the 8 GHz wide dFFTS4G backend. The total on-source integration time was 6.25 h, with PWV in the range 0.25–0.5 mm. The receiver was tuned to 649.403 GHz in the lower sideband, and used the wobbler with an amplitude of 40″ and frequency of 0.5 Hz in 2016/2017 and 30″ with 1.5 Hz in 2019. We reduced the data using the standard procedures in CLASS, using a first order fit for baseline removal. We used antenna gains of Jy K−1 = 120 for the DSB receiver and Jy K−1 = 70 for the 2SB receiver, based on carefully pointed and focused observations of Uranus close in time to the science target. We re-observed the 8.3σ DSB line with the 2SB receiver because after the DSB observations, a construction problem with the SEPIA band 9 selection mirror (NMF3 in Figs. 1 and 3 of Belitsky et al. 2018) was identified (and soon thereafter corrected). As a result, there were time-variable optical losses, which are difficult to characterize for long integrations like ours. This effect cannot be separated from other changes between 2016 and 2019, including a move of SEPIA to the central position in the Nasmyth A cabin using a new set of mirrors, the change of receiver from DSB to 2SB, the change of the telescope subreflector, and the replacement of the main telescope mirror surface panels to improve the overall telescope efficiency. However, thanks to the frequent planet calibration observations, we were able to combine these changes into the 70% improvement in the antenna gain, leading to consistent [OIII] 88 μm line detections of 108 ± 13 and 82 ± 16 Jy km s−1 for the DSB and 2SB receivers, respectively. For our analysis, we use the weigthed average flux of 108 ± 11 Jy km s−1.

2.2. Additional data

In addition to the two FSL reported above, we also include three other lines in our analysis. Gullberg et al. (2015) report a bright [CII] 158 μm line detected with APEX/FLASH, reproduced next to the [OIII] and [OI] lines in Fig. 1. Bothwell et al. (2017) report the detection of the [CI] 609 μm line from ALMA Cycle 3 spectroscopy. The [NII] 205 μm line was detected as part of an Atacama Compact Array (ACA) survey of 41 SPT sources (Cunningham et al. 2019). The [NII] 122 μm line was detected in a similar continuation of this ACA project (Cunningham et al., in prep.). Table 1 lists the main parameters from these observations.

thumbnail Fig. 1.

Spectra of [OIII] 88 μm, [NII] 122 μm, [OI] 145 μm, [CII] 158 μm (Gullberg et al. 2015), and [NII] 205 μm (Cunningham et al. 2019), fine structure lines in SPT 0418−47. The channel widths are 71, 60, 116, 70 and 50 km s−1, respectively. All velocities are given relative to z = 4.2248.

Table 1.

Observational parameters of the fine structure lines in SPT 0418−47 (no gravitational amplification corrections applied).

The dust continuum SED has been sampled by seven photometric points (Table 2), yielding a dust temperature Td = 45 ± 2 K (Strandet et al. 2016) and intrinsic far-IR (8–1000 μm) luminosity LFIR = (2.8 ± 0.3) × 1012 L⊙ (Bothwell et al. 2017). The lensing amplification factor μ = 32.7 ± 2.7 was derived by Spilker et al. (2016) based on 0 . ″ $ {{\overset{\prime\prime}{.}}} $5 resolution ALMA imaging. The lens model also shows that the lensed galaxy is very compact with an effective radius of 0 . ″ $ \overset{\prime \prime }{.} $092 ± 0.008 (∼0.6 kpc).

Table 2.

Continuum photometry of SPT 0418−47, taken from Strandet (2018) and Reuter et al., in prep.

3. Results

Figure 1 shows our detections of both the [OIII] 88 μm and [OI] 145 μm lines at the ∼5σ level. We fitted an unconstrained single Gaussian to the observed profiles and list the redshift, velocity width ΔV and integrated line luminosity in Table 1. For consistency, we have used the same fitting procedure for all lines shown in Fig. 1. We note that the [NII] 122 μm line width is less than half that of the other lines, but this may be due to the low signal to noise ratio (S/N) of our detection. Integrating over −300 to +350km s−1, that is the full velocity width of the bright [CII] 158 μm line yields an integrated [NII] 122 μm line flux of 7.1 ± 2.5 Jy km s−1. This is fully consistent with our unconstrained Gaussian fit. Deeper data are required to obtain a more reliable line flux, but the uncertainties are not to a level that they would change the conclusions in this paper.

The line luminosities are calculated following Solomon et al. (1997):

L line = 1.04 × 10 − 3 S line Δ V ν rest ( 1 + z ) − 1 D L 2 , $$ \begin{aligned} L_{\rm line}=1.04\times 10^{-3}S_{\rm line}\Delta V\nu _{\rm rest}(1+z)^{-1}D^{2}_{\mathrm{L}},\end{aligned} $$

where SlineΔV is the velocity-integrated line flux in Jy km s−1 and νrest is the rest frequency in GHz. The redshifts and velocity widths are consistent with each other within the uncertainties for all FSL, with the exception of the low S/N [NII] 122 μm line. The spectral profiles of the five well-detected FSL are also consistent with the CO(2-1) and CO(4-3) lines, suggesting that differential lensing is not important in this source (Gullberg et al. 2015; Aravena et al. 2016). Moreover there is no evidence for any companion galaxies in the lens model of SPT 0418−47 (Spilker et al. 2016), so we can safely assume that the integrated line fluxes are tracing a single galaxy.

We first compare our new FSL detections with other high-z galaxies. The only FSL that have been studied in significant samples of high-z galaxies are [CII] and [CI] (Stacey et al. 2010; Walter et al. 2011; Alaghband-Zadeh et al. 2013; Gullberg et al. 2015; Bothwell et al. 2017; Lagache et al. 2018). In both of these lines, the luminosity of SPT 0418−47 falls near the average of the SPT DSFGs. A first survey of the [NII] 205 μm lines has been observed in 41 galaxies from our SPT sample (Cunningham et al. 2019); SPT 0418−47 has been observed as part of this survey.

Detections of the [OIII] 88 μm line have now been reported in a dozen other high-z objects. Figure 2 compares the ratio of those [OIII] 88 μm lines with [CII] 158 μm as a function of LFIR. Compared to eleven other high-z DSFGs, SPT 0418−47 has one of the brightest relative [OIII] 88 μm emission, emitting about twice as much power as the [CII] 158 μm line. Only the z = 7.2 Lyman-α emitter SXDF-NB1006-2 has a ratio that is at least an order of magnitude higher, which has been argued as evidence for a very low metallicity (Inoue et al. 2016). Interestingly, SPT 0418−47 falls close to the starburst galaxy rather than the quasar in the double system CFHQS J2100 (Walter et al. 2018). The only low-z sample that reaches such similarly high LFIR is the Great Observatories All-sky LIRG Survey (GOALS; Sanders et al. 2003; Díaz-Santos et al. 2017). Compared to this sample, SPT 0418−47 has a high [OIII] 88 μm/[CII] 158 μm ratio, but there appears to be an increasing trend (and scatter) towards the highest LFIR. SPT 0418−47 appears to be more consistent with both the nearby dwarf galaxies (Cormier et al. 2015) and Seyferts (Spinoglio et al. 2015). Increasing the number statistics is clearly required to perform a better comparison, but our bright [OIII] 88 μm detection suggests that detecting this line may be easier than suggested by the previously reported upper limits.

thumbnail Fig. 2.

Ratio of the [OIII] 88 μm and [CII] 158 μm line luminosities as a function of (intrinsic) far-IR luminosity. The comparison low-z samples are the Herschel Dwarf Galaxies (green circles; Cormier et al. 2015), nearby Seyferts (orange diamonds; Spinoglio et al. 2015), and the GOALS sample (grey squares, assuming a 20% uncertainty in LFIR; Sanders et al. 2003; Díaz-Santos et al. 2017). The high-z sources marked in blue are SDP.81 (Valtchanov et al. 2011), the Eyelash (George et al. 2013), the Cloverleaf (Uzgil et al. 2016), SXDF-NB1006-2 (Inoue et al. 2016), SPT 0311−58 (Marrone et al. 2018), CFHQS J2100 (Walter et al. 2018), COSMOS-AzTEC1 (Tadaki et al. 2019), B14-65666 (Hashimoto et al. 2019), and J1342+0928 (Novak et al. 2019).

Thus far, only one other [OI] 145 μm detection has been reported at z >  1 (Yang et al. 2019). Although there are a few detections of the 63 μm line especially in stacked Herschel spectra (Uzgil et al. 2016; Wardlow et al. 2017; Wilson et al. 2017; Zhang et al. 2018), this is not the case in the fainter 145 μm line. To check if our detection is expected, we scale from the [OIII] 88 μm and [CII] 158 μm detections in the stacked spectrum of DSFGs (Zhang et al. 2018), which predicts a 63 μm luminosity of 1.6–5.5 × 1010 L⊙. To convert this to an [OI] 145 μm luminosity, we use the observed 145 μm to 63 μm [OI] line luminosity ratio between 2% and 35% (Brauher et al. 2008; Vasta et al. 2010), which yields a predicted range of 0.02–1.9 × 1010 L⊙. Our detection at the 2.1 × 1010 L⊙ level (Table 1) is at the upper bound of this range. However, we warn that the ratio of the two [OI] lines is subject to extinction and self-absorption effects (e.g. Vasta et al. 2010), and is therefore quite uncertain. We therefore also directly compare our observed [OI] 145 μm detection with low-z samples. Figure 3 plots the [CII] 158 μm and [OIII] 88 μm both normalized by the [OI] 145 μm luminosity. While the [OIII] 88 μm/[OI] 145 μm ratio is consistent with most normal galaxies (except dwarfs), the [CII] 158 μm in SPT 0418−47 appears to be underluminous relative to the [OI] 145 μm. In the next section, we combine all FSL results to determine the physical conditions of the ISM in SPT 0418−47.

thumbnail Fig. 3.

Line luminosity ratio of [CII] 158 μm and [OI] 145 μm against [OIII] 88 μm and [OI] 145 μm. The low-z samples are drawn from Fernández-Ontiveros et al. (2016), Brauher et al. (2008), Herrera-Camus et al. (2018), Cormier et al. (2015) and Spinoglio et al. (2015). We note that relative to [OI] 145 μm, SPT 0418−47 (red star) has an [OIII] 88 μm luminosity that is rather consistent with normal (black, brown and blue) rather than dwarf galaxies (green), but a weaker [CII] 158 μm. Arrows indicate which physical processes dominate these line ratios (see text).

4. Discussion

4.1. The nature of SPT 0418−47

Before discussing the physical parameters of its ISM, we first compare SPT 0418−47 with the main sequence galaxy population at z ∼ 4. Unfortunately, a direct stellar mass determination is rather complex because the foreground lensing galaxy outshines the background DSFG (Vieira et al. 2013). Using the CIGALE SED modelling described in §4.4 and the deblending technique of Ma et al. (2015), we derive an intrinsic M* = (9.5 ± 3.0) × 109 M⊙ (Phadke et al., in prep.). This value can be compared with the stellar masses 2–8 × 1010 M⊙ in four other SPT DSFGs derived using HST and Spitzer imaging (Ma et al. 2015). Using the gas masses derived from low-J CO lines (Aravena et al. 2016) in these four DSFGs, this implies stellar to gas mass ratios in the range 0.3–1.0. The relatively low gas mass in SPT 0418−47 of 0.5 × 1010 M⊙ then suggests the stellar mass in the range 0.5–1.7 × 1010 M⊙, fully consistent with the SED-derived M* above. Detailed JWST imaging of SPT0418-47 scheduled as part of the Early Release Science observations will provide a much more reliable stellar mass determination.

This rather low stellar mass contrasts with the high SFR = 280 ± 30 M⊙ yr−1 derived as 10−10 × LFIR. This places SPT 0418 −47 well above the galaxy main sequence at z ∼ 4, which predicts a SFR ∼ 25 M⊙yr−1 for a 1010M⊙ galaxy (Schreiber et al. 2017).

4.2. Contributions of PDR and HII regions

As mentioned in the introduction, the FSL we observed originate from both the PDR and HII regions. It is important to keep in mind that we are only looking at galaxy-integrated line luminosities rather than ratios within a single Giant Molecular Cloud (GMC). What we are comparing are thus the luminosity-weighted combinations of all GMCs within these galaxies.

The line that has the most uncertain relative contribution of PDR and HII regions is [CII] 158 μm. To estimate the expected contribution of HII regions to the [CII] 158 μm line flux, we can use the [NII] 205 μm line, which has the same critical density and similar ionization range. Following Díaz-Santos et al. (2017), we assume L [ CII ] ion / L [ NII ] ≃ 3.0 ± 0.5 $ L_{{[\mathrm{CII}]}}^{\mathrm{ion}}/L_{[\mathrm{NII}]} \simeq 3.0 \pm 0.5 $, which covers the usual full range of HII region physical conditions. Comparing this to the observed ratio of 24.5 ± 2.0 (Table 1), we derive an HII fraction of 12 ± 2%. While there is some remaining model uncertainty on this number, it seems clear that the [CII] line is dominated by PDRs. Given that this contribution is rather small, we do not correct for it in the following; moreover such a correction cannot be calculated for all comparison samples due to the absence of systematic [NII] 205 μm observations.

Another line ratio that can be used to determine the ratio of the filling factors of warm ionized and dense neutral gas is [OIII] 88 μm/[OI] 63 μm (Cormier et al. 2015, 2019; Díaz-Santos et al. 2017). This ratio provides a measure that is completely independent of abundance, and only very weakly on optical depth. The ratio could depend on temperature, though it would in this case measure the variation in the ratio of temperatures in the PDR and HII phases between galaxies. Díaz-Santos et al. (2017, their Fig. 7) do not find such a trend in the local luminious infrared galaxies from GOALS. The only galaxies that show clearly higher [OIII] 88 μm/[OI] 63 μm ratios are dwarf galaxies, which has been interpreted as indicating the presence of low density channels allowing the hard UV photons to reach the outer parts of galaxies (Cormier et al. 2015, 2019). As a result, these low metallicity galaxies are likely to be filled with relatively low density (below n crit e − ( [ OIII ] 88 μ $ n^{e^-}_{\mathrm{crit}}(\rm [OIII]\,88\,\mu $m) ≃ 500 cm−3) gas rather than PDRs. One uncertainty in this ratio remains the optical depth and self-absorption of the [OI] 63 μm line. We here use instead the [OI] 145 μm line, which is not subject to these uncertainties. Figure 3 confirms that the dwarf galaxies (green points) are clearly offset from the other samples, with the normal galaxies and Seyferts of Brauher et al. (2008) (orange points) lying in between. SPT 0418−47 instead lies at the lower end of the [OIII] 88 μm/[OI] 145 μm distribution, which implies that the galaxy-integrated line emission is dominated by PDR. This confirms our conclusion above from the [CII] 158 μm/[NII] 205 μm ratio.

4.3. Gas density

To determine the gas density, one should ideally use the ratio of two FSL of the same ion or neutral atom such as [OIII] 88 μm/52 μm, [NII] 122 μm/205 μm, or [OI] 145 μm/63 μm (e.g. Spinoglio et al. 2015; Herrera-Camus et al. 2016). Of these, we only have both [NII] lines, which allow to probe the density range 1 cm−3 <  ne <  103 cm−3. Our observed [NII] 122 μm/205 μm luminosity ratio of 1.40 ± 0.43 corresponds to an electron density ne = 50 ± 20 cm−3 following the relation shown in for example Fig. 2 of Herrera-Camus et al. (2016). This places SPT 0418−47 slightly above the average density found in nearby spirals and ultra-luminous galaxies from the “Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel” (KINGFISH), “Beyond the Peak” (BtP), and GOALS samples (Herrera-Camus et al. 2016; Díaz-Santos et al. 2017). We are currently performing a survey of this [NII] 122 μm/205 μm ratio in the SPT sample, which will be the first systematic survey in high-redshift galaxies, and refer to a forthcoming publication for a more in depth discussion (Cunningham et al., in prep.).

It is important to keep in mind that the [NII] lines only trace the diffuse HII regions, but not the denser PDR. Bothwell et al. (2017) modelled the observed ratio of four PDR-dominated lines (CO(2-1), CO(4-3), [CI] 609 μm, and [CII] 158 μm) with the 3D-PDR code (Bisbas et al. 2012), which also includes the effect of cosmic rays. They derive a consistent model with a PDR density of n ∼ 2.5 × 105 cm−3, and a far-UV field strength G0 ∼ 3.2 × 105. This places SPT 0418−47 at the high end of the observed PDR gas densities, while the G0 is the highest one in sample of 13 SPT DSFGs modelled by Bothwell et al. (2017), that is firmly in the region occupied by ULIRGs. As mentioned by Bothwell et al. (2017), the densities from the 3D-PDR model are almost an order of magnitude higher than the ones derived using the PDR toolbox2 (Kaufman et al. 2006; Pound & Wolfire 2008). The [OI] 145 μm line adds a powerful new diagnostic to the PDR modelling, as the [OI] 145 μm to [CII] 158 μm ratio is particularly sensitive to density, especially in the high G0 range. This ratio has the advantage of being mostly immune to extinction effects, but one has to keep in mind that it is also subject to abundance variations and HII contributions to [CII] 158 μm. From the PDR toolbox modelling of all five PDR-dominated lines, we derive a G0 ∼ 103 and n ∼ 2 × 104 cm−3. We defer a more detailed PDR modelling using 3D-PDR to a future publication.

4.4. AGN contributions

Optical line ratios are often used as an AGN diagnostic (the classical “BPT diagram”; Baldwin et al. 1981). In lensed DSFGs, this technique can unfortunately not be used as the optical spectra are dominated by the foreground lensing galaxy, while the background DSFG is almost completely dust obscured. At far-IR wavelengths, the situation is reversed: our sub-arcsecond 850 μm ALMA imaging (Vieira et al. 2013) shows that the lens is undetected while the background DSFG is strongly detected. We can thus use the hot dust continuum shape and/or the mid- and far-IR FSL to look for AGN contributions.

Turning first to the dust continuum SED (Fig. 4), we find that it is consistent with the starburst template of Magdis et al. (2012), suggesting that an AGN is not required. To obtain a quantitative measure of the AGN fraction, we use the method described by Spilker et al. (2018), using the CIGALE SED modelling code (Noll et al. 2009; Boquien et al. 2019) including the dust models of Fritz et al. (2006). We derive an AGN fraction < 5%. This is consitent with results in other DSFGs from the same SPT sample, where no significant AGN fractions were found (Ma et al. 2016; Spilker et al. 2018; Apostolovski et al. 2019).

thumbnail Fig. 4.

Continuum Spectral Energy Distribution of SPT 0418−47 adjusted with the template of Magdis et al. (2012) and a simple quadratic interpolation.

The alternative method we examine is to use the FSL. Several authors have used mid-IR and far-IR FSL to determine AGN contributions (e.g. Fernández-Ontiveros et al. 2016; Herrera-Camus et al. 2018). These tracers are mostly based on the ratio of very high and lower ionization lines of the same species. Unfortunately, most of these species are in the mid-IR (e.g. [OIV] 25.9 μm, [NeIII] 15.6 μm, [NeII] 12.8 μm, [SIV] 10.5 μm, [SIII] 18.7 μm), and are inaccessible from the ground, even at high redshift. However, one ratio that does appear to separate AGN regions is the [OI] 63 μm/[CII] 158 μm line (Fig. 12 of Herrera-Camus et al. 2018). While we did not observe the [OI] 63 μm line in SPT 0418−47, in Fig. 3, we compare the [CII] 158 μm/[OI] 145 μm ratio in a variety of galaxies. The advantage of the 145 μm line is that it is not affected by extinction and self-absorption. Cormier et al. (2019) also showed that adding X-ray dominated regions to their CLOUDY models leads to a slight increase by a factor 1.2 in the [OI] 145 μm/63 μm luminosity ratio (their Fig. 8). A similar effect is seen when increasing the cosmic ray flux. As AGN would produce both stronger X-rays and cosmic ray fluxes, the 145 μm line is therefore at least as good an AGN tracer, if not a better one. In addition, AGN are expected to exite the higher [OI] level, that is increase the [OI] 145 μm/63 μm ratio.

Figure 3 shows that the local AGN (blue symbols) and Seyferts (brown and a number of black symbols) cover the full range of [CII] 158 μm/[OI] 145 μm, but they are also the only ones with [CII] 158 μm/[OI] 145 μm ≲ 5. Interestingly, SPT 0418−47 falls among these lowest ratios, suggesting the possible presence of an AGN. However, we warn that this low [CII] 158 μm/[OI] 145 μm ratio may also be a direct consequence of the relatively high Tdust = 45 ± 2 K. Such higher dust and gas temperatures may be due to a higher Far-UV flux G0, which will also increase [OI] faster than [CII] (assuming the 145 μm line behaves like the 63 μm line, Malhotra et al. 2001). Surprisingly, the AGN galaxies do not show higher [OIII] 88 μm/[OI] 145 μm ratios. This may be because we are comparing only galaxy-integrated line luminosities, and AGN emission is known to be very anisotropic due to an obscuring torus; the ionized gas is thus only observed in ionization cones and not throughout the galaxy. As such, the fact that SPT 0418−47 has a ratio within the region occupied by both AGN and normal galaxies does not provide any information on the presence of an AGN. We conclude that while the SED seems to suggest a negligible AGN contribution in SPT 0418−47, the FSL ratio method is more inconclusive.

4.5. Metallicity

As argued by Nagao et al. (2011) and Pereira-Santaella et al. (2017), the ratio of [OIII] 52 μm + [OIII] 88 μm and [NIII] 57 μm is the most promising tracer of the gas metallicity using FSL. While these lines have now been observed with Herschel in a significant number of nearby objects, they are beyond the reach of ALMA for all but the very highest redshift objects known to date. More accessible alternatives are the ratios of [OIII] 88 μm with the [NII] 122 μm and/or [NII] 205 μm lines (Pereira-Santaella et al. 2017; Tadaki et al. 2019). The main disadvantage of using the [NII] lines is that they have a much lower ionization potential of 14.53 eV compared to 35.12 eV for [OIII] 88 μm, creating a degeneracy between the ionization parameter U and metallicity Z. An additional, though slightly less critical degeneracy exists with the gas electron density ne. To solve this degeneracy between U and Z, Rigopoulou et al. (2018) use the ratio of the restframe 88/122 μm continuum flux measurements. This assumes that the gas and dust are well mixed; the dust temperature Td will then increase with the ionizing flux, creating a correlation between the continuum flux ratio and U (see also Díaz-Santos et al. 2017). For PDR-dominated lines, a similar correlation between the Far-UV flux G0 and the 60/100 μm continuum flux ratio has indeed been observed in ISO data (e.g. Malhotra et al. 2001; Brauher et al. 2008), suggesting the gas and dust temperatures are indeed correlated.

To allow a direct comparison with the method of Rigopoulou et al. (2018), we first derive the rest frame 88 and 122 μm continuum flux ratio C(88−122) by fitting the continuum photometry of SPT 0418−47 (Table 2) with the template of Magdis et al. (2012) and a simple quadratic interpolation. Figure 4 shows that the results do not depend strongly on the adopted function, with C(88−122) = 1.26 ± 0.01 for the template and C(88−122) = 1.13 ± 0.14 for the quadratic interpolation. We conservatively adopt the latter, which yields −3.2 <  logU <  −2.0 using Fig. 4 of Rigopoulou et al. (2018).

In Fig. 5, we reproduce the CLOUDY models of Pereira-Santaella et al. (2017), where the coloured bands indicate different ionization parameters, and their spread different gas electron densities. The logU constraint limits the range of models to the orange, yellow and green shaded areas in Fig. 5. We note that there is very little dependence on the density, as the [OIII] 88 μm and [NII] 122 μm lines have very similar n crit e − $ n_{{\rm{crit}}}^{{e^ - }} $ = 510 and 310 cm−3, respectively. The observed [OIII] 88 μm to [NII] 122 μm ratio thus suggests 0.3 <  Z/Z⊙ <  1.3. This is consistent with the [CII] 158 μm/[NII] 205 μm ratio, which has also been used as a rough metallicity indicator (Nagao et al. 2012). The ratio in SPT 0418−47 is lower than the one in ALESS 73.1 (De Breuck et al. 2014), which also suggests a near solar metallicity.

thumbnail Fig. 5.

[OIII] 88 μm over [NII] 122 μm line ratio as a function of metallicity. The colour shaded areas mark the range of ratios predicted by the CLOUDY models of Pereira-Santaella et al. (2017) with the same ionization parameter but different densities (circles n = 10 cm−3 and squares n = 104 cm−3). The colour coding (grey, red, orange, yellow, and green) indicates the ionization parameter (log U = −4.0, −3.5, −3.0, −2.5, and −2.0, respectively). The solid horizontal purple line marks the observed ratio in SPT 0418−47 with the uncertainty range delimited by dashed lines. The thick part of the horizontal line indicates the allowed range of log U as derived from the continuum SED (Fig. 4).

Our new metallicity determination allows us to compare SPT 0418−47 to the mass-metallicity relationship studied at lower redshifts, and verify if this relation evolves over cosmic time. Although there is a rather large uncertainty, we derived a rather low stellar mass of ∼1010 M⊙ (see Sect. 4.1). The almost solar metallicity in SPT 0418−47 (corresponding to 8.18 <  12 + log(O/H) < 8.81) places it a bit higher on the mass-metallicity relation compared to other DSFGs at slightly lower redshifts (Fig. 7 of Rigopoulou et al. 2018), near the track for z = 2.2 from Maiolino et al. (2008) and Mannucci et al. (2010).

5. Conclusions

We have detected five FSL in SPT 0418−47, including the second detection of [OI] 145 μm at z >  1. Most of the analytic power of the FSL comes from line ratios. We draw the following conclusions from our study: Firstly, the ratio of [OIII] 88 μm to [OI] 145 μm suggests that the emission from the ISM in SPT 0418−47 is dominated by PDR rather than HII regions. The [CII] 158 μm line is almost 90% dominated by PDR. Secondly, the ratios of [CII] 158 μm to [OI] 145 μm and of [NII] 205 μm to [NII] 122 μm suggest that both the PDR and HII regions have high gas densities. Thirdly, using a new diagram of [CII] 158 μm to [OI] 145 μm against [OIII] 88 μm to [OI] 145 μm, we find that SPT 0418−47 lies in a region occupied by nearby AGN in low redshift samples. While this suggests that SPT 0418−47 may also harbour an AGN, we warn that the observed line ratios may also be a direct consequence of the relatively warm Td = 45.3 ± 2.3 K created by an extremely powerful starburst. Fourthly, following the method of Rigopoulou et al. (2018), we derive an ionization parameter −3.2 <  logU <  −2.0 from the dust continuum flux ratios. Combined with the observed [OIII] 88 μm to [NII] 122 μm line ratio, we derive a gas metallicity of 0.3 <  Z/Z⊙ <  1.3 in SPT 0418−47.

Taken together, these results imply that SPT 0418−47 has a dense, roughly solar metallicty ISM that is dominated by PDR. The presence of an AGN cannot be excluded, but it is equally possible that the extreme strength of the starburst is heating both the dust and gas in the entire galaxy to higher than average values.

The high gravitational amplification (μ = 32.7 ± 2.7) of SPT 0418−47 has allowed us to detect even the weaker [OI] 145 μm line with APEX. This line has been often neglected at high redshift, probably because it was also difficult to study in nearby samples before the advent of Herschel. It is a powerful alternative to its brighter counterpart at 63 μm, being at a significantly easier frequency to observe from the ground at z >  1.17, rather than z >  4 for the 63 μm line. Moreover, the 145 μm line is not subject to extinction and self-absorption affects, which may reduce the observed flux of the 63 μm line. We expect the [OI] 145 μm line to become a powerful new tracer of PDR in high-z galaxies.

1

This assumes the [CII] 158 μm is optically thin; in the SPT sample, the line appears to have moderate optical depth (Gullberg et al. 2015). For optically thick [CII] 158 μm, the critical densities are lower by an order of magnitude, see Fig. 3 of Lagache et al. (2018).

Acknowledgments

We thank the anonymous referee for useful suggestions, Myha Vuong for stimulating discussions on fine structure lines, and the APEX staff for their excellent support in obtaining the data. This publication is based on data acquired with the Atacama Pathfinder Experiment (APEX) under programme IDs 098.A-0513 (ESO), 0101.F-9522, and 093.F-9512 (both Max Planck). APEX is a collaboration between the Max-Planck-Institut fur Radioastronomie, the European Southern Observatory, and the Onsala Space Observatory. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2016.1.00133.T and ADS/JAO.ALMA#2018.1.01670.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 SPT is supported by the NSF through grant PLR-1248097, with partial support through PHY-1125897, the Kavli Foundation and the Gordon and Betty Moore Foundation grant GBMF 947. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc.

References

  1. Abel, N. P., Ferland, G. J., Shaw, G., & van Hoof, P. A. M. 2005, ApJS, 161, 65 [NASA ADS] [CrossRef] [Google Scholar]
  2. Alaghband-Zadeh, S., Chapman, S. C., Swinbank, A. M., et al. 2013, MNRAS, 435, 1493 [Google Scholar]
  3. Apostolovski, Y., Aravena, M., Anguita, T., et al. 2019, A&A, 628, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Aravena, M., Spilker, J. S., Bethermin, M., et al. 2016, MNRAS, 457, 4406 [NASA ADS] [CrossRef] [Google Scholar]
  5. Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Belitsky, V., Lapkin, I., Fredrixon, M., et al. 2018, A&A, 612, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Béthermin, M., De Breuck, C., Gullberg, B., et al. 2016, A&A, 586, L7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Bisbas, T. G., Bell, T. A., Viti, S., et al. 2012, MNRAS, 427, 2100 [NASA ADS] [CrossRef] [Google Scholar]
  9. Boquien, M., Burgarella, D., Roehlly, Y., et al. 2019, A&A, 622, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Bothwell, M. S., Aguirre, J. E., Aravena, M., et al. 2017, MNRAS, 466, 2825 [Google Scholar]
  11. Brauher, J. R., Dale, D. A., & Helou, G. 2008, ApJS, 178, 280 [NASA ADS] [CrossRef] [Google Scholar]
  12. Carniani, S., Maiolino, R., Pallottini, A., et al. 2017, A&A, 605, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Cigan, P., Young, L., Cormier, D., et al. 2016, AJ, 151, 14 [NASA ADS] [CrossRef] [Google Scholar]
  14. Combes, F., Rex, M., Rawle, T. D., et al. 2012, A&A, 538, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Coppin, K. E. K., Danielson, A. L. R., Geach, J. E., et al. 2012, MNRAS, 427, 520 [NASA ADS] [CrossRef] [Google Scholar]
  16. Cormier, D., Madden, S. C., Lebouteiller, V., et al. 2015, A&A, 578, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Cormier, D., Abel, N. P., Hony, S., et al. 2019, A&A, 626, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Croxall, K. V., Smith, J. D., Pellegrini, E., et al. 2017, ApJ, 845, 96 [NASA ADS] [CrossRef] [Google Scholar]
  19. Cunningham, D. J. M., Chapman, S. C., Aravena, M., et al. 2019, MNRAS, submitted [arXiv:1906.02293] [Google Scholar]
  20. De Breuck, C., Röttgering, H., Miley, G., van Breugel, W., & Best, P. 2000, A&A, 362, 519 [NASA ADS] [Google Scholar]
  21. De Breuck, C., Williams, R. J., Swinbank, M., et al. 2014, A&A, 565, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Decarli, R., Walter, F., Neri, R., et al. 2012, ApJ, 752, 2 [NASA ADS] [CrossRef] [Google Scholar]
  23. Decarli, R., Walter, F., Carilli, C., et al. 2014, ApJ, 782, L17 [NASA ADS] [CrossRef] [Google Scholar]
  24. Díaz-Santos, T., Armus, L., Charmandaris, V., et al. 2017, ApJ, 846, 32 [NASA ADS] [CrossRef] [Google Scholar]
  25. Farrah, D., Lebouteiller, V., Spoon, H. W. W., et al. 2013, ApJ, 776, 38 [NASA ADS] [CrossRef] [Google Scholar]
  26. Ferkinhoff, C., Hailey-Dunsheath, S., Nikola, T., et al. 2010, ApJ, 714, L147 [NASA ADS] [CrossRef] [Google Scholar]
  27. Ferkinhoff, C., Brisbin, D., Nikola, T., et al. 2011, ApJ, 740, L29 [NASA ADS] [CrossRef] [Google Scholar]
  28. Ferkinhoff, C., Brisbin, D., Nikola, T., et al. 2015, ApJ, 806, 260 [NASA ADS] [CrossRef] [Google Scholar]
  29. Fernández-Ontiveros, J. A., Spinoglio, L., Pereira-Santaella, M., et al. 2016, ApJS, 226, 19 [NASA ADS] [CrossRef] [Google Scholar]
  30. Fritz, J., Franceschini, A., & Hatziminaoglou, E. 2006, MNRAS, 366, 767 [Google Scholar]
  31. George, R. D., Ivison, R. J., Hopwood, R., et al. 2013, MNRAS, 436, L99 [NASA ADS] [CrossRef] [Google Scholar]
  32. Graciá-Carpio, J., Sturm, E., Hailey-Dunsheath, S., et al. 2011, ApJ, 728, L7 [NASA ADS] [CrossRef] [Google Scholar]
  33. Gullberg, B., De Breuck, C., Vieira, J. D., et al. 2015, MNRAS, 449, 2883 [NASA ADS] [CrossRef] [Google Scholar]
  34. Güsten, R., Nyman, L. Å., Schilke, P., et al. 2006, A&A, 454, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Hashimoto, T., Laporte, N., Mawatari, K., et al. 2018, Nature, 557, 392 [NASA ADS] [CrossRef] [Google Scholar]
  36. Hashimoto, T., Inoue, A. K., Mawatari, K., et al. 2019, PASJ, 70, S8 [Google Scholar]
  37. Herrera-Camus, R., Bolatto, A., Smith, J. D., et al. 2016, ApJ, 826, 175 [Google Scholar]
  38. Herrera-Camus, R., Sturm, E., Graciá-Carpio, J., et al. 2018, ApJ, 861, 94 [NASA ADS] [CrossRef] [Google Scholar]
  39. Inoue, A. K., Tamura, Y., Matsuo, H., et al. 2016, Science, 352, 1559 [NASA ADS] [CrossRef] [Google Scholar]
  40. Kamenetzky, J., Rangwala, N., Glenn, J., Maloney, P. R., & Conley, A. 2014, ApJ, 795, 174 [NASA ADS] [CrossRef] [Google Scholar]
  41. Kaufman, M. J., Wolfire, M. G., & Hollenbach, D. J. 2006, ApJ, 644, 283 [NASA ADS] [CrossRef] [Google Scholar]
  42. Kennicutt, R. C., & Evans, N. J. 2012, ARA&A, 50, 531 [NASA ADS] [CrossRef] [Google Scholar]
  43. Klein, T., Ciechanowicz, M., Leinz, C., et al. 2014, IEEE Trans. Terahertz Sci. Technol., 4, 588 [NASA ADS] [CrossRef] [Google Scholar]
  44. Lagache, G., Cousin, M., & Chatzikos, M. 2018, A&A, 609, A130 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Lu, N., Zhao, Y., Díaz-Santos, T., et al. 2017, ApJ, 842, L16 [NASA ADS] [CrossRef] [Google Scholar]
  46. Ma, J., Gonzalez, A. H., Spilker, J. S., et al. 2015, ApJ, 812, 88 [NASA ADS] [CrossRef] [Google Scholar]
  47. Ma, J., Gonzalez, A. H., Vieira, J. D., et al. 2016, ApJ, 832, 114 [CrossRef] [Google Scholar]
  48. Magdis, G. E., Daddi, E., Béthermin, M., et al. 2012, ApJ, 760, 6 [NASA ADS] [CrossRef] [Google Scholar]
  49. Maiolino, R., Nagao, T., Grazian, A., et al. 2008, A&A, 488, 463 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Malhotra, S., Kaufman, M. J., Hollenbach, D., et al. 2001, ApJ, 561, 766 [Google Scholar]
  51. Mannucci, F., Cresci, G., Maiolino, R., et al. 2010, MNRAS, 408, 2115 [NASA ADS] [CrossRef] [Google Scholar]
  52. Marrone, D. P., Spilker, J. S., Hayward, C. C., et al. 2018, Nature, 553, 51 [NASA ADS] [CrossRef] [Google Scholar]
  53. Moriwaki, K., Yoshida, N., Shimizu, I., et al. 2018, MNRAS, 481, L84 [NASA ADS] [CrossRef] [Google Scholar]
  54. Nagao, T., Maiolino, R., Marconi, A., & Matsuhara, H. 2011, A&A, 526, A149 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Nagao, T., Maiolino, R., De Breuck, C., et al. 2012, A&A, 542, L34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Noll, S., Burgarella, D., Giovannoli, E., et al. 2009, A&A, 507, 1793 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Novak, M., Bañados, E., Decarli, R., et al. 2019, ApJ, 881, 63 [NASA ADS] [CrossRef] [Google Scholar]
  58. Palay, E., Nahar, S. N., Pradhan, A. K., et al. 2012, MNRAS, 423, L35 [NASA ADS] [CrossRef] [Google Scholar]
  59. Pavesi, R., Riechers, D. A., Capak, P. L., et al. 2016, ApJ, 832, 151 [NASA ADS] [CrossRef] [Google Scholar]
  60. Pereira-Santaella, M., Rigopoulou, D., Farrah, D., Lebouteiller, V., & Li, J. 2017, MNRAS, 470, 1218 [NASA ADS] [CrossRef] [Google Scholar]
  61. Pety, J. 2005, SF2A-2005: Semaine de l’Astrophysique Francaise, 721 [Google Scholar]
  62. Planck Collaboration XIII. 2016, A&A, 594, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Pound, M. W., & Wolfire, M. G. 2008, in Astronomical Data Analysis Software and Systems XVII, eds. R. W. Argyle, P. S. Bunclark, & J. R. Lewis, ASP Conf. Ser., 394, 654 [NASA ADS] [Google Scholar]
  64. Rangwala, N., Maloney, P. R., Glenn, J., et al. 2011, ApJ, 743, 94 [NASA ADS] [CrossRef] [Google Scholar]
  65. Rigopoulou, D., Pereira-Santaella, M., Magdis, G. E., et al. 2018, MNRAS, 473, 20 [NASA ADS] [CrossRef] [Google Scholar]
  66. Rosenberg, M. J. F., van der Werf, P. P., Aalto, S., et al. 2015, ApJ, 801, 72 [NASA ADS] [CrossRef] [Google Scholar]
  67. Sanders, D. B., Mazzarella, J. M., Kim, D.-C., et al. 2003, AJ, 126, 1607 [Google Scholar]
  68. Sargsyan, L., Samsonyan, A., Lebouteiller, V., et al. 2014, ApJ, 790, 15 [NASA ADS] [CrossRef] [Google Scholar]
  69. Schreiber, C., Pannella, M., Leiton, R., et al. 2017, A&A, 599, A134 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Solomon, P. M., Downes, D., Radford, S. J. E., et al. 1997, ApJ, 478, 144 [NASA ADS] [CrossRef] [Google Scholar]
  71. Spilker, J. S., Marrone, D. P., Aravena, M., et al. 2016, ApJ, 826, 112 [NASA ADS] [CrossRef] [Google Scholar]
  72. Spilker, J. S., Aravena, M., Béthermin, M., et al. 2018, Science, 361, 1016 [NASA ADS] [CrossRef] [Google Scholar]
  73. Spinoglio, L., Pereira-Santaella, M., Dasyra, K. M., et al. 2015, ApJ, 799, 21 [NASA ADS] [CrossRef] [Google Scholar]
  74. Stacey, G. J., Hailey-Dunsheath, S., Ferkinhoff, C., et al. 2010, ApJ, 724, 957 [NASA ADS] [CrossRef] [Google Scholar]
  75. Strandet, M. 2018, PhD Thesis, Rheinischen Friedrich-Wilhelms-Universität Bonn [Google Scholar]
  76. Strandet, M., Weiss, A., Vieira, J., et al. 2016, ApJ, 822, 80 [NASA ADS] [CrossRef] [Google Scholar]
  77. Tadaki, K.-I., Iono, D., Hatsukade, B., et al. 2019, ApJ, 876, 1 [NASA ADS] [CrossRef] [Google Scholar]
  78. Tamura, Y., Mawatari, K., Hashimoto, T., et al. 2019, ApJ, 874, 27 [NASA ADS] [CrossRef] [Google Scholar]
  79. Uzgil, B. D., Bradford, C. M., Hailey-Dunsheath, S., Maloney, P. R., & Aguirre, J. E. 2016, ApJ, 832, 209 [NASA ADS] [CrossRef] [Google Scholar]
  80. Valtchanov, I., Virdee, J., Ivison, R. J., et al. 2011, MNRAS, 415, 3473 [NASA ADS] [CrossRef] [Google Scholar]
  81. Vasta, M., Barlow, M. J., Viti, S., et al. 2010, MNRAS, 404, 1910 [NASA ADS] [Google Scholar]
  82. Vieira, J. D., Marrone, D. P., Chapman, S. C., et al. 2013, Nature, 495, 344 [NASA ADS] [CrossRef] [Google Scholar]
  83. Villar-Martin, M., Tadhunter, C., & Clark, N. 1997, A&A, 323, 21 [NASA ADS] [Google Scholar]
  84. Vishwas, A., Ferkinhoff, C., Nikola, T., et al. 2018, ApJ, 856, 174 [NASA ADS] [CrossRef] [Google Scholar]
  85. Walter, F., Weiß, A., Downes, D., Decarli, R., & Henkel, C. 2011, ApJ, 730, 18 [NASA ADS] [CrossRef] [Google Scholar]
  86. Walter, F., Riechers, D., Novak, M., et al. 2018, ApJ, 869, L22 [NASA ADS] [CrossRef] [Google Scholar]
  87. Wardlow, J. L., Cooray, A., Osage, W., et al. 2017, ApJ, 837, 12 [NASA ADS] [CrossRef] [Google Scholar]
  88. Weiß, A., De Breuck, C., Marrone, D. P., et al. 2013, ApJ, 767, 88 [NASA ADS] [CrossRef] [Google Scholar]
  89. Wilson, D., Cooray, A., Nayyeri, H., et al. 2017, ApJ, 848, 30 [NASA ADS] [CrossRef] [Google Scholar]
  90. Yang, J., Venemans, B., Wang, F., et al. 2019, ApJ, 880, 153 [Google Scholar]
  91. Zhang, Z.-Y., Ivison, R. J., George, R. D., et al. 2018, MNRAS, 481, 59 [NASA ADS] [CrossRef] [Google Scholar]
  92. Zhao, Y., Lu, N., Xu, C. K., et al. 2016, ApJ, 819, 69 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1.

Observational parameters of the fine structure lines in SPT 0418−47 (no gravitational amplification corrections applied).

In the text
Table 2.

Continuum photometry of SPT 0418−47, taken from Strandet (2018) and Reuter et al., in prep.

In the text

All Figures

thumbnail Fig. 1.

Spectra of [OIII] 88 μm, [NII] 122 μm, [OI] 145 μm, [CII] 158 μm (Gullberg et al. 2015), and [NII] 205 μm (Cunningham et al. 2019), fine structure lines in SPT 0418−47. The channel widths are 71, 60, 116, 70 and 50 km s−1, respectively. All velocities are given relative to z = 4.2248.
In the text
thumbnail Fig. 2.

Ratio of the [OIII] 88 μm and [CII] 158 μm line luminosities as a function of (intrinsic) far-IR luminosity. The comparison low-z samples are the Herschel Dwarf Galaxies (green circles; Cormier et al. 2015), nearby Seyferts (orange diamonds; Spinoglio et al. 2015), and the GOALS sample (grey squares, assuming a 20% uncertainty in LFIR; Sanders et al. 2003; Díaz-Santos et al. 2017). The high-z sources marked in blue are SDP.81 (Valtchanov et al. 2011), the Eyelash (George et al. 2013), the Cloverleaf (Uzgil et al. 2016), SXDF-NB1006-2 (Inoue et al. 2016), SPT 0311−58 (Marrone et al. 2018), CFHQS J2100 (Walter et al. 2018), COSMOS-AzTEC1 (Tadaki et al. 2019), B14-65666 (Hashimoto et al. 2019), and J1342+0928 (Novak et al. 2019).
In the text
thumbnail Fig. 3.

Line luminosity ratio of [CII] 158 μm and [OI] 145 μm against [OIII] 88 μm and [OI] 145 μm. The low-z samples are drawn from Fernández-Ontiveros et al. (2016), Brauher et al. (2008), Herrera-Camus et al. (2018), Cormier et al. (2015) and Spinoglio et al. (2015). We note that relative to [OI] 145 μm, SPT 0418−47 (red star) has an [OIII] 88 μm luminosity that is rather consistent with normal (black, brown and blue) rather than dwarf galaxies (green), but a weaker [CII] 158 μm. Arrows indicate which physical processes dominate these line ratios (see text).
In the text
thumbnail Fig. 4.

Continuum Spectral Energy Distribution of SPT 0418−47 adjusted with the template of Magdis et al. (2012) and a simple quadratic interpolation.
In the text
thumbnail Fig. 5.

[OIII] 88 μm over [NII] 122 μm line ratio as a function of metallicity. The colour shaded areas mark the range of ratios predicted by the CLOUDY models of Pereira-Santaella et al. (2017) with the same ionization parameter but different densities (circles n = 10 cm−3 and squares n = 104 cm−3). The colour coding (grey, red, orange, yellow, and green) indicates the ionization parameter (log U = −4.0, −3.5, −3.0, −2.5, and −2.0, respectively). The solid horizontal purple line marks the observed ratio in SPT 0418−47 with the uncertainty range delimited by dashed lines. The thick part of the horizontal line indicates the allowed range of log U as derived from the continuum SED (Fig. 4).
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.