Performance of the VLA during the Next Semester
Resolution
Resolution
The VLA's resolution is generally diffraction-limited, and thus is set by the array configuration and the observing frequency. Like all synthesis arrays, the VLA is sensitive only to structures on a range of angular scales between the diffraction limit (the smallest angular scale detectable) and a "Largest Angular Scale" (which depends on the fringe spacing formed by the shortest baselines in the configuration). For emission structures smaller than the diffraction limit (θ ∼ λ/Bmax), the VLA acts like a single-dish instrument—the resulting image is smoothed to the resolution of the array. For emission structures larger than the detectable range, the VLA is simply blind to the emission; this is a limitation unique to interferometers. No subsequent processing can fully recover the missing information from these large scales. It can only be obtained by observing in a more compact VLA array configuration or with data from an instrument that is sensitive to the missing angular scales, such as a large single dish or a compact array of smaller antennas.
Table 3.1.1 displays the VLA's resolution and the scale at which severe attenuation of large-scale structure occurs. This table shows the maximum and minimum antenna separations, the approximate synthesized beam size (full width at half-power; the resolution element) for the central frequency for each band, and the largest angular scale of detectable emission.
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) | 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 |
Band | Largest Angular Scale θLAS(arcsec)1,4 | |||
74 MHz (4) | 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.
- Notes:
- 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 angular scale structure visible to the array.
- 2. The listed resolutions are appropriate for sources with declinations between −15 and +75 degrees.
- 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 may slightly 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.
The following figure is a graphical representation of the synthesized beamwidths for natural and robust weighting for the four main array configurations between 1 and 50 GHz. Also available are synthesized beamwidth figures for the low frequency (1–12 GHz) and the high frequency (12–50 GHz) receiver bands.
Sensitivity
Sensitivity
The theoretical thermal noise expected for an image using natural weighting of the visibility data is given by:
\[\Delta I_m = \frac{SEFD}{\eta_{\rm c}\sqrt{n_{\rm pol}N(N-1)t_{\rm int}\Delta\nu}}\] | (eq. 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 (~0.93 with the use of the 8-bit samplers).
- - 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 3.2.1 shows the SEFDs as a function of frequency used in the VLA exposure calculator for those Cassegrain bands currently installed on VLA antennas, and include the contribution to Tsys from atmospheric emission at the zenith. Figure 3.2.2 shows the SEFDs as a function of frequency for the P-band; these measurements are based on imaging of a field far from the galactic plane. Table 3.2.1 gives the SEFD at some fiducial VLA frequencies.
Figure 3.2.1: SEFD used in the Exposure Calculator for the VLA. Left: The system equivalent flux density as a function of frequency for the L, S, C and X-band receivers. Right: The system equivalent flux density as a function of frequency for the Ku, K, Ka, and Q-band receivers. SEFDs at Ku, K, Ka, and Q bands include contributions from Earth's atmosphere and were determined under good conditions.
Figure 3.2.2: The SEFD used in the VLA Exposure Calculator as a function of frequency for the P-band receiver
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 3.2.1 shows the confusion noise in D configuration assuming robust weighting (see Condon et al. 2012, ApJ, 758), 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 3.2.1; (ii) confusion from the sidelobes of uncleaned sources lying outside the image, often from sources in the sidelobes of the primary beam. This confusion 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–15% loss in sensitivity above that expected—i.e., the efficiency factor ηc = 0.78 to 0.83.
Frequency | SEFD (Jy) |
RMS confusion level |
||
---|---|---|---|---|
0.39 GHz (P) | 2790 | 5330 | ||
1.5 GHz (L) | 420 | 74 | ||
3.0 GHz (S) | 370 | 12 | ||
6.0 GHz (C) | 310 | 2 | ||
10.0 GHz (X) | 250 | negligible | ||
15 GHz (Ku) | 320 | negligible | ||
20 GHz (K) | 500 | negligible | ||
33 GHz (Ka) | 600 | negligible | ||
45 GHz (Q) | 1300 | negligible |
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. Observing 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 (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—which are currently unavailable but will be implemented at some time in the future—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 complex gain calibrator (amplitude and phase) and target source are observed.
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 items:
- 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.
Both CASA and AIPS allow the application of elevation-dependent gains and an estimated opacity generated from ground-based weather through the CASA tasks gencal and plotweather, and AIPS task INDXR.
- 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 11 mph (5 m/s). Wind speeds near the stow limit 45 mph (20 m/s) will have a similar negative effect at 8 and 15 GHz.
To achieve increased pointing accuracy, referenced pointing is recommended where a nearby calibrator is observed in interferometric pointing mode every hour or so. The local pointing corrections 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–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 online guide to the exposure calculator; the exposure calculator provides a graphical user interface to these equations.
Special caveats apply for P-band (230–470 MHz) observing. The SEFD's in Figure 3.2.2 or that listed in table 3.1.2 are from an observation taken far from the Galactic plane, where the sky brightness is about 30K. At P-band, Galactic synchrotron emission is very bright in directions near the Galactic plane. The system temperature increase due to Galactic emission will degrade sensitivity by factors of two to three for observations in the plane, and by a factor of five or more at or near the Galactic center. Additionally, the antenna efficiency (currently about 0.31 for 300 MHz) will decline with both increasing and decreasing frequencies from the center of P-band.
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:
\[T_{\rm b} = \frac{S \lambda^2}{2k\Omega} = F \cdot S\] | (eq. 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 parameter F depends on the synthesized beam, therefore on the array configuration, and has the approximate value of 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, then 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 6 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):
\[M_{\rm HI} = 2.36 \times 10^5 D^2 \sum S \Delta V ~~~~M_\odot\] | (eq. 3) |
where D is the distance to the galaxy in Mpc, and SΔV is the HI line area in units of Jy km/s.
VLA Frequency Bands and Tunability
-
- 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 the NRAO Helpdesk for further information.
- 2. The 4-band system is currently under development. The default frequency range maximizes sensitivity of the system and provides a nominal bandwidth of 12 MHz and a channel resolution of 64 kHz.
- 3. The default setup for P-band will provide 16 subbands from the A0/C0 IF pair, each 16 MHz wide, to cover the frequency range 224–480 MHz. The channel resolution is 125 kHz.
- 4. The default frequency setup 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.
Tuning Restrictions
In general, for all frequency bands except Ka, if the total span of the two independent IF pairs of the 8-bit system (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 A0/C0 IF pair must be greater than 32.0 GHz. There is no restriction on the B0/D0 frequency, unless the B0/D0 band overlaps the A0/C0 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 B0/D0 frequency must be lower than the A0/C0 frequency.
For the 3-bit system, the maximum frequency span permitted for the A1/C1 and A2/C2 IF pairs is about 5000 MHz. The same restriction applies to B1/D1 and B2/D2. The tuning restrictions given above for the separation and location of the 8-bit pairs A0/C0 and B0/D0 also apply to the 3-bit pairs, with A0/C0 replaced by A1/C1 and A2/C2, and B0/D0 replaced by B1/D1 and B2/D2.
VLA Samplers
The VLA is equipped with two different types of samplers, 8-bit with 1GHz bandwidth, and 3-bit with 2GHz bandwidth. The choice depends on your science goals and on technicalities described below.
The 8-bit 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. Some restrictions apply, depending on band, as described in the section on Frequency Bands and Tunability.
The 3-bit 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. The tuning restrictions are described in the section on Frequency Bands and Tunability. Taken together, the eight 3-bit samplers offer a maximum of 8192 MHz coverage with full polarization.
Which set to use?
- S, L, and 4/P-band observations, whether line or continuum, should use the 8-bit sampler set.
- C and X-band continuum observations should use 3-bit samplers in order to exploit the full 4 GHz bandwidth: in spite of the 15% reduction in sensitivity that comes with 3-bit (at equal bandwidth to the 8-bit samplers—see below for details) and the reduced effective bandwidth after removing RFI, this still provides superior overall sensitivity. For more details we refer to EVLA memo 166.
- Note: C-band is impacted by strong RFI caused by microwave links near 6 GHz in the A and B configurations. As a result, 3-bit data obtained with the standard setup are corrupted. We advise observers to use mixed 3-bit and 8-bit samplers. For more details, refer to the VLA Observing Guide.
- Ku, K, Ka, and Q-band continuum observations should use the 3-bit samplers for maximum bandwidth.
- Wide-band spectral line searches requiring more than 2 GHz span should use the 3-bit samplers.
- Spectral-line observations which fit within two, possibly disjoint, 1 GHz bands should use the 8-bit set.
- Simultaneous continuum and high resolution spectral line observation can use mixed 3-bit and 8-bit samplers. The 3-bit samplers in this case will be set up to deliver the continuum data, while the 8-bit samplers will be for the spectral line data. This mix mode can be used in C-band and higher.
Major Characteristics of each Set
The 8-bit samplers are warranted for observations at 4/P, L, and S-bands. The full analog bandwidth from the receivers fits within the 2048 MHz span covered by the samplers.
For the 3-bit samplers, users need to be aware of the following issues:
- Sensitivity: compared to the 8-bit system, the sensitivity of the 3-bit samplers is worse by ~15% (at equal bandwidth). Alternatively, a given continuum noise level requiring on-source integration time T with the 8-bit (two bands of 1GHz), requires 0.33T with the 3-bit (4 bands of 2GHz, assuming the bandwidth is available from the front end).
- Resonances: each of the eight 3-bit samplers on an antenna has a resonance about 3 MHz wide. Each resonance is independent of all others, so there is no correlated signal between antennas. The resonance degrades the spectrum in its narrow frequency range, but has little effect on continuum observing. Bandpass solutions will be affected, but can be interpolated over. Spectral-line calibration and images at the affected frequencies will show significant loss in sensitivity. The resonances are easily seen in autocorrelation spectra, and it is recommended that users, especially spectral-line users, utilize these to locate the compromised frequencies.
- Amplitude Calibration: The traditional method for both 8- and 3-bit systems is to observe a flux-density calibrator, use self-cal to determine the antenna amplitude calibration factors (gains), and transfer the gains to the phase calibrator and target. For 3-bit samplers this procedure gives results good to 5% between elevations of 20–70 degrees. (Expect worse at the upper edge of Q-band and/or during bad weather.) The switched power data can be used to correct for system gain variations and works well for the 8-bit samplers. For 3-bit samplers, the Pdif depends on the Psum, i.e., Pdif is non-linear and its application will bias the resulting visibilities by 5–10%. The origin of this effect is understood, but we have not yet determined how best to compensate for it. Because of this, we do not recommend use of the Psum and Pdif data to calibrate visibilities from the 3-bit samplers. We do, however, recommend that the requantizer gains in the switched power data be applied to remove gain changes. For more information about the switched power, Psum, and Pdif, see EVLA memo 145.
Setting up the 8-bit or 3-bit Samplers
Either set requires an initial scan for each individual LO (frequency) tuning, during which power levels are optimized.
For the 8-bit system, a dummy scan of 1 minute duration is sufficient for each tuning. This is usually done while the antennas are slewing at the start of an observing file, as the pointing direction of the antennas is not critical.
For the 3-bit system, the requirements are more demanding, see the section on 3-bit setup within the Guide to Observing with the VLA. The minimum setup time is 1 minute for each tuning to adjust the power levels and bandpass slopes across the 2GHz samplers. These values are retained and applied if the tuning is re-encountered in the same observation. Additionally, every time the LO setup is changed—whether or not it is new (e.g., changing from 8-bit X-band reference pointing back to target)—a scan of 30 seconds is needed to reset the subband gains (requantizers) in the correlator. For better amplitude calibration at high frequencies, the 3-bit initial setup should be near the elevation of the target, so do it after the first 8-bit setup described above. For 3-bit observing without 8-bit (e.g., C or X-band without reference pointing), the power variation with elevation is small, so the 3-bit setup can be done at any elevation.
For settings that use a mix of 3-bit and 8-bit samplers, the guidelines to set up the 3-bit samplers should be followed.
Other issues
The overhead for setup of 3-bit samplers can eat into observing time, especially for projects with many different LO settings, and/or sources all over the sky accompanied by band change, reference pointing, and requantizer reset for each direction. The impact is most severe for short scheduling blocks.
Polarization testing conducted so far indicates no degradation of performance by using the 3-bit samplers.
Time Resolution and Data Rates
The default integration times for the various array configurations and frequency bands are as follows:
Configurations | Observing Bands | Default integration time |
---|---|---|
A, B, C, D | 4 P | 2 seconds |
A | L S C X Ku K Ka Q | 2 seconds |
B | L S C X Ku K Ka Q | 3 seconds |
C, D | X Ku K Ka Q | 3 seconds |
C, D | L S C | 5 seconds |
Observations with the 3-bit (wideband) samplers, when applicable, should 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:
Proposal type |
Minimum integration time | Maximum data rate |
---|---|---|
General Observing (GO) | 50 msec | up to 60 MB/s (216 GB/hr) |
Shared Risk Observing (SRO) | 50 msec | > 60 MB/s (216 GB/hour) and up to 100 MB/s (360 GB/hour) |
Resident Shared Risk Observing (RSRO) | < 50 msec | > 100 MB/s (360 GB/hr) |
Note that integration times as short as 5 msec and data rates as high as 300 MB/s can be supported for some observing, though any such observing is considered Resident Shared Risk Observing. For these short integration times and high data rates there will be limits on bandwidth and/or number of antennas involved in the observation. Those desiring to utilize such short integration times and high data rates should consult with NRAO staff.
The maximum recommended integration time for any VLA observing is 10 seconds.
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) × Nant × (Nant − 1)/(27×26) / (Δt/1 sec) |
~ 160 GB/hr × (Nchpol/16384) x Nant × (Nant − 1)/(27×26) / (Δt/1 sec) | |
~ 3.7 TB/day × (Nchpol/16384) × Nant × (Nant − 1)/(27×26) / (Δt/1 sec) |
Here Nchpol is the sum over all subbands of spectral channels times polarization products:
Nchpol = Σi Nchan,i × 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 and 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 ~2 seconds for standard observations, and 0.8 seconds for shared risk observations.
We note that frequency averaging in the correlator will reduce the total number of channels. Therefore, the data rate and the data volume will be reduced by the same channel averaging factor. See the Chromatic Aberration section for more details on the frequency averaging in the correlator and to assess its impact on your science.
These data rates are challenging for transfer and analysis. Data may either be downloaded via ftp over the Internet, or shipped on hard drives for large data sets or for those with slow Internet connections (please review the data shipping policy). For users whose science permits, the Archive Access Tool allows some level of frequency averaging in order to decrease data set sizes before ftp; note that the full spectral resolution will be retained in the NRAO archive for all observations.
†Note: The data rate formula given above does not account for the auto-correlations delivered by WIDAR. Precise data rate values can be obtained through the use of the Resource Catalog Tool for proposing (RCT-proposing).
Radio Frequency Interference
The very wide bandwidths of the VLA mean that radio frequency interference (RFI) will be present in a far larger fraction of current observations than in observations made with the old systems. Considerable effort has gone into making the VLA's 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 (C, S, L, and the low-band system) bands, and is most serious in 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 protected for radio astronomy are: 73.0-74.6 MHz, 322.0-328.6 MHz, 1400–1427 MHz, 1660.6–1670.0 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 significant external interference should occur within these bands.
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. Users concerned about the precise frequencies of strong RFI, and the likelihood of being affected, are encouraged to peruse these plots. To access these plots, or for more information on RFI, including the impact of satellite transmissions, please see the RFI section in the Guide to Observing with the VLA.
Subarrays
The continuum subarray option offers two 1 GHz baseband pairs with the 8-bit samplers in up to 3 subarrays or four 2 GHz baseband pairs with the 3-bit samplers, 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 8×4 = 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 3×3 = 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.
Table 3.8.1 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.
Number of antennas correlated using each Baseline Board column | Number of antennas | |||||||
---|---|---|---|---|---|---|---|---|
C7 | C6 | C5 | C4 | C3 | C2 | C1 | C0 | |
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 |
For more information on subarrays, please see the Subarray section in the Guide to Observing with the VLA.
Positional Accuracy & Astrometry
The position of a target 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 must be known to within the antenna primary beam, except perhaps for mosaicking observations. In the special case of using the phased VLA as a VLBI element, the a priori position must be accurate to within the synthesized beam of the array.
In post-processing, target positions are typically determined from an image made after phase calibration, i.e., correcting the antenna and atmospheric phases as determined on the reference source. The accuracy of the calibration determines the accuracy of the positions in the image. Note that phase self-calibration imposes the assumed position of the model, i.e., makes the position indeterminate. Therefore, an absolute position cannot be determined after self-calibration, but relative positions between features within a self-calibrated image are valid.
It may help to think of astrometry as two methods, narrow-field and wide-field.
Narrow-field astrometry
In narrow-field astrometry, the target is close to the phase tracking center and the antennas nod every few minutes between the target and a calibrator. If no special calibration provisions are taken, under typical conditions, an astrometric accuracy of ~10% of the synthesized beam FWHM can often be obtained. For example, an observation in Ka-band (~33 GHz) in A-configuration might reach an astrometric accuracy of ~10 milliarcseconds (mas). When care is taken (special calibrations and ideal observing conditions), the accuracy can approach 1–2% of the synthesized beam, with a floor of ~2 mas. If such accuracies are needed, we strongly recommend obtaining advice from VLA staff in setting up the observations.
Astrometric calibrators are marked J2000 A in the VLA calibrator list, and have a positional 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.
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.
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 ~12 GHz. Dynamic scheduling uses the API data to run a project under suitable conditions specified by the user. VLBI projects using the phased VLA will typically be fixed date and not dynamically scheduled.
Wide-field astrometry
Wide-field astrometry is used 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. With no special effort, the errors build up to roughly one synthesized beam at a separation of ~104 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 wide-band, wide-field case is a research area. These effects are handled somewhat differently in the post-processing packages. Check with VLA staff for more details via the NRAO Helpdesk.
Calibrating the Flux Density Scale
Normal calibration of the flux density scale for VLA observations is effected by including a scan on a source of presumed known flux density in each Scheduling Block (SB). Using that known flux density source, the flux density of the complex gain calibrator(s) can be determined and then transferred to your target source(s). Historically, 3C48 and 3C286 have been the standard sources for which NRAO has assumed flux densities are known as a function of frequency, and which have been recommended as flux density scale calibrator sources for the VLA. Restrictions on baseline length as a function of VLA configuration and observing band were supplied which, if followed, allowed relatively accurate flux density scale calibration. We have recently improved the ability to calibrate the flux density scale by providing sky brightness models for these sources in CASA and AIPS, which loosens the restrictions on configurations and bands. We have also added the sources 3C138** and 3C147 to the list of calibrators that have models. However, 3C48*, 3C138**, and 3C147 have spectral flux densities that vary with time (3C286, along with 3C295 and 3C196, are constant), so some care should be taken if the most accurate flux density scale calibration is desired.
Note: While accurate models are available in both AIPS and CASA for various frequency bands for the calibrators 3C286, 3C48*, 3C147, and 3C138**, neither 3C295 nor 3C196 has such models in CASA. Therefore, the VLA CASA calibration pipeline will fail if these two calibrators are used. Furthermore, 3C295 and 3C196 may not be suitable for all VLA configurations and frequencies even if one chooses to not use the pipeline.
A single observation of a few minutes of one of the above-mentioned flux density scale calibrators will suffice for most observers. If possible, the flux density scale calibrator should be observed at a time when it is nearly at the same elevation as the complex gain calibrator, especially for the highest four bands (Ku-, K-, Ka-, and Q-band). This is not always possible because of timing and geometry of sources, and that it is not typically known when an SB will be executed (so elevations versus time are uncertain). Flux density scale calibration accuracy in this case should be of order 10% at 4- and P-bands, 5% at L- through Ku-bands, and 10-15% for the three higher bands. If more accuracy is needed, a more careful strategy should be adopted, potentially using multiple flux density scale calibrators. The fundamental accuracy of the scale is ~5% at 4- and P-bands, 3% at L- through Ku-bands, increasing to 5% at Q-band. See Perley and Butler (2017) for more details on how the spectral flux densities of 3C48*, 3C138**, 3C147, and 3C286 (and many other sources) have been determined across the frequency range from 50 MHz to 50 GHz and how they vary versus time, along with information on the fundamental accuracy of the flux density scale when using these sources.
If less accuracy is needed in the flux density scale calibration, an observation of one of these standard sources need not necessarily be included in an SB. As an example, for a short triggered observation where a simple detection is desired, the time spent slewing back and forth to the flux density scale calibrator can make the SB significantly longer than it could otherwise be. In this scenario, the switched power measurement can be used to calibrate the flux density scale; see EVLA Memo 120 for some background. This technique, which should only be used with the 8-bit samplers, is not a standard path of calibration, but it is possible. The flux density scale accuracy in this case is ~10% for L- through Ku-bands, increasing to ~20% at Q-band; not nearly as good as using the "standard" method of flux density scale calibration, but it may be sufficient for some observers.
For reference, the polynomial expression for the spectral flux density for 3C286 determined in Perley and Butler (2017) is: \[\log(S) = 1.2481 - 0.4507 \log(f) - 0.1798 \log^2(f) + 0.0357 \log^3(f)\] where S is the flux density in Jy, and f is the frequency in GHz.
The tables below show flux densities determined using the polynomial coefficients for a few sources at a single frequency within each of the VLA bands.
Source | 75 MHz | 350 MHz | 1500 MHz | 3000 MHz | 6000 MHz | 10000 MHz | 15000 MHz | 22000 MHz | 33000 MHz | 45000 MHz |
---|---|---|---|---|---|---|---|---|---|---|
3C48* = J0137+3309 | 72.8 | 42.2 | 15.4 | 8.44 | 4.42 | 2.68 | 1.79 | 1.22 | 0.815 | 0.601 |
3C138** = J0521+1638 | 26.5 | 16.1 | 8.25 | 5.44 | 3.39 | 2.33 | 1.72 | 1.28 | 0.949 | 0.761 |
3C147 = J0542+4951 | 58.0 | 52.3 | 21.0 | 12.0 | 6.45 | 3.99 | 2.73 | 1.93 | 1.39 | 1.13 |
3C196 = J0813+4813 | 129 | 44.4 | 13.6 | 6.98 | 3.38 | 1.91 | 1.20 | 0.763 | 0.473 | 0.329 |
3C286 = J1331+3030 | 30.0 | 25.9 | 14.6 | 9.91 | 6.39 | 4.50 | 3.37 | 2.54 | 1.88 | 1.49 |
3C295 = J1411+5212 | 124 | 58.4 | 21.2 | 11.0 | 5.06 | 2.70 | 1.60 | 0.970 | 0.571 | 0.385 |
Source | 328 MHz | 1465 MHz | 2565 MHz | 4885 MHz | 6680 MHz | 11320 MHz | 16564 MHz | 25564 MHz | 32064 MHz | 48064 MHz |
---|---|---|---|---|---|---|---|---|---|---|
3C48* = J0137+3309 | 43.9 | 15.6 | 9.82 | 5.48 | 4.12 | 2.56 | 1.86 | 1.33 | 1.11 | 0.816 |
3C138** = J0521+1638 | 15.9 | 8.26 | 6.00 | 4.00 | 3.23 | 2.24 | 1.69 | 1.25 | 1.06 | 0.821 |
3C147 = J0542+4951 | 53.9 | 21.4 | 13.8 | 7.88 | 5.91 | 3.67 | 2.61 | 1.82 | 1.53 | 1.14 |
3C196 = J0813+4813 | 46.5 | 13.8 | 8.14 | 4.22 | 3.00 | 1.67 | 1.08 | 0.656 | 0.508 | 0.313 |
3C286 = J1331+3030 | 25.8 | 14.6 | 10.9 | 7.33 | 5.97 | 4.12 | 3.15 | 2.30 | 1.92 | 1.44 |
3C295 = J1411+5212 | 61.1 | 21.6 | 12.8 | 6.42 | 4.45 | 2.33 | 1.43 | 0.819 | 0.596 | 0.405 |
We refer the reader to the VLA Observing Guide for the practical considerations (e.g., observing frequency and array configuration, as well as post processing) regarding the choice of the flux density scale calibrator in their scheduling blocks.
* The flux density scale calibrator 3C48 has been undergoing a flare since January 2018 or so. While we have not fully characterized this with the VLA, other instruments have measured it at some frequencies. At Ku-band the magnitude of the flare is of order 10%. The effect will be smaller at lower frequencies (of order 5% at L-band), and might be larger at higher frequencies (of order 20% at Q-band). If you care about the flux density scale of your observations at that level, you may want to re-calibrate your data once new time-variable values have been put into CASA and AIPS.
** The flux density scale calibrator 3C138 is currently undergoing a flare. From VLA calibration pipeline results, we have noticed that 3C138 is deviating from the model. The amount of this deviation is still being investigated by NRAO staff, but does seem to effect frequencies of 10 GHz and higher. At K and Ka-bands the magnitude of the flare is currently of order 40-50% compared to Perley-Butler 2017 flux scale. If you care about the flux density scale of your observations above 10 GHz, monitoring datasets are publicly available in the archive under project code TCAL0009, from which you may find an updated flux density ratio to use for your data.
Complex Gain Calibration
General Guidelines for Complex Gain Calibration
Adequate complex gain calibration (tracking amplitude and phase fluctuations as a function of time) is a complicated function of source-calibrator separation, frequency, array scale (configuration), and weather. Since what defines adequate for some experiments is completely inadequate for others, it is difficult to define simple guidelines to ensure adequate phase calibration. However, some general statements remain valid most of the time. These are given below.
- Under decent conditions with no thunderstorms or ionospheric storms, tropospheric effects dominate at frequencies higher than about 4 GHz; ionospheric effects dominate at frequencies lower than about 4 GHz.
- 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. For these cases, calibration need only be done to track instrumental changes—three or four times per hour is usually sufficient for tracking the system gains.
- If your target object has sufficient flux density to permit phase self-calibration, there is no need to calibrate more than a couple of times per hour at low frequencies or 15 minutes at high frequencies in order to track pointing or other effects that might influence the amplitude scale. The enhanced sensitivity of the VLA guarantees, for full-band continuum observations, that every field will have enough background sources to enable phase self-calibration at L and S-bands. At higher frequencies, the background sky is not sufficient, and only the flux of the target source itself will be available.
- In principle, the smaller the source-calibrator angular separation, the better (even if the closer calibrator is weaker). However, the final choice will depend on the observing frequency. If deciding between a nearby calibrator with an S code in the calibrator database, and a more distant calibrator with a P code, for low frequencies (L-band and below) the nearby calibrator is usually the better choice (low frequency strategy). For higher frequencies it is advisable to use the further away but higher quality P-code calibrator (high frequency strategy). A detailed description of calibrator codes is available in the calibrator list.
- Phase stability often deteriorates dramatically after about 10AM due to small-scale convective cells set up by solar heating, even in clear and calm conditions, especially in the summer. Observers should consider a more rapid calibration cycle for observations between this time and a couple of hours after sundown.
- At high frequencies, and longer configurations, rapid switching between the source and nearby calibrator is needed to track tropospheric phase fluctuations if the target cannot be self-calibrated. See Rapid Phase Calibration and the Atmospheric Phase Interferometer (API) (below).
Rapid Phase Calibration and the Atmospheric Phase Interferometer
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 or a single observation of an unpolarized source. The complex gain 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 15–50 degrees during a several 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 (schedule four scans, if possible, in case one is lost); if at all possible, it is strongly recommended that five or more observations covering 100 degrees (or more) of parallactic angle in roughly uniform steps be run. 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 small objects as 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 can be used at frequencies of ~3 GHz and higher, or 3C147 at frequencies over ~10 GHz. Note that 3C48 and 3C138 are variable—the polarization properties are known to be changing significantly over time, most notably at the higher frequencies (for details see Perley and Butler (2013b)).
More information on polarization calibration strategy can be found in the Polarimetry section in the Guide to Observing with the VLA.
VLBI Observations
The VLA can participate in VLBI observations with the VLBA. This is possible by either using the VLA as a phased array (Y27) or using one of its dishes (single dish: Y1). We note that currently P-band cannot be phased. For more details see the VLBI at the VLA documentation. In phased array mode, the program TelCal derives the antenna-based delay and phase corrections needed for antenna phasing in real time. This correction is applied to the antenna signals before they are summed, requantized to 2-bits, and recorded in VDIF format on the Mark5C disk 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.
Standard VLA data, i.e., correlations between VLA antennas, are also archived in the NRAO science data archive. By default, a maximum of 512 MHz dual polarization is phased/recorded and sent to the archive. This leaves a large portion of the possible VLA bandwidth unused (up to 2 or 8 GHz total depending on sampler choice). Increasing the bandwidth to the maximum that the VLA can provide has obvious benefits for most continuum observations. Pulsar gating can be performed commensally with the phased VLA, avoiding the need for extra single dish observations. Line VLBI observations that run through the WIDAR correlator currently produce each requested VLBI subband bandwidth divided into full polarization products with 64 frequency channels, regardless of the requested spectral resolution obtained in the VLBI correlation. By adding unused baseline boards to the VLBI-specified line subbands, a closer match in spectral resolution can be offered for the standalone VLA data. If you think your science could benefit from these capabilities, please review our documentation on VLBI at the VLA or contact the helpdesk for further guidance.
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 and longer wavelengths in the D configuration. For instance, at 20 cm, 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. The deconvolution of the data 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://sundog.stsci.edu/top.html) (B configuration) or the NRAO VLA Sky Survey (NVSS, http://www.cv.nrao.edu/nvss/) (D configuration, all-sky).
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 when using the default inputs with CASA tasks importasdm and importevla. AIPS task UVFLG can be used to flag VLA data based on shadowing, although it will only flag based on antennas in the dataset, and is ignorant of antennas in other subarrays. The CASA task flagdata can also be used to flag data based on shadowing. For more information on shadowing, please see the Antenna Shadowing section in the Guide to Observing with the VLA.
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. The affected data would show time-variable high-amplitude points.
Combining Configurations and Mosaicking
Any single VLA configuration will allow accurate imaging of a range of spatial scales determined by the shortest and longest baselines. For extended and structured objects, it may be required to obtain observations in multiple array configurations. It is advisable that the frequencies used be the same for all configurations to be combined. The ideal combination of arrays results in a uv-plane with all cells equally filled by uv-points. To first order, this can be achieved by using the beam sizes of the individual arrays to inversely scale the on-source integration time. This approach is equivalent to achieving the same surface brightness sensitivity for all arrays on all scales. For the VLA, observations in the different configurations generate beam sizes that decrease by factors of three, i.e., C configuration generates a three times smaller beam than D configuration, B three times smaller than C, and A three times smaller than B. Thus, on-source integrations would increase by about an order of magnitude between each array. Such a drastic increase is very expensive and, in fact, not necessary since some spatial scales are common to more than a single array, which is equivalent to some uv-cells being filled more than others. The best way to fill the uv-plane depends on many factors such as declination of the source, LST time of the observation, and bandwidth.
Experience shows for the VLA that a factor of about three in on-source integration time for the different array configurations works well for most experiments. For example, a 20min on-source time in D, 1hr in C, 3hrs in B, and 9hrs in A should produce a decent map. Using large bandwidths and multi-frequency synthesis will broaden all uv tracks radially and one may need even fewer array configurations or shorter integration times between the different arrays.
Objects larger than the primary antenna pattern may be mapped through the technique of interferometric mosaicking. The VLA has no limit on the number of pointings for each mosaic. Typically hexagonal, rectangular, or individual pointing patterns are used and the overlap regions will result in an improved rms over each individual pointing. Given the many, potentially short observations, it is important to obey the data rate limits outlined in the Time Resolution and Data Rates Section. In addition to discrete or pointed mosaics, on-the-fly (OTF) mosaics (i.e. dumping the data while moving the telescopes across the source) are also available.
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 of the Documentation for more information.
Guidelines for mosaicking with the VLA are given in the Guide to Observing with the VLA.
Pulsar Observing
The VLA can be used for several kinds of pulsar observing: phase-binning using the WIDAR correlator, using the phased-array for single-beam pulsar processing in either search or fold modes, or simply standard imaging mode with fast integrations. Both phase-binning and phased-array (YUPPI) modes are available under General Observing (GO). The only exception is the 4-band YUPPI which is a Resident Shared Risk Observing (RSRO) capability. For any questions not addressed here regarding the capabilities of these observing modes, please contact the NRAO Helpdesk.
Phased-array pulsar processing
The "Y" Ultimate Pulsar Processing Instrument (YUPPI) is a software suite that runs in the correlator backend (CBE) computer cluster and can process a single-beam phased/summed-array data stream for pulsar observations in real time, into either folded profiles or search mode (filterbank) output. Coherent dedispersion can be optionally applied in either mode.
In the phased-array pulsar processing mode, the voltage data streams from each antenna are divided into a number of frequency subbands within the correlator, then summed and requantized before being output to the cluster for pulsar processing. The limitation on bandwidth comes primarily from the available network connections between the correlator and cluster. In all cases, a maximum of 64 subbands total can be processed. Depending on the number of bits chosen, this results in the following total bandwidth constraints:
Subband bandwidth |
Subband quantization |
Max total bandwidth |
Samplers |
---|---|---|---|
32 MHz | 8 bits | 2048 MHz | 8-bit |
64 MHz | 4 bits | 4096 MHz | 3-bit |
128 MHz | 2 bits | 8192 MHz | 3-bit |
<32 MHz | 8 bits | 64*BWsub | 8-bit |
As described in the VLA Frequency Bands and Tunability section, the 8-bit samplers provide two independently tunable 1 GHz IFs, while the 3-bit samplers provide four tunable 2 GHz IFs.
The pulsar-specific processing is done in real time using the DSPSR software package and, in principle, any processing option supported by DSPSR can be used; this will be constrained by the real-time computing power available in the cluster. In general, each subband can be divided into an arbitrary (2n) number of channels; 1 (summed), 2 or 4 detected polarization products can be output; and coherent dedispersion can be enabled or not.
Fold mode
In fold mode, the data are averaged modulo a known pulsar ephemeris (provided via a standard TEMPO/TEMPO2 "par file") into pulse profiles. The data can also be folded at a constant topocentric period, for example at 10 Hz to detect the injected noise cal signal. Fold integration times as short as 1 second have been tested. Up to 16384 profile bins can be used. The data are recorded in PSRFITS format using the standard 16-bit data encoding. This means the final output data rate is given by:
Data rate = 2 bytes × Nsubband × Nchannel × Nbin × Npoln / Tint
If the desired data rate exceeds ~25 MB/s, additional testing ahead of time may be required.
Search mode
In search mode, the data are simply detected and averaged over a specified amount of time before being output to disk, resulting in a filterbank data array (power vs time and frequency). Coherent dedispersion at a known DM can optionally be enabled for this. Data can be recorded using 2, 4, 8, 16 or 32 bits, resulting in a final data rate of:
Data rate = (Nbit/8) bytes × Nsubband × Nchannel × Npoln / Tint
The maximum sustained output rate in this mode should be kept less than ~400 MB/s.
Subarrays
It is possible to use any of the phased-array pulsar modes listed here in a subarray observation, following the guidelines described in the Subarrays section. In addition, the above constraints on the pulsar processing apply to the total of all simultaneously-used subarrays, rather than each subarray separately. For example, the total number of subbands in use across all subarrays must not be greater than 64; the total (not per-subarray) data rate must meet the above constraints, etc. It is possible to use different parameters such as subband bandwidth, number of bits, or processing mode (fold versus seach) in the different subarrays.
VLBI
It is possible to use phased-array pulsar processing as part of a VLBI experiment; see the VLBI Observations section, and links therein, for additional information about VLBI at the VLA. The main constraint on this type of observation is that a subband that is being recorded for VLBI can not be sent to the pulsar processing system. However, since VLBI recording typically only requires a small number of subbands (2 through 8), any additional subbands produced by WIDAR can be sent to the pulsar system, following the constraints above. This provides a high time resolution data stream covering wider bandwidth than the VLBI data. One typical use case is to detect a pulsar using the VLA data stream and determine a short-term timing ephemeris covering the observation. This can then be used to gate the VLBI correlation, reducing uncertainties associated with extrapolating existing timing solutions or obtaining time on other telescopes for these purposes.
Gated or binned visibilities
The WIDAR correlator has the capability to internally integrate (fold) visibilities into 1 or more pulse phase bins, following a standard TEMPO-compatible pulsar ephemeris. This mode can be used to image the emission from a pulsar of known period anywhere in the telescope field of view. This provides both higher signal-to-noise ratio on the pulsar than a standard image, and allows the pulsed emission to be separated from continuous emission from other sources in the field.
The following constraints apply to binning-mode observations:
- Binning is limited to the case where the pulse period is divided evenly into bins covering the full pulse period. "Gating" style observations (common in VLBI) where a single on-pulse bin is used are not supported.
- The maximum number of pulse phase bins is 1000.
- There is a tradeoff between the total bandwidth and the minimum bin width (pulse period divided by number of bins):
- With 4 subbands (up to 512 MHz total), the minimum allowed bin width is 12.5 μs.
- With 16 subbands (up to 2048 MHz total), the minimum allowed bin width is 50 μs.
- With 64 subbands (up to 8192 MHz total), the minimum allowed bin width is 200 μs.
- The number of channels per subband is currently limited to 128 maximum. Combining recirculation and binning is not allowed.
- Integration (dump) time must be an integer number of pulse periods.
- The data rate produced in this mode is the standard VLA data rate (see the Data Rate section) multiplied by the number of bins. The data rate must be kept less than 60 MB/s.
It should also be noted that there is currently very limited support for binned observations in standard data processing software (e.g., CASA). Development of data analysis procedures is ongoing and users of this mode should be aware that this will likely involve some advanced/low-level manipulation of raw VLA data sets.
Fast-dump visibilities
While not specifically a pulsar mode, standard visibility data can be dumped as fast as 5 ms, which may be sufficient for imaging of slow pulsars. See the Time Resolution and Data Rates section for more details.
Solar Observing
Observations of the Sun can currently be performed in five frequency bands: 1-2 GHz (L band), 2-4 GHz (S band), 4-8 GHz (C band), 8-12 GHz (X band), and 12-18 GHz (Ku band). To observe the Sun, any 8-bit resource using these bands can be used in Solar mode (see below). As the Sun is moving across the sky, the source position should be defined in the source catalog.
If the center of the solar disk needs to be tracked, select the Sun using the drop-down for "Solar System Body with Internal Ephemeris". Otherwise, to allow for specifying the location of interest on or near the solar disk and the differential rotation model to be used for tracking the source, solar observations require selection of a source position type of "Solar System Body with Uploaded Ephemeris". An ephemeris file can be generated using the guidelines given in the Moving Object section or, for example, by using the ALMA Solar Ephemeris Generator.
The user must select the “Solar” scan mode when specifying scan details under Observation Preparation in the OPT. This scan mode ensures that the necessary attenuators are switched into the signal path when observing the Sun and that special noise cal sources are employed in the switched power system; this allows users to flux calibrate their data. Two observing modes are currently supported:
- For quiet Sun observations, when no strong active regions or flares expected, users should select a scan mode of “Solar Attenuators with Low Noise Internal Cal”.
- If the science objective involves observing a strong active region and/or flare activity, users should select a scan mode of “Solar Attenuators with High Noise Internal Cal”.
These selections ensure that cal noise levels of these scans are of order a few percent of the anticipated system temperature. (A third scan mode, “Noise Reverse Coupler Setup (L Band only)”, is not yet supported.)
Solar observations rely on the switched power system to flux calibrate the data. This, in turn, requires the use of 8-bit sampling modes in the frequency bands that support solar observing. Solar observing is such that input power levels may change from scan to scan, particularly for active Sun targets. This being the case each source scan requires a setup scan to trigger the system to reset the requantizer gain.
Considerations for the calibration of solar observations are otherwise quite similar to standard gain calibration. Phase calibration proceeds in the same manner for solar observations as it does for general observing. However, it is important to be aware of two factors: First, when observing a calibrator source, the attenuation used when observing the Sun is removed from the signal path; the RFI that is suppressed, to a large degree, when observing the Sun is fully present in observations of calibrators. Second, when the attenuators are removed, solar radiation may be present in the sidelobes of the primary beam if the phase calibrator is too close to the Sun, particularly if the Sun is in an active phase. Therefore, it is advisable to observe phase calibrators that are further in angular distance from the source than would normally be used.
VLA+LWA (eLWA) Observing
The VLA can be used to record individual antenna signals (voltages) to VLBI Data Interchange Format (VDIF) files. Similar to Very Long Baseline Interformetry, using the phased-VLA or individual antennas, this specifically allows for joint correlation of VLA antennas with the Long Wavelength Array (LWA) stations in New Mexico in the overlapping 4m-band range. Joint correlation is performed offline using the LWA software correlator, which is located inside the VLA control building and produces FITS-IDI compatible data output.
For shared risk observing, this mode is available for a single subband with a center frequency of 76 MHz, a bandwidth of 8 MHz, and 4-bit VDIF output. Other possible modes or center frequencies (limited by the VLA MJP-dipole response) are available through resident-shared risk observing. Use of any RSRO modes should also be brought to the attention of the LWA Director in order to verify that correlation is possible.
eLWA proposals that are granted VLA time are automatically granted time for the LWA stations, for more information on how to select the eLWA resource for your proposal, refer to the PST manual. Observations are prepared and scheduled like any other VLA observations through the OPT and the LWA stations are automatically triggered to follow the VLA pointing. Specific instructions for eLWA instrument setups are provided in the OPT manual.
The two LWA stations currently available in this mode are: LWA1 (close to the center of the VLA) and LWA-Sevilleta (on Sevilleta National Wildlife Refuge, ~80 km North-East of the VLA).
Current limitations
- VLA+LWA (eLWA) observing cannot be mixed with other observing modes in a single scheduling block/observing script. This limitation is due to the disabling of the array geometric delay model.
- Data quality is highly susceptible to the low-frequency interference environment. Known sources of interference can be the active sun, powerline arcing, or self-generated interference from AC power components or digital electronics. The environment is regularly monitored and mitigation measures are taken, which in rare cases can significantly delay or prevent execution of observations.
- Correlated data products are available through the LWA archive only and will eventually also be available through the regular NRAO archive. The PI of a proposal will be notified when correlated products are available.
- Data calibration is recommended to be performed using AIPS.
Connect with NRAO