VLA > VLA Observational Status Summary 2018B

# VLA Observational Status Summary 2018B

VLA capabilities October 2018 - January 2019

# 1. Introduction

### Purpose of Document, Older Versions of the OSS

This Observational Status Summary (OSS) summarizes the instrumental status of the Karl G. Jansky Very Large Array (VLA) for the C-configuration, for the observing period 4 October 2018 through 7 January 2019 (Semester 2018B), and should be used when preparing proposals for the 1 February 2018 deadline. Multi-configuration proposals that include this configuration may also be submitted. Additionally, proposals requesting only configurations that will fall in semester 2019A (or later) may be submitted if the Principal Investigator is a graduate student. NRAO offers this service to provide scientific and technical feedback for students, and to provide them with an opportunity to re-submit their proposals for their principal semester with this information in hand.

The OSS is intended as a ready reference for those contemplating use of the VLA for their astronomical research. The information is in summary form; those requiring greater detail should use the NRAO Helpdesk, or refer to the manuals and documentation listed in Documentation. Most of the information contained here, and much more, is available through the VLA science web page and the companion VLBA science web page. For capabilities offered in previous semesters, we refer to our overview of all OSS versions available online.

The VLA is a large and complex modern instrument. Some familiarity with the principles and practices of its operation is necessary for its efficient use. Although the NRAO strives to make using the VLA as simple as possible, users must be aware that proper selection of observing mode and calibration technique is often crucial to the success of an observing program. Inexperienced and first–time users are encouraged to enlist the assistance of an experienced colleague or NRAO staff member for advice on, or direct participation in, an observing program. For more details, refer to the Visiting the DSOC and VLA page. The VLA is an extremely flexible instrument, and we are always interested in imaginative and innovative ways of using it.

### Overview of the VLA

The Karl G. Jansky Very Large Array (VLA) is a 27–element interferometric array, arranged along the arms of an upside-down Y, which produces images of the radio sky at a wide range of frequencies and resolutions. The VLA is located at an elevation of 2100 meters on the Plains of San Agustin in southwestern New Mexico, and is managed from the Pete V. Domenici Science Operations Center (DSOC) in Socorro, New Mexico.

The basic data produced by the VLA are the visibilities, or measures of the spatial coherence function, formed by correlation of signals from the array's elements. The most common mode of operation will use these data, suitably calibrated, to form images of the radio sky as a function of sky position and frequency. Another mode of observing, commonly called phased array, allows operation of the array as a single element through coherent summation of the individual antenna signals. This mode is most commonly used for Very Long Baseline Interferometry (VLBI) observing and for observations of rapidly varying objects, such as pulsars.

The VLA can vary its resolution over a range exceeding a factor of ∼50 through movement of its component antennas. There are four basic antenna arrangements, called configurations, whose scales vary by the ratios 1 : 3.28 : 10.8 : 35.5 from smallest to largest. These configurations are denoted D, C, B, and A, respectively. For details about antenna positions in the various configurations we refer to the stations position file (pdf version).

The VLA completes one cycle through all four configurations in an approximately 16 month period. Consult the Configuration Plans and Proposal Deadlines page or recent NRAO and AAS newsletters for current and up-to-date configuration schedules and associated proposal deadlines. Refer to the Guide to Proposing for the VLA for information on how to submit an observing proposal.

Observing projects on the VLA will vary in duration from as short as 1/2 hour to as long as several weeks. Most observing runs have durations of a few to 24 hours with only one or, perhaps, a few target sources. However, since the VLA is a two-dimensional array, images can be made with data durations of less than one minute. This mode, commonly called snapshot mode, is well suited to surveys of relatively strong, isolated objects. See the section on Snapshots for more detail.

All VLA antennas are outfitted with eight receivers providing continuous frequency coverage from 1 to 50 GHz. These receivers cover the frequency ranges of 1–2 GHz (L-band), 2–4 GHz (S-band), 4–8 GHz (C-band), 8–12 GHz (X-band), 12–18 GHz (Ku-band), 18–26.5 GHz (K-band), 26.5–40 GHz (Ka-band), and 40–50 GHz (Q-band). Additionally, all antennas of the VLA have receivers for lower frequencies, enabling observations at P-band (230–470 MHz). These low frequency receivers also work at 4-band (54–86 MHz), and new feeds have been deployed on a number of VLA antennas to observe at this frequency range.

The VLA correlator is both powerful and flexible. Details of the correlator configurations being offered for VLA science are described in the WIDAR Section of this document. It is important to realize that the VLA correlator is fundamentally a spectral line correlator and that even continuum observations are done in a wide-band mode with many channels.

# 2. Offered VLA Capabilities during the Next Semester

### The Call for Proposals

The most recent Call for Proposals summarizes the General Observing (GO) capabilities being offered for the Karl G. Jansky Very Large Array (VLA).

In addition to these general capabilities, NRAO continues to offer shared risk observing options for those who would like to push the capabilities of the VLA beyond those offered for general use. These are the Shared Risk Observing (SRO) and Resident Shared Risk Observing (RSRO) programs.

Details about what is being offered for each program is given below. If you have any questions or problems with any link or tool, please submit a ticket through the NRAO Helpdesk.

Considering the lack of future hybrid configurations after semester 2016A, guidelines on how to substitute such configurations with the use of principal array configurations are presented in the Array Configurations section of the Guide to Proposing for the VLA.

### General Observing (GO) and Shared-Risk Observing (SRO)

#### Summary of Capabilities

As described in the Call for Proposals, the VLA offers continuous frequency coverage from 1–50 GHz in the following observing bands: 1–2 GHz (L-band); 2–4 GHz (S-band); 4–8 GHz (C-band); 8–12 GHz (X-band); 12–18 GHz (Ku-band); 18–26.5 GHz (K-band); 26.5–40 GHz (Ka-band); and 40–50 GHz (Q-band). Both single pointing and mosaics with discrete, multiple field centers will be supported under General Observing (GO). P-band (224–480 MHz) Stokes I continuum and spectral line observations are also offered under GO; see the Low Frequency Observing section below for more details. Data rates of up to 60 MB/s (216 GB/hour) will be available to all users, combined with correlator integration time limits per band and per configuration, as described in the Time Resolution and Data Rates section. However, data rates in excess of 25 MB/s (90 GB/hour) require additional justification. Limitations on frequency settings and tuning ranges are described in the Frequency Bands and Tunability section.

The GO capabilities being offered are:

 Capability Description 8-bit samplers Standard default setups for: 2 GHz bandwidth continuum observations at S/C/X/Ku/K/Ka/Q bands (16 × 128 MHz subbands) 1 GHz bandwidth continuum observations at L band (16 × 64 MHz subbands) 256 MHz bandwidth continuum observations at P band (16 × 16 MHz subbands) Flexible setups for spectroscopy using two, independently tunable, 1 GHz baseband pairs, each of which can be split into up to 16, flexibly tunable, subbands Single, dual, and full polarization products Number of channels summed over all polarization products up to 16384 (no recirculation) or up to 1,048,576 (with recirculation; subject to the data rate limits and the chosen subband bandwidth) 3-bit samplers Standard default setups for: 8 GHz bandwidth continuum observations at K/Ka/Q bands 6 GHz bandwidth at Ku band 4 GHz bandwidth at C/X bands Flexible setups for spectroscopy using four, independently tunable, 2 GHz baseband pairs, each of which can be split into up to 16, flexibly tunable, subbands Single, dual, and full polarization products Number of channels summed over all polarization products up to 16384 (no recirculation) or up to 1,048,576 (with recirculation; subject to the data rate limits and the chosen subband bandwidth) Mixed 3-bit and 8-bit samplers Allows more flexibility for simultaneous continuum and high-resolution spectral line observing Subarrays Up to 3 independent subarrays using standard 8-bit continuum setups Phased-array for VLBI VLA Phased Array (Y27) for VLBI (see the VLBA call for proposals and the VLBA Observational Status Summary for a detailed description of the VLA Phased Array capabilities being offered)

SRO capabilities can be set up via the Observing Preparation Tool (OPT) and run through the dynamic scheduler without intervention, but are not as well tested as GO capabilities. A summary of the SRO capabilities being offered are:

• On-the-Fly (OTF) mosaicing (used when each pointing on the sky is no more than a few seconds);
• Up to 32 subbands per baseband with the 8-bit samplers;
• Solar observing;
• P-band (224–480 MHz) spectroscopy; and
• Frequency averaging in the correlator: Averaging channels by either a factor of two or four in the correlator is offered to reduce the data volumes for all wide-band continuum science projects that utilize the 3-bit samplers (applicable to the VLA bands C through Q). This capability is only available for single subarray and non-OTF observations.

‡Note: OTF mosaic observations have certain limitations on the length and number of scans. Contact the VLA staff through the NRAO Helpdesk to check if the OTF mosaic observations to be proposed are technically feasible ahead of submitting a proposal.

We expect that most SRO programs will have no or only minor problems that can be corrected quickly. If an SRO program fails, however, and it becomes clear that detailed testing with additional expertise is needed, then the project must make an experienced member from their team available to help troubleshoot the problem. In some cases, this may require the presence of that experienced member in Socorro. If adequate support from the project is not given, then the time on the telescope will be forfeited. The additional effort is to be determined based on discussions with the NRAO staff and management and the project team.

The guidelines for General and Shared Risk observing proposals, along with information about tools and other advice, can be found in the VLA Proposal Submission Guidelines.

### Resident Shared Risk Observing (RSRO)

#### Summary of Capabilities

The RSRO program provides access to extended capabilities of the VLA that require additional testing in exchange for a period of residence to help commission those capabilities. Capabilities that would fall under the RSRO program include:

• Correlator dump times shorter than 50 msec, including integration times as short as 5 msec for transient detection;
• Pulsar observations;
• Data rates above 60 MB/s;
• Recirculation beyond a factor of 64 in the correlator;
• P-band system (224–480 MHz) polarimetry;
• 4-band system (54–86 MHz; see the Low Frequency Observing section below);
• More than 3 subarrays or subarrays with the 3-bit system; and
• Complex phased array observations (e.g., pulsar and complex VLBI observing modes).

The guidelines for Resident Shared Risk Observing proposing, along with requirements and considerations, can be found in the VLA Proposal Submission Guidelines.

### Low Frequency Observing (P-band and 4-band)

The new low frequency receiver system developed in collaboration with the Naval Research Laboratory is available for observations at P-band (224–480 MHz) through the GO and SRO program. Use of the P-band system for polarimetry is through the RSRO program as noted above. The new receivers also work at 4-band (54–86 MHz), and new feeds are expected to be deployed on all VLA antennas by the time of the 2018 semester. Both 4-band and P-band can be observed simultaneously, but access to the 4-band system is only available through the RSRO program.

### Commensal Observing Systems at the VLA

There are two commensal systems on the VLA that may take data at the same time as your proposed observation. The first is the VLITE system, which will take data at P-band during regular observations that use bands other than P-band. The VLITE system is deployed on fifteen VLA antennas. Observers wishing to gain access to the commensal VLITE data taken during their VLA observations should follow the instructions on the VLITE web page for doing so. The second is the realfast system, which takes data at very fast dump rates in an effort to detect Fast Radio Bursts (FRBs). This system is still in development, but may be observing during the 2018B semester.

To report errors or problems encountered in any link or while using any NRAO tool listed here, please submit a ticket through the NRAO Helpdesk.

# 3.1. 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.

Table 3.1.1: Configuration Properties
ConfigurationABCD
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.

# 3.2. Sensitivity

### Sensitivity

The theoretical thermal noise expected for an image using natural weighting of the visibility data is given by:

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.62TsysA, 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.

Table 3.2.1: SEFDs and D-Configuration Confusion Limits
Frequency SEFD
(Jy)

RMS confusion level
in D config (µJy/beam)
Robust Weighting

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) 350 negligible
22 GHz (K) 560 negligible
33 GHz (Ka) 710 negligible
45 GHz (Q) 1260 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:

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):

where D is the distance to the galaxy in Mpc, and SΔV is the HI line area in units of Jy km/s.

# 3.3. VLA Frequency Bands and Tunability

### Bands

For observations taken with the 8-bit samplers, each receiver can tune to two different frequencies, each 1024 MHz 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 (X and Y). 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 (or Y when applicable), IFs C and D provide LCP (or X when applicable). 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. In 8-bit mode, a maximum of 1024 MHz can be correlated for each IF pair (see the WIDAR Section), for a total maximum bandwidth of 2048 MHz. To distinguish this 8-bit system from the 3-bit system, these IF pairs are denoted A0/C0 and B0/D0.

More options are available with the 3-bit samplers. This system provides four (R,L) polarization pairs, each 2048 MHz wide. The A/C IF pair provides two sampled pairs, labelled A1/C1 and A2/C2, and the B/D IF pair provides two sampled pairs, labelled B1/D1 and B2/D2.

For more details on the 8-bit and 3-bit samplers see the VLA Samplers section.

The tuning ranges, along with default frequencies for continuum applications, are given in Table 3.3.1 below.

Table 3.3.1: Default frequencies for continuum applications
BandRange18-bit continuum applications (GHz)3-bit continuum applications (GHz)

(GHz)IF pair A0/C0IF pair B0/D0IF pair A1/C1IF pair A2/C2IF pair B1/D1IF pair B2/D2
4 m (4) 0.058 – 0.0842 .054 – .086
90 cm (P) 0.22 – 0.503 0.224 – 0.4803
20 cm (L) 1.0 – 2.04 1.0 – 1.54 1.5 – 2.04
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 4.0 – 6.0 6.0 – 8.0
3 cm (X) 8.0 – 12.0 8.0 – 9.0 9.0 – 10.0 8.0 – 10.0 10.0 – 12.0
2 cm (Ku) 12.0 – 18.0 13.0 – 14.0 14.0 – 15.0 12.0 – 14.0 14.0 – 16.0 16.0 – 18.0
1.3 cm (K) 18.0 – 26.5 20.2 – 21.2 21.2 – 22.2 22.0 – 24.0 24.0 – 26.0 18.0 – 20.0 20.0 – 22.0
1 cm (Ka) 26.5 – 40.0 32.0 – 33.0 31.0 – 32.0 33.0 – 35.0 35.0 – 37.0 29.0 – 31.0 31.0 – 33.0
0.7 cm (Q) 40.0 – 50.0 40.0 – 41.0 41.0 – 42.0 44.0 – 46.0 46.0 – 48.0 40.0 – 42.0 42.0 – 44.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 the NRAO Helpdesk for further information.
2. The 4-band system is currently under development. Observing time may be requested through the RSRO program.
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.

# 3.4. 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 2048 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 4096 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.
• 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.

# 3.5. Field of View

### 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 = 42/νGHz for frequencies between 1 and 50 GHz (L- through Q-band). At P-band the approximate value is θPB = 50/νGHz. New precise measurements of the primary beam shape have been reported in EVLA Memo 195; these allow for the correction of the primary beam attenuation in wide-field images. Both AIPS and CASA (5.0 and later versions) have these new parameters incorporated.

With the wide-bandwidths of the VLA it is necessary to account for the variation of the primary beam with frequency in order to achieve high-dynamic range images. For this and other imaging details we refer to the Limitations on Imaging Performance section of the OSS.

Objects larger than approximately half this angle cannot be directly observed by the array. However, a technique known as mosaicking, 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 mosaicking with the VLA are given in the Guide to Observing with the VLA.

### 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 (θ0HPBW). Table 3.5.1 shows the decrease in peak response and the increase in apparent radial width as a function of this parameter and should be used to determine how much spectral averaging can be tolerated when imaging a particular field.

Table 3.5.1: Reduction in Peak Response Due to Bandwidth Smearing
(Δν/ν0)*(θ0HPBW) 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; θ0HPBW is the source offset from the phase tracking center in units of the synthesized beam.

Note: The VLA correlator supports frequency averaging for single subarray and non-OTF observations under the Shared Risk Observing program. Currently this capability is limited to an averaging by a factor of 2 or 4  and only for wide-band continuum science projects that utilize the 3-bit samplers (applicable for C-band through Q-band observations). Observers interested in this capability should consult the EVLA memo 199 to assess the suitability of the frequency averaging in the correlator for their observations, because the extent of the bandwidth smearing is heavily dependent on the frequency averaging factor, the array configuration, and the observing frequency.

### 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 3.5.2, 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.

Table 3.5.2: Averaging Time for a Given Amplitude Loss
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) results 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.

Note: For both the chromatic aberration and the time-averaging loss, the issue is not a simple reduction in amplitude for sources far from the phase center, but a convolution to the extent that a point source far from the phase center will become resolved due to bandwidth and/or time smearing. Furthermore, the description given above for the bandwidth smearing is based on the assumption that the radiation is monochromatic to parameterize the smearing, and does not take into account the consequences of having wide-bandwidths as is the case for the VLA. Therefore, while proposing and planning for VLA observations, and depending on the objectives of the science and the location of the sources of interest within the field, including confusing sources which may be far outside the science field, the above noted guidelines need to be used to conservatively estimate the proper channel width and the correlator integration time in order to minimize the effects of the bandwidth smearing and the time smearing, respectively.

### 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 VLA 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 cm 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 task IMAGR and CASA task clean. Refer to the package help files for these tasks, or consult with the NRAO Helpdesk 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 for more details).

# 3.6. Time Resolution and Data Rates

The default integration times for the various array configurations and frequency bands are as follows:

Table 3.6.1: Default Integration Times
ConfigurationsObserving
Bands
Default
integration time
A, B, C, D 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:

Table 3.6.2: Minimum integration times and maximum data rates
Proposal type

Minimum integration time

Maximum data rate
General Observing (GO) and
Shared Risk Observing (SRO)
50 msec 25 MB/s (90 GB/hr), or up to
60 MB/s (216 GB/hr) with additional justification
Resident Shared Risk Observing (RSRO) < 50 msec > 60 MB/s (216 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, which is offered under the Shared Risk Observing program, 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 General Observing Setup Tool (GOST), or the Resource Catalog Tool (RCT) portion of the Observation Preparation Tool (OPT).

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.

# 3.8. Subarrays

The continuum subarray option offers two 1 GHz baseband pairs with the 8-bit samplers 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 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.

Table 3.8.1: Some Possible Subarray Options
Number of antennas correlated using each Baseline Board columnNumber
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.

# 3.9. 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.

# 3.10. Limitations on Imaging Performance

### Imaging Fidelity

Image fidelity is a measure of the accuracy of the reconstructed sky brightness distribution. A related metric, dynamic range, is a measure of the degree to which imaging artifacts around strong sources are suppressed, which in turn implies a higher fidelity of the on-source reconstruction.

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 effective phase stability of the telescope due to atmospheric/ionospheric fluctuations, although pointing errors and changes in atmospheric opacity may also be a limiting factor. If a good model of the sky brightness distribution exist (e.g., use of compact structures of sufficient strength—though a good model of resolved sources in the field of view may also be used), standard self-calibration can be counted on to improve the images. At low frequencies, where the dominant phase error is due to ionospheric plasma density fluctuations, more advanced techniques may be required to account for change of ionospheric phase across the field of view. Depending on the underlying nature of the errors, dynamic ranges in the thousands to hundreds of thousands can be achieved using these techniques. With the new WIDAR correlator, with its much greater bandwidths and much higher sensitivities, self-calibration methods can be extended to observations of sources with much lower flux densities than ever possible with the old VLA.

The choice of image reconstruction algorithm also affects the correctness of the on-source brightness distribution. The CLEAN algorithm is most appropriate for predominantly point-source dominated fields. Extended structure is better reconstructed with multi-resolution and multi-scale algorithms. For high dynamic ranges with wide bandwidths, algorithms that model the sky spectrum as well as the average intensity can yield more accurate reconstructions.

### Invisible Structures at Large Angular Scales

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. Image reconstruction algorithms such as multi-resolution and multi-scale CLEAN may help reduce these negative bowls, but care must be taken in choosing appropriate scale sizes to work with.

Table 3.1.1 in the Resolution section gives the largest scale visible to each configuration/band combination.

### 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.

### Sidelobes from non-Deconvolved 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. Sidelobe confusion is a strong function of observing band—affecting most strongly L and P-band observations—and is rarely a significant problem for observations at frequencies above 4 GHz.

### 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.

### Wide-band Imaging

The very wide bandpasses provided by the VLA enable imaging over 2:1 bandwidth ratios: at L, S, and C-bands, the upper frequency is twice that of the lower frequency. It is this wide bandwidth which enables sub-microJy sensitivity.

In many cases, where the observation goal is a simple detection and there are no strong sources near to the region of interest, standard imaging methods that combine the data from all frequencies into one single image (multi-frequency-synthesis) may suffice. A rough rule of thumb is that—provided a strong source is not adjacent to the target zone—if the necessary dynamic range in the image is less than 1000:1 (i.e., the strongest source in the beam is less than 1000 times higher than the noise), a simple wide-band map may suffice.

For higher dynamic ranges, complications arise from the fact that the brightness in the field of view dramatically changes as a function of frequency, both due to differing structures in the actual sources in the field of view and due to the attenuation of the sources by the primary beam. One symptom of such problems is the appearance of radial spokes around bright sources, visible above the noise floor, when imaged as described above.

The simplest solution is to simply make a number of maps (e.g., one for each subband) which can then be suitably combined after correction for the primary beam shape. But with up to 64 subbands available with the VLA's new correlator, this is not always the optimal approach, both from a practical standpoint and frequency-dependent image reconstruction uncertainties. Further, images at all bands must be smoothed to the angular resolution at the lowest frequency before any spectral information can be extracted. With a 2:1 bandwidth, the difference in angular resolution across the band will be significant.

A better approach is to process all subbands simultaneously, utilizing software which takes into account the possibility of spatially variant spectral index and curvature and knows the instrumentally-imposed attenuation due to the primary beam. Such wideband (plus wide-field) imaging algorithms are now available within CASA as part of the clean and tclean tasks

### Wide-field Imaging

Wide-field observing refers to both the non-coplanar nature of the VLA when observing in non-snapshot mode and the effect of the primary beam patterns.

At high angular resolutions and low frequencies, standard imaging methods will produce artifacts around sources away from the phase center due to the non-coplanar nature of the VLA . Faceted imaging (AIPS, CASA) and W-projection (CASA) techniques can be used to solve this problem.

Wide field imaging also requires the accurate representation of primary beam patterns and their use during imaging. This is relevant only for very high dynamic ranges ( > 10000:1 ) or when there are very strong, confusing sources at and beyond the half-power point of the primary beam. This problem is worse with a wide-band instrument simply because the size of the primary beam—and the radius at which the half-power point occurs—varies with frequency while there is increased sensitivity out to a wider field of view. In CASA, the tclean task contains algorithms to deal with these effects by modeling and correcting for polarization squint, frequency-dependent and rotating primary beams, per antenna, during imaging. Work is under way to develop full-polarization primary beam corrections.  Please note, however, that most advanced methods will lead to a significant increase in processing time and may not always be required. Therefore, in the interest of practicality, they should be used only if there is evidence of artifacts without these methods. At low frequencies, wide-field imaging typically includes both non-coplanar effects plus primary beams whereas at high frequencies it is primarily primary beam effects.

Mosaicking is another form of wide-field imaging in which data from multiple pointings are combined either during imaging as a joint mosaic (CASA) or after imaging by stitching several images together (CASA, AIPS).

# 3.11. 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.

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 is not a standard path of calibration, and is not easily supported in CASA or AIPS, 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 table below shows flux densities determined using the polynomial coefficients in Perley and Butler (2017) for a few sources at a single frequency within each of the VLA bands. Since several of these sources are variable with time, these values are appropriate for January 2016.

Flux densities (Jy) of Standard Calibrators for January 2016
Source75 MHz350 MHz1500 MHz3000 MHz6000 MHz10000 MHz15000 MHz22000 MHz33000 MHz45000 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

# 3.12. 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.
• The smaller the source-calibrator angular separation, the better. However, when 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. 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

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 for very small source-calibrator separations. Observing efficiency declines, however, 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. This is represented in the Observation Preparation Tool as 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. Short source-calibration scans have been demonstrated to result in images of faint sources with diffraction-limited spatial resolution on the longest baselines. Under average weather conditions, and using a 120 second cycle time, the residual phase at 43 GHz should be reduced to ≤ 30 degrees. Note that at a typical wind velocity in the compact D-configuration, this effective baseline length is the same as—or larger than—the longest baseline in the array and it is not worth the increased overhead of short cycle times. Under these circumstances, it is sufficient to calibrate every 5-10 minutes to track the instrumental changes. The fast switching technique will not work in bad weather (such as rain showers or when there are well-developed convection cells (thunderstorms)). It is important to correctly specify the required tropospheric phase stability as measured by the Atmospheric Phase Interferometer at observe time (see below).

Further details can be found in VLA Scientific Memos # 169 and 173. These memos, and other useful information, can be obtained from References 9 and 10 in Documentation.  Also see the High Frequency Strategy guide for additional recommendations on observing at high frequencies.

An Atmospheric Phase Interferometer (API) is used to continuously measure the tropospheric contribution to the interferometric phase. The API uses an interferometer of two, 1.5 meter antennas, separated by 300 meters, observing an 11.7 GHz beacon from a geostationary satellite. The API data are heavily used for the dynamic scheduling of the VLA.

Characteristic seasonal averages are shown in Table 3.12.1 below:

Table 3.12.1: Seasonal API/wind values at the VLA
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

Day indicates sunrise to sunset values; night indicates sunset to sunrise values.

### Other Issues that Affect Complex Gain (Amplitude)

There are other instrumental effects that cause fluctuations in gain over time for VLA observations - we describe three of them here.  Note that all of these effects are to the gain amplitude, not phase.

#### Gain Curves

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 per antenna and per band, and the necessary corrections can be applied to the visibility data in both AIPS and CASA. Additionally, atmospheric opacity will also cause an elevation-dependent gain which is also particularly notable at these four highest frequency bands. We currently do not have an atmospheric opacity monitoring procedure; users should utilize the appropriate tasks available in both AIPS and CASA to estimate and correct for the opacity using ground-based weather data. If your complex gain calibrator is near your target source, and the flux density scale calibrator is also observed at a similar elevation, then most of this elevation-based gain will be calibrated correctly during normal calibration. Note also that a good procedure for removing elevation-based 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 incorrect antenna gain curves. Use of this procedure is only practical if your observations span a wide range in elevation.

#### Antenna Pointing Offsets

Another important gain variation effect at the four highest frequency bands is that due to antenna pointing offsets. Daytime observations on sunny days can suffer pointing errors, primarily in elevation, of up to one arcminute. This effect can be largely removed by utilizing the referenced pointing procedure which determines the pointing offset of a nearby calibrator. This offset is then applied to subsequent target source observations. It is recommended that this local offset be determined at least hourly, more often during sunrise or sunset, utilizing an object within 15 degrees of the target source—preferentially at an earlier right ascension. Studies show that the maximum pointing error will be reduced to about 7 arcseconds or better if proper referenced pointing is utilized.

#### Electronic Gain Variations

The VLA's post-amplifiers are not temperature stabilized and exhibit significant gain (amplitude) changes between night and day, particularly at the four highest frequency bands. Changes as large as 30% have been seen between night and day in calm, clear conditions. These gain changes, and others caused by possible changes in attenuator settings, are monitored and can be removed by application of the internal calibration signal, whose results are recorded in the switched power table in both AIPS and CASA. These corrections are not applied by the calibration pipeline—users who wish to correct for these gain changes must utilize the appropriate tasks in AIPS or CASA. For the most accurate flux density bootstrapping, this table must be applied to the visibility data before calibration.

# 3.13. 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 3C138 is 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.

# 3.14. VLBI Observations

The VLA can participate in VLBI observations with the VLBA. Presently, this is only allowed in phased array mode (single dish is only available through the VLBA Resident Shared Risk Observing program) with restricted WIDAR correlator configurations. Also 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.

# 3.15. 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).

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.

# 3.17. 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 less 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) is 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.

# 3.18. 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. Any of these types of pulsar observing are considered Resident Shared Risk Observing (RSRO), and participants in this program are expected to work closely with NRAO staff. 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 suite of software 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 this 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:

Table 3.18.1: Pulsar Observing 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, although 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, 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 also 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.

Use of this mode will likely require testing ahead of the observations to verify the desired setup works correctly. Some constraints on what we expect to be possible include:

• 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 (512 MHz total), the minimum allowed bin width is 12 μs.
• With 16 subbands (2048 MHz total), the minimum allowed bin width is 50 μs.
• With 64 subbands (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.

# 4. WIDAR Correlator

## Introduction

The correlator configurations offered for general observing may be divided into three basic modes: wideband, spectral line, and subarrays. 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 and in the Resources section of the Proposal Submission Tool (PST), which must be used to define the correlator configuration for General Observing (GO) and Shared Risk Observing (SRO) proposals.

Note that for general observing phased array configurations are only allowed as part of VLBI experiments (see the section on VLBI Observations); other uses such as pulsar observations must be done as Resident Shared Risk Observations.

Wideband and spectral line observing modes with the WIDAR correlator are described below. For the subarray mode, we refer to the Subarrays section of the OSS.

For technical details about the WIDAR correlator, refer to References 14 in Documentation.

## WIDAR Correlator: Wideband Observing

The wideband observing setups provide the widest possible bandwidth for a given observing band, with channel spacing depending on the number of polarization products as listed in the following table 4.1.1:

Table 4.1.1: Wideband & Subarray Correlator Options
(all but P and L-bands)
Polarization productsChannel spacing
Full (RR, RL, LR, LL) 2 MHz
Dual (RR and LL) 1 MHz
Single (RR or LL) 0.5 MHz

8-bit wideband setups are available for all observing bands, providing a total of 2 GHz of bandwidth per polarization (1 GHz per polarization at L-band, and 256 MHz per polarization at P-band). 3-bit setups are available for all bands above S-band, providing total bandwidths per polarization of 4 GHz (C/X-bands), 6 GHz (Ku-band), or 8 GHz (K/Ka/Q-bands). In all cases, except for P and L-band, each of the subbands is 128 MHz wide. At L-band the default is 64 MHz/subband, yielding channels twice as narrow as those listed in the table above, while at P-band the default is 16 MHz/subband, resulting in 125 kHz channel spacing.

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, or choose to have a mix 8-bit and 3-bit samplers. The tuning restrictions are described in the section on VLA Frequency Bands and Tunability, and the 8-bit and 3-bit samplers are described in the section on VLA Samplers.

## WIDAR Correlator: Spectral Line Observing

### Basebands and Subbands

Observers have access to very flexible correlator configurations using up to 64 subbands in up to 4 basebands sampled with the 8-bit and/or the 3-bit samplers. These capabilities may be summarized as follows:

• Two 1 GHz baseband pairs using the 8-bit samplers, or four 2 GHz baseband pairs using the 3-bit samplers, independently tunable within the limits outlined in the section on VLA Frequency Bands and Tunability. The 8-bit baseband pairs are referred to as A0/C0 and B0/D0, while the 3-bit samplers are A1/C1, A2/C2, B1/D1, and B2/D2. The AC/BD nomenclature corresponds to that of the IF pairs in the pre-expansion VLA.
• Up to 16 subband pairs (spectral windows) in each 3-bit baseband pair, and up to 32 subbands in each 8-bit baseband pair, for a total of up to 64 subbands in any correlator configuration:
• 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 16384 (without recirculation):
• These may be spread flexibly over subbands and polarization products, in multiples of 64: 64, 128, 192, 256, 384, …, 16384 cross-correlation products;
• Recirculation may be used to increase the number of channels per subband for subbands narrower than 128 MHz. Baseline Board stacking may be used to increase the number of channels per subband for setups requiring less than 64 subbands;
• Assigning many channels to a given subband may reduce the total bandwidth and/or 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 (, …, 15) (any integer between 0 and 7 for 8-bit, between 0 and 15 for 3-bit); νsbLow the lower edge of the subband (the subband center frequency minus half the subband bandwidth); νsbHigh the upper edge of the subband (the subband center frequency plus half the subband bandwidth).

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. Note in particular that the center of a baseband is a boundary and no line should be observed at the baseband center.

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 and 128 MHz subbands are thus constrained to cover a region between two of those boundaries, with no finer tuning being 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.)

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 the upper sideband, and at the highest sky frequency in the baseband when using the lower sideband.

### Subband Bandwidths and 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 percent within a few percent of the subband edge. The precise filter shape and noise increase is a complex but predictable function of the subband bandwidth (sbBW) 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, without recirculation enabled, can have a different number of channels and polarization products, subject to two limitations:

1. For the ithsubband, 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 and LL)
• 256 nBlBP,i with a single polarization product (RR or LL)
Here nBlBP,i= 1, 2, 3, 4, 5, …, 64 is the number of Baseline Board Pairs (BlBPs) assigned to that subband.
2. The sum over all subbands of nBlBP,i must be less than or equal to 64, the number of Baseline Board pairs in the correlator. Equivalently, the sum over all subbands of spectral channels times polarization products is limited to 64 × 256 = 16,384 (without recirculation).

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. Use of more than one pair per subband (i.e., nBlBP>1) is known as Baseline Board stacking; see additional details about this below.

Limitation #1 corresponds to table 4.2.1 of the options for subband bandwidth and spectral resolution when using nBlBP Baseline Board pairs for a subband:

Table 4.2.1: Subband Bandwidth and Spectral Resolution Options (without recirculation)
Subband bandwidth &
total velocity coverage
Full polarization products
(RR, RL, LR, LL)
64nBlBP spectral channels

Channel spacing:
Dual polarization products
(RR and 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/nBlBPGHz km/s 1000/nBlBP kHz 300/nBlBPGHz km/s 500/nBlBP kHz 150/nBlBPGHz 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.0325 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.

• nBlBP is the number of Baseline Board Pairs assigned to the subband.
• Each subband may have a different number of spectral channels and polarization products, and each may be tuned independently.
• There can be at most 16 subbands per baseband, and nBlBP must be an integer: 1, 2, 3, 4, 5, …, 64.
• The sum of nBlBP over all subbands must be less than or equal to 64.
• Use of more than one BlBP for a subband may further restrict the number of subbands available in one or more of the basebands; see text for details.

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'nProducts Spectralchannels 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-5 RR,RL,LR,LL 64 6 x 1 A0/C0 sb6 RR, LL 3840 1 x 30 A0/C0 sb7 RR 768 1 x 3 A0/C0 sb8 RR,RL,LR,LL 192 1 x 3 B0/D0 sb0-2 RR,RL,LR,LL 64 3 x 1 B0/D0 sb3 LL 768 1 x 3 B0/D0 sb4 RR, LL 2048 1 x 16

### Recirculation

When the subband bandwidth is less than the maximum 128 MHz, the spare clock cycles that become available in the correlator hardware can be re-purposed to compute additional lags using a single baseline board. This increases the spectral resolution within the subband, and is known as recirculation. Each factor of two reduction in subband bandwidth results in an additional factor of two maximum lags; therefore for subbands of 128 MHz / N the possible spectral resolution (in units of frequency) can be increased by a factor of N2.

#### Recirculation vs. Baseline Board Stacking

When faced with the choice between recirculation and Baseline Board stacking to increase the number of channels in a subband, we recommend recirculation for subbands narrower than 128 MHz; this is supported in observatory software (GOST, OPT). Using recirculation rather than stacking frees up more Baseline Boards for other user; alternatively the experiment becomes less dependent on all Baseline Board pairs being available/working at the time of observation. For subbands of 128 MHz, recirculation is not possible, and Baseline Board stacking must be utilized to increase the number of channels.

The present implementation of recirculation is that, for each halving of the subband bandwidth, the number of channels in the subband may be doubled without having to use additional correlator hardware. The maximum recirculation factor for a subband is 128/(subband bandwidth in MHz) and, of course, subject to other configuration restrictions such as data rate.

### Baseline Board Stacking

As opposed to recirculation, which increases the number of channels in a subband by exploiting otherwise unused CPU resources, Baseline Board stacking adds more channels to a subband by adding correlator hardware resources, i.e., using up more Baseline Board pairs. Using Baseline Board stacking may therefore 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. The following section is for those who are attempting complex line experiments and who find the GOST or the RCT restricting the number of subbands and/or channels they can use in unexpected ways. Most observers can skip the following section.

#### Baseline Board Stacking and Correlator Use

First let us consider how the correlator hardware is organized. The cross-multiplications in the WIDAR correlator are spread across 64 Baseline Board pairs (BlBP), arranged into 4 quadrants of 16 BlBP each. Each baseband is connected directly to one of those quadrants. In the simplest mode, each of the 16 BlBP of a quadrant handles the correlations for one of the 16 subbands of the corresponding baseband. Four basebands and four quadrants are required to handle the full 8 GHz of bandwidth per polarization provided by the 3-bit (wideband) samplers: 8 GHz is split into four basebands of 2 GHz each, with each baseband fed into a different BlBP quadrant. Each BlBP in that quadrant handles a subband of maximum bandwidth 128 MHz, so 16 BlBP handles 16 subbands for a total of 16×128 MHz = 2048 MHz.

A single BlBP produces 256 cross-correlations per baseline for a single subband, which can be used for a single polarization product (e.g., RR or LL with 256 spectral channels), or two (RR and LL with 128 spectral channels each), or four (RR, RL, LR, and LL with 64 spectral channels each).

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 8×128 MHz = 1024 MHz of each baseband. Three-quarters of the correlator BlBP hardware remain unused.

The spectral line mode allows access to these extra correlator resources through Baseline Board stacking: using multiple BlBPs 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 BlBPs. Those BlBPs can then be used to produce more spectral channels for that subband, with n BlBPs producing 256×n cross-correlations per baseline. The limit on the total number of cross-correlations (16384) stems from the total number of BlBPs (64): 64×256 = 16384.

Unfortunately, completely flexible crossbar switches are expensive and could not be implemented in the VLA's correlator. This means that one cannot route a given subband to a randomly-chosen BlBP. The routings which are possible, are as follows:

1. A subband in a baseband can be routed to any BlBP within the corresponding quadrant.
2. Data coming into a given BlBP in one quadrant, can be routed to the corresponding BlBP in any other quadrant.

Routing option #1 means that one could use all the BlBPs within a quadrant to correlate a single subband, yielding 16×256 = 4096 cross-correlations for that subband:

Routing option #2 means that one could use the BlBPs in all 4 quadrants to correlate a single subband. One simple case would use 4 BlBPs to correlate each of the 16 subbands in a single baseband, yielding 4×256 = 1024 cross-correlations for each of those subbands. Note that in this case, no BlBPs are left to correlate any data from the second baseband.

Using routing option #2 does come with a subtle cost: assigning a BlBP in quadrant X to correlate a subband corresponding to quadrant Y removes that BlBP from use in the baseband corresponding to quadrant X…and therefore 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 BlBPs), 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 31×128 MHz = 3968 MHz per polarization. Actually, however, there are only 15+15 subbands available, or 30×128 MHz = 3840 MHz per polarization, because the spectral line subband has eaten one BlBP corresponding to the other baseband:

If the same spectral line required twice as many channels, this will result in the loss of two subbands in both of the basebands:

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:

Understanding these confusing constraints can help observers set up the VLA more effectively to achieve their scientific goals. For instance, in 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 above examples all use BlBP pair stacking in powers of 2, but this is not required. To give some idea of more complex possibilities, the following tables (4.2.3 and 4.2.4) give two examples of other possible configurations. The RCT display shows how the Baseline Boards are used to process the individual subbands. The cyan boxes (shaded when printed out in black and white) 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.

Table 4.2.3: Complex Configuration Example #1
BasebandSubbandPol'n productsSpectral channelsnBlBP
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
RCT display:

Table 4.2.4: Complex Configuration Example #2
BasebandSubbandPol'n productsSpectral channelsnBlBP
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
RCT display:

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 a bandwidth of 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.

Use of the 3-bit samplers further extends the possibilities. Here is one example:

Table 4.2.5: 3-bit Complex Configuration Example #1
A1/C1 sb0-8 RR, LL, RL, LR 9 x 64 9 x 1 Q1: 0–8
A1/C1 sb9 RR, LL 1 x 1152 1 x 9 Q1 & Q3: 9–11, 14 / Q4: 9
A1/C1 sb10 RR 1 x 1792 1 x 7 Q1 & Q3 & Q4 : 12,13 / Q2: 13
A1/C1 sb11 RR, LL 1 x 384 1 x 3 Q1 & Q2 & Q3: 15
A2/C2 sb0-11 RR, LL, RL, LR 12 x 64 12 x 1 Q2: 0–11
A2/C2 sb12 LL 1 x 768 1 x 3 Q2: 12, 14 / Q4: 14
B1/D1 sb0-3 RR, LL, RL, LR 4 x 64 4 x 1 Q3: 0–3
B1/D1 sb4 RR, LL, RL, LR 1 x 320 1 x 5 Q3: 4–8
B2/D2 sb0-6 RR, LL, RL, LR 7 x 64 7 x 1 Q4: 0–6
B2/D2 sb7 RR, LL 1 x 640 1 x 5 Q4: 7, 8, 10, 11, 15
RCT display:

Once again, the GOST implements all of these constraints, and is generally smart enough to figure out the routing scheme that works best for your particular request.

# 5. Documentation

### Documentation

Documentation for VLA data reduction, image making, observing preparation, etc., can be found in various manuals. Current manuals are available on-line. Those manuals marked by an asterisk (*) can be mailed out upon request, or are available for downloading from the NRAO website. Direct your requests for mailed hardcopy to Lori Appel. Many other documents of interest to the VLA user, not listed here, are available from our website.

1. PROCEEDINGS FROM THE 1988 SYNTHESIS IMAGING WORKSHOP: Synthesis theory, technical information and observing strategies can be found in: "Synthesis Imaging in Radio Astronomy." This collection of lectures given in Socorro in June 1988 has been published by the Astronomical Society of the Pacific as Volume 6 of their Conference Series. The lectures of the 2014 workshop are available at the 14th Synthesis Imaging Workshop web site.
2. PROCEEDINGS FROM THE 1998 SYNTHESIS IMAGING WORKSHOP: This is an updated and expanded version of Reference 1, taken from the 1998 Synthesis Imaging Summer School, held in Socorro in June, 1998. These proceedings are published as Volume 180 of the ASP Conference Series.
3. GUIDE TO OBSERVING WITH THE VLA: Describes details of how to observe with the VLA once you have been allocated time on the VLA (https://science.nrao.edu/facilities/vla/docs/manuals/obsguide). Including special observing modes such as:
1. OBSERVING WITH THE 8-BIT (up to 2 GHz bandwidth) & 3-BIT (up to 8 GHz bandwidth) SAMPLER SYSTEMS (https://science.nrao.edu/facilities/vla/docs/manuals/obsguide/modes/set-up);
2. SPECTRAL LINE OBSERVING (https://science.nrao.edu/facilities/vla/docs/manuals/obsguide/modes/line);
3. HIGH FREQUENCY OBSERVING (https://science.nrao.edu/facilities/vla/docs/manuals/obsguide/topical-guides/hifreq);
4. LOW FREQUENCY OBSERVING (https://science.nrao.edu/facilities/vla/docs/manuals/obsguide/topical-guides/lofreq);
5. VERY LOW FREQUENCY OBSERVING (< 500 MHz) (https://science.nrao.edu/facilities/vla/docs/manuals/obsguide/topical-guides/vlofreq);
6. POLARIMETRY (https://science.nrao.edu/facilities/vla/docs/manuals/obsguide/modes/pol);
7. MOSAICKING (https://science.nrao.edu/facilities/vla/docs/manuals/obsguide/modes/mosaicking);
9. MOVING OBJECTS (https://science.nrao.edu/facilities/vla/docs/manuals/obsguide/modes/moving);
10. VLBI AT THE VLA (https://science.nrao.edu/facilities/vla/docs/manuals/obsguide/modes/vlbi).
4. *CASA COOKBOOK: The CASA Cookbook for use of the package for data reduction of VLA (& ALMA) data is available, along with other documentation, from the CASA home page (http://casa.nrao.edu). See (http://casa.nrao.edu/docs/cookbook/)
5. VLA CASA Guides: Tutorials and data reduction examples of VLA data in CASA (https://casaguides.nrao.edu/index.php/Karl_G._Jansky_VLA_Tutorials)
6. *AIPS COOKBOOK: The Astronomical Image Processing System (AIPS) software is able to fully calibrate VLA data and do most imaging operations. The exception is the wide-band (bandwidth synthesis) deconvolution which is being developed in CASA only. ALMA data may also be reduced in AIPS although the package is not fully qualified to calibrate data from the ALMA linearly-polarized feeds. The Cookbook description for calibration and imaging under the AIPS system can be found near all public workstations in the SOC. The latest version has expanded descriptions of data calibration imaging, cleaning, self-calibration, spectral line reduction, and VLBI reductions. See (http://www.aips.nrao.edu/cook.html)
7. *GOING AIPS: This is a two-volume programmers manual for those wishing to write programs under AIPS. It is now somewhat out of date. See (http://www.aips.nrao.edu/goaips.html)
8. *VLA CALIBRATOR LIST: This page contains the list of VLA Calibrators in both 1950 and J2000 epoch. See (https://science.nrao.edu/facilities/vla/observing/callist)
9. *The Very Large Array: Design and Performance of a Modern Synthesis Radio Telescope, Napier, Thompson, and Ekers, Proc. of IEEE, 71, 295, 1983.
10. *HISTORICAL VLA MEMO SERIES: archive memo series from the early days of the VLA. See (http://library.nrao.edu/vlam.shtml)
11. *RECENT VLA MEMO SERIES: the memo series relating to the expanded capabilities of the VLA. See (http://library.nrao.edu/evla.shtml)
12. *The VLA Expansion Project: Construction Project Book. The Expanded VLA Project Books contains the technical details of the VLA Expansion construction project. It is available online at http://www.aoc.nrao.edu/evla/pbook.shtml.
13. INTRODUCTION TO THE NRAO VERY LARGE ARRAY (Green Book): This manual has general introductory information on the VLA. Topics include theory of interferometry, hardware descriptions, observing preparation, data reduction, image making and display. Major sections of this 1983 manual are now out of date, but it nevertheless remains a useful source of information on much of the VLA. There are a few hard copies at the VLA and in the DSOC. Much of this document is now available for download (https://science.nrao.edu/facilities/vla/obsolete/green-book). Note: it does not include any information about the hardware and software specific to the expanded Karl G. Jansky VLA.
14. WIDAR: The DRAO design and development documents of the WIDAR correlator of the VLA are available at http://www.aoc.nrao.edu/widar/docs/.

### Online Tools & Important Links

The NRAO User Portal. (https://my.nrao.edu) This is a gateway to the NRAO interactive services that include the Proposal Submission Tool (PST).

The NRAO Proposal Submission Tool (PST) online manual. (https://science.nrao.edu/facilities/vla/docs/manuals/proposal-guide/pst)

The VLA Exposure Calculator Tool (ECT) online manual. (https://science.nrao.edu/facilities/vla/docs/manuals/propvla/determining)

The VLA Exposure Calculator Tool (ECT). (

The VLA GO Setup Tool (GOST). (https://science.nrao.edu/facilities/vla/docs/manuals/propvla/gost)

### Acknowledgements

Many thanks to all the VLA staff and our RSRO participants who have worked long and hard to commission these capabilities and who have helped to create this extensively updated set of documentation.

NRAO is grateful to Professor Rob Ivison for supporting the upgrade of some of the 3-bit samplers on the VLA via a grant from the European Research Council. For observations using the 3-bit samplers between May 2015 and March 2018 we encourage users to include the following text in the Acknowledgments section of their publications:

"We acknowledge funding towards the 3-bit samplers used in this work from ERC Advanced Grant 321302, COSMICISM."

# 6. Contact Information

Please go to the People page for information on key personnel at NRAO-Socorro.

Please direct queries to the NRAO Helpdesk; you can expect a response within one to two business days.

# 7. Editor's Notes

This Observational Status Summary for the Karl G. Jansky (expanded) VLA is based substantially on its predecessor, the VLA Observational Status Summary. Over the VLA history of almost 30 years, many individuals contributed to that document by writing sections, editing previous versions, commenting on draft material, and implementing the capabilities described herein. We thank all these contributors for their efforts. For questions on the content, or suggestions that would enhance the clarity of this guide, we recommend contacting the NRAO Helpdesk.