====== Calibration ======
For heterodyne observations, there are three types of calibrations: temperature (chopper) calibration, telescope efficiency, and (if necessary) AOS frequency calibration.
It is best to calibrate regularly.
* [[#Temperature Calibration]] determines //T//A*, the antenna temperature corrected for atmospheric absorption and ambient temperature telescope losses, and measures the system temperature. \\ ''__UIP__> **cal**''
* [[#Telescope efficiency]]. After pointing, measure the brightness of planets to determine the telescope efficiency. \\ ''__UIP__> **planet //planet// **''
* Then with beam switching (secondary mirror wobbling) \\ ''__UIP__> **chop 2**''
* or with position switching using an ''//offset//'' in arcsec \\ ''__UIP__> **oo 1 /step //offset// **'' \\
* [[#AOS Frequency Calibration|Frequency calibration]] is necessary for [[aos5|AOS5]]. It is done automatically each time AOS5 is restarted. Frequency calibration is not necessary for the FFT spectrometers.
//**REMEMBER**// Observers must decide how important the calibration is to them and how to deal with it.
Several authors (Tom Phillips, Ken [Taco] Young, and Todd Groesbeck) contributed to this description of the calibration procedure.
===== Temperature Calibration =====
Heterodyne observations are calibrated with the chopper method first discussed for millimeter astronomy by Penzias and Burrus (1973 ARA&A 11, 51 [[http://adsabs.harvard.edu/abs/1973ARA%26A..11...51P|ADS]]). Two measurements are made: one of the sky and one of an ambient temperature absorber (hot load or chopper) placed into the beam at a point between the receiver and the secondary mirror. As shown below, these measurements provide the information necessary to determine //T//A*, the antenna temperature corrected for atmospheric absorption and ambient temperature telescope losses, including hot spillover and blockage.
Temperature calibration is performed with the ''Cal'' command. This is done automatically at the start of the next scan after changing sources, retuning the receiver, etc. The results of the most recent appropriate ''Cal'' measurement are automatically applied so the recorded data are calibrated on the //T//A* scale.
=== Theory ===
Define the following:
| //α// | = hot spillover efficiency ||
| | = 1 - (fraction of power falling on ground, etc.) ||
| //β// | = cold spillover efficiency ||
| | = 1 - (fraction of power falling on sky, but not forming part of the beam) ||
| //γ// |= source coupling efficiency ||
| //T//RX | = receiver noise temperature | [K] |
| //T//h | = hot load temperature (= ground temperature = air temperature) | [K] |
| //T//c | = cold load temperature | [K] |
| //T//s,sig | = source temperature, signal sideband | [K] |
| //T//s,img | = source temperature, image sideband | [K] |
| //T//A* | = corrected antenna temperature, SSB | [K] |
| //V//h | = hot load signal | [V] |
| //V//c | = cold load signal | [V] |
| //Y// | = //V//h / //V//c ||
| //V//sky | = sky (off position) signal | [V] |
| //V//s | = source and sky (on position) signal | [V] |
| //G// | = //G//sig + //G//img = overall system gain | [V K-1] |
| //G//sig | = system gain, signal sideband | [V K-1] |
| //G//img | = system gain, image sideband | [V K-1] |
| //τ//z | = zenith atmospheric optical depth at the observing frequency ||
| //τ//sig | = zenith atmospheric optical depth at the signal frequency ||
| //τ//img | = zenith atmospheric optical depth at the image frequency ||
| //C//SB | = Sideband atmospheric correction ratio ||
| //A// | = secant (zenith angle) = airmass ||
Several assumptions simplify the theory. The major ones are:
- The hot load, ground, and air are the same temperature, //T//h = //T//ground = //T//air. Often //T//air may be less than //T//h by about 20 K, but this is not too important. The observatory software uses //T//h = 280 K.
- The receiver sideband gains are equal, //G//sig = //G//sig.
- All temperatures are large compared to the observing frequency, i. e., //kT// >> //hυ//. The lowest temperature is likely to be //T//c, which is about 70 K. The Rayleigh-Jeans approximation is valid, therefore, over the submillimeter range.
At the CSO, all the receivers are double sideband (DSB). Continuum (broad band) sources, such as calibrators, spillover, the atmosphere, and planets, contribute to the measured signal in both sidebands. Spectral lines, however, only contribute in one sideband.
Looking at the calibrator (hot load or chopper), the measured signal includes contributions from the receiver and the load,
* //V//h = (//G//sig+//G//img) [//T//RX + //T//h ].
Looking at the sky (off position) and at the source and sky (on position), the measured signal include contributions from the receiver, the atmosphere, the ground (hot) spillover, the sky (cold) spillover, and the source,
| //V//sky | = (//G//sig+//G//img) | × [//T//RX | + (1 - //α//)//T//h ] | + //G//sig //αβ// (1 - //e//-//τ//sig//A// )//T//h | + //G//img //αβ// (1 - //e//-//τ//img//A// )//T//h | + //G//sig //α//(1 - //β//)(1 - //e//-//τ//sig//A// )//T//h | + //G//img //α//(1 - //β//)(1 - //e//-//τ//img//A// )//T//h | , |
| //V//s | = (//G//sig+//G//img) | × [//T//RX | + (1 - //α//)//T//h ] | + //G//sig //αβ// (1 - //e//-//τ//sig//A// )//T//h | + //G//img //αβ// (1 - //e//-//τ//img//A// )//T//h | + //G//sig //α//(1 - //β//)(1 - //e//-//τ//sig//A// )//T//h | + //G//img //α//(1 - //β//)(1 - //e//-//τ//img//A// )//T//h | + //αβγ// [ //G//sig //T//s,sig //e//-//τ//sig//A// + //G//img //T//s,img //e//-//τ//img//A// ]. |
| || receiver | hot spillover | atmosphere || cold spillover || source |
or
| //V//sky | = (//G//sig+//G//img) //T//RX | + //G//sig (1 - //αe//-//τ//sig//A// )//T//h | + //G//img (1 - //αe//-//τ//img//A// )//T//h | , |
| //V//s | = (//G//sig+//G//img) //T//RX | + //G//sig (1 - //αe//-//τ//sig//A// )//T//h | + //G//img (1 - //αe//-//τ//img//A// )//T//h | + //αβγ// [ //G//sig //T//s,img //e//-//τ//sig//A// + //G//img //T//s,img //e//-//τ//img//A// ]. |
Observations are made with some switching scheme, i. e., position switching. The observing procedure measures
* ''OO'' = (//V//s - //V//sky)///V//sky = //αβγ// [ //G//sig //T//s,sig //e//-//τ//sig//A// + //G//img //T//s,img //e//-//τ//img//A// ] / //V//sky ,
the calibration procedure measures
* ''Cal'' = (//V//h - //V//sky)///V//sky = //α// //T//h [ //G//sig //e//-//τ//sig//A// + //G//img //e//-//τ//img//A// ] / //V//sky .
When the sideband gains are equal, //G//sig = //G//sig, the ratio ''OO''/''Cal'' is independent of the hot spillover efficiency, //α//, and of the system gain, //G//sig + //G//img. The data are recorded as SSB antenna temperature, corrected for atmospheric absorption and hot spillover,
* //T//A* = 2//T//h × ( ''OO'' / ''Cal'' ) = 2//βγ// (//T//s,sig //e//-//τ//sig//A// + //T//s,img //e//-//τ//img//A// ) / ( //e//-//τ//sig//A// + //e//-//τ//img//A// ).
=== Sideband correction ===
== Lines ==
Individual spectral lines appear only in one sideband, so //T//s,img = 0 and
* //βγ// //T//s,sig [line] = //C//SB //T//A* ,
where the sideband atmospheric correction ratio
* //C//SB = {[1 + //e//(//τ//sig-//τ//img)//A//]/2}.
The correction ratio may be calculated with an atmospheric [[http://cso.caltech.edu/atm/tuning.html|model]].
Note the correction ratio depends on the observing zenith angle.
If the atmospheric transmission in the two sidebands is equal, //τ//sig = //τ//img and //C//SB = 1, independent of zenith angle, and the recorded data are properly calibrated.
If the atmospheric transmission in the two sidebands is unequal, //τ//sig ≠ //τ//img, then for proper calibration the recorded data should be multiplied by //C//SB ≠ 1.
If the atmospheric transmission in the image sideband is better than in the signal sideband, //τ//sig > //τ//img and //C// > 1.
In this case the atmospheric contribution to the system noise in the image sideband is less than the contribution in the signal sideband.
== Continuum sources ==
For continuum sources, the source temperature is the same in both sidebands, //T//s,sig = //T//s,img = //T//s, so
* //βγ// //T//s [continuum] = //T//A* / 2 .
This is independent of the relative atmospheric transmission in the two sidebands.
=== System Temperature ===
The system temperature, //T//sys, indicates the system noise level corrected for atmospheric absorption and hot spillover. Among other things, it is used for weighting scans when summing them. The system temperature includes contributions in both sidebands from the receiver, from the atmosphere, and from hot spillover. For line (SSB) observations,
* //T//sys = (//e////τ//sig//A// / //α//) {[//T//rx + //α//(1 - //e//-//τ//sig//A//)//T//atm + (1 - //α//)//T//spill ] + [//T//rx + //α//(1 - //e//-//τ//img//A//)//T//atm + (1 - //α//)//T//spill]} .
If the air temperature and spillover temperature are the same, //T//atm = //T//spill = //T//h, then
* //T//sys = (//e////τ//sig//A// / //α//) [2//T//RX + (2 - //α////e//-//τ//sig//A// - //α////e//-//τ//img//A//)//T//h] .
The system temperature can be determined from a calibration measurement,
* ''Cal'' = [//α// (//e//-//τ//sig//A// + //e//-//τ//img//A//) //T//h] / [2//T//RX + (2 - //α////e//-//τ//sig//A// - //α////e//-//τ//img//A//)//T//h] ,
so
* //T//sys = {[1 + //e//(//τ//sig-//τ//img)//A//]/2} × 2//T//h / ''Cal'' = //C//SB 2//T//h / ''Cal'' .
The observatory software assumes the atmospheric transmission in the two sidebands is equal, //τ//sig = //τ//img, i. e., //C//SB = 1, and reports
* //T//sys = 2//T//h / ''Cal'' .
=== Receiver Temperature and Y factor ===
The receiver temperature, //T//RX, and //Y// factor may be determined from measurements of two loads, hot (ambient) and cold (LN2).
* //V//h = (//G//sig+//G//img) [//T//RX + //T//h ] ,
* //V//c = (//G//sig+//G//img) [//T//RX + //T//c ] ,
* //Y// = //V//h / //V//c = [//T//RX + //T//h ] / [//T//RX + //T//c ] , and
* //T//RX = (//T//h - //Y// //T//c) / (//Y// - 1) .
They may also be determined by substituting a cold load for the sky when making a (manual) calibration.
* ''Cal'' = (//V//h - //V//c)///V//c ,
* //T//sys = 2//T//h / ''Cal'' ,
* //Y// = //V//h / //V//c = 1 + ''Cal'' = 1 + 2 //T//h / //T//sys , and
* //T//RX = (//T//h - //Y// //T//c) / (//Y// - 1) .
For a piece of Eccosorb drenched in LN2, it is usually assumed that //T//c = 80 K.
===== Telescope efficiency =====
To correct for cold spillover and for the coupling of the beam to the source, it is necessary to measure the telescope efficiency, //βγ//. Then the source temperature may be determined,
* //T//s = //T//A*///βγ// .
The actual value of //βγ//, and what you may wish to call it, depends on the size of the source. There are three typical cases: sources large, similar, or small compared to the beam.
- For large sources, those a few times the beam size, determine //βγ// from measurements of the full moon, which fills the main beam and the near side lobes, \\ //βγ// (Moon) = //T//A*(Moon)/2//T//Moon , \\ where //T//A*(Moon) is the observed (SSB) lunar brightness and //T//Moon is the intrinsic lunar temperature. Probably you will still want to call data calibrated this way //T//A*, but note that they are calibrated with respect to the Moon.\\ \\
- For sources comparable to the beam size, determine //βγ// from measurements of a planet. Usually this called the main beam temperature scale, //T//mb. The best planets for this are Jupiter and Mars, with Venus possibly OK. Note the planet temperatures depend on observing frequency (cf. Butler 2012 ALMA [[https://science.nrao.edu/facilities/alma/aboutALMA/Technology/ALMA_Memo_Series/alma594/abs594|memo 594]]). Mars is the most reliable but it does have absorption lines in its atmosphere. For critical observations, measure the telescope efficiency at different zenith angles. For a planet of uniform temperature, calculate \\ //βγ// (Main beam) = (//T//A*(Planet)/2//T//Planet) × [1 - exp(-(//D/////Θ//)2 ln(2))]-1 , \\ where //T//A*(Planet) is the observed (SSB) planetary brightness and //T//Planet is the intrinsic planetary temperature, //D// is the planet diameter, and //Θ// the half power full beamwidth. This calculation assumes that the main beam is gaussian. The {{beam_eff.astro}} and {{beam_eff.class}} scripts perform this calculation in ''Astro'' and ''Class'', respectively.\\ \\
- Small sources can be calibrated using the main beam efficiency or, more traditionally, using the aperture efficiency, which is very similar except that it is assumed that the beam is the diffraction function of a uniformly illuminated, perfect telescope, rather than a gaussian. In that case, slightly higher temperatures will be obtained and the result is called brightness temperature, //T//b. This is not a good unit because, in its original definition, it requires a correction for the Planck factor also. This makes sense for very hot radio sources, but not for low temperature molecular sources. To avoid this difficulty, we should either invent yet another temperature unit or use Janskys. I suggest you use Janskys.
For line observations : \\
| //T//A* | = 2//T//h × ( ''OO'' / ''Cal'' ) , | recorded data |
| //T//A*(Moon-corrected) | = //T//A* / //βγ// (Moon), | large source |
| //T//MB | = //T//A* / //βγ// (Main Beam) | smallish source |
The efficiencies vary with the receiver and, strictly speaking, with frequency for each receiver. It is best to measure the efficiency yourself, but approximate values are:
| Receiver | 230 GHz | 345 GHz | 492 GHz | 650 GHz | 850 GHz |
| //βγ// (Main beam) | 76% | 65% | 53% | 32% ?? | |
===== Skydip Measurements =====
=== Theory ===
We will use many of the same quantities and some of the results of the previous section. In addition, let us define the following quantities:
| //T//spillover | = spillover temperature ≡ (1 - //η//hot) //T//h |
| //τ//z | = atmospheric opacity at zenith |
| //A// | = airmass at a given zenith angle = 1/cos(''ZA'') (assuming a plane-parallel atmosphere) |
| //Y// | = //V//h///V//c = Receiver //Y// factor |
[//η//hot) ≡ //α//.]
As noted, we assume a plane-parallel model for the atmosphere, so for a given zenith angle we have\\
//τ// = //τ//z //A//.
We will be comparing measurements at different zenith angles rather than on- and off-source measurements. The entire passband of the receiver is treated as a single channel (the total power output is used), and we assume that the contribution of any source to the total power is negligible in comparison with the contributions of the receiver, spillover, and the atmosphere.
From the previous section, we have \\
| //V//h | = //G// [//T//RX + //T//h ] |
| //V//c | = //G// [//T//RX + //T//c ] |
and\\
//V//sky = //G// [//T//RX + (1 - //η//hot//e//-//τ// )//T//h ]. \\
If we now form the quantity \\
| //S// | = ln [(//V//h - //V//c)/(//V//h - //V//sky)] |
| | = ln [(//T//h - //T//c)/ //η//hot //e//-//τ// //T//h ] |
| | = //τ// + ln [(//T//h - //T//c)/ //η//hot //T//h ] |
| | = //τ//z//A// + ln [(//T//h - //T//c)/ //η//hot //T//h ] , |
we see that //A// is a linear function of airmass with slope //τ//Z. With //T//h and //T//c known, //η//hot can be found from the intercept as \\
//η//hot = (1 - //T//h///T//c) e-(intercept). \\
//T//spillover can then be calculated from its definition. Note that for any measurement of //V//sky we may define a //T//equiv by requiring \\
//V//sky = //G// [ //T//RX + //T//equiv ].\\
//T//spillover can then also be thought of as the //T//equiv that would be measured for an airmass of zero.
If //V//c (or //T//c) is unknown, //τ//z can be found by defining \\
| //S'// | = ln [//V//h/(//V//h - //V//sky)] |
| | = ln [(//T//h + //T//RX)/ //η//hot //e//-//τ// //T//h ] |
| | = //τ//z//A// + ln [(//T//h + //T//RX) / //η//hot //T//h ] , |
where now //S'// is linear in airmass with a slope of //τ//z. If //T//RX is known, this intercept can be used to find //η//hot. However, it is more likely that //T//RX will be determined from hot and cold load measurements of //V//h and //V//c at known temperatures as follows: \\
//Y// = //V//h///V//c = (//T//RX + //T//h)/(//T//RX + //T//c) \\
and \\
//T//RX = (//T//h - Y //T//c)/(//Y// - 1) .\\
(This is of course unrelated to doing a skydip except insofar as measurements are typically made of //V//h and //V//c.)
Whether //S// or //S'// is used, it is necessary to make more than one measurement in order to find the slope and the intercept. The usual way of doing this is outlined below.
=== Practice ===
//N. B. This section is obsolete. The ''UIP'' ''SKYDIP'' macro will position the telescope to successive airmasses.//
A skydip is usually performed at the CSO using the ''1M'', ''1.5M'', etc., macros under the ''UIP''. Issuing one of these macros causes the telescope to go to the corresponding zenith angle for 1, 1.5, 2, 3, 4, or 5 airmasses. The total power must be recorded at each of these positions (some day it may be automated, but not yet) and for the hot and cold loads.
A utility program called TAUPLOT exists at the CSO which will do all of the necessary calculations and produce pretty plots. This may be run by simply typing ''tauplot'' at either the ''UIP>'' prompt or the ''$'' prompt. The program will then prompt you for an airmass and voltage, and will repeat the prompt after the values are entered. A simple carriage return will cause the program to go on, and you are prompted for //V//h and then //V//c. Default values for //T//h and //V//c will be used unless other values are specified at the next prompt.
IMPORTANT NOTE - The program is set up under the assumption that //G// is positive (it checks that //V//sky is always less than //V//h). This means that although the total power voltage may be displayed as negative, the absolute value of this should be used. Furthermore, it is often more convenient to use integral values though this is not necessary. Thus, if the meter reads -0.142, a value of 142 should be entered.
Once the data entry is complete, the program calculates //τ//z, the intercept, //η//hot, //T//spillover, //Y//, and //T//RX, as well as equivalent temperatures //T//equiv for all of the measured //V//sky. These are displayed on the screen, and you are asked if you wish to plot the results. The default is to produce first a screen plot and then hardcopy, but you may specify otherwise if you desire. All of the calculated values are also displayed on the plot.
If for some reason, //V//c is not measured, the program may be run using the //A//' formalism outlined above by not entering a value at the //V//c prompt. In this case, //τ//z is the only value calculated and is displayed with a message to that effect.
===== AOS Frequency Calibration =====
//N. B. This section is obsolete.//
Frequency calibration of the AOS backends at the CSO is done using a multiplied crystal oscillator to generate a comb with peaks separated by 100 MHz for the 500 MHz AOS or 10 MHz for the 50 MHz AOS. These allow one to determine the relationship between the IF frequency //υ//IF and channel number of the AOS. (This relationship is assumed to be linear and experience shows this to be so.)
When the ''LO'' command is issued, the LO frequency //υ//LO is set so that the desired line frequency will correspond to the center of the passband, or approximately channel 512. This will take into account velocity and frequency offsets as well as the correct //υ//IF. A default relationship between channel number and frequency is then written into the observation headers. However, because the AOS exhibits some variation in response to temperature changes, the calibration comb should be used to determine the precise relationship as follows:
* Start by using the ''FCAL'' macro to record a scan showing the response of the AOS to the comb. Then use CLASS to look at this scan using units of channel numbers. \\ Determine the channel numbers corresponding to the peaks of the comb. As a peak may cover more than one channel, it may be necessary to estimate where the center of the peak occurs. Note that as displayed in CLASS, channel //N// extends from //N//-0.5 to //N//+0.5 and that for the purposes of these calculations, non-integral channel numbers should be used, e. g., a peak may occur at channel number 308.7.
* Next use the number of channels separating the peaks to find the width of each channel in MHz. These values should be roughly the total backend width (500 MHz or 50 MHz) divided by the total number of channels (1024), or about 0.49 MHz or 0.049 MHz. The channel width in velocity is determined from the channel width in frequency and is calculated for the desired line frequency.
* Finally determine the channel number which actually corresponds to the desired line frequency. For the 500 MHz AOS, this is the channel location of the central peak. For the 50 MHz AOS, this is the channel location halfway between the two central peaks. This should not vary by more than a few channels from the default channel location of 511.5 for both of the backends.
If after determining the correct values for the reference channel number (corresponding to the desired line frequency) and the width of the channels you desire to change from the default values to the correct values, you can do this within CLASS for data taken in the IRAM format. At the present time, there is no provision made for changing the default values that will be written as you observe, so corrections must be made after the fact. As you will want to write out the scans after they have been corrected, this means you must first open a new data file for writing. See the CLASS manual for a description of the ''FILE OPEN'' command if you are not sure how to do this.
Variables exist corresponding to the reference channel location and the channel spacing in both frequency and velocity which are used by CLASS in plotting, etc. The names of these variables are ''REFERENCE'', ''FREQ_STEP'', and ''VELO_STEP''. The command ''EXAMINE variable_name'' can be used to display the current values, while the command ''LET variable_name = value'' can be used to change the values. Unfortunately, these new values do not become effective until the scan is written out to a data file and the new scan then read back in. Alternatively, the command ''MODIFY'' can be used (with conveniently different variable names, of course). To change the reference channel, type ''MODIFY RECENTER value''; to change the channel frequency spacing, type ''MODIFY WIDTH value''. Using ''MODIFY'' has some advantages: the values become effective immediately, and changing the frequency spacing automatically causes the velocity spacing to be changed correctly. However, with either ''LET'' or ''MODIFY'' //the scan must be written to a data file or the changes will be lost.// This can be done using either ''WRITE'' or ''UPDATE''. Again, see the CLASS manual for more details on these commands.
Some conditions for which having the correct values is essential would include: adding together spectra where a line has been observed in different parts of the passband; observing weak or narrow lines which might be lost if any smearing were to occur because of AOS drifts; observing lines near the edge of the passband, or very wide lines. Note that because all values are calculated from the reference channel at around 512, errors in the channel width value become progressively more important as you move toward the edge of the passband.
At the present time, the frequency calibration comb is switched in only upon the user's request (or if the antenna is sent on a very long slew). The ''FCAL'' macro in the ''UIP'' program does this; to do this in a command file you must give the actual series of commands which correspond to the ''FCAL'' macro. Type ''DEFINE/LIST FCAL'' to see this series of commands.