SCUBA chop throw dependence
The dependence of SCUBA sensitivity and image shape on chop throw |
|
Introduction
SCUBA's flux sensitivity, measured by the Flux Calibration Factor (FCF, see the Calibration Cookbook) and image characteristics both depend on the chop throw used for the observation. During every observation made with SCUBA, the secondary mirror unit (SMU) chops between the "on" position and an "off" position at a constant rate of 7.8125 Hz. Since this is a real physical movement, there is some portion of every cycle during which the SMU is at neither position but is moving between the two. This leads to a smearing out of the image along the chop axis. The effect is strongest with the largest chop throws since the constant rate means that less time per cycle is spent in the "on" and "off" positions. This leads to lower flux in the "on" position, and therefore higher FCFs. This can be seen in the sequence of images below, which shows jiggle-maps of CRL618 with three different chop throws:
| 850 micron (50, 120 and 180" chops respectively, contours same in each image at 0.002 (white - may be difficult to see), 0.0015 0.001, 0.0009, 0.0008, 0.0007, 0.0006, 0.0005, 0.0004, 0.0003 units): |
 |
 |
 |
(Note about blanking times: When performing beamswitching heterodyne observations with the JCMT, a "blanking time", dependent on the chop throw, is applied to the signal. This effectively cuts out the signal from the periods where the SMU is known to be travelling between "on" and "off" positions. Plans were made to do the same for SCUBA observations, however when starting to probe this issue we discovered that the blanking time was in fact never activated for SCUBA. As a result the entire signal, incuding travelling time, is included in a SCUBA image.)
Hence we strongly recommend that you make calibration observations using the same chop parameters as your science observations! However, there are times when this is not possible and so we have attempted to evaluate the relationships between FCFs at different chop throws. At 850 microns this turns out to be a relatively simple task; at 450 microns the data is not so friendly! Also below is a brief discussion of the image shape characteristics at different chop throws.
850 micron FCF relations
The FCF relations at 850 microns are relatively easy to pick out. We reduced all Uranus jiggle-map observations in the archive that were a) made in good observing conditions, b) made between 8pm and 4am, to avoid thermal effects on the dish, c) made during times when we know SCUBA's performance was not atypical, using the monthly FCF summaries, d) non-polarimetric observations.
The table below summarises these results (note that the aperture FCFs are calculated on the basis of our standard 40" aperture):
| Chop throw (") |
Number of data points |
Beam FCF (Jy/beam/volt) |
Aperture FCF (Jy/square arcsec/volt) |
| 30 |
34 |
196 ± 8 |
0.878 ± 0.03 |
| 35 |
6 |
198 ± 10 |
0.864 ± 0.04 |
| 44 |
20 |
204 ± 12 |
0.856 ± 0.04 |
| 45 |
6 |
201 ± 10 |
0.850 ± 0.04 |
| 60 |
259 |
210 ± 26 |
0.859 ± 0.11 |
| 68 |
16 |
206 ± 7 |
0.839 ± 0.04 |
| 90 |
5 |
212 ± 6 |
0.859 ± 0.03 |
| 120 |
73 |
227 ± 15 |
0.897 ± 0.05 |
| 150 |
4 |
240 ± 17 |
0.951 ± 0.05 |
| 180 |
8 |
255 ± 21 |
0.975 ± 0.07 |
Some notes on these results:
- The data were reduced using the in-house version of ORACDR, using CSO tau fits for opacity correction and FLUXES for FCF value calculation.
- The aperture flux calculation uses a 40" diameter aperture only if this is greater than the chop throw; otherwise, an aperture equal to the chop throw is used. Hence the 30" and 35" aperture FCFs are calculated using smaller apertures than the remainder of the data set. This effect, in combination with the interference of the negative beams within the aperture, causes the upturn in aperture FCFs with low chop throws.
The trends shown in the above table are plotted below, along with lines of best fit. For the Beam FCF dependence, a quadratic fit to all data points is given, whereas for the Aperture FCF dependence, the fit was constrained only for data points ≥ 60" (see point 2 above). The error bars shown in these plots are reduced by the factor of √N from the values in the table above.
| Trend |
Beam FCF |
Aperture FCF |
| |
|
|
| Best fit formula (dashed line, chop in ") |
0.001*(chop)^2 + 0.17*chop + 191 |
0.000004*(chop)^2 - 0.0001*chop + 0.83 |
450 micron FCF relations
When we tried to repeat the same exercise with 450 micron data, we found that, even taking the steps described above to minimise effects such as dish temperature and known SCUBA problems, we were unable to find smooth trends in the averaged data. 450 micron observations with SCUBA are extremely sensitive to small changes in focus and opacity as well as to the night-to-night state of SCUBA, and so these sources of noise in the datasets drown out the trend due to changing chop throw.
To minimise the night-to-night changes, we instead chose to only make intra-night comparisons to produce typical ratios of FCFs, rather than worry about the absolute values. Unfortunately, nights when more than three different chop throws are used are rare. However, by combining all the useful data for Uranus and CRL618 (Mars was not used since it can become too large) we were able to get some consistency. These results are summarised in the table below. Note that the ratios are always expressed in the order shown i.e. FCF(smaller chop throw)/FCF (larger chop throw).
| Chop throw ratio |
Number of data points |
Beam FCF average ratio |
Aperture FCF average ratio |
| 30/35 |
1 |
0.964 |
1.004 |
| 30/40 |
1 |
1.102 |
1.202 |
| 30/44 |
51 |
1.096 ± 0.012 |
1.092 ± 0.008 |
| 30/45 |
1 |
1.095 |
1.133 |
| 30/60 |
4 |
1.000 ± 0.061 |
1.124 ± 0.030 |
| 30/68 |
52 |
1.042 ± 0.016 |
1.049 ± 0.013 |
| 30/90 |
2 |
1.065 ± 0.048 |
1.101 ± 0.033 |
| 30/100 |
3 |
1.019 ± 0.036 |
1.058 ± 0.003 |
| 30/120 |
15 |
0.890 ± 0.017 |
0.973 ± 0.015 |
| 30/150 |
2 |
0.726 ± 0.108 |
0.829 ± 0.083 |
| 30/180 |
4 |
0.985 ± 0.062 |
1.107 ± 0.062 |
| 35/100 |
1 |
0.899 |
1.010 |
| 35/120 |
3 |
0.785 ± 0.044 |
0.905 ± 0.038 |
| 40/50 |
1 |
0.982 |
1.014 |
| 44/68 |
54 |
0.960 ± 0.009 |
0.977 ± 0.007 |
| 44/90 |
1 |
1.052 |
0.962 |
| 44/100 |
2 |
0.931 ± 0.006 |
0.948 ± 0.003 |
| 44/120 |
9 |
0.839 ± 0.036 |
0.862 ± 0.040 |
| 44/180 |
1 |
0.626 |
0.847 |
| 45/60 |
2 |
1.068 ± 0.045 |
0.934 ± 0.024 |
| 45/68 |
1 |
0.974 |
1.029 |
| 45/120 |
5 |
0.813 ± 0.019 |
0.858 ± 0.021 |
| 45/150 |
1 |
0.876 |
0.812 |
| 50/60 |
1 |
0.989 |
0.976 |
| 60/90 |
2 |
0.928 ± 0.002 |
0.948 ± 0.012 |
| 60/100 |
1 |
1.029 |
0.964 |
| 60/120 |
3 |
0.889 ± 0.110 |
0.908 ± 0.016 |
| 60/180 |
1 |
0.630 |
0.764 |
| 68/90 |
1 |
1.095 |
1.004 |
| 68/100 |
3 |
0.964 ± 0.002 |
0.946 ± 0.030 |
| 68/120 |
8 |
0.901 ± 0.024 |
0.963 ± 0.050 |
| 90/120 |
4 |
0.943 ± 0.017 |
0.938 ± 0.010 |
| 90/180 |
2 |
0.711 ± 0.596 |
0.838 ± 0.647 |
| 100/120 |
3 |
0.889 ± 0.043 |
0.973 ± 0.052 |
| 120/150 |
2 |
0.724 ± 0.602 |
0.865 ± 0.658 |
| 120/160 |
1 |
0.791 |
0.781 |
| 120/180 |
6 |
0.735 ± 0.021 |
0.915 ± 0.023 |
| 150/180 |
1 |
0.945 |
1.051 |
As is obvious from the above table, there is a wide variety in the number of data points available for each chop throw ratio. The most represented chop throws are 30", 44" and 68", and these often come together, since these are the three chop throws used for SHADES observations. Other areas of "chop throw space" are markedly underrepresented including, unfortunately, > 120", which is the region where the largest changes are expected.
For the five most represented chop throws in the above table, we have produced plots to show the trend of their ratios. In the plots below, the constant chop throw is always regarded as the numerator; e.g. most of the ratios involving 120" have been inverted for inclusion in the 120" plot. The aperture ratios are shown by a dotted circle whilst the beam ratios are shown by a cross. Only those data points with errors ≤ 0.1 have been included.
| 30"/x |
44"/x |
68"/x |
100"/x |
120"/x |
 |
 |
 |
 |
 |
Some notes on these results:
- There is some degree of circularity in these results and plots, which was inevitable given the small data set. This is especially the case with the plotted data, so the trend seen in one plot may be similar to that in another because they are using the same results - i.e. beware of overinterpreting these data!
- It seems unlikely that the apparent upturn in the 30" ratio data set for 180" is real - a downwards trend with increasing denominator makes more sense and is seen in every other plot.
- We have not attempted to fit to these data plots because of the high degree of circularity in the results, especially the 44, 68 and 100" ratio data sets.
- The error on each ratio is typically lower for the Aperture FCF than for the Beam FCF. This is as expected since the Aperture FCF is least sensitive to changes in focus. Hence, we would recommend that, if you are going to use the ratios presented above to make an FCF conversion, you use the Aperture calibration technique whenever possible (see the Calibration Cookbook for further discussion).
Image characteristics
The change in FCF with chop throw reflects a change in the image shape, since the actual signal is conserved across the whole image. We have evaluated this effect using the KAAPA command PSF, which fits a Sersic profile, with D = A exp (-0.5 * [r/σ]^ γ), where r is the radial distance from the fitted centre of the source, σ is the profile width and γ is the radial fall-off parameter. If γ = 2, the fitted profile is a perfect Gaussian. The results of running PSF on the above data sets are shown below.
Notes:
- There seem to be quite clear trends in the variation of FWHM at both wavelengths, with a relatively large increase over the lower chop throws, then a settling down to a fairly consistent value of ~ 14.8"e; at 850 µm and ~ 7.8" at 450 µm. The largest increase takes place over the range ~ 30 - 70" at at 850 µm and over the range ~ 45 - 100" at 450 µm.
- The trends are not quite so consistent in the variation of gamma. At 850 µm the trend appears to be relatively flat at around γ ~ 2.16, with two much lower values at the lowest chop throws, 30 and 40". At 450 µm, there is a very clear turnover from a chop throw of ~ 60", where γ ~ 2.08, to a chop throw of 180", where γ ~ 1.86. This is despite the apparent consistency in the FWHM between these chop throws. Perhaps this indicates that at high chop throws, the combination of the off-source smeared out portion of the data and the on-source beam-shaped portion looks most like two separate but combined components, whereas at lower chop throws the smearing is lower and the appearance (to PSF at least) is more of a single component.
|