Performance of the VLA
VLA capabilities January 2013 - April 2014
1. Introduction
This section contains details of the EVLA's resolution, expected sensitivity, tuning range, dynamic range, pointing accuracy, and modes of operation. Detailed discussions of most of the observing limitations are found elsewhere. In particular, see References 1 and 2, listed in Documentation].
2. Resolution
VLA Resolution
The VLA's resolution is generally diffraction-limited, and thus is set by the array configuration and frequency of observation. It is important to be aware that a synthesis array is "blind" to structures on angular scales both smaller and larger than the range of fringe spacings given by the antenna distribution. For the former limitation, the VLA acts like any single antenna - structures smaller than the diffraction limit (θ ∼ λ/D) are broadened to the resolution of the antenna. The latter limitation is unique to interferometers; it means that structures on angular scales significantly larger than the fringe spacing formed by the shortest baseline are not measured. No subsequent processing can fully recover this missing information, which can only be obtained by observing in a smaller array configuration, by using the mosaicing method, or by utilizing data from an instrument (such as a large single antenna or an array comprising smaller antennas) which provides this information.
Table 5 summarizes the relevant information. This table shows the maximum and minimum antenna separations, the approximate synthesized beam size (full width at half-power) for the central frequency for each band, and the scale at which severe attenuation of large scale structure occurs.
Configuration | A | B | C | D |
---|---|---|---|---|
Bmax (km1) | 36.4 | 11.1 | 3.4 | 1.03 |
Bmin (km1) | 0.68 | 0.21 | 0.0355 | 0.035 |
Band |
Synthesized Beamwidth θHPBW(arcsec)1,2,3 | |||
74 MHz (4 band) | 24 | 80 | 260 | 850 |
350 MHz (P) | 5.6 | 18.5 | 60 | 200 |
1.5 GHz (L) | 1.3 | 4.3 | 14 | 46 |
3.0 GHz (S) | 0.65 | 2.1 | 7.0 | 23 |
6.0 GHz (C) | 0.33 | 1.0 | 3.5 | 12 |
10 GHz (X) | 0.20 | 0.60 | 2.1 | 7.2 |
15 GHz (Ku) | 0.13 | 0.42 | 1.4 | 4.6 |
22 GHz (K) | 0.089 | 0.28 | 0.95 | 3.1 |
33 GHz (Ka) | 0.059 | 0.19 | 0.63 | 2.1 |
45 GHz (Q) | 0.043 | 0.14 | 0.47 | 1.5 |
Largest Angular Scale θLAS(arcsec)1,4 | ||||
74 MHz (4 band) | 800 | 2200 | 20000 | 20000 |
350 MHz (P) | 155 | 515 | 4150 | 4150 |
1.5 GHz (L) | 36 | 120 | 970 | 970 |
3.0 GHz (S) | 18 | 58 | 490 | 490 |
6.0 GHz (C) | 8.9 | 29 | 240 | 240 |
10 GHz (X) | 5.3 | 17 | 145 | 145 |
15 GHz (Ku) | 3.6 | 12 | 97 | 97 |
22 GHz (K) | 2.4 | 7.9 | 66 | 66 |
33 GHz (Ka) | 1.6 | 5.3 | 44 | 44 |
45 GHz (Q) | 1.2 | 3.9 | 32 | 32 |
- These estimates of the synthesized beamwidth are for a uniformly weighted, untapered map produced from a full 12 hour synthesis observation of a source which passes near the zenith.
- Footnotes:
- 1. Bmax is the maximum antenna separation, Bmin is the minimum antenna separation, θHPBW is the synthesized beam width (FWHM), and θLAS is the largest scale structure "visible" to the array.
- 2. The listed resolutions are appropriate for sources with declinations between −15 and 75 degrees. For sources outside this range, the extended north arm hybrid configurations (DnC, CnB, BnA) should be used, and will provide resolutions similar to the smaller configuration of the hybrid, except for declinations south of −30. No double-extended north arm hybrid configuration (e.g., DnB, or CnA) is provided.
- 3. The approximate resolution for a naturally weighted map is about 1.5 times the numbers listed for θHPBW. The values for snapshots are about 1.3 times the listed values.
- 4. The largest angular scale structure is that which can be imaged reasonably well in full synthesis observations. For single snapshot observations the quoted numbers should be divided by two.
- 5. For the C configuration an antenna from the middle of the north arm is moved to the central pad "N1". This results in improved imaging for extended objects, but will degrade snapshot performance. Note that although the minimum spacing is the same as in D configuration, the surface brightness sensitivity and image fidelity to extended structure is considerably inferior to that of the D configuration.
A project with the goal of doubling the longest baseline available in the A configuration by establishing a real-time fiber optic link between the VLA and the VLBA antenna at Pie Town was established in the late 1990s, and used through 2005. This link is no longer operational; there is a goal (unfunded, at present) of implementing a new digital Pie Town link after the EVLA construction project has been completed.
3. Sensitivity
The theoretical thermal noise expected for an image using natural weighting of the visibility data is given by:
[display]\Delta I_m = \frac{SEFD}{\eta_{\rm c}\sqrt{n_{\rm pol}N(N-1)t_{\rm int}\Delta\nu}}[/display] | (1) |
where:
- - SEFD is the "system equivalent flux density" (Jy), defined as the flux density of a radio source that doubles the system temperature. Lower values of the SEFD indicate more sensitive performance. For the VLA's 25-meter paraboloids, the SEFD is given by the equation SEFD = 5.62Tsys/ηA, where Tsys is the total system temperature (receiver plus antenna plus sky), and ηA is the antenna aperture efficiency in the given band.
- - ηc is the correlator efficiency (at least 0.92 for the VLA).
- - npol is the number of polarization products included in the image; npol = 2 for images in Stokes I, Q, U, or V, and npol = 1 for images in 'RCP' or 'LCP'.
- - N is the number of antennas.
- - tint is the total on-source integration time in seconds.
- - Δν is the bandwidth in Hz.
Figure 2 shows the SEFD as a function of frequency for those bands currently installed on VLA antennas, and include the contribution to Tsys from atmospheric emission at the zenith. Table 6 gives the SEFD at some fiducial VLA frequencies.
Frequency | SEFD | RMS confusion level |
---|---|---|
(Jy) | in D config (µJy/beam) | |
1.5 GHz (L) | 420 | 89 |
3.0 GHz (S) | 370 | 14 |
6.0 GHz (C) | 310 | 2.3 |
10.0 GHz (X) | 250 | negligible |
15 GHz (Ku) | 350 | negligible |
22 GHz (K) | 560 | negligible |
33 GHz (Ka) | 730 | negligible |
45 GHz (Q) | 1400 | negligible |
- Note: SEFDs at Ku, K, Ka, and Q bands include contributions from Earth's atmosphere, and were determined under good conditions.
Note that the theoretical rms noise calculated using equation 1 is the best limit possible. There are several factors that will tend to increase the noise compared with theoretical:
- For the more commonly-used "robust" weighting scheme, intermediate between pure natural and pure uniform weightings (available in the AIPS task IMAGR and CASA task clean), typical parameters will result in the sensitivity being a factor of about 1.2 worse than the listed values.
- Confusion. There are two types of confusion: (i) that due to confusing sources within the synthesized beam, which affects low resolution observations the most. Table 6 shows the confusion noise in D configuration (see Condon 2002, ASP Conf. 278, 155), which should be added in quadrature to the thermal noise in estimating expected sensitivities. The confusion limits in C configuration are approximately a factor of 10 less than those in Table 6; (ii) confusion from the sidelobes of uncleaned sources lying outside the image, often from sources in the sidelobes of the primary beam. This primarily affects low frequency observations.
- Weather. The sky and ground temperature contributions to the total system temperature increase with decreasing elevation. This effect is very strong at high frequencies, but is relatively unimportant at the other bands. The extra noise comes directly from atmospheric emission, primarily from water vapor at K-band, and from water vapor and the broad wings of the strong 60 GHz O2 transitions at Q-band.
- Losses from the 3-bit samplers. The VLA's 3-bit samplers incur an additional 10 to 15% loss in sensitivity above that expected -- i.e., the efficiency factor ηc = 0.8 to 0.85.
- In general, the zenith atmospheric opacity to microwave radiation is very low - typically less than 0.01 at L, C and X-bands, 0.05 to 0.2 at K-band, and 0.05 to 0.1 at the lower half of Q-band, rising to 0.3 by 49 GHz. The opacity at K-band displays strong variations with time of day and season, primarily due to the 22 GHz water vapor line. Conditions are best at night, and in the winter. Q-band opacity, dominated by atmospheric O2, is considerably less variable.
- Observers should remember that clouds, especially clouds with large water droplets (read, thunderstorms!), can add appreciable noise to the system temperature. Significant increases in system temperature can, in the worst conditions, be seen at frequencies as low as 5 GHz.
- Tipping scans can be used for deriving the zenith opacity during an observation. In general, tipping scans should only be needed if the calibrator used to set the flux density scale is observed at a significantly different elevation than the range of elevations over which the phase calibrator and target source are observed. However, the antenna tip capability is currently unavailable for the upgraded VLA -- it is hoped that this will be again be available in the next year.
- When the flux density calibrator observations are within the elevation range spanned by the science observing, elevation dependent effects (including both atmospheric opacity and antenna gain dependencies) can be accounted for by fitting an elevation-dependent gain term. See the following item.
- Antenna elevation-dependent gains. The antenna figure degrades at low elevations, leading to diminished forward gain at the shorter wavelengths. The gain-elevation effect is negligible at frequencies below 8 GHz. The antenna gains can be determined by direct measurement of the relative system gain using the AIPS task ELINT on data from a strong calibrator which has been observed over a wide range of elevation. If this is not possible, care should be taken to observe a primary flux calibrator at the same elevation as the target.
The AIPS task INDXR applies standard elevation-dependent gains and an estimated opacity, generated from ground-based weather, in CL table version 1. The CASA calibration tasks (e.g. gaincal, bandpass) also use the standard gain curves.
- Pointing. The SEFD quoted above assumes good pointing. Under calm nighttime conditions, the antenna blind pointing is about 10 arcsec rms. The pointing accuracy in daytime can be much worse -- occasionally exceeding 1 arcminte, due to the effects of solar heating of the antenna structures. Moderate winds have a very strong effect on both pointing and antenna figure. The maximum wind speed recommended for high frequency observing is 15 mph (7 m/s). Wind speeds near the stow limit (45 mph) will have a similar negative effect at 8 and 15 GHz.
- To achieve better pointing, "referenced pointing" is recommended, where a nearby calibrator is observed in interferometer pointing mode every hour or so. The local pointing corrections thus measured can then be applied to subsequent target observations. This reduces rms pointing errors to as little as 2 - 3 arcseconds (but more typically 5 to 7 arcseconds) if the reference source is within about 15 degrees (in azimuth and elevation) of the target source, and the source elevation is less than 70 degrees. At source elevations greater than 80 degrees (zenith angle < 10 degrees), source tracking becomes difficult; it is recommended to avoid such source elevations during the observation preparation setup.
- Use of referenced pointing is highly recommended for all Ku, K, Ka, and Q-band observations, and for lower frequency observations of objects whose total extent is a significant fraction of the antenna primary beam. It is usually recommended that the referenced pointing measurement be made at 8 GHz (X-band), regardless of what band your target observing is at, since X-band is the most sensitive, and the closest calibrator is likely to be weak. Proximity of the reference calibrator to the target source is of paramount importance; ideally the pointing sources should precede the target by 20 or 30 minutes in Right Ascension (RA). The calibrator should have at least 0.3 Jy flux density at X-band and be unresolved on all baselines to ensure an accurate solution.
To aid VLA proposers there is an exposure tool calculator on-line at http://science.nrao.edu/facilities/evla/tools/exposure/evlaExpoCalc.jnlp that provides a graphical user interface to these equations.
The beam-averaged brightness temperature measured by a given array depends on the synthesized beam, and is related to the flux density per beam by:
[display]T_{\rm b} = \frac{S \lambda^2}{2k\Omega} = F \cdot S[/display]
-
-
-
-
- ...equation (2)
-
-
-
where Tb is the brightness temperature (Kelvins) and Ω is the beam solid angle. For natural weighting (where the angular size of the approximately Gaussian beam is ∼ 1.5λ/Bmax), and S in mJy per beam, the constant F depends only upon array configuration and has the approximate value F = 190, 18, 1.7, 0.16 for A, B, C, and D configurations, respectively. The brightness temperature sensitivity can be obtained by substituting the rms noise, ΔIm, for S. Note that Equation 2 is a beam-averaged surface brightness; if a source size can be measured the source size and integrated flux density should be used in Equation 2, and the appropriate value of F calculated. In general the surface brightness sensitivity is also a function of the source structure and how much emission may be filtered out due to the sampling of the interferometer. A more detailed description of the relation between flux density and surface brightness is given in Chapter 7 of Reference 1, listed in Documentation.
For observers interested in HI in galaxies, a number of interest is the sensitivity of the observation to the HI mass. This is given by van Gorkom et al. (1986; AJ, 91, 791):
[display]M_{\rm HI} = 2.36 \times 10^5 D^2 \sum S \Delta V ~M_\odot[/display]
where D is the distance to the galaxy in Mpc, and SΔV is the HI line area in units of Jy km/s.
Also, look at Figure 1.
4. VLA Frequency Bands and Tunability
For observations taken with the 8-bit samplers under the general observing program, each receiver can tune to two different frequencies, each 1.024 GHz wide, within the same frequency band. Right-hand circular (RCP) and left-hand circular (LCP) polarizations are received for both frequencies, except for the low-band receiver (50 -- 500 MHz), which provides linear polarization. Each of these four data streams follows the VLA nomenclature, and are known as IF (for "Intermediate Frequency" channel) "A", "B", "C", and "D". IFs A and B provide RCP, IFs C and D provide LCP. IFs A and C are always at the same frequency, as are IFs B and D (but note that the A and C IFs frequency is usually different from the B and D frequency). We normally refer to these two independent data streams as "IF pairs" -- i.e., the 'A/C' pair and the 'B/D' pair. Currently, a maximum of 1.024 GHz can be correlated for each IF pair (see Correlator Configurations), for a total maximum bandwidth of 2.048 GHz. To distinguish this 8-bit system from the 3-bit system (described below), these IF pairs are denoted A0/C0 and B0/D0.
The tuning ranges, along with default frequencies for continuum applications, are given in Table 7 below.
Band | Range1 | Default frequencies for continuum applications (GHz) | |
---|---|---|---|
(GHz) | IF pair A0/C0 | IF pair B0/D0 | |
4 m (4) | 0.058-0.084 | TBD | TBD |
90 cm (P) | 0.23-0.47 | TBD | TBD |
20 cm (L) | 1.0-2.0 | 1.0 -- 1.52 | 1.5 -- 2.02 |
13 cm (S) | 2.0-4.0 | 2.0 -- 3.0 |
3.0 -- 4.0 |
6 cm (C) | 4.0-8.0 | 4.5 -- 5.5 | 5.5 -- 6.5 |
3 cm (X) | 8.0-12.0 | 8.0 -- 9.0 | 9.0 -- 10.0 |
2 cm (Ku) | 12.0-18.0 | 13.0 -- 14.0 |
14.0 -- 15.0 |
1.3 cm (K) | 18.0-26.5 | 20.2 -- 21.2 |
21.2 -- 22.2 |
1 cm (Ka) | 26.5-40.0 | 32.0 -- 33.0 |
31.0 -- 32.0 |
0.7 cm (Q) | 40.0-50.0 | 40.0 -- 41.0 |
41.0 -- 42.0 |
- Notes:
- 1. Listed here are the nominal band edges. For all bands, the receivers can be tuned to frequencies outside this range, but at the cost of diminished performance. Contact VLA staff for further information.
- 2. The default frequency set-up for L-band comprises two 512 MHz IF pairs (each comprising 8 contiguous subbands of 64 MHz) to cover the entire 1-2 GHz of the L-band receiver.
In general, for all frequency bands except Ka, if the total span of the two independent IFs (defined as the frequency difference between the lower edge of one IF pair and the upper edge of the other) is less than 8.0 GHz, there are no restrictions on the frequency placements of the two IF pairs. For K, Ka and Q bands (the only bands where a span greater than 8 GHz is possible), there are special rules:
- At Ka band, the low frequency edge of the AC IF must be greater than 32.0 GHz. There is no restriction on the BD frequency, unless the BD band overlaps the AC band when the latter is tuned at or near the 32.0 GHz limit. In this case, the Observation Preparation Tool (OPT) may not allow the requested frequency setups. Users wanting to use such a frequency setup are encouraged to contact the NRAO Helpdesk for possible tuning options.
- At K and Q bands, if the frequency span is greater than 8.0 GHz, the BD frequency must be lower than the AC frequency.
With the 3-bit samplers, more options are available. This system provides four (R,L) polarization pairs, each 2048 MHz wide. The A/C IF pair provides two sampled pairs, labelled A1C1 and A2C2, and the B/D IF pair provides two sampled pairs, labelled B1D1 and B2D2. The maximum frequency span permitted for the A1C1 and A2C2 pairs is about 5000 MHz. The same restriction applies to B1D1 and B2D2. The tuning restrictions given above for the separation and location of the 8-bit pairs A0C0 and B0D0 also apply to the 3-bit pairs.
5. VLA Samplers
The VLA is now equipped with two different types of samplers. Which set you should use depends primarily on your science goals.
A) 8-bit Set
This set consists of four 8-bit samplers running at 2.048 GSamp/sec. The four samplers are arranged in two pairs, each pair providing 1024 MHz bandwidth in both polarizations. The two pairs are denoted A0/C0 and B0/D0. Taken together, the four samplers offer a maximum of 2048 MHz coverage with full polarization. The frequency spans sampled by the two pairs need not be adjacent, but some restrictions apply, depending on band. Tuning restrictions for these two pairs are described in section 3.4 (VLA Frequency Bands and Tunability) of this document.
B) 3-bit Set
This set consists of eight 3-bit samplers running at 4.096 GSamp/sec. The eight samplers are arranged as four pairs, each pair providing 2048 MHz bandwidth in both polarizations. Two of these pairs, denoted A1/C1 and A2/C2 cannot span more than 5000 MHz (lower edge of one to the higher edge of the other). The same limitation applies to the second pair, denoted B1/D1 and B2/D2. A more complete description of the tuning restrictions is given in Section 3.4 (VLA Frequency Bands and Tunability) of this document.
C) Major Characteristics of each Set
The 8-bit samplers have been in use for many years, and their characteristics are well understood. Use of this sampler set is warranted for low-band, L-band, and S-band observations, where the full analog bandwidth provided by the receivers is less than or equal to the 2048 MHz span covered by the samplers.
The 3-bit samplers, at the time of writing, are still being implemented on the array. Their characteristics are currently being explored, so users should refer to other sections of this document for updated and current instructions for optimum use. Major issues users need to be aware of include:
- The sensitivity of the array is reduced by 10 to 15% with the 3-bit system, compared to the ideal analog system. For comparison, note that the loss is less than 3% for the 8-bit system.
- Each of the eight 3-bit samplers has a resonance of approximate width 3 MHz. Each sampler's resonance is independent of, and different from, all others, so there is no correlated signal between antennas. The resonance, in effect, greatly reduces the SNR within the narrow frequency span of that resonance. This will degrade bandpass solutions within the narrow resonance window, and images made at frequencies within any one of the resonances will show significant loss in sensitivity. The resonance effect is easily seen in autocorrelation spectra, and it is recommended that users utilize these to inform themselves of the compromised frequencies. As not all samplers are (as of this writing) outfitted on the array, we can not yet provide a table of the compromised frequencies.
- The SY (switched power) values utilized to correct for system gain variations are sensitive to the total power from each antenna. Application of the switched power values will bias the resulting visibilities by a value of 5 to 10%. This origin of this effect is well understood, but we have not yet determined how best to compensate for it. Because of this, we do not recommend use of the SY table data for data taken with the 3-bit samplers.
D) Setting up the 3-bit or 8-bit Samplers
Setting up either set of samplers requires a short initial observation for each individual LO (frequency) tuning. For the 8-bit system, a scan of 1 minute duration is sufficient for each tuning. The pointing direction of the antennas is not critical.
For the 3-bit system, the requirements are more demanding. The recommended duration is two minutes per frequency tuning. In addition, we recommend that the antennas be pointed at an elevation near to 30 degrees. A slew time of up to three minutes will likely be required to position the antennas at this elevation.
E) Which set to use?
Observations at S, L, and low bands should use the 8-bit sampler set in all cases. There is no benefit from use of the 3-bit samplers at these bands.
Continuum-application observations at Ku, K, Ka, and Q bands, where the best sensitivity is critical, should use the 3-bit sampler sets. Similarly, wide-band spectral line searches requiring more than the 2 GHz span of the 8-bit set should use the 3-bit sampler sets.
Spectral-line application observations whose frequency span lies within the two 1024 MHz-wide spans offered by the 8-bit set should use the 8-bit set.
It is not known at this time if full-band continuum observations at C and X bands, where the maximum bandwidth is 4 GHz, should use the 8-bit or 3-bit sets. In principle, the 3-bit system should offer ~40% better continuum sensitivity, but this will be eroded by:
- the reduction of 10 to 15% noted above,
- the increased setup time needed for 3-bit observations,
- the extra RFI present at these bands which can be avoided by tuning the 8-bit system. This is especially notable at C-band.
F) Other issues to be aware of.
The 8-bit system provides much better protection against strong RFI. However, we do not yet know if the 3-bit system will be
significantly compromised by RFI at any band where its use is warranted. This is most likely to be an issue at C-band, where
strong RFI from broadcast satellites and local microwave links are known to be present. Our current understanding of the RFI spectrum at higher frequency bands suggests there will be no problems with the 3-bit system, except possibly when pointing near geostationary satellites.
Polarization testing conducted so far indicates no degradation of performance by using the 3-bit samplers.
6. Field of View
6.1. Introduction
At least four different effects will limit the field of view. These are: primary beam; chromatic aberration; time-averaging; and non-coplanar baselines. We discuss each briefly:
6.2. Primary Beam
The ultimate factor limiting the field of view is the diffraction-limited response of the individual antennas. An approximate formula for the full width at half power in arcminutes is: θPB = 45/νGHz. More precise measurements of the primary beam shape have been derived and are incorporated in AIPS (task PBCOR) and CASA (clean task and the imaging toolkit) to allow for correction of the primary beam attenuation in wide-field images. Objects larger than approximately half this angle cannot be directly observed by the array. However, a technique known as "mosaicing," in which many different pointings are taken, can be used to construct images of larger fields. Refer to References 1 and 2 in Documentation for details.
Guidelines for mosaicing with the EVLA are given at http://evlaguides.nrao.edu/index.php?title=Mosaic_Observing
6.3. Chromatic Aberration (Bandwidth Smearing)
The principles upon which synthesis imaging are based are strictly valid only for monochromatic radiation. When visibilities from a finite bandwidth are gridded as if monochromatic, aberrations in the image will result. These take the form of radial smearing which worsens with increased distance from the delay-tracking center. The peak response to a point source simultaneously declines in a way that keeps the integrated flux density constant. The net effect is a radial degradation in the resolution and sensitivity of the array.
These effects can be parameterized by the product of the fractional bandwidth (Δν/ν0) with the source offset in synthesized beamwidths (θ0/θHPBW). Table 8 shows the decrease in peak response and the increase in apparent radial width as a function of this parameter. Table 8 should be used to determine how much spectral averaging can be tolerated when imaging a particular field.
(Δν/ν0)*(θ0/θHPBW) | Peak | Width |
---|---|---|
0.0 | 1.00 | 1.00 |
0.50 | 0.95 | 1.05 |
0.75 | 0.90 | 1.11 |
1.0 | 0.80 | 1.25 |
2.0 | 0.50 | 2.00 |
- Note: The reduction in peak response and increase in width of an object due to bandwidth smearing (chromatic aberration). Δν/ν0 is the fractional bandwidth; θ0/θHPBW is the source offset from the phase tracking center in units of the synthesized beam.
6.4. Time-Averaging Loss
The sampled coherence function (visibility) for objects not located at the phase-tracking center is slowly time-variable due to the motion of the source through the interferometer coherence pattern, so that averaging the samples in time will cause a loss of amplitude. Unlike the bandwidth loss effect described above, the losses due to time averaging cannot be simply parametrized, except for observations at δ = 90°. In this case, the effects are identical to the bandwidth effect except they operate in the azimuthal, rather than the radial, direction. The functional dependence is the same as for chromatic aberration with Δν/ν0 replaced by ωeΔtint, where ωe is the Earth's angular rotation rate, and Δtint is the averaging interval.
For other declinations, the effects are more complicated and approximate methods of analysis must be employed. Chapter 13 of Reference 1 (in Documentation) considers the average reduction in image amplitude due to finite time averaging. The results are summarized in Table 9, showing the time averaging in seconds which results in 1%, 5% and 10% loss in the amplitude of a point source located at the first null of the primary beam. These results can be extended to objects at other distances from the phase tracking center by noting that the loss in amplitude scales with (θΔtint)2, where θ is the distance from the phase center and Δtint is the averaging time. We recommend that observers reduce the effect of time-average smearing by using integration times as short as 1 or 2 seconds (also see the section on Time Resolution and Data Rates) in the A and B configurations.
Amplitude loss | |||
---|---|---|---|
Configuration | 1.0% | 5.0% | 10.0% |
A | 2.1 | 4.8 | 6.7 |
B | 6.8 | 15.0 | 21.0 |
C | 21.0 | 48.0 | 67.0 |
D | 68.0 | 150.0 | 210.0 |
- Note: The averaging time (in seconds) resulting in the listed amplitude losses for a point source at the antenna first null. Multiply the tabulated averaging times by 2.4 to get the amplitude loss at the half-power point of the primary beam. Divide the tabulated values by 4 if interested in the amplitude loss at the first null for the longest baselines.
6.5. Non-Coplanar Baselines
The procedures by which nearly all images are made in Fourier synthesis imaging are based on the assumption that all the coherence measurements are made in a plane. This is strictly true for E-W interferometers, but is false for the EVLA, with the single exception of snapshots. Analysis of the problem shows that the errors associated with the assumption of a planar array increase quadratically with angle from the phase-tracking center. Serious errors result if the product of the angular offset in radians times the angular offset in synthesized beams exceeds unity: θ > λB/D2, where B is the baseline length, D is the antenna diameter, and λ is the wavelength, all in the same units. This effect is most noticeable at λ90 and λ20 cm in the larger configurations, but will be notable in wide-field, high fidelity imaging for other bands and configurations.
Solutions to the problem of imaging wide-field data taken with non-coplanar arrays are well known, and have been implemented in AIPS (IMAGR) and CASA (clean). Refer to the package help files for these tasks, or consult with Rick Perley, Frazer Owen, or Sanjay Bhatnagar for advice. More computationally efficient imaging with non-coplanar baselines is being investigated, such as the "W-projection" method available in CASA; see EVLA Memo 67 (http://www.aoc.nrao.edu/evla/geninfo/memoseries/evlamemo67.pdf) for more details.
7. Time Resolution and Data Rates
The default integration times for the various array configurations are as follows:
Configuration | Observing Bands | Default integration time |
---|---|---|
D, C |
L S C |
5 seconds |
D, C |
X Ku K Ka Q |
3 seconds |
B | all | 3 seconds |
A | all | 1 second |
Observations with the 3-bit (wideband) samplers must use these integration times. Observations with the 8-bit samplers may use shorter integration times, but these must be requested and justified explicitly in the proposal, and obey the following restrictions:
Minimum Integration Times and Maximum Data Rates (2013A) |
||
---|---|---|
Proposal type |
Minimum integration time |
Maximum data rate |
Standard | 1 sec |
20 MB/s ( ~70 GB/hr) |
Shared risk | 50 msec | 60 MB/s (~210 GB/hr) |
The maximum recommended integration time for any EVLA observing is 60 seconds. For high frequency observations with short scans (e.g., fast switching, as described in Rapid Phase Calibration and the Atmospheric Phase Interferometer (API)), shorter integration times may be preferable.
Observers should bear in mind the data rate of the VLA when planning their observations. For Nant antennas and integration time Δt, the data rate is:
Data rate | = 45 MB/sec × (Nchpol/16384) x Nant × (Nant1)/(27×26) / (Δt/1 sec) |
= 160 GB/hr × (Nchpol/16384) x Nant × (Nant1)/(27×26) / (Δt/1 sec) | |
= 3.7 TB/day × (Nchpol/16384) x Nant × (Nant1)/(27×26) / (Δt/1 sec) | |
...equation (4) |
Here Nchpol is the sum over all subbands of spectral channels times polarization products:
Nchpol = Σsb Nchan,i x Npolprod,i |
where Nchan,i is the number of spectral channels in subband i, and Npolprod,i is the number of polarization products for subband i (1 for single polarization [RR or LL], 2 for dual polarization [RR & LL], 4 for full polarization products [RR, RL, LR, LL]). This formula, combined with the maximum data rates given above, imply that observations using the maximum number of channels currently available (16384) will be limited to minimum integration times of ~3 seconds for standard observations, and 0.8 seconds for shared risk observations.
These data rates are challenging for data transfer, as well as data analysis. Currently data may either be downloaded via ftp within the Science Operations Centers, or mailed on hard drives for those not in the same building as the archive. The Archive Access Tool (https://science.nrao.edu/facilities/evla/data-archive/evla) allows some level of frequency averaging to decrease data set sizes before ftp, for users whose science permits; note that the full spectral resolution will be retained in the NRAO archive for all observations.
Higher time resolutions and data rates are possible in principle but will be considered only through the resident shared risk program.
8. Radio-Frequency Interference
The very wide bandwidths of the upgraded Very Large Array mean that RFI (radio-frequency interference) will be present in a far larger fraction of VLA observations than in observations made with the old systems. Considerable effort has gone into making the VLA's new electronics as linear as possible, so that the effects of any RFI will remain limited to the actual frequencies at which the RFI exists. Non-linear effects, such as receiver saturation, should occur only for those very unlikely, and usually very brief, times when the emitter is within the antenna primary beam.
RFI is primarily a problem within the low frequency bands (C, S, L, and the low-band system), and is most serious to the D configuration. With increasing frequency and increasing resolution comes an increasing fringe rate, which is often very effective in reducing interference to tolerable levels.
The bands within the tuning range of the VLA which are allocated exclusively to radio astronomy are 1400-1427 MHz, 1660-1670 MHz, 2690-2700 MHz, 4990-5000 MHz, 10.68-10.7 GHz, 15.35-15.4 GHz, 22.21-22.5 GHz, 23.6-24.0 GHz, 31.3-31.8 GHz, and 42.5-43.5 GHz. No external interference should occur within these bands.
RFI seen in VLA data can be internal or external. Great effort has been expended to eliminate all internally-generated RFI. Nevertheless, some internal RFI remains, which we are working hard to eliminate. Nearly all such internally-generated signals are at multiples of 128 MHz. So far as we know, all such internal signals are unresolved in frequency, and hence will affect only a single channel.
Radio frequency interference of external origin will be an increasing problem to astronomical observations. Table 10 lists some of the sources of external RFI at the VLA site that might be observed within the VLA's expanded tuning range within L and S bands. Figure 3 shows a raw power cross-power spectrum at L-band. Figure 4 shows a similar plot for the lower half of S-band.
Frequency (MHz) | Source | Comments |
---|---|---|
1025-1150 | Aircraft navigation | Very strong |
1200.0 | VLA modem | |
1217-1237 | GPS L2 | Very strong |
1243-1251 | GLONASS L2 | |
1254 | Aeronautical radar | |
1263 | Aeronautical radar | |
1268 | COMPASS E6 | |
1310 | Aeronautical radar | |
1317 | Aeronautical radar | |
1330 | Aeronautical radar | |
1337 | Aeronautical radar | |
1376-1386 | GPS L3 | Intermittent |
1525-1564 | INMARSAT satellites | |
1564-1584 | GPS L1 | Very strong |
1598-1609 | GLONASS L1 | |
1618-1627 | IRIDIUM satellites | |
1642 | 2nd harmonic VLA radios | Sporadic |
1683-1687 | GOES weather satellite | |
1689-1693 | GOES weather satellite | |
1700-1702 | NOAA weather satellite | |
1705-1709 | NOAA weather satellite | |
1930-1990 | PCS cell phone base stations | |
2178-2195 | Satellite Downlink | very strong |
2320-2350 | Satellite radio | very strong |
VLA staff periodically observes the entire radio spectrum, with the VLA, from 1.0 through 50.0 GHz with 125 kHz channel resolution to monitor the ever-changing RFI spectrum. Plots from this program, accompanied with tables of identified sources are available online, at http://science.nrao.edu/facilities/evla/observing/RFI/index. Users concerned about the precise frequencies of strong RFI, and the likelihood of being affected, are encouraged to peruse these plots.
Although most of the stronger sources of RFI are always present, it is very difficult to reliably predict their effect on observations. Besides the already noted dependence on frequency and array configuration, there is another significant dependency on sky location for those satellites in geostationary orbit. For these transmitters, (for example, the frequency range from 3.8 to 4.2 GHz), the effect on observing varies dramatically on the declination of the target source. Sources near zero declination will be very strongly affected, while observations north of the zenith may well be nearly unaffected, especially at the highest resolutions.
Also available are total-power plots of all RFI observations made by the interference protection group, from 1993 onwards at http://www.vla.nrao.edu/cgi-bin/rfi.cgi. For general information about the RFI environment, contact the head of the IPG (Interference Protection Group) by sending e-mail to nrao-rfi@nrao.edu.
The VLA electronics (including the WIDAR correlator) have been designed to minimize gain compression due to very strong RFI signals, so that in general it is possible to observe in spectral windows containing RFI, provided the spectra are well sampled to constrain Gibbs ringing, and spectral smoothing (such as Hanning) is applied. Both AIPS and CASA provide useful tasks which automatically detect and flag spectral channels/times which contain strong RFI.
Extracting astronomy data from frequency channels in which the RFI is present is much more difficult. Testing of algorithms which can distinguish and subtract RFI signals from interferometer data is ongoing.
The new 3-bit samplers will be more susceptible to RFI signals than the 8-bit samplers, since the latter have more 'levels' within which these strong signals can be accommodated. However, since the RFI power at the bands where the 3-bit samplers will most commonly be used (X, Ku, K, Ka, and Q) is nearly always less than the total noise power, we do not expect problems when wide-band 3-bit observing is done in these bands. At C-band however, there are strong TV and microwave communications signals which may saturate the 3-bit samplers. Until VLA staff can better measure this effect, we recommend use of the 8-bit samplers at C-band.
Calibration of VLA data when strong RFI is present within a subband can be difficult. Careful editing of the data, using newly available programs within CASA and AIPS, will be necessary before sensible calibration can be done. The use of spectral smoothing (typically, Hanning), prior to editing and calibration, is strongly recommended when RFI is present within a subband.
Identification and removal of RFI is always more effective when the spectral and temporal resolutions are high. However, the cost of higher spectral and temporal resolution is in database size and, especially, in computing time. A good strategy is to observe with high resolution, then average down in time and frequency once the editing is completed.
9. Subarrays
The separation of the VLA into up to three subarrays will be supported for general observing in semester 2013A, Each subarray is independent for the entire duration of an observation, and is set up using an independent scheduling block (SB). All SBs must be the same length, and will be scheduled simultaneously. The antennas comprising each subarray will remain in the same subarray for the duration of the observation. There are some limitations on the number of antennas that can comprise each subarray, due to the architecture of the correlator. The correlator operates on 8 groups of 4 antennas (up to 32 antennas total), and each group can only contain antennas from one subarray. Thus, for three subarrays and 27 antennas, antennas may be distributed as 8, 9, 10 antennas, or 11, 11, 5 antennas, but not 9, 9, 9 antennas since each subarray of 9 antennas requires 3 groups, and the three subarrays together would require 3x3=9 groups, where only eight are available. For two subarrays the grouping requirement gives no restrictions, and the antennas may be flexibly allocated between the two subarrays.
Subarrays will only be available for general observing with the 8-bit samplers and 2 GHz bandwidth. Access to the 3-bit samplers and/or more flexible correlator set-ups will be offered under the shared risk observing programs.
10. Positional Accuracy
Summary: A target's position can be determined to a small fraction of the synthesized beam, limited by atmospheric phase stability, the proximity of an astrometric calibrator, the calibrator-source cycle time, and the SNR on target.
In preparation for observing, the a-priori position of a target must be known to within the antenna primary beam, except perhaps for mosaicing. In post-processing, target positions are typically determined from an image made after phase calibration, i.e. zeroing of antenna phases on the reference source. (Note that phase self-calibration makes the position indeterminate.)
It may help to think of astrometry in 2 steps, narrow and wide-field. In the former, the target is close to the phase tracking center and the antennas nod every few minutes between the target and a calibrator of known position. Under good conditions of phase stability, accurate antenna positions, (so-called 'baselines'), a strong target, a close calibrator of known position, and rapid switching, the accuracy can approach 1-2% of the synthesized beam, with a floor of ~2 mas. Under more typical conditions, 10% of the beam is readily achieved.
Astrometric calibrators are marked 'A' in the VLA calibrator list, and have an accuracy of ~2 mas. Other catalogs from the USNO and the VLBA are also useful, but offsets may exist between the VLA and VLBA centroids, arising from extended structure in the particular source, and the different resolutions of the arrays. Phase stability can be assessed in real time from the Atmospheric Phase Interferometer (API) at the VLA site, which uses observations of a geostationary satellite at ~12GHz.
For studies of proper motion and parallax, the absolute accuracy of a calibrator may be less important than its stability over time. Close or in-beam calibrators with poor a-priori positions can be used, and tied to the ICRF reference frame in the same or separate observations.
The widefield case is to determine the positions of targets within the primary beam, referenced to a calibrator within the beam or close by. In addition to the previous effects, there are distortions as a function of position in the field, from small errors in the Earth orientation parameters (EOP) used at correlation time, differential aberration, and phase gradients across the primary beam. These effects are handled somewhat differently in the various reduction packages. With no special effort, the errors build up to roughly ~1 synthesized beam at a separation of ~10^4 beams from the phase tracking center. Not all these errors are fully understood, and accurate recovery of positions over the full primary beam in the wideband, widefield case is a research area.
11. Limitations on Imaging Performance
11.1. Introduction
Imaging performance can be limited in many different ways. Some of the most common are listed in the following subsections.
11.2. Image Fidelity
With conventional external calibration methods, even under the best observing conditions, the achieved dynamic range will rarely exceed a few hundred. The limiting factor is most often the atmospheric phase stability, although pointing errors and changes in atmospheric opacity can also be a limiting factor. If the target source contains compact structures of sufficient strength (depending on the band, bandwidth, atmospheric coherence time, and source complexity), self-calibration can be counted on to improve the images. Dynamic ranges in the thousands to hundreds of thousands can be achieved using these techniques, depending on the underlying nature of the errors. With the new WIDAR correlator and its much greater bandwidths and much higher sensitivities, self-calibration methods can be extended to observations of sources with much lower flux densities than very possible with the old VLA.
11.3. Invisible Structures
An interferometric array acts as a spatial filter, so that for any given configuration, structures on a scale larger than the fringe spacing of the shortest baseline will be completely absent. Diagnostics of this effect include negative bowls around extended objects, and large-scale stripes in the image. Table 5 gives the largest scale visible to each configuration/band combination.
11.4. Poorly Sampled Fourier Plane
Unmeasured Fourier components are assigned values by the deconvolution algorithm. While this often works well, sometimes it fails noticeably. The symptoms depend upon the actual deconvolution algorithm used. For the CLEAN algorithm, the tell-tale sign is a fine mottling on the scale of the synthesized beam, which sometimes even organizes itself into coherent stripes. Further details are to be found in Reference 1 in Documentation.
11.5. Sidelobes from Confusing Sources
At the lower frequencies, large numbers of detectable background sources are located throughout the primary antenna beam, and into its first sidelobe. Sidelobes from those sources which have not been deconvolved will lower the image quality of the target source. Although bandwidth smearing and time-averaging will tend to reduce the effects of these sources, the very best images will require careful imaging of all significant background sources. The deconvolution tasks in AIPS (IMAGR) and CASA (clean) are well suited to this task.
11.6. Sidelobes from Strong Sources
An extension of the previous section is to very strong sources located anywhere in the sky, such as the Sun (especially when a flare is active), or when observing with a few tens of degrees of the very strong sources Cygnus A and Casseopeia A. Image degradation is especially notable at lower frequencies, shorter configurations, and when using narrow-bandwidth observations (especially in spectral line work) where chromatic aberration cannot be utilized to reduce the disturbances. In general, the only relief is to include the disturbing sources in the imaging, or to observe when these objects are not in the viewable hemisphere.
12. Calibration and Flux Density Scale
The VLA Calibrator List contains information on 1860 sources sufficiently unresolved and bright to permit their use as calibrators. The list is available within the Observation Preparation Tool and may be accessed on the Web at http://www.vla.nrao.edu/astro/calib/manual/.
Accurate flux densities can be obtained by observing one of 3C286, 3C147, 3C48 or 3C138 during the observing run. Not all of these are suitable for every observing band and configuration - consult the VLA Calibrator Manual for advice. Over the last several years, we have implemented accurate source models directly in AIPS and CASA for much improved calibration of the amplitude scales. Models are available for 3C48, 3C138, 3C147, and 3C286 for L, C, X, Ku, K, and Q bands. At Ka band either of the K or Q band models works reasonably well. For S-band, use the L or C band models.
Since the standard source flux densities are slowly variable, we monitor their flux densities when the array is in its D configuration. As the VLA cannot accurately measure absolute flux densities, the values obtained must be referenced to assumed or calculated standards, as described in the next paragraph. Table 11 shows the flux densities of these sources in January 2012 at the standard VLA bands. The flux density scale for the VLA, from 1 through 50 GHz, is based on emission models of the planet Mars, which is then calibrated to the CMB dipole using WMAP (Wilkinson Microwave Anisotropy Probe) observations. The source 3C286 (=J1331+3030) is known to be non-variable, and has thus been adopted as the prime flux density calibrator source for the VLA. The adopted polynomial expression for the spectral flux density for 3C286 is:
log(S) = 1.2553 - 0.4689 log(f) - 0.1597 log2(f) + 0.0286 log3(f)
where S is the flux density in Jy, and f is the frequency in GHz.
The absolute accuracy of our flux density scale is estimated to be about 3%. However, with care, the internal accuracy in flux density bootstrapping is better than 1% at all bands except Q-band, where pointing errors limit bootstrap accuracy to perhaps 3%. Note that such high internal accuracies are only possible in long-duration observations where the antenna gains curves and atmospheric opacity can be directly measured, and where there is good elevation overlap between the target source(s) and the flux density standard calibrator.
Source | 1465 | 2565 | 4885 | 8435 | 14965 | 22460 | 36435 | 43340 |
---|---|---|---|---|---|---|---|---|
3C48 = J0137+3309 | 15.56 | 9.80 | 5.39 | 3.14 | 1.77 | 1.19 | 0.73 | 0.63 |
3C138 = J0521+1638 | 8.71 | 6.17 | 4.02 | 2.78 | 1.89 | 1.46 | 1.03 | 0.92 |
3C147 = J0542+4951 | 21.85 | 13.75 | 7.59 | 4.49 | 2.59 | 1.77 | 1.10 | 0.94 |
3C286 = J1331+3030 | 14.90 | 10.03 | 7.34 | 5.09 | 3.39 | 2.52 | 1.75 | 1.53 |
3C295 = J1411+5212 | 22.15 | 12.95 | 6.41 | 3.34 | 1.62 | 0.957 | 0.507 | 0.403 |
NGC7027 | 1.62 | 3.59 | 5.38 | 5.79 | 5.62 | 5.42 | 5.18 | 5.04 |
The sources 3C48, 3C147, and 3C138 are all slowly variable. VLA staff monitor these variations on timescale of a year or two, and suitable polynomial coefficients are determined for them which should allow accurate flux density bootstrapping. These coefficients are updated approximately every other year, and are used in the AIPS task SETJY and in the CASA task setjy.
The VLA antennas have elevation-dependent gain variations which are important to account for at the four highest-frequency bands. Gain curves are determined by VLA staff approximately every other year, and the necessary corrections are applied to the visibility data when these data are downloaded from the archive. In addition to this, atmospheric opacity will also cause an elevation-dependent gain which is particularly notable at these four highest frequency bands. At the current time, we do not have an atmospheric opacity monitoring procedure, so users should utilize the appropriate tasks available in both AIPS and CASA to estimate and correct for the opacity using ground-based weather data. Correction of these gain dependencies, plus regular calibration using a nearby phase calibrator, should enable good amplitude gain calibration for most users. Note that extraordinary attenuation by clouds can only be (approximately) corrected for by regular observation of a nearby calibrator.
A better procedure for removing elevation gain dependencies uses the AIPS task ELINT. This task will generate a 2nd order polynomial gain correction utilizing your own calibrator observations. This will remove both the antenna and opacity gain variations, and has the decided advantage of not utilizing opacity models or possibly outdated antenna gain curves. Use of this procedure is only practical if your observations span a wide range in elevation.
By far the most important gain variation effect is that due to pointing. Daytime observations on sunny days can suffer pointing errors of up to one arcminute (primarily in elevation). This effect can be largely removed by utilizing the 'referenced pointing' procedure. This determines the pointing offset of a nearby calibrator, which is then applied to subsequent target source observations. It is recommended that this local offset be determined at least hourly, utilizing an object within 15 degrees of the target source -- preferentially at an earlier HA. Studies show that the maximum pointing error will be reduced to about 7 arcseconds, or better. VLA staff continue to work on improving this essential methodology.
Although the VLA's new electronics are very stable, gain variations in the system will occur, either due to changes in amplifier gains, or due to changes in internal attenuator settings made by the observing system to ensure that correct voltage levels are provided to the samplers. These changes are monitored by an internal calibration signal, the results of which are stored in the switched power table (SY table). For the most accurate flux density bootstrapping, this table must be applied to the visibility data before calibration. Gain bootstrapping better than 1% can be accomplished for the 8-bit sampler system after application of the SY table data. For the 3-bit system there is an additional complication, as the values of the SY data are sensitive to the total power, as well as the system gain. VLA staff are currently working on a methodology to remove the total power dependency. Not applying the SY table data will reduce bootstrapping accuracy to perhaps 5 to 10%.
13. Complex Gain Calibration
13.1. General Guidelines for Gain Calibration
Adequate gain calibration is a complicated function of source-calibrator separation, frequency, array scale, and weather. And, since what defines adequate for some experiments is completely inadequate for others, it is impossible to define any simple guidelines to ensure adequate phase calibration in general. However, some general statements remain valid most of the time. These are given below.
- Tropospheric effects dominate at wavelengths shorter than 20 cm, ionospheric effects dominate at wavelengths longer than 20 cm.
- Atmospheric (troposphere and ionosphere) effects are nearly always unimportant in the C and D configurations at L and S bands, and in the D configuration at X and C bands. Hence, for these cases, calibration need only be done to track instrumental changes - once per hour is generally sufficient.
- If your target object has sufficient flux density to permit phase self-calibration, there is no need to calibrate more than once hourly at low frequencies (L/S/C bands) or 15 minutes at high frequencies (K/Ka/Q bands) in order to track pointing or other effects that might influence the amplitude scale. The newly-enhanced sensitivity of the VLA now guarantees, for full-band continuum observations, that every field will have enough background sources to enable phase self-calibration at L and S bands, and probably also at C-band. At higher frequencies, the background sky is not sufficient, and only the flux of the target source itself will be available.
- The smaller the source-calibrator angular separation, the better. In deciding between a nearby calibrator with an "S" code in the calibrator database, and a more distant calibrator with a "P" code, the nearby calibrator is usually the better choice (see http://www.vla.nrao.edu/astro/calib/manual/key.html for a description of calibrator codes).
- In clear and calm conditions, day-time phase stability often deteriorates dramatically after about 10AM, due to small-scale convective cells set up by solar heating. Observers should consider a more rapid calibration cycle for observations between this time and a couple hours after sundown.
- At high frequencies, and longer configurations, rapid switching between the source and nearby calibrator is often helpful. See Rapid Phase Calibration and the Atmospheric Phase Interferometer (API).
13.2. Rapid Phase Calibration and the Atmospheric Phase Interferometer (API)
For some objects, and under suitable weather conditions, the phase calibration can be considerably improved by rapidly switching between the source and calibrator. Source-Calibrator observing cycles as short as 40 seconds can be used. However, observing efficiency declines for very short cycle times, so it is important to balance this loss against a realistic estimate of the possible gain. Experience has shown that cycle times of 100 to 150 seconds at high frequencies have been effective for source-calibrator separations of less than 10 degrees. For the old VLA this was known as "fast-switching." For the upgraded VLA it is just a loop of source-calibrator scans with short scan length. This technique "stops" tropospheric phase variations at an effective baseline length of ∼vat/2 where va is the atmospheric wind velocity aloft (typically 10 to 15 m/sec), and t is the total switching time. It has been demonstrated to result in images of faint sources with diffraction-limited spatial resolution on the longest EVLA baselines. Under average weather conditions, and using a 120 second cycle time, the residual phase at 43 GHz should be reduced to ≤ 30 degrees. Further details can be found in VLA Scientific Memos # 169 and 173 (http://www.vla.nrao.edu/memos/sci/). These memos, and other useful information, can be obtained from Reference 12 in Documentation.
Note that the fast switching technique will not work in bad weather (such as rain showers, or when there are well-developed convection cells - most notably, thunderstorms).
An Atmospheric Phase Interferometer (API) is used to continuously measure the tropospheric contribution to the interferometric phase using an interferometer comprising two 1.5 meter antennas separated by 300 meters, observing an 11.7 GHz beacon from a geostationary satellite. The API data can be used to estimate the required calibration cycle times when using fast switching phase calibration, and in the worst case, to indicate to the observer that high frequency observing may not be possible with current weather conditions. A detailed description of the API may be found at http://www.vla.nrao.edu/astro/guides/api/.
Plots of current/historical data can be found at: https://www.aoc.nrao.edu/cgi-bin/weather/apipg.cgi
Characteristic seasonal averages are represented below:
Month | API (night) [deg] | API (median) [deg] | API (day) [deg] | Wind (night) [m/s] | Wind (median) [m/s] | Wind (day) [m/s] |
---|---|---|---|---|---|---|
January | 2.3 | 2.8 | 3.6 | 1.6 | 1.9 | 2.3 |
February | 2.9 | 3.4 | 4.5 | 4.0 | 4.3 | 4.5 |
March | 2.8 | 3.7 | 5.5 | 3.4 | 3.9 | 4.7 |
April | 3.3 | 4.5 | 6.2 | 5.3 | 5.5 | 5.8 |
May | 2.9 | 4.6 | 6.7 | 2.6 | 3.2 | 3.7 |
June | 3.8 | 5.5 | 7.4 | 2.5 | 3.9 | 6.3 |
July | 6.2 | 8.3 | 10.5 | 2.9 | 2.9 | 3.0 |
August | 5.4 | 7.1 | 11.3 | 1.7 | 2.3 | 3.0 |
September | 5.2 | 6.6 | 8.8 | 2.3 | 3.0 | 3.6 |
October | 4.2 | 5.3 | 7.4 | 2.3 | 2.9 | 3.7 |
November | 2.6 | 3.0 | 4.0 | 1.2 | 2.5 | 1.6 |
December | 2.8 | 3.2 | 4.1 | 1.2 | 1.6 | 2.7 |
Click on the Month links above to see plots of phase and wind speed versus time.
- Note: day indicates sunrise to sunset values; night indicates sunset to sunrise values.
14. Polarization
For projects requiring imaging in Stokes Q and U, the instrumental polarization should be determined through observations of a bright calibrator source spread over a range in parallactic angle. The phase calibrator chosen for the observations can also double as a polarization calibrator provided it is at a declination where it moves through enough parallactic angle during the observation (roughly Dec 15deg to 50deg for a few hour track). The minimum condition that will enable accurate polarization calibration from a polarized source (in particular with unknown polarization) is three observations of a bright source spanning at least 60 degrees in parallactic angle (if possible schedule four scans in case one is lost). If a bright unpolarized unresolved source is available (and known to have very low polarization) then a single scan will suffice to determine the leakage terms. The accuracy of polarization calibration is generally better than 0.5% for objects small compared to the antenna beam size. At least one observation of 3C286 or 3C138 is required to fix the absolute position angle of polarized emission. 3C48 also can be used to fix the position angle at wavelengths of 6 cm or shorter. The results of a careful monitoring program of these and other polarization calibrators can be found at http://www.aoc.nrao.edu/~smyers/evlapolcal/polcal_master.html
High sensitivity linear polarization imaging may be limited by time dependent instrumental polarization, which can add low levels of spurious polarization near features seen in total intensity and can scatter flux throughout the polarization image, potentially limiting the dynamic range. Preliminary investigation of the EVLA's new polarizers indicates that these are extremely stable over the duration of any single observation, strongly suggesting that high quality polarimetry over the full bandwidth will be possible.
The accuracy of wide field linear polarization imaging will be limited, likely at the level of a few percent at the antenna half-power width, by angular variations in the antenna polarization response. Algorithms to enable removable of this angle-dependent polarization are being tested, and observations to determine the antenna polarizations have begun. Circular polarization measurements will be limited by the beam squint, due to the offset secondary focus feeds, which separates the RCP and LCP beams by a few percent of the FWHM. The same algorithms noted above to correct for antenna-induced linear polarization can be applied to correct for the circular beam squint. Measurement of the beam squints, and testing of the algorithms, is ongoing.
Ionospheric Faraday rotation of the astronomical signal is always notable at 20 cm. The typical daily maximum rotation measure under quiet solar conditions is 1 or 2 radians/m2, so the ionospherically-induced rotation of the plane of polarization at these bands is not excessive - 5 degrees at 20 cm. However, under active conditions, this rotation can be many times larger, sufficiently large that polarimetry is impossible at 20 cm with corrrection for this effect. The AIPS program TECOR has been shown to be quite effective in removing large-scale ionospherically induced Faraday Rotation. It uses currently-available data in IONEX format. Please consult the TECOR help file for detailed information. In addition, the interim EVLA receivers generally have poor polarization performance outside the frequency range previously covered by the VLA (e.g., outside the 4.5-5.0 GHz frequency range for C band, and outside 1.3-1.7 GHz for L-band), and the wider frequency bands of these interim receivers may be useful only for total intensity measurements.
For more details see: https://science.nrao.edu/facilities/evla/early-science/polarimetry
15. Correlator Configurations
15.1. Introduction
The VLA correlator is very flexible, and can provide data in many ways. In the 2013A call for proposals we are offering a greatly expanded set of general capabilities, as well as some extensions as (non-resident) shared risk observing.
15.2. General Capabilities
15.2.1. Introduction
The correlator configurations offered for general observing in the 2013A semester may be divided into three basic modes: wideband, spectral line, and subarrays. Note that the possible setups are also subject to the integration time and data rate restrictions outlined in the section on Time Resolution and Data Rates. The possibilities and restrictions are embodied in the General Observing Setup Tool (GOST), which must be used to define the correlator configuration for general observing proposals for the 2013A semester. We are in the process of updating the Resource Configuration Tool (RCT) to reflect the new possibilities.
Note that phased array configurations are only allowed as part of VLBI experiments (see the section on VLBI Observations) or as Resident Shared Risk observations.
15.2.2. Wideband Observing
The wideband observing setup provides the widest possible bandwidth for a given observing band, with channel spacing depending on the number of polarization products.
Wideband & Subarray Correlator Options | |
---|---|
Polarization products | Channel spacing |
Full (RR, RL, LR, LL) | 2 MHz |
Dual (RR, LL) | 1 MHz |
Single (RR or LL) | 0.5 MHz |
For the low-frequency bands (L, S, C, X, and Ku bands) this uses the 8-bit samplers to provide a total of 2 GHz of bandwidth per polarization. The three highest-frequency bands (K, Ka, and Q bands) employ the new 3-bit (wideband) samplers, providing a total of 8 GHz of bandwidth; in this case all subbands are 128 MHz wide. In many frequency bands the total processed bandwidth is less than that delivered by the front-end. In those cases the observer may independently tune two 1 GHz baseband pairs when using the 8-bit samplers, or four 2 GHz baseband pairs when using the 3-bit samplers. The (fairly mild) tuning restrictions are described in the section on VLA Frequency Bands and Tunability.
15.2.3. Spectral Line Configurations
Starting with semester 2013A observers have access to very flexible subbands within two 1 GHz baseband pairs (using the 8-bit samplers). These capabilities may be summarized as follows:
- Two 1 GHz baseband pairs (using the 8-bit samplers), independently tunable within the limits outlined in the section on VLA Frequency Bands and Tunability. These baseband pairs are referred to as A0/C0 and B0/D0, and correspond to the IF pairs in the pre-expansion VLA.
- Up to 16 subband pairs (spectral windows) in each baseband pair
- Tuning, bandwidth, number of polarization products, and number of channels can be selected independently for each subband
- All subbands must share the same integration time
- No part of a subband can cross a 128 MHz boundary
- Subband bandwidths can be 128, 64, 32, ..., 0.03125 MHz (128 / 2n, n=0, 1, ..., 12)
- The sum over subbands of channels times polarization products is limited to 16,384.
- These may be spread flexibly over subbands and polarization products, in powers of two: 64, 128, 256, 512, ..., 16384 cross-correlation products.
- Assigning many channels to a given subband may reduce the total number of subbands available.
The remainder of this section discusses the various limitations in more detail, including some examples to show how they come up in practice.
Subband tuning restrictions
Each subband may be placed anywhere within a baseband, with the caveat that no subband may cross a 128 MHz boundary. Mathematically:
νBB0 + n*128 MHz <= νsbLow <= νsbHigh <= νBB0 + (n+1)*128 MHz |
where:
νBB0 | the lower frequency edge of the baseband; | |
n= 0, 1, ..., 7 | (i.e., any integer between 0 and 7); | |
νsbLow | the lower edge of the subband (i.e., the subband center frequency minus half the subband bandwidth); | |
νsbHigh | the upper edge of the subband (i.e., the subband center frequency plus half the subband bandwidth). |
So for example, if the baseband were tuned to cover 10000-11024 MHz, one could place a 64 MHz subband to cover 10570-10634 MHz, but not to cover 10600-10664 MHz (because that would cross the 128 MHz boundary at 10640 MHz).
The figure below illustrates these restrictions:
The black curve shows the analog filter response for an 8-bit baseband covering 1024 MHz, starting at νBB0. The dashed blue vertical lines show the 128 MHz boundaries; no subband can cross those boundaries. 128 MHz subbands are thus constrained to cover a region between two of those boundaries, and no finer tuning is possible. Narrower subbands, like the 64 MHz subband shown here in red, can be shifted around arbitrarily within one of the 128 MHz "slots", but cannot cross any of these boundaries. (The dotted vertical red lines show the boundaries of the 64 MHz subband, while the solid curve shows an illustrative line within the subband.) Observers may explore the impact of these 128 MHz boundaries on their particular observations by using the TUNE tool.
The analog filter shape defining the baseband rolls off severely at one edge of the baseband, so the 128 MHz "slot" at that edge has reduced sensitivity. The baseband edge is at the lowest sky frequency in the baseband when using upper sideband, and at the highest sky frequency in the baseband when using lower sideband.
Subband bandwidths & the digital filter response
The bandwidth for each subband may be selected independently, and can be any of 128/2n MHz, for n= 0, 1, ..., 12: 128, 64, 32, 16, 8, 4, 2, or 1 MHz, or 500, 250, 125, 62.5, or 31.25 kHz.
The usable portion of the subband is set by three effects. First, as discussed above, the analog filters which define the baseband are not perfect, leading to lower sensitivity in the 128 MHz near the baseband edge for the 8-bit samplers
Second, because the digital filters are not infinitely sharp, the rejected sideband leaks in at both edges of the subband. This leads to additional (aliased) noise, with a factor ~2 increase in the noise at the subband edges, dropping to a few per cent within a few per cent of the subband edge. The precise filter shape and noise increase is a complex but predictable function of the subband bandwidth and the subband tuning.
The third effect stems from the offset frequencies used for sideband rejection in the WIDAR correlator. The local oscillators at the individual antennas are tuned to slightly different frequencies, with those offsets taken out in the correlator. This means that each antenna observes a slightly different sky frequency, and thus some baselines will not give an interesting correlation near one edge of the subband. The maximum frequency shift is currently set to 32*f0, with the fundamental f0 being set to f0= max(25.6 kHz*sbBW/128 MHz, 100 Hz). Here sbBW is the smallest subband bandwidth within the baseband. For the wider subband bandwidths the maximum frequency shift corresponds to <1% of that bandwidth, but for narrower subbands the effect can be severe. For instance, a 31.25 kHz subband has f0= 100 Hz, and a maximum frequency shift of 3.2 kHz -- 10% of the subband may be lost on some baselines.
Spectral channels and polarization products
Each subband can have a different number of channels and polarization products, subject to two limitations:
- For the ith subband, the number of spectral channels can be:
- 64 nBlBP,i with full polarization products (RR,RL,LR,LL)
- 128 nBlBP,i with dual polarization products (RR,LL)
- 256 nBlBP,i with a single polarization product (RR or LL)
- The sum over all subbands of nBlBP,i must be less than or equal to 64. Equivalently, the sum over all subbands of spectral channels times polarization products is limited to 64x256= 16,384.
Baseline Boards are the boards in the WIDAR correlator where the actual cross-multiplications are done. There are 128 Baseline Boards arranged as 64 Baseline Board pairs (BlBPs). The limitations given here correspond to the capabilities of the individual boards, and the finite number of boards the correlator has.
Limitation #1 corresponds to the following table of the options for subband bandwidth and spectral resolution when using nBlBP Baseline Board pairs for a subband:
Subband Bandwidth and Spectral Resolution Options | |||||||
---|---|---|---|---|---|---|---|
Subband bandwidth & total velocity coverage |
Full polarization products (RR, RL, LR, LL) 64nBlBP spectral channels Channel spacing: |
Dual polarization products (RR, LL) 128nBlBP spectral channels Channel spacing: |
Single polarization product (RR or LL) 256nBlBP spectral channels Channel spacing: |
||||
128 MHz | 38400/νGHz km/s | 2000/nBlBP kHz | 600/nBlBP/νGHz km/s | 1000/nBlBP kHz | 300/nBlBP/νGHz km/s | 500/nBlBP kHz | 150/nBlBP/νGHz km/s |
64 | 19200 | 1000 / nBlBP | 300 / nBlBP | 500 / nBlBP | 150 / nBlBP | 250 / nBlBP | 75 / nBlBP |
32 | 9600 | 500 / nBlBP | 150 / nBlBP | 250 / nBlBP | 75 / nBlBP | 125 / nBlBP | 37.5 / nBlBP |
16 | 4800 | 250 / nBlBP | 75 / nBlBP | 125 / nBlBP | 37.5 / nBlBP | 62.5 / nBlBP | 18.75 / nBlBP |
8 | 2400 | 125 / nBlBP | 37.5 / nBlBP | 62.5 / nBlBP | 18.75 / nBlBP | 31.25 / nBlBP | 9.375 / nBlBP |
4 | 1200 | 62.5 / nBlBP | 18.75 / nBlBP | 31.25 / nBlBP | 9.375 / nBlBP | 15.625/nBlBP | 4.687 /n BlBP |
2 | 600 | 31.25 / nBlBP | 9.375 / nBlBP | 15.625/nBlBP | 4.687 / nBlBP | 7.8125 / nBlBP | 2.344 / nBlBP |
1 | 300 | 15.625/nBlBP | 4.687 / nBlBP | 7.8125 / nBlBP | 2.344 / nBlBP | 3.906 / nBlBP | 1.172 / nBlBP |
0.5 | 150 | 7.8125 / nBlBP | 2.344 / nBlBP | 3.906 / nBlBP | 1.172 / nBlBP | 1.953 / nBlBP | 0.586 / nBlBP |
0.25 | 75 | 3.906 / nBlBP | 1.172 / nBlBP | 1.953 / nBlBP | 0.586 / nBlBP | 0.977 / nBlBP | 0.293 / nBlBP |
0.125 | 37.5 | 1.953 / nBlBP | 0.586 / nBlBP | 0.977 / nBlBP | 0.293 / nBlBP | 0.488 / nBlBP | 0.146 / nBlBP |
0.0625 |
18.75 | 0.977 / nBlBP | 0.293 / nBlBP | 0.488 / nBlBP | 0.146 / nBlBP | 0.244 / nBlBP | 0.073 / nBlBP |
0.03.25 | 9.375 | 0.488 / nBlBP | 0.146 / nBlBP | 0.244 / nBlBP | 0.073 / nBlBP | 0.122 / nBlBP | 0.037 / nBlBP |
Subband bandwidth and spectral resolution options. Note that the table entries refer to the spacing between spectral channels -- that spacing is before any frequency smoothing, so these channels are not independent.
|
Here are four examples of allowed general observing setups which use all 64 BlBPs to produce the maximum number of channels times polarization products:
Baseband | Subband |
Pol'n Products |
Spectral channels |
nBlBP | |
Example 1 | A0/C0 | sb0 | RR | 16384 | 64 |
Example 2 | A0/C0 | sb0 | RR | 8192 | 32 |
A0/C0 | sb1 | RR, LL | 1024 | 8 | |
A0/C0 | sb2 | RR, LL | 512 | 4 | |
B0/D0 | sb0 | RR, LL | 2048 | 16 | |
B0/D0 | sb1 | RR,RL,LR,LL | 256 | 4 | |
Example 3 | A0/C0 | sb0 | RR | 8192 | 32 |
A0/C0 | sb1 | LL | 1024 | 4 | |
A0/C0 | sb2 | RR, LL | 1024 | 8 | |
A0/C0 | sb3 | RR,RL,LR,LL | 1024 | 16 | |
A0/C0 | sb4 | RR,RL,LR,LL | 256 | 4 | |
Example 4 | A0/C0 | sb0-7 | RR,RL,LR,LL | 64 | 8 x 1 |
A0/C0 | sb8 | RR, LL | 4096 | 1 x 32 | |
B0/D0 | sb0-7 | RR,RL,LR,LL | 64 | 8 x 1 | |
B0/D0 | sb8 | RR, LL | 2048 | 1 x 16 |
Number of subbands and Baseline Board stacking
The number of subbands available in a baseband depends on the correlator resources available. This means that using Baseline Board stacking to increase the number of channels in one subband, may limit the number of subbands available in one or more of the basebands. Understanding how this works requires understanding some of the details of the correlator hardware. That understanding is built into the GOST, and observers may simply use that tool to find out whether their particular setup will, in fact, work. But the results can be confusing without some understanding of the hardware constraints from which they arise. These hardware constraints are complex, and most observers will not need to understand these details. This section is for those few who are attempting complex line experiments, and who find the GOST or SRCT restricting the number of subbands and/or channels they can use in unexpected ways.
First let us consider how the correlator hardware is organized. The cross-multiplications in the WIDAR correlator are spread across 64 Baseline Board (BlB) pairs, arranged into 4 quadrants of 16 BlB pairs each. Each baseband is connected directly to one of those quadrants. In the simplest mode, each of the 16 BlB pairs of a quadrant handles the correlations for one of the 16 subbands of the corresponding baseband. Four basebands and four quadrants are required, in order to handle the full 8 GHz of bandwidth per polarization provided by the 3-bit (wideband) samplers: that 8 GHz is split into four basebands of 2 GHz each, with each baseband fed into a different BlB quadrant. Each BlB pair in that quadrant handles a subband of maximum bandwidth 128 MHz, so 16 BlB pairs handles 16 subbands for a total of 16x128 MHz= 2048 MHz.
A single BlB pair produces 256 cross-correlations per baseline for a single subband, which can be used for a single polarization product (e.g., RR with 256 spectral channels), or two (RR+LL with 128 spectral channels each), or four (RR,RL,LR,LL, with 64 spectral channels each).
Baseline Board (BlB) pair usage for a simple 4 x 2 GHz (3-bit sampler) experiment. The different colors correspond to different baseband pairs. The BlB quadrants are indicated along the left-hand side, while the numbers across the top represent the BlB pairs within those quadrants. Numerals within the colored boxes label the subbands within the basebands. In this setup each subband is correlated using a single BlB pair, using all of the 64 available BlB pairs. |
When using the 8-bit samplers, the total bandwidth is only 2 GHz per polarization, split into two basebands of 1 GHz each. The simplest continuum setup uses only two quadrants, since there are only two basebands; and only 8 subbands are required to span the 8x128 MHz= 1024 MHz of each baseband. Three-quarters of the correlator BlB hardware remain unused.
Baseline Board (BlB) pair usage for a simple 2 x 1 GHz (8-bit sampler) experiment. The cyan boxes correspond to the first baseband pair (IF A0/C0), fed directly into Q1; the yellow boxes to the second (IF B0/D0), fed directly into Q3. The BlB quadrants are indicated along the left-hand side, while the numbers across the top represent the BlB pairs within those quadrants. Numerals within the colored boxes label the subbands within the basebands. In this setup each subband is correlated using a single BlB pair, using only 16 of the 64 available BlB pairs. |
The spectral line mode allows access to these `extra' correlator resources through Baseline Board stacking: using multiple BlB pairs to process the same subband and produce more cross-correlations for that subband. This is done using crossbar switches which make the data for a single subband available to several BlB pairs. Those BlB pairs can then be used to produce more spectral channels for that subband, with n BlB pairs producing 256*n cross-correlations per baseline. The limit on the total number of cross-correlations (16384) stems from the total number of BlB pairs (64): 64x256= 16384.
Unfortunately completely flexible crossbar switches are expensive, and could not be implemented in the VLA's new correlator. This means that one cannot route a given subband to a randomly-chosen BlB pair. The routings which are possible, are as follows:
- A subband in a baseband can be routed to any BlB pair within the corresponding quadrant.
- Data coming into a given BlB pair in one quadrant, can be routed to the corresponding BlB pair in any other quadrant.
Routing option #1 means that one could use all the BlB pairs within a quadrant to correlate a single subband, yielding 16x256= 4096 cross-correlations for that subband:
Baseline Board (BlB) pair usage when assigning all BlB pairs within a quadrant to a single subband within the corresponding baseband. The cyan boxes correspond to the first baseband pair (IF A0/C0), fed directly into Q1; the yellow boxes to the second (IF B0/D0), fed directly into Q3. The BlB quadrants are indicated along the left-hand side, while the numbers across the top represent the BlB pairs within those quadrants. Numerals within the colored boxes label the subbands within the basebands. In this setup each subband is correlated using 16 BlB pairs, yielding 16x256= 4096 cross-correlations for those subbands, and using 32 of the 64 available BlB pairs. |
Routing option #2 means that one could use the BlB pairs in all 4 quadrants to correlate a single subband. One simple case would use 4 BlB pairs to correlate each of the 16 subbands in a single baseband, yielding 4x256= 1024 cross-correlations for each of those subbands. Note that in this case, no BlB pairs are left to correlate any data from the second baseband.
Baseline Board (BlB) pair usage when assigning the BlB pairs in all 4 quadrants to correlate the corresponding subbands in a single baseband. The BlB quadrants are indicated along the left-hand side, while the numbers across the top represent the BlB pairs within those quadrants. Numerals within the colored boxes label the subbands within the baseband. In this setup each subband is correlated using 4 BlB pairs, yielding 4x256= 1024 cross-correlations for those subbands, and using all of the 64 available BlB pairs. |
Using routing option #2 does come with a subtle cost: assigning a BlB pair in quadrant X to correlate a subband corresponding to quadrant Y, removes that BlB pair from use in the baseband corresponding to quadrant X...and hence also removes the corresponding subband in that baseband. So getting more channels for a subband in one baseband, may prevent the use of a subband in a different baseband. To take a simple example, consider an experiment where one wishes to observe a single line in dual polarization with 512 channels (requiring 4 BlB pairs), plus as much continuum bandwidth as possible. Naively one would say there are 16 subbands in each baseband; one is used for the spectral line, so that leaves 16+15= 31 subbands, and with the widest subband bandwidth (128 MHz) the total available continuum should be 31x128 MHz= 3968 MHz per polarization. Actually however there are only 15+15 subbands available, or 30x128 MHz= 3840 MHz per polarization, because the spectral line subband has "eaten" one BlB pair corresponding to the other baseband:
Baseline Board (BlB) pair usage for a continuum plus single line experiment. The cyan boxes correspond to the first baseband pair (IF A0/C0), fed directly into Q1; the yellow boxes to the second (IF B0/D0), fed directly into Q3. The BlB quadrants are indicated along the left-hand side, while the numbers across the top represent the BlB pairs within those quadrants. Numerals within the colored boxes label the subbands within the basebands. In this setup subband #15 of A0/C0 is correlated using 4 BlB pairs, yielding 4x256= 1024 cross-correlations for that subband, while the other 15 subbands in each baseband are correlated using a single BlB pair each. The use of 4 BlB pairs for a subband in baseband A0/C0 has eliminated the corresponding subband of the other baseband (B0/D0). Note that only 30+4= 34 of the 64 available BlB pairs have been used -- there is plenty of correlator hardware available, but the signal for the spectral line subband cannot be routed to them. |
If the same spectral line required twice as many channels, this result in the loss of two subbands in both of the basebands:
Baseline Board (BlB) pair usage for a continuum plus single line experiment. The cyan boxes correspond to the first baseband pair (IF A0/C0), fed directly into Q1; the yellow boxes to the second (IF B0/D0), fed directly into Q3. The BlB quadrants are indicated along the left-hand side, while the numbers across the top represent the BlB pairs within those quadrants. Numerals within the colored boxes label the subbands within the basebands. In this setup subband #14 of A0/C0 is correlated using 8 BlB pairs, yielding 8x256= 2048 cross-correlations for that subband, while the other 14 subbands in each baseband are correlated using a single BlB pair each. The use of 8 BlB pairs for a subband in baseband A0/C0 has eliminated one additional subband in baseband A0/C0, and two subbands in B0/D0. Note that only 28+8= 36 of the 64 available BlB pairs have been used -- there is no way to route the signals from the "missing" subbands to that hardware. |
In some cases one may want to use a different routing, to use up subbands in one baseband in preference to another. For instance, the same spectral line setup (2048 cross-correlations for a single spectral line subband, plus as much continuum as possible) could be set up to allow 13 "continuum" subbands in the A0/C0 baseband, and the full 16 "continuum" subbands in B0/D0:
Baseline Board (BlB) pair usage for a continuum plus single line experiment. The cyan boxes correspond to the first baseband pair (IF A0/C0), fed directly into Q1; the yellow boxes to the second (IF B0/D0), fed directly into Q3. The BlB quadrants are indicated along the left-hand side, while the numbers across the top represent the BlB pairs within those quadrants. Numerals within the colored boxes label the subbands within the basebands. In this setup subband #13 of A0/C0 is correlated using 8 BlB pairs, yielding 8x256= 2048 cross-correlations for that subband, while the other 13 subbands in A0/C0 and 16 subbands in B0/D0 are correlated using a single BlB pair each. In this setup the use of 8 BlB pairs for a subband in baseband A0/C0 has eliminated two additional subbands in baseband A0/C0, while retaining all 16 subbands in B0/D0. Note that only 29+8= 37 of the 64 available BlB pairs have been used. [Shared risk observations can assign a number of BlB pairs which is not a power of two, and hence could assign one additional BlBP (Q4-BlBP 15) to the spectral line subband.] |
Understanding these confusing constraints can help observers set up the VLA more effectively to achieve their scientific goals. For instance, for a mixed line+continuum experiment it works best to use the resource tools to set up the baseband tunings and subband channelization for the most important lines first, then add the desired continuum, and then see what correlator resources remain for any lines of secondary interest. The Shared Risk Observing (SRO) program provides even more options and hence benefits even more from a clear understanding of the correlator hardware.
Once again, both the GOST and the Shared Risk Resource Catalog Tool (SRCT) implement all of these constraints, and are generally smart enough to figure out the routing scheme that works best for your particular request.
15.2.4. Continuum Subarrays
The continuum subarray option offers two 1 GHz baseband pairs (two 8-bit baseband pairs, each with 1 GHz in RCP and 1 GHz in LCP) in up to 3 subarrays, with the same spectral channel and polarization product options as are available for wideband observing. The setup for each subarray is completely independent, in terms of observing frequency, polarization products, and integration times.
When using three subarrays, there are some restrictions on the number of antennas in each subarray. The Baseline Board in the correlator treats each set of 4 antennas independently, using a separate column of correlator chips. With 8 such columns, the correlator can handle up to 8x4= 32 antennas. The correlator configuration software requires that a given column not be split across subarrays. This does not matter when using only two subarrays, but forces some subtle restrictions when using three. For instance, one cannot observe with 9 antennas in each of 3 subarrays, because 9 antennas requires three columns (two with 4 antennas each, and one with 1 antenna); three subarrays of 9 antennas each would require 3x3= 9 columns, one more than are actually available. Splitting the array into 10, 9, and 8 antennas is allowed, since the first two subarrays use 3 columns each, while the third uses only two.
The following table gives four examples of how correlator resources can be split into multiple subarrays. Antennas in each subarray are color-coded: red for subarray 1, green for subarray 2 (if present), and blue for subarray 3 (if present). The last column gives the number of antennas in each subarray (e.g., in the setup shown in the first row, subarray 1 has 10 antennas; subarray 2 has 9 antennas; and subarray 3 has 8 antennas). In all cases a total of 27 antennas are used. [The columns are numbered in reverse order (C7 to C0) to match the numbering scheme used on the actual Baseline Boards.]
Some Possible Subarray Options | ||||||||
---|---|---|---|---|---|---|---|---|
Number of antennas correlated using each Baseline Board column | Number | |||||||
C7 | C6 | C5 | C4 | C3 | C2 | C1 | C0 | of antennas |
4 | 4 | 2 | 4 | 4 | 1 | 4 | 4 | 10 + 9 + 8 |
4 | 4 | 4 | 2 | 4 | 4 | 4 | 1 | 14 + 13 |
4 | 4 | 3 | 4 | 4 | 3 | 3 | 2 | 11 + 11 + 5 |
4 | 4 | 4 | 4 | 4 | 4 | 3 | -- | 27 |
15.3. Shared Risk Observing Program
The Shared Risk Observing (SRO) program allows observers access to capabilities that can be set up via the Observation Preparation Tool and run through the dynamic scheduler (without intervention), but are not well tested. In semester 2013A this program provides a number of new correlator capabilities:
- Data rates as high as 60 MB/s and integration times as short as 50msec.
- Use of the 3-bit samplers in any frequency band, with subbands which can be narrower than 128 MHz. All subbands within a baseband must be contiguous, and must have the same bandwidth. Note that it makes little sense to use the 3-bit samplers for L and S bands, where the 8-bit system cover the entire available observing band with higher dynamic range and better sensitivity.
- Mixing 3- and 8-bit basebands (3-bit AC and 8-bit BD, or vice versa) in the same correlator configuration.
- A simple spectral sweep mode using recirculation to provide 512 channels in dual polarization (RR, LL) for each of 32 64-MHz subbands spanning two 1 GHz basebands.
- Recirculation is a technique whereby narrowband data are fed back into a cross-multiplier chip multiple times to obtain more lags, and hence more spectral channels.
- This particular recirculation configuration provides 2 x 32 x 512= 32768 correlations per baseline.
Those proposing for Shared Risk Observating must use the Shared Risk version of the Resource Catalog Tool (SRCT) to define their correlator set-up(s) instead of the GOST. Note however that the SRCT does not fully enforce the SRO restrictions outlined here, particularly for 3-bit observations, so some configurations allowed by the SRCT may correspond to Resident Shared Risk observations. Those in doubt are encouraged to get in touch with NRAO staff well before the proposal deadline.
Shared risk observations: spectral line configurations
Arbitrary Baseline Board stacking using the 8-bit samplers is allowed under the Shared Risk Observing (SRO) program for semester 2013A. The basic capabilities are given in the table above (Subband Bandwidth and Spectral Resolution Options), the difference being that the number of Baseline Board pairs assigned to a subband (nBlBP) can be any integer up to 64, rather than only a power of two. Similarly the interaction between BlBP stacking the number of independent subbands is as described above, but this interaction is more subtle and potentially confusing in the shared risk case.
To give some idea of the possibilities, the following tables give two examples of possible shared risk configurations. The SRCT display shows how the Baseline Boards are used to process the individual subbands. The cyan boxes show the Baseline Boards used to process data from baseband A0/C0, while the yellow boxes show Baseline Boards used to process data from baseband B0/D0.
Shared Risk Configuration Example #1 |
||||
---|---|---|---|---|
Baseband | Subband | Pol'n products | Spectral channels | nBlBP |
A0/C0 | sb0 | RR | 10240 | 40 |
A0/C0 | sb1 | LL | 768 | 3 |
A0/C0 | sb2 | RR,LL | 2176 | 17 |
B0/D0 | sb0 | RR | 256 | 1 |
B0/D0 | sb1 | RR,LL | 384 | 3 |
SRCT display: |
Shared Risk Configuration Example #2 |
||||
---|---|---|---|---|
Baseband | Subband | Pol'n products | Spectral channels | nBlBP |
A0/C0 | sb0 | RR | 4352 | 17 |
A0/C0 | sb1 | RR, LL | 1152 | 9 |
B0/D0 | sb0 | RR,RL,LR,LL | 192 | 3 |
B0/D0 | sb1 | RR, LL |
4480 | 35 |
SRCT display: |
Note that the individual subbands can have different bandwidths, and those bandwidths may be chosen completely independently of the number of spectral channels in each subband. So for instance a subband with bandwidth 2 MHz and 1152 spectral channels would have a channel separation of 2 MHz/1152= 1.736 kHz; but the observer could equally well choose a bandwidth of 64 MHz for that subband, leading to a channel separation of 64 MHz/1152= 55.56 kHz.
Again, while the SRCT understands the various hardware constraints on correlator configuration, it does not distinguish between configurations which are available as Shared Risk Observing (SRO), and those which are only available through the Resident Shared Risk Observing (RSRO) program. If in doubt, please consult with NRAO staff.
15.4. Resident Shared Risk Observing Program
Even more sophisticated correlator configurations are available under the Resident Shared Risk Observing (RSRO) program. Examples include:
- More extensive use of recirculation to provide many times more spectral channels, flexibly distributed among the subbands.
- Access to more than 16 subbands per 8-bit baseband.
- Data rates beyond 60 MB/s and integration times below 50 msec.
- Pulsar phase bins and real-time de-dispersion.
- Multiple phase/delay centers.
- Access to wideband phased array output for high time resolution observations.
- Real-time RFI excision and transient detection.
- 7-bit correlations.
- More subarrays, and more complex correlator configurations for each subarray.
- Multiple integration times within a single subarray (e.g., slower integrations for spectral line than for continuum subbands).
Those interested in these capabilities should consult with NRAO staff, and be prepared to make a significant commitment to aid in commissioning the desired mode(s).
Correlator configurations and observing details for RSRO programs must be fully explained and justified in the Scientific and Technical Justification section of the proposal.
16. VLBI Observations
The VLA can now participate in VLBI observations with the VLBA. Currently, this is only allowed in phased array mode (i.e. single dish is not allowed) with highly restricted WIDAR correlator configurations. For more details see the VLBI at the VLA documentation. In phased array mode the program TelCal derives the antenna-based delay and/or phase corrections needed for antenna phasing in real time. This delay correction is then applied to the antenna signals in the correlator and the visibilities are then summed. The resulting phased array sum is recorded in VDIF format on the Mark5C recorder at the VLA site. The disk(s) are then transported to Socorro, NM and correlated on the DiFX correlator with other VLBI stations which participated in the observation. The standard WIDAR data are also archived in the NRAO science data archive.
17. Snapshots
The two-dimensional geometry of the VLA allows a snapshot mode whereby short observations can be used to image relatively bright unconfused sources. This mode is ideal for survey work where the sensitivity requirements are modest.
Single snapshots with good phase stability of strong sources should give dynamic ranges of a few hundred. Note that because the snapshot synthesized beam contains high sidelobes, the effects of background confusing sources are much worse than for full syntheses, especially at 20 cm in the D configuration, for which a single snapshot will give a limiting noise of about 0.2 mJy. This level can be reduced by taking multiple snapshots separated by at least one hour. Use of the AIPS program 'IMAGR' or CASA task 'clean' is necessary to remove the effects of background sources. Before considering snapshot observations at 20 cm, users should first determine if the goals desired can be achieved with the existing "Faint Images of the Radio Sky at Twenty-centimeters survey(FIRST, http://www.cv.nrao.edu/first/)" (B configuration) or the "Co-Ordinated Radio 'N' Infrared Survey for High-mass star formation (CORNISH, http://www.ast.leeds.ac.uk/Cornish/public/index.php)" (B configuration), or the NRAO VLA Sky Survey (NVSS, http://www.cv.nrao.edu/nvss/) (D configuration, all-sky) surveys.
18. Shadowing and Cross-Talk
Observations at low elevation in the C and D configurations will commonly be affected by shadowing. It is strongly recommended that all data from a shadowed antenna be discarded. This will automatically be done during filling (CASA task importasdm) when using the default inputs (Note: For archival VLA data, the AIPS task FILLM can also flag based on shadowing). AIPS task UVFLG can be used to flag modern VLA data based on shadowing as well, although it will only flag based on antennas in the dataset, and is ignorant of antennas in other sub-arrays.
Cross-talk is an effect in which signals from one antenna are picked up by an adjacent antenna, causing an erroneous correlation. This effect is important at low frequencies in compact configurations. Careful examination of the visibilities is necessary to identify and remove this form of interference.
19. Combining Configurations and Mosaicing
Any single EVLA configuration will allow accurate imaging up to a scale approximately 30 times the synthesized beam. Objects larger than this will require multiple configuration observations. It is advisable that the frequencies used be the same for all configurations to be combined. Objects larger than the primary antenna pattern may be mapped through the technique of interferometric mosaicing. Time-variable structures (such as the nuclei of radio galaxies and quasars) cause special, but manageable, problems. See the article by Mark Holdaway in reference 2 (Documentation) for more information.
Guidelines for mosaicing with the EVLA are given at http://evlaguides.nrao.edu/index.php?title=Mosaic_Observing
20. Pulsar Observing
While the WIDAR correlator will eventually allow flexible pulsar phase-binning and real-time correction for dispersion, those aspects of the correlator have not yet been commissioned. Similarly the use of the phased array for high time-resolution observations is not yet possible for general use.