6.1 Specifications

Table of Contents

Return to the Table of Contents for this section at any time by selecting Return to Table of Contents. Users may also navigate through the entire Observer's Handbook by using the complete Table of Contents menu to the right.

6.1   Specifications

6.1.1   Instrument Overview

The German REceiver for Astronomy at Terahertz Frequencies (GREAT) is a dual channel heterodyne instrument that will provide high resolution spectra (up to R = 108) in several frequency windows in the 0.490–4.74477749 THz range.

The front-end unit consists of two independent dewars that operate simultaneously. Each dewar contains one of the following channels:

upGREAT Low Frequency Array (LFA)

The LFA is a 14-pixel array (two polarizations and times; 7 pixels) arranged in a hexagonal geometry with a central beam. The current tuning range of the LFA is approximately 1.810–2.070 THz.
 
upGREAT High Frequency Array (HFA)
The HFA is a 7-pixel array arranged in a hexagonal geometry with a central beam. The tuning range of the HFA is limited to the [OI] frequency of 4.74477749 THz (+/- 100 km/s LSR velocity).
 
4GREAT
4GREAT consists of four co-aligned (to within a few arcsec) pixels at four different frequencies. The tuning ranges of these pixels are 490–635 GHz (Herschel/HIFI band 1), 890–1100 GHz (Herschel/HIFI band 4), 1200–1500 GHz (GREAT L1 channel), and 2490–2590 GHz (GREAT M channel).
 

The GREAT instrument uses Fast Fourier Transform (FFT) spectrometers as backends. Each XFFTS spectrometer has a bandwidth of 2.5 GHz and 64,000 channels, providing a resolution of 44 kHz. The beam size is close to the diffraction limit—about 14 arcsec at 160 μm.

A detailed description of the GREAT instrument and its performance during Basic Science can be found in the GREAT special issue (Heyminck et al. 2012, A&A, 542, L1). The LFA is described in Risacher et al. 2016 (A&A, 595, A34). Because the instrument is regularly being upgraded and its performance improved, always check the GREAT website for the latest information.

Return to Table of Contents

6.1.1.1   Design

Heterodyne receivers work by mixing the signal from a source at a given frequency νs with that from a local oscillator (LO) at a specified (and precisely controlled) frequency νLO and amplifying the result. The mixing results in two frequency bands, called the signal and the image bands, located symmetrically on either side of νLO and separated from νLO by the intermediate frequency νIF = |ν– νLO|. GREAT operates in double sideband mode, i.e. both the image and signal bands are equally sensitive to incoming radiation. By definition the spectral line of interest is always placed in the signal band, which can be chosen to be either above (Upper Sideband, USB) or below the LO frequency (Lower Sideband, LSB); see Figure 6-1. For sources rich in spectral lines, care has to be taken so that a spectral line in the image band does not overlap or blend with the line in the signal band.

Figure 6-1.
Diagram showing how heterodyne receiver works

Figure 6-1. Schematics of one GREAT receiver channel. The second mixer-stage is needed to match the ope­­ra­tional frequencies of the first mixer-element to that of the microwave spectrometer.

Return to Table of Contents

6.1.1.2   Configurations

The channels that are currently operational are listed in Table 6-1. Not every frequency in each tuning range has been checked, so there may be gaps where the LOs do not provide enough power to pump the mixer.

 
Table 6-1: GREAT Configurations
Front-End Frequencies (GHz) Lines of Interest DSB5 Receiver Temperatures (K)
upGREAT LFA 1810–2070 [OI], [CΙΙ], CO, OH2π1/2 1000
upGREAT HFA 4744.77749 [OI] 1100
4GREAT 490–6351 NH3, [CI], CO, CH 120
890–11002 CO, CS 350
1200–15003 [NII], CO, OD, HCN, SH, H2D+ 800
2490–2590 OH2π3/2 18OH2π3/2  1500
1Herschel/HIFI Band 1
2Herschel/HIFI Band 4
3GREAT L1 channel
4GREAT M channel
5DSB = Dual Sideband
 

In Cycle 8, GREAT offers two configurations:

upGREAT LFA with upGREAT HFA
This configuration employs the Low Frequency Array (LFA) in one dewar and the High Frequency Array (HFA) in the second dewar. The LFA array consists of 2 x 7 pixels (H and V polarizations) and the HFA consists of seven pixels. At the frequency of the CII line, 1900.536 GHz, the frequency coverage of both LFA polarizations overlap and all 14 pixels are in operation. At higher and lower frequencies, only one LFA polarization will be available and only seven LFA pixels can be operated at a given time.
 
4GREAT with upGREAT HFA
This configuration implements 4GREAT in one dewar and the HFA in the second dewar. The 4GREAT pixels are aligned on the sky to within several arcseconds.
 

Each of the two installed mixers is provided with XFFTS Fast Fourier Transform spectrometers. The XFFTS has a 2.5 GHz bandwidth and 64,000 channels, providing a resolution of 44 kHz.

The usable instantaneous bandwidth is generally less than that covered by the XFFTS, generally about 1-1.5 GHz. The tuning range of the HFA is limited because the LO (a Quantum Cascade Laser) can only be tuned in discrete steps.

Return to Table of Contents

6.1.2   Performance

6.1.2.1   Imaging Sensitivities

The GREAT sensitivities and integration times are calculated with the SOFIA Instrument Time Estimator (SITE). Presented here are the background of these calculations and some worked out examples. Because of the way a heterodyne receiver is calibrated (by measuring the receiver temperature, Trx, with a hot and a cold load), the logical intensity unit for a heterodyne observation is temperature, expressed in Kelvin (K). Either the antenna temperature, TA* (the asterisk refers to values after correction for sky transmission, telescope losses and rearward spillover, see e.g. Kutner and Ulich, ApJ, 250, 341 (1981)) or the main beam brightness temperature, Tmb, are used. Similarly, the noise is expressed in temperature units as well, ΔTA * or ΔTmb, and the sensitivity (or signal-to-noise ratio) of the observations is given by the ratio of the source temperature and the rms temperature of the spectrum. In order to calculate these quantities, we first must estimate the single sideband (SSB) system temperature, Tsys, which also includes losses from the atmosphere and the telescope. Tsys is given by:

(Eq.6-1)
Tsys = 2 x [Trx + ηtelJ x Tsky + ηtel x Ts + Ttel ] /( ηtel x ηsky)

where
Trx is the dual sideband (DSB) receiver temperature
Tsky is the radiation temperature of the sky
Ts is the continuum temperature of the source; This term is usually completely negligible (unless observing an extremely bright object such as Jupiter, Saturn, or Sgr B2) and is therefore set to 0 in the GREAT time estimator.
Ttel is the Rayleigh-Jeans telescope temperature
ηsky is the fraction of radiation transmitted through the atmosphere
ηtel is the efficiency of the telescope, which includes ohmic losses and spillover
 

The factor 2 in Equation 6-1 assumes that the noise temperature is the same in both signal and image band, which is true for the HEB mixers used by GREAT. The transmission of the atmosphere, ηsky, at the altitude, observing frequency, and airmass that we plan to observe at can be estimated using the atmospheric transmission code ATRAN. The GREAT time estimator calls ATRAN directly, so those estimating the integration time needed by using the time estimator have no need to run ATRAN separately. Tsky depends on ηsky and the ambient temperature of the sky (Tamb) where the signal is absorbed and can be derived from Equation 6-2.

(Eq. 6-2)
Tsky = J(Tamb) x (1 - ηsky)

where J(Tamb) is the mean Rayleigh-Jeans (R-J) temperature of the atmosphere, which we assume to have a physical temperature of 220 K at 41,000 ft, resulting in J(Tamb) = 177.5 K at 1.9 THz. Likewise the telescope temperature, Ttel, is related to ηtel by Equation 6-3:

(Eq. 6-3)
Ttel = Jtel x (1 - ηtel)

where Jtel is the radiation temperature of the telescope, with a physical temperature ~ 230 K (Jtel = 187.4 K at 1.9 THz). If we assume an ηtel of 0.92, then Ttel = 14.8 K.

As an example, let us calculate the system temperature at the [CII] fine structure line at 157.74 μm (1.9005369 THz). In this example we calculate what we would have at the beginning of a flight, when we are still at low altitude. We therefore assume that we fly at an altitude of 39,000 ft and observe at an elevation of 30 degree. For a standard atmospheric model this corresponds to a transmission of ~ 76%, which gives Tsky = 42.6 K. For a receiver temperature Trc = 1100 K, Equation 6-1 therefore predicts a single sideband system temperature Tsys = 3301 K when observing the sky.

Now we are ready to calculate the sensitivity. The rms antenna temperature, (corrected for the atmospheric absorption, and telescope losses), ΔTA*, for both position switching
and beam switching is given by Equation 6-4

(Eq. 6-4)
ΔTA* = (2 x Tsys x κ ) x (t x Δν)-0.5

where k is the backend degradation factor, t is the total integration time of the number of on and off pairs that we plan to take, and Δν is the frequency resolution of our spectra. Strictly speaking, Δν is the noise bandwidth, which can be slightly different than the frequency resolution, depending on the design of the spectrometer. For our example we expect the Full Width Half Maximum (FWHM) of the line to be a few km/s, and we will therefore calculate the rms for a velocity resolution of 1 km/s, corresponding to a frequency resolution of 6.3 MHz. Since the GREAT backends have much higher resolution, this is not a problem. We can easily bin the spectrum to our desired velocity resolution. For an observation with three pairs of 40 seconds in each beam, or t = 4 min, and assuming the backend degradation factor k=1, we then find ΔTA* = 0.17 K, which is the one sigma rms antenna temperature.

To convert antenna temperature to brightness temperature TR*, we have to make one more correction as shown in Equation 6-5:

(Eq. 6-5)
TR* = TA* / ηfss

where ηfss is the forward scattering efficiency, usually measured for a very extended source (like the Moon). For GREAT ηfss = 0.97. Therefore our brightness rms temperature, ΔTR* = 0.18 K5. Note: The GREAT time estimator assumes the line temperature in TR*, and not in main beam brightness temperature, Tmb = TR*MB , which is for a source that just fills the main beam.

If we want to express our results in flux density, Sn, rather than brightness temperature, we can convert antenna temperature, TA, to flux density, Sν, using the standard relation given in Equation 6-6:

(Eq. 6-6)
Sν = 2 x κ x ηfss x TA* / Aeff  

where κ is the Boltzmann constant, and Aeff is the effective area of the telescope. Aeff is related to the geometrical surface area of the telescope, Ag, by the aperture efficiency, ηa,
i.e. Aeff = ηa x Ag. For the measured main beam efficiency in early April 2013 (0.67) and a Half Power Beam Width (HPBW) of ~ 14.1 arcsec (+/- 0.3 arcsec) an aperture efficiency6 of 55 +/- 2 % is derived. Equation 6-7 yields the following simple form for the 2.5 m SOFIA telescope:

(Eq. 6-7)
S (Jy) = 971 x TA* (K) or within errors ~ 1000 x TA* (K)

Normally we use Jy only for spatially unresolved sources, but we can also use Equation 6-7 to convert line intensities into W/m2, which maybe a more familiar unit for the far infrared community. If we assume that the [CII] line we are observing is a Gaussian with a Full With Half Maximum (FWHM) of = 5 km/s, i.e. 31.8 MHz, the integrated line intensity is given by 1.065 x Tpeak x Δν, where Δν = 31.8 MHz. If we take Tpeak equal to our rms antenna temperature, we find using Equation 6-7 that our four minute integration therefore corresponds to a one sigma brightness limit of ~ 6.1 10-17 W/m2 for a 5 km/s wide line observed with 1 km/s resolution. If we only aim for a detection, we can probably degrade the resolution to 2 km/s. In this case we gain a square root of 2, and therefore our one sigma detection limit is 4.3⋅10-17 W/m2.

This is the reverse of what is typical when writing a proposal, when the proposer has an estimate how wide and how bright the line is expected to be and knows what signal to noise is needed for the analysis. Assume we want to observe the [CII] 158 μm line in T Tauri, a young low-mass star. Podio et al. (2012, A&A, 545, A44) find a line intensity of 7.5 ⋅10-16 W m-2 for the Herschel PACS observations, which are unresolved in velocity (the PACS velocity resolution is ~ 240 km/s for [CII]). Here we want to velocity resolve the line to see if it is outflow dominated or whether it is emitted from a circumstellar disk or both. We therefore need a velocity resolution of 1 km/s or better. If we assume that the line is outflow dominated with a FWHM of say 20 km/s (127.2 MHz; little or no contribution from the circumstellar disk) we get a peak antenna temperature (using Equation 6-7) of 0.52 K or a radiation temperature TR* = 0.55 K. In this case we want a SNR of at least 10 and a velocity resolution of 1 km/s or better. Let’s check whether it is feasible. If we plug in the values we have in the GREAT time estimator (assume 40 degrees elevation, standard atmosphere, and we fly at 41,000 ft) or we can estimate it from the equations given above.

With these assumptions ATRAN gives us an atmospheric transmission of 0.86 integrated over the receiver band-pass. The receiver temperature is 1100 K (DSB). Using Equation 6-2, we find that the sky only adds 24.9 K to the system temperature and from Equation 6-1 we therefore get Tsys = 2881 K. Since we want to reach a signal to noise of 10, the rms antenna temperature ΔTA* = 0.052 K. We can now solve for the integration time using Equation 6-4, where we set Δν = 6.338 MHz (1 km/s resolution). In this case t = 1937 sec or 32.3 min. The PACS observations show the emission to be compact, so we can do the observations in Dual Beam Switching Mode (DBS; see Section 6.2.1.1), with a chop throw of 60 arcsecond. Both DBS and Total Power (TP) modes are currently estimated to have an overhead of 100% and a setup time for tuning and calibration of two minutes (which get added when entering the observations in USPOT, the SOFIA proposal tool). Our observation would therefore take 60 minutes, which is completely feasible. SITE gives t = 1930 sec. The difference is negligible.

Sensitivity calculations for an On-the-Fly (OTF) map, when data are taken while the telescope is scanning), are done a bit differently. For example, for an OTF map of the [CII] line (Half Power Beam Width, HPBW ~ 14 arcsec), we need to sample the beam about every 7 arcsec. If we read out the average once per second, for example, this means that we scan with a rate of 7 arcsec/second. To do a 3 arcmin scan will therefore take 26 sec, resulting in 26 map points, let us make it 27, to get an odd number of points. We therefore need to spend (1 second) * (√27) = 5.2 seconds on the reference position in total power mode. We ignore the time it takes the telescope to slew to the next row and any time needed for calibration. For a 3 x 3 arcmin2 map, i.e. 27 x 27 positions with a cell size of 7 x 7 arcsec2. The integration time for each row is therefore 5.2 +27 seconds or 32.2 seconds/row. The total integration time for the map is therefore 14.5 minutes. We definitely want to do one repeat, so the total integration time is therefore 29 minutes. For all observing modes, we assume a 100% overhead and 2-minute setup for tuning and calibration. The total duration of the two maps is therefore 61.9 minutes. Thus, a 3 x 3 arcmin2 map in the [C II] line is entirely feasible.

Our 3 by 3 arcmin2 map with one repeat has an integration time of 2 seconds per map point. For typical observing conditions (41,000 feet, 30 degrees elevation) the GREAT time estimator (settings: TP OTF map, Non=27, Classic OTF; see Section 6.2.3) gives us an rms temperature/map point of 0.9 K for a velocity resolution of 1 km/s.

Return to Table of Contents