{'fraction_non_alphanumeric': 0.05523471896620733, 'fraction_numerical': 0.05008862713705758, 'mean_word_length': 4.16996748448123, 'pattern_counts': {'":': 0, '<': 5, '<?xml version=': 0, '>': 6, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
ALMA reveals a chemically evolved submillimeter galaxy at z=4.76

The chemical properties of high-z galaxies provide important information to constrain galaxy evolutionary scenarios. However, widely-used metallicity diagnostics based on rest-frame optical emission lines are not usable for heavily dust-enshrouded galaxies (such as Sub-Millimeter Galaxies; SMGs), especially at z>3. Here we focus on the flux ratio of the far-infrared fine-structure emission lines [NII]205um and [CII]158um to assess the metallicity of high-z SMGs. Through ALMA cycle 0 observations, we have detected the [NII]205um emission in a strongly [CII]-emitting SMG, LESS J033229.4-275619 at z=4.76. The velocity-integrated [NII]/[CII] flux ratio is 0.043 +/- 0.008. This is the first measurement of the [NII]/[CII] flux ratio in high-z galaxies, and the inferred flux ratio is similar to the ratio observed in the nearby universe (~0.02-0.07). The velocity-integrated flux ratio and photoionization models suggest that the metallicity in this SMG is consistent with solar, implying the chemical evolution has progressed very rapidly in this system at z=4.76. We also obtain a tight upper limit on the CO(12-11) transition, which translates into CO(12-11)/CO(2-1)<3.8 (3 sigma). This suggests that the molecular gas clouds in LESS J033229.4-275619 are not affected significantly by the radiation field emitted by the AGN in this system.


Introduction
Investigating the metal content in galaxies is a powerful diagnostic for testing galaxy evolutionary models, because the metallicity is determined by the past star-formation history, gas inflow, and outflow phenomena. The gas metallicity in galaxies has been investigated in galaxies up to z∼3 (Maiolino et al. 2008;Mannucci et al. 2010). However, rest-frame optical metallicity diagnostics (e.g., Nagao et al. 2006a) are not applicable for galaxies at z > 4, where the optical emission lines required to measure the metallicity are redshifted out of the near-infrared atmospheric window (see Maiolino et al. 2008). Moreover, high-z young galaxies with vigorous star formation are often obscured by dust and recognized as, e.g., ultraluminous infrared galaxies (ULIRGs) and submillimeter galaxies (SMGs). In these galaxies, the rest-frame optical lines are significantly affected by the dust extinction and accordingly the optical metallicity diagnostics may lead to large systematic errors. For instance, it has been claimed that ULIRGs show systematically lower metallicities than expected from the mass-metallicity relation (Rupke et al. 2008;Caputi et al. 2008); but this may simply be a consequence of the fact that in these heavily obscured systems the optical metallicity tracers only probe the outer, less enriched regions. Santini et al. (2010) reported that the high dust mass measured in these galaxies disagrees with the metallicity inferred from the optical metallicity diagnostics, suggesting that the optical spectral indices are likely unreliable in these dusty systems.
Metallicity diagnostics exploiting far-infrared fine-structure emission lines are powerful alternatives ). Here we focus on the two strongest fine-structure lines at λ rest > 150 µm, i.e., [C ii] 158 µm and [N ii] 205 µm. The [C ii] 158 µm emission (ν rest = 1900.539 GHz) is one of the strongest lines among the whole electromagnetic spectrum and is consequently observed for some high-z objects even up to z ∼ 6−7 (Maiolino et al. 2005;Walter et al. 2009a;Venemans et al. 2012). This emission arises mainly in H ii regions and photodissociation regions (PDRs), and the relative contribution from those two regions depends on the physical conditions of gas clouds (e.g., Abel et al. 2005). The [N ii] 205 µm emission (ν rest = 1461.132 GHz; Brown et al. 1994) arises in H ii regions (see Nagao et al. 2011). In H ii regions, the [N ii]/[C ii] flux ratio depends mostly on the N/C elemental abundance ratio, while it is relatively insensitive to other gas physical conditions such as the gas density, since their critical densities are similar (n cr = 44 cm −3 and 46 cm −3 for [N ii] 205 µm and [C ii] 158 µm at 8000 K, respectively). Note that the N/C ratio is proportional to Z gas at least for log Z gas /Z ⊙ > −0.5 as observed in Galactic H ii regions (e.g., van Zee et al. 1998) because nitrogen is a secondary element (see also Hamann & Ferland 1999). Since the volume ratio of H ii regions and PDRs also depends on cloud properties such as the gas density and ionization structure (hence the ionization parameter), the observed [N ii]/[C ii] flux ratio accordingly depends on those parameters. Recently, [N ii] 205 µm detections have been obtained for three strongly lensed galaxies at z∼4-5 (Decarli et al. 2012;Combes et al. 2012); however a [C ii] 158 µm detection has not been reported in any of these galaxies.
Here we show that the [N ii]/[C ii] flux ratio is a good metallicity indicator based on photoionization models, and apply this new method to a high-z luminous SMG, LESS J033229.4-275619 (hereafter LESS J0332). We selected this SMG be-

Observations and results
We observed LESS J0332 at z = 4.76, with the band 6 receiver on the Atacama Large Millimeter/submillimeter Array (ALMA) in the dual-polarization setup, as a cycle 0 program. The observation was carried out in three separated runs; on 2011 October 1, 2012 January 12, and 2012 January 13. Each observing run consisted of 1.2 hours and the total observing time is 3.6 hours, including overheads. The receiver was tuned to 253.96001 GHz to cover the redshifted [N ii] 205 µm emission with the upper sideband, whose redshift is based on our previous [C ii] 158 µm detection (De Breuck et al. 2011). The lower sideband was used to cover the redshifted CO(12-11) emission (ν rest = 1381.995 GHz). The correlator was used in the frequency domain mode with a bandwidth of 1875 MHz (488.28 kHz × 3840 channels). The observation was performed with 18 antennas in the compact configuration, but the data obtained with one antenna was flagged out in the last run due to its high system temperature. Callisto was also observed as a flux calibrator. The bandpass and phase were calibrated with J0522-364 and J0403-360, respectively. The atmospheric condition was PWV = 3.0-5.0 mm in the first run and PWV = 1.0-2.0 mm for the last two runs.
The data were processed with Common Astronomy Software Applications (   The best-fit result is parameterized by the peak frequency of 253.868±0.007 GHz, peak intensity of 2.58 ± 0.21 mJy, FWHM of 230 ± 22 km s −1 , and the velocity-integrated flux of 630 ± 78 mJy km s −1 . The redshift inferred from the observed [N ii] 205 µm frequency is z [NII] = 4.7555±0.0002, which is significantly blueshifted from the Lyα redshift. This blueshift is probably caused by the Lyα absorption by the intergalactic neutral hydrogen. The redshift is consistent with the CO(2-1) redshift, though there is a velocity offset, significant at ∼2σ, between the [N ii] 205 µm and [C ii] 158 µm lines. This velocity offset may be a consequence of the modest signal-to-noise ratio in the two spectra, especially in the [C ii] spectrum. Forthcoming [C ii] ALMA observations of the same objects will either confirm or reconcile the velocity discrepancy. However, if this velocity offset is confirmed, this results in different [N ii]/[C ii] flux ratios in the blue and red components (see §3.2).
The redshifted CO(12-11) emission is not detected in our ALMA data. The inferred 3σ upper limit on the velocityintegrated flux is 344 mJy km s −1 (adopting the same aperture size as adopted for the [N ii] 205 µm line and assuming the same velocity profile).
The continuum emission at λ obs ∼ 1.2 mm (i.e., λ rest ∼ 210 µm) is clearly detected with high accuracy. Its flux is 3.5 ± 0.1 mJy with the same aperture as adopted for the [N ii] 205 µm measurement. These measurements were made by combining the upper and lower sideband data and excluding the channels affected by the [N ii] 205 µm emission and by a relatively high noise level at the edges of the sidebands.  Gilli et al. (2011) reported the presence of a Compton-thick active galactic nucleus (AGN) in LESS J0332. If the AGN contributes significantly to the observed far-IR emission lines, its interpretation would become accordingly more complex. Within this context the CO spectral line energy distribution can help, since it is sensitive to the heating energy source. CO lines at high excitation levels are significantly stronger when the molecular gas clouds are affected by the X-ray emission from AGNs than in cases without AGNs (see, e.g., Spaans & Meijerink 2008). The nearby ULIRG-Quasar Mrk 231 shows strong high-J CO lines up to J = 13 − 12, which are properly accounted for by introducing X-ray dominated regions (XDRs) into models (van der Werf et al. 2010). At high-z the gravitationally magnified quasar at z = 3.91, APM 08279+5255 also shows strong high-J CO lines (Weiß et al. 2007;Bradford et al. 2011), which are also well described by XDR models. The CO spectral line energy distribution of other star-formation-dominated high-z galaxies, such as SMM J16359+6612 at z = 2.5 (Weiß et al. 2005) and IRAS F10214+4724 at z = 2.3 (Ao et al. 2008), is completely different from that of APM 08279+5255, showing weaker high-J CO lines (see, e.g., Fig. 14 in Weiß et al. 2007).

Possible AGN contribution
By combining our upper limit on the CO(12-11) flux and the previous measurement of the CO(2-1) flux (Coppin et al. 2010), we obtain a 3σ upper limit on the flux ratio of CO(12-11)/CO(2-1) of 3.8. This upper limit is inconsistent with the CO spectral line energy distribution of the quasar APM 08279+5255, but is fully consistent with other star-formation-dominated high-z objects (see Fig. 14 in Weiß et al. 2007). This suggests that the molecular clouds in LESS J0332 are not described by XDR models, i.e., the AGN contribution to the heating and excitation of the ISM in LESS J0332 is not significant. This is consistent with our earlier study on LESS J0332 (De Breuck et al. 2011), where we estimated that the XDR contribution to the [C ii] 158 µm is ∼1.3%, based on the absorption-corrected X-ray luminosity of L 2−10keV = 2.5 × 10 44 erg s −1 (Gilli et al. 2011) and a scaling relation of L [CII],AGN = 2 × 10 −3 L 2−10keV (Stacey et al. 2010).

Fig. 3. Observed [N ii]/[C ii]
flux ratios compared with model predictions. The green hatched range denotes the observed range for low-z galaxies. The horizontal dashed line shows our ALMA result on LESS J0332, where the emission-line fluxes are measured by integrating the best-fit Gaussian function. Dotted magenta lines denote the flux ratio at the red part (from -350 km s −1 to -150 km s −1 ) and blue part (from -550 km s −1 to -350 km s −1 ) of the lines. The red and blue lines show Cloudy model results as a function of Z gas with log n HII = 1.5 and 3.0 respectively, while solid and dashed lines denote the models with log U HII = -2.5 and -3.5, respectively.  Rangwala et al. 2011). Therefore the [N ii]/[C ii] flux ratio of LESS J0332 is similar to the observed ratios reported for nearby galaxies, suggesting similar Z gas .

Gas metallicity
At high-z [N ii] 205 µm has been detected in HLS J091828.6+514223 at z = 5.24 (Combes et al. 2012), APM 08279+5255 at z = 3.91, and MM 18423+5938 at z = 3.93 (Decarli et al. 2012). However, their [C ii] 158 µm line has not been observed and consequently their [N ii]/[C ii] ratio is unknown (see Walter et al. 2009b and references therein). Therefore our [N ii] 205 µm detection allows us to infer the first measurement of the diagnostic [N ii]/[C ii] flux ratio at highz. Note that there are many [N ii] 122 µm detections in nearby galaxies (e.g., Graciá-Carpio et al. 2011) and also in a few highz galaxies (Ferkinhoff et al. 2011). Although there are attempts to infer the [N ii] 205 µm flux from the [N ii] 122 µm emission (e.g., Decarli et al. 2012), this method may introduce a large systematic error because the flux ratio of [N ii] 122 µm and [N ii] 205 µm varies by a factor of ∼10, i.e., it is strongly dependent on the gas density (see Oberst et al. 2006).
To explore the gas metallicity in LESS J0332 more quantitatively, we carried out model calculations using Cloudy (Ferland et al. 1998) version 08.00. Since the [C ii] line arises in both H ii regions and PDRs, a consistent treatment to connect those two regions is required to investigate the [N ii]/[C ii] flux ratio. We assumed a pressure-equilibrium gas cloud that is characterized by certain gas densities and ionization parameters at the illuminated face (n HII and U HII ) for each model run. Here we examine gas clouds with log n HII = 1.5 and 3.0, and log U HII = -3.5 and -2.5, for Z gas /Z ⊙ = 0.05 -3.0. Note that the gas density in PDRs is higher than log n HII under the pressure-equilibrium assumption. The input continua are the starburst99 constant starformation model spectra (Leitherer et al. 1999) with an age of 1 Myr and a stellar metallicity equal to the gas metallicity. We did not take AGN effects into account in the models since the expected AGN contribution is small (see §3.1). The relative chemical composition of gas clouds is scaled from the solar elemental abundance ratios except for nitrogen, which scales as Z 2 gas at Z gas > 0.23Z ⊙ (see Kewley & Dopita 2002). Orion-type graphite and silicate grains are included. Calculations are stopped at the depth of A V = 100 mag to cover the whole [C ii] 158 µm emitting regions, following Abel et al. (2005). Note that the resulting [N ii]/[C ii] flux ratios are not sensitive to the adopted stopping criterion. The calculation details are given in Nagao et al. (2011).
The [N ii]/[C ii] flux ratio predicted by the model runs are compared with observations in Fig. 3. The models show that the [N ii]/[C ii] flux ratio increases monotonically with Z gas . The moderate dependences on parameters other than Z gas , such as n HII and U HII , prevent us from determining the accurate Z gas value. The observed [N ii]/[C ii] flux ratio (∼ 0.043) infers log(Z gas /Z ⊙ ) ∼ 0.0 ± 0.3; i.e., consistent to the solar metallicity. Although the constraint on the gas metallicity is currently fairly loose, by observing additional lines in the future (such as [N ii] 122 µm and [O i] 145 µm) it will be possible to constrain the gas properties (density and ionization parameter) and therefore determine the gas metallicity much more accurately.
However, we note that the velocity profiles of the [N ii] 205 µm and [C ii] 158 µm emission are different (Fig. 2). In the redder part of the line the [N ii]/[C ii] flux ratio is significantly lower than that in the bluer part. If the different [N ii]/[C ii] ratios for different velocities are caused by metallicity variations, this implies that the system is not chemically homogeneous, but probably resulting from a merging systems where gas with different metallicities has not fully mixed yet. More specifically, the redder component with a high [N ii]/[C ii] flux ratio may be associated with a chemically enriched galaxy that is merging with a metal poor galaxy, associated with the bluer component characterized by the low [N ii]/[C ii] ratio. To investigate this scenario more quantitatively, we measured the [N ii]/[C ii] flux ratio on the redder part (-350 km s −1 < V < -150 km s −1 ) and on the bluer part (-550 km s −1 < V < -350 km s −1 ) of the lines, resulting in a flux ratio of 0.086 and 0.020, respectively. These values roughly correspond to the highest and lowest flux ratios seen in nearby galaxies (Fig. 3), which is indicative of solar-supersolar metallicity in the former case and sub-solar metallicity in the later case. Note that the two distinct [N ii] 205 µm velocity components are seen in another high-z [N ii] 205 µm emitter, HLS J091828.6+514223 (Combes et al. 2012), which may be another example of a chemically inhomogeneous system at high-z.
The low-metallicity component in LESS J0332 was already identified in our previous work on this SMG, where we argued (based at that time solely on the [C ii] 158 µm, FIR and CO properties) that the gas metallicity in LESS J0332 is low (De Breuck et al. 2011). However, the most interesting result obtained here is that this system does also host a significant metal-rich component (with solar-supersolar metallicity), indicating that this system is chemically highly evolved already at z = 4.76 (the cosmic age of 1.27 Gyr).
This result is consistent with past metallicity studies on highz AGNs. The diagnostic flux ratios of some UV metallic permitted lines of type-1 quasars show no redshift evolution of the quasar broad-line region up to z ∼ 6 (e.g., Nagao et al. 2006c;Juarez et al. 2009). The metallicity in narrow-line regions of type-2 AGNs (which trace larger spatial scales and are therefore more closely related to the host galaxy properties than to the broad-line region) also show no redshift evolution (e.g., Nagao et al. 2006b) even up to z ∼ 5.2 . Note that the high-z AGNs investigated in those studies are hosted mostly by massive galaxies (see, e.g., Seymour et al. 2007), similar to SMGs (e.g., Borys et al. 2005). These results on the metallicity of AGNs and LESS J0332 suggest that the chemical evolution of massive systems has progressed very rapidly in the early epoch, z > 5, which is qualitatively consistent with the so-called downsizing evolution seen in galaxies and AGNs (e.g., De Lucia et al. 2006;Ikeda et al. 2011).