{'fraction_non_alphanumeric': 0.05162573185550687, 'fraction_numerical': 0.06142062819899105, 'mean_word_length': 3.9904446894524073, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 3, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Observations of Cyanopolyynes toward Four High-Mass Star-Forming Regions Containing Hot Cores

We carried out line survey observations at the 26-30 GHz band toward the four high-mass star-forming regions containing hot cores, G10.30-0.15, G12.89+0.49, G16.86-2.16, and G28.28-0.36, with the Robert C. Byrd Green Bank Telescope. We have detected HC5N from all of the sources, and HC7N from the three sources, except for G10.30-0.15. We further conducted observations of HC5N at the 42-46 GHz and 82-103 GHz bands toward the three sources, G12.89+0.49, G16.86-2.16, and G28.28-0.36, with the Nobeyama 45 m radio telescope. The rotational lines of HC5N with the high excitation energies (Eu/k=63-100 K), which are hardly excited in the cold dark clouds, have been detected from the three sources. The rotational temperatures of HC5N are found to be 13-20 K in the three sources. The detection of the lines with the high excitation energies and the derived rotational temperatures indicate that HC5N exists in the warm gas within 0.07-0.1 pc radii around massive young stellar objects. The column densities of HC5N in the three sources are derived to be (2.0+-2.8)*10^(13) cm-2 . We compare the ratios between N(HC5N) the column density of HC5N and W(CH3OH) the integrated intensity of the thermal CH3OH emission line among the three high-mass star-forming regions. We found a possibility of the chemical differentiation in the three high-mass star-forming regions; G28.28-0.36 shows the largest N (HC5N)/W (CH3OH) ratio of>8.0*10^(14) in units of (K km s-1 )-1 cm-2 , while G12.89+0.49 and G16.86-2.16 show the smaller values (~2*10^(13) ).


INTRODUCTION
As the star formation process progresses, the chemical composition of molecular gas changes. The star formation starts with the cold dark molecular clouds (T ∼ 10 K, n ∼ 10 4 cm −3 ). In the dark clouds, carbon-chain molecules such as CCS and HC 5 N are abundant (Suzuki et al. 1992;Hirota et al. 2009). After formation of protostars, the cloud temperature rises. Such physical changes significantly affect the chemical composition, leading to hot core chemistry in the high-mass star-forming regions and hot corino chemistry in the low-mass star-forming regions. In hot core and hot corino regions (T ≥ 100 K, n ∼ 10 6 cm −3 ), while carbon-chain molecules except for HC 3 N are deficient, saturated complex organic molecules (COMs) such as CH 3 OH and CH 3 CN are rich. These COMs are formed by the combination of the grain-surface reactions and the following gas-phase reactions (e.g. Garrod & Herbst 2006;Herbst & van Dishoeck 2009).
Conversely, two low-mass star-forming cores with a variety of carbon-chain molecules have been discovered; L1527 in Taurus (Sakai et al. 2008) and IRAS 15398-3359 in Lupus (Sakai et al. 2009a). The rotational temperatures of carbon-chain molecules in L1527 were derived to be ∼ 12 − 15 K (Sakai et al. 2009b), which are higher than those in the cold dark clouds. Aikawa et al. (2008) and Hassel et al. (2008) showed that CH 4 , which is evaporated from grain mantles in the warm regions (T ∼ 20 K) around the protostar, could form carbon chains in the gas phase. In fact, carbon-chain molecules were found to be enhanced in the warm regions where the temperature becomes higher than 20−30 K (Sakai et al. 2010). Such a chemical mechanism was named warm carbon chain chemistry (WCCC, Sakai et al. 2008).
While there have been many studies about carbon-chain molecules in the low-mass star-forming regions, our knowledge of carbon-chain molecules in the high-mass star-forming regions is still poor. As with the low-mass star-forming cores, understanding chemical mechanisms of carbon-chain molecules in the high-mass star-forming cores is one of the keys to revealing the initial conditions of massive star formation. Cyanoacetylene (HC 3 N), the shortest member of cyanopolyynes, is one of the hot core tracers (a dense gas tracer; n H2 ∼ 10 5 − 10 6 cm −3 , Bergin et al. 1996), whereas longer cyanopolyynes are rarely detected. Chapman et al. (2009) demonstrated their chemical model calculation, and suggested that cyanopolyynes (HC 2n+1 N; n = 1−4) could be efficiently formed by the neutral-neutral reactions between CN and C 2n H 2 . The latter is formed from C 2 H 2 evaporated from grain mantles in hot core regions. Taniguchi et al. (2016b) proposed that the main formation pathway of HC 3 N in the high-mass star-forming region G28.28-0.36 is the neutral-neutral reaction between C 2 H 2 and CN, which is the same reaction suggested by Chapman et al. (2009), based on their observations of the 13 C isotopic fractionation with the Nobeyama 45 m radio telescope. Motivated by the chemical model calculation of Chapman et al. (2009), Green et al. (2014) conducted survey observations of HC 5 N using the J = 12 − 11 rotational line toward 79 hot cores with the Tidbinbilla 34 m radio telescope, and detected HC 5 N toward 35 sources. They could not derive the excitation temperatures due to the detection of the single rotational line which can be excited even in the cold dark clouds. Thus, the detected emission lines may come from the cold envelopes around the protostars. We still cannot confirm that long cyanopolyynes are formed and exist in hot cores.
In this paper, we carried out observations at the 26−30 GHz band toward four high-mass star-forming regions containing hot cores, G10.30-0.15, G12.89+0.49, G16.86-2.16, and G28.28-0.36, with the Robert C. Byrd Green Bank Telescope (GBT) 1 . Two rotational transitions of HC 5 N (J = 10 − 9 and 11−10) and three transitions of HC 7 N (J = 24 − 23, 25−24, and 26−25) are in the frequency band. HC 7 N has been detected in several star-forming regions (Sakai et al. 2008;Cordiner et al. 2012;Friesen et al. 2013;Feng et al. 2015), which is considered as an indicator of cyanopolyyne-rich sources. We also observed the rotational lines of HC 5 N at the 42−46 GHz and 82−103 GHz bands toward the three sources, except for G10.30-0.15, with the Nobeyama 45 m radio telescope 2 . Several lines of HC 5 N with the high excitation energies (E u /k ∼ 63 − 100 K) are in the 82−103 GHz band. Such lines are hardly excited in the cold environments, and the detection of these lines means that HC 5 N exists in the warm regions, as discussed later. From the results with the two radio telescopes, we derive the rotational temperatures and the column densities of HC 5 N in the three sources using the rotational diagram. We confirm that HC 5 N exists in the warm gas in the high-mass star-forming regions from the rotational temperatures. The thermal emission line of CH 3 OH, which is one of the hot core tracers and representative COMs, is in the frequency band of the GBT observations. The origin of CH 3 OH is considered to be different from that of carbon-chain species. We then investigate the relationship between HC 5 N and CH 3 OH in the three high-mass star-forming regions in Section 4.2.1. We found a possibility of the chemical differentiation among the three sources.

OBSERVATIONS
We describe the observational details of each frequency band in this section. Table 1 summarizes the observing parameters of each frequency band.
The observed positions, distances, systemic velocities, and related objects (ultracompact H II regions and outflows) of our four targets are summarized in Table 2. We selected these targets from the list of sources detected in HC 5 N (Green et al. 2014) by applying the following three criteria: 1. the source declination is above -21 • , 2. the distance (D) is within 3 kpc, and 3. CH 3 CN was detected (Purcell et al. 2006).
The first criterion enables us to observe with the GBT and the Nobeyama 45 m telescope. The second criterion is required for preventing from heavy beam dilution. The last one ensures that target sources contain hot cores. CH 3 CN is one of the hot core tracers, and it is always detected in the hot gas (> 100 K, Bisschop et al. 2007). Eight sources in the list meet the three criteria. Since our purpose is to derive the accurate rotational temperatures, we chose three sources among the eight selected sources which show the highest peak intensities of HC 5 N (G12. 89+0.49, G16.86-2.16, and G28.28-0.36). In order to investigate different conditions, we selected one source which shows the low peak intensity of HC 5 N (G10.30-0.15) with the high peak intensities of CH 3 CN.

Observations with the Robert C. Byrd Green Bank Telescope (GBT)
We carried out observations with the Robert C. Byrd Green Bank Telescope (GBT) of the National Radio Astronomy Observatory in 2016 April and July. We used the Ka-band receiver to cover the 26.1−29.7 GHz range, except for 28.7−28.9 GHz. The beam size (HPBW), the aperture efficiency (η A ), and the main beam efficiency (η mb ) of the Ka-band receiver were approximately 27 ′′ , 67%, and 77%, respectively. The typical system temperatures were from 70 to 100 K in April, and from 90 to 120 K in July, depending on the weather conditions and elevations. We used the VEGAS spectrometer whose bandwidth and frequency resolution were 1080 MHz and 66 kHz, respectively. The frequency resolution of 66 kHz corresponds to the velocity resolution of 0.73 km s −1 at 27 GHz. The center frequencies of four spectral windows were set at 26.525,27.425,28.325,and 29.275 GHz,respectively. We employed the position-switching mode. The integration time was 60 seconds per on-source and off-source positions. The observed positions of the four target sources are summarized in Table 2. The off-source positions were set to be +15 ′ away in the declination.

Observations with the Nobeyama 45 m radio telescope
The observations were carried out between 2017 January and March. We observed toward the three sources G12.89+0.49, G16.86-2.16, and G28.28-0.36. The on-source and off-source positions were set at the same ones with the GBT observations. The position-switching mode was employed, and the integration time was 20 seconds per scan.
We used the Z45 receiver (Nakamura et al. 2015) for the 42−46 GHz band observations. The system temperatures of the Z45 receiver were 120−150 K, depending on the weather conditions and elevations. The main beam efficiency (η mb ) and the beam size (HPBW) at 43GHz were 71% and 37 ′′ , respectively. The TZ receiver (Nakajima et al. 2013) was used for the 82−103 GHz band observations. The system temperatures were from 120 to 200 K. The main beam efficiency and the beam size at 86 GHz were 54% and 18 ′′ , respectively.
We used the SAM45 FX-type digital correlator (Kamazaki et al. 2012) in frequency settings whose bandwidths and resolutions are 500 MHz and 122.07 kHz for observations of the Z45 receiver, and 1000 MHz and 244.14 kHz for those of the TZ receiver, respectively. The frequency resolutions correspond to the velocity resolution of ∼ 0.85 km s −1 .

Observational results with the GBT
We conducted data reduction 3 using the GBTIDL 4 . The overviews of obtained spectra toward the four sources are shown in Figures 1 and 2. The rms noise levels are 14, 10, 17, and 11 mK in T mb for G10.30-0.15, G12.89+0.49, G16.86-2.16, and G28.28-0.36, respectively. We fitted the spectra with a Gaussian profile and obtained the spectral line parameters, as summarized in Table 3. The lower detection limit is set to be 4σ.
Two and three rotational lines of HC 5 N (J = 10 − 9 and 11 − 10) and HC 7 N (J = 24 − 23, 25 − 24, and 26 − 25) are in the observed frequency range, respectively. All of the observed rotational lines of HC 5 N were detected from all of the four sources. We also detected HC 7 N from the three sources, except for G10.30-0.15. Detection of HC 7 N has been reported in several star-forming regions (Sakai et al. 2008;Cordiner et al. 2012;Friesen et al. 2013;Feng et al. 2015). The detection of HC 7 N from the three high-mass star-forming regions suggests that our target cores contain the long cyanopolyynes abundantly.
In addition, several other lines were detected. The hydrogen recombination lines are prominent only in G10.30-0.15, which seems to be caused by a few ultracompact H II (UCH II) regions associated with G10.30-0.15 (Thompson et al. 2006). Although a UCH II region is associated with G28.28-0.36 (Walsh et al. 2003), we did not detect any hydrogen recombination lines in G28.28-0.36. The rotational lines of HC 3 N were detected toward all of the four sources. HC 3 N is considered as one of the hot core tracers (a dense gas tracer), and it is usually detected in hot core regions, as mentioned in Section 1. The thermal CH 3 OH emission lines (J K = 4 0 − 3 1 E) with the low-excitation energy (E u /k = 36.3 K) were detected except for G28.28-0.36, and its lines with high-excitation energies (E u /k ∼ 120 − 340 K) were detected only toward G12.89+0.49. Metastable inversion transition lines of NH 3 with very high-excitation energies, i.e., (J, K) = (8, 8) at 26.51898 GHz (E u /k ∼ 686 K) and (J, K) = (9, 9) at 27.47794 GHz (E u /k ∼ 852 K), were detected toward G12.89+0.49 and G16.86-2.16. The results indicate that extremely hot gas is contained in these regions (Goddi et al. 2011). The H 2 CO (J Ka,Kc = 3 1,2 − 3 1,3 ) emission lines were detected toward G12.89+0.49 and G16.86-2.16.

Observational results with the Nobeyama 45 m radio telescope
We conducted data reduction using the Java Newstar, which is an astronomical data analyzing system of Nobeyama data 5 . Figures 3 and 4 show the spectra of HC 5 N at the 42−46 GHz and 82−103 GHz bands, respectively. The rms noise levels are from 6 to 14 mK in T * A at the 42−46 GHz band and from 3 to 6 mK at the 82−103 GHz band. Table  4 summarizes the line parameters obtained from the Gaussian fitting. The main beam efficiencies which we used in the analysis are summarized in Table 4.
The signal-to-noise ratios were between 4 and 23. The two rotational lines of HC 5 N at the 42−46 GHz band (J = 16 − 15 and 17−16) were detected from all of the three sources. Moreover, we detected six rotational lines of HC 5 N with the high excitation energies (E u /k ∼ 63 − 100 K) from all of the sources. The V LSR values of the detected lines are summarized in Table 4. The V LSR values of the HC 5 N lines are consistent with the systemic velocities of each source (Table 2).

Analyses
We derived the rotational temperatures and column densities of HC 5 N in the three sources (G12.89+0.49, G16.86-2.16, and G28.28-0.36) from the rotational diagram analysis, using the following formula (Goldsmith & Langer 1999); where k is the Boltzmann constant, S is the line strength, µ is the permanent electric dipole moment, N is the column density, and Q(T rot ) is the partition function. We used 4.330 D (Alexander et al. 1976) for µ of HC 5 N. Figure 5 shows the fitting results of the rotational diagram for the three sources. The errors include the Gaussian fitting errors, the uncertainties from the main beam efficiency of 10%, the chopper-wheel method of 10%, and from the other factors such as calibration and filling factor of 30%. Since we do not know the spatial distributions of HC 5 N, we assume the filling factor of unity. In the fitting, the values of the J = 16 − 15 and 17−16 transitions are systematically lower than the fitting lines. This suggests a smaller filling factor than unity, because the beam size at the 42−46 GHz band with the Nobeyama 45 m telescope is significantly larger than those at the 26−30 GHz band with the GBT and at the 82−103 GHz band with the Nobeyama 45 m telescope. The ranges of the derived rotational temperatures and column densities in the three sources are 17−23 K and (7.5−11)×10 12 cm −2 , as summarized in Table 5. In the case of G10.30-0.15, we detected only two lines with the similar energies using the GBT, and we cannot conduct the rotational diagram analysis. We derived the column density of HC 5 N in G10.30-0.15 to be 8.5 +2.0 −1.7 × 10 12 cm −2 , using the average value of the rotational temperatures in the three sources (21 ± 11 K). The error of the rotational temperature in G10.30-0.15 27.29429 2.6 1.14 (5) 7.9 (4) 9.6 (6) 1.60 (11)     was estimated from the relative errors of the rotational temperatures in the three sources, and the error of the column density corresponds to the change in the rotational temperature. Assuming a small beam filling factor, we multiplied the integrated intensities of the GBT data by ( 27 ′′ 18 ′′ ) 2 and the 45GHz band data by ( 37 ′′ 18 ′′ ) 2 for the correction of the different beam sizes. Figure 6 shows the rotational diagram with the beam-size correction. All of the data, including the 42−46 GHz band, are better fitted. From the above  considerations of beam sizes and a small beam filling factor, the emission region sizes are limited to below 18 ′′ , which correspond to 0.07−0.1 pc radii at the target source distances (1.67−3.0 kpc, Table 2). Table 5 summarizes the derived values of the rotational temperatures and column densities. The rotational temperatures are derived to be 13.8−18 K, and the corresponding column densities are (2.05−2.78) ×10 13 cm −2 . The derived rotational temperatures with and without the beam-size correction are consistent within their 1σ errors. The rotational temperature and the column density in G10.30-0.15 were derived by the same manner as mentioned in the previous paragraph (T rot = 16 ± 4 K, N = 7.5 +0.8 −0.6 × 10 12 cm −2 ). We use the results with the beam-size correction in the following sections, because the fitting results are better. We compare the results in the high-mass star-forming regions with that in the cyanopolyyne peak in Taurus Molecular Cloud-1 (TMC-1 CP, d = 140 pc), which is a well studied dark cloud and carbon-chain-species-rich source. The derived excitation temperature of HC 5 N in TMC-1 CP is 6.5 ± 0.2 K (Taniguchi et al. 2016a). The rotational temperatures of HC 5 N in our three target high-mass star-forming regions except for G10.30-0.15, ∼ 13 − 20 K, are significantly higher than the excitation temperature in the dark cloud TMC-1 CP.
The rotational temperature of HC 5 N in the low-mass star-forming region L1527 (d = 140 pc) is known to be 14.7±5.3 K (3σ, Sakai et al. 2009b). The rotational temperatures of HC 5 N in our target high-mass star-forming regions are comparable with or slightly higher than that in L1527. L1527 is one of the WCCC sources, in which carbon-chain  species are formed in the warm gas around the protostar with the temperature of ∼ 20 − 30 K (Sakai & Yamamoto 2013). Therefore, HC 5 N seems to exist in the warm gas in our target high-mass star-forming regions. The detection of the emission lines with the high excitation energies also allows us to confirm the warm components of HC 5 N by itself. The detected emission lines of HC 5 N with the high-excitation energies (E u /k ∼ 90 − 100 K) should not come from the cold environments with the gas kinetic temperature of 10 K. We assume that the column density of HC 5 N is comparable with that in TMC-1 CP (N (HC 5 N) = (6.2 ± 0.3) ×10 13 cm −2 , Taniguchi et al. 2016a) and the rotational temperature is a typical gas kinetic temperature in cold dark cloud (10 K). The peak intensities of the J = 38 − 37 and 39−38 lines are calculated to be 0.8 mK in T * A , even though the column density of the particular carbon-chain-rich source is considered. Thus, the detection of the emission lines of HC 5 N with the energies of ∼ 90−100 K implies that HC 5 N exists in the warm gas. There is a possibility that the emission lines with the energies of < 20 K come from not only the warm gas but also the cold gas. In that case, the derived rotational temperatures will be the lower limits due to the mixing of the warm gas and the cold gas in the single-dish beams. To summarize, we can conclude that HC 5 N exists in the warm gas within 0.07−0.1 pc radii around the massive young stellar objects. We investigate the relationship between HC 5 N and CH 3 OH in the high-mass star-forming cores. CH 3 OH is one of the representative hot core tracers and COMs, and its origin is considered to be different from carbon-chain molecules. Figure 7 shows the comparison of the ratios between N (HC 5 N) the column density of HC 5 N and W (CH 3 OH) the integrated intensity of the thermal CH 3 OH emission line summarized in Table 3. We cannot derive the excitation temperature and column density of CH 3 OH due to the detection of only single line with low-excitation energy (E u < 50 K). The other detected lines with high-excitation energies in G12.89+0.49 could possibly be masers (Müller et al. 2004). Therefore, we cannot use these lines to derive the excitation temperatures.
Since the column density of HC 5 N in G10.30-0.15 was derived using its average rotational temperature of the other three high-mass star-forming regions, we do not discuss the result of G10.30-0.15 in detail due to its large uncertainty. We can divide three high-mass star-forming regions into two groups; G28.28-0.36 shows the largest N (HC 5 N)/W (CH 3 OH) ratio of > 8.0 × 10 14 in units of (K km s −1 ) −1 cm −2 , while those of G12.89+0.49 and G16.86-2.16 are only (1.9 ± 0.3) ×10 13 and (2.6 ± 0.5) ×10 13 , respectively. The ratio in G28.28-0.36 is significantly larger than those in G12.89+0.49 and G16.86-2.16 by one order of magnitude, even though there is a possibility that the chemical and physical structures in the beams cause the observed large difference. The N (HC 5 N) values are almost comparable among the three sources. Therefore, the differences in the N (HC 5 N)/W (CH 3 OH) ratio come from the W (CH 3 OH) values.
The integrated intensity is proportional to both the column density and the rotational temperature. When we assume the same N (HC 5 N)/N (CH 3 OH) ratios in all of the three sources, the differences in the N (HC 5 N)/W (CH 3 OH) ratios depend on only the rotational temperature of CH 3 OH. If the rotational temperature of CH 3 OH in G28.28-0.36 is the typical value in hot cores (100−250 K, Bisschop et al. 2007), the large difference in the integrated intensities of CH 3 OH by an order of magnitude (e.g., for the case of G28.28-0.36 and G16.86-2.16) cannot be explained given any rotational temperature assumed in G16.86-2.16. Alternatively, if the rotational temperatures of CH 3 OH in G12.89+0.49 and G16.86-2.16 are 100−250 K, G28.28-0.36 will have an unlikely high rotational temperature of CH 3 OH, 1200−2600 K. Thus, the differences in the N (HC 5 N)/W (CH 3 OH) ratio among the three sources cannot be explained only by changes in the rotational temperature of CH 3 OH, but the column density of CH 3 OH in G28.28-0.36 should be lower than those in G12.89+0.49 and G16.86-2.16. In summary, the observed variation in the N (HC 5 N)/W (CH 3 OH) ratio reflects the chemical differentiation in the single-dish beams.

Fractional abundances of HC5N in the four high-mass star-forming regions
We also derived the fractional abundances of HC 5 N defined as X(HC 5 N) = N (HC 5 N)/N (H 2 ) toward the four sources. We calculated the column densities of H 2 from the SCUBA, installed on the James Clerk Maxwell Telescope, 850 µm continuum data using the following formula (Kauffmann 2005 The continuum data were obtained from the Canadian Astronomy Data Center, JCMT Science Archive 6 . We used the continuum fluxes (S beam ν ) toward the observed positions and the beam size (18 ′′ ) as summarized in Table 6 7 . We used 0.0182 cm 2 g −1 for κ ν (Kauffmann 2005), and 18 ′′ for θ HPBW . T is the dust temperature obtained by SED fitting (Purcell et al. 2009), as summarized in Table 6. In G28.28-0.36, the dust temperature was not derived, and we used the average value of 58 K, determined from 65 good fits (Purcell et al. 2009). Table 6 summarizes the column densities of molecular hydrogen N (H 2 ) and the X(HC 5 N) values. The errors of N (H 2 ) are derived assuming that uncertainty of the dust temperatures is 20 K. Figure 8 shows the comparison of the X(HC 5 N) values among the four sources. The X(HC 5 N) values show the same tendency as the N (HC 5 N)/W (CH 3 OH) ratios. G28.28-0.36 shows the larges X(HC 5 N) value (4.2 +2.6 −2.0 × 10 −9 ) in the four sources, and G12.89+0.49 and G16.86-2.16 show the similar values, which is smaller than that in G28.28-0.36. As mentioned before, the column density of HC 5 N in G10.30-0.15 was derived using the average rotational temperature of the three sources, and so we do not discuss in detail. Note-The errors represent the standard deviation.
a Taken from online material of Purcell et al. (2009).
b The errors were derived from the uncertainties of the dust temperatures of 20 K. c The average value derived from 65 good fits (Purcell et al. 2009).
Both Figures 7 and 8 imply a possible signature of chemical differentiation among the high-mass star-forming cores, although we selected only four biased samples associated with HC 5 N emission and the 6.7 GHz methanol masers. We cannot exclude a possibility of the chemical and physical structures in the beam, and one possible explanation for such chemical differentiation is the different beam filling factor of hot core regions (Hatchell et al. 1998). Observations to reveal the spatial distributions of carbon-chain molecules and COMs with the interferometry are essential for revealing the origin of the chemical differentiation.
The fractional abundance of HC 5 N in L1527 was derived to be ∼ 2 × 10 −10 (N (HC 5 N) = (6.8 ± 1.4)×10 12 cm −2 and N (H 2 ) = 2.8 × 10 22 cm −2 were derived by Sakai et al. (2009b) and Jørgensen et al. (2002), respectively). Even though there are uncertainties in column densities of H 2 and HC 5 N due to assumed parameters (e.g. temperature, dust opacity, source size), the derive fractional abundances of HC 5 N in the three high-mass star-forming regions (except for G10.30-0.15) are comparable to or slightly higher than that in L1527 by a factor of 5−20. Such high HC 5 N abundances in our target sources suggest that the warm gas around these massive young stellar objects is rich in HC 5 N. Moreover, the high HC 5 N abundances would suggest efficient formation mechanisms of carbon chains in these high-mass star-forming regions (Aikawa et al. (2008) and Hassel et al. (2008)  (∼ 2 × 10 13 ). Such chemical differentiation may originate from the complex chemical and physical structures in the beam.
We deeply appreciate the anonymous referee for constructive suggestions, which are very useful to develop our discussion. We would like to express our great thanks to the staff members of the National Radio Astronomy Observatory. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. K. T. deeply appreciates Dr. Nichol Cunningham (NRAO) for her kind help to prepare and carry out observations. We would like to express our thanks to the staff members of the Nobeyama Radio Observatory. The Nobeyama Radio Observatory is a branch of the National Astronomical Observatory of Japan, National Institutes of Natural Sciences.
Facilities: Robert C. Byrd Green Bank Telescope, Nobeyama 45 m radio telescope Software: GBTIDL, Java Newstar