Free Access
Issue
A&A
Volume 564, April 2014
Article Number A101
Number of page(s) 6
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201323281
Published online 11 April 2014

© ESO, 2014

1. Introduction

The star formation rate in the Galaxy may be related to the pressure of the interstellar medium (ISM), which itself is a function of the interplay of dynamical processes and energy sources on the interstellar gas. It has also been suggested that ISM pressure plays a role in the formation of giant molecular clouds (Blitz & Rosolowsky 2004, 2006). Thus, the vertical (z) distribution of the various ISM gas components is an important parameter for understanding these dynamical processes throughout the Milky Way. In the Galaxy, interstellar clouds are distributed in a thin disk about the mid-plane at b = 0° and their scale height depends, in hydrostatic equilibrium, on a number of factors, including thermal pressure, the random motion of the clouds, magnetic pressure, ionization pressure, and the gravitational force of the stars and gas in the disk. Thus, determining the z-distribution, which may be different for various ISM components, provides information about these parameters. The scale height of the diffuse atomic hydrogen clouds is known from extensive maps of the H i 21-cm line (cf. Boulares & Cox 1990; Dickey & Lockman 1990) and that for the dense molecular hydrogen clouds from large scale maps of the 12CO J = 1 → 0 rotational line (e.g. Sanders et al. 1984; Dame et al. 1987; Bronfman et al. 1988; Clemens et al. 1988; Malhotra 1994; Dame et al. 2001; Jackson et al. 2006). The scale height for H i and H2 (as traced by CO) clouds are different and each varies by a factor of ~3 across the inner Galaxy, and increases in the outer Galaxy (Narayan & Jog 2002). In addition to tracing gas that can be observed in H i and CO, the 1.9 THz emission from ionized carbon, [C ii], traces the H2 gas where carbon is ionized but little, or no, CO or neutral carbon is found (the CO-dark H2 gas), and also traces the warm ionized medium (WIM).

The scale height for clouds traced by [C ii] is not well established because the necessary spectral line surveys of its 158-μm line have not, until recently been available. The COBE FIRAS instrument made the only large-scale survey of spectrally unresolved [C ii] (Wright et al. 1991; Bennett et al. 1994), however, COBE with its 7° beam and ~1000 km s-1 velocity resolution, is unable to resolve the latitudinal distribution. There are two moderate-scale Galactic surveys of spectrally unresolved [C ii], the Far-Infrared Line Mapper (FILM) onboard the Infrared Telescope in Space (IRTS; Shibai et al. 1994; Makiuti et al. 2002) and the Balloon-borne Infrared Carbon Explorer (BICE; Nakagawa et al. 1998), and an earlier small-scale survey with the Balloon-borne Infrared Telescope (BIRT; Shibai et al. 1991). However, there is only one spectrally resolved survey, the Herschel open time key program, Galactic Observations of Terahertz C+, hereafter GOT C+ (see Langer et al. 2010; Pineda et al. 2013; Langer et al. 2014). FILM and BICE had an angular resolution of order ten to fifteen arcminutes, sufficient to determine the latitudinal b distribution of [C ii], but the velocity resolution was, at best, ~175 km s-1 for BICE (for FILM it was ~750 km s-1 and for BIRT 143 km s-1). Only the GOT C+ survey had the spectral resolution (<1 km s-1) sufficient to resolve the velocity structure of individual gas clouds and thus locate their Galactic radius using a position–velocity rotation curve. However, GOT C+ surveyed [C ii] sparsely in longitude, l, and latitude, b.

Here we determine the average scale height distribution of [C ii] in z by combining the high spectral resolution [C ii] radial distribution from GOT C+ at b = 0° with the BICE latitudinal angular distribution in b by adopting a hydrostatic model result for the distribution in z. We use the BICE survey to determine the scale height of [C ii] emission because it had better coverage in both longitude and latitude in the Galactic disk than the BIRT and FILM surveys. BICE had an angular resolution of 15, and spectral resolution of ~175 km s-1.

The GOT C+ survey contains several hundred lines-of-sight of spectrally resolved [C ii] emission throughout the Galactic disk from l = 0° to 360° and b = 0°, ±0.5°, and ±1.0°. However, because GOT C+ is a sparse survey it does not have sufficient coverage in latitude b to derive a smooth continuous distribution in the vertical distribution z. In contrast, the BICE survey had insufficient spectral resolution, but had much better coverage in b, but only observed up to latitudes b = ± 3° and only within longitude 350° ≤ l ≤ 25°. GOT C+ has a 3σ sensitivity ~0.24 K km s-1 (Langer et al. 2014) which, over the bandwidth of the velocity resolution of BICE, corresponds to 7.4×10-6 erg s-1 cm-2 sr-1. BICE has a 3σ detection limit ~2×10-5 erg s-1 cm-2 sr-1 (Nakagawa et al. 1998). Thus GOT C+ is almost three times more sensitive than BICE and capable of detecting the [C ii] emission seen by BICE.

We also use results from FILM (Shibai et al. 1996; Makiuti et al. 2002), which observed [C ii] at higher latitudes than BICE, to recalibrate the results of Nakagawa et al. (1998) for | b | = 3° to . FILM observed [C  ii] along a great circle crossing the plane at l = 50° (inner Galaxy) and 230° (outer Galaxy), but measured the [C ii] intensity for larger latitudes than BICE. However, Shibai et al. (1996) and Makiuti et al. (2002) smoothed their data to 1° to improve their sensitivity, insufficient to use at low latitudes to determine the scale height in the disk.

Nakagawa et al. (1998) compared the latitudinal distribution of [C ii] from BICE with those of H i, 12CO, and far-infrared dust emission, finding that H i had the largest distribution in b followed by [C ii] and then far-infrared dust emission, and that CO had the smallest distribution in b. However, without knowing where the [C ii] emission came from they could not assign a spatial scale height to the gas traced by the [C ii] 158-μm line. We begin with a summary of the BICE and GOT C+ [C ii] distributions and then derive an approximate relationship between the radial and b distributions that allow us to deconvolve the [C ii] distribution in z in the disk. Finally, we compare the [C ii] scale height with those of CO and H i and discuss its implications for understanding the sources of [C ii] emission.

2. The z-distribution of [C ii]

To determine the latitudinal distribution of [C ii] Nakagawa et al. (1998) averaged the [C ii] line intensity (measured in erg s-1 cm-2 sr-1) only over longitudes 5° < l < 25° in order to avoid the Galactic center, and calculated a relative intensity as a function of b such that the peak of this distribution is unity. Here we define the relative intensity as I(b) /I(bc), where I(b) is the intensity along b, and I(bc) is the intensity at the peak of the distribution located at bc (i.e. this term accounts for any offset in the peak from b = 0°). However, Nakagawa et al. (1998) set the normalized intensity to zero at b = ±4°, thus suppressing potential contributions from | b | = 3° to . The relative intensity distribution from BICE is plotted in their Fig. 9, and recreated here in Fig. 1. Nakagawa et al. (1998) also plotted the corresponding relative intensity of other ISM gas and dust tracers: H i, 12CO(1 → 0), and far-infrared continuum from IRAS. They derived the full width at half maximum (FWHM) in b, for these four tracers and their results are summarized here in Table 1. In Fig. 1 it can be seen that the peak in the [C ii] relative intensity is shifted slightly below the plane to bc ~ −0.2°; the far-IR emission also peaks there (see their Fig. 9), but their plot of the CO and H i peak at b = 0°.

thumbnail Fig. 1

Distribution of the normalized intensity of [C ii], I(b) /I(bc), as a function of latitude as derived by Nakagawa et al. (1998) from the BICE survey (solid line). The modified distribution using the FILM latitudinal data to calibrate the normalized intensity at b ± 4° (dashed line).

Table 1

b(FWHM) of gas and dust tracers: Nakagawa et al. (1998).

The assumption that I(b = ±4°) /I(bc) = 0 is in conflict with the FILM observations of [C ii] and could introduce a slight error in determining the FWHM. Shibai et al. (1996) observed that [C ii] emission in the FILM latitudinal scans extended out to perhaps ±50° (see their Fig. 1, which plots the [C ii] intensity distribution in latitude smoothed over a 1° beam). They interpreted the latitudinal [C ii] distribution as having two components: (1) concentrated emission in the disk; and, (2) a weaker component, decreasing slowly with b beyond ~10°. Makiuti et al. (2002) analyzed the distribution at large b and conclude that it mainly traces the WIM at high latitudes. We have fit the wings of the FILM inner Galaxy scan from to 10° in order to recalibrate the BICE distribution such that the BICE relative intensity corresponds to the relative intensity ratio observed by FILM at ±4°. We do not use the FILM fit beyond ± 4° because the FILM observing path results in b being a strong function of longitude l and it does not cover the regions needed to compare GOT C+ with the BICE survey, and because at large b FILM is surveying the Galaxy near the solar radius. The revised BICE distribution is shown in Fig. 1 and labeled BICE + FILM. This recalibration makes very little difference to the bulk of the [C ii] distribution, but, as will be seen later, provides a much better fit to the wings at | b | ≥ 2°.

The GOT C+ survey observed spectrally resolved [C ii] along 151 lines-of-sight covering 0° < l < 360° at b = 0°. The details of the observing mode and data reduction are discussed in Pineda et al. (2013), and representative spectra can be viewed there and in Langer et al. (2014). In Pineda et al. (2013) the [C ii] spatial-velocity maps were used along with a Galactic rotation curve to assign the intensity as a function of Galactic radius. These intensities were then summed in rings about the Galactic center and then used to calculate the azimuthally averaged emissivity for [C ii], ϵ[ CII ](Rgal) as a function of Galactic radius, Rgal, in units of K km s-1 kpc-1, except for the innermost 0.5 kpc, because this region is undersampled. In Fig. 2 we reproduce the radial distribution for ϵ[ CII ](Rgal) from Fig. 7 in Pineda et al. (2013). In Fig. 2 it can be seen that the [C ii] emissivity peaks in the molecular ring at about 5 kpc and most of the emission comes from Rgal = 3 to 7 kpc.

Our approach assumes that the GOT C+ radial distribution sees all the [C ii] components observed by BICE including the WIM, which is believed to more prominent at higher latitudes than in the plane. In the Galactic plane for b = 0° the contribution of the WIM component in the radial profile was estimated by Pineda et al. (2013) to be only ~4% using the electron abundance of low density ionized gas in the plane as given by the NE2001 code (Cordes & Lazio 2002). Note that this value is only an estimate and is not directly observed in the GOT C+ data. However, Velusamy et al. (2012) used the GOT C+ data to detect emission from the compressed WIM along selected lines of sight corresponding to the spiral arm tangencies. Combined with its higher sensitivity and adequate sampling in the plane, the distribution of [C ii] emission in the GOT C+ radial profile contains all [C ii] features as seen by BICE, including the emission from the WIM. Therefore we believe that the combined effects of all [C ii] components will be fully represented in our solution for the z-scale derived below.

thumbnail Fig. 2

Galactic radial emissivity of [C ii], ϵ[ CII ](Rgal), derived by Pineda et al. (2013) using the GOT C+ [C ii] survey at b = 0°.

We can determine the average z distribution in the inner Galactic disk by combining the GOT C+ radial distribution of ϵ[ CII ](Rgal) for b = 0° with the modified BICE + FILM latitudinal distribution if we know the functional form of the emissivity in z, and if we assume that the distribution in z is independent of location in the inner Galactic disk. For the first assumption, models of hydrostatic equilibrium of gas in the ISM suggest a Gaussian distribution for many ISM components, however, the second assumption can only be correct on average, because the pressure in the ISM varies across the Galaxy, and the observed scale heights for CO and H i vary with Galactic radius.

Adopting these two assumptions the averaged [C ii] intensity observed by BICE + FILM along a line-of-sight is just the integral of the emissivity along a path length s from the solar system, as illustrated in Fig. 3 for l = 0°. If we know the value of the emissivity at z = 0 we can determine its value at height z for a given scale height of the Gaussian distribution. Thus the intensity measured by BICE + FILM as a function of b is related to the [C ii] emissivity, ϵ[ CII ](Rgal,l,x,z), as a function of the distance to the source in the plane, x, the height above the plane, z, along a line-of-sight longitude l, and latitude, b. Most of the observed [C ii] intensity measured by BICE comes from the near side of the inner Galaxy for | b | ≥ 1° because at higher latitudes the line of sight passes far above the plane on the far side of the Galaxy. For example, at b = 1° the line of sight, s, passes ~150 pc above the Galactic center, while for b = 3°, it is ~450 pc above the plane, and on the far side of the molecular ring it is >600 pc above b = 0°.

thumbnail Fig. 3

Schematic showing two representative lines of sight (dashed lines) from the solar system through the Galactic disk (GC labels the Galactic center and b labels the latitude). The intensity of [C ii] is the integral of the emissivity, ϵ, along the path length s (see text). The actual lines of sight used to calculate the intensity cover l = 5° to 25°.

We can write the intensity as a function of latitude, I(b), along a given line of sight in terms of the integration of the emissivity as follows, I[CII](b)=0smaxϵ[CII](s)ds\begin{equation} I_{\rm [CII]}(b)=\int_0^{s_{\rm max}} \epsilon_{\rm [CII]}(s){\rm d}s \end{equation}(1)where the integral of the emissivity ϵ[ CII ](s) is along a line s from the solar system, smin = 0 to a maximum distance smax, where the emissivity is small. As illustrated in Fig. 3 we can relate the intensity of emission at each location along s to the emission along the plane at b = 0°, by changing the integration along s to x in the plane. Substituting x = scos(b) in Eq. (1) we get, I[CII](l,b)=0xmaxϵ[CII](x)f(z)(cos(b))-1dx\begin{equation} I_{\rm [CII]}(l,b)=\int_0^{x_{\rm max}} \epsilon_{\rm [CII]}(x )f(z)(\cos\,(b))^{-1} {\rm d}x \end{equation}(2)where f(z) is the emissivity distribution in z.

Table 2

Tracer fitting parameters.

We have chosen to integrate up to xmax = 28 kpc, which is on the far side of the Galaxy and where the GOT C+ [C ii] emissivity is very small. However, for all practical purposes, as discussed above, most of the contribution to the intensity comes from the near side of the Galaxy, except for the very lowest values of b. The usual form for f(z) derived from equations of hydrostatic equilibrium balancing the ISM pressure with the gravitational force of the stars and gas in the Galaxy is a Gaussian (cf. Spitzer 1978), although other dynamical processes (e.g. supernova, outflows) can lead to non-Gaussian terms. Here we assume a Gaussian distribution for the [C ii] emissivity of the form, f(z)=f(zc)e0.5((zzc)/z0)2\begin{equation} f(z) = f(z_{\rm c}){\rm e}^{-0.5((z-z_{\rm c})/z_0)^2} \end{equation}(3)where z0 is the scale height, zc accounts for an offset in the peak of the distribution, and f(zc) normalizes the distribution to unity. The FWHM([C ii]) is equal to 2(2ln2)z00.5\hbox{$^{0.5}z_0$}. We substitute this form in Eq. (2) and, as can be seen in Fig. 3, rewrite the integral using z = xcos(b). The relative intensity, F(b) = I(b) /I(bc) is given by, F(b)=0xmaxϵ[CII](x)e0.5((xsinbzc)/z0)2cos(b)-1dx0xmaxϵ[CII](x)e0.5(zc/z0)2dx·\begin{equation} \label{eq:Equation_F(b)} F(b)=\frac{\int_0^{x_{\rm max}} \epsilon_{\rm [CII]}(x){\rm e}^{-0.5((xsinb - z_{\rm c})/z_0)^2}\cos(b)^{-1}{\rm d}x}{\int_0^{x_{\rm max}} \epsilon_{\rm [CII]}(x){\rm e}^{-0.5(z_{\rm c}/z_0)^2}{\rm d}x}\cdot \end{equation}(4)where cos(b = 0°) is unity in the denominator.

As mentioned above, Nakagawa et al. (1998) calculated the relative intensity as a function of b by averaging the BICE observations over longitudes from 5° to 25°, to avoid the Galactic center, thus omitting [C ii] over a region with Galactic radius ~0.75 kpc. We solve for the average scale height covering longitudes 5° to 25°, where we use the GOT C+ radial profiles to calculate the intensity along l for b = 0°. We also assume that [C  ii] emission is optically thin because most of the GOT C+ [C ii] spectra have a main beam temperature much less than the kinetic temperature (c.f. discussion of C+ excitation and [C ii] radiative transfer in Goldsmith et al. 2012).

We iterated on the two parameters, z0 and zc in Eq. (4) to minimize the rms deviation of the model averaged over five longitudes compared to the modified BICE + FILM [C ii] distribution. The best fit is given by z0 = 73 pc and zc = −28 pc, and is listed in Table 2, along with the corresponding FWHM([C ii]) = 172 pc. We also fit the original BICE distribution without the FILM correction and find the same offset and a slightly different FWHM = 170 pc. The fitting parameters to the BICE and the BICE + FILM distributions do not differ significantly as the central ± 2° dominates the quality of the fit. However, the fit is much better for the modified BICE + FILM distribution for | b | > 2°. In Fig. 4 we show a plot of the relative intensity for the BICE + FILM combination and for just the original BICE distribution, compared to the fitted distribution as a function of b. It can be seen that our fit is very good over the complete range b = −4° to +, differing by at most ~15% from the modified BICE + FILM distribution, except at the bumps in the shoulders (b ~ ±2°) where the difference is as much as 30%.

thumbnail Fig. 4

Model fit (dotted line) of the normalized intensity using the GOT C+ observations at b = 0° compared to the normalized intensity from the modified BICE + FILM profile (solid line) versus latitude b. The unmodified BICE profiles are shown for comparison (dashed line).

Nakagawa et al. (1998) calculated the total flux measured by BICE over the entire mapped angular region | b | ≤ 3° and 350°l ≤ 25° to be 6 × 10-6 erg s-1 cm-2. The BICE calibration errors are ±35% and include an estimated uniform offset of 2×10-5 erg s-1 cm-2 sr-1, comparable to the calibration uncertainty (see their Sect. 3.1). FILM observations are better calibrated than those of BICE and Makiuti et al. (2002) compared FILM with BICE and found that the [C ii] line intensity for FILM is ~65% of BICE for I([C ii]) >10-4 erg s-1 cm-2 sr-1 and ~85% for I([C ii]) < 10-4 erg s-1 cm-2 sr-1. To confirm that GOT C+ and BICE are sensitive to the same ISM gas components we calculated the total flux that would be detected by GOT C+ over this same angular area, | b | ≤ 3° and 350°l25°, using the emissivity in the plane (b = 0°) in Pineda et al. (2013) and the Gaussian distribution in z with a FWHM = 172 pc. We calculate from GOT C+ a total flux ~4.5 × 10-6 erg s-1 cm-2. This value is about 25% less than calculated in the BICE survey, however, it is in good agreement if we recalibrate BICE using the calibration from FILM. The agreement between GOT C+ and BICE is well within the errors of the BICE result, and supports our assumption that GOT C+ and BICE trace the same ISM gas components.

thumbnail Fig. 5

Gaussian fits for the normalized intensities of [C ii], CO, and H i as a function of z. The fits include the measured offsets zoffset for CO and [C ii], and an inferred offset for H i equal to that for CO. Towards the inner Galaxy CO has, on average, the narrowest distribution in the disk, followed by [C ii], and H i has the thickest distribution.

The FWHM for CO has been estimated by several authors using different Galactic surveys, including Sanders et al. (1984), Dame et al. (1987, 2001), Bronfman et al. (1988), Clemens et al. (1988), Malhotra (1994), and Jackson et al. (2006). They all find about the same result that CO emission peaks between 4 and 7 kpc, and in the Galactic center, and that the scale height increases with increasing Galactic radius beyond 2 kpc and is offset below b = 0°. Sanders et al. (1984) found that the FWHM (or 2z1 /2 in their notation) ranged from 60 to 140 pc for Rgal increasing from 3 to 8 kpc. To compare the CO distribution with the average derived for [C ii] we calculated the average FWHM(CO) using their results from 3 to 8 kpc where the [C ii] is largest. We find FWHM ⟩ = ~ 110 pc and zc ~ −25 pc, and these values are listed in Table 2. Their CO offset below the plane is very close to what we derived for [C ii], however, the CO FWHM is noticeably smaller than that for [C ii].

The H i distribution in the Galaxy has been discussed by many authors (c.f. Boulares & Cox 1990; Dickey & Lockman 1990; Narayan & Jog 2002; Kalberla & Dedes 2008; Kalberla & Kerp 2009). Here we will use the parameters from Dickey & Lockman (1990) for fitting the distribution of H i relative intensity with z. H i is roughly constant from Rgal = 4 to 8 kpc, but the distribution in z is complicated by having several components, some of which extend into the halo. Dickey & Lockman (1990) find that the best estimate of the z mean density distribution, n(z), is given by two Gaussians with peak mean densities of n1(0) = 0.395 and n2(0) = 0.107 cm-3 and FWHM1 = of 212 and FWHM2 = 530 pc, plus an exponential with a peak mean density, n3(0) = 0.064 cm-3 and scale height 403 pc. In Table 2 we list their two Gaussian parameters and one exponential parameter, but to calculate a relative H i intensity we normalize their contributions to the total intensity using the corresponding central mean densities to weight the contributions, where the total F(z) = F1(z) + F2(z) + F3(z), with Fi(bc) = ni(bc) /ntot(bc), and ntot(bc) = n1(bc) + n2(bc) + n3(bc). We have included an offset zc =−25 pc for H i in Table 2, to facilitate comparison of the relative distributions of the three gas tracers. The spatial resolution of H i Galactic surveys is not high enough to resolve clearly an offset zc(HI) of order 25 pc, but there is some indication that the solar system is offset from the warped Galactic plane as seen in the Leiden-Argentine-Bonn 21-cm survey (see Fig. 3 in Kalberla & Kerp 2009).

In Fig. 5 we plot the Gaussian distributions of the relative intensities as a function of z for [C ii], CO, and H i using the average FWHM for each of the components listed in Table 2. As can be seen in Fig. 5 the average distribution of 12CO as a function of z in the Galactic disk is narrower than [C ii], and both are narrower than that for H i. The distribution for [C ii] consists mainly of a disk component confined to ±200 pc.

3. Discussion

There is a very simple explanation why the scale height for [C ii] is more than that of CO but less than that of H i. While H i traces mainly the atomic hydrogen clouds in the Galaxy and CO traces the dense molecular clouds, [C ii] traces both of these regions, as well as the diffuse molecular clouds that have [C ii] but no CO emission (CO-dark H2 clouds), and the WIM. The solution for the density distribution of clouds in hydrostatic equilibrium is a Gaussian function, exp (− 0.5(z/z0)2), with the scale factor, z0 proportional to the velocity dispersion, vi20.5\hbox{$\langle v_i^2\rangle^{0.5}$}, where i labels the ISM cloud component. The diffuse atomic clouds with their lower densities have higher velocity dispersions, vH20.5\hbox{$\langle v_{\rm H}^2\rangle^{0.5}$}~10 km s-1, while those for the denser CO clouds have vCO20.5\hbox{$\langle v_{\rm CO}^2\rangle^{0.5}$}~5 km s-1 (see discussion and models in Narayan & Jog 2002). Thus in hydrostatic equilibrium, for an equal ISM pressure, the H i clouds have a larger scale height than those of the denser CO clouds.

Pineda et al. (2013) found that the PDRs of dense molecular clouds emit ~43% of the total [C ii] throughout the plane at b = 0°, diffuse atomic hydrogen clouds ~23%, diffuse molecular clouds (CO-dark H2 clouds) ~30%, and for the WIM estimated ~4%. In the inner Galaxy, Rgal< 9 kpc, where we evaluate the [C ii] distribution, these percentages are only slightly different. Therefore, it is no surprise that the distribution of [C ii], which arises from all ISM components, would have a distribution intermediate between that of the dense CO clouds and less dense atomic H i clouds and the WIM. Thus [C ii] traces a mixture of clouds of different mass, density, and velocity dispersion, some of which are also traced by H i and CO, and the WIM. The different scale heights then depend on the different physical properties and energetics of the clouds that enter into the hydrostatic mechanisms responsible for the distribution of gas in the plane.

As seen in Fig. 5, the derived z distribution in [C ii] does not follow that of H i for z greater than about ~200 pc. At this height there are few dense molecular clouds as traced by CO and likely very few diffuse molecular clouds (CO-dark H2 clouds), so any [C ii] would have to come from the H i clouds and/or the WIM. Makiuti et al. (2002) compared the distribution of [C ii] and H i at high latitudes and conclude that the [C ii] emission comes primarily from the WIM (see their Fig. 3). However, the FILM results at high latitude are limited to a region near the solar radius and cannot be extrapolated across the Galaxy. The combined GOT C+ and BICE data suggest that this conclusion also holds for the inner Galaxy as well, because H i clouds high above the plane have low densities and are not likely to emit [C ii] efficiently. This low emissivity is due to the difference in the excitation conditions for H i and [C ii]. The intensity of H i is proportional to its column density in the optically thin regime, I([HI])=5.5×10-19N(HI),\begin{equation} I([{\rm H\,I}]) = 5.5\times 10^{-19}N({\rm H\,I}), \end{equation}(5)in units of (K km s-1), and is relatively insensitive to density and kinetic temperature. In contrast the [C ii] emission is very sensitive to kinetic temperature, Tkin because the energy of the upper level P3/22\hbox{$^{\rm 2}P_{\rm 3/2}$}, Eu/k = 91.25 K is typically higher than the gas temperature in neutral clouds, and density where the atomic, n(H), and molecular, n(H2), hydrogen densities are much lower than the critical densities for thermalizing the C+, ncr(H) ~ 3000 cm-3 and the recently revised value ncr(H2) ~ 4500 cm-3 (Wiesenfeld & Goldsmith 2014). The radiative transfer equation for the [C ii] intensity for optically thin emission is given in Goldsmith et al. (2012) and Langer et al. (2014), and, in the diffuse clouds, such that the intensity can be written as, Ij([CII])=1.73×10-16eΔE/kTn(j)ncr(j)Nj(C+),\begin{equation} I_j([{\rm C\,II}]) = 1.73\times10^{-16} {\rm e}^{{-\Delta}E/kT} \frac{ n(j)}{n_{\rm cr}(j)} N_j({\rm C^+}) , \end{equation}(6)where I is in units of K km s-1 and the index j labels H i or H2. Therefore, while the H i intensity depends only on the column density of atomic hydrogen, the [C ii] intensity also depends very sensitively on the density and will be much smaller in low density atomic hydrogen clouds above the plane.

The scale height for [C ii] derived here depends on the radial distribution derived from the GOT C+ sampling at b = 0°, which as noted above is a sparse sample. The premise of the GOT C+ survey was that a well designed unbiased sampling in longitude would represent statistically the distribution of [C ii] in the Galactic plane. Therefore, the fact that GOT C+, along with our model of the z distribution, reproduces the total flux observed by BICE (rescaled to the FILM calibration) supports this approach.

Another potential uncertainty is the adoption of a Galactic rotation curve in Pineda et al. (2013) to locate the source of the [C ii] emission. In Langer et al. (2014) we adopted a rotation curve based on gas-flow hydrodynamical models to assign a distance based on velocity. We found that it made a difference mainly in the inner Galaxy, | l | ≤ 6°, but this region is mostly excluded from the BICE analysis (see above). There are also regions where clouds have peculiar velocities due to Galactic dynamics, where the rotation curve may assign the wrong distance. For example, Zhang et al. (2014) find non-rotational cloud motions at the end of the Galactic bar from parallax observations of masers at l ~ 30°. While 30° is outside of the longitudinal range observed by BICE, this region contributes to our GOT C+ data set. We cannot quantitatively assess the error introduced into our radial distribution but note that, because we average lines of sight from all across the Galaxy, the edges of the bar contribute a small fraction of the emission in any given ring.

We cannot calculate the radial dependence for the [C ii] FWHM from the GOT C+ survey without a better sampling in b. However, to gain some insight on the effect of a variable FWHM(Rgal), we assume that it varies similar to that for CO. Sanders et al. (1984) and Clemens et al. (1988) found that the CO scale height varied roughly as Rgal0.5\hbox{$R_{\rm gal}^{0.5}$} between 3 and 9 kpc. We replaced zo in Eq. (3) with one that varied Rgal0.5\hbox{$ R_{\rm gal}^{0.5}$} for Rgal> 3 kpc at the value for 3 kpc. We solved for the scale factor that best fit the BICE distribution in b, similar to what was done to determine an average scale factor. We find that the best fit is given by, FWHM(Rgal) = 172(Rgal/ 4.7)0.5 pc. Thus the average value for FWHM for [C ii] of 172 pc corresponds to the radial solution at ~4 to 5 kpc, essentially in the molecular ring. For the assumed radial dependence, the FWHM ranges from ~140 pc to ~230 pc over Galactic radii 3 kpc to 8 kpc.

4. Summary

We have combined the GOT C+ spectrally resolved [C ii] survey in the Galactic plane at b = 0° with the latitudinal distribution derived from the BICE survey of spectrally unresolved [C ii] to derive, for the first time, the average scale height of [C ii] over the inner Galactic plane. GOT C+ is slightly more sensitive than BICE and the total flux measured by GOT C+ is close to that of BICE given the uncertainties of the BICE calibration. Therefore

these two surveys are likely tracing the same ISM gas components. The average distribution in the inner Galactic disk is well fit by a single Gaussian with FWHM([C ii]) = 172 pc and an offset − 28 pc below the plane (b = 0°).

In this paper we find that the [C ii] distribution is larger in z than that of CO, but smaller than H i. The origin of the [C ii] emission has been attributed to different sources by various authors based on the spectrally unresolved surveys. However, the result here suggests to us a more complicated picture with [C ii] tracing a mix of ISM cloud categories. The GOT C+ data for b ≠ 0° may be able to give some insights on the distribution of the different ISM components, but to determine completely the distribution of [C ii] for the separate ISM components as a function of Galactic radius and z, we will need more detailed spectrally resolved latitudinal maps across the Galaxy and with finer steps in b. We also need to extend the spectral line observations to higher values of b than observed in the GOT C+ survey to understand the contributions of the warm ionized medium and low density high latitude H i clouds to the [C ii] emissivity above the disk.

Acknowledgments

We thank the referee for a careful reading of the manuscript and several suggestions that improved the discussion. We also thank P. F. Goldsmith for constructive comments and edits. This work was performed at the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.

References

  1. Beichman, C. A., Neugebauer, G., Habing, H. J., Clegg, P. E., & Chester, T. J. 1988, Infrared astronomical satellite (IRAS) catalogs and atlases, Explanatory supplement, 1 [Google Scholar]
  2. Bennett, C. L., Fixsen, D. J., Hinshaw, G., et al. 1994, ApJ, 434, 587 [CrossRef] [Google Scholar]
  3. Blitz, L., & Rosolowsky, E. 2004, ApJ, 612, L29 [NASA ADS] [CrossRef] [Google Scholar]
  4. Blitz, L., & Rosolowsky, E. 2006, ApJ, 650, 933 [NASA ADS] [CrossRef] [Google Scholar]
  5. Boulares, A., & Cox, D. P. 1990, ApJ, 365, 544 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bronfman, L., Cohen, R. S., Alvarez, H., May, J., & Thaddeus, P. 1988, ApJ, 324, 248 [NASA ADS] [CrossRef] [Google Scholar]
  7. Clemens, D. P., Sanders, D. B., & Scoville, N. Z. 1988, ApJ, 327, 139 [NASA ADS] [CrossRef] [Google Scholar]
  8. Cordes, J. M., & Lazio, T. J. W. 2002, unpublished [arXiv:astro-ph/0207156] [Google Scholar]
  9. Dame, T. M., Ungerechts, H., Cohen, R. S., et al. 1987, ApJ, 322, 706 [NASA ADS] [CrossRef] [Google Scholar]
  10. Dame, T. M., Hartmann, D., & Thaddeus, P. 2001, ApJ, 547, 792 [NASA ADS] [CrossRef] [Google Scholar]
  11. Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  12. Goldsmith, P. F., Langer, W. D., Pineda, J. L., & Velusamy, T. 2012, ApJS, 203, 13 [Google Scholar]
  13. Hartmann, D., & Burton, W. B. 1997, Atlas of Galactic Neutral Hydrogen (Cambridge, UK: Cambridge University Press) [Google Scholar]
  14. Jackson, J. M., Rathborne, J. M., Shah, R. Y., et al. 2006, ApJS, 163, 145 [NASA ADS] [CrossRef] [Google Scholar]
  15. Kalberla, P. M. W., & Dedes, L. 2008, A&A, 487, 951 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Kalberla, P. M. W., & Kerp, J. 2009, ARA&A, 47, 27 [NASA ADS] [CrossRef] [Google Scholar]
  17. Langer, W. D., Velusamy, T., Pineda, J. L., et al. 2010, A&A, 521, L17 [Google Scholar]
  18. Langer, W. D., Velusamy, T., Pineda, J. L., Willacy, K., & Goldsmith, P. F. 2014, A&A, 561, A122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Makiuti, S., Shibai, H., Nakagawa, T., et al. 2002, A&A, 382, 600 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Malhotra, S. 1994, ApJ, 433, 687 [NASA ADS] [CrossRef] [Google Scholar]
  21. Nakagawa, T., Yui, Y. Y., Doi, Y., et al. 1998, ApJS, 115, 259 [NASA ADS] [CrossRef] [Google Scholar]
  22. Narayan, C. A., & Jog, C. J. 2002, A&A, 390, L35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Pineda, J. L., Langer, W. D., Velusamy, T., & Goldsmith, P. F. 2013, A&A, 554, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Sanders, D. B., Solomon, P. M., & Scoville, N. Z. 1984, ApJ, 276, 182 [NASA ADS] [CrossRef] [Google Scholar]
  25. Shibai, H., Okuda, H., Nakagawa, T., et al. 1991, ApJ, 374, 522 [NASA ADS] [CrossRef] [Google Scholar]
  26. Shibai, H., Yui, M., Matsuhara, H., et al. 1994, ApJ, 428, 377 [NASA ADS] [CrossRef] [Google Scholar]
  27. Shibai, H., Okuda, H., Nakagawa, T., et al. 1996, PASJ, 48, L127 [NASA ADS] [Google Scholar]
  28. Spitzer, L. 1978, Physical processes in the interstellar medium (New York: Wiley-Interscience Publication) [Google Scholar]
  29. Velusamy, T., Langer, W. D., Pineda, J. L., & Goldsmith, P. F. 2012, A&A, 541, L10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Wiesenfeld, L., & Goldsmith, P. F. 2014, ApJ, 780, 183 [NASA ADS] [CrossRef] [Google Scholar]
  31. Wright, E. L., Mather, J. C., Bennett, C. L., et al. 1991, ApJ, 381, 200 [NASA ADS] [CrossRef] [Google Scholar]
  32. Zhang, B., Moscadelli, L., Sato, M., et al. 2014, ApJ, 781, 89 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

b(FWHM) of gas and dust tracers: Nakagawa et al. (1998).

Table 2

Tracer fitting parameters.

All Figures

thumbnail Fig. 1

Distribution of the normalized intensity of [C ii], I(b) /I(bc), as a function of latitude as derived by Nakagawa et al. (1998) from the BICE survey (solid line). The modified distribution using the FILM latitudinal data to calibrate the normalized intensity at b ± 4° (dashed line).

In the text
thumbnail Fig. 2

Galactic radial emissivity of [C ii], ϵ[ CII ](Rgal), derived by Pineda et al. (2013) using the GOT C+ [C ii] survey at b = 0°.

In the text
thumbnail Fig. 3

Schematic showing two representative lines of sight (dashed lines) from the solar system through the Galactic disk (GC labels the Galactic center and b labels the latitude). The intensity of [C ii] is the integral of the emissivity, ϵ, along the path length s (see text). The actual lines of sight used to calculate the intensity cover l = 5° to 25°.

In the text
thumbnail Fig. 4

Model fit (dotted line) of the normalized intensity using the GOT C+ observations at b = 0° compared to the normalized intensity from the modified BICE + FILM profile (solid line) versus latitude b. The unmodified BICE profiles are shown for comparison (dashed line).

In the text
thumbnail Fig. 5

Gaussian fits for the normalized intensities of [C ii], CO, and H i as a function of z. The fits include the measured offsets zoffset for CO and [C ii], and an inferred offset for H i equal to that for CO. Towards the inner Galaxy CO has, on average, the narrowest distribution in the disk, followed by [C ii], and H i has the thickest distribution.

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.