Radiometer Equation Applied to Telescopes
From AstroBaki
[edit] Prerequisites
[edit] Short Topical Videos
[edit] Reference Material
- * Synthesis Imaging in Radio Astronomy II, ed. Taylor, Carilli, Perley, Ch. 9 (Wrobel & Walker)
- Much of the preceding is based on the Radio 101 single-dish lecture, Single Dish Basics.
- Handwritten notes used to make the video by Chat Hull
- The NRAO Course on Radiometers
Contents |
The Radiometer Equation
The radiometer equation is really all about signal to noise. Here’s the answer:
where:
- T_{src} is the signal of the source you’re observing
- T_{sys} is your system temperature
- T_{rms} is the noise in your system, or the RMS fluctuations in your system temperature
- Δν is the bandwidth of your correlator (in Hz)
- τ is integration time (seconds)
What’s the deal with all the temperatures?
All of the temperatures in the above equation have their origin in the Rayleigh-Jeans limit of the blackbody equation:
Radio astronomers like to refer to the above temperature T as the "brightness temperature," T_{B}, of a source. T_{B} is the temperature of the blackbody needed to produce the observed specific intensity at that frequency:
Now, back to specific intensity. I_{ν} is defined as:
where:
- Δt is the exposure time
- Δν is a small interval in frequency
- ΔE is the energy emitted in bandwidth Δν over time Δt
- ΔA is the area of the telescope
- ΔΩ is the observed solid angle on the source, or your beam size (steradians)
So,
It is common to report I_{ν} in Jy/beam, where the beam has units of Sr^{ − 1}, and Janskys (Jy) are the units of flux density S_{ν}, defined below.
To arrive at other quantities, we must integrate over I_{ν}:
- Flux Density: . The fundamental unit of flux density is the Jansky.
- Power received: .
Source Temperature, T_{src}
The source temperature is defined as the brightness temperature associated with the power received by the telescope from the source you’re observing:
The last step comes from the Antenna Theorem, which states that dAdΩ = λ^{2}.
K / Jy, or "forward gain"
Now that we’ve related the power P_{ν} received by the antenna at a given frequency, to the source temperature, T_{src}, and the observed intensity, I_{ν}, we can arrive at a conversion between flux density and brightness temperature for an unresolved source:
Therefore, the conversion between K and Jy, also known as the "forward gain" of an antenna, is just
Thus, brightness temperature is just another measure of the brightness of a source. The forward gain is a physical property of a telescope, which dictates the telescope receivers’ response to a given increase in Janskys.
The System Temperature
While the source temperature T_{src} describes the energy received from the source you are interested in, the system temperature, T_{sys} describes the actual power received due to both the sky (T_{sky}) and the receivers (T_{Rx}):
The main component is the receiver temperature, T_{Rx}, which comes from the thermal (or Johnson) noise from the receiver electronics. Depending on the observation waveband, receivers are sometimes cooled to reduce T_{Rx}.
T_{sky} is the power generated by everything you don’t want to be looking at on the sky: background sources, water vapor in the atmosphere (especially in mm-wave astronomy), galactic backgrounds (especially in long-wavelength astronomy), etc.
Effective Area A_{e}
Above we have written the area of an antenna as A_{e}, or the "effective area."
where (η_{a} is the aperture efficiency and A_{p} is the projected area of the telescope). η_{a} is a number less than one, and signifies that not all radiation incident upon the telescope actually makes it to the receiver because of the dish’s finite reflectivity. This loss of signal requires that the signal later be amplified. There are two scenarios that can be encountered:
- T_{sys} is dominated by T_{Rx}: in this case aperture efficiency matters, because in amplifying the signal, noise from the receiver is also amplified.
- T_{sys} is dominated by T_{sky}: in this case, which might happen on a cloudy day at a millimeter telescope where the noise from the atmosphere is dominant, aperture efficiency doesn’t matter. This is because the dominant sky noise is first cut down by a low aperture efficiency, and then multiplied back up to its original level by the amplifiers.
Detecting the signal...
Typically T_{sys} > > T_{src}, so how do we detect the source we are interested in? Beat down the noise!
In order to detect the source, we need T_{src} > T_{rms}, where the observation-dependent quantity T_{rms} is the noise in our measurement of the observation-independent quantity T_{sys}:
where N is the number of independent data points.
For a telescope, the number of independent samples is , where Δν is the bandwidth (Hz) and τ is the integration time (seconds). With a bandwidth Δν, the signal is statistically independent over a time interval 1 / (2Δν) (the factor of 2 comes from the Nyquist sampling rate), so that the number of independent samples is just τ divided by 1 / (2Δν), so N = 2τΔν. Then, as explained in the NRAO course, the fluctuations in a given measurement of T_{sys} are , and the factor of two disappears, giving us the final result:
It is important to note that the "noise" (or RMS variations) in any given measurement is proportional to the uncertainty in each measurement, divided by . The fact that is not because T_{sys} represents an uncertainty in each measurement! It is because, in the Rayleigh-Jeans limit of the blackbody equation, the uncertainty in a measurement in a given mode is proportional to the number of photons N_{γ} in that mode. And since , it follows that . Even though T_{sys} is sometimes called "noise temperature," we emphasize here that T_{sys} is a real signal; it does not represent real "noise," or fluctuations in measurements.
The Radiometer Equation
Now we can write down an expression for the signal to noise ratio (the radiometer equation):
Typical values might be Δν = 10 MHz, τ = 1 sec so that . Typical values for T_{sys} are 40 - 200 K.
SEFD
The SEFD is the ‘system equivalent flux density’, which is the flux density equivalent of T_{sys}:
The SEFD is a useful way to compare the sensitivity of two different systems since it folds in both T_{sys} and A_{e}. This also greatly simplifies the sensitivity calculation: if you know the flux in Jy of the source you want to detect, and you know the SEFD, then you can easily calculate the integration time you need to make a given S/N detection (for an unresolved source):
After substituting in the temperature-based Radiometer Equation for S, we arrive at an intuitive expression for the RMS variations in flux density S_{ν,rms}:
In the temperature-based Radiometer Equation, signal-to-noise increases with increased bandwidth and integration time. In the flux-density-based Radiometer equation, RMS flux density varations decrease with increased bandwidth and integration time.