VLBA Observational Status Summary 2017B
VLBA capabilities August 2017 - January 2018
1. Introduction, Older Versions of the OSS
VLBA Observational Status Summary — 2017B
This document summarizes the current observational capabilities of the Very Long Baseline Array (VLBA) instrument. It is intended specifically to accompany the Call for Proposals for observing semester 2017B, with a submission deadline of 2017 February 1, but is also the best source for current information on VLBA instrumentation. For capabilities prior to that date, we refer to our overview of previous OSS versions available online.
The VLBA is an array of ten 25-m diameter antennas at stations distributed over United States territory (Napier et al. 1994; Napier 1995). It is the first astronomical array dedicated to observations using the technique of Very Long Baseline Interferometry (VLBI), which was pioneered in the 1960s. The VLBA offers (1) in absentia, year-round station and correlator operation; (2) station locations selected to optimize u-v plane coverage; (3) ten observing bands at wavelengths ranging from 90 cm to 3 mm (two stations are not equipped at 3 mm); (4) rapid, automated selection of receivers and of frequencies within a given receiver; and (5) integrated data flow from acquisition to correlation to post-processing. VLBA observations can acquire simultaneous dual circular polarizations from any single receiver, from widely separated frequencies within the 6-cm band, and from receiver pairs at 13/4 cm or 90/50 cm. The VLBA is operated remotely from the Pete V. Domenici Science Operations Center (SOC) in Socorro, New Mexico.
Broad overviews of the range of astronomical research possible with the VLBA are presented in the VLBA 10th anniversary meeting proceedings (Romney & Reid 2005), and the conference proceedings edited by Zensus, Taylor, & Wrobel (1998). Recommended reading for users new to the VLBA includes a "Guide to Using the VLBA", and a short VLBI overview (Walker 1999).
This document's primary intent is to provide, in concise form, the minimal information needed to formulate technically sound proposals requesting VLBA resources. A secondary aim is to describe some of the subtleties of data reduction and telescope scheduling. It is updated synchronously with the VLBA calls for proposals, or more often when required by major changes.
Requests for information beyond the scope of this document should be directed to the NRAO Helpdesk.
2. New Developments
The VLBA Sensitivity Upgrade project is now complete, with the VLBA's original observing system having been decommissioned at the end of 2013.
The Sensitivity Upgrade project traced its origins to a two-year internal NRAO study of the scientific goals (Romney 1999; VLBA Sensitivity Upgrade Memo # 2) and the derived technical requirements (Romney 2000; VLBA Sensitivity Upgrade Memo # 3), for an upgrade of the VLBA. These concepts were subsequently cited among the initial steps recommended in Phase II of the VLBI Future Committee report "Mapping the Future of VLBI Science in the U.S." (Taylor, Lonsdale, et al. 2004; VLBA Sensitivity Upgrade Memo # 4), upon which the initial Sensitivity Upgrade proposal (Walker et al. 2007; VLBA Sensitivity Upgrade Memo # 15) was based. Funding to begin implementation of these goals became available in 2007.
The Sensitivity Upgrade instrumentation now in full operational use addresses, finally, many of the VLBI Future Committee recommendations. The primary discrepancy lies in the achieved maximum recording rate being limited to 2 Gigabits/sec (Gbps), rather than the desired minimum 4-Gbps capability. It should be noted, however, that the Sensitivity Upgrade equipment was designed to support a 4-Gbps rate if sufficient resources can be found to buy a second recorder for each station, and to double the pool of recording media. The 2-Gbps rate, corresponding to a continuum bandwidth of 256 MHz per polarization, is nevertheless a 16-fold increase over the standard "sustainable rate" at which the VLBA operated for many early years.
Overviews of the new instrumentation are presented in this document in the sections on the ROACH Digital Backend and Mark 5C Recorder. More detailed information is available in the Sensitivity Upgrade memo series. Questions can be submitted via the NRAO Helpdesk.
3. Antenna Sites
Table 1 gives the geographic locations of the ten stations comprising the VLBA, ordered from East through West. All locations are based on the WGS84 ellipsoid used by the GPS system, with Earth radius [inline]a=6378.137\: \rm{km}[/inline] and flattening [inline]1/f=298.257223563[/inline]. See Napier (1995) for further site information.
North | West | |||
Latitude | Longitude | Elevation | Code | |
Location | [° ′ ″] | [ ° ′ ″] | [m] | |
Saint Croix, VI | 17:45:23.68 | 64:35:01.07 | -15 | SC |
Hancock, NH | 42:56:00.99 | 71:59:11.69 | 296 | HN |
North Liberty, IA | 41:46:17.13 | 91:34:26.88 | 222 | NL |
Fort Davis, TX | 30:38:06.11 | 103:56:41.34 | 1606 | FD |
Los Alamos, NM | 35:46:30.45 | 106:14:44.15 | 1962 | LA |
Pie Town, NM | 34:18:03.61 | 108:07:09.06 | 2365 | PT |
Kitt Peak, AZ | 31:57:22.70 | 111:36:44.72 | 1902 | KP |
Owens Valley, CA | 37:13:53.95 | 118:16:37.37 | 1196 | OV |
Brewster, WA | 48:07:52.42 | 119:40:59.80 | 250 | BR |
Mauna Kea, HI | 19:48:04.97 | 155:27:19.81 | 3763 | MK |
Several other radio telescopes often participate in VLBI observing in conjunction with the VLBA. A High Sensitivity Array (HSA) currently comprises the VLBA, Karl G. Jansky Very Large Array (VLA), Robert C. Byrd Green Bank Telescope (GBT), Arecibo, and Effelsberg telescopes; locations of the HSA telescopes are listed in Table 2. The HSA is described in detail in Section 14.2 of this document.
4. Frequency Bands & Performance
The nominal frequency ranges for VLBA receiver systems are shown in Table 3. Actual frequency ranges are broader; consult the measurements reported by Hronek & Walker (1996) for details. Updates on frequency-dependent performance across VLBA bands are available at http://www.vlba.nrao.edu/cgi-bin/wbd_dir.pl. These actual ranges may be especially important for avoiding radio frequency interference (RFI), and for programs involving extragalactic spectral lines, rotation measures (Cotton 1995b; Kemball 1999), and multi-frequency synthesis (Conway & Sault 1995; Sault & Conway 1999).
(c) Different ranges within the 20-cm receiver.
(f) Different ranges within the 1 cm receiver. Continuum performance better at 23.8 GHz, away from water line.
(g) See Table 4 below for individual station details.
(j) Data rate 256 Mbps. (k) Data rate 32 Mbps.
(l) 8-station array; 4-hour integration.
Also shown in Table 3 are parameters characterizing the performance of a typical VLBA station for the various frequency bands. Columns [3] and [5] give typical VLBA system-equivalent-flux-density (SEFD) values at zenith and opacity-corrected peak gains, respectively. These are means over measurements in both polarization at all ten antennas, at the frequencies in column [4].
The typical zenith SEFD can be combined with the aggregate recorded data rate and appropriate integration times to estimate the root-mean-square (RMS) noise level on a single VLBA baseline, and in a VLBA image. Characteristic values tabulated in columns [6] and [7] are computed assuming, for most cases, the VLBA's upgraded 2-Gbps recording rate for continuum observations; a typical fringe-fit interval of 2 minutes; and a total on-source integration time of 8 hours. Exceptions, indicated in the table notes, apply to the fringe-fit intervals at the lowest and highest frequency bands, where shorter intervals are often required; for the recording rate limits imposed by the available RF bandwidth at the lowest frequency bands; and for most parameters at the extreme 3-mm band. Performance may be worse than the tabulated estimates on some baselines due to poor primary or subreflector surfaces or poor atmospheric conditions.
Antenna | Nominal | Typical | Typical | Typical | Baseline |
Frequency |
Zenith | Peak | Zenith | Sensitivity |
|
Range | SEFD | Gain | Tsys |
ΔS2048,30s |
|
[GHz] | [Jy] | [K Jy-1] | [K] | [mJy] | |
BR | 80 - 90 | 3500 | 0.039 | 135 | 28. |
NL | 80 - 96 | 4900 | 0.055 | 270 | 33. |
FD | 80 - 96 | 3600 | 0.034 | 120 | 28. |
LA | 80 - 90 | 3100 | 0.051 | 160 | … |
PT | 80 - 96 | 4100 | 0.024 | 100 | 28. |
KP | 80 - 96 | 4600 | 0.025 | 110 | 30. |
OV | 80 - 96 | 5800 | 0.020 | 100 | 33. |
MK | 80 - 96 | 4100 | 0.023 | 100 | 28. |
The 3 mm band extends beyond the design specification for the VLBA antenna, and is challenging for the panel-setting accuracy of the primary reflectors, the figure of the subreflectors, and the pointing of the antennas. In addition, performance in this band is highly dependent on weather conditions. Table 4 gives a snapshot of performance at 86 GHz for each antenna, as well as the RMS noise in 30 seconds (at 2 Gbps) on a baseline to LA, which is one of the most sensitive3 mm antennas.
5. Station Signal Processing
5.1. Primary Signal Path
This sub-section describes the instrumentation that collects and amplifies radio-frequency (RF) radiation from a source, converts, and transmits it to the station control building. Napier et al. (1994) includes further information on most of the following aspects.
The antenna brings the RF signals to a focus at one of ten feeds. The main reflector is a 25-m diameter shaped figure of revolution with a focal-length-to-diameter ratio of 0.354. A 3.5-m diameter Cassegrain subreflector with a shaped asymmetric figure is used at observing wavelengths shorter than 30 cm, while the prime focus is used at longer wavelengths. The antenna features a wheel-and-track mount, with an advanced-design reflector support structure. Antenna motions are relatively rapid for an antenna of this size, to facilitate rapid source changes: 30° per minute in elevation and 90° per minute in azimuth.
The feed couples free-space electromagnetic waves into waveguides for transmission to the receiver system. Feeds at observing wavelengths shorter than 30 cm are located on a ring at the offset Cassegrain focus, and are selected by rotation of the subreflector with a maximum transition time of about 20 seconds. A frequency-selective dichroic system enables simultaneous 13/4-cm observations. The 90- and 50-cm feeds are crossed dipoles mounted on the subreflector near prime focus; simultaneous 90/50-cm observations are possible.
The polarizer extracts orthogonal circularly-polarized signals, which are routed separately to dual receivers. In receiver wavelength bands shorter than 30 cm, the polarizer is cooled to cryogenic temperatures.
The receiver amplifies the signal. Most VLBA receivers use HFET (Heterostructure Field Effect Transistor) amplifiers at a physical temperature of 15 K, but the 90- and 50-cm receivers use GaAsFETs (Gallium Arsenide FETs) at room temperature. All receivers produce dual-polarization outputs, in opposite hands of circular polarization.
The IF converter mixes the receiver output signals with the first LO generated by a front end synthesizer. Two IF bands between 512 and 1024 MHz are output by each converter, one in each sense of circular polarization. The same LO signal is used for mixing with both polarizations in most cases. However, the 4 cm IF converter has a special mode that allows both output signals to be connected to the RCP output of the receiver and to use separate LO signals, thereby allowing the use of spanned bandwidths exceeding 512 MHz. Also, the 90 cm and 50 cm signals are combined and transmitted on the same IFs. The 50 cm signals are not frequency converted, while the 90 cm signals are upconverted to 827 MHz before output.
Four IF cables, designated A, B, C, D, carry the IF signals from the antenna vertex room to the station control building. Normally only two IFs are in use at a time, with the signals from each IF converter transmitted via A and C, or B and D; by convention, RCP is normally carried by IFs A and B, and LCP by C and D. However, several 4-IF configurations are available for special cases. These include dual-polarization observations at two arbitrary frequencies anywhere within the 4-8 GHz range of the new 6-cm receiver, and combinations of dual IF outputs from both the 13- and 4-cm receivers (using the dichroic system described in the paragraph on feeds above). Either of these configurations can be activated simply by specifying the desired frequencies and polarizations. The 4-cm receiver also supports a 3-IF configuration, with a second first LO available for RCP signals; this mode requires that the SCHED setup include the 'dualx' parameter.
5.2. Frequency and Time Standard
Essential auxiliary instrumentation, required to make simultaneous observations feasible at VLBA stations separated by thousands of kilometers, is described in this sub-section.
A hydrogen maser provides an ultra-stable frequency reference at each VLBA station. Its standard signals, at 100 MHz and 5 MHz, and multiplied versions thereof, are used throughout the station electronics, both in the antenna and in the station building.
The front end synthesizer generates reference signals to convert the receiver output from RF to IF, with lock points at (n× 500) ± 100 MHz (for n an integer). Output frequencies range between 2.1 and 15.9 GHz. There are 3 such synthesizers, each of which is locked to the maser. One synthesizer is used for most frequency bands, but two are used at 1 cm, at 7 mm, 3 mm, and for the 4 cm wideband mode.
5.3. Calibration Signals
VLBA stations support several different types of calibration measurements.
Two calibration signals are injected near the beginning of the primary data path, and detected elsewhere in the VLBA system, with derived corrections applied in data analysis:
The switched-noise system injects well calibrated, broadband noise, switched on-off in a 50% duty cycle at 80 Hz. This noise signal is synchronously detected in the RDBE, to provide a time-tagged system-temperature table that is delivered with the primary fringe visibility data. Application of these measurements for amplitude calibration is discussed separately.
The pulse-cal system injects a series of pulses at intervals of 1.0 or 0.2 microseconds, to generate monochromatic, phase-stable tones, spaced at multiples of 1 MHz or 5 MHz. The tones are detected in the VLBA DiFX correlator. Application of these measurements for phase calibration is discussed separately.
There is also a round trip cable calibration scheme that monitors the length of the signal cables, to enable corrections for temperature and pointing induced variations.
5.4. Roach Digital Backend (RDBE)
The RDBE replaces much of the VLBA's original analog signal processing in the station control building. The baseband converters, in particular, are eliminated by sampling directly from the IF outputs of each station's receivers, with 8-bit precision. All subsequent processing is performed digitally. For clarity in the following descriptions, two items of essential VLBA terminology are defined here:
An "IF" refers to one of a maximum of four 512-MHz wide intermediate-frequency analog signals transmitted from the receiver(s) to the RDBE. Most receivers provide two IFs, in opposite circular polarizations. However, four IFs are available to support specialized observing modes at some wavelengths: two dual-polarization pairs, at arbitrary frequencies within the full range of the new 6-cm receiver; or from different receivers in 13/4-cm or 90/50-cm dual-receiver operation.
A "channel" refers to a single contiguous frequency range (of any bandwidth), observed in a single polarization, that is sampled, filtered, and recorded as a separate entity. This approach is essential for the VLBA, where capabilities are fundamentally limited by the overall data-transmission bandwidth.
'RDBE' is an acronym for "ROACH Digital Backend''. ROACH, in turn, refers to the FPGA-based central signal processing board ("Reconfigurable Open Architecture Computing Hardware'') that was developed in a collaboration among NRAO, the South African KAT project, and the Collaboration for Astronomy Signal Processing and Electronics Research (CASPER) at UC Berkeley. In addition to the ROACH, the RDBE includes an input analog level control module, a sampler developed by CASPER, and a synthesizer board which generates the 1024-MHz sample clock. RDBEs accept two 512-1024 MHz IF inputs, and deliver packetized output via a 10G Ethernet interface. Each VLBA station is equipped with two RDBE units.
Currently, two separate "observing systems" are available within the VLBA's RDBE. Inputs to either data system can come from any of the four VLBA IFs. Some suggestions for choosing between the observing systems follow the functional outlines below.
PFB: The polyphase filterbank digital signal-processing algorithm produces sixteen fixed-bandwidth 32-MHz channels within a single RDBE unit. Channels can be selected flexibly between two input IFs, and can be placed at 32-MHz steps along the entire IF frequency range. Some typical selection modes include [a] a compact dual-polarization configuration of eight contiguous 32-MHz channels at matching frequencies in each polarization; [b] a spanned-band dual-polarization configuration, with eight 32-MHz channel pairs spaced every 64 MHz; and [c] a single-polarization configuration of 16 channels, contiguous across the entire width of one IF. (In case [c], one end channel will not lie within the IF band, and does not produce usable data.) The selected channels are requantized at two bits per Nyquist sample and output in a packetized stream at a total data rate of 2048 Mbps (referred to subsequently as "2 Gbps'').
DDC: The digital downconverter algorithm supports a wide range of bandwidths. A total of 1, 2, or 4 channels can be processed within a single RDBE unit; 4 or 8 channels are available using both RDBEs. Available bandwidths range downward from 128 MHz to 1 MHz in binary steps; recording rate limitations restrict the 128-MHz bandwidth to a maximum of 4 channels. All channels must use the same bandwidth within an observing scan. Channels can be selected flexibly among up to four input IFs, and in either sideband. Tuning of individual channels can be set in steps of 15.625 kHz, although 250-kHz steps are recommended when compatibility with legacy systems is required. Channels may not cross IF zone boundaries at 640 and 896 MHz. Each channel is requantized at two bits per Nyquist sample and output in a packetized stream, at a total data rate ranging from 4 to 2048 Mbps (subsequently "2 Gbps").
Suggestions for Observing System Selection: Wideband science will be possible using either the PFB observing system, at its fixed 2048 Mbps data rate, or the DDC system at 2048 Mbps or lower rates. Both systems provide output at two bits per Nyquist sample. The primary instrumental differences are in the numbers and bandwidths of channels, and in the channel passbands. The PFB's many narrower channels may be advantageous in avoiding spectral ranges impacted by interference, particularly in the 18-cm band. On the other hand, the smaller number of wider-band channels available in the DDC may simplify data analysis in some cases. Digital logic capacity of the RDBE limits the PFB's signal processing to fewer filter taps for each of its 16 channels than for the 4-channel DDC system, so that the DDC's passbands cut off significantly more sharply.
Spectroscopic and other narrow-band observations will generally be best supported by the DDC system, which incorporates scientifically equivalent counterparts for all modes of the VLBA legacy system, and extends these to wider bandwidths. Even extremely narrow bands can be accommodated by observing at 1 MHz bandwidth and selecting a narrower range using the DiFX correlator's spectral zoom mode.
Most VLBA receivers produce only two IFs, in opposite polarizations, but some receivers support four-IF modes, such as dual-polarization dual-frequency. The four-IF capability of the DDC allows these modes to be exploited.
Conversion of Legacy Schedules to RDBE/DDC: A separate web page describes the relatively straightforward conversion of SCHED “keyin” files applicable to the VLBA's legacy data system, to use the DDC system instead. Explanations of the "scientifically equivalent" modes referred to above are included. It is designed primarily for users with some VLBA experience.
5.5. Programmable Network Switch
A software-based network switch, purchased from XCube Research and Development, is an essential element of the data path at each VLBA station. Its primary functions are to merge the packet streams from each of the two RDBE units into a single stream that is sent to the Mark 5C recorder, and to regulate the timing of these packets so as not to overflow the Mark 5C's input buffer.
Other switching and real-time data analysis functions may be added as part of future developments. A phase-cal detection capability is currently under test.
5.6. Mark 5C Recorder
The VLBA's data transmission system comprises the recorder units at the stations, playback units at the correlator, and the magnetic disk modules that are shipped between those units. The Mark 5C recording system was developed jointly by NRAO, Haystack Observatory, and Conduant Corporation. It closely resembles the Mark 5A version used previously by the VLBA, and the Mark 5B used at some other observatories. In particular, identical disk modules are used. However, Mark 5C is a packet-based system, which allows a more straightforward functionality than its predecessors. It simply records the payload of each 10G Ethernet packet originating in the RDBE, without imposing any special recording format. All formatting of the observed data — most essentially, the precision time tag — is internal to the packet payloads, which are transmitted directly from recorder to playback by the Mark 5C system. Mark 5B formatting is used internally by the RDBE's PFB personality, but the newer VLBI Data Interchange Format (VDIF) is used in the DDC personality.
Each Mark 5C unit accommodates two removable modules, each in turn comprising eight commercial disk drives. As used on the VLBA, these modules are recorded sequentially at a maximum rate of 2 Gbps, matching the current maximum RDBE output rate. Available recording media are expected to support the highest data rates for approximately half of all observing hours.
Further information on the Mark 5C system is available in the Sensitivity Upgrade memo series, and in the Haystack Mark 5 series.
6. Correlator
6.1. Introduction
The correlator is situated in the SOC, at the end of the data path. Its role is to reproduce the signals recorded at the VLBA stations and any others involved in the observation, and to combine them in two-station baseline pairs, to yield the visibility function which is the fundamental measurement produced by the VLBA. VLBA observations are processed using the DiFX software correlator. DiFX was developed at Swinburne University in Melbourne, Australia (Deller et al. 2007), and adapted to the VLBA operational environment by NRAO staff (Brisken 2008). Subsequent references to "DiFX" apply specifically only to this VLBA implementation.
We encourage users to include the following text in the Acknowledgments section of any publication arising from VLBA observations made since December 2009:
This work made use of the DiFX software correlator developed at Swinburne University of Technology as part of the Australian Major National Research Facilities program.
... and to cite the following paper by the developers: Deller, et al. 2011, PASP, 123, 275.
Software correlation has become feasible in recent years, and is especially well suited to applications like VLBI with bandwidth-limited data-transmission systems and non-realtime processing. Among its several advantageous aspects are: (1) flexible allocation of processing resources to support correlation of varying numbers of stations, frequency and time resolution, and various special processing modes, with no fundamental fixed limits other than the finite performance of the processing cluster; (2) optimization of resource usage to minimize processing time; (3) integration of control and processing functions; (4) continuously scalable, incremental upgrade paths; and (5) relatively straightforward implementation of special modes and tests. These and other virtues of software correlation are discussed in more detail by Deller et al. (2007).
Despite the absence of fixed limits cited in item (1) above, guidelines have been established for the extremes of spectral resolution, integration period, and output rate, for routine DiFX processing, as specified in the appropriate sections below. Exceptions will be considered for proposals including a sufficiently compelling scientific justification.
The VLBA DiFX correlator is not configured to process data from a single antenna, nor is a multi-station autocorrelation-only mode available.
Operation of DiFX is governed primarily by an observation description in VEX format (currently vex1.5). This format is used for both station and correlator control functions in a number of VLBI arrays, and VLBA program SCHED (Walker 2011) has been producing it for many years.
DiFX only accepts data on Mark5 disk modules, as recorded by a Mark5A, Mark5B, Mark5B+ or Mark5C recorder. It can process data in a variety of formats including VLBA, Mark4, and Mark5B. Support for VDIF format is currently incomplete but includes those versions created by the VLBA RDBE and the VLA WIDAR correlator. Modes recorded at EVN stations are also thought to be fully supported.
Correlator output is written according to the FITS Interferometry Data Interchange Convention (Greisen 2009). In addition to the fundamental visibility function measurements and associated meta-data, the FITS files include amplitude and phase calibration measurements, weather data, and editing flags, all derived from data logged at the observing stations. An up-to-date release of AIPS is required to handle DiFX data properly.
Conversion of DiFX correlator output to the Mark 4 format that is used primarily in analysis of geodetic observations is also available. To enable this additional output, a SCHED parameter CORDFMT=MARK4 should be specified.
6.2. Spectral Resolution
DiFX allows quite flexible selection of the desired number of "spectral points" spanning each individual data channel. Any number that can be factored as 2n · 5m can be specified, subject to these limitation:
- A maximum of 4096 points per channel, for routine DiFX processing.
- A total of 132,096, summed over all channels and polarization products, for compatibility with AIPS.
- A minimum spectral resolution of 2 Hz.
The number of spectral points must be the same for all data channels at any given time, although multiple passes are possible with different sets of channels. The actual spectral resolution obtained, and statistical independence of the spectral points, depends on subsequent smoothing and other processing.
DiFX also supports "spectral zooming'', selection of a subset of correlated spectral points from any or all data channels. Only the selected spectral points are included in the output dataset. This capability is of value mainly in maser studies, where a recorded data channel may be much wider than the maser emission in two main categories of observations: (1) Maser astrometry with in-beam continuum calibrators. Wideband observing is required for maximum sensitivity on the calibrators, while zooming allows high spectral resolution at the frequencies where maser emission appears. (2) Multiple maser transitions. When wideband data channels are used to cover a large number of widely separated maser transitions, spectral zooming allows the empty portions of high-resolution spectrum to be discarded.
Spectral zooming does not work with mixed sideband observing, this can happen for HSA or global observations where some telescopes require upper sideband and others require lower sideband. In proposing observations that will use spectral zooming, the required number of spectral points before zooming should be specified in the Proposal Submission Tool. Currently, the location and width of the "zoom" bands must be communicated directly to VLBA operations before correlation.
6.3. Integration Period
DiFX accommodates a nearly continuous range of correlator integration periods over the range of practical interest. Individual integrations are quantized in multiples of the indivisible internal FFT interval, which is equal to the number of spectral points requested, divided by the data channel bandwidth.
For most cases, with low to moderate spectral resolution, and/or wide data channels, the FFT intervals are fairly short, and it is straightforward to find an integration period in any desired range that is an optimal integral multiple of the FFT interval. ("Optimal'' refers here to the performance of DiFX.) Extreme cases of very high spectral resolution (many spectral points across a narrow data channel - resolution of less than about 100 Hz) imply FFT intervals long enough that only limited choices of integral multiples are available.
For flexibility in these situations (although the option exists in all cases), integration periods other than an integral multiple of the FFT interval can be approximated, in a long-term mean, by an appropriate sequence of nearby optimal integral multiples. In this case, output records are time-tagged as if correlated with exactly the requested period.
SCHED accepts an additional parameter so that users can indicate that the requested integration period is to be implemented exactly, as described above. Otherwise, the nearest optimal integral multiple of the FFT interval is passed to the correlator.
6.4. Pulsar Gating and Binning
DiFX supports the following three time-based selection modes to facilitate pulsar observations. In all cases, a pulsar spin ephemeris must be provided by the user. Except for certain applications of mode 3, the ephemeris must be capable of predicting the absolute rotation phase of the pulsar. Pulsar modes incur a minimum correlation-time penalty of about 50%. High output data rates may require greater correlator resource allocations.
- Binary Gating: A simple pulse-phase driven on-off accumulation window can be specified, with "on" and "off" phases. Such gating increases the signal to noise ratio of pulsar observations by a factor of typically 3 to 6, and can also be used to search for off-pulse emission.
- Matched-filter Gating: If the pulse profile at the observation frequency is well understood and the pulse phase is very well predicted by the provided pulse ephemeris, additional signal to noise over binary gating can be attained by appropriately scaling the correlation coefficients as a function of pulse phase. Depending on the pulse shape, additional gains of up to 50% in sensitivity over binary gating can be realized.
- Pulsar Binning: This mode entails generating a separate visibility spectrum for each requested range of pulse phase. There are no explicit limits to the number of pulse phase bins that are supported, however, data rates can become increasingly large. Currently AIPS does not support databases with multiple phase bins. Until post-processing support is available, a separate FITS file will be produced for each pulsar phase bin.
Details of pulsar observing, including practical aspects of using the pulsar modes, and limitations imposed by operations, are documented by Brisken & Deller (2010).
6.5. Multiple Phase Centers
The field of view in VLBI observations is very small, around 10-4 of the primary antenna beam area. This restricted interferometer beam arises in the correlation process from smearing at positions away from the correlation phase center, due to averaging in time (with, typically, a 2-second period) and/or across bandwidth ("chromatic aberration'' over, typically, 0.5 MHz spectral resolution). Thus, imaging of targets that are widely spaced in the primary beam requires multiple processing passes in typical correlator implementations. If the visibilities are maintained at high time and frequency resolution, it is possible to perform a u-v shift after correlation, essentially repointing the correlated dataset to a new phase center. However, this approach would require prohibitively large visibility datasets.
DiFX implements multiple u-v shifts inside the correlator, to generate as many phase centers as are necessary, in a single correlation pass. The output consists of one dataset of normal size for each phase center. This mode consumes around three times the correlator resources of a normal continuum correlation, due to the need for finer frequency resolution before the u-v shift, but the additional cost is only weakly dependent on the number of phase centers. For reasonable spectral and temporal resolution requirements (for example, adequate for smearing < 10% at the 50% contour of the VLBA primary beam), 200 phase centers require only 20% more correlator time than 2 phase centers. Extremely high spectral and/or temporal resolution (e.g. for shifts even closer to the edge of the primary beam) carry a higher overhead per additional phase center. This mode thus should be requested only for imaging of three or more sources within any single antenna pointing. The correlator output rate expands proportionally to the number of phase centers.
Correlator memory limits the product of baselines, spectral points, and phase centers for one correlator pass. The current limit is approximately 600 phase centers for the 10 element VLBA at 2 Gbps record rate (512 MHz polarization-summed bandwidth). An unlimited number of phase centers can ultimately be achieved in multiple correlation passes, however.
Multiple phase-center correlation is requested in the NRAO Proposal Submission Tool by setting the "Number of Fields'' item in the resource section to the maximum number of phase centers required for any antenna pointing specified in a given resource. The requested spectral resolution and integration time should correspond to the desired initial number of spectral points per data channel (required to minimize bandwidth smearing) and the desired integration between u-v shifts (to minimize time smearing). A resulting expanded output data rate that exceeds the current limit, as well as any required multiple passes, must be justified specifically in the proposal.
SCHED includes facilities to support specification of the actual phase center locations to be used in correlation.
For more details on wide-field imaging techniques, see Bridle & Schwab (1999), and Garrett et al. (1999).
6.6. Output Rate
Correlation parameters should result in an output rate less than 10 MBytes per second (of observing time) for routine DiFX processing; higher rates may be considered if required and adequately justified. Observers should ensure that their data-analysis facilities can handle the dataset volumes that will result from the correlation parameters they specify.
An approximate parametrization of the output rate is given by
[display]R = 4 \cdot \frac{N_{\rm stn} \cdot (N_{\rm stn}+1) \cdot N_{\rm sbb} \cdot N_{\rm spc}}{T_{\rm int}}\cdot N_{\rm ppb} \cdot N_{\rm phc} \cdot p[/display]
where the rate [inline]R[/inline] is in Byte/s;
[inline] N_{\rm stn} , \; N_{\rm sbb} , \; N_{\rm spc} [/inline] are the numbers of observing stations, data channels, and spectral points per data channel, respectively;
[inline]T_{\rm int}[/inline] is the correlator integration period;
[inline]N_{\rm ppb}[/inline] is the number of pulsar phase bins; and
[inline]N_{\rm phc}[/inline] is the number of phase centers.
The polarization factor [inline]p=1[/inline] for single-polar, or dual-polar parallel-hand output; or [inline]p=2[/inline] for cross-polar, four-Stokes processing.
Output data rates are also estimated by SCHED.
7. Angular Resolution & u-v Coverage
Table 5 gives the maximum lengths ([inline]B^{\rm km}_{\rm max}[/inline]) for each of the VLBA's 45 internal baselines as well as the baselines to HSA telescopes. A measure of the corresponding resolution ([inline]\theta_{\rm HPBW}[/inline]) in milliarcseconds (mas) is
[display]\theta_{\rm HPBW} \sim 2063 \times \frac{\lambda^{\rm cm}}{B^{\rm km}_{\rm max}} \; {\rm mas}[/display]
where [inline]\lambda^{\rm cm}[/inline] is the receiver wavelength in cm (Wrobel 1995). A uniformly weighted image made from a long u-v plane track will have a synthesized beam with a slightly narrower minor axis.
SC | HN | NL | FD | LA | PT | KP | OV | BR | MK | EB | AR | GB | Y27 | |
SC | ... | 2853 | 3645 | 4143 | 4458 | 4579 | 4839 | 5460 | 5767 | 8611 | 6822 | 238 | 2708 | 4532 |
HN | 2853 | ... | 1611 | 3105 | 3006 | 3226 | 3623 | 3885 | 3657 | 7502 | 5602 | 2748 | 829 | 3198 |
NL | 3645 | 1611 | ... | 1654 | 1432 | 1663 | 2075 | 2328 | 2300 | 6156 | 6734 | 3461 | 1064 | 1640 |
FD | 4143 | 3105 | 1654 | ... | 608 | 564 | 744 | 1508 | 2345 | 5134 | 8084 | 3922 | 2354 | 515 |
LA | 4458 | 3006 | 1432 | 608 | ... | 236 | 652 | 1088 | 1757 | 4970 | 7831 | 4246 | 2344 | 226 |
PT | 4579 | 3226 | 1663 | 564 | 236 | ... | 417 | 973 | 1806 | 4795 | 8014 | 4365 | 2551 | 52 |
KP | 4839 | 3623 | 2075 | 744 | 652 | 417 | ... | 845 | 1913 | 4466 | 8321 | 4623 | 2939 | 441 |
OV | 5460 | 3885 | 2328 | 1508 | 1088 | 973 | 845 | ... | 1214 | 4015 | 8203 | 5255 | 3323 | 1025 |
BR | 5767 | 3657 | 2300 | 2345 | 1757 | 1806 | 1913 | 1214 | ... | 4398 | 7441 | 5585 | 3326 | 1849 |
MK | 8611 | 7502 | 6156 | 5134 | 4970 | 4795 | 4466 | 4015 | 4398 | ... | 10328 | 8434 | 7028 | 4835 |
EB | 6822 | 5602 | 6734 | 8084 | 7831 | 8014 | 8321 | 8203 | 7441 | 10328 | ... | 6911 | 6335 | 8008 |
AR | 238 | 2748 | 3461 | 3922 | 4246 | 4365 | 4623 | 5255 | 5585 | 8434 | 6911 | ... | 2545 | 4317 |
GB | 2708 | 829 | 1064 | 2354 | 2344 | 2551 | 2939 | 3323 | 3326 | 7028 | 6335 | 2545 | ... | 2516 |
Y27 | 4532 | 3198 | 1640 | 515 | 226 | 52 | 441 | 1025 | 1849 | 4835 | 8008 | 4317 | 2516 | ... |
Values of [inline]\theta_{\rm HPBW}[/inline] for the longest VLBA baseline, at the center frequencies of the VLBA's observing bands (Table 3), are shown in Table 6. The longest VLBA baseline at 3 mm is currently that between MK and NL, which is about 30% shorter than the longest baseline at lower frequencies.
Customized plots of the u-v plane coverage with the VLBA and/or other VLBI stations can be generated by VLBA program SCHED (Walker 2011).
8. Baseline Sensitivity
Baseline sensitivity is the RMS thermal noise (ΔS) in the visibility amplitude in a single polarization on a single baseline. Adequate baseline sensitivity is required for VLBI fringe fitting. Baseline sensitivities between VLBA antennas, for typical observing parameters, are listed in column [6] of the Receiver Frequency Ranges & Performance table.
Alternatively, the baseline sensitivity for two identical antennas, in the weak source limit, can be calculated using the formula (Walker 1995a; Wrobel & Walker 1999):
[display]\Delta S = {\rm SEFD} / [\eta_s \cdot ( 2 \cdot \Delta \nu \cdot \tau_{\rm ff} ) ^{1/2} ] \; {\rm Jy}[/display]
SEFD or "system equivalent flux density" is the system noise expressed in Janskys. [inline]\eta_s \le 1 \ \ [/inline] accounts for the VLBI system inefficiency (primarily quantization in the data recording). Kogan (1995b) provides the combination of scaling factors and inefficiencies appropriate for VLBA visibility data. The bandwidth in Hz is [inline]\Delta\nu[/inline]. For a continuum target, use the channel bandwidth or the full recorded bandwidth, depending on the fringe-fitting mode; for a line target, use the channel bandwidth divided by the number of spectral points that span the channel. [inline]\tau_{\rm ff}[/inline] is the fringe-fit interval in seconds, which should be less than or about equal to the coherence time.
Moran & Dhawan (1995) discuss expected coherence times. The actual coherence time appropriate for a given observation can be estimated using observed fringe amplitude data on an appropriately strong and compact source.
For non-identical antennas 1 and 2, SEFD can be replaced by the geometric mean [inline]\sqrt{{\rm SEFD}_1 \times {\rm SEFD}_2}[/inline].
9. Image Sensitivity
Image sensitivity is the RMS thermal noise (ΔIm) expected in a single-polarization image. Image sensitivities for the 10-station VLBA, for typical observing parameters, are listed in column [7] of the Receiver Frequency Ranges & Performance table.
Alternatively, the image sensitivity for a homogeneous array with natural weighting can be calculated using the following formula (Wrobel 1995; Wrobel & Walker 1999).
[display]\Delta I_m = {\rm SEFD} / [\eta_s \cdot ( N \cdot (N-1) \cdot \Delta \nu \cdot t_{\rm int} ) ^{1/2} ] \; \rm{Jy\; beam^{-1}}[/display]
Parameters SEFD, [inline]\eta_s[/inline], and [inline]\Delta\nu[/inline] are the same as those used in computing baseline sensitivity, [inline] N [/inline] is the number of observing stations, and [inline]t_{\rm int}[/inline] is the total integration time on source in seconds.
The expression for image noise becomes rather more complicated for a heterogeneous array such as the HSA, and may depend quite strongly on the data weighting that is chosen in imaging. The EVN sensitivity calculator provides a convenient estimate. For example, the RMS noise at 22 GHz for the 10-station VLBA in a 1-hr integration is reduced by a factor between 4 and 5 by adding the GBT and the phased VLA.
If simultaneous dual polarization data are available with the above value of ΔIm per polarization, then for an image of Stokes I, Q, U, or V,
[display]\Delta I = \Delta Q = \Delta U = \Delta V = \frac{\Delta I_m}{\sqrt{2}}[/display]
For a polarized intensity image of [inline]P = \sqrt{Q^2 + U^2}[/inline]
[display]\Delta P = 0.655 \times \Delta Q = 0.655 \times \Delta U[/display]
It is sometimes useful to express [inline]\Delta I_m[/inline] in terms of an RMS brightness temperature in Kelvins ([inline]\Delta T_B[/inline]) measured within the synthesized beam. An approximate formula for a single-polarization image is
[display]\Delta T_b \sim 320 \times \Delta I_m \times (B^{\rm km}_{\rm max})^2 \; {\rm K}[/display]
where [inline]B^{\rm km}_{\rm max}[/inline] is as in Table 5.
10. Calibration Data
Data necessary to perform accurate calibration for the VLBA are supplied as part of the correlator output files, and will appear as extension tables within the AIPS datasets created by task FITLD. These tables include GC (gain), TY (system temperature), and WX (weather) tables for amplitude calibration, PC (pulse-cal) tables for system phase calibration, and FG (flag) tables for editing.
For non-VLBA stations, some or all of these tables may be missing, since relevant measurements are not available at the time of correlation. For example, for the HSA, GC and TY information are available for most stations, except that calibration of the phased VLA requires additional information about the flux density of at least one source. Flag (FG) tables for non-VLBA stations generally are absent or only partially complete, lacking information about antenna off-source times. However, the "flag" file that is written by program SCHED (Walker 2011) is quite good at predicting the on-source times for the HSA stations. In using this file as an input to AIPS task UVFLG, it is recommended that all entries for the ten VLBA stations be deleted. The FG table supplied with the correlator output files includes the actual on-source times for these antennas, obtained directly from VLBA monitor data. For further information on applying calibrations, see Appendix C of the AIPS Cookbook (NRAO staff, 2006) or the relevant AIPS HELP files.
11. Amplitude Calibration
Traditional calibration of VLBI fringe amplitudes for continuum sources requires knowing the on-source system noise in Jy (SEFD; Moran & Dhawan 1995). System temperatures in Kelvin ([inline]T_{\rm sys}[/inline]) are measured continuously during observations at VLBA stations, with mean values tabulated at least once per source/frequency combination or once every user-specified interval (default 2 minutes), whichever is shorter. These [inline]T_{\rm sys}[/inline] values are used in fringe amplitude calibration by AIPS task APCAL, which converts [inline]T_{\rm sys}[/inline] to SEFD by dividing by the VLBA antenna gains in [inline]\rm K\: Jy^{-1}[/inline], expressed as a peak gain multiplied by a normalized "gain curve". The latter data are based on regular monitoring of all receiver and feed combinations. [inline]T_{\rm sys}[/inline] and gain values for VLBA antennas are delivered in TY and GC tables, respectively. Single-station spectra can be used for amplitude calibration of spectral line observations.
Additional amplitude adjustments may be necessary to correct for the atmospheric opacity above an antenna, which can be significant at high frequencies (Moran & Dhawan 1995). Leppänen (1993) describes a method for opacity adjustments. AIPS task APCAL uses weather data from the WX table and the system temperature data to carry out such adjustments.
Further corrections are usually applied to observations taken with 2-bit (4-level) sampling, for the effects of non-optimal setting of the quantizer voltage thresholds (Kogan 1995a). These adjustments are usually relatively minor but can induce systematic effects. Sampling-based calibration adjustments are determined by AIPS task ACCOR. The combination of the antenna and quantizer calibrations may be found and applied in AIPS using the procedure VLBACALA.
Although experience with VLBA calibration shows that it probably yields fringe amplitudes accurate to 5% or less at the standard frequencies in the 1-10 GHz range, it is recommended that users observe a few amplitude calibration check sources during their VLBA program. Such sources can be used (1) to assess the relative gains of VLBA antennas plus gain differences among channels at each station; (2) to test for non-closing amplitude and phase errors; and (3) to check the correlation coefficient adjustments, provided contemporaneous source flux densities are available independent of the VLBA observations. These calibrations are particularly important if non-VLBA stations are included in an observation, since their a priori gains and/or measured system temperatures may be much less accurate than for the well-monitored VLBA stations. The recommended technique for this situation is to restrict the gain normalization in self-calibration to a subset of trusted stations (generally some of the VLBA stations), and to high elevations. AIPS task CALIB can do both.
The VLBA gains are measured at the center frequencies appearing in column [4] of the Receiver Frequency Ranges & Performance table; users observing at other frequencies may be able to improve their amplitude calibration by including brief observations, usually of their amplitude check sources, at the appropriate frequencies. Amplitude check sources should be point-like on inner VLBA baselines. Some popular choices in the range 13 cm to 2 cm are J0555+3948=DA193, J0854+2006=OJ287, and J1310+3220. Other check sources may be selected from various VLBI surveys. It might be prudent to avoid sources known to have exhibited extreme scattering events (e.g., Fiedler et al. 1994a, b).
12. Phase Calibration & Imaging
12.1. Fringe Finders
VLBI fringe phases are much more difficult to deal with than fringe amplitudes. If the a priori correlator model assumed for VLBI correlation is particularly poor, then the fringe phase can wind so rapidly in both time (the fringe rate) and in frequency (the delay) that no fringes will be found within the finite fringe rate and delay windows examined during correlation. Reasons for a poor a priori correlator model include source position and station location errors, atmospheric (tropospheric and ionospheric) propagation effects, and the behavior of the independent clocks at each station. Users observing sources with poorly known positions should plan to refine the positions first on another instrument. To allow accurate location of any previously unknown antennas and to allow VLBA staff to conduct periodic monitoring of clock drifts, each user should include one or more "fringe finder" sources which are strong, compact, and have accurately known positions. Consult Markowitz & Wurnig (1998) to select a fringe finder for observations between between 20 cm and 7 mm; your choice will depend on your wavelengths but J0555+3948=DA193, J0927+3902=4C39.25, J1642+3948=3C345, and J2253+1608=3C454.3 are generally reliable in the range 13 cm to 2 cm. In addition, at 90 and 50 cm we recommend either J1331+3030=3C286 or J2253+1608=3C454.3. Fringe-finder positions, used by default by VLBA program SCHED (Walker 2011) and the VLBA correlator, are given in the standard source catalog available as an ancillary file with SCHED.
12.2. The Pulse Cal System
Fringe phases should be coherent across the entire set of channels produced by each RDBE. Correction of phase offsets between the two RDBEs at each station, and/or between the oppositely polarized signal channels, can be determined using the "phase cal'' or "pulse cal'' system (Thompson 1995). In conjunction with the LO cable length measuring system, this system can also be used to measure changes in the delays through the cables and electronics which must be removed for accurate geodetic and astrometric observations.
The phase cal system consists of a pulse generator and a sine-wave detector. The interval between the pulses can be either 0.2 or 1 microsecond. They are injected into the signal path at the receivers and serve to define the delay reference point for astrometry. The pulses appear in the spectrum as a "comb'' of very narrow, weak spectral lines at integral multiples of 1 or 5 MHz. The phases of one or more of these lines are measured by the detector, logged as a function of time, and delivered in a PC table.
AIPS tasks can load and apply the PC data. However, some VLBA observers may still want to use a strong compact source to do a "manual'' phase cal if necessary (Diamond 1995). Spectral line users will not want the pulse cal comb to appear in their observations, and should ensure that their observing schedules both disable the pulse cal generators and include observations suitable for a manual phase cal. Manual phase calibration also is likely to be necessary for non-VLBA stations that have no tone generators or detectors, and in VLBA observations at 3 mm, where the VLBA receivers have no pulse calibration tones.
12.3. Fringe Fitting
After correlation and application of the pulse calibration, the phases on a VLBA target source still can exhibit high residual fringe rates and delays. Before imaging, these residuals should be removed to permit data averaging in time and, for a continuum source, in frequency. The process of finding these residuals is referred to as fringe fitting. Before fringe fitting, it is recommended to edit the data based on the a priori edit information provided for VLBA stations. Such editing data are delivered in the FG table. The old baseline-based fringe search methods have been replaced by more powerful global fringe search techniques (Cotton 1995a; Diamond 1995). Global fringe fitting is simply a generalization of the phase self-calibration technique, as during a global fringe fit the difference between model phases and measured phases are minimized by solving for the station-based instrumental phase, its time slope (the fringe rate), and its frequency slope (the delay). Global fringe fitting in AIPS is done with the program FRING or associated procedures. If the VLBA target source is a spectral line source or is too weak to fringe fit on itself, then residual fringe rates and delays can be found on an adjacent strong continuum source and applied to the VLBA target source in a phase-referencing technique.
VLBA delays do tend to be very stable (1 to a few ns) during observations at higher frequencies where the ionospheric variations are limited. Thus one, or a small number, of high SNR fringe fits on strong fringe finders may provide superior results for delay, over trying to fit weaker sources. The phases will still need to be corrected, but that can be done with self-calibration or a phase-only fringe fit.
12.4. Editing
After fringe-fitting and averaging, VLBA visibility amplitudes should be inspected and obviously discrepant points removed (Diamond 1995; Walker 1995b). Usually such editing is done interactively using tasks in AIPS or the Caltech program Difmap (Shepherd 1997). VLBA correlator output data also includes flags derived from monitor data output in an FG table, containing information such as off-source flags for the stations during slews to another source.
12.5. Self-Calibration, Imaging, and Deconvolution
Even after global fringe fitting, averaging, and editing, the phases on a VLBA target source can still vary rapidly with time. Most of these variations are due to inadequate removal of station-based atmospheric phases, but some variations also can be caused by an inadequate model of the source structure during fringe fitting. If the VLBA target source is sufficiently strong and if absolute positional information is not needed, then it is possible to reduce these phase fluctuations by looping through cycles of Fourier transform imaging and deconvolution, combined with phase self-calibration in a time interval shorter than that used for the fringe fit (Cornwell 1995; Walker 1995b; Cornwell & Fomalont 1999). Fourier transform imaging is straightforward (Briggs, Schwab, & Sramek 1999), and done with AIPS task IMAGR or the Caltech program Difmap (Shepherd 1997). The resulting VLBI images are deconvolved to rid them of substantial sidelobes arising from relatively sparse sampling of the u-v plane (Cornwell, Braun, & Briggs 1999). Such deconvolution is achieved with AIPS tasks based on the CLEAN or Maximum Entropy methods or with the Caltech program Difmap.
Phase self-calibration just involves minimizing the difference between observed phases and model phases based on a trial image, by solving for station-based instrumental phases (Cornwell 1995; Walker 1995b; Cornwell & Fomalont 1999). After removal of these instrumental phases, the improved visibilities are used to generate an improved set of model phases, usually based on a new deconvolved trial image. This process is iterated several times until the phase variations are substantially reduced. The method is then generalized to allow estimation and removal of complex instrumental antenna gains, leading to further image improvement. Both phase and complex self-calibration can be accomplished using AIPS task CALIB or with the Caltech program Difmap. Self-calibration should only be done if the VLBA target source is detected with sufficient signal-to-noise in the self-calibration time interval (otherwise, fake sources can be generated!), and if absolute positional information is not needed.
The useful field of view in VLBI images can be limited by finite bandwidth, integration time, and non-coplanar baselines (Wrobel 1995; Cotton 1999b; Bridle & Schwab 1999; Perley 1999b). Measures of image correctness -- image fidelity and dynamic range -- are discussed by Walker (1995a) and Perley (1999a).
12.6. Phase Referencing
If the VLBA target source is not sufficiently strong for self-calibration or if absolute positional information is needed but geodetic techniques are not used, then VLBA phase referenced observations must be employed (Beasley & Conway 1995). Currently, 63% of all VLBA observations employ phase referencing. Wrobel et al. (2000) recommend strategies for phase referencing with the VLBA, covering the proposal, observation, and correlation stages. A VLBA phase reference source should be observed frequently and be within a few degrees of the VLBA target region, otherwise differential atmospheric (tropospheric and ionospheric) propagation effects will prevent accurate phase transfer. VLBA users can draw candidate phase calibrators from the source catalog distributed with SCHED, which contains over 7000 sources. Easy searching for the nearest calibrators is available online through the VLBA Calibrator Survey (Beasley et al. 2002). Most of these candidate phase calibrators now have positional uncertainties below 1 mas.
Calibration of atmospheric effects for either imaging or astrometric observations can be improved by the use of multiple phase calibrators that enable multi-parameter solutions for phase effects in the atmosphere. See AIPS Memos 110 (task DELZN, Mioduszewski 2004) and 111 (task ATMCA, Fomalont & Kogan 2005), available from the AIPS home page, for further information.
Walker & Chatterjee (1999) have investigated ionospheric corrections using GPS based ionospheric models. Such corrections can be of significant benefit for even the highest frequencies on the VLBA. These corrections may be made with the AIPS task TECOR, as described in AIPS Cookbook Appendix C (NRAO 2006), or the procedure VLBATECR. In addition, it is strongly recommended that the most accurate Earth-Orientation values be applied to the calibration, since correlation may have taken place before final values were available; this may be done with AIPS task CLCOR or more easily with the AIPS procedure VLBAEOPS.
The rapid motion of VLBA antennas often can lead to very short time intervals for the slew between target source and phase reference source, and some data may be associated with the wrong source. Users should be alert for visibility points of very low amplitude at the beginnings of scans.
13. Specialized Observational Techniques
13.1. Polarimetry
In VLBA polarimetric observations, channels are assigned in pairs to opposite hands of circular polarization at each frequency. Typical "impurities" of the antenna feeds are about 3% for the center of most VLBA bands and degrade toward the band edges and away from the pointing center in the image plane. Without any polarization calibration, an unpolarized source will appear to be polarized at the 2% level. Furthermore, without calibration of the RCP-LCP phase difference, the polarization angle is undetermined. With a modest investment of time spent on calibrators and some increased effort in the calibration process, the instrumental polarization can be reduced to less than 0.5%.
To permit calibration of the feed impurities (sometime also called "leakage" or "D-terms"), VLBA users should include observations of a strong (≈ 1 Jy) calibration source, preferably one with little structure. This source should be observed during at least 5 scans covering a wide range (> 100 degrees) of parallactic angle, with each scan lasting for several minutes. The electric vector polarization angle (EVPA) of the calibrator will appear to rotate in the sky with parallactic angle while the instrumental contribution stays constant. Some popular calibrator choices are J0555+3948=DA193 and J1407+2827=OQ208, although either or both may be inappropriate for a given frequency or an assigned observing time. Fortunately, many calibrators satisfying the above criteria are available.
A viable alternative approach to measuring polarization leakage is to use an unpolarized calibrator source. This can be done with a single scan.
The wide channel bandwidths now available may make it necessary to apply frequency-dependent instrumental polarization corrections. Most VLBI calibrators are resolved, and the usual AIPS tool for solving for instrumental polarization using resolved sources, LPCAL, does not handle this frequency dependence. Procedures are being tested for dealing with this issue. They are likely to involve making polarization images based on the best LPCAL results, using them to divide the data, then using the frequency dependence capability in PCAL to do the rest.
To set the absolute EVPA on the sky, it is necessary to determine the phase difference between RCP and LCP. For VLBA users at frequencies of 5 GHz and above, the best method for EVPA calibration is to observe one or two of the compact sources that are being monitored with the VLA; see the VLA/VLBA Polarization Calibration Page (Taylor & Myers 2000). At 1.6 GHz it may be preferable to observe a source with a stable, long-lived jet component with known polarization properties. At frequencies of 5 GHz and below one can use J0521+1638=3C138 (Cotton et al. 1997a), J1331+3030=3C286 (Cotton et al. 1997b), J1829+4844=3C380 (Taylor 1998), or J1902+3159=3C395 (Taylor 2000). At 8 GHz and above one may use J1256-0547=3C279 (Taylor 1998) or J2136+0041=2134+004 (Taylor 2000), although beware that some of these jet components do change on timescales of months to years. It will be necessary to image the EVPA calibrator in Stokes I, Q, and U, and to determine the appropriate correction to apply. Thus it is recommended to obtain 2 to 4 scans, each scan lasting at least 3 minutes, over as wide a range in hour angle as is practical.
To permit calibration of the RCP-LCP delays, VLBA users should include a 2-minute observation of a very strong (≈ 10 Jy) calibration source. While 3C279 is a good choice for this delay calibration, any very strong fringe-finder will suffice.
Post-processing steps include amplitude calibration; fringe-fitting; solving for the RCP-LCP delay; self-calibration and Stokes I image formation; instrumental polarization calibration; setting the absolute position angle of electric vectors on the sky; and correction for ionospheric Faraday rotation, if necessary (Cotton 1995b, 1999a; Kemball 1999). All these post-processing steps can currently be done in AIPS, as can the polarization self-calibration technique described by Leppänen, Zensus, & Diamond (1995).
13.2. Spectral Line Observations
Diamond (1995) and Reid (1995, 1999) describe the special requirements for data acquisition, correlation, and post-processing of spectroscopic VLBI observations. The transition rest frequency, approximate velocity, and velocity width for the line target must be known in order to set the observing frequency and bandwidth correctly. The schedule should include observations of one or more strong continuum sources to be used for fringe-finding, "manual" phase calibration, and bandpass calibration. In addition, scans of a continuum source reasonably close to the line target should be scheduled, for use in delay and fringe-rate calibration. The pulse cal generators should be disabled.
Post-processing steps include performing Doppler corrections for the Earth's orbital motion (a correction for Earth rotation is not necessary for VLBA observations since station-based fringe rotation is applied in the correlator); amplitude calibration using single-antenna spectra; fringe fitting the continuum calibrators and applying the results to the line target; referencing phases to a strong spectral feature in the line source itself; deciding whether to use normal synthesis imaging or fringe rate mapping; and then forming a spectral line cube. All these post-processing steps can be done in AIPS.
Data reduction techniques for VLBI spectral line polarimetry are discussed by Kemball, Diamond, & Cotton (1995) and Kemball (1999).
13.3. Pulsar Observations
All special processing required for pulsar observations is supported within the correlator. Details of the available gating and binning options, and their impact on the output data rate, are presented in the pulsar modes and output rate subsections of the DiFX section, respectively.
14. Extended "VLBA-Plus" Arrays
14.1. Introduction
In the interest of enhancing VLBI sensitivity and/or angular resolution, several inter-observatory agreements have been established to support combinations of the VLBA with other radio telescopes. This section describes the observing capabilities of each of these combinations.
14.2. The High Sensitivity Array (HSA)
The HSA comprises the VLBA, phased Very Large Array (VLA), Green Bank Telescope (GBT), Effelsberg, and Arecibo telescopes, and subsets thereof. All of these are equipped with instrumentation compatible with the VLBA observing capabilites described in the Station Signal Processing section. VLBI observations combining the VLBA with any one or more of the other four HSA stations can be requested in a single HSA proposal. Proposal deadlines for the HSA coincide with those for the VLBA alone, as described in Section 15. Further information on "Observing with the High Sensitivity Array" is available in a separate document.
• The VLA is available as a single phased array ("Y27"; no subarrays) with two independently-tunable subband pairs, one polarization pair (RCP+LCP) in the A0/C0 basebands and the other (RCP+LCP) in the B0/D0 basebands. Setups also matching the VLBA PFB and DDC (4- and 8-channel) observing systems are available on the VLA. The VLA must be set up to match the VLBA; mixed modes are not allowed.
Bandwidths must be uniform across the entire VLBI array, and throughout the entire duration of the observation. In particular, VLA phasing and VLBI observing must be carried out at the same bandwidth. Subband bandwidths of 16 MHz and wider are available as a general capability. Bandwidths narrower than 16 MHz may work if the source is strong enough, but are expected to be of limited use, have not been tested, and are available only on a shared-risk basis.
Two adjunct documents: VLBI at the VLA and Guide to Observing with the VLBA: HSA/GMVA/Global VLBI, discuss the available phased-VLA capabilities, and provide instructions for their use.
• The GBT, like the VLBA, transitioned to a new partnership arrangement as of semester 2017A. Time available for VLBI on the GBT has been reduced compared to earlier observing semesters due to this change. Proposers should only include the GBT in the proposal if it is essential for the science and if it is clearly justified in the text.
The GBT is frequency agile (with some limitations) at all its bands.
The GBT's 6 cm receiver is similar to the VLBA's new system, but does differ in converting to circular polarization at ambient temperature. Tests have seen substantial polarization leakage between the RCP and LCP channels. Proposals to use this receiver will be considered only for total-intensity observations. Such proposals should request full dual-polarization modes for both observation and correlation, and careful calibration of the leakage terms should be included in the data analysis.
Proposals including the GBT in an HSA observation must include time to set up the telescope (pointing, focus, etc.) prior to the start of the observation. This can take 0.5-1 hour depending on the frequency.
Further information on use of the GBT may be found in Section 5.7 and Chapter 7 of the GBT Proposer's Guide, in the GBT Observer's Guide, and in the VLBI at the GBT page.
• The Effelsberg telescope supports both of the VLBA observing systems, and is frequency agile at 5 GHz and above.
Further information on use of Effelsberg is available at the Effelsberg HSA page.
• The Arecibo 305-m telescope is currently available only with the PFB observing system. An RSRO project is encouraged to help qualify the 4-channel DDC observing system at Arecibo.
The telescope operates at frequencies up to 10 GHz, and can observe between declinations of -1 and +38 degrees, with transits ranging up to 2.75 hours at 18.5 degrees.
Further information on use of Arecibo is available in The VLBI Guide at Arecibo and The New Users Guide.
14.3. The Global 3mm VLBI Array (GMVA)
The GMVA observes at 3mm wavelength, using 8 VLBA stations, GBT, Effelsberg, Pico Veleta, Onsala, Metsaehovi, Yebes, and Korean VLBI Network (KVN) telescopes. (The VLBA stations in high-humidity locations at HN and SC are not equipped with 3mm receivers.)
The European part of the GMVA is coordinated by the Max-Planck-Institut für Radioastronomie (MPIfR), which hosts the main GMVA website. Proposal deadlines for the GMVA coincide with those for the VLBA alone, as described in Section 15. Observations are block-scheduled in two sessions per year, typically in April-May and September-October, and correlated at the MPIfR in Bonn.
In GMVA observations, all telescopes operate at their highest possible bandwidth. Currently, this is a 2 Gbps data rate in all cases except at the KVN telescopes, which operate in a compatible 1 Gbps mode.
14.4. The European VLBI Network (EVN) and Global cm VLBI
The EVN is a VLBI network of stations operated by an international consortium of institutes (Schilizzi 1995). The EVN home page provides access to the EVN User's Guide. Included in the guide is an EVN Status Table, giving details of current observing capabilities of all EVN stations; and the EVN Call for Proposals, which specifies EVN session dates and the wavelengths to be observed.
The EVN provides proposal, review, and scheduling mechanisms for such programs, and conducts regular sessions of several weeks throughout the year to carry out these observations. Unlike the VLBA, HSA, and GMVA, the EVN operates on a trimester cycle, with proposal deadlines on February 1, June 1, and October 1.
Proposals requesting the EVN in combination with the VLBA or other affiliates are classified as "Global cm VLBI". EVN and Global cm VLBI proposals must be prepared and submitted to the EVN using the EVN's NorthStar Tool. Approved observations will be carried out during EVN sessions.
15. Proposal Preparation, Submission, & Review
Observing time on the VLBA is scheduled on a semester basis, with each semester lasting six months. Proposal deadlines are February 1 and August 1, with the February 1 proposal deadline nominally covering time to be scheduled during the following August through January, and the August 1 deadline covering time to be scheduled from February through July. Time can be requested over multiple semesters if scientifically justified.
A Call for Proposals (see the latest VLBA proposal preparation and submission document for more information) is published approximately four weeks in advance of each semester submission deadline. Proposals must be prepared and submitted using the NRAO Proposal Submission Tool (PST), available via NRAO Interactive Services. Observing proposals may specify the VLBA alone, or in the various extended arrays described in Section 14.
Submitted proposals are reviewed by a Science Review Panel (SRP) in relevant subdisciplines (e.g., solar system, stellar, Galactic, extragalactic, etc.). The SRP's comments and rating are strongly advisory to the Time Allocation Committee (TAC), and the comments of both groups are passed on to the proposers soon after each semi-annual meeting of the TAC, and prior to the next proposal submission deadline. A detailed description of the time allocation process is available.
Approved programs are scheduled by the VLBA scheduling officers, who may be contacted at schedvlb@lbo.us. A "Guide to Using the VLBA" is available, aimed specifically at inexperienced users but also useful to fill in knowledge gaps for more experienced observe
The VLBA SCHED program (Walker 2011) can be used to determine the Greenwich Sidereal Time range during which the VLBI target sources are visible at various stations. This program can also be used to evaluate the u-v plane coverage and synthesized beams provided by the selected array.
An accurate source position service is available to obtain positions accurate enough for correlation. This should be requested simultaneously with the proposal, if not earlier. Requirements for source position accuracy in correlation are discussed in the "Guide to Using the VLBA".
16. Observation Preparation & Execution
Approximately 65% of all VLBA observations are scheduled dynamically, based on array and weather conditions predicted 1-2 days in advance. Users allocated VLBA observing time, either on fixed dates or on a dynamically-scheduled basis, will be sent instructions for preparing observing schedules.
Most VLBA observations are scheduled using the VLBA SCHED program (Walker 2011). SCHED includes a variety of extremely useful facilities cross-linking and error-checking the many options available for VLBA observing, and substantially simplifies the overall scheduling process. A comprehensive SCHED User Manual includes instructions for obtaining and installing the software.
Help with preparing VLBA observations is available through the NRAO Helpdesk.
Each VLBA program is run remotely from the SOC by VLBA operations. No observing assistance by a VLBA user is expected, although VLBA operations should be able to reach the observer by telephone during the program, at a number that should be specified in the schedule file. As the program progresses, the array operator monitors the status of the antennas and the station data path. Logging, calibration, and flagging data are automatically recorded by the monitor and control system. If necessary, the array operator can request local assistance from a site technician at each VLBA station. Recorded media are automatically shipped from each VLBA station to the correlator specified by the observer.
17. Data Archive and Distribution
All output from the VLBA is maintained in the NRAO data archive. The user(s) who proposed the observations retains a proprietary right to access the archived data for an interval of 12 months following the end of correlation of the last observations requested in the original proposal, or a direct extension of that proposal. Thereafter, the archived data are available to any person on request. Data can be obtained from the archive either as multiple correlator output files, or as large FITS files with default calibrations attached.
Distributed data conform to the FITS Interferometry Data Interchange Convention (Greisen 2009), which is read by AIPS task FITLD.
18. Post-Processing Software
AIPS: The Astronomical Image Processing System, a set of programs for the analysis of continuum and line observations, is widely used with VLBA and VLBI data. These programs are supported on a variety of computer operating systems, including Linux and Mac-OS/X. The AIPS Cookbook (NRAO staff, 2014) includes an entire chapter on reducing VLBI data, including discussion of VLBA calibration, space VLBI, polarimetry, and phase referencing; Appendix C provides a step-by-step guide to calibrating many types of VLBA data in AIPS, including VLBA+VLA datasets.
A new ''frozen" version of AIPS is produced each year, and a newer version is updated and made available throughout the calendar year. Observers are encouraged to use a very recent version of AIPS, in order to keep up with ongoing developments in VLBA instrumentation. AIPS home page: http://aips.nrao.edu.
CASA: Common Astronomy Software Applications is the new data reduction package for the Jansky VLA and ALMA. It does not yet offer an end-to-end reduction path for VLBA data. However, CASA does contain imaging and calibration tools that may be of use for VLBI data. CASA home page: http://casa.nrao.edu.
Difmap: Difmap (Shepherd 1997), developed as part of the Caltech VLBI Analysis Programs. provides editing, imaging, self-calibration, and pipelining capabilities in an interactive package. Development has been frozen, and continued support is limited primarily to assistance in installation. Difmap download site: ftp://ftp.astro.caltech.edu/pub/difmap/difmap.html. Contact: M.C. Shepherd, mcs@astro.caltech.edu.
19. References
A list of papers and articles referred to in this document follows. Numerous citations are made to two particular conference proceedings volumes:
- VLBI & the VLBA = Very Long Baseline Interferometry and the VLBA, Astronomical Society of the Pacific Conference Series, Volume 82, 1995, eds. J.A. Zensus, P.J. Diamond, & P.J. Napier.
This conference, held 1993 June 23-20, was the inaugural meeting held as the VLBA entered operational service. An overview of the entire proceedings volume is available. Links in the citations below will retrieve preprints of the relevant chapter, in gzipped PostScript format.
- Synthesis II = Synthesis Imaging in Radio Astronomy II, Astronomical Society of the Pacific Conference Series, Volume 180, 1999, eds. G.B. Taylor, C.L. Carilli, & R.A. Perley.
- Beasley, A.J., & Conway, J.E. 1995, in VLBI & the VLBA, p. 327
- Beasley, A.J., Gordon, D., Peck, A.B., Petrov, L., MacMillan, D.S., Fomalont, E.B., & Ma, C. 2002, Astrophysical Journal Supplement, 141, 13
- Bridle, A.H., & Schwab, F.R. 1999, in Synthesis II, p. 371.
- Briggs, D.S., Schwab, F.R., & Sramek, R.A. 1999, in Synthesis II, p. 127.
- Brisken, W. 2008, VLBA Sensitivity Upgrade Memo # 23.
- Brisken, W., & Deller, A. 2010, VLBA Sensitivity Upgrade Memo # 32.
- Conway, J.E., & Sault, R.J. 1995, in VLBI & the VLBA, p. 309.
- Cornwell, T.J. 1995, in VLBI & the VLBA, p. 39.
- Cornwell, T.J., Braun, R., & Briggs, D.S. 1999, in Synthesis II, p. 151.
- Cornwell, T.J., & Fomalont, E.B. 1999, in Synthesis II, p. 187.
- Cotton, W.D. 1995a, in VLBI & the VLBA, p. 189.
- Cotton, W.D. 1995b, in VLBI & the VLBA, p. 289.
- Cotton, W.D. 1999a, in Synthesis II, p. 111.
- Cotton, W.D. 1999b, in Synthesis II, p. 357.
- Cotton, W.D., Dallacasa, D., Fanti, C., Fanti, R., Foley, A.R., Schilizzi, R.T., & Spencer, R. E. 1997a, Astronomy & Astrophysics, 325, 493.
- Cotton, W.D., Fanti, C., Fanti, R., Dallacasa, D., Foley, A.R., Schilizzi, R.T., & Spencer, R. E. 1997b, Astronomy & Astrophysics, 325, 479.
- Deller, A.T., Tingay, S.J., Bailes, M., & West, C. 2007, Publications of the Astronomical Society of the Pacific, 119, 318.
- Deller, A.T., Brisken, W.F., Phillips, C.J., Morgan, J., Alef, W. Cappallo, R., Middleberg, E., Romney, J., Rotmann, H., Tingay, S. J., & Wayth, R., 2011, Publications of the Astronomical Society of the Pacific, 123, 275.
- Diamond, P.J. 1995, in VLBI & the VLBA, p. 227.
- Fiedler, R., Dennison, B., Johnston, K.J., Waltman, E.B., & Simon, R.S. 1994a, Astrophysical Journal, 430, 581.
- Fiedler, R., Pauls, T., Johnston, K.J., & Dennison, B. 1994b, Astrophysical Journal, 430, 595.
- Fomalont, E., & Kogan, L. 2005, AIPS Memo No. 111.
- Garrett, M. A., Porcas, R. W., Pedlar, A., Muxlow, T. W. B., & Garrington, S. T. 1999, New Astronomy Reviews, 43, 519.
- Greisen, E. W. 2009, AIPS Memo 114.
- Hronek, A., & Walker, R.C. 1996, VLBA Test Memo No. 51.
- Kemball, A.J. 1999, in Synthesis II, p. 499.
- Kemball, A.J., Diamond, P.J., & Cotton, W.D. 1995, Astronomy & Astrophysics Supplement Series, 110, 383.
- Kogan, L. 1995a, VLBA Scientific Memo No. 9.
- Kogan, L. 1995b, VLBA Scientific Memo No. 12.
- Leppänen, K.J. 1993, VLBA Scientific Memo No. 1.
- Leppänen, K.J., Zensus, J.A., & Diamond, P.J. 1995, Astronomical Journal, 110, 2479.
- Markowitz, A., & Wurnig, J. 1998, VLBA Test Memo No. 60.
- Mioduszewski, A. 2004, AIPS Memo No. 110.
- Moran, J.M., & Dhawan, V. 1995, in VLBI & the VLBA, p. 161.
- Napier, P.J. 1995, in VLBI & the VLBA, p. 59.
- Napier, P.J., Bagri, D.S., Clark, B.G., Rogers, A.E.E., Romney, J.D., Thompson, A.R., & Walker, R.C. 1994, Proc. IEEE, 82, 658.
- NRAO staff, 2012, AIPS Cookbook.
- Perley, R.A. 1999a, in Synthesis II, p. 275.
- Perley, R.A. 1999b, in Synthesis II, p. 383.
- Reid, M.J. 1995, in VLBI & the VLBA, p. 209.
- Reid, M.J. 1999, in Synthesis II, p. 481.
- Rogers, A.E.E. 1995, in VLBI & the VLBA, p. 93.
- Romney, J.D. 1999, internal NRAO document reissued as VLBA Sensitivity Upgrade Memo # 2.
- Romney, J.D. 2000, internal NRAO document reissued as VLBA Sensitivity Upgrade Memo # 3.
- Romney, J.D., & Reid, M.J. 2005, Future Directions in High Resolution Astronomy, Astronomical Society of the Pacific Conference Series, Volume 340, 2005.
- Romney, J.D. 2007, VLBA Sensitivity Upgrade Memo # 1.
- Sault, R.J., & Conway, J.E. 1999, in Synthesis II, p. 419.
- Schilizzi, R.T. 1995, in VLBI & the VLBA, p. 397.
- Shepherd, M.C. 1997, ADASS IV, Astronomical Society of the Pacific Conference Series, Volume 125, eds. G. Hunt & H.E. Payne, p. 77. http://www.nrao.edu/meetings/proceedings.shtml
- Taylor, G.B. 1998, Astrophysical Journal, 506, 637.
- Taylor, G.B. 2000, Astrophysical Journal, 533, 95.
- Taylor, G.B., Lonsdale, C.J., et al. 2004, Mapping the Future of VLBI Science in the U.S.
- Taylor, G.B., & Myers, S.T. 2000, VLBA Scientific Memo No. 26.
- Thompson, A.R. 1995, in VLBI & the VLBA, p. 73.
- Walker, R.C. 1995a, in VLBI & the VLBA, p. 133.
- Walker, R.C. 1995b, in VLBI & the VLBA, p. 247.
- Walker, R.C. 1999, in Synthesis II, p. 433.
- Walker, R.C. & Chatterjee, S. 1999, VLBA Scientific Memo No. 23.
- Walker, C., Durand, S., Kutz, C., & Hayward, R. 2007a, VLBA Sensitivity Upgrade Memo # 10.
- Walker, C., Romney, J., Brisken, W., & Durand, S. 2007b, VLBA Sensitivity Upgrade Memo # 15.
- Walker, C., Durand, S., Kutz, C., & Hayward, R. 2008, VLBA Sensitivity Upgrade Memo # 21.
- Walker, R.C. 2011 et seq., The SCHED User Manual. http://www.aoc.nrao.edu/software/sched/
- Wrobel, J.M. 1995, in VLBI & the VLBA, p. 411.
- Wrobel, J.M., & Walker, R.C. 1999, in Synthesis II, p. 171.
- Wrobel, J.M., Walker, R.C., Benson, J.M., & Beasley, A.J. 2000, VLBA Scientific Memo No. 24.
- Zensus, J.A., Taylor, G.B., & Wrobel, J.M. 1998, IAU Colloquium 164: Radio Emission from Galactic and Extragalactic Compact Sources, Astronomical Society of the Pacific Conference Series, Volume 144. http://www.nrao.edu/meetings/past.shtml
20. Editor's Notes
The VLBA Observational Status Summary is the collective work of innumerable individuals who wrote and edited the text, commented on draft material, and implemented the capabilities described herein, during the 23 years since the VLBA's dedication in 1993. We thank these many colleagues for their contributions.
Editors:
1992 - 2005 | J. M. Wrobel |
2005 - 2009 | J. S. Ulvestad |
2009 - 2016 | J. D. Romney |
The current editor is responsible for the most recent revisions, and thus is the best contact for readers who may have questions on the material, or suggestions that would enhance the usefulness of this guide.
Versions issued prior to July 2012 were based on a LaTeX source, from which PS and PDF versions were produced. A hypertext version was generated using the LaTeX2HTML translator. The LaTeX version of January 2012 was ported to Plone in April 2012.