next previous
Up: Molecular gas in

3. Observations and data reduction

 

The observations were carried out between the beginning of August 1993 and the end of August 1994. On Cerro Tololo, the Galactic center is above an elevation of 20tex2html_wrap5415 for 11 hours per day. Because the telescope cannot track sources at elevations higher than 78tex2html_wrap5415, the total daily observing time was 9 hours.

The tex2html_wrap_inline5241 line was placed in the lower sideband of the receiver because the line frequency is close to the lower limit of the tuning range. With this, the image frequency of the system was 112.556165 GHz. The double-sideband receiver temperature, tex2html_wrap_inline5433, determined with a standard "Hot-Cold-Test'' using the filterbank backend, was always in the range of 180 to 220 K. The radial velocity of the center channel of the AOS was set to 0  with respect to the Local Standard of Rest (LSR).

A detailed description of the observing procedure is given in Appendix A (click here). This includes the scan integration with weighted-OFFs position switching, the methods of scan calibration, and the telescope pointing.

3.1. Survey plan

Because the FWHM beam width of the telescope is 9tex2html_wrap54172 at the tex2html_wrap_inline5241 line frequency, a fully sampled map requires data spaced at 4tex2html_wrap54176. As a compromise to cover a larger region of the Galactic center, we have used a grid spacing of 9tex2html_wrap5305 (= 0tex2html_wrap538115), that is half-sampling. We chose the grid axes to follow Galactic coordinates centered on (tex2html_wrap_inline5449). Since the tex2html_wrap_inline5225 emission is weak and the AOS was available only for a limited time, we decided to concentrate on the inner region of the Galactic center and to extend the map only to the Clump 2 region (Bania 1977) in the range of l = +3tex2html_wrap53810 to +3tex2html_wrap53815.

3.2. Data acquisition

Each position was observed with a series of 10-minute-scans. Typically, a single scan had an rms noise temperature per channel of about 0.12 to 0.15 K, with extremes from 0.08 to 0.25 K, on a scale, as obtained with the AOS (see Sect. 3.4.2 (click here)). The scans were averaged until the rms noise temperature per channel in the average spectrum was less than 0.05 K. Typically, twelve 10-minute-scans per position were needed to reach this limit, and on average 0.04 K was reached. Thus, two to three positions per day could be completely measured.

This long total integration time was necessary to reach a reasonable signal-to-noise ratio in the spectra. The most intense tex2html_wrap_inline5241 emission, at SgrB2, is only 0.45 K, i.e. for this source, the maximum signal-to-noise ratio per channel was about 10. Away from this peak the emission decreases quickly to about 0.2 K peak temperature and outside the inner region to below the tex2html_wrap_inline5473 detection limit of 0.15 K. At many grid positions tex2html_wrap_inline5241 emission is only detected by integrating over a number of velocity channels.

The OFF positions used, verified to be free of tex2html_wrap_inline5237 emission stronger than 0.04 K, are given in Appendix B (click here). Since these OFF positions are appropriate for tex2html_wrap_inline5237 observations, they are very well suited for measurements of the much weaker tex2html_wrap_inline5241 emission.

3.3. Data reduction

 

Each spectrum was examined individually, and bad spectrometer channels, difficult to avoid on a CCD chip with 1499 channels, were removed and replaced by interpolating the neighbouring channels. After this cleaning, all spectra of the same position were averaged, the regions for the baseline subtraction were selected and a polynomial baseline was subtracted. By inspecting every single spectrum individually we checked for artefacts which might mimic emission in the averaged spectrum. Such features occur rarely. These were not considered to be emission but rather as instrumental baseline effects contributing to the uncertainties.

  figure660
Figure 2: The averaged spectrum of 13 individual 10-minute-spectra toward the position l = 0tex2html_wrap538175, b = -0tex2html_wrap538115 of the tex2html_wrap_inline5241 Galactic Center Survey. The line windows and the tex2html_wrap_inline5493 order baseline are also displayed. The rms of the fit of the baseline to the spectrum is 0.047 K

The reduction procedure was as follows: If the noise in all spectra was roughly the same, as was the case for most positions, the spectra were first averaged, weighted by time, i.e. equally. Then, the regions for the baseline subtraction were selected. The average spectrum was used for this selection, because here weak line emission could be recognized more easily and, therefore, the window regions, where the channel content could not be used for the baseline fit, could be set accordingly. In cases of uncertainty these windows were determined by comparison with tex2html_wrap_inline5237 spectra. If the noise of individual spectra for one position was significantly different, first, the line windows were set in a provisional average of the spectra. Then, baselines were subtracted from the individual spectra and these spectra were averaged, weighted by the rms per channel, determined from the baseline fit. Finally, a baseline of order 0 was subtracted from the average spectrum. We have checked that averaging weighted by rms did not have any effect on the intensity and line shape compared to averaging weighted by time. Such an average yielded lower rms per channel only if the noise of the single spectra contributing to one position was significantly different.

The polynomial baseline subtracted was typically of order 6 with long "periods'' (tex2html_wrap_inline5497 ) and low amplitudes (tex2html_wrap_inline5499). This baseline was caused by standing waves from reflections in the telescope. The FWZI extent of the tex2html_wrap_inline5225 emission is typically 50 to 250 , which is less than in tex2html_wrap_inline5325. In addition, the FWHP width of individual tex2html_wrap_inline5225 emission features are typically about 30 to 60  which is one order of magnitude less than the period of the subtracted baseline. Therefore, we found that the baselines and tex2html_wrap_inline5225 emission were clearly distinguishable and we could ascertain that the tex2html_wrap_inline5493 order baselines do not affect the intensities of the 18 emission systematically nor introduce artefacts. In Fig. 2 (click here), we show a typical example of an averaged spectrum with the chosen emission regions and the subtracted baseline of tex2html_wrap_inline5493 order.

  figure678
Figure 3: The integrated intensity of the innermost tex2html_wrap_inline5513- and b = -0tex2html_wrap538115-positions of the Galactic center region which were measured twice. The open squares represent the results of August 1993, the crosses the results of March 1994. The positions at b = -0tex2html_wrap538115 are plotted with their l value reduced by 0tex2html_wrap538105

To check that the observations were reproduceable the innermost positions of the Galactic center region (l = -0tex2html_wrap538115 to 1tex2html_wrap538135, b = -0tex2html_wrap538115 to 0tex2html_wrap53810), which were first observed in August 1993, were remeasured in March 1994. In Fig. 3 (click here), the integrated intensities, tex2html_wrap_inline5539, from both observing periods are compared for the velocity ranges from -373.7 to -225.0, from -225.0 to +225.0, and from +225.0 to +373.7  which represent the emission of the tex2html_wrap_inline5371) line, the emission of the tex2html_wrap_inline5241 line, and an emission-free region, respectively. These plots show that the observations were reproducible within a certain scatter and that no time-dependent systematic changes were present in the data. A good measure for the scatter is the mean deviation of the integrated intensity of the two periods which is given by:
equation698
where Ii is the integrated intensity of position i. The average value of this mean deviation, tex2html_wrap_inline5561, is given in Table 1 (click here) and compared to the expected value, tex2html_wrap_inline5563, from the average rms per channel, 0.04 K, determined by:
 equation708
where tex2html_wrap_inline5565 is the number of channels covered by the integration area and tex2html_wrap_inline5567 is the channel width. From this comparison, it becomes clear that in the emission-free area the scatter is within the expected limit whereas in the areas containing emission it is larger. This indicates that the uncertainties of the integrated intensities in areas of emission is mostly caused by uncertainties in the baselines and that the rms noise per channel plays a minor role. The reason for the rather strong influence of the baseline choice on the integrated intensity is that the observed emission lines are rather weak but have extended wings. Therefore, a small error in the baseline causes a relatively large error in the integrated intensity. This effect is even larger for the tex2html_wrap_inline5371) line due to additional uncertainties caused by this line being at the edge of the frontend and backend bands. It illuminates the necessity of having not only a sensitive system but also a stable system with flat baselines when broad lines are observed (see also Radford et al. 1996 for an analysis of baseline problems when observing broad, in their case extragalactic, lines).

   

Velocity range Line tex2html_wrap_inline5571 tex2html_wrap_inline5573
K K
-373.7 - -225.0 HNCO 7.10 0.84
-225.0 - +225.0 C18O 5.42 1.45
+225.0 - +373.7 -- 0.25 0.84
Table 1: The averaged relative deviation of the integrated intensity of the innermost Galactic center region of the two observation shifts in August 1993 and March 1994 and the expected deviation

3.4. Data calibration

 

The calibration of the scans follows the standard chopper wheel method first described by Penzias & Burrus (1973) and applied at many other millimeter telescopes. See Appendix A.2 (click here) for a detailed description of this method.

3.4.1. Calibration stability

In addition to daily antenna tippings, the calibration stability was monitored by taking a 10-minute-spectrum toward the position l = 0tex2html_wrap5381625, b = 0tex2html_wrap53810 (SgrB2) once or twice per day. With a peak temperature of 1.14 K, the tex2html_wrap_inline5371) line is considerably stronger than the tex2html_wrap_inline5241 line (0.45 K) toward this position and its intensity, integrated over a velocity range from -322.0 to -245.0  (with respect to the tex2html_wrap_inline5241 line frequency) is a good indicator for the calibration stability.

The stability of the calibration, throughout the course of the survey, proved to be excellent (better than 1% on different days), and the mean uncertainty in the calibration of a single scan is tex2html_wrap_inline5629. For 12 scans, the typical number of scans per position, the calibration uncertainty just from statistics is tex2html_wrap_inline5631 and probably better, because a large part of the calibration uncertainty is caused by the uncertainty in the determination of the baseline which was much better determined in the averaged spectrum (see also Sect. 3.3 (click here)).

3.4.2. Calibration scaling of the AOS spectra

 

The chopper wheel method calibrates the signal coming out of the backend to the antenna temperature, . This calibration bases on the scaling factor tex2html_wrap_inline5637 (see Appendix A.2 (click here)), which, for instrumental reasons, could be determined only with the filterbank. Because the signal path to the AOS differed from that to the filterbank, the response of the two backends was not equal. Therefore, as a last step in the calibration process, the correct scaling of the spectra taken with the AOS to , the AOS "efficiency'', had to be determined. The effect of different response could be seen most clearly toward positions where the tex2html_wrap_inline5241 line emission was narrow enough and close enough to = 0  to be also observed with the filterbank. A comparison of more than 150 spectra, taken simultaneously with the AOS and the filterbank toward 13 different positions, yielded a ratio of peak antenna temperatures between the AOS and the filterbank, resampled to the resolution of the AOS, of tex2html_wrap_inline5641, and a ratio of tex2html_wrap_inline5643 for integrated intensities (see Dahmen 1995 for more details).

However, to finally determine the different scaling factor of the AOS it is necessary to compare an emission feature of the AOS, which is well resolved and which lies well above the detection limit, with the result of an observation toward the same position with the same beamwidth and with comparable frequency resolution the calibration of which is well known. Obviously, SgrB2 (l=0tex2html_wrap5381625, b=0tex2html_wrap53810) is a suitable position for this purpose. SgrB2 was observed in tex2html_wrap_inline5241 with the 0.5 MHz-resolution-filterbank at the 1.2m NMWT during November 1994. This telescope has the same beamwidth as the SMWT and the calibration is well established (Cohen et al. 1986). Because the scaling factor from to , tex2html_wrap_inline5655, is the same for both telescopes (see Appendix A.2 (click here)) the values could be compared directly. This comparison was done using peak intensities instead of integrated intensities because a significant part of the integrated intensity originates in the extended wings of the broad line. Thus, small errors in the baseline give rise to larger errors in the integrated intensity, since the integrated intensity is the product of the velocity window and the average intensity in this window. An identical absolute baseline error is translated to a larger relative error if related to the average intensity than if related to the peak intensity. In both spectra, the peak intensity was fitted, with a Gaussian curve. With these fits, the "efficiency'' of the AOS calibration at the tex2html_wrap_inline5241 line emission, tex2html_wrap_inline5659, is:
equation774
This value is, within the noise, the same as the one determined by comparing spectra obtained with the AOS with spectra obtained with the FB at the 1.2m SMWT. Thus, the main (and probably single) cause for scaling differences between spectra taken with the SMWT and the 1.2m NMWT is the integration of the AOS into the system of the SMWT. The calibration of the SMWT filterbank is identical to within the noise with the calibration obtained at the NMWT, and the following final conversions for and can be established:
 eqnarray779


next previous
Up: Molecular gas in

Copyright by the European Southern Observatory (ESO)
web@ed-phys.fr