Free Access
Issue
A&A
Volume 650, June 2021
Article Number A86
Number of page(s) 24
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/202039295
Published online 10 June 2021

© ESO 2021

1 Introduction

The spatial locations of photon-dominated regions (PDRs) – which are defined as regions where the far-ultraviolet field (FUV) has photon energies ranging between 6.0 eV < < 13.6 eV (Röllig et al. 2013) - and H II regions surrounding stellar nurseries with photon energies 13.6 eV ≤ (Draine 2011) are closely related. As a consequence, and given the low ionization potential of neutral carbon of ~11.2 eV, a fraction of the observed [CII] 158 μm emission does not originate within PDRs but arises from their adjacent H II regions instead, where ionized carbon can be collisionally excited by electrons (Goldsmith et al. 2012). Moreover, Abel (2006) showed that between 10% and 60% of the observed [CII] 158 μm emission can arise from H II regions. This can significantly alter the derived energy balance in PDRs by attributing more energy cooling to them via the [CII] 158 μm line than what they truly emit. More recently, Goldsmith et al. (2015) showed that, along several lines of sight (LOS) in the Galactic plane, velocity-resolved [NII] 205 μm line emission is highly correlated with the [CII] 158 μm line emission in the sense that any significant detection of the former is always accompanied by a detection of the latter, though relative line intensities can vary dramatically. This strongly suggests a physical connection between the observed emission of both species on large scales.

The contribution of a H II region to the observed [CII] emission toward PDRs can be traced by [NII] line observations at 205 μm, as shown by semi-theoretical work in the literature (Heiles 1994; Abel 2006). Given the large nitrogen ionization potential (~14.6 eV), [NII] emission can only originate from within H II regions, removing the origin ambiguity present in [CII] line observations. On the other hand, the correlation between the [CII] 2P3∕22P1∕2 and [NII] 3P13P0 fine structure line transitions within H II regions allows us to estimate the [CII] emission contribution of H II regions to the total observed emission, assuming both transitions emit at the same radial velocities. The intensity of the [CII] line inside a H II region is mainly influenced by the electron density ne, which controls the excitation of ionized carbon, the size of the H II region with respect to the PDR’s size constraining the volume extent from within which [CII] emission can arise, and the intensity of the stellar continuum controlling the [CIII] → [CII] line transition (Abel 2006). Therefore, a semiempirical relationship between I([CII]) and I([NII]) can be established to predict the [CII] emission contribution from H II regions to PDRs across massive star forming regions (Abel 2006). The total observed [CII] emission toward PDRs can be then corrected accordingly, yielding a more accurate energy balance determination in such regions. In this sense, the interplay between the extremely massive and luminous Arches Cluster (Cotera et al. 1996) and the large scale structures known as the Arched Filaments (AF; Pauls et al. 1976) offers a unique possibility to test the validity of such an approach to determine the relative [CII] emission contribution between PDRs and H II regions in the extremely harsh and high-metallicity environment of the Galactic center (GC). In Fig. 1, the 20 cm continuum emission view of the region containing the Arches Cluster and the AF is shown (Yusef-Zadeh & Morris 1987). Further discussion on the figure is featured later in the text.

The Arches Cluster, first discovered by Cotera et al. (1996) and located at a projected distance of ~27 pc from Sgr A, is tightly packed, young, contains a mass of ~2 × 104 M within a radius of ~0.19 pc, and has an estimated age between 2–2.5 Myr (Clark et al. 2018). The bolometric luminosity of its ~160 OB stars is Lbol ~ 108 L, producing a Lyman continuum flux QLyc ~ 1–4 × 1051 s−1 (Figer 2008; Simpson et al. 2007). The measured local standard of rest (LSR) velocity of the cluster is +95 ± 8 km s−1 (Figer et al. 2002), while its proper motion is 172 ± 15 km s−1 (Clarkson et al. 2012), yielding a 3D velocity ~232 ± 30 km s−1 (Stolte et al. 2008). The large 3D velocity of the cluster suggests that it must be in a noncircular orbit if the cluster is to be found within the X2 orbits (Stolte et al. 2008). Given the short orbital time (≲1 Myr) for a stellar cluster to orbit the GC, it is likely that the cluster is no longer spatially related to its parent molecular cloud (Simpson et al. 2007). The interaction of the stellar radiation field of the Arches Cluster with the surrounding ISM has given rise to arc-like structures known as the AF (Simpson et al. 2007), shown in Fig. 1. Their radio continuum emission originates mainly in a thermal plasma of electron temperature Te ~ 7000 K, which is derived from hydrogen recombination lines and radio continuum observations (Pauls et al. 1976; Lang et al. 2001). García et al. (2016) showed that the brightest [CII] and [NII] emission around this region is found within the AF, at negative radial velocities. These structures are closely associated to a system of dense molecular clouds traced by CS(2–1) emission (Bally et al. 1987, 1988). Previous work by Genzel et al. (1990) used low spectral (67 km s−1) and spatial (55″) [CII] line observations in combination with [OI] at 63 μm, mid-, far-infrared (FIR), and radio continuum observations to model the emission from the massive (~2 × 104 M) G0.07+0.04 region within the AF in the context of PDRs. Their results were consistent with a PDR nature of the emission and they interpreted the [CII] line observations as coming from partially ionized interfaces between fully molecular and fully ionized gas.

The main goal of the present work is to assess whether the semi-theoretical predictions in Abel (2006) for the [CII] emission contribution from H II regions to nearby PDRs, derived for other star forming regions across the Galactic disk such as the Orion Nebula, are applicable to the AF region in the GC, where the ISM physical conditions are much more extreme (Mezger et al. 1996; Güsten & Philipp 2004). For this purpose, the behavior of the observed [CII] versus [NII] emission in the AF (García et al. 2016) is analyzed. In Sect. 2, a summary of the observations used in the present work is given. In Sect. 3, we discuss the criteria used to overcome the intrinsic difficulties involved in the [CII] line profile interpretation, mainly due to the presence of emission and absorption features along the LOS and opacity uncertainties. In Sect. 4, the applicability of the Abel (2006) analysis under the extreme physical conditions in the GC is investigated. We explore how local variations on the carbon-to-nitrogen elemental abundance ratio [C/N] can explain the overall behavior of the emission. In Sect. 5, we present the derived [CII] emission contribution from the H II region component to the observed [CII] emission in the AF at the selected positions in our analysis. Using these results, in Sect. 6 we derive the average [C/N] elemental abundance ratio for the AF combining independent estimates of electron temperatures Te and electron volume densities ne with our observations. Also, by assuming a carbon elemental abundance [C/H] consistent with values found toward the GC, an average nitrogen elemental abundance [N/H] is estimated for the AF region and compared with other values in the literature found across the GC. In Sect. 7, our main results are summarized.

2 Summary of the observations

We made use of the [NII] 3P13P0 at 205 μm (1462.1 GHz) and [CII] 2P3∕22P1∕2 at 158 μm (1900.5 GHz) fine structure line observations from García et al. (2016) covering the AF region. These observations were part of the Sgr A complex observations obtained with the HIFI heterodyne receiver (Roelfsema et al. 2012) onboard the Herschel Satellite (Pilbratt et al. 2010) as part of the Herschel guaranteed time HEXGAL key program (PI: R. Güsten, MPIfR) in which the [CI] 3P13P0, [CI] 3P23P1, [NII] 3P13P0, and [CII] 2P3∕22P1∕2 submillimeter lines were measured. The HIFI receiver was a double sideband (DSB) single pixel receiver covering the frequency range from 480 to 1250 GHz and from 1410 up to 1910 GHz. The [NII] 3P13P0 and [CII] 2P3∕22P1∕2 lines (hereinafter referred to as [NII] and [CII]) were covered by the receiver bands 6a and 7b, respectively.Two wide-band acousto-optical spectrometers (WBS) were used as backends, measuring vertical andhorizontal polarizations. The Herschel beam sizes at the [NII] and [CII] frequencies are 14.5″ and 11.2″, respectively, with corresponding spectral resolutions of 0.103 and 0.079 km s−1. We used the on-the-fly (OTF) map with load chop and position-switch reference observing mode of Herschel-HIFI to measure the [NII] and [CII] line intensities. This observing mode is very efficient in spending a large fraction of the OTF scanning time ON source, with a total OTF scanning time shorter than the Allan stability time for these high-frequency lines (Roelfsema et al. 2012). The downside is that it produces under-sampled maps in the OTF scan direction as shown in Fig. 3.5 by García (2015). The pointing accuracy of the telescope was estimated to be ~2″, and the calibration error in the measured intensities for bands 6a and 7b is ~25% (Roelfsema et al. 2012; García et al. 2016). For the present work, all data sets have been resampled to a common 1 km s−1 spectral resolution and spatially smoothed to 46″ spatial resolution. This allowed us to further increase the signal-to-noise ratio (S/N) in the [NII] and [CII] data sets, and to match the lowest spatial resolution of the data in García et al. (2016) obtained for the [CI] 3P13P0 line observations (hereinafter referred to as [CI](1–0)). Given the observing mode used to measure the [NII] and [CII] line transitions, a considerable amount of emission is interpolated in this process. We used the TA*$T_{\textrm{A}}^{*}$ antenna temperature scale for the reported spectra, which is the best approximation of the true convolved antenna temperature measured in the region given the error beam pickup level in the HIFI 6a and 7b bands, and the large spatial extension of the emission (García et al. 2016; García 2015).

In order to investigate the [CII] line profile and the number of radial velocity components along the LOS at the targeted positions in this work, we included the [CI](1–0) and [CI](2–1) line observations from García et al. (2016) in our analysis. For the same purpose, complementary C18 O(1–0) observations obtained with the EMIR1 receiver (Carter et al. 2012) at the IRAM 30m telescope under project number 020-16 are also used. These data sets were smoothed to the common 1 km s−1 and 46″ spectral and spatial resolutions shared by the [CI](1–0), [NII], and [CII] line observations. The C18 O(1–0) data sets are part of a different and ongoing work, so the corresponding technical information will be published in more detail in a future paper. The [CII] and [NII] line intensities are weak compared to what is observed toward other massive star forming regions in the Galactic disk, where [CII] peak line intensities are well above 50 K (Ossenkopf et al. 2013). For our selected positions in the AF, the [CII] line is not brighter than ~8 K, despite the fact that the AF are the brightest structures in [CII] emission within the Sgr A complex. The low-level emission is particularly troublesome for the [NII] line observations for which the resampling to a 46″ spatial resolution increases the S/N only to yield detections with ≳3σ significance level, while for the [CII] line observations, the S/N increase allows us to study the line profile in more detail. Despite the aforementioned limitations, the [NII] and [CII] data sets in García et al. (2016) constitute the most extended velocity-resolved ones available for these transitions in the Sgr A complex at the time of publication.

thumbnail Fig. 1

AF region outlined in 20 cm continuum emission from Yusef-Zadeh & Morris (1987). Selected positions for the I([CII]) versus I([NII]) integrated intensity analysis in this work: E1, E2-N, E2-S, W1, W2, G0.07+0.04, G0.10+0.02, P1, and P2, are shown as open squares in the García et al. (2016) data grid. The 46″ spatial resolution of their data is represented by open circles. Filled triangles depict either the closest continuumemission peak to each position or one of the two CS(2–1) peaks (P1 or P2) in the observations by Serabyn & Güsten (1987). Panel insets contain the [CII] (blue), [NII] (red), and scaled C18 O(1–0) (filled-gray) spectra at the positions selected for this work, showing very complex line profiles with large variations in line shape and brightness across the region. Black crosses show the positions of several H II regions labeled H1–H13 spanning a region south of the filaments (the so-called H-Region). The ⋆ symbol shows the position of the Arches Cluster, ☆ while (text here emptystar) symbols indicate the locations of some of the massive (O, B, and WR) field stars scattered across the region (Dong et al. 2012, 2015).

3 Complexity of line profiles

The overall view of the area containing the Arches Cluster and the AF traced by 20 cm continuum observations (Yusef-Zadeh & Morris 1987) is shown in Fig. 1. The region contains several small H II regions (labeled H1–H13; Lang et al. 2001; Yusef-Zadeh & Morris 1987; Zhao et al. 1993), and O, B, and WR field stars (Dong et al. 2012, 2015). The selected positions (open circles) for our analysis in Sect. 5 are named after their association with the northern (N) or southern (S) parts of the E1, E2, W1, and W2 AF segments (e.g., E1, E2-N, E2-S, W1, and W2), other known sources (e.g., G0.10+0.02 and G0.07+0.04), or high-density gas identified in CS rotational line observations (e.g., P1 and P2 Lang et al. 2001; Timmermann et al. 1996; Serabyn & Güsten 1987). The E1, E2-N, E2-S, and G0.10+0.0 positions trace the two sub-filaments forming the eastern (E) filament, while the W1, W2, and G0.07+0.04 positions trace the western (W) filament. The G0.07+0.04 region is particularly interesting since the gas there is thought to be interacting with the Northern Thread filament (NTF), a structure identified in high-resolution 6 cm continuum observations (Lang et al. 1999b). The local magnetic field is aligned with the thread, indicating the possible interaction of the ionized, neutral atomic and molecular material with the pervasive large-scale, dipolar magnetic field in the GC. Most of the selected positions are seen as local maxima in the 3.6 cm continuum observations by Lang et al. (2001), while P1 and P2 are located at the edges of the E and W filaments. The criteria for selecting these positions is based on the [CII] and [NII] emission distribution shown in García et al. (2016), which is most intense along the AF and closely follows the 20 cm continuum emission in the region (Yusef-Zadeh & Morris 1987), and on the location of the high-density gas, where PDR-like emission is mostly expected if the stellar FUV field is the main excitation source of the gas (Draine 2011). The measured [CII], [NII], and C18 O(1–0) spectra for all positions are displayed in the corresponding panel insets, while their equatorial coordinates are listed in Table 1.

When observing the GC, the LOS crosses a large fraction of the Galactic plane. Hence, related structures such as spiral arms appear in the measured spectra as either deep absorption or narrow bright emission features, depending on the observed species. Typical LSR radial velocities of spiral arms in the LOS toward the GC are approximately −55 km s−1 for the 3-kpc expanding arm, approximately −30 km s−1 for the 4-kpc arm, and approximately −5 km s−1 for the local arm (Oka et al. 1998; Jones et al. 2012; García et al. 2016). On top of this, the emission originating from within the GC volume is extremely broad and is usually detected between ±200 km s−1 at Galactic longitude l = 0° (García et al. 2016), making the identification of individual gas components in the observed spectra particularly troublesome due to source blending, varying opacity effects across the spectra, different gas excitation conditions, etcetera. Moreover, in the case of the [CII] emission this effect is worsened as C+ can exist in a plethora of environments given its low ionization potential, while its 1.9 THz fine structure transition can be excited in different phases of the ISM, from diffuse clouds to PDRs (Langer et al. 2010, 2016; Pineda et al. 2014), making the identification of its exact origin along the LOS in broad, blended, and complex line profiles a very challenging task. In order to illustrate the complexity of the emission detected in the AF, Fig. 2 shows the [CII], [CI](1–0), and C18 O(1–0) spectra for all positions in Table 1 tracing the molecular, atomic, and ionized content of the ISM in the GC. From the figure, bright and narrow components at negative radial velocities seen in the optically thin C18 O(1–0) emission (marked by vertical-dotted lines) are superimposed on a much broader emission distribution along the radial velocities.

In order to identify the physical components in the GC traced by the [CII] emission, a natural choice would be to target the optically thin [13CII] F = 2 →1, 1 →0, 1 →1 hyperfine triplet. The selection of this isotope would drastically decrease the confusion between different components along the LOS, foreground emission and absorption features, etcetera, as the conditions for sufficient line excitation (mainly warm, high-column-density gas) are not easily met across the entire Galactic plane. Unfortunately, these diagnostic lines are not easily detectable in the GC for two main reasons. Firstly, the small triplet frequency separation, equivalent to +11.2 km s−1, −65.2 km s−1, and +63.2 in radial velocities with regard to the rest frequency of the main isotope at 0 km s−1 (with relative intensities 0.625, 0,250, and 0.125, respectively) are located too close to the rest frequency of the main isotope to avoid strongline blending, given the extremely broad [CII] line profiles shown in Fig. 2. Secondly, even if there were a position where the [CII] line became narrow enough, the peak [CII] line intensities in Fig. 1 are quite low (<10 K) when compared with measurements in other massive star forming regions in the Galactic disk (Guevara et al. 2020). The detection of even the strongest [13CII] hyperfine component would require prohibitively long integration times on the SOFIA Airborne Observatory, which is currently the only astronomical facility capable of observing the [CII] line. On the other hand, [13 CII] observationsof massive star forming regions outside the GC such as M43, Mon R2, M17 SW, Orion Bar, and NGC 2024 have shown that opacity effects are a major issue when interpreting [CII] line profiles, sometimes with up to 90% of the line intensity absorbed away in the measured spectra (Graf et al. 2012; Ossenkopf et al. 2013; Guevara et al. 2020). This is a major problem for the [CII] data analysis, potentially causing a strong underestimation of the gas cooling rate and the misidentification of radial velocity components along the LOS, which causes misleading model results. Therefore, in dealing with any [CII] data interpretation, this issue must be addressed. Fortunately, this problem can be overcome in the GC given the extreme ISM physical conditions and the low brightness of the [CII] line, in contrast to what it is observed in the Galactic disk. In the absence of [13 CII] observations for the AF, other complementary data can be used to assess the impact of opacity effects on the observed [CII] spectra. Experimentally, the C18 O(1–0) line seems to closely follow the [13CII] line profile in massive star forming regions such as M43, Mon R2, and M17 SW in the Galactic disk (Guevara et al. 2020). Therefore, this optically thin transition could potentially be used as a template to determine the number of [CII] physical components contained in each spectrum, when no [13CII] observations are available. Figure 2 shows a remarkable similarity between the [CI](1–0) and C18 O(1–0) line profiles, except for some of the very narrow spectral features, most likely tracing cold gas in the foreground, where the [CI](1–0) transition is not very bright. Such a similarity implies that opacity effects do not affect the observed [CI](1–0) transition very strongly. At the same time, the [CII] line shape at W2, W1, P2, G0.10+0.02, E2S, and E1 closely follows that of the more optically thin C18O(1–0) and [CI](1–0) lines. This provides empirical support to the idea that opacity effects in the [CII] spectra at the selected positions can not be the dominant factor responsible for the observed line shapes. Moreover, the fact that all these positions, which are several parsecs away from each other, exhibit the same spectral behavior in these transitions suggests this as a common feature across the entire AF region.

At positions P1, G0.07+0.04, and E2N, the [CII] line shapes are also similar to those of C18 O(1–0) and [CI](1–0), but they are largely (>5 km s−1) blueshifted toward more negative radial velocities. The difference in central radial velocities (obtained from Gaussian fits in Sect. 3 and displayed in Appendix A) between the [CII] and C18 O(1–0) radial velocity components identified in the spectra are approximately −13 (P1-46 km s−1), −6 (P1-13 km s−1), −8 (G007-41 km s−1), and −9 km s−1 (E2N-18 km s−1). On the other hand, at position E2S, although a visual inspection of Fig. 2 shows a similar emission distribution of these two lines, the central radial velocity of the [CII] line shows a large red-shift of approximately +7 km s−1 (E2S-17 km s−1) with respect to what is obtained for the C18O(1–0) line. For therest of the radial velocity components identified at all positions, the difference in central radial velocities between both lines is <5 km s−1. A large shift in central radial velocities between ionized and neutral species could be attributed to geometry effects or local interactions with the magnetic field in the GC, rather than local opacity effects. For instance, the G0.07+0.04 source is known to exhibit a large radial velocity blueshift of the emission produced by ions (e.g., H+ and C+) with regardto the line profiles produced by neutral species such as C, CO, and C18 O (Genzel et al. 1990), which is consistent with what is shown in Fig. 2. Two interpretations for the blueshifted emission at G0.07+0.04 have been given in the literature. On the one hand, the radial velocity shift between the H92α line, at approximately −44 km s−1 (Lang et al. 2001), and the molecular emission could be a consequence of the location of the AF ionized gas at the near side of a molecular cloud (Serabyn & Güsten 1987; Lang et al. 2001). Genzel et al. (1990) found that their velocity-unresolved observations of [CII] and H110α recombination lines are blueshifted by 10–15 km s−1 with respect to their CO and CS observations, indicating streaming motions between ionized and molecular gas. In that context, the increasing magnitude of the radial velocity blueshift, from the emission originating in neutral molecular/atomic gas to the [CII] emission and further to the H92α recombination line, would reflect a layered spatial distribution along the LOS of the gas at this position. On the other hand, it has been proposed that G0.07+0.04 is interacting with the nonthermal NTF, whose intrinsic magnetic field is aligned with the filament (Yusef-Zadeh & Morris 1988; Lang et al. 1999b). Such interpretation is mainly based on the brightness discontinuity in 20 cm continuum emission at the location where both regions intersect, as shown in Fig. 23 in Lang et al. (1999b). The NTF is expected to have a much stronger magnetic field than its surrounding diffuse ISM (Yusef-Zadeh et al. 2007), so the interaction between the ionized gas at G0.07+0.04 and the NTF magnetic field could cause the observed radial velocity shift. Despite of the observational evidence supporting this scenario, there is also observational evidence against it. Using H92α and 8.3 GHz continuum observations, Lang et al. (2001) found that the intersection region between G0.07+0.04 and the NTF shows no discontinuities concerning the rest of the W1 filament in terms of the gas velocity field, line width, or line-to-continuum ratio. Hence, the exact nature of the interaction between both regions is still unclear, if it exists at all (Lang et al. 1999b, 2001).

Another strong argument against large opacity effects in the observed [CII] line shape is given by Goldsmith et al. (2012). It is based on the fact that, for [CII] line intensities ≲8 K, the emission falls into their so-called effective optically thin (EOT) limit where the I([CII]) integrated intensity is proportional to the C+ column density along the LOS, irrespectively of the optical depth of the line. This is indeed the case for all positions in the present work, as can be seen from the antenna temperature scale of the spectra in Fig. 1. The optically thin behavior of the [CII] line in the GC resembles what is observed in the CO(1–0) transition when used to trace the large-scale mass distribution in the Milky Way’s Galactic disk. Despite the fact that locally the CO(1–0) line is known to be optically thick, at large scales it behaves as if it were optically thin, tracing most of the molecular mass in giant molecular clouds (GMCs; Dame et al. 1986; García et al. 2014). The most common interpretation for this feature is the presence of a large collection of unresolved small CO(1–0) clumps within the telescope’s beam, moving at differential radial velocities, which when observed yield an effectively optically thin transition. Based on all the aforementioned evidence, in the following we assume that spectral features in the [CII] data are mostly due to the existence of different radial velocity components along the LOS, and opacity effects play a less important role in shaping the observed line profiles compared to what is indeed the case in other regions of the Galactic disk.

In order to estimate integrated intensities of the [CII] and [NII] lines, we follow two different approaches depending on the data analysis. In Sect. 4, integrated intensities are obtained by simply adding the measured intensities in the observed spectra (in units of K km s−1) over different LSR radial velocity ranges, which were then converted into units of ergs s−1 cm−2. In this way, the overall behavior of the I([CII]) versus I([NII]) relationship in the AF region can be explored. For the data analysis in Sect. 5, we follow the discussion in Sect. 3, focusing on identifying the radial velocity components traced by the [CII] spectra, using the optically thin C18 O(1–0) line, and to a minor extent the [CI](1–0) line, as a template. The procedure is outlined here:

  • 1.

    We identify all velocity components traced by the optically thin C18 O(1–0) transition. Broad (ΔV (FWHM) ≥ 10 km s−1) and high SN (>5) detected components are assumed to originate within the GC, given the large column densities and turbulence of the gas in the region (Mezger et al. 1996; Güsten & Philipp 2004). Narrow features with radial velocities close to those of gas associated with spiral arms along the LOS are assumed to be in the foreground and are excluded from the analysis. With this information, Gaussian profiles (assumed for simplicity) are fit to obtain the central radial velocity and velocity width of each component.

  • 2.

    The radial velocity information is used to perform Gaussian fits to identify the same components in the [CII] spectrum. Due to the complexity of the spectra, usually multi-Gaussian profiles are fit, as the figures in Appendix A show. Given that the C18 O(1–0) and [CII] lines trace gas under very different physical conditions, it is not always the case that each detected component (narrow or broad) in either line has a corresponding counterpart in the other. Therefore, in order to minimize misidentifications, the [CI](1–0) and [CI](2–1) data in García et al. (2016) are used to check the consistency of the identified components in the [CII] spectra. From this approach, three I([CII]) integrated intensities are obtained for each [CII] component: (1) the integral over the fit Gaussian profile I([CII])g; (2) the sum of measured intensities (in K) multiplied by the channel width (in km s−1) over the relevant spectral channel range listed in Table 2 where the component is found I([CII])s and indicated by the vertical black dotted lines in the figures of Appendix A; and (3) the integrated intensity in the LSR radial velocity −90 km s−1 to +25 km s−1 range where most of the emission in the AF region is found I([CII])i (Simpson et al. 2007; Lang et al. 2001; García et al. 2016), calculated in the same way as I([CII])s. A graphic description of the approach followed to derive I([CII])g, I([CII])s, and I([CII])i is shown in Fig. A.1. After these definitions, the greatest differences between I([CII])g and I([CII])s appear when radial velocity components are largely blended.

  • 3.

    Given the different nature of the [NII] line originating exclusively inside the H II region (i.e., no emission contamination of spiral arms along the LOS is expected) and the poorer S/N of the [NII] data set, only I([NII])s integrated intensities are calculated over the same radial velocity ranges in which I([CII])s intensities were obtained (vertical dotted black lines for figures in Appendix A) as the Gaussian fitting approach in this case is likely not reliable. Finally, the integrated intensity I([NII])i is obtained in the same way as for I([CII])i.

All derived Ig, Is, and Ii integrated intensities were then converted into units of ergs s−1 cm−2. A total of twelve radial velocity components are identified and listed in Table 2 within the nine selected positions in Table 1. The name of each radial velocity component consists of the position’s name followed by the radial velocity obtained from the Gaussian fit to the [CII] spectrum. All details pertaining to deriving integrated intensities used for the data analysis are given in Appendix A.

Table 1

Selected positions in the AF region shown in Fig. 1.

thumbnail Fig. 2

Comparison between [CII] (filled light gray), [CI](1–0) (gray), and C18O(1–0) (black) line profiles for all positions in Table 1. Relative intensities among the spectra have been scaled for display purposes. The observations trace different phases of the ISM from ionized, neutral, and molecular gas. Narrow and strong emission and absorption features in the C18 O(1–0) optically thin line common to all positions are marked with vertical dotted lines as they are most likely foreground components located in the Galactic disk somewhere along the LOS and outside the GC.

Table 2

I([CII]) and I([NII]) integrated intensities in units of ergs s−1 cm−2 for the velocity components identified at the positions in Table 1.

4 Observational and theoretical I([CII]) versus I([NII]) relationship

The spatial distribution of the I([CII])i/I([NII])i integrated intensity ratio in the AF region is shown in Fig. 3. For each ratio, only integrated intensities with detections above the 3σ significance level were considered, where σ is the propagated noise in the radial velocity integration interval (σ = TA,rms*$T_{\textrm{A,rms}}^{*}$ × ΔV×N$\Delta V\,{\times}\,\sqrt{N}$, with ΔV the radial velocity resolution and N the number of spectral channels). Overplotted on Fig. 3 is the 20 cm continuum emission of the region and the selected positions from Table 1.

There are several interesting features in Fig. 3. First, the spatial distribution of the integrated intensity ratios is not homogeneous, with the bulk of the values in the range 0 < I([CII])i/I([NII])i < 10. These ratios are enclosed between two prominent lanes of larger ratios with I([CII])i/I([NII])i > 10 that are close (at least in projection) to the southern edge of the radio arc and the northern edge of the Sgr A-east supernova remnant (SNR). Second, at G0.07+0.04 there is a local enhancement of the integrated intensity ratios with regard to its surrounding area. While there is no definitive evidence for the interaction of G0.07+0.04 with the NTF tracing the local magnetic field in the GC (Lang et al. 1999b), such a localized enhancement suggests that the gas is undergoing a different physical process at this location than in the rest of the region. In the opposite sense, the northern part of the W1 filament demonstrates a notorious decrease in the I([CII])i/I([NII])i ratio, where W1 intersects the large ratio lane south of the radio arc. At this location, there are no evident structures tracing the local magnetic field, either in the 20 or 6 cm continuum observations (Lang et al. 1999b), contrary to what is observed at the G0.07+0.04 location. Therefore, it is likely that the physical process affecting the I([CII])i/I([NII])i ratio at the intersection is not strongly related to a local enhancement of the magnetic field. Third, the same lane of high I([CII])i/I([NII])i values intersects the E1, E2S, and E2N filaments without altering the spatial continuity of such large ratios. This could be indicative of a geometric effect in which the plasma responsible for the continuum emission (Pauls et al. 1976; Lang et al. 1999b) is spatially decoupled from the gas emitting in these submillimeter spectral lines. The prominent lane of I([CII])i/I([NII])i ratios north of the Sgr A-east SNR seems to follow the edge of its continuum emission over the entire map, with a decrease in the ratio toward a local 20 cm continuum maximum at (+50″, +300″) and the location of the H3 region (shown in Fig. 1), similarly to what is seen toward W1. The upper panel in Fig. 3 shows the histogram of the I([CII])i/I([NII])i ratios in the map. A Gaussian profile was fit to characterize the distribution, from which the mean ratio Rc = 5.1 ± 0.1 and its standard deviation Rσ = 2.0 ± 0.1 were derived. The vast majority of the I([CII])i/I([NII])i ratios seem to be relatively consistent, with a Gaussian-like distribution suggesting that there are no dramatic changes in the physical conditions of the bulk of the gas. Nonetheless, a more detailed interpretation of the I([CII])i/I([NII])i ratio spatialdistribution in Fig. 3 is challenging as the PDR-H II contribution to the observed [CII] emission (or from foreground gas along the LOS) is not disentangled and the [NII] emission traces only the H II region.

Models for the relative contribution of H II regions to the observed [CII] emission toward PDRs have been presented previously in the literature (Heiles 1994; Abel 2006). In the models of Abel (2006), the [CII] 2P3∕22P1∕2 emission coming from H II regions is traced by [NII] 3P13P0 line observations (Heiles 1994), as the power-law fit in their Eq. (1) shows. This reflects the amount of [CII] emission expected to arise from an H II region, while the rest of the emission above the prediction is attributed to the adjacent PDR, under the assumption of pressure equilibrium between both regions. To investigate the relationship between the [NII] and [CII] emission in the AF, Fig. 4 shows the evolution of the I([CII]) versus I([NII]) integrated intensities as a function of radial velocity integration interval, increasing from the upper-left to the bottom-right panel, where I([CII])i versus I([NII])i values are plotted. Given the complexity of the line profiles, the selection of any particular radial velocity interval to explore the behavior of I([CII]) versus I([NII]) integrated intensities could potentially affect the true scatter relationship in the GC, as different contributions to the [CII] emission (e.g., from foreground sources such as spiral arms)could be included in (or excluded from) the derived values. In Fig. 4, the Abel (2006) model predictions of the amount of [CII] emission expected to originate within H II regions for metallicities Z = 0.1–1.0 Z are overplotted as dotted and dashed lines, respectively. The vast majority of the I([CII]) versus I([NII]) data points move along the Abel (2006) predictions for Z = 0.1 Z. The few positions where the [NII] emission decreases with increasing radial velocity range are dominated by the rms noise in the spectra, artificially altering I([NII]). The overall behavior of the I([CII]) versus I([NII]) scatter relationships in Fig. 4 shows that the emission from both species is correlated (independently of the radial velocity range) as predicted by the Abel (2006)models, but it does not seem to be “calibrated” for the GC in the sense that most of the data points fall below or along solar and subsolar metallicities, respectively, while in the GC super-solar metallicities from Z = 1–3 Z are measured (Mezger et al. 1979, 1996; Cox & Laureijs 1989; Shields & Ferland 1994; Yusef-Zadeh et al. 2007). In the following, we argue against some of the reasons that could be responsable for the observed discrepancy between model and observations such as:

thumbnail Fig. 3

Bottom: spatial distribution of the I([CII])i/I([NII])i integrated intensity ratio calculated for the radial velocity range −90 km s−1 to +25 km s−1. The AF traced by 20 cm continuum observations and the selected positions in Table 1 are shown as black contours and white circles, respectively. The white ⋆ symbol represents the location of the Arches Cluster. Top: histogram of the spatial distribution of I([CII])i/I([NII])i integrated intensity ratios shown in the map. A Gaussian fit was performed to the histogram to characterize the central ratio (Rc) and standard deviation (Rσ) of the distribution.

Choice of antenna temperature scale

García et al. (2016) showed that the most reasonable antenna temperature scale to apply to their data is TA*$T_{\textrm{A}}^{*}$, in which the calibrated data have been divided by the telescope’s forward efficiency (Feff). Nonetheless, this represents a lower limit to the true convolved antenna temperatures. Since the main beam efficiencies (Beff) for the [CII] and [NII] line transitions observed with the Herschel-HIFI satellite are virtually identical around ~0.58 (Müller & Jellema 2014), scaling the integrated intensities of both lines by Beff /Feff to change the antenna temperature scale to main beam (Tmb), which is an upper limit to the true convolved antenna temperatures, has no impact on the overall trend in Fig. 4.

Underestimation of I([CII]) integrated intensities

The [CII] emission at the G0.07+0.04 position was measured by Genzel et al. (1990) as ~1.2 × 10−3 ergs s−1 cm−2 within a 55″ beam (relative calibration error 30%). This is somewhat lower than what we estimate by integrating the spectrum in the velocity range −90 km s−1 < Vlsr < +25 km s−1 in a 46″ beam, yielding ~1.9 × 10−3 ergs s−1 cm−2 (relative calibration error of 25%, García 2015). Nonetheless, and despite the slightly different beam sizes, both estimates are consistent within the error uncertainties, which rules out a systematic underestimation of I([CII]) integrated intensities.

Overestimation of I([NII]) integrated intensities

In the AF, the [NII] emission is much weaker than the [CII] emission, yielding a lower S/N in the [NII] data than that in the [CII] line observations. For example, at G0.07+0.04 the rms noise of the [NII] spectrum (~0.34 K) yields SN ~ 5, while for the [CII] spectrum, the rms noise (~0.42 K) yields SN ~ 18. Since integrated intensities are calculated mostly for positions where SN > 3 for each line, overestimating I([NII]) is rather unlikely.

C++ instead of C+

Simpson et al. (2007) showed that strong shocks from the Quintuplet Cluster Wolf-Rayet stars can ionize oxygen to [OIV] (ionization potential ~54.9 eV, Draine 2011) and destroy dust grains. If this were also the case for carbon in the AF and most of the ionized carbon were in the form of C++ (ionization potential ~24.4 eV, Draine 2011) rather than in C+ due to a strong ultra-violet (UV) field from the Arches Cluster, the low [CII] line intensities with respect to the predicted ones from the [NII] line intensities could be explained. Nonetheless, this would not explain why this is still the case at positions far away from the Arches’ location such as W1 and W2 where, given their large projected distances to the cluster (>10 pc), the UV radiation field should be significantly reduced. Therefore, this possibility seems rather unlikely.

With none of the aforementioned reasons able to fully account for the behavior of the I([CII]) versus I([NII]) scatter relationship in Fig. 4, we explore the validity of the main underlying assumptions in the Abel (2006) theoretical work to account for it. There, the simplified case of an H II region and PDR in pressure equilibrium, for a single massive star, elemental abundances for the Orion Nebula, and a constant [C/N] ratio, are assumed. By the time Abel (2006) was published, their model was quite advanced in the sense that no other theoretical code could handle both H II and PDR regions. In terms of their assumed elemental abundances, variations in the [C/N] ratio were not investigated. Following Heiles (1994), it is possible to write I([CII])i/I([NII])i as a function of the [C/N] abundance ratio, assumed to be ~[N(C+)/N(N+)], and the stellar temperature T as follows: I([CII])iI([NII])i4.8×[N(C+)N(N+)]×(T104)0.15.\begin{equation*}\frac{I({[CII]})_{i}}{I({[NII]})_{i}} \approx 4.8\,{\times}\,\Bigg[\frac{N(C^{&#x002B;})}{N(N^{&#x002B;})}\Bigg]\,{\times}\,\Bigg(\frac{T_{\star}}{10^{4}}\Bigg){}^{0.15}. \end{equation*}(1)

From Eq. (1), it is clear that the [C/N] abundance ratio is the dominant factor controlling the expected I([CII])i/I([NII])i ratio, while it depends only marginally on the assumed stellar temperature. Moreover, Abel (2006) showed that the intensity ratio is relatively insensitive to ionization parameter U (defined as the dimensionless ratio of hydrogen-ionizing flux to hydrogen density), and as a consequence to the assumed stellar continuum, which is a function of U and T.

The abundance ratio dependence in Eq. (1) is specially interesting since it suggests that decoupled carbon and nitrogen elemental abundances in the GC could be driving the observed scatter relationship in Fig. 4. For instance, the carbon abundance is thought to be between 3 (Arimoto et al. 1996) and 10 times the solar value (Sodroski et al. 1995; Oka et al. 2005) in the GC, while the nitrogen elemental abundance is poorly constrained (Rodríguez-Fernández, & Martín-Pintado 2005; Simpson et al. 1995). At the same time, secondary nitrogen production has been found in other Galactic nuclei (Liang et al. 2006; Pérez-Montero et al. 2011). Liang et al. (2006) examined the oxygen abundance and [N/O] abundance ratio in nearly 40.000 metal-rich (8.4 < 12 + log (O/H) < 9.3 or equivalently 0.6 < Z < 4.5) star-forming galaxies and their nuclei. They reported an increasing [N/O] abundance ratio with increasing 12 + log (O/H) values, interpreting such a trend as being consistent with a combination of primary nitrogen production (from massive stars forming in primordial gas) and secondary (from intermediate to low-mass stars), the latter being the dominant process in high-metallicity environments. Pérez-Montero et al. (2011) found that nitrogen is overabundant in blue compact dwarf galaxies, which is something that cannot be explained solely by the presence of WR stars, pointing to a secondary production of nitrogen. Florido et al. (2015) found observational evidence that bars in galaxies enhance the [N/O] abundance ratio, with log (N/O) ~ 0.09 dex larger than in unbarred galaxies. The presence of a bar also tends to enhance dust concentrations, star formation, and electron densities in the central regions of galaxies. In their work, they suggest that secondary nitrogen production could explain, at least qualitatively, the higher [N/O] abundance ratio observed in bulges less massive than ~1010 M, which puts the central bar of the Milky Way into this category (Mezger et al. 1996).

Since the Abel (2006) calculations did not investigate the variation in [C/N] abundance ratio and did not include SEDs consistent with stellar clusters, we broadened their calculations to investigate the model sensitivities to these parameters. The revised calculations use a Starburst99 SED (Leitherer et al. 2010), which is more appropriate for gas irradiated by a stellar cluster, and they use a range of ionization parameters varying from log U = −4 to −1, in increments of 1 dex. Our choice of SED is a 4 Myr instantaneous starburst. Our calculations include two metallicities (1 Z and 2 Z) and two different hydrogen densities (2.5 and 3.5 in log n(H) units). Elemental abundances in the model range from [C/H] = (0.13−1.00) × 10−3 and [N/H] = (0.40−3.16) × 10−4 which include the corresponding Orion Nebula (ON) values in Simón-Díaz & Stasińska (2011) [C/H]ON = 2.334 × 10−4 and [N/H]ON = 0.832 × 10−4, which are taken (in this work) as the standard elemental abundances for the Galactic disk. Thus, the model [C/N] ratios vary between 0.84 and 12.65, whereas [C/N]ON = 2.81 is obtained forthe Galactic disk. Carbon and nitrogen abundances for each model are scaled between 0.5 and 2, in increments of 0.25. Scaling them in this way implies that our models vary the [C/N] abundance ratio by a factor of 64. Our calculations include only the H II region, which we consider to end when the hydrogen ionization fraction falls below 0.01. Our results are shown in Fig. 5, where the model [C/N] abundance ratios for the AF ([C/N]AF) have been normalized to the Galactic disk value. For comparison, the (extrapolated) Heiles (1994) linear correlation between I([CII]) versus I([NII]) derived for the extended low-density warm ionized medium (ELDWIM; Eq. (12b) in their work) is shown. A direct comparison between the measurements obtained toward the AF and the results of Heiles (1994) is difficult as the latter were derived from measurements in the −4.70 < log I([CII]) < −3.74 range, much lower than the overall integrated intensities in the AF. The scatter in Fig. 5 can be largely explained by abundance variations in the [C/N] ratio from our model, which mimic the scatter in the I([CII])i versus I([NII])i relationship reasonably well. Only a small amount of scatter can be attributed to variations in the hardness of the radiation field due to changes in the ionization parameter U, as the I([CII])i/I([NII])i ratio is largely insensitive to it over the parameter space covered by our calculations (Abel 2006). The models most consistent with the general behavior of the measurements have metallicities between 1 Z and 2 Z, a volume density of log n(H) = 3.5, and an ionization parameter of log U = −1 to −2. These values are fully consistent with average conditions of the gas in the GC (Mezger et al. 1996; Güsten & Philipp 2004; Rodríguez-Fernández, & Martín-Pintado 2005). The vast majority of data points are well reproduced by models that satisfy [C/N]AF < [C/N]ON. Given the large carbon elemental abundances found toward the GC ranging from 3 to 10 times the solar value (Arimoto et al. 1996; Sodroski et al. 1995; Oka et al. 2005), the model results generally suggest the presence of a nitrogen overabundance in the GC with regard to the value in the Galactic disk (Simón-Díaz & Stasińska 2011) that keeps [C/N]AF < [C/N]ON in the AF region. Moreover, these results suggest that [N/H] does not scale in the same way as [C/H] when moving from the Galactic disk toward the GC, but that it must rise faster in order to keep [C/N]AF < [C/N]ON. If this is indeed the case, any combination of physical processes responsible for this would have to satisfy at least these two conditions: (1) to enrich the GC ISM with nitrogen more efficiently than with carbon, and (2) to transport this nitrogen-enriched gas toward the inner parts of the GC, where the AF are thought to be located.

For the first condition, it would be necessary for carbon and nitrogen production to be decoupled. Garnett et al. (1995) found that, in irregular galaxies, the trend of the [C/N] abundance ratio increases with [O/H] abundance, while for solar abundance stars and H II regions, it abruptly decreases, suggesting that the bulk production of nitrogen might indeed be decoupled from the carbon production in those galaxies. In terms of nitrogen formation in the CNO cycle, the precursor agent for 14 N is the 13 C carbon isotope via the 13C + H → 14 N + γ reaction. The relative amount of 13C in the GC is about three times that available in the solar neighborhood with [12 C/13C] ~ 25 (Güsten & Philipp 2004), so stars forming from non-primordial gas in the GC will already contain a larger amount of 13 C available forthe formation of nitrogen than in other parts of the Galaxy. In this case, the secondary production of nitrogen from low- to intermediate-mass stars in the GC bulge could be a plausible origin of such an elemental abundance increase. Moreover, secondary nitrogen production has been reported for other Galactic nuclei similar to our GC (Liang et al. 2006; Pérez-Montero et al. 2011; Florido et al. 2015), and in our own Galaxy the synthesis of nitrogen via secondary production has been found to be a very important contributor to the present-day [N/H], especially for H II regions in the inner 5 kpc of the Galactic disk (Martín-Hernández et al. 2002). Concerning the transport of nitrogen-enriched gas, the mass supply mechanism to the GC has to be investigated. Galactic bar models suggest that ejected gas from stars in the bulge can replenish the gas in the outer X1 orbits within 1 to 2 orbital periods of the stellar bar rotation (Jenkins & Binney 1994; Mezger et al. 1996; Martel et al. 2018). This gas would then migrate to the inner X2 orbits via loss of angular momentum, mixing with the gas already in these orbits, keeping [N/H] high with regard to the Galactic disk value. As a result, a large fraction of the metals in the GC (carbon, nitrogen, and oxygen mainly) would not have originated within the GC. They would have been produced along the entire length of the Galactic bulge (~few kpc) and then transported to the GC by gas flows (Martel et al. 2013). If this scenario is correct, we expect elemental abundances in the GC to be consistent with values found at somewhat larger galactocentric radii (RGAL). Given that secondary nitrogen production in galactic nuclei of barred galaxies with physical conditions similar to the Milky Way’s GC has been claimed to be responsible for the enhancement of the [N/H] abundance, we suggest that the scatter in Fig. 5 is consistent with the same [N/H] enhancement mechanism in the Milky Way’s GC. If the gas supply mechanism for the GC is mainly due to the mass loss from low- to intermediate-mass stars in the Galactic bulge, this would not only impact the [N/H] abundance via secondary nitrogen production, but it would alter the elemental abundances of other species such as the 17 O and 18 O oxygen isotopes as well, lowering the [18 O/17O] ratio with regard to the Galactic value ≿4 (Wouterloot et al. 2008). Consistently, lower values around [18 O/17O] ~ 3 have been reported in a few places within the GC volume (Penzias 1981; Zhang et al. 2015). Moreover, our ongoing work on the [18 O/17O] abundance ratio in the AF region (Steinke et al., in prep.) has preliminary results consistent with the values previously reported for the GC.

thumbnail Fig. 4

Scatter plot of I([CII]) versus I([NII]) integrated intensities calculated for increasingly larger radial velocity ranges from −40 to −25 km s−1, −50 to −15 km s−1, −60 to −5 km s−1, −70 to +5 km s−1, −80 to +15 km s−1, and −90 to +25 km s−1 (corresponding to I([CII])i versus I([NII])i). The I([CII]) versus I([NII]) relationship in Abel (2006) for metallicities 1.0 Z and 0.1 Z are shown as dashed and dotted straight lines in each panel. As the radial velocity range increases, the bulk of the data points move along the Abel (2006) prediction scale for the [CII] emission contribution of H II regions.

thumbnail Fig. 5

Scatter plot of the I([CII])i versus I([NII])i integrated intensities. Data points are shown as open black circles. The Abel (2006) predictions for the [CII] emission contribution from H II regions, for metallicities 1.0 Z and 0.1 Z (Eqs. (1) and (2) in their work), and the Heiles (1994) linear correlation between I([CII]) and I([NII]) for the extended low-density warm ionized medium (ELDWIM; Eq. (12b) in their work) are shown as dashed, dotted, and dotted-dashed lines, respectively. The model predictions in this work, using different [C/N] abundance ratios normalized to the Galactic disk value in Simón-Díaz & Stasińska (2011) are displayed in colors. Metallicities 1.0 Z and 2.0 Z are shown asfilled and open ■ and ▴ symbols for log n(H) = 3.5 and 2.5, respectively. The symbol size is proportional to the ionization parameter log10 U = −1, −2,−3, and −4.

5 [CII] emission from H II regions

The [CII] and [NII] emission distributions in the AF are spatially and spectrally very closely related (García et al. 2016). Therefore, given their different physical origin, it is expected that part of the observed [CII] emission does not originate from PDRs, but instead from the H II region component. In order to disentangle the latter from the total observed [CII] emission, we re-calibrated the I([CII]) versus I([NII]) relationship in H II regions (Abel 2006) for the physical conditions of the gas in the AF and applied this to the radial velocity components listed in Table 2. This is shown in Fig. 6. In the figure, the model data points satisfy 0.3 < [C/N]AF/[C/N]ON < 0.5, while measurement points represent I([CII])g versus I([NII])s integrated intensities. A least-square fit to the model data points yields logI(CIIHII)=1.068×logI([NII])+0.645.\begin{equation*}\log I(\CII_{\textrm{HII}}) = 1.068\,{\times}\, \log I({[NII]}) &#x002B; 0.645. \end{equation*}(2)

Equation (2) allows us to predict the [CII] emission contribution from the H II region I([CII]HII) in the AF based on our [NII] line observations for each of the identified radial velocity components. Then, by subtracting [CII]HII from the total [CII] emission, the amount of [CII] emission originating only in the PDR component I([CII]PDR), was derived. Following our integrated intensity definitions in Sect. 3, I([NII])s and I([NII])i were used in Eq. (2), yielding the I([CII]HII) values Is and Ii listed in Table 3, respectively. Then, I([CII]PDR) was obtained by subtracting the I([CII]HII) estimations from the I([CII])g, I([CII])s, and I([CII])i integrated intensities. These results are represented by the values Dg,s = I([CII])gIs, Ds,s = I([CII])sIs, and Di,i = I([CII])iIi in same table, with subindexes g, s, and i referring to the integrated intensities defined in Sect. 3. Finally, the relative [CII] emission contribution from PDRs to the total observed emission contained in I([CII])g, I([CII])s, and I([CII])i, is estimated as Rgs,g = Dg,s/I([CII])g, Rss,s = Ds,s /I([CII])s, and Rii,i = Di,i /I([CII])i, respectively.

To obtain the results in Table 3, two main assumptions were made. First, we assumed that all the features present in the observed [CII] spectra were the result of the superposition of emission and absorption features arising fromforeground sources in the Galactic disk (e.g., spiral arms whose emission contribution was properly removed by fitting Gaussian profiles to the observed spectra), gas associated with an H II region component, and gas associated with a PDR component. We did not include the possibility of [CII] emission arising from diffuse CO−dark H2 clouds, which has been claimed to be the case for Galactic disk sources (Langer et al. 2014, 2017). This is justified by the large average volume densities (~104 cm−3) in the GC (Güsten & Philipp 2004), which make such a contribution small, if any at all. Second, the results shown in Table 3 rely entirely on Eq. (2) being a reasonable model to account for the [CII] emission contribution from H II regions to the observed [CII] emission in the AF. Models satisfying 0.3 < [C/N]AF/[C/N]ON < 0.5 can explain the I([CII])i versus I([NII])i scatter relationship in Fig. 5, in the same way shown in Fig. 6 of Abel (2006) for the Orion Nebula. For [C/N]ON = 2.81 (Simón-Díaz & Stasińska 2011), our results imply that the [C/N]AF abundance ratio has to satisfy 0.84<[C/N]AF<1.41.\begin{equation*}0.84 < \textrm{[C/N]}_{\textrm{AF}} < 1.41. \end{equation*}(3)

We defer the discussion on the feasibility for the gas in the AF to have [C/N] abundance ratios in the range defined by the condition in Eq. (3) to Sect. 6.

Calculating integrated intensities for both species in three different ways (Ig, Is, and Ii) allows us to check the consistency between different [CII]HII and [CII]PDR estimates, since the complexity of the line profiles shown in Appendix A could easily yield misleading results. The relative [CII] emission contributions from PDRs to the observed total emission Rgs,g, Rss,s, and Rii,i are remarkably similar. This is particularly true for Rgs,g and Rss,s, which are fully consistent within the error uncertainties. In the case of W1-24km s−1, the derived Rii,i value is much lower than both Rgs,g and Rss,s, outside error estimates. From the W1-24 km s−1 spectra in Appendix A, it seems that [NII] emission unrelated to the source might be included in the −90 km s−1 to +25 km s−1 radial velocity range, artificially increasing the I([CII]HII) prediction via Eq. (2), decreasing the PDR contribution for this component, and consequently yielding an artificially low Rii,i value. The Rgs,g values in Table 3 are the closest to the true fraction of [CII] emission originating from PDRs in the AF as they are derived from the I([CII])g and I([NII])s integrated intensities, which are the most accurate approximation of the true emission of the radial velocity components in both spectral lines. The fraction of [CII] emission in PDRs varies significantly among the different radial velocity components, from ~0.20 up to ~0.75, with no apparent spatial trend with regard to the position of the Arches Cluster, shown in Fig. 1. Our results are in agreement with the predictions in Abel (2006) showing that between 10% and 60% of the observed [CII] emission can be expected to arise from H II regions in the Orion Nebula. Langer et al. (2017) used [CII] line observations from the GOT C+ Survey (Langer et al. 2010; Pineda et al. 2013; Langer et al. 2014) and [NII] line observations from a Herschel-HIFI survey of ten GOT C+ lines (Goldsmith et al. 2015; Langer et al. 2016) to disentangle the [CII] emission contribution from the fully ionized and neutral (defined as the combination of PDR and CO−dark H2) gas, toward the LOS (l, b) = (0,0), and within the GC volume. They reported that the fraction of [CII] emission arising from the fully ionized and neutral gas toward this LOS is ~0.30 and ~0.70 of the total observed [CII] emission, respectively. Our results in Table 3 show that at least 5 out of the 12 radial velocity components are fully consistent with their results, within error uncertainties, despite the fact that they are found toward different LOS. Since Langer et al. (2017) showed that the [CII] emission in the central molecular zone (CMZ) arises primarily from PDRs and highly ionized gas, we have excluded the contribution of any CO−dark H2 gas in the AF from the comparison. For the other components in Table 3, the fraction of [CII] emission found toward H II regions and PDRs changes substantially. This reinforces the idea that extrapolating these results to different LOS, or even to other locations within the GC volume, has to be taken with caution. Moreover, such large variations could have a strong impact when trying to model the [CII] emission from the AF, as [CII] emission is a standard probe of the FUV fields in PDRs, dominating the cooling of atomic gas at low FUV fields (G° < 103) and low densities (ncrit([CII]) ~3 × 103 cm−3) (Hollenbach, & Tielens 1999). Hence, any PDR modeling attempt that does not properly disentangle the [CII] emission contributions from H II regions and PDRs from the total [CII] emission could give misleading results.

thumbnail Fig. 6

Scatter plot of I([CII])g versus I([NII])s integrated intensities for radial velocity components listed in Table 2. Data points are shown as filled black circles, while model points satisfying 0.3 < [C/N]AF/[C/N]ON < 0.5 are shown in gray, with symbol shape and size depending on the gas physical parameters, as listed in the upper left corner of the figure. A least-squares fit to the model predictions yields log I([CII]HII) = 1.068 × log I([NII]) + 0.645, shown by the gray straight line. For comparison, the expected [CII] emission contribution from H II regions in the Abel (2006) models for 1.0 Z and 0.1 Z are shown as dashed and dotted black straight lines, respectively.

6 [C/N] elemental abundance ratio in the AF

The results in Sect. 5 suggest that a lower [C/N] abundance ratio than that found toward the Galactic disk (Simón-Díaz & Stasińska 2011) is needed to explain the I([CII])i versus I([NII])i observationsin Fig. 5 seen toward the AF in the GC. This can be tested in a relatively independent way by combining our [NII] line observations and the I([CII]) integrated intensity model predictions in Eq. (2) with estimates for the electron density (ne) and electron temperature(Te) in each radial velocity component. Together, these observational/empirical constraints allow us to derive the [C/N] abundance ratio directly, under a few reasonable assumptions, so the consistency of the [C/N] abundance ratio range in Eq. (3) estimated for the AF region can be tested.

Following Röllig et al. (2016), under the assumption of local thermodynamic equilibrium (LTE) and optically thin [NII] emission (Goldsmith et al. 2015), the ionized nitrogen column density N(N+) for any given radial velocity component in Table 4 is derived as N(N+)=I([NII])c(ne,Te)[cm2],\begin{equation*}N(N^{&#x002B;}) = \frac{I({[NII]})}{c(n_{\textrm{e}},T_{\textrm{e}})} [\textrm{cm}^{-2}], \end{equation*}(4)

where I([NII]) is the integrated intensity I([NII])s or I([NII])i (in units of K km s−1) listed in Table 2. The c(ne,Te) coefficients are scaling factors between I([NII]) and N(N+) as a function of the adopted ne and Te, that take into account the relative population of the three fine-structure levels of the nitrogen ion (3 P0, 3 P1, and 3 P2) due to collisions with electrons and spontaneous decay. They are calculated for the [NII] 205 and 122 μm transitions in Röllig et al. (2016). These coefficients are almost constant for 5000 K < Te < 8000 K and vary around 15% for 100 cm−3 < ne < 500 cm−3 (Röllig et al. 2016), which are the relevant ranges for the gas physical conditions in the AF region and for the observed [NII] line transition at 205 μm used in this work. The optically thin assumption for the [NII] line emission is well justified as Goldsmith et al. (2015) showed that, for regions with ne ≥ 100 cm−3, column densities up to N(N+) ~ 1018 cm−2 still yield optically thin emission. This is in agreement with our nitrogen column density results in Table 4, which are larger than the average in the Galactic plane N(N+) ~ 4.7 × 1016 cm−2 (Goldsmith et al. 2015), which is expected given the [NII] emission concentration toward the GC.

The ionized carbon column density N(C+) is calculatedfollowing Eq. (1) in Langer et al. (2017), under the assumption of the EOT limit for the [CII] line transition (Goldsmith et al. 2012): N(C+)=2.92×1015[1+eΔE/kTkin2(1+ncr,ene)]I([CII])[cm2],\begin{equation*}N(C^{&#x002B;}) = 2.92\,{\times}\, 10^{15}\Bigg[1 &#x002B; \frac{e^{\Delta E/kT_{\textrm{kin}}}}{2} \Bigg(1&#x002B;\frac{n_{\textrm{cr,e}}}{n_{e}}\Bigg) \Bigg] I({[CII]})[\textrm{cm}^{-2}], \end{equation*}(5)

where I([CII]) is the [CII]HII integrated intensity (in units of K km s−1) listed in Table 3 for each radial velocity component. In Eq. (5), a critical electron density ncr,e ~ 50 cm−3 at kinetic temperatureTkin ~ 8000 K is adopted. These values are typical for H II regions and are close to the average Te ~ 7000 K found toward the AF (Pauls et al. 1976; Goldsmith et al. 2012; Langer et al. 2015). The fact that collision partners other than electrons are not considered in Eq. (5) is based on the assumption that electrons dominate collisions over H and H2 in regions where the [NII] line is detected (Langer et al. 2015). For the [CII] line transition at 158 μm, the term ΔE/k = 91.21 K. The factor 2.92 × 1015 is obtained from a combination of fundamental constants (e.g., h, c, k), the frequency of the [CII] line transition, and the adaptation of Eq. (5) to K km s−1 units for the integrated intensities (Langer et al. 2015). Finally, we approximate Tkin ~ Te in Eq. (5) for each radial velocity component in Table 4.

Under the aforementioned assumptions, Eqs. (4) and (5) are basically a function of I([CII]) and I([NII]) integrated intensities, and the gas physical parameters ne and Te. The latter two were derived by Lang et al. (2001) for the entire extension of the AF using high-resolution 8.3 GHz continuum and H92α recombination line observations. We assigned ne and Te estimates listed in their Tables 4 and 7 to our positions W1, W2, G010, G007, E2N, E2S, and E1 based on spatial coincidence and assuming a single value for each parameter along the LOS. For position E1, we adopted the electron density derived for their source D, which is the spatially closest source to E1 in their data. At the location of P1 and P2, just outside the AF as shown in Fig. 1, no prominent 8.3 GHz continuum emission was detected in Lang et al. (2001). Given that their measured I([NII]) integrated intensities in Table 2are not particularly lower than the ones detected toward the other locations within the AF, the lack of prominent continuum emission might indicate an LOS distance effect. Therefore, we adopted the upper limit ne ~ 100 cm−3 derived for molecular clouds located far from thermal continuum sources in the GC derived by Rodríguez-Fernández, & Martín-Pintado (2005) and the lowest electron temperature found by Lang et al. (2001) toward the AF with Te ~ 5000 K. The c(ne, Te) coefficients for each position were derived by interpolating in ne parameter space the values in Röllig et al. (2016), as c(ne, Te) is constant over the relevant Te range in the present work.

From Eqs. (4) and (5), the [C/N] abundance ratio is obtained by first assuming that the gas is fully single ionized (Langer et al. 2017), and secondly, that the ratio of column densities reflects the ratios of the carbon and nitrogen elemental abundances, similarly to our model assumptions in Sect. 4. Then, it follows: [N(C+)N(N+)]~[CN].\begin{equation*}\Bigg[\frac{ N(C^{&#x002B;})}{N(N^{&#x002B;})} \Bigg] \sim \Bigg[\frac{C}{N} \Bigg]. \end{equation*}(6)

The first assumption in Eq. (6) is analogous to assuming that the carbon and nitrogen ionization correction factors (ICFs), defined as semiempirical correction factors that scale any given ionic elemental abundance to derive the corresponding elemental abundance of the neutral species, and usually used in observations of FIR to optical lines from H II regions (Simpson et al. 1995; Martín-Hernández et al. 2002; Rodríguez-Fernández, & Martín-Pintado 2005), are around unity for the [CII] line transition at 158 μm and [NII] line transition at 205 μm. The second assumption in Eq. (6) implies that the depletion of carbon and nitrogen atoms onto dust grains is low. There is some observational evidence suggesting that, in Nebulae, a significant fraction of carbon atoms might be trapped in carbonaceous grains and in polycyclic aromatic hydrocarbons (PAHs), although it is rather unknown how severe these effects can be (Draine 2003; Jenkins 2009; Simón-Díaz & Stasińska 2011). Given that [C/H] elemental abundances in the GC are rather large in relation to Galactic disk values (Arimoto et al. 1996; Sodroski et al. 1995; Oka et al. 2005; Simón-Díaz & Stasińska 2011), and that there is mounting evidence of the presence of smaller dust grain sizes in the GC compared to the Galactic disk (Nishiyama et al. 2006, 2008, 2009), probably pointing to efficient large-scale dust destruction mechanisms in the GC, it seems reasonable to expect that most of the carbon within the harsh environmental conditions in the GC is found in the gas phase rather than in the solid phase. In the nitrogen case, Jenkins (2009) conducted an extensive study toward several LOS in the Milky Way characterizing the depletion of C, O, and N, among other elements, onto dust grains, finding no evidence for nitrogen depletion. A possible cause for this could be related to the formation of N2, which is very stable, and that has been suggested to inhibit the depletion of nitrogen into any solid compound in the ISM (Gail & Sedlmayr 1986). The only exception for nitrogen to be found in the solid phase could be related to the presence of ices (Hily-Blant et al. 2010; Simón-Díaz & Stasińska 2011), which are not expected to live long within the GC environment containing the ionized gas. To summarize, we consider the assumptions yielding Eq. (6) a reasonable first-order approximation to describe the relationship between ionized and neutral carbon and nitrogen gas in the AF. The values used in Eqs. (4) and (5) to derive the column densities N(N+) and N(C+), and the [C/N] abundance ratios in the AF via Eq. (6), are summarized in Table 4. For the calculations, integrated intensities I([NII])s and I([NII])i in Table 2 and I([CII])HII values in Table 3 were used to derive [C/N]s,s and [C/N]i,i, where subindexes identify the corresponding N+ and C+ column densities involved in the [C/N] abundance ratio calculation. Error uncertainties were obtained by propagating the experimental errors in Is and Ii, through Eqs. (4)–(6). The derived [C/N]s,s and [C/N]i,i abundance ratios are fully consistent within error uncertainties. Therefore, in the following, we focus our discussion only on the [C/N]s,s results.

The [C/N] abundance ratios found toward the AF are surprisingly homogeneous, with all estimates lying within a narrow 0.95 ≤ [C/N] ≤ 1.25 range. This result is in full agreement with the [C/N] range shown in the condition defined in Eq. (3). This shows that using the derived [CII] and observed [NII] emission from the H II region to derive the [C/N] abundance ratios in the AF yields consistent results with two almost independent approaches. The only caveat in the direct comparison between both [C/N] estimates is that the N(C+) and N(N+) column densities in Eqs. (4) and (5) are not fully independent as the I([CII])HII integrated intensities used in Eq. (5) are derived from I([NII]) integrated intensities via Eq. (2). Nonetheless,when considering all the variables involved in the [C/N] calculations, the results in the present work strongly suggest thatthe [C/N] abundance ratio in the AF differs substantially from that measured in the Galactic disk (Simón-Díaz & Stasińska 2011), pointing to a different ISM enrichment process in the GC. Using the average [C/N]AF = 1.13 ± 0.09 in Table 4 and the abundance ratio derived for the Orion Nebula [C/N]ON = 2.81 ± 0.01 (Simón-Díaz & Stasińska 2011), we obtain [C/N]AF/[C/N]ON = 0.40 ± 0.03, in full agreement with what is expected from our model results in Sect. 5.

The results in Table 4 allow us to predict the average nitrogen elemental abundance toward the AF [N/H]AF, as long as the average carbon elemental abundance [C/H]AF in the GC is provided. A simple expression to estimate [N/H]AF as a function of our results and the abundance derived in the Galactic disk for the Orion Nebula (Simón-Díaz & Stasińska 2011) is given by [NH]AF=α×([C/N]ON[C/N]AF)×[NH]ON,  \begin{equation*}\Bigg[\frac{N}{H} \Bigg]_{\textrm{AF}} = \alpha \times \Bigg(\frac{[C/N]_{\textrm{ON}}}{ [C/N]_{\textrm{AF}}} \Bigg) \times \Bigg[\frac{N}{H} \Bigg]_{\textrm{ON}}, \\ \end{equation*}(7)

where α > 1 is a scaling factor between the carbon elemental abundance in the GC and the one found in the Galactic disk so that [C/H]AF = α × [C/H]ON, with α ranging between 3 and 10 in the GC (Arimoto et al. 1996; Sodroski et al. 1995; Oka et al. 2005). In a conservative estimate with α = 3, adopting [C/N]ON = 2.81 (Simón-Díaz & Stasińska 2011), [C/N]AF = 1.13 from Table 4, and [N/H]ON = 8.32 × 10−5 for the Orion Nebula (Simón-Díaz & Stasińska 2011), Eq. (7) yields an average nitrogen elemental abundance in the AF of [N/H]AF = 6.21 × 10−4. This prediction can be compared with the values derived by Simpson et al. (1995) for G0.10+0.02 ([N/H]S95 = 2.66 ± 0.64 × 10−4) and G1.13−0.11 ([N/H]S95 = 4.39 ± 1.00 × 10−4) also located in the GC. The average [N/H]AF value is a factor of 2.3 and 1.4 greater than those derived by Simpson et al. (1995), respectively. Considering that their results rely heavily on a rather uncertain ICF(N+2) value for their [NIII] observations at 57 μm, such a discrepancy is not surprising. We also notice that the same factor of 2.3 is obtained for G0.10+0.02 (G010-17 km s−1 in Table 4) if its individual [C/N]s,s and [C/N]i,i values are used to calculate the [C/N] average in the GC instead of using [C/N]AF in Eq. (7).

The production of nitrogen and its evolution in Galaxies is a highly nontrivial process compared to other species such as oxygen and carbon. Nitrogen can be produced by high-, intermediate-, and low-mass stars, enriching the ISM on different timescales, while oxygen is produced mostly by massive, short-lived stars. Hence, the [N/O] abundance ratio gives important information on the primary or secondary nature of nitrogen. Esteban & García-Rojas (2018) derived the radial distribution of the 12 + log(N/H) elemental abundance and the [N/O] abundance ratio as a function of RGAL down to 5.7 kpc. In their Fig. 4, a steadily increasing [N/H] with decreasing RGAL can be seen, yielding a least−squares fit 12 + log(N/H) = 8.21 − 0.059 × RGAL. Very recently, Arellano-Córdova et al. (2020) reassessed these results taking advantage of the large improvement in the accuracy of heliocentric distances for H II regions in the Galactic disk obtained from Gaia parallaxes of the second data release (DR2; Gaia Collaboration 2018), confirming the general trend of the nitrogen elemental abundance radial distribution found by Esteban & García-Rojas (2018). Nonetheless, whether the [N/H] radial distribution holds for smaller RGAL is unclear. An extrapolation of the Esteban & García-Rojas (2018) results toward the GC yields [N/H] ~ 1.6 × 10−4 which is a factor of ~4 lower than our [N/H]AF prediction. Esteban & García-Rojas (2018) and Arellano-Córdova et al. (2020) also showed that the [N/O] ratio is almost flat for RGAL > 8 kpc, suggesting that the bulk of nitrogen in the Galaxy is not formed by standard secondary processes above this RGAL, but for the H II region S2-61 located at RGAL ~ 5.7 kpc, the largest [N/O] ratio in their sample is found, showing a large enhancement in the inner parts of the Galactic disk, and suggesting that the secondary production of nitrogen plays an important role for small RGAL. Stanghellini & Haywood (2018) showed than nitrogen can be highly enhanced at the expense of carbon in high-mass stars, suggesting that a nitrogen-rich shell is ejected by hot-bottom burning (HBB) stars, transforming most of their carbon into nitrogen. Colzi et al. (2018) studied the fractionation of nitrogen in the Galaxy by measuring the isotopic [14 N/15N] abundance ratio as a function of RGAL. Both isotopes have different origins, with massive stars in their He-burning phase as the primary contributors of 14 N to the ISM, while 15N is produced mostly in the hot CNO cycle of novae outbursts by low-mass stars. They reported a decrease in the [14 N/15N] ratio toward the GC, sampled down to RGAL ~ 2 kpc, implying that the secondary production of nitrogen dominates over its primary production.

Our results suggest that the nitrogen elemental abundance in the GC might be underestimated, and that the nitrogen enrichment of the ISM increases much faster toward smaller RGAL than the carbon elemental abundance. A plausible explanation for this behavior is the secondary production of nitrogen from low- to intermediate-mass stars in the Galactic bulge, whose mass loss would act as an important mass supply source to the inner parts of the GC where the AF are located.

Table 3

[CII] emission contribution from H II regions ([CII]HII) and PDRs ([CII]PDR) in the Arched Filaments (in units of ergs s−1 cm−2) for the radialvelocity components listed in Table 2.

Table 4

[C/N] abundance ratios derived from [CII] and [NII] line observations of the AF.

7 Summary and conclusions

In this work, we studied the relative contributions of H II regions and PDRs to the total [CII] emission observed toward the AF and nearby positions, and their implications for the ISM composition in the GC. Our analysis is based on [CII] and [NII] line observations at 158 μm and 205 μm, respectively, published in García et al. (2016) (publicly available), complemented with observations of the optically thin C18O(1–0) transition to be published in detail elsewhere. The [CII] and [NII] data sets in García et al. (2016) represent the most spatially extended velocity-resolved observations available for the entire Sgr A complex (containing the AF) at the time of publication. Our main results can be summarized as follows:

  • The spatial distribution of the I([CII])/I([NII]) integrated intensity ratios (in units of ergs s−1 cm−2) toward the AF is relatively homogeneous. The bulk of the ratios is confined between two prominent lanes of high I([CII])/I([NII]) values, with one located south of the radio arc and the other found to the north of the Sgr A-east SNR, outlined by 20 cm radio continuum observations (Yusef-Zadeh & Morris 1987). Some striking features in the I([CII])/I([NII]) spatial distribution are: (1) a strong decrease in the high ratio lanes at the intersections (at least in projection) with W1 (south of the radio arc) and the location of the H3 H II region (north of Sgr A-east), while no decrease is seen at the intersection with the E1, E2S, and E2N (south of the radio arc); and (2) a local enhancement of the I([CII])/I([NII]) ratio at the location of G0.07+0.04, where it has been suggested that the gas is interacting with the NTF, tracing the local magnetic field in the GC (Lang et al. 1999b). The features in the I([CII])/I([NII]) spatial distribution strongly suggest that geometry effects in the form of the distance component along the LOS must play a role in the coupling between the thermal plasma responsible for the radio continuum emission and the ionized gas responsible for the observed I([CII])/I([NII]) ratios in the AF, while localized physical phenomena might be driving the high ratios close to major nonthermal sources such as the radio arc and the Sgr A-east SNR. A histogram of the I([CII])/I([NII]) ratios in the AF can be well described by a Gaussian profile, with a mean of Rc = 5.1 ± 0.1 and standard deviation Rσ = 2.0 ± 0.1.

  • The I([CII]) versus I([NII]) relationship in the AF region was studied by expanding the models in Abel (2006) for the relative [CII] emission contribution from H II regions to the total observed [CII] emission toward PDRs. Different values of the [C/N] abundance ratio, SEDs consistent with stellar clusters, and a parameter grid for the hydrogen density n(H), metallicity Z, and ionization parameter U, suitable for the physical conditions of the ISM at the GC, were included. Our models show that the large scatter in the I([CII]) versus I([NII]) relationship can be largely explained by variations in the [C/N] abundance ratio in the AF. The models satisfying 0.3 < [C/N]AF/[C/N]ON < 0.5, with the abundance ratio [C/N]AF predicted for the AF and [C/N]ON derived for the Orion Nebula in the Galactic disk (Simón-Díaz & Stasińska 2011) best disentangle the [CII] emission contribution from H II regions and PDRs toward the AF in the GC, similarly to the predictions in Abel (2006) for the Orion Nebula. We interpret the fact that carbon-to-nitrogen abundance ratios in the AF and in the Galactic disk satisfy [C/N]AF/[C/N]ON < 1 as a strong indication of a large nitrogen enrichment of the ISM when moving from the Galactic disk to the GC. Carbon elemental abundances have already been estimated at between 3 and 10 times larger in the GC than the values measured in the Galactic disk (Arimoto et al. 1996; Sodroski et al. 1995; Oka et al. 2005).

  • We derive the [CII] emission contribution from H II regions traced by [NII] line observations toward nine selected positions in and around the AF: seven positions trace the filaments’ segments E1, E2N, E2S, G0.10+0.02, W1, W2, and G0.07+0.04, and two positions, P1 and P2, trace nearby high-density gas outside them (Serabyn & Güsten 1987). By combining the optically thin [CI](1–0) and C18O(1–0) transitions, we identify a total of twelve radial velocity components within the nine selected positions in [CII] for which we derive the H II region emission contribution. For this purpose, we select the models satisfying 0.3 < [C/N]AF/[C/N]ON < 0.5, with model parameter ranges consistent with the ISM characteristics in the GC, that is, metallicities from 1 Z to 2 Z, hydrogen volume density log n(H) = 3.5, and ionization parameters log U = −1 to −2. A least-squares fit to the model data points yields log I([CII]) = 1.068 × log I([NII]) + 0.645. Then, assuming that the total observed [CII] emission toward the AF originates either in the H II region or PDR only, we show that the fraction of [CII] emission in PDRs can vary dramatically between locations in the AF, with a fraction ranging from ~0.20 to ~0.75. Therefore, any [CII] modeling attempt in the region has to properly disentangle the contribution from H II regions and PDRs in order to avoid the production of misleading results.

  • The results for the carbon-to-nitrogen abundance ratio in the AF are fully consistent with an (almost) independent derivation of the same ratio by calculating N(C+) and N(N+) columndensities, and assuming that they directly mirror the values of [C/N]. For the twelve radial velocity components identified in this work, we obtain abundance ratios in the range 0.95 ≤ [C/N]AF ≤ 1.25 with an average value of [C/N]AF = 1.13 ± 0.09 for the AF, in full agreement with the range previously derived by our models. By combining this result with the carbon-to-nitrogen abundance ratio measured in the Galactic disk [C/N]ON ~ 2.81 (Simón-Díaz & Stasińska 2011) and adopting a three-times-higher carbon elemental abundance in the GC than in the solar vicinity, we predict an average nitrogen elemental abundance for the AF of [N/H]AF = 6.21 × 10−4, which is a factor of 2.3 higher than the one derived by Simpson et al. (1995) for G0.10+0.02, relying on a rather poorly constrained ionization correction factor ICF(N+2) for their [NIII] observations at 57 μm. Nonetheless, given the vastly different approaches between the present work and the one in Simpson et al. (1995), this relatively small discrepancy is remarkable.

  • We propose that our [N/H]AF prediction is close to the real average value in the GC. If this is indeed the case, it would imply that the nitrogen elemental abundance in the GC is a factor of ~7.5 larger than that of the Galactic disk. We suggest that a plausible explanation for such a scenario could be the nitrogen enrichment of the ISM in the GC via the secondary production of nitrogen from low- to intermediate-mass stars in the Galactic bulge. If their mass loss through stellar winds acts as the main gas supply source for the gas in the inner orbits of the GC, a large nitrogen elemental abundance in the AFlocated within the GC is expected. Overall, our results show that tight constraints on the [C/N] abundance ratio for the GC are needed, significantly tighter than previous abundance measurements have discerned.

Our results have shown that a large fraction of the observed [CII] 158 μm emission in the AF region does not originate within PDRs but in H II regions. This can have a considerable impact on PDR model results used to derived the gas physical conditions such as impinging FUV fields, hydrogen densities, hydrogen ionization rates due to cosmic-rays, and so on, because the gas cooling rate traced by this ionized carbon fine-structure line can be largely overestimated and give misleading results. Moreover, these results show that the [C/N] abundance ratio in the AF is considerably smaller than the Galactic plane average value pointing to a nitrogen-enriched ISM in the GC and suggesting a secondary origin of the nitrogen production. The impact of these results on the PDR modeling of the AF region will be investigated in a follow-up publication.

Acknowledgements

This work was sponsored (in part) by the Chinese Academy of Sciences (CAS), through a grant to the CAS South America Center for Astronomy (CASSACA) in Santiago, Chile, by the “Comisión Nacional de Ciencia y Tecnología (CONICYT)” via Project FONDECYT de Iniciación 11170551, and by the Collaborative Research Centre 956, sub-project A5, funded by the Deutsche Forschungsgemeinschaft (DFG). Robert Simon acknowledges support by the french Agence National de Recherche (ANR) and the german Deutsche Forschungsgemeinschaft (DFG) through the project “GENESIS” (project numbers ANR-16-CE92-0035-01/DFG1591/2-1). The data used in this work are based partially on observations carried out under project ID 020-16 with the IRAM 30m telescope. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain).

Appendix A Estimation of I([CII]) and I([NII]) integrated intensities

In this appendix, the procedure to obtain I([CII]) and I([NII]) integrated intensities for each position in Table 1 is outlined. The lack of detection of the optically thin [13 CII] transition, the low brightness of both [CII] and [NII] lines, and the poor S/N in the [NII] data force us to use complementary observations acting as templatesto identify the radial velocity components in both transitions, later used for the data analysis. Many of the existing [13 CII] observationsin massive star-forming regions have shown that its line profile is closely followed by the optically thin C18O(1–0) transition (Guevara et al. 2020). Goldsmith et al. (2012) showed that the [CII] line falls into the so-called EOT limit for peak antenna temperatures ≤8 K, where opacity effects are not the dominant factor shaping the line profile. This can also be seen in the comparison between the [CII], [CI](1–0), and C18O(1–0) observationsin Fig. 2 and in the [CI](2–1) spectra shown in the figures of this appendix. Despite the very different gas physical conditions the aforementioned transitions trace, they show relatively similar broad spectral features (assumed to originate in the GC). Therefore, though opacity effects may be present, they are unlikely to have a significant effect on shaping the observed line profiles. Instead, they are more likely dominated by the presence of different radial velocity components in the spectra. In this sense, by carrying out a multiwavelength comparison, we can estimate the corresponding [CII] emission contribution to all radial velocity components by fitting multi-Gaussian profiles (assumed for simplicity) in the proper radial velocity ranges, where the components are found. Then, for the [NII] line observations, the corresponding integrated intensities are calculated by simply adding the emission, multiplied by the channel spectral width, in the radial velocity range where the [CII] emission for each component is found.

thumbnail Fig. A.1

Graphical description of the approach adopted to derived Ig, Is, and Ii integrated intensities for the [CII] and [NII] line observations, as described in Sect. 3. In the example, the spectra in the upper, middle, and lower panels correspond to the [CII] line at the W2 position. The spectra were normalized to the peak antenna temperature for display purposes.

Figures A.2A.6 show the [CII] (open blue), [NII] (filled black), [CI](2–1) (filled dark gray), [CI](1–0) (filled gray), and C18O(1–0) (open red) color-coded spectra. Antenna temperatures have been scaled by the factor next to the line’s name for display purposes. The red velocity components in the multi-Gaussian fits to the C18O(1–0) lines are used to recover their counterpart radial velocity components in the [CII] spectra multi-Gaussian fits, shown in blue. Dashed-line curves are the radial velocity components identified in both [CII] and C18O(1–0) lines that are used for the analysis in this work, while dotted-line Gaussian components are needed to reproduce other relevant features in the spectra besides the GC-related emission. The vertical dotted lines mark the radial velocity range for which the I([NII]) integrated intensity, associated with each GC [CII] radial velocity component, was calculated. With this approach, five different integrated intensities (listed in Table 2) are calculated. Firstly, we have the Ig values derived from the multi-Gaussian fits to the [CII] spectra. Then we have the Is values for both [CII] and [NII] line transitions obtained by directly adding the emission (multiplied by the spectral channel width) in the radial velocity channels marked by the vertical dotted black lines. Finally, we have the Ii values obtained in the same way as Is, but for the radial velocity range of −90 km s−1 to +25 km s−1, which is used in the data analysis described in Sect. 4. The approach to derive Ig, Is, Ii is graphically shown in Fig. A.1. Given the complexity of the [CII] spectra, there are several emission features that are not used in our analysis as they do not satisfy our requirements to be classified as a radial velocity component of the AF in the GC. Narrow radial velocity components obtained from the Gaussian fits are excluded from the analysis, as they most likely originate outside the GC, given their small line width and central radial velocities coincident with those of spiral arms along the LOS. Broad radial velocity components with positive central radial velocities (VLSR > 0 km s−1) might be part of the GC but are most likely not associated with the AF because most of their gas is found toward negative radial velocities, which is shown in their I([CII]) and I([NII]) integrated intensity channel maps (García et al. 2016). [CII] and [NII] emission features at very negative radial velocities are seen at positions P1, W1, W2, G0.10+0.02, and G0.07+0.04. Since they have neither [CI](1–0) nor C18O(1–0) evident counterparts detected above the rms noise level of the spectra, they are excluded from our analysis as no assessment can be made about opacity effects, origin (inside or outside the GC), or the number of radial velocity components that describe them. In the following, individual notes on the multi-Gaussian fitting procedure are given for each position in Table 1:

  • P1: The C18O(1–0) multi-Gaussian fit identifies two broad components tracing the AF. They are blended with two narrow features assumed to be foreground components. The [CI](1–0) spectrum shows all features, while the [CI](2–1) transition only shows the broad components, reflecting different excitation conditions between the foreground and the GC gas. The multi-Gaussian [CII] fit is forced to reproduce the two broad components, which are blueshifted with regard to the radial velocity fit in the C18O(1–0) spectrum. The narrow foreground feature at −35 km s−1 is detected at the same radial velocity as the one seen in the optically thin C18O(1–0) transition. The [NII] spectrum has most of the emission toward the brightest [CII] component.

  • P2: The C18O(1–0) multi-Gaussian fit identifies one broad component tracing the AF, blended with four narrow features assumed to be foreground components. The [CI](1–0) spectrum shows all the complex features in the C18O(1–0) spectrum, while the [CI](2–1) transition shows mostly the broad component, reflecting different excitation conditions between the foreground and the GC gas. The multi-Gaussian [CII] fit is forced to reproduce the broad component at the center of the line, which is slightly blueshifted with regard to the radial velocity fit in the C18O(1–0) spectrum. The emission in the [NII] spectrum is distributed between −60 km s−1 and +30 km s−1 in nearly three major components, roughly correlated with the emission in the [CII] spectrum.

  • W1: The C18O(1–0) multi-Gaussian fit identifies one broad component tracing the AF, blended with four narrow features assumed to be foreground components. The [CI](1–0) shows all features in the C18O(1–0) spectrum, while the [CI](2–1) transition shows only the broad component and the two narrow features at −6 km s−1 and +2 km s−1, reflecting the changing excitation conditions of the gas along the same LOS. The multi-Gaussian [CII] fit is forced to reproduce the broad component, which is slightly redshifted in relation to the radial velocity fit in the C18O(1–0) spectrum, as well as several less bright features. The [NII] spectrum has most of the emission in the same spectral range as the main [CII] component.

  • W2: The C18O(1–0) multi-Gaussian fit identifies one broad component tracing the AF, blended with two narrow features assumed to be foreground components, and two other narrow features around radial velocity 0 km s−1, probably associated with cold gas along the LOS to the GC. The [CI](1–0) and [CI](2–1) transitions mainly show the broad component and the prominent feature near 0 km s−1, which are seen in the C18O(1–0) spectrum. The multi-Gaussian [CII] fit is forced to reproduce the broad component at the center of the line, which is slightly redshifted with regard to the radial velocity fit in the C18O(1–0) spectrum. The [NII] spectrum has most of the emission toward the main [CII] component.

  • G0.10+0.02: The C18O(1–0) multi-Gaussian fit identifies one broad component tracing the AF, with four other narrow features in the spectrum assumed to be foreground components. The [CI](1–0) spectrum shows almost all the features of the C18O(1–0) spectrum, while the [CI](2–1) transition mostly shows the broad component and the feature near 0 km s−1, reflecting different excitation conditions of the gas along the LOS. The multi-Gaussian [CII] fit is forced to reproduce the broad component at the center of the line, which is slightly redshifted in relation to the radial velocity fit in the C18O(1–0) spectrum. The emission in the [NII] spectrum is mostly distributed between −75 km s−1 and +20 km s−1, and it is roughly correlated with the emission in the [CII] spectrum.

  • G0.07+0.04: The C18O(1–0) multi-Gaussian fit identifies two broad components tracing the AF. They are blended with one small and two narrow components. The [CI](1–0) and [CI](2–1) transitions show all features seen in the C18O(1–0) spectrum. The multi-Gaussian [CII] fit is forced to reproduce the two broad components, which are blueshifted with regard to the radial velocity fit in the C18O(1–0) spectrum. The narrow foreground feature at −35 km s−1 is detected at the same radial velocity as the one seen in the optically thin C18O(1–0) transition. The [NII] spectrum has most of the emission toward the brightest [CII] component, although a considerable fraction of the emission is found in a radial velocity range with no C18O(1–0) counterpart emission, probably tracing gas deeper into the H II region.

  • E2N: The C18O(1–0) multi-Gaussian fit identifies one broad component tracing the AF, blended with three narrow features assumed to be foreground components. The [CI](1–0) shows all features in the C18O(1–0) spectrum, while the [CI](2–1) transition shows only the broad component, reflecting the changing excitation conditions of the gas along the same LOS. The multi-Gaussian [CII] fit is forced to reproduce the broad component, which is blueshifted in relation to the radial velocity fit in the C18O(1–0) spectrum. The narrow feature in the center of the very broad velocity component is interpreted as a different source, and its emission contribution, though small, is fit away from the main component. The [NII] spectrum has most of the emission in the same spectral range as the main [CII] component.

  • E2S: The C18O(1–0) multi-Gaussian fit identifies one very broad component tracing the AF, with four other narrow features assumed to be foreground components. The [CI](1–0) and the C18O(1–0) transitions mostly trace the broad component identified in the C18O(1–0) spectrum, while the narrow feature at +3 km s−1 is detected only in the [CI](1–0) transition, reflecting different excitation conditions of the gas along the LOS. The [CII] spectrum is extremely complex, with several narrow features at different radial velocities blended with a much broader emission component. Hence, the multi-Gaussian [CII] fit is forced to reproduce what we interpret as an underlying broad component at the center of the line, which is consistent with the one observed in the C18O(1–0) line. This component is redshifted with regard to the radial velocity fit in the C18O(1–0) spectrum. The emission in the [NII] spectrum is mostly distributed between −75 km s−1 and +30 km s−1, and it is roughly correlated with the emission of the underlying [CII] broad component.

  • E1: The C18O(1–0) multi-Gaussian fit identifies two broad components tracing the AF. They are blended with three narrow features assumed to be foreground components. The [CI](1–0) spectrum shows emission at all features except for the narrow component at −35 km s−1, while the [CI](2–1) transition only shows the broad components, reflecting different excitation conditions between the foreground and the GC gas. The multi-Gaussian [CII] fit is forced to reproduce the two broad components, which are slightly blue-shifted with regard to the radial velocity fit in the C18O(1–0) spectrum. The narrow feature at −42 km s−1 is fit to improve the overall fit of the broad components. The [NII] spectrum has most of the emission toward the identified [CII] broad components, with the rest found toward positive radial velocities, where a counterpart in the [CII] line is also detected.

thumbnail Fig. A.2

Comparison of [CII] (open blue), [NII] (filled black), [CI](2–1) (filled dark gray), [CI](1–0) (filled gray), and C18 O(1–0) (open red) spectra for positions P1 and P2 in Table 1. Antenna temperatures have been scaled by the factor shown next to the line’s name for display purposes. Gaussian component fits to the [CII] and C18 O(1–0) spectra are shown in blue and red, respectively. Only dashed-line Gaussian components in the [CII] spectra are used for the analysis in this work as only for them can a counterpart in the C18 O(1–0) optically thin emission be identified. Vertical dotted black lines show the radial velocity limits for which integrated intensities were calculated for the [NII] line by directly adding the observed emission, multiplied by the spectral channel width, to that LSR velocity range.

thumbnail Fig. A.3

Comparison of [CII] (open blue), [NII] (filled black), [CI](2–1) (filled dark gray), [CI](1–0) (filled gray), and C18 O(1–0) (open red) spectra for positions W1 and W2 in Table 1. Antenna temperatures have been scaled by the factor shown next to the line’s name for display purposes. Gaussian component fits to the [CII] and C18 O(1–0) spectra are shown in blue and red, respectively. Only dashed-line Gaussian components in the [CII] spectra are used for the analysis in this work as only for them can a counterpart in the C18 O(1–0) optically thin emission be identified. Vertical dotted black lines show the radial velocity limits for which integrated intensities were calculated for the [NII] line by directly adding the observed emission, multiplied by the spectral channel width, to that LSR velocity range.

thumbnail Fig. A.4

Comparison of [CII] (open blue), [NII] (filled black), [CI](2–1) (filled dark gray), [CI](1–0) (filled gray), and C18 O(1–0) (open red) spectra for positions G0.10+0.02 and G0.07+0.04 in Table 1. Antenna temperatures have been scaled by the factor shown next to the line’s name for display purposes. Gaussian component fits to the [CII] and C18 O(1–0) spectra are shown in blue and red, respectively. Only dashed-line Gaussian components in the [CII] spectra are used for the analysis in this work as only for them can a counterpart in the C18 O(1–0) optically thin emission be identified. Vertical dotted black lines show the radial velocity limits for which integrated intensities were calculated for the [NII] line by directly adding the observed emission, multiplied by the spectral channel width, to that LSR velocity range.

thumbnail Fig. A.5

Comparison of [CII] (open blue), [NII] (filled black), [CI](2–1) (filled dark gray), [CI](1–0) (filled gray), and C18 O(1–0) (open red) spectra for positions E2N and E2S in Table 1. Antenna temperatures have been scaled by the factor shown next to the line’s name for display purposes. Gaussian component fits to the [CII] and C18 O(1–0) spectra are shown in blue and red, respectively. Only dashed-line Gaussian components in the [CII] spectra are used for the analysis in this work as only for them can a counterpart in the C18 O(1–0) optically thin emission be identified. Vertical dotted black lines show the radial velocity limits for which integrated intensities were calculated for the [NII] line by directly adding the observed emission, multiplied by the spectral channel width, to that LSR velocity range.

thumbnail Fig. A.6

Comparison of [CII] (open blue), [NII] (filled black), [CI](2–1) (filled dark gray), [CI](1–0) (filled gray), and C18 O(1–0) (open red)spectra for position E1 in Table 1. Antenna temperatures have been scaled by the factor shown next to the line’s name for display purposes. Gaussian component fits to the [CII] and C18 O(1–0) spectra are shown in blue and red, respectively. Only dashed-line Gaussian components in the [CII] spectra are used for the analysis in this work as only for them can a counterpart in the C18 O(1–0) optically thin emission be identified. Vertical dotted black lines show the radial velocity limits for which integrated intensities were calculated for the [NII] line by directly adding the observed emission, multiplied by the spectral channel width, to that LSR velocity range.

References

  1. Abel, N. P. 2006, MNRAS, 368, 1949 [NASA ADS] [CrossRef] [Google Scholar]
  2. Arellano-Córdova, K. Z., Esteban, C., García-Rojas, J., et al. 2020, MNRAS, 496, 1051 [Google Scholar]
  3. Arimoto, N., Sofue, Y., & Tsujimoto, T. 1996, PASJ, 48, 275 [NASA ADS] [CrossRef] [Google Scholar]
  4. Bally, J., Stark, A. A., Wilson, R. W., & Henkel, C. 1987, ApJS, 65, 13 [NASA ADS] [CrossRef] [Google Scholar]
  5. Bally, J., Stark, A. A., Wilson, R. W., & Henkel, C. 1988, ApJ, 324, 223 [NASA ADS] [CrossRef] [Google Scholar]
  6. Carter, M., Lazareff, B., Maier, D., et al. 2012, A&A, 538, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Clark, J. S., Lohr, M. E., Najarro, F., Dong, H., & Martins, F. 2018, A&A, 617, A65 [CrossRef] [EDP Sciences] [Google Scholar]
  8. Clarkson, W. I., Ghez, A. M., Morris, M. R., et al. 2012, ApJ, 753, 92 [Google Scholar]
  9. Colzi, L., Fontani, F., Rivilla, V. M., et al. 2018, MNRAS, 478, 3693 [NASA ADS] [CrossRef] [Google Scholar]
  10. Cotera, A. S., Erickson, E. F., Colgan, S. W. J., et al. 1996, ApJ, 461, 750 [NASA ADS] [CrossRef] [Google Scholar]
  11. Cox, P., & Laureijs, R. 1989, IAU Symp., 136, 121 [Google Scholar]
  12. Dame, T. M., Elmegreen, B. G., Cohen, R. S., & Thaddeus, P. 1986, ApJ, 305, 892 [NASA ADS] [CrossRef] [Google Scholar]
  13. Dong, H., Wang, Q. D., & Morris, M. R. 2012, MNRAS, 425, 884 [NASA ADS] [CrossRef] [Google Scholar]
  14. Dong, H., Mauerhan, J., Morris, M. R., Wang, Q. D., & Cotera, A. 2015, MNRAS, 446, 842 [NASA ADS] [CrossRef] [Google Scholar]
  15. Draine, B. T. 2011, Physics of the Interstellar and Intergalactic Medium (Princeton: Princeton University Press) [Google Scholar]
  16. Draine, B. T. 2003, ARA&A, 41, 241 [NASA ADS] [CrossRef] [Google Scholar]
  17. Esteban, C. & García-Rojas, J. 2018, MNRAS, 478, 2315 [NASA ADS] [CrossRef] [Google Scholar]
  18. Figer, D. F. 2008, IAU Symp., 250, 247 [Google Scholar]
  19. Figer, D. F., Najarro, F., Gilmore, D., et al. 2002, ApJ, 581, 258 [NASA ADS] [CrossRef] [Google Scholar]
  20. Florido, E., Zurita, A., Perez, I., et al. 2015, A&A, 584, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Gail, H.-P., & Sedlmayr, E. 1986, A&A, 166, 225 [NASA ADS] [Google Scholar]
  23. García, P. 2015, PhD thesis, University of Cologne, Germany [Google Scholar]
  24. García, P., Bronfman, L., Nyman, L.-Å., Dame, T. M., & Luna, A. 2014, ApJS, 212, 2 [NASA ADS] [CrossRef] [Google Scholar]
  25. García, P., Simon, R., Stutzki, J., et al. 2016, A&A, 588, A131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Garnett, D. R., Skillman, E. D., Dufour, R. J., et al. 1995, ApJ, 443, 64 [NASA ADS] [CrossRef] [Google Scholar]
  27. Genzel, R., Stacey, G. J., Harris, A. I., et al. 1990, ApJ, 356, 160 [NASA ADS] [CrossRef] [Google Scholar]
  28. Graf, U. U., Simon, R., Stutzki, J., et al. 2012, A&A, 542, L16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Goldsmith, P. F., Langer, W. D., Pineda, J. L., & Velusamy, T. 2012, ApJS, 203, 13 [Google Scholar]
  30. Goldsmith, P. F., Yıldız, U. A., Langer, W. D., et al. 2015, ApJ, 814, 133 [NASA ADS] [CrossRef] [Google Scholar]
  31. Güsten, R., & Philipp, S. D. 2004, The Dense Interstellar Medium in Galaxies (Berlin: Springer), 253 [Google Scholar]
  32. Guevara, C., Stutzki, J., Ossenkopf-Okada, V., et al. 2020, A&A, 636, A16 [CrossRef] [EDP Sciences] [Google Scholar]
  33. Heiles, C. 1994, ApJ, 436, 720 [NASA ADS] [CrossRef] [Google Scholar]
  34. Hily-Blant, P., Walmsley, M., Pineau Des Forêts, G., et al. 2010, A&A, 513, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Hollenbach, D. J., & Tielens, A. G. G. M. 1999, Rev. Mod. Phys., 71, 173 [Google Scholar]
  36. Jenkins, E. B. 2009, ApJ, 700, 1299 [NASA ADS] [CrossRef] [Google Scholar]
  37. Jenkins, A., & Binney, J. 1994, MNRAS, 270, 703 [Google Scholar]
  38. Jones, P. A., Burton, M. G., Cunningham, M. R., et al. 2012, MNRAS, 419, 2961 [NASA ADS] [CrossRef] [Google Scholar]
  39. Lang, C. C., Morris, M., & Echevarria, L. 1999, ApJ, 526, 727 [NASA ADS] [CrossRef] [Google Scholar]
  40. Lang, C. C., Goss, W. M., & Morris, M. 2001, AJ, 121, 2681 [NASA ADS] [CrossRef] [Google Scholar]
  41. Langer, W. D., Velusamy, T., Pineda, J. L., et al. 2010, A&A, 521, L17 [Google Scholar]
  42. Langer, W. D., Velusamy, T., Pineda, J. L., et al. 2014, A&A, 561, A122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Langer, W. D., Goldsmith, P. F., Pineda, J. L., et al. 2015, A&A, 576, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Langer, W. D., Goldsmith, P. F., & Pineda, J. L. 2016, A&A, 590, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Langer, W. D., Velusamy, T., Morris, M. R., Goldsmith, P. F., & Pineda, J. L. 2017, A&A, 599, A136 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Leitherer, C., Ortiz Otálvaro, P. A., Bresolin, F., et al. 2010, ApJS, 189, 309 [Google Scholar]
  47. Liang, Y. C., Yin, S. Y., Hammer, F., et al. 2006, ApJ, 652, 257 [NASA ADS] [CrossRef] [Google Scholar]
  48. Martel, H., Carles, C., Robichaud, F., Ellison, S. L., & Williamson, D. J. 2018, MNRAS, 477, 5367 [Google Scholar]
  49. Martel, H., Kawata, D., & Ellison, S. L. 2013, MNRAS, 431, 2560 [NASA ADS] [CrossRef] [Google Scholar]
  50. Martín-Hernández, N. L., Peeters, E., Morisset, C., et al. 2002, A&A, 381, 606 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Mezger, P. G., Pankonin, V., Schmid-Burgk, J., Thum, C., & Wink, J. 1979, A&A, 80, L3 [NASA ADS] [Google Scholar]
  52. Mezger, P. G., Duschl, W. J., & Zylka, R. 1996, A&ARv, 7, 289 [NASA ADS] [CrossRef] [Google Scholar]
  53. Müller, M., & Jellema, W., “The HIFI Beam: Release #1, Release Note for Astronomers”, 19 Sept 2014 [Google Scholar]
  54. Nishiyama, S., Nagata, T., Kusakabe, N., et al. 2006, ApJ, 638, 839 [Google Scholar]
  55. Nishiyama, S., Nagata, T., Tamura, M., et al. 2008, ApJ, 680, 1174 [NASA ADS] [CrossRef] [Google Scholar]
  56. Nishiyama, S., Tamura, M., Hatano, H., et al. 2009, ApJ, 696, 1407 [NASA ADS] [CrossRef] [Google Scholar]
  57. Oka, T., Hasegawa, T., Sato, F., Tsuboi, M., & Miyazaki, A. 1998, ApJS, 118, 455 [NASA ADS] [CrossRef] [Google Scholar]
  58. Oka, T., Geballe, T. R., Goto, M., Usuda, T., & McCall, B. J. 2005, ApJ, 632, 882 [NASA ADS] [CrossRef] [Google Scholar]
  59. Ossenkopf, V., Röllig, M., Neufeld, D. A., et al. 2013, A&A, 550, A57. [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Pauls, T., Downes, D., Mezger, P. G., & Churchwell, E. 1976, A&A, 46, 407 [NASA ADS] [Google Scholar]
  61. Pérez-Montero, E., Vílchez, J. M., Cedrés, B., et al. 2011, A&A, 532, A141 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Penzias, A. A. 1981, ApJ, 249, 518 [NASA ADS] [CrossRef] [Google Scholar]
  63. Pilbratt, G. L., Riedinger, J. R., Passvogel, T., et al. 2010, A&A, 518, LL1 [Google Scholar]
  64. Pineda, J. L., Langer, W. D., Velusamy, T., et al. 2013, A&A, 554, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Pineda, J. L., Langer, W. D., & Goldsmith, P. F. 2014, A&A, 570, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Rodríguez-Fernández, N. J., & Martín-Pintado, J. 2005, A&A, 429, 923 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. Roelfsema, P. R., Helmich, F. P., Teyssier, D., et al. 2012, A&A, 537, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Röllig, M., Simon, R., Güsten, R., et al. 2016, A&A, 591, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Röllig, M., Szczerba, R., Ossenkopf, V., & Glück, C. 2013, A&A, 549, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Serabyn, E., & Guesten, R. 1987, A&A, 184, 133 [Google Scholar]
  71. Shields, J. C., & Ferland, G. J. 1994, ApJ, 430, 236 [NASA ADS] [CrossRef] [Google Scholar]
  72. Simón-Díaz, S., & Stasińska, G. 2011, A&A, 526, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  73. Simpson, J. P., Colgan, S. W. J., Rubin, R. H., et al. 1995, ApJ, 444, 721 [NASA ADS] [CrossRef] [Google Scholar]
  74. Simpson, J. P., Colgan, S. W. J., Cotera, A. S., et al. 2007, ApJ, 670, 1115 [NASA ADS] [CrossRef] [Google Scholar]
  75. Sodroski, T. J., Odegard, N., Dwek, E., et al. 1995, ApJ, 452, 262 [NASA ADS] [CrossRef] [Google Scholar]
  76. Stanghellini, L., & Haywood, M. 2018, ApJ, 862, 45 [Google Scholar]
  77. Stolte, A., Ghez, A. M., Morris, M., et al. 2008, ApJ, 675, 1278 [NASA ADS] [CrossRef] [Google Scholar]
  78. Timmermann, R., Genzel, R., Poglitsch, A., et al. 1996, ApJ, 466, 242 [NASA ADS] [CrossRef] [Google Scholar]
  79. Wouterloot, J. G. A., Henkel, C., Brand, J., et al. 2008, A&A, 487, 237 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  80. Yusef-Zadeh, F., & Morris, M. 1987, ApJ, 320, 545 [NASA ADS] [CrossRef] [Google Scholar]
  81. Yusef-Zadeh, F., & Morris, M. 1988, ApJ, 329, 729 [NASA ADS] [CrossRef] [Google Scholar]
  82. Yusef-Zadeh, F., Muno, M., Wardle, M., & Lis, D. C. 2007, ApJ, 656, 847 [NASA ADS] [CrossRef] [Google Scholar]
  83. Zhang, J. S., Sun, L. L., Riquelme, D., et al. 2015, ApJS, 219, 28 [NASA ADS] [CrossRef] [Google Scholar]
  84. Zhao, J.-H., Desai, K., Goss, W. M., & Yusef-Zadeh, F. 1993, ApJ, 418, 235 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Selected positions in the AF region shown in Fig. 1.

Table 2

I([CII]) and I([NII]) integrated intensities in units of ergs s−1 cm−2 for the velocity components identified at the positions in Table 1.

Table 3

[CII] emission contribution from H II regions ([CII]HII) and PDRs ([CII]PDR) in the Arched Filaments (in units of ergs s−1 cm−2) for the radialvelocity components listed in Table 2.

Table 4

[C/N] abundance ratios derived from [CII] and [NII] line observations of the AF.

All Figures

thumbnail Fig. 1

AF region outlined in 20 cm continuum emission from Yusef-Zadeh & Morris (1987). Selected positions for the I([CII]) versus I([NII]) integrated intensity analysis in this work: E1, E2-N, E2-S, W1, W2, G0.07+0.04, G0.10+0.02, P1, and P2, are shown as open squares in the García et al. (2016) data grid. The 46″ spatial resolution of their data is represented by open circles. Filled triangles depict either the closest continuumemission peak to each position or one of the two CS(2–1) peaks (P1 or P2) in the observations by Serabyn & Güsten (1987). Panel insets contain the [CII] (blue), [NII] (red), and scaled C18 O(1–0) (filled-gray) spectra at the positions selected for this work, showing very complex line profiles with large variations in line shape and brightness across the region. Black crosses show the positions of several H II regions labeled H1–H13 spanning a region south of the filaments (the so-called H-Region). The ⋆ symbol shows the position of the Arches Cluster, ☆ while (text here emptystar) symbols indicate the locations of some of the massive (O, B, and WR) field stars scattered across the region (Dong et al. 2012, 2015).

In the text
thumbnail Fig. 2

Comparison between [CII] (filled light gray), [CI](1–0) (gray), and C18O(1–0) (black) line profiles for all positions in Table 1. Relative intensities among the spectra have been scaled for display purposes. The observations trace different phases of the ISM from ionized, neutral, and molecular gas. Narrow and strong emission and absorption features in the C18 O(1–0) optically thin line common to all positions are marked with vertical dotted lines as they are most likely foreground components located in the Galactic disk somewhere along the LOS and outside the GC.

In the text
thumbnail Fig. 3

Bottom: spatial distribution of the I([CII])i/I([NII])i integrated intensity ratio calculated for the radial velocity range −90 km s−1 to +25 km s−1. The AF traced by 20 cm continuum observations and the selected positions in Table 1 are shown as black contours and white circles, respectively. The white ⋆ symbol represents the location of the Arches Cluster. Top: histogram of the spatial distribution of I([CII])i/I([NII])i integrated intensity ratios shown in the map. A Gaussian fit was performed to the histogram to characterize the central ratio (Rc) and standard deviation (Rσ) of the distribution.

In the text
thumbnail Fig. 4

Scatter plot of I([CII]) versus I([NII]) integrated intensities calculated for increasingly larger radial velocity ranges from −40 to −25 km s−1, −50 to −15 km s−1, −60 to −5 km s−1, −70 to +5 km s−1, −80 to +15 km s−1, and −90 to +25 km s−1 (corresponding to I([CII])i versus I([NII])i). The I([CII]) versus I([NII]) relationship in Abel (2006) for metallicities 1.0 Z and 0.1 Z are shown as dashed and dotted straight lines in each panel. As the radial velocity range increases, the bulk of the data points move along the Abel (2006) prediction scale for the [CII] emission contribution of H II regions.

In the text
thumbnail Fig. 5

Scatter plot of the I([CII])i versus I([NII])i integrated intensities. Data points are shown as open black circles. The Abel (2006) predictions for the [CII] emission contribution from H II regions, for metallicities 1.0 Z and 0.1 Z (Eqs. (1) and (2) in their work), and the Heiles (1994) linear correlation between I([CII]) and I([NII]) for the extended low-density warm ionized medium (ELDWIM; Eq. (12b) in their work) are shown as dashed, dotted, and dotted-dashed lines, respectively. The model predictions in this work, using different [C/N] abundance ratios normalized to the Galactic disk value in Simón-Díaz & Stasińska (2011) are displayed in colors. Metallicities 1.0 Z and 2.0 Z are shown asfilled and open ■ and ▴ symbols for log n(H) = 3.5 and 2.5, respectively. The symbol size is proportional to the ionization parameter log10 U = −1, −2,−3, and −4.

In the text
thumbnail Fig. 6

Scatter plot of I([CII])g versus I([NII])s integrated intensities for radial velocity components listed in Table 2. Data points are shown as filled black circles, while model points satisfying 0.3 < [C/N]AF/[C/N]ON < 0.5 are shown in gray, with symbol shape and size depending on the gas physical parameters, as listed in the upper left corner of the figure. A least-squares fit to the model predictions yields log I([CII]HII) = 1.068 × log I([NII]) + 0.645, shown by the gray straight line. For comparison, the expected [CII] emission contribution from H II regions in the Abel (2006) models for 1.0 Z and 0.1 Z are shown as dashed and dotted black straight lines, respectively.

In the text
thumbnail Fig. A.1

Graphical description of the approach adopted to derived Ig, Is, and Ii integrated intensities for the [CII] and [NII] line observations, as described in Sect. 3. In the example, the spectra in the upper, middle, and lower panels correspond to the [CII] line at the W2 position. The spectra were normalized to the peak antenna temperature for display purposes.

In the text
thumbnail Fig. A.2

Comparison of [CII] (open blue), [NII] (filled black), [CI](2–1) (filled dark gray), [CI](1–0) (filled gray), and C18 O(1–0) (open red) spectra for positions P1 and P2 in Table 1. Antenna temperatures have been scaled by the factor shown next to the line’s name for display purposes. Gaussian component fits to the [CII] and C18 O(1–0) spectra are shown in blue and red, respectively. Only dashed-line Gaussian components in the [CII] spectra are used for the analysis in this work as only for them can a counterpart in the C18 O(1–0) optically thin emission be identified. Vertical dotted black lines show the radial velocity limits for which integrated intensities were calculated for the [NII] line by directly adding the observed emission, multiplied by the spectral channel width, to that LSR velocity range.

In the text
thumbnail Fig. A.3

Comparison of [CII] (open blue), [NII] (filled black), [CI](2–1) (filled dark gray), [CI](1–0) (filled gray), and C18 O(1–0) (open red) spectra for positions W1 and W2 in Table 1. Antenna temperatures have been scaled by the factor shown next to the line’s name for display purposes. Gaussian component fits to the [CII] and C18 O(1–0) spectra are shown in blue and red, respectively. Only dashed-line Gaussian components in the [CII] spectra are used for the analysis in this work as only for them can a counterpart in the C18 O(1–0) optically thin emission be identified. Vertical dotted black lines show the radial velocity limits for which integrated intensities were calculated for the [NII] line by directly adding the observed emission, multiplied by the spectral channel width, to that LSR velocity range.

In the text
thumbnail Fig. A.4

Comparison of [CII] (open blue), [NII] (filled black), [CI](2–1) (filled dark gray), [CI](1–0) (filled gray), and C18 O(1–0) (open red) spectra for positions G0.10+0.02 and G0.07+0.04 in Table 1. Antenna temperatures have been scaled by the factor shown next to the line’s name for display purposes. Gaussian component fits to the [CII] and C18 O(1–0) spectra are shown in blue and red, respectively. Only dashed-line Gaussian components in the [CII] spectra are used for the analysis in this work as only for them can a counterpart in the C18 O(1–0) optically thin emission be identified. Vertical dotted black lines show the radial velocity limits for which integrated intensities were calculated for the [NII] line by directly adding the observed emission, multiplied by the spectral channel width, to that LSR velocity range.

In the text
thumbnail Fig. A.5

Comparison of [CII] (open blue), [NII] (filled black), [CI](2–1) (filled dark gray), [CI](1–0) (filled gray), and C18 O(1–0) (open red) spectra for positions E2N and E2S in Table 1. Antenna temperatures have been scaled by the factor shown next to the line’s name for display purposes. Gaussian component fits to the [CII] and C18 O(1–0) spectra are shown in blue and red, respectively. Only dashed-line Gaussian components in the [CII] spectra are used for the analysis in this work as only for them can a counterpart in the C18 O(1–0) optically thin emission be identified. Vertical dotted black lines show the radial velocity limits for which integrated intensities were calculated for the [NII] line by directly adding the observed emission, multiplied by the spectral channel width, to that LSR velocity range.

In the text
thumbnail Fig. A.6

Comparison of [CII] (open blue), [NII] (filled black), [CI](2–1) (filled dark gray), [CI](1–0) (filled gray), and C18 O(1–0) (open red)spectra for position E1 in Table 1. Antenna temperatures have been scaled by the factor shown next to the line’s name for display purposes. Gaussian component fits to the [CII] and C18 O(1–0) spectra are shown in blue and red, respectively. Only dashed-line Gaussian components in the [CII] spectra are used for the analysis in this work as only for them can a counterpart in the C18 O(1–0) optically thin emission be identified. Vertical dotted black lines show the radial velocity limits for which integrated intensities were calculated for the [NII] line by directly adding the observed emission, multiplied by the spectral channel width, to that LSR velocity range.

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.