We characterize the accuracy of linear-polarization mosaics made using the Atacama Large Millimeter/submillimeter Array (ALMA). First, we observed the bright, highly linearly polarized blazar 3C 279 at Bands 3, 5, 6, and 7 (3 mm, 1.6 mm, 1.3 mm, and 0.87 mm, respectively). At each band, we measured the blazar's polarization on an 11$\times$11 grid of evenly-spaced offset pointings covering the full-width at half-maximum (FWHM) area of the primary beam. After applying calibration solutions derived from the on-axis pointing of 3C 279 to all of the on- and off-axis data, we find that the residual polarization errors across the primary beam are similar at all frequencies: the residual errors in linear polarization fraction $P_\textrm{frac}$ and polarization position angle $\chi$ are $\lesssim$0.001 ($\lesssim$0.1% of Stokes $I$) and $\lesssim$1$^\circ$ near the center of the primary beam; the errors increase to $\sim$0.003-0.005 ($\sim$0.3-0.5% of Stokes $I$) and $\sim$1-5$^\circ$ near the FWHM as a result of the asymmetric beam patterns in the (linearly polarized) $Q$ and $U$ maps. We see the expected double-lobed "beam squint" pattern in the circular polarization (Stokes $V$) maps. Second, to test the polarization accuracy in a typical ALMA project, we performed observations of continuum linear polarization toward the Kleinmann-Low nebula in Orion (Orion-KL) using several mosaic patterns at Bands 3 and 6. We show that after mosaicking, the residual off-axis errors decrease as a result of overlapping multiple pointings. Finally, we compare the ALMA mosaics with an archival 1.3 mm CARMA polarization mosaic of Orion-KL and find good consistency in the polarization patterns.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

In 2017, an Atacama Large Millimeter/submillimeter Array (ALMA) high-frequency long baseline campaign was organized to test image capabilities with baselines up to 16 km at submillimeter (submm) wavelengths. We investigated image qualities using ALMA receiver Bands 7, 8, 9, and 10 (285-875 GHz) by adopting band-to-band (B2B) phase referencing in which a phase calibrator is tracked at a lower frequency. For B2B phase referencing, it is expected that a closer phase calibrator to a target can be used, comparing to standard in-band phase referencing. In the first step, it is ensured that an instrumental phase offset difference between low- and high-frequency Bands can be removed using a differential gain calibration in which a phase calibrator is certainly detected while frequency switching. In the next step, comparative experiments are arranged to investigate the image quality between B2B and in-band phase referencing with phase calibrators at various separation angles. In the final step, we conducted long baseline imaging tests for a quasar at 289 GHz in Band 7 and 405 GHz in Band 8 and complex structure sources of HL Tau and VY CMa at ~670 GHz in Band 9. The B2B phase referencing was successfully applied, allowing us to achieve an angular resolution of 14x11 and 10x8 mas for HL Tau and VY CMa, respectively. There is a high probability of finding a low-frequency calibrator within 5.4 deg in B2B phase referencing, bright enough to use an 8 s scan length combined with a 7.5 GHz bandwidth.

Using the Atacama Large Millimeter/submillimeter Array (ALMA), we observed the young Herbig star HD 100546, host to a prominent disk with a deep, wide gap in the dust. The high-resolution 1.3 mm continuum observation reveals fine radial and azimuthal substructures in the form of a complex maze of ridges and trenches sculpting a dust ring. The $^{12}$CO(2-1) channel maps are modulated by wiggles or kinks that deviate from Keplerian kinematics particularly over the continuum ring, where deviations span 90$^\circ$ in azimuth, covering 5 km s$^{-1}$. The most pronounced wiggle resembles the imprint of an embedded massive planet of at least 5 M$_{\rm Jup}$ predicted from previous hydrodynamical simulations (Perez, Casassus, & Benitez-Llambay 2018). Such planet is expected to open a deep gap in both gas and dust density fields within a few orbital timescales, yet the kinematic wiggles lie near ridges in the continuum. The lesser strength of the wiggles in the $^{13}$CO and C$^{18}$O isotopologues show that the kinematic signature weakens at lower disk heights, and suggests qualitatively that it is due to vertical flows in the disk surface. Within the gap, the velocity field transitions from Keplerian to strongly non-Keplerian via a twist in position angle, suggesting the presence of another perturber and/or an inner warp. We also present VLT/SPHERE sparse aperture masking data which recovers scattered light emission from the gap's edges but shows no evidence for signal within the gap, discarding a stellar binary origin for its opening.

We present a detailed description of the special procedures for calibration and quality assurance of Atacama Large Millimeter/submillimeter Array (ALMA) observations in Very Long Baseline Interferometry (VLBI) mode. These procedures are required to turn the phased ALMA array into a fully calibrated VLBI station. As an illustration of these methodologies, we present full-polarization observations carried out with ALMA as a phased array at 3mm (Band 3) and 1.3mm (Band 6) as part of Cycle-4. These are the first VLBI science observations conducted with ALMA and were obtained during a 2017 VLBI campaign in concert with other telescopes worldwide as part of the Global mm-VLBI Array (GMVA, April 1-3) and the Event Horizon Telescope (EHT, April 5-11) in ALMA Bands 3 and 6, respectively.

We present a new package for joint deconvolution of ALMA 12m, 7m, and Total Power (TP) data, dubbed ``Total Power Map to Visibilities (TP2VIS)". It converts a TP (single-dish) map into visibilities on the CASA platform, which can be input into deconvolvers (e.g., CLEAN) along with 12m and 7m visibilities. A manual is presented in the Github repository (https://github.com/tp2vis/distribute). Combining data from the different ALMA arrays is a driver for a number of science topics, namely those that probe size scales of extended and compact structures simultaneously. We test TP2VIS using model images, one with a single Gaussian and another that mimics the internal structures of giant molecular clouds. The result shows that the better uv coverage with TP2VIS visibilities helps the deconvolution process and reproduces the model image within errors of only 5% over two orders of magnitude in flux.

Using the $\approx$15km ALMA long baselines, we imaged the Stokes $I$ emission and linearly polarized intensity ($PI$) in the 1.1-mm continuum band of a very young intermediate-mass protostellar source, MMS 6, in the Orion Molecular Cloud-3. The achieved angular resolution, $0''.02{\times}0''.03$ ($\approx$10 AU), shows for the first time a wealth of data on the dust emission polarization in the central 200 AU of a protostar. The $PI$ peak is offset to the south-west (SW) by $\approx$20 AU with respect to the Stokes $I$ peak. Its polarization degree is 11 % with its $E$-vector orientation of P.A.${\approx}135^{\circ}$. A partial ring-like structure with a radius of $\approx$80 AU is detected in $PI$ but not in the Stokes $I$. NW (north-west) and SE (south-east) parts of the ring are bright with a high polarization degree of $\gtrsim$10 %, and their $E$-vector orientations are roughly orthogonal to those observed near the center. We also detected arm-like polarized structures, extending to 1000 AU scale to the north, with the $E$-vectors aligned along the minor axis of the structures. We explored possible origins of the polarized emission comparing with magnetohydrodynamical (MHD) simulations of the toroidal wrapping of the magnetic field. The simulations are consistent with the $PI$ emission in the ring-like and the extended arm-like structures observed with ALMA. However, the current simulations do not completely reproduce observed polarization characteristics in the central 50 AU. Although the self-scattering model can explain the polarization pattern and positional offset between the Stokes $I$ and $PI$, this model is not able to reproduce the observed high degree of polarization.

We present ALMA images of the sub-mm continuum polarisation and spectral index of the protoplanetary ringed disk HD163296. The polarisation fraction at 870{\mu}m is measured to be ~0.9% in the central core and generally increases with radius along the disk major axis. It peaks in the gaps between the dust rings, and the largest value (~4%) is found between rings 1 and 2. The polarisation vectors are aligned with the disk minor axis in the central core, but become more azimuthal in the gaps, twisting by up to +/-9degrees in the gap between rings 1 and 2. These general characteristics are consistent with a model of self-scattered radiation in the ringed structure, without requiring an additional dust alignment mechanism. The 870/1300{\mu}m dust spectral index exhibits minima in the centre and the inner rings, suggesting these regions have high optical depths. However, further refinement of the dust or the disk model at higher resolution is needed to reproduce simultaneously the observed degree of polarisation and the low spectral index.

In this lecture, we describe a number of advanced gain calibration techniques. In particular, self-calibration is an important tool in interferometric imaging at all wavelengths. It allows the observer to determine and remove residual phase and amplitude errors that remain in the data after normal calibration while simultaneously producing a more accurate and more sensitive image of the target source. We describe the basic technique of self-calibration and attempt to dispel some common misconceptions. We proceed to give a range of useful tips, and provide continuum, spectral line, and mosaic self-calibration examples using ALMA data. We also discuss fast switching and radiometric phase correction along with advanced phase transfer techniques that can become advantageous or even essential at high frequency where the density of sufficiently bright phase calibrators becomes sparse.

The Atacama Large millimetre/submillimetre Array (ALMA) makes use of water vapour radiometers (WVR), which monitor the atmospheric water vapour line at 183 GHz along the line of sight above each antenna to correct for phase delays introduced by the wet component of the troposphere. The application of WVR derived phase corrections improve the image quality and facilitate successful observations in weather conditions that were classically marginal or poor. We present work to indicate that a scaling factor applied to the WVR solutions can act to further improve the phase stability and image quality of ALMA data. We find reduced phase noise statistics for 62 out of 75 datasets from the long-baseline science verification campaign after a WVR scaling factor is applied. The improvement of phase noise translates to an expected coherence improvement in 39 datasets. When imaging the bandpass source, we find 33 of the 39 datasets show an improvement in the signal-to-noise ratio (S/N) between a few to ~30 percent. There are 23 datasets where the S/N of the science image is improved: 6 by <1%, 11 between 1 and 5%, and 6 above 5%. The higher frequencies studied (band 6 and band 7) are those most improved, specifically datasets with low precipitable water vapour (PWV), <1mm, where the dominance of the wet component is reduced. Although these improvements are not profound, phase stability improvements via the WVR scaling factor come into play for the higher frequency (>450 GHz) and long-baseline (>5km) observations. These inherently have poorer phase stability and are taken in low PWV (<1mm) conditions for which we find the scaling to be most effective. A promising explanation for the scaling factor is the mixing of dry and wet air components, although other origins are discussed. We have produced a python code to allow ALMA users to undertake WVR scaling tests and make improvements to their data.

This paper presents the first detailed investigation of the characteristics of mm/submm phase fluctuation and phase correction methods obtained using ALMA with baseline lengths up to ~15 km. Most of the spatial structure functions (SSFs) show that the phase fluctuation increases as a function of baseline length, with a power-law slope of ~0.6. In many cases, we find that the slope becomes shallower (average of ~0.2-0.3) at baseline lengths longer than ~1 km, namely showing a turn-over in SSF. The phase correction method using water vapor radiometers (WVRs) works well, especially for the cases where PWV >1 mm, which reduces the degree of phase fluctuations by a factor of two in many cases. However, phase fluctuations still remain after the WVR phase correction, suggesting the existence of other turbulent constituent that cause the phase fluctuation. This is supported by occasional SSFs that do not exhibit any turn-over; these are only seen when the PWV is low or after WVR phase correction. This means that the phase fluctuation caused by this turbulent constituent is inherently smaller than that caused by water vapor. Since there is no turn-over in the SSF up to the maximum baseline length of ~15 km, this turbulent constituent must have scale height of 10 km or more, and thus cannot be water vapor, whose scale height is around 1 km. This large scale height turbulent constituent is likely to be water ice or a dry component. Excess path length fluctuation after the WVR phase correction at a baseline length of 10 km is large (>200 micron), which is significant for high frequency (>450 GHz or <700 micron) observations. These results suggest the need for an additional phase correction method, such as fast switching, in addition to the WVR phase correction. We simulated the fast switching, and the result suggests that it works well, with shorter cycle times linearly improving the coherence.

In a radio interferometer, the determination of geometrical antenna positions relies on accurate calibration of the dry and wet delay of the atmosphere above each antenna. For the Atacama Large Millimeter/Submillimeter Array (ALMA), which has baseline lengths up to 16 kilometers, the geography of the site forces the height above mean sea level of the more distant antenna pads to be significantly lower than the central array. Thus, both the ground level meteorological values and the total water column can be quite different between antennas in the extended configurations. During 2015, a network of six additional weather stations was installed to monitor pressure, temperature, relative humidity and wind velocity, in order to test whether inclusion of these parameters could improve the repeatability of antenna position determinations in these configurations. We present an analysis of the data obtained during the ALMA Long Baseline Campaign of Oct. through Nov. 2015. The repeatability of antenna position measurements typically degrades as a function of antenna distance. Also, the scatter is more than three times worse in the vertical direction than in the local tangent plane, suggesting that a systematic effect is limiting the measurements. So far we have explored correcting the delay model for deviations from hydrostatic equilibrium in the measured air pressure and separating the partial pressure of water from the total pressure using water vapor radiometer (WVR) data. Correcting for these combined effects still does not provide a good match to the residual position errors in the vertical direction. One hypothesis is that the current model of water vapor may be too simple to fully remove the day-to-day variations in the wet delay. We describe possible avenues of improvement, including measuring and applying more accurate values of the sky coupling efficiency of the WVRs.

We present full-polarization observations of the compact, steep-spectrum radio quasar 3C~286 made with the ALMA at 1.3~mm. These are the first full-polarization ALMA observations, which were obtained in the framework of Science Verification. A bright core and a south-west component are detected in the total intensity image, similar to previous centimeter images. Polarized emission is also detected toward both components. The fractional polarization of the core is about 17\%, this is higher than the fractional polarization at centimeter wavelengths, suggesting that the magnetic field is even more ordered in the millimeter radio core than it is further downstream in the jet. The observed polarization position angle (or EVPA) in the core is $\sim$\,$39^{\circ}$, which confirms the trend that the EVPA slowly increases from centimeter to millimeter wavelengths. With the aid of multi-frequency VLBI observations, we argue that this EVPA change is associated with the frequency-dependent core position. We also report a serendipitous detection of a sub-mJy source in the field of view, which is likely to be a submillimeter galaxy.

Magnetized plasmas traversed by linearly polarized light will reveal their presence by the frequency dependent Faraday rotation of the angle of polarization. The regions surrounding the black holes powering the jets in AGNs are expected to have dense magnetized plasmas, possibly giving rise to very large Faraday rotations. Compact steep spectrum (CSS) sources are good candidates to search for very large Faraday rotated components as they contain compact emission from close to the black hole and many are strongly depolarized at centimeter wavelengths as expected from strong Faraday effects. We present data on several CSS sources (3C48, 3C138 and 3C147) observed with the VLA at frequencies between 20 and 48 GHz in the most extended configuration. Large, but not excessive rotation measures are reported.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations from the 2014 Long Baseline Campaign in dust continuum and spectral line emission from the HL Tau region. The continuum images at wavelengths of 2.9, 1.3, and 0.87 mm have unprecedented angular resolutions of 0.075 arcseconds (10 AU) to 0.025 arcseconds (3.5 AU), revealing an astonishing level of detail in the circumstellar disk surrounding the young solar analogue HL Tau, with a pattern of bright and dark rings observed at all wavelengths. By fitting ellipses to the most distinct rings, we measure precise values for the disk inclination (46.72pm0.05 degrees) and position angle (+138.02pm0.07 degrees). We obtain a high-fidelity image of the 1.0 mm spectral index ($\alpha$), which ranges from $\alpha\sim2.0$ in the optically-thick central peak and two brightest rings, increasing to 2.3-3.0 in the dark rings. The dark rings are not devoid of emission, we estimate a grain emissivity index of 0.8 for the innermost dark ring and lower for subsequent dark rings, consistent with some degree of grain growth and evolution. Additional clues that the rings arise from planet formation include an increase in their central offsets with radius and the presence of numerous orbital resonances. At a resolution of 35 AU, we resolve the molecular component of the disk in HCO+ (1-0) which exhibits a pattern over LSR velocities from 2-12 km/s consistent with Keplerian motion around a ~1.3 solar mass star, although complicated by absorption at low blue-shifted velocities. We also serendipitously detect and resolve the nearby protostars XZ Tau (A/B) and LkHa358 at 2.9 mm.

We present Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm continuum images of the asteroid 3 Juno obtained with an angular resolution of 0.042 arcseconds (60 km at 1.97 AU). The data were obtained over a single 4.4 hr interval, which covers 60% of the 7.2 hr rotation period, approximately centered on local transit. A sequence of ten consecutive images reveals continuous changes in the asteroid's profile and apparent shape, in good agreement with the sky projection of the three-dimensional model of the Database of Asteroid Models from Inversion Techniques. We measure a geometric mean diameter of 259pm4 km, in good agreement with past estimates from a variety of techniques and wavelengths. Due to the viewing angle and inclination of the rotational pole, the southern hemisphere dominates all of the images. The median peak brightness temperature is 215pm13 K, while the median over the whole surface is 197pm15 K. With the unprecedented resolution of ALMA, we find that the brightness temperature varies across the surface with higher values correlated to the subsolar point and afternoon areas, and lower values beyond the evening terminator. The dominance of the subsolar point is accentuated in the final four images, suggesting a reduction in the thermal inertia of the regolith at the corresponding longitudes, which are possibly correlated to the location of the putative large impact crater. These results demonstrate ALMA's potential to resolve thermal emission from the surface of main belt asteroids, and to measure accurately their position, geometric shape, rotational period, and soil characteristics.

We present initial results of very high resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations of the $z$=3.042 gravitationally lensed galaxy HATLAS J090311.6+003906 (SDP.81). These observations were carried out using a very extended configuration as part of Science Verification for the 2014 ALMA Long Baseline Campaign, with baselines of up to 15 km. We present continuum imaging at 151, 236 and 290 GHz, at unprecedented angular resolutions as fine as 23 milliarcseconds (mas), corresponding to an un-magnified spatial scale of ~180 pc at z=3.042. The ALMA images clearly show two main gravitational arc components of an Einstein ring, with emission tracing a radius of ~1.5". We also present imaging of CO(10-9), CO(8-7), CO(5-4) and H2O line emission. The CO emission, at an angular resolution of ~170 mas, is found to broadly trace the gravitational arc structures but with differing morphologies between the CO transitions and compared to the dust continuum. Our detection of H2O line emission, using only the shortest baselines, provides the most resolved detection to date of thermal H2O emission in an extragalactic source. The ALMA continuum and spectral line fluxes are consistent with previous Plateau de Bure Interferometer and Submillimeter Array observations despite the impressive increase in angular resolution. Finally, we detect weak unresolved continuum emission from a position that is spatially coincident with the center of the lens, with a spectral index that is consistent with emission from the core of the foreground lensing galaxy.

We report ALMA Band 7 (350 GHz) imaging at 0.4 - 0.6arcsec resolution and Band 9 (696 GHz) at ~0.25arcsec resolution of the luminous IR galaxies Arp 220 and NGC 6240. The long wavelength dust continuum is used to estimate ISM masses for Arp 220 East, West and NGC 6240 of 1.9, 4.2 and 1.6x10^9 msun within radii of 69, 65 and 190 pc. The HCN emission was modeled to derive the emissivity distribution as a function of radius and the kinematics of each nuclear disk, yielding dynamical masses consistent with the masses and sizes derived from the dust emission. In Arp 220, the major dust and gas concentrations are at radii less than 50 pc in both counter-rotating nuclear disks. The thickness of the disks in Arp 220estimated from the velocity dispersion and rotation velocities are 10-20 pc and the mean gas densities are n_H2 ~10^5 cm^-3 at R < 50 pc. We develop an analytic treatment for the molecular excitation (including photon trapping), yielding volume densities for both the HCN and CS emission with n_H2 ~2x10^5 cm^-3. The agreement of the mean density from the total mass and size with that required for excitation suggests that the volume is essentially filled with dense gas, i.e. it is not cloudy or like swiss cheese.

We present the phase characteristics study of the Atacama Large Millimeter/submillimeter Array (ALMA) long (up to 3 km) baseline, which is the longest baseline tested so far using ALMA. The data consist of long time-scale (10 - 20 minutes) measurements on a strong point source (i.e., bright quasar) at various frequency bands (bands 3, 6, and 7, which correspond to the frequencies of about 88 GHz, 232 GHz, and 336 GHz). Water vapor radiometer (WVR) phase correction works well even at long baselines, and the efficiency is better at higher PWV (>1 mm) condition, consistent with the past studies. We calculate the spatial structure function of phase fluctuation, and display that the phase fluctuation (i.e., rms phase) increases as a function of baseline length, and some data sets show turn-over around several hundred meters to 1 km and being almost constant at longer baselines. This is the first millimeter/submillimeter structure function at this long baseline length, and to show the turn-over of the structure function. Furthermore, the observation of the turn-over indicates that even if the ALMA baseline length extends to the planned longest baseline of 15 km, fringes will be detected at a similar rms phase fluctuation as that at a few km baseline lengths. We also calculate the coherence time using the 3 km baseline data, and the results indicate that the coherence time for band 3 is longer than 400 seconds in most of the data (both in the raw and WVR-corrected data). For bands 6 and 7, WVR-corrected data have about twice longer coherence time, but it is better to use fast switching method to avoid the coherence loss.

Planetary ephemerides have been developed and improved over centuries. They are a fundamental tool for understanding solar system dynamics, and essential for planetary and small body mass determinations, occultation predictions, high-precision tests of general relativity, pulsar timing, and interplanetary spacecraft navigation. This paper presents recent results from a continuing program of high-precision astrometric very-long-baseline interferometry (VLBI) observations of the Cassini spacecraft orbiting Saturn, using the Very Long Baseline Array (VLBA). We have previously shown that VLBA measurements can be combined with spacecraft orbit determinations from Doppler and range tracking and VLBI links to the inertial extragalactic reference frame (ICRF) to provide the most accurate barycentric positions currently available for Saturn. Here we report an additional five years of VLBA observations along with improved phase reference source positions, resulting in an improvement in residuals with respect to the Jet Propulsion Laboratory's dynamical ephemeris.

Radio source counts constrain galaxy populations and evolution, as well as the global star formation history. However, there is considerable disagreement among the published 1.4-GHz source counts below 100 microJy. Here we present a statistical method for estimating the microJy and even sub-microJy source count using new deep wide-band 3-GHz data in the Lockman Hole from the Karl G. Jansky Very Large Array (VLA). We analyzed the confusion amplitude distribution P(D), which provides a fresh approach in the form of a more robust model, with a comprehensive error analysis. We tested this method on a large-scale simulation, incorporating clustering and finite source sizes. We discuss in detail our statistical methods for fitting using Monte Carlo Markov chains, handling correlations, and systematic errors from the use of wide-band radio interferometric data. We demonstrated that the source count can be constrained down to 50 nJy, a factor of 20 below the rms confusion. We found the differential source count near 10 microJy to have a slope of -1.7, decreasing to about -1.4 at fainter flux densities. At 3GHz the rms confusion in an 8arcsec FWHM beam is ~ 1.2 microJy/beam, and a radio background temperature ~ 14mK. Our counts are broadly consistent with published evolutionary models. With these results we were also able to constrain the peak of the Euclidean normalized differential source count of any possible new radio populations that would contribute to the cosmic radio background down to 50 nJy.

Gaseous giant planet formation is thought to occur in the first few million years following stellar birth. Models predict that giant planet formation carves a deep gap in the dust component (shallower in the gas). Infrared observations of the disk around the young star HD142527, at ~140pc, found an inner disk ~10AU in radius, surrounded by a particularly large gap, with a disrupted outer disk beyond 140AU, indicative of a perturbing planetary-mass body at ~90 AU. From radio observations, the bulk mass is molecular and lies in the outer disk, whose continuum emission has a horseshoe morphology. The vigorous stellar accretion rate would deplete the inner disk in less than a year, so in order to sustain the observed accretion, matter must flow from the outer-disk into the cavity and cross the gap. In dynamical models, the putative protoplanets channel outer-disk material into gap-crossing bridges that feed stellar accretion through the inner disk. Here we report observations with the Atacama Large Millimetre Array (ALMA) that reveal diffuse CO gas inside the gap, with denser HCO+ gas along gap-crossing filaments, and that confirm the horseshoe morphology of the outer disk. The estimated flow rate of the gas is in the range 7E-9 to 2E-7 Msun/yr, which is sufficient to maintain accretion onto the star at the present rate.

In the past decade, our understanding of galaxy evolution has been revolutionized by the discovery that luminous, dusty, starburst galaxies were 1,000 times more abundant in the early Universe than at present. It has, however, been difficult to measure the complete redshift 2 distribution of these objects, especially at the highest redshifts (z > 4). Here we report a redshift survey at a wavelength of three millimeters, targeting carbon monoxide line emission from the star-forming molecular gas in the direction of extraordinarily bright millimetrewave-selected sources. High-resolution imaging demonstrates that these sources are strongly gravitationally lensed by foreground galaxies. We detect spectral lines in 23 out of 26 sources and multiple lines in 12 of those 23 sources, from which we obtain robust, unambiguous redshifts. At least 10 of the sources are found to lie at z > 4, indicating that the fraction of dusty starburst galaxies at high redshifts is greater than previously thought. Models of lens geometries in the sample indicate that the background objects are ultra-luminous infrared galaxies, powered by extreme bursts of star formation.

We present Atacama Large Millimeter/submillimeter Array (ALMA) 860 micrometer imaging of four high-redshift (z=2.8-5.7) dusty sources that were detected using the South Pole Telescope (SPT) at 1.4 mm and are not seen in existing radio to far-infrared catalogs. At 1.5 arcsec resolution, the ALMA data reveal multiple images of each submillimeter source, separated by 1-3 arcsec, consistent with strong lensing by intervening galaxies visible in near-IR imaging of these sources. We describe a gravitational lens modeling procedure that operates on the measured visibilities and incorporates self-calibration-like antenna phase corrections as part of the model optimization, which we use to interpret the source structure. Lens models indicate that SPT0346-52, located at z=5.7, is one of the most luminous and intensely star-forming sources in the universe with a lensing corrected FIR luminosity of 3.7 X 10^13 L_sun and star formation surface density of 4200 M_sun yr^-1 kpc^-2. We find magnification factors of 5 to 22, with lens Einstein radii of 1.1-2.0 arcsec and Einstein enclosed masses of 1.6-7.2x10^11 M_sun. These observations confirm the lensing origin of these objects, allow us to measure the their intrinsic sizes and luminosities, and demonstrate the important role that ALMA will play in the interpretation of lensed submillimeter sources.

Deep radio observations at 1.4GHz for the Extended Chandra Deep Field South were performed in June through September of 2007 and presented in a first data release (Miller et al. 2008). The survey was made using six separate pointings of the Very Large Array (VLA) with over 40 hours of observation per pointing. In the current paper, we improve on the data reduction to produce a second data release (DR2) mosaic image. This DR2 image covers an area of about a third of a square degree and reaches a best rms sensitivity of 6 uJy and has a typical sensitivity of 7.4 uJy per 2.8" by 1.6" beam. We also present a more comprehensive catalog, including sources down to peak flux densities of five or more times the local rms noise along with information on source sizes and relevant pointing data. We discuss in some detail the consideration of whether sources are resolved under the complication of a radio image created as a mosaic of separate pointings each suffering some degree of bandwidth smearing, and the accurate evaluation of the flux densities of such sources. Finally, the radio morphologies and optical/near-IR counterpart identifications (Bonzini et al. 2012) are used to identify 17 likely multiple-component sources and arrive at a catalog of 883 radio sources, which is roughly double the number of sources contained in the first data release.

Atacama Large Millimeter/submillimeter Array (ALMA) will be the world largest mm/submm interferometer, and currently the Early Science is ongoing, together with the commissioning and science verification (CSV). Here we present a study of the temporal phase stability of the entire ALMA system from antennas to the correlator. We verified the temporal phase stability of ALMA using data, taken during the last two years of CSV activities. The data consist of integrations on strong point sources (i.e., bright quasars) at various frequency bands, and at various baseline lengths (up to 600 m). From the observations of strong quasars for a long time (from a few tens of minutes, up to an hour), we derived the 2-point Allan Standard Deviation after the atmospheric phase correction using the 183 GHz Water Vapor Radiometer (WVR) installed in each 12 m antenna, and confirmed that the phase stability of all the baselines reached the ALMA specification. Since we applied the WVR phase correction to all the data mentioned above, we also studied the effectiveness of the WVR phase correction at various frequencies, baseline lengths, and weather conditions. The phase stability often improves a factor of 2 - 3 after the correction, and sometimes a factor of 7 improvement can be obtained. However, the corrected data still displays an increasing phase fluctuation as a function of baseline length, suggesting that the dry component (e.g., N2 and O2) in the atmosphere also contributes the phase fluctuation in the data, although the imperfection of the WVR phase correction cannot be ruled out at this moment.

We used the Karl G. Jansky Very Large Array (VLA) to image one primary beam area at 3 GHz with 8 arcsec FWHM resolution and 1.0 microJy/beam rms noise near the pointing center. The P(D) distribution from the central 10 arcmin of this confusion-limited image constrains the count of discrete sources in the 1 < S(microJy/beam) < 10 range. At this level the brightness-weighted differential count S^2 n(S) is converging rapidly, as predicted by evolutionary models in which the faintest radio sources are star-forming galaxies; and ~96$% of the background originating in galaxies has been resolved into discrete sources. About 63% of the radio background is produced by AGNs, and the remaining 37% comes from star-forming galaxies that obey the far-infrared (FIR) / radio correlation and account for most of the FIR background at lambda = 160 microns. Our new data confirm that radio sources powered by AGNs and star formation evolve at about the same rate, a result consistent with AGN feedback and the rough correlation of black hole and bulge stellar masses. The confusion at centimeter wavelengths is low enough that neither the planned SKA nor its pathfinder ASKAP EMU survey should be confusion limited, and the ultimate source detection limit imposed by "natural" confusion is < 0.01 microJy at 1.4 GHz. If discrete sources dominate the bright extragalactic background reported by ARCADE2 at 3.3 GHz, they cannot be located in or near galaxies and most are < 0.03 microJy at 1.4 GHz.

This paper presents accurate absolute positions from a 24 GHz Very Long Baseline Array (VLBA) search for compact extragalactic sources in an area where the density of known calibrators with precise coordinates is low. The goals were to identify additional sources suitable for use as phase calibrators for galactic sources, determine their precise positions, and produce radio images. In order to achieve these goals, we developed a new software package, PIMA, for determining group delays from wide-band data with much lower detection limit. With the use of PIMA we have detected 327 sources out of 487 targets observed in three 24 hour VLBA experiments. Among the 327 detected objects, 176 are within 10 degrees of the Galactic plane. This VGaPS catalogue of source positions, plots of correlated flux density versus projected baseline length, contour plots, as well as weighted CLEAN images and calibrated visibility data in FITS format, are available on the Web at http://astrogeo.org/vgaps. Approximately one half of objects from the 24 GHz catalogue were observed at dual band 8.6 GHz and 2.3 GHz experiments. Position differences at 24 GHz versus 8.6/2.3 GHz for all but two objects on average are strictly within reported uncertainties. We found that for two objects with complex structure positions at different frequencies correspond to different components of a source.

Four close radio sources in the International Celestial Reference Frame (ICRF) catalog were observed using phase referencing with the VLBA at 43, 23 and 8.6 GHz, and with VERA at 23 GHz over a one year period. The goal was to determine the stability of the radio cores, and to assess structure effects associated with positions in the ICRF. Although the four sources were compact at 8.6 GHz, the VLBA images at 43 GHz with 0.3-mas resolution showed that all were composed of several components. A component in each source was identified as the radio core using some or all of the following emission properties: compactness, spectral index, location at the end of the extended emission region, and stationary in the sky. Over the observing period, the relative positions between the four radio cores were constant to 0.02 mas---the phase referencing positional accuracy obtained at 23 and 43 GHz among the sources---suggesting that once a radio core is identified, it remains stationary in the sky to this accuracy. Other radio components in two of the four sources had detectable motion in the radio jet direction. Comparison of the 23 and 43 GHz VLBA images with the VLBA 8.6 GHz images and the ICRF positions suggests that some ICRF positions are dominated by a moving jet component; hence, they can be displaced up to 0.5 mas from the radio core, and may also reflect the motion of the jet component. Future astrometric efforts to determine a more accurate quasar reference frame at 23 and 43 GHz and from the VLBI2010 project are discussed, and supporting VLBA or EVN observations of ICRF sources at 43 GHz are recommended in order to determine the internal structure of the sources. A future collaboration between the radio (ICRF) and the optical frame of Gaia is discussed.

The planetary ephemeris is an essential tool for interplanetary spacecraft navigation, studies of solar system dynamics (including, for example, barycenter corrections for pulsar timing ephemeredes), the prediction of occultations, and tests of general relativity. We are carrying out a series of astrometric VLBI observations of the Cassini spacecraft currently in orbit around Saturn, using the Very Long Baseline Array (VLBA). These observations provide positions for the center of mass of Saturn in the International Celestial Reference Frame (ICRF) with accuracies ~0.3 milli-arcsecond (1.5 nrad), or about 2 km at the average distance of Saturn. This paper reports results from eight observing epochs between 2006 October and 2009 April. These data are combined with two VLBA observations by other investigators in 2004 and a Cassini-based gravitational deflection measurement by Fomalont et al. in 2009 to constrain a new ephemeris (DE 422). The DE 422 post-fit residuals for Saturn with respect to the VLBA data are generally 0.2 mas, but additional observations are needed to improve the positions of all of our phase reference sources to this level. Over time we expect to be able to improve the accuracy of all three coordinates in the Saturn ephemeris (latitude, longitude, and range) by a factor of at least three. This will represent a significant improvement not just in the Saturn ephemeris but also in the link between the inner and outer solar system ephemeredes and in the link to the inertial ICRF.

The next-generation space VLBI mission, VSOP-2, is expected to provide unprecedented spatial resolutions at 8.4, 22, and 43GHz. In this report, phase referencing with VSOP-2 is examined in detail based on a simulation tool called ARIS. The criterion for successful phase referencing was to keep the phase errors below one radian. Simulations with ARIS reveal that phase referencing achieves good performance at 8.4GHz, even under poor tropospheric conditions. At 22 and 43GHz, it is recommended to conduct phase referencing observations under good or typical tropospheric conditions. The satellite is required to have an attitude-switching capability with a one-minute or shorter cycle, and an orbit determination accuracy higher than about 10cm at apogee; the phase referencing calibrators are required to have a signal-to-noise ratio larger than four for a single scan. The probability to find a suitable phase referencing calibrator was estimated by using VLBI surveys. From the viewpoint of calibrator availability, VSOP-2 phase referencing at 8.4GHz is promising. However, the change of finding suitable calibrators at 22 and 43GHz is significantly reduced; it is important to conduct specific investigations for each target at those frequencies.

We have measured the sub-milli-arcsecond structure of 274 extragalactic sources at 24 and 43 GHz in order to assess their astrometric suitability for use in a high frequency celestial reference frame (CRF). Ten sessions of observations with the Very Long Baseline Array have been conducted over the course of $\sim$5 years, with a total of 1339 images produced for the 274 sources. There are several quantities that can be used to characterize the impact of intrinsic source structure on astrometric observations including the source flux density, the flux density variability, the source structure index, the source compactness, and the compactness variability. A detailed analysis of these imaging quantities shows that (1) our selection of compact sources from 8.4 GHz catalogs yielded sources with flux densities, averaged over the sessions in which each source was observed, of about 1 Jy at both 24 and 43 GHz, (2) on average the source flux densities at 24 GHz varied by 20%-25% relative to their mean values, with variations in the session-to-session flux density scale being less than 10%, (3) sources were found to be more compact with less intrinsic structure at higher frequencies, and (4) variations of the core radio emission relative to the total flux density of the source are less than 8% on average at 24 GHz. We conclude that the reduction in the effects due to source structure gained by observing at higher frequencies will result in an improved CRF and a pool of high-quality fiducial reference points for use in spacecraft navigation over the next decade.

Since VLBI techniques give microarcsecond position accuracy of celestial objects, tests of GR using radio sources as probes of a gravitational field have been made. We present the results from two recent tests using the VLBA: In 2005, the measurement of the classical solar deflection; and in 2002, the measurement of the retarded gravitational deflection associated with Jupiter. The deflection experiment measured PPN-gamma to an accuracy of 0.0003; the Jupiter experiment measured the retarded term to 20% accuracy. The controversy over the interpretation of the retarded term is summarized.

We report on recent VLBA/VERA/IVS observational tests of General Relativity. First, we will summarize the results from the 2005 VLBA experiment that determined gamma with an accuracy of 0.0003 by measuring the deflection of four compact radio sources by the solar gravitational field. We discuss the limits of precision that can be obtained with VLBA experiments in the future. We describe recent experiments using the three global arrays to measure the aberration of gravity when Jupiter and Saturn passed within a few arcmin of bright radio sources. These reductions are still in progress, but the anticipated positional accuracy of the VLBA experiment may be about 0.01 mas.

We have used the Very Long Baseline Array (VLBA) at 43, 23 and 15 GHz to measure the solar gravitational deflection of radio waves among four radio sources during an 18-day period in October 2005. Using phase-referenced radio interferometry to fit the measured phase delay to the propagation equation of the parameterized post-Newtonian (PPN) formalism, we have determined the deflection parameter gamma = 0.9998 +/- 0.0003$ (68% confidence level), in agreement with General Relativity. The results come mainly from 43 GHz observations where the refraction effects of the solar corona were negligible beyond 3 degrees from the sun. The purpose of this experiment is three-fold: to improve on the previous results in the gravitational bending experiments near the solar limb; to examine and evaluate the accuracy limits of terrestrial VLBI techniques; and to determine the prospects and outcomes of future experiments. Our conclusion is that a series of improved designed experiments with the VLBA could increase the presented accuracy by at least a factor of 4.

Astrometry can bring powerful constraints to bear on a variety of scientific questions about neutron stars, including their origins, astrophysics, evolution, and environments. Using phase-referenced observations at the VLBA, in conjunction with pulsar gating and in-beam calibration, we have measured the parallaxes and proper motions for 14 pulsars. The smallest measured parallax in our sample is 0.13+-0.02 mas for PSR B1541+09, which has a most probable distance of 7.2+1.3-1.1 kpc. We detail our methods, including initial VLA surveys to select candidates and find in-beam calibrators, VLBA phase-referencing, pulsar gating, calibration, and data reduction. The use of the bootstrap method to estimate astrometric uncertainties in the presence of unmodeled systematic errors is also described. Based on our new model-independent estimates for distance and transverse velocity, we investigate the kinematics and birth sites of the pulsars and revisit models of the Galactic electron density distribution. We find that young pulsars are moving away from the Galactic plane, as expected, and that age estimates from kinematics and pulsar spindown are generally in agreement, with certain notable exceptions. Given its present trajectory, the pulsar B2045-16 was plausibly born in the open cluster NGC 6604. For several high-latitude pulsars, the NE2001 electron density model underestimates the parallax distances by a factor of two, while in others the estimates agree with or are larger than the parallax distances, suggesting that the interstellar medium is irregular on relevant length scales. The VLBA astrometric results for the recycled pulsar J1713+0747 are consistent with two independent estimates from pulse timing, enabling a consistency check between the different reference frames.

We present the multiwavelength properties of 266 cataloged radio sources identified with 20 and 6 cm VLA deep observations of the CDFS at a flux density limit of 42 \mu Jy at the field centre at 1.4 GHz. These new observations probe the faint end of both the star formation and radio galaxy/AGN population. X-ray data, including upper limits, turn out to be a key factor in establishing the nature of faint radio sources. We find that, while the well-known flattening of the radio number counts below 1 mJy is mostly due to star forming galaxies, these sources and AGN make up an approximately equal fraction of the sub--millijansky sky, contrary to some previous results. We have also uncovered a population of distant AGN systematically missing from many previous studies of sub-millijansky radio source identifications. The AGN include radio galaxies, mostly of the low-power, Fanaroff-Riley I type, and a significant radio-quiet component, which amounts to approximately one fifth of the total sample. We also find that radio detected, X-ray AGN are not more heavily obscured than the X-ray detected AGN. This argues against the use of radio surveys as an efficient way to search for the missing population of strongly absorbed AGN.

We discuss the X-ray properties of the radio sources detected in a deep 1.4 and 5 GHz VLA Radio survey of the Extended Chandra Deep Field South (E-CDFS). Among the 266 radio sources detected, we find 89 sources (1/3 of the total) with X-ray counterparts in the catalog of the 1Ms exposure of the central 0.08 deg^2 (Giacconi et al. 2002; Alexander et al. 2003) or in the catalog of the 250 ks exposure of the 0.3 deg^2 E-CDFS field (Lehmer et al. 2005). For 76 (85%) of these sources we have spectroscopic or photometric redshifts, and therefore we are able to derive their intrinsic properties from X-ray spectral analysis, namely intrinsic absorption and total X-ray luminosities. We find that the population of submillijansky radio sources with X-ray counterparts is composed of a mix of roughly 1/3 star forming galaxies and 2/3 AGN. The distribution of intrinsic absorption among X-ray detected radio sources is different from that of the X-ray selected sample. Namely, the fraction of low absorption sources is at least two times larger than that of X-ray selected sources in the CDFS. This is mostly due to the larger fraction of star forming galaxies present among the X-ray detected radio sources. If we investigate the distribution of intrinsic absorption among sources with L_X>10^42 erg s^-1 in the hard 2-10 keV band (therefore in the AGN luminosity regime), we find agreement between the X-ray population with and without radio emission. In general, radio detected X-ray AGN are not more heavily obscured than the non radio detected AGN. This argues against the use of radio surveys as an efficient way to search for the missing population of strongly absorbed AGN.

In February 1997, the Japanese radio astronomy satellite HALCA was launched to provide the space-bourne element for the VLBI Space Observatory Programme (VSOP) mission. A significant fraction of the mission time was to be dedicated to the VSOP Survey of bright compact Active Galactic Nuclei (AGN) at 5 GHz, which was lead by ISAS. The VSOP Survey Sources are an unbiased dataset of 294 targets, of which 82% were successfully observed. These are now undergoing statistical analysis to tease out the characteristics of typical AGN sources. We present here the summary of the imaging and conclusions we have reached.

We present a detailed analysis of 256 radio sources from our deep (flux density limit of 42 microJy at the field centre at 1.4 GHz) Chandra Deep Field South 1.4 and 5 GHz VLA survey. The radio population is studied by using a wealth of multi-wavelength information in the radio, optical, and X-ray bands. The availability of redshifts for ~ 80% of the sources in our complete sample allows us to derive reliable luminosity estimates for the majority of the objects. X-ray data, including upper limits, for all our sources turn out to be a key factor in establishing the nature of faint radio sources. Due to the faint optical levels probed by this study, we have uncovered a population of distant Active Galactic Nuclei (AGN) systematically missing from many previous studies of sub-millijansky radio source identifications. We find that, while the well-known flattening of the radio number counts below 1 mJy is mostly due to star forming galaxies, these sources and AGN make up an approximately equal fraction of the sub-millijansky sky, contrary to some previous results. The AGN include radio galaxies, mostly of the low-power, Fanaroff-Riley I type, and a significant radio-quiet component, which amounts to approximately one fifth of the total sample. The ratio of radio to optical luminosity depends more on radio luminosity, rather than being due to optical absorption.

This paper presents the sixth part to the Very Long Baseline Array (VLBA) Calibrator Survey. It contains the positions and maps of 264 sources of which 169 were not previously observed with very long baseline interferometry (VLBI). This survey, based on two 24 hour VLBA observing sessions, was focused on 1) improving positions of 95 sources from previous VLBA Calibrator surveys that were observed either with very large a priori position errors or were observed not long enough to get reliable positions and 2) observing remaining new flat-spectrum sources with predicted correlated flux density in the range 100-200 mJy that were not observed in previous surveys. Source positions were derived from astrometric analysis of group delays determined at the 2.3 and 8.6 GHz frequency bands using the Calc/Solve software package. The VCS6 catalogue of source positions, plots of correlated flux density versus projected baseline length, contour plots and fits files of naturally weighted CLEAN images, as well as calibrated visibility function files are available on the Web at http://vlbi.gsfc.nasa.gov/vcs6

We have obtained 8.4 GHz VLBA observations of a 31-GHz complete sample of ~100 sources between 10 and 100 mJy. The main goals of these observations are: to determine the angular size, radio spectra and identification for a weak sample of high frequency sources; to find the fraction of sources which have sufficiently compact emission for use as calibrators for VLBI observations; and for design considerations of the proposed DSN Array. We find that a large fraction of observed sources have VLBI detections. A majority of these sources have most of their emission in a compact <1 mas radio core, with remaining sources having steep radio spectra. The source list was provided from GBT observations to remove discrete sources in the CBI fields.

We report 20 and 6 cm VLA deep observations of the CDF-S including the Extended CDF-S. We discuss the radio properties of 266 cataloged radio sources, of which 198 are above a 20 cm completeness level reaching down to 43 microJy at the center of the field. Survey observations made at 6 cm over a more limited region covers the original CDF-S to a comparable level of sensitivity as the 20 cm observations. Of 266 cataloged radio sources, 52 have X-ray counterparts in the CDF-S and a further 37 in the E-CDF-S area not covered by the 1 Megasecond exposure. Using a wide range of material, we have found optical or infrared counterparts for 254 radio sources, of which 186 have either spectroscopic or photometric redshifts (Paper II). Three radio sources have no apparent counterpart at any other wavelength. Measurements of the 20 cm radio flux density at the position of each CDF-S X-ray source detected a further 30 radio sources above a conservative 3-sigma detection limit. X-ray and sub-mm observations have been traditionally used as a measure of AGN and star formation activity, respectively. These new observations probe the faint end of both the star formation and radio galaxy/AGN population, as well as the connection between the formation and evolution of stars and SMBHs. Both of the corresponding gravitational and nuclear fusion driven energy sources can lead to radio synchrotron emission. AGN and radio galaxies dominate at high flux densities. Although emission from star formation becomes more prominent at the microjansky levels reached by deep radio surveys, even for the weakest sources, we still find an apparent significant contribution from low luminosity AGN as well as from star formation.

We present the optical and infrared identifications of the 266 radio sources detected at 20 cm with the Very Large Array in the Chandra Deep Field South (Kellermann et al. 2008). Using deep i-band Advanced Camera for Surveys, R-band Wide Field Imager, K-band SOFI/NTT, K-band ISAAC/VLT and Spitzer imaging data, we are able to find reliable counterparts for 254 (~95%) VLA sources. Twelve radio sources remain unidentified and three of them are ``empty fields''. Using literature and our own data we are able to assign redshifts to 186 (~70%) radio sources: 108 are spectroscopic redshifts and 78 reliable photometric redshifts. Based on the rest frame colors and morphological distributions of the host galaxies we find evidences for a change in the submillijansky radio source population: a) above ~ 0.08 mJy early-type galaxies are dominating; b) at flux densities below ~0.08 mJy, starburst galaxies become dominant.

We have observed the Extended Chandra Deep Field South (E-CDF-S) using a mosaic of six deep Very Large Array (VLA) pointings at 1.4GHz. In this paper, we present the survey strategy, description of the observations, and the first data release. The observations were performed during June through September of 2007 and included from 15 to 17 ``classic'' VLA antennas and 6 to 11 that had been retrofitted for the Expanded VLA (EVLA). The first data release consists of a 34.1' x 34.1' image and the attendant source catalog. The image achieves an rms sensitivity of 6.4 uJy per 2.8" x 1.6" beam in its deepest regions, with a typical sensitivity of 8 uJy. The catalog is conservative in that it only lists sources with peak flux densities greater than seven times the local rms noise, yet it still contains 464 sources. Nineteen of these are complex sources consisting of multiple components. Cross matching of the catalog to prior surveys of the E-CDF-S confirms the linearity of the flux density calibration, albeit with a slight possible offset (a few percent) in scale. Improvements to the data reduction and source catalog are ongoing, and we intend to produce a second data release in January 2009.

In February 1997, the Japanese radio astronomy satellite HALCA was launched to provide the space-bourne element for the VLBI Space Observatory Programme (VSOP) mission. Approximately twenty-five percent of the mission time was dedicated to the VSOP Survey of bright compact Active Galactic Nuclei (AGN) at 5 GHz. This paper, the fifth in the series, presents images and models for the remaining 140 sources not included in Paper III, which contained 102 sources. For most sources, the plots of the uv-coverage, the visibility amplitude versus uv-distance, and the high resolution image are presented. Model fit parameters to the major radio components are determined, and the brightness temperature of the core component for each source is calculated. The brightness temperature distributions for all of the sources in the VSOP AGN survey are discussed.

Experimental verification of the existence of gravimagnetic fields generated by currents of matter is important for a complete understanding and formulation of gravitational physics. Although the rotational `intrinsic' gravimagnetic field has been extensively studied and is now being measured by the Gravity Probe B, the `extrinsic' gravimagnetic field generated by the translational current of matter is less well studied. The present paper uses the post-Newtonian parametrized Einstein and light geodesics equations to show that the `extrinsic' gravimagnetic field generated by the translational current of matter can be measured by observing the relativistic time delay and/or light deflection caused by the moving mass. We prove that the `extrinsic' gravimagnetic field is generated by the relativistic effect of the aberration of the gravity force caused by the Lorentz transformation of the metric tensor and the Levi-Civita connection. We show that the Lorentz transformation of the gravity field variables is equivalent to the technique of the retarded Lienard-Wiechert gravitational potentials predicting that a light particle is deflected by gravitational field of a moving body from its retarded position so that both general-relativistic phenomena -- the aberration and the retardation of gravity -- are tightly connected and observing the aberration of gravity proves that gravity has a causal nature. We explain in this framework the 2002 deflection experiment of a quasar by Jupiter where the aberration of gravity from its orbital motion was measured with accuracy 20%. We describe a theory of VLBI experiment to measure the gravitational deflection of radio waves from a quasar by the Sun, as viewed by a moving observer from the geocentric frame, to improve the measurement accuracy of the aberration of gravity to a few percent.

We analyse the 267 radio sources from our deep (flux limit of 42 microJy at the field center at 1.4 GHz) Chandra Deep Field South 1.4 and 5 GHz VLA survey. The radio population is studied by using a wealth of multi-wavelength information, including morphology and spectral types, in the radio, optical, and X-ray bands. The availability of redshifts for ~ 70% of our sources allows us to derive reliable luminosity estimates for the majority of the objects. Contrary to some previous results, we find that star-forming galaxies make up only a minority (~ 1/3) of sub-mJy sources, the bulk of which are faint radio galaxies, mostly of the Fanaroff-Riley I type.

In February 1997 the Japanese radio astronomy satellite HALCA was launched to provide the space-borne element for the VSOP mission. HALCA provided linear baselines three-times greater than that of ground arrays, thus providing higher resolution and higher AGN brightness temperature measurements and limits. Twenty-five percent of the scientific time of the mission was devoted to the ``VSOP survey'' of bright, compact, extra-galactic radio sources at 5 GHz. A complete list of 294 survey targets were selected from pre-launch surveys, 91% of which were observed during the satellite's lifetime. The major goals of the VSOP Survey are statistical in nature: to determine the brightness temperature and approximate structure, to provide a source list for use with future space VLBI missions, and to compare radio properties with other data throughout the electromagnetic spectrum. All the data collected have now been analysed and is being prepared for the final image Survey paper. In these papers we present details of the mission, images of the sources, and some statistics of the contributions and completeness plus some statistics of the images and brightness temperatures. Full resolution images are to be found on http://www.vsop.isas.jaxa.jp/survey/publications/evn.ps.gz

This paper presents the fifth part of the Very Long Baseline Array (VLBA) Calibrator Survey (VCS), containing 569 sources not observed previously with very long baseline interferometry in geodetic or absolute astrometry programs. This campaign has two goals: (i) to observe additional sources which, together with previous survey results, form a complete sample, (ii) to find new strong sources suitable as phase calibrators. This VCS extension was based on three 24-hour VLBA observing sessions in 2005. It detected almost all extragalactic flat-spectrum sources with correlated flux density greater than 200 mJy at 8.6 GHz above declination -30 degrees which were not observed previously. Source positions with milliarcsecond accuracy were derived from astrometric analysis of ionosphere-free combinations of group delays determined from the 2.3 GHz and 8.6 GHz frequency bands. The VCS5 catalog of source positions, plots of correlated flux density versus projected baseline length, contour plots and FITS files of naturally weighted CLEAN images, as well as calibrated visibility function files are available on the Web at http://vlbi.gsfc.nasa.gov/vcs5

We present a 1.4-GHz catalog of 810 radio sources (560 sources in the complete sample) with 1.8" resolution found within a 17' radius in the SSA13 field (RA=13h12m,DEC=42d38'). The radio image from the VLA has an rms noise level of 4.82 microJy/beam at the field center, and Subaru optical images in r-band (6300A) and z-band (9200A) have a three-sigma detection magnitude of 26.1 and 24.9, respectively. 88% of the radio sources are identified with an optical counterpart, and there is significantly more reddening for objects fainter than 24-mag. The radio and optical parameters are tabulated, and source morphologies are displayed by radio contours overlaying optical false-color images. The radio structures show a wealth of complexity and these are classified into a small number of categories. About one-third of the radio sources are larger than 1.2" and their orientation is often similar to that of the associated galaxy or binary-galaxy system. The density of sources in the SSA13 field above 75 microJy is 0.40 per square arcmin, with a slope of -2.43 in the differential counts. The radio spectral index may steepen for sources below 75 microJy. We estimate that at most 40% of the microJansky radio sources are dominated by AGN processes.

We describe our explicit Lorentz-invariant solution of the Einstein and null geodesic equations for the deflection experiment of 2002 September 8 when a massive moving body, Jupiter, passed within 3.7' of a line-of-sight to a distant quasar. We develop a general relativistic framework which shows that our measurement of the retarded position of a moving light-ray deflecting body (Jupiter) by making use of the gravitational time delay of quasar's radio wave is equivalent to comparison of the relativistic laws of the Lorentz transformation for gravity and light. Because, according to Einstein, the Lorentz transformation of gravity field variables must depend on a fundamental speed $c$, its measurement through the retarded position of Jupiter in the gravitational time delay allows us to study the causal nature of gravity and to set an upper limit on the speed of propagation of gravity in the near zone of the solar system as contrasted to the speed of the radio waves. We discuss the misconceptions which have inhibited the acceptance of this interpretation of the experiment. We also comment on other interpretations of this experiment by Asada, Will, Samuel, Pascual-Sanchez, and Carlip and show that their `speed of light' interpretations confuse the Lorentz transformation for gravity with that for light, and the fundamental speed of gravity with the physical speed of light from the quasar.

Recent papers by Samuel declared that the linearized post-Newtonian v/c effects are too small to have been measured in the recent experiment involving Jupiter and quasar J0842+1845 that was used to measure the ultimate speed of gravity defined as a fundamental constant entering in front of each time derivative of the metric tensor in the Einstein gravity field equations. We describe our Lorentz-invariant formulation of the Jovian deflection experiment and confirm that v/c effects are do observed, as contrasted to the erroneous claim by Samuel, and that they vanish if and only if the speed of gravity is infinite.

This paper presents the fourth extension to the Very Long Baseline Array (VLBA) Calibrator Survey, containing 258 new sources not previously observed with very long baseline interferometry (VLBI). This survey, based on three 24 hour VLBA observing sessions, fills remaining areas on the sky above declination -40 degrees where the calibrator density is less than one source within a 4 degree radius disk at any given direction. The share of these area was reduced from 4.6% to 1.9%. Source positions were derived from astrometric analysis of group delays determined at 2.3 and 8.6 GHz frequency bands using the Calc/Solve software package. The VCS4 catalogue of source positions, plots of correlated flux density versus projected baseline length, contour plots and fits files of naturally weighted CLEAN images, as well as calibrated visibility function files are available on the Web at http://gemini.gsfc.nasa.gov/vcs4 .

The highest velocity neutron stars establish stringent constraints on natal kicks, asymmetries in supernova core collapse, and the evolution of close binary systems. Here we present the first results of a long-term pulsar astrometry program using the VLBA. We measure a proper motion and parallax for the pulsar B1508+55, leading to model-independent estimates of its distance (2.37+0.23-0.20 kpc) and transverse velocity (1083+103-90 km/s), the highest velocity directly measured for a neutron star. We trace the pulsar back from its present Galactic latitude of 52.3 degrees to a birth site in the Galactic plane near the Cyg OB associations, and find that it will inevitably escape the Galaxy. Binary disruption alone is insufficient to impart the required birth velocity, and a natal kick is indicated. A composite scenario including a large kick along with binary disruption can plausibly account for the high velocity.

Determination of pulsar parallaxes and proper motions addresses fundamental astrophysical questions. We have recently finished a VLBI astrometry project to determine the proper motions and parallaxes of 27 pulsars, thereby doubling the total number of pulsar parallaxes. Here we summarise our astrometric technique and present the discovery of a pulsar moving in excess of 1000 km/s. As an example of the application of high precision pulsar astrometry we also infer the identification of 2 pulsars originating from a disrupted binary in the Cygnus Superbubble.

We present full polarization Very Long Baseline Array (VLBA) observations at 5 GHz and 15 GHz of 24 compact active galactic nuclei (AGN). These sources were observed as part of a pilot project to demonstrate the feasibility of conducting a large VLBI survey to further our understanding of the physical properties and temporal evolution of AGN jets. The sample is drawn from the Cosmic Lens All-Sky Survey (CLASS) where it overlaps with the Sloan Digital Sky Survey at declinations north of 15 degrees. There are 2100 CLASS sources brighter than 50 mJy at 8.4 GHz, of which we have chosen 24 for this pilot study. All 24 sources were detected and imaged at 5 GHz with a typical dynamic range of 500:1, and 21 of 24 sources were detected and imaged at 15 GHz. Linear polarization was detected in 8 sources at both 5 and 15 GHz, allowing for the creation of Faraday rotation measure (RM) images. The core RMs for the sample were found to have an average absolute value of 390 +/- 100 rad/m^2. We also present the discovery of a new Compact Symmetric Object, J08553+5751. All data were processed automatically using pipelines created or adapted for the survey.

Eighteen days of MERLIN data and 42 hours of A-array VLA data at 1.4 GHz have been combined to image a 10-arcmin field centred on the Hubble Deep and Flanking Fields (HDF and HFF). A complete sample of 92 radio sources with 1.4-GHz flux densities above 40 microJy has been imaged using MERLIN+VLA. The images are amongst the most sensitive yet made at 1.4 GHz, with rms noise levels of 3.3 microJy/beam in the 0.2-arcsec images. Virtually all the sources are resolved, with angular sizes in the range 0.2 to 3 arcsec. No additional sources were detected down to 23 microJy in the central 3 arcmin, indicating that sources fainter than 40 microJy are heavily resolved with MERLIN and must have typical angular sizes greater than 0.5 arcsec. Compact radio sources were used to align the optical data to the ICRF, to <50 mas in the HDF. We find a statistical association of very faint (2 microJy and above) radio sources with optically bright HDF galaxies down to about 23 mag. Of the 92 radio sources above 40 microJy, about 85 percent are identified with galaxies brighter than about I = 25 mag; the remaining 15 percent are associated with optically faint systems. We identify several very red, optically faint systems including the the strongest sub-mm source in the HDF, HDF850.1. 72 percent of the radio sources are starburst or AGN-type systems; the remainder are unclassified. The proportion of starburst systems increases with decreasing flux density; below 100 microJy 70 percent of the sources are starburst-type systems in the redshift range 0.3 -- 1.3. Chandra detections are associated with 55 of the 92 radio sources but their X-ray flux densities do not appear to be correlated with the radio flux densities or morphologies.

This paper presents the third extension to the Very Large Baseline Array (VLBA) Calibrator Survey, containing 360 new sources not previously observed with very long baseline interferometry (VLBI). The survey, based on three 24 hour VLBA observing sessions, fills the areas on the sky above declination -45 degrees where the calibrator density is less than one source within a 4 degrees radius disk at any given direction. The positions were derived from astrometric analysis of the group delays determined at 2.3 and 8.6 GHz frequency bands using the Calc/Solve software package. The VCS3 catalogue of source positions, plots of correlated flux density versus the length of projected baseline, contour plots and fits files of naturally weighted CLEAN images as well as calibrated visibility function files are available on the Web at http://gemini.gsfc.nasa.gov/vcs3

The sensitivity and versatility of SKA will provide microarcsec astrometric precision and high quality milliarcsec-resolution images by simultaneously detecting calibrator sources near the target source. To reach these goals, we suggest that the long-baseline component of SKA contains at least 25% of the total collecting area in a region between 1000 to 5000 km from the core SKA. We also suggest a minimum of 60 elements in the long-baseline component of SKA to provide the necessary (u-v) coverage. For simultaneous all-sky observations, which provide absolute astrometric and geodetic parameters, we suggest using ten independent subarrays each composed of at least six long-baseline elements correlated with the core SKA. We discuss many anticipated SKA long-baseline astrometric experiments: determination of distance, proper motion and orbital motion of thousands of stellar objects; planetary motion detections; mass determination of degenerate stars using their kinetics; calibration of the universal distance scale from 10 pc to 10 Mpc; the core and inner-jet interactions of AGN. With an increase by a factor of 10 in absolute astrometric accuracy using simultaneous all sky observations, the fundamental quasar reference frame can be defined to <10 microarcsec and tied to the solar-system dynamic frame to this accuracy. Parameters associated with the earth rotation and orientation, nutation, and geophysical parameters, can be accurately monitored. Tests of fundamental physics include: solar and Jovian deflection experiments, the sky frame accuracy needed to interpret the gravity wave/pulsar-timing experiment, accurate monitoring of spacecraft orbits which impacts solar system dynamics.

The VSOP (VLBI Space Observatory Programme) mission is a Japanese-led project to study radio sources with sub-milliarcsec angular resolution, using an orbiting 8-m telescope on board the satellite HALCA with a global earth-based array of telescopes. A major program is the 5 GHz VSOP Survey Program, which we supplement here with VLBA observations to produce a complete and flux-density limited sample. Using statistical methods of analysis of the observed visibility amplitude versus projected (u,v) spacing, we have determined the angular size and brightness temperature distribution of bright AGN radio cores. On average, the cores have a diameter (full-width, half-power) of 0.20 mas which contains about 20% of the total source emission, and (14+/-6)% of the cores are <0.04 mas in size. About (20+/-5)% of the radio cores have a source frame brightness temperature 10^{13}K, and (3+/-2)% have 10^{14}K. A model of the high brightness temperature tail suggests that the radio cores have a brightness temperatures approx 10^{12}K, and are beamed toward the observer with an average bulk motion of v/c=0.993 +/- 0.004.

The VSOP mission is a Japanese-led project to study radio sources with sub-milliarcsecond angular resolution using an orbiting 8-m telescope, HALCA and global arrays of Earth-based telescopes. Approximately 25% of the observing time has been devoted to a survey of compact AGN at 5 GHz which are stronger than 1 Jy -- the VSOP AGN Survey. This paper, the second in a series, describes the data calibration, source detection, self-calibration, imaging and modeling, and gives examples illustrating the problems specific to space VLBI. The VSOP Survey web-site which contains all results and calibrated data is described.

The VSOP mission is a Japanese-led project to study radio sources with sub-milliarcsec resolution using an orbiting 8 m telescope, HALCA, along with global arrays of Earth-based telescopes. Approximately 25% of the observing time is devoted to a survey of compact AGN which are stronger than 1 Jy at 5 GHz-the VSOP AGN Survey. This paper, the third in the series, presents the results from the analysis of the first 102 Survey sources. We present high resolution images and plots of visibility amplitude versus projected baseline length. In addition, model-fit parameters to the primary radio components are listed, and from these the angular size and brightness temperature for the radio cores are calculated. For those sources for which we were able to determine the source frame core brightness temperature, a significant fraction (53 out of 98) have a source frame core brightness temperature in excess of 10^12 K. The maximum source frame core brightness temperature we observed was 1.2 X 10^13 K. Explaining a brightness temperature this high requires an extreme amount of relativistic Doppler beaming. Since the maximum brightness temperature one is able to determine using only ground-based arrays is of the order of 10^12 K, our results confirm the necessity of using space VLBI to explore the extremely high brightness temperature regime.

Reference systems and frames are crucial for high precision absolute astrometric work, and their foundations must be well-defined. The current frame, the International Celestial Reference Frame, will be discussed: its history, the use of the group delay as the measured quantity, the positional accuracy of 0.3 milliarcsec, and possible future improvements. On the other hand, for the determination of the motion of celestial objects, accuracies approaching 0.01 milliarcsec can be obtained by measuring the differential position between the target object and nearby stationary sources. This astrometric technique uses phase referencing, and the current techniques and limitations are discussed, using the results from four experiments. Brief comments are included on the interpretation of the Jupiter gravity deflection experiment of September 2002.

In February 1997 the Japanese radio astronomy satellite Halca was launched to provide the space-bourne element for the VSOP mission. Approximately twenty-five percent of the mission time has been dedicated to the VSOP Survey, a 5 GHz survey of bright, compact AGN. Both the final, calibrated, high resolution images and plots of visibility amplitude versus uv distance for the first 102 of the sources have been prepared and has been submitted. Papers on the methods and the models from fitting the cumulative {\em uv} amplitudes will also be submitted. The analysis of the second half is well underway.

We use redshift observations of two deep 1.4GHz fields to probe the evolution of the bright end of the radio galaxy luminosity function to z=1.5. We show that the number of galaxies with radio power that would correspond locally to an ultraluminous infrared galaxy (ULIG) evolves rapidly over this redshift range. The optical spectra and X-ray luminosities are used to refine this result by separating the sources with clear active galactic nucleus (AGN) signatures from those that may be dominated by star formation. Both populations show extremely rapid evolution over this redshift range. We find that the number of sources with ULIG radio power and no clear AGN signatures evolves as (1+z)^(7).

Deep VLT optical spectroscopy, HST+ACS (GOODS) imaging and VLA observations are used to unveil the nature of a complete sample of 47 EROs with R-Ks>5 and Ks<20. The spectroscopic redshift completeness is 62%. Morphological classification was derived for each ERO through visual inspection and surface brightness profile fitting. Three main ERO morphological types are found: E/S0 galaxies (30-37%), spiral-like (24-46%) and irregular systems (17-39%). The only ERO detected in the radio is likely to host an obscured AGN. The average radio luminosity of the star-forming EROs undetected in the radio implies star formation rates of the order of 33 Msun yr^-1. The colors, redshifts and masses of the E/S0 galaxy subsample imply a minimum formation redshift z_f \sim 2. With this z_f there is enough time to have old and massive stellar spheroids already assembled at z \sim 1. We verify that the R-Ks vs. J-Ks color diagram is efficient in segregating old and dusty-star-forming EROs.

Positions and maps of 1608 new compact sources were obtained in twelve sessions observed during 1994--2002 at the VLBA network at 8.4/2.3 GHz. These sources are recommended for use as calibrators for phase reference imaging and as geodetic sources for astrometric/geodetic VLBI applications.

We have determined the relativistic light deflection of the quasar J0842+1835 as Jupiter passed within 3.7' on 2002 September 8, by measuring the time delay using the VLBA and Effelsberg radio telescopes at 8.4 GHz. At closest approach, General Relativity (GR) predicts a radial (static) deflection of 1190 microarcsec, and tangential (retarded) deflection in the direction of Jupiter's motion of 51 microarcsec. Our experiment achieved an rms position error of <10 microarcsec, and measured this retarded deflection to be 0.98 +/- 0.19 (rms error) times that predicted by GR. The increase positional accuracy for this VLBI phase referencing experiment was achieved by using two calibrator sources. Comments on the interpretation of this experiment are given.

(abridged) We present multiwavelength observations for a large sample of microjansky radio sources detected in ultradeep 1.4GHz maps centered on the Hubble Deep Field-North (HDF-N) and the Hawaii Survey Fields SSA13 and SSA22. Our spectroscopic redshifts for 169 radio sources reveal a flat median redshift distribution, and these sources are hosted by similarly luminous optical L* galaxies, regardless of redshift. This is a serious concern for radio estimates of the local star formation rate density, as a substantial fraction of the ultraviolet luminosity density is generated by sub-L* galaxies at low redshifts. From our submillimeter measurements for 278 radio sources, we find error-weighted mean 850micron fluxes of 1.72$\pm$0.09 mJy for the total sample, 2.37$\pm$0.13 mJy for the optically-faint (I>23.5) subsample, and 1.04$\pm$0.13 mJy for the optically-bright (I<23.5) subsample. We significantly (>3\sigma) detect in the submillimeter 50 of the radio sources, 38 with I>23.5. Spectroscopic redshifts for three of the I<23.5 submillimeter-detected radio sources are in the range z=1.0-3.4, and all show AGN signatures. Using only the submillimeter mapped regions we find that 69\pm9% of the submillimeter-detected radio population are at I>23.5. We also find that 66\pm7% of the S850>5 mJy (>4\sigma) sources are radio-identified. We find that millimetric redshift estimates at low redshifts are best made with a FIR template intermediate between a Milky Way type galaxy and a starburst galaxy, and at high redshifts with an Arp220 template.

A relativistic sub-picosecond model of gravitational time delay in radio astronomical observations is worked out and a new experimental test of general relativity is discussed in which the effect of retardation of gravity associated with its finite speed can be observed. As a consequence, the speed of gravity can be measured by differential VLBI observations. Retardation in propagation of gravity is a central part of the Einstein theory of general relativity which has not been tested directly so far. The idea of the proposed gravitational experiment is based on the fact that gravity in general relativity propagates with finite speed so that the deflection of light caused by the body must be sensitive to the ratio of the body's velocity to the speed of gravity. The interferometric experiment can be performed, for example, during the very close angular passage of a quasar by Jupiter. Due to the finite speed of gravity and orbital motion of Jupiter, the variation in its gravitational field reaches observer on Earth not instantaneously but at the retarded instant of time and should appear as a velocity-dependent excess time delay in addition to the well-known Shapiro delay, caused by the static part of the Jupiter's gravitational field. Such Jupiter-QSO encounter events take place once in a decade. The next such event will occur on September 8, 2002 when Jupiter will pass by quasar J0842+1835 at the angular distance 3.7 arcminutes. If radio interferometric measurement of the quasar coordinates in the sky are done with the precision of a few picoseconds (about 5 microarcseconds) the effect of retardation of gravity and its speed of propagation may be measured with an accuracy about 10%.

We present the results from two radio integrations at 8.4 GHz using the VLA. One of the fields, at 13h,+43d (SA13 field), has an rms noise level of 1.49 microJy and is the deepest radio image yet made. Thirty-four sources in a complete sample were detected above 7.5 microJy and 25 are optically identified to a limit of I=25.8, using our deep HST and ground-based images. The radio sources are usually located within 0.5" (typically 5 kpc) of a galaxy nucleus, and generally have a diameter less than 2.5". The second field at 17h, +50d (Hercules Field) has an rms noise of 35 microJy and contains 10 sources. We have also analyzed a complete flux density-limited sample at 8.4 GHz of 89 sources from five deep radio surveys, including the Hubble deep field. Half of all the optical counterparts are with galaxies brighter than I=23 mag, but 20% are fainter than I=25.5 mag. We confirm the tendency for the micro-Jansky radio sources to prefer multi-galaxy systems. The distribution of the radio spectral index between 1.4 and 8.4 GHz peaks at alpha = -0.75~ with a median value of -0.6. The average spectral index becomes steeper (lower values) for sources below 35 microJy, and for sources identified with optical counterparts fainter than I=25.5 mag. The differential radio count between 7.5 and 1000 microJy has a slope of -2.11 +/-0.13 and a surface density of 0.64 sources per square-arcmin with flux density greater than $7.5 microJy.

A catalog containing milliarcsecond--accurate positions of 1332 extragalactic radio sources distributed over the northern sky is presented - the Very Long Baseline Array Calibrator Survey (VCS1). The positions have been derived from astrometric analysis of dual-frequency 2.3 and 8.4 GHz VLBA snapshot observations; in a majority of cases, images of the sources are also available. These radio sources are suitable for use in geodetic and astrometric experiments, and as phase-reference calibrators in high-sensitivity astronomical imaging. The VCS1 is the largest high-resolution radio survey ever undertaken, and triples the number of sources available to the radio astronomy community for VLBI applications. In addition to the astrometric role, this survey can be used in active galactic nuclei, Galactic, gravitational lens and cosmological studies. The VCS1 catalog will soon be available at www.nrao.edu/vlba/VCS1 .

We present the results of a multi-frequency study of the inversion symmetric radio galaxy NGC326 based on Very Large Array observations at 1.4, 1.6, 4.8, 8.5 and 14.9 GHz. The morphological, spectral and polarization properties of this peculiar object are studied at different levels of spatial resolutions. The interpretation of the data will be discussed in forthcoming papers.

W75N is a star-forming region containing various ultracompact HII regions and OH, water, and methanol maser emission. Our VLBA map shows that the OH masers are located in a thin disk rotating around an O-star which is the exciting star of the ultracompact HII region VLA1. A separate set of maser spots is connected with the ultracompact HII region VLA2. The radial velocity of OH maser spots varies across the disk from 3.7 km/s to 10.9 km/s. The diameter of the disk is 4000 A.U. All maser spots are strongly polarized. This are the first OH masers showing nearly 100 per cent linear polarization in several spots. Two maser spots seem to be Zeeman pairs corresponding to a magnetic field of 5.2 mgauss and 7.7 mgauss, and in one case we tentatively found a Zeeman pair consisting of two linearly polarized components. The linearly polarized maser spots are shown to be sigma-components which is the case of the magnetic field being perpendicular to the line of sight. The direction of the magnetic field as determined from linearly polarized spots is perpendicular to the plane of the disk, although the galactic Faraday rotation may significantly affect this conclusion.

We present a comprehensive classification of all observed QPOs within the framework of the transition layer model using a large set of Rossi X-ray Timing Explorer (RXTE) data for Sco X-1. The model assumes an optically thin material along the observer's line of sight in the horizontal branch and an increasingly optically thick material while in the other two branches that is consistent with X-ray and radio observations and the disk transition layer model of QPOs. We identify the ~ 6 Hz frequencies in the normal branch as acoustic oscillations of a spherical shell around the neutron star (NS) that is formed after radiation pressure near the Eddington accretion rate destroys the disk. The size of the shell is on the order of one NS radii from the NS. We also estimate the upper limit of Sco X-1's magnetic field to be 0.7 x 10^6 G at about one NS radii above the NS surface while in the horizontal X-ray branch.

The radio components associated with the LMXB Sco X-1 have been monitored with extensive VLBI imaging at 1.7 and 5.0 GHz over four years, including a 56-hour continuous VLBI observation in 1999 June. We often detected one strong and one weak compact radio component, moving in opposite directions from the radio core. We suggest that the moving components are radio lobes generated by the disruption of energy flow in a twin-beam from the binary system. The average lifetime of a lobe-pair, the space motion of the lobes and the measured energy flow in the beam are discussed in arXiv:astro-ph/0104325. The lobe has a flux density that is variable over a time-scale of one hour, a measured minimum size of 1 mas (2.8 au), and is extended perpendicular to its motion. This short electron radiative lifetime may be caused by synchrotron losses if the lobe magnetic field is 300 G, or by adiabatic expansion of the electrons as soon as they are produced at the working surface. The lobes also show periods of slow expansion and a steepening radio spectrum, perhaps related to the characteristics of the beam energy flow. The radio morphology for Sco X-1 is more simple than for most other Galactic jet sources. The lobes of Sco X-1 are similar to hot-spots found in many extragalactic double sources. Scaling the phenomena observed in Sco X-1 to extragalactic sources implies hot-spot variability time-scales of 10^4 yr and hot-spot lifetimes of 10^5 yr. The recurrent formation of lobes in Sco X-1 probably does not occur for extragalactic radio sources.

The evolution of the radio emission from Sco X-1 is determined from a 56-hour continuous VLBI observation and from shorter observations over a four-year period. The radio source consists of a variable core near the binary, and two variable compact radio lobes which form near the core, move diametrically outward, then fade away. Subsequently, a new lobe-pair form near the core and the behavior repeats. The differences in the radio properties of the two lobes are consistent with the delay and Doppler-boosting associated with an average space velocity of 0.45c at 44 deg to the line of sight. Four lobe speeds, between 0.32c and 0.57c, were measured for several lobe-pairs on different days. The speed during each epoch remained constant over many hours. The direction of motion of the lobes over all epochs remained constant to a few degrees. Two core flares are contemporaneous with two lobe flares after removal of the delay associated with an energy burst moving with speed >0.95c in a twin-beam from the core to each lobe. This is the first direct measurement of the speed of energy flow within an astrophysical jet. The similarity of the core and lobe flares suggests that the twin-beam flow is symmetric and that the core is located near the base of the beam.

We present preliminary results of a multi-frequency study of the inversion symmetric radio galaxy NGC326 based on VLA observations at 1.4, 1.6, 4.8, 8.5, and 14.9 GHz. These data allow us to investigate in detail the morphological, spectral and polarization properties of this peculiar object at different levels of spatial resolution.

Results are presented from a long-term astrometry program on PSR B0919+06 using the NRAO Very Long Baseline Array. With ten observations (seven epochs) between 1994--2000, we measure a proper motion of 18.35 +/- 0.06 mas/yr in RA, 86.56 +/- 0.12 mas/yr in Dec, and a parallax of 0.83 +/- 0.13 mas (68% confidence intervals). This yields a pulsar distance of 1.21 +/- 0.19 kpc, making PSR B0919+06 the farthest pulsar for which a trigonometric parallax has been obtained, and the implied pulsar transverse speed is 505 +/- 80 km/s. Combining the distance estimate with interstellar scintillation data spanning 20 years, we infer the existence of a patchy or clumpy scattering screen along the line of sight in addition to the distributed electron density predicted by models for the Galaxy, and constrain the location of this scattering region to within about 250 parsecs of the Sun. Comparison with the lines of sight towards other pulsars in the same quadrant of the Galaxy permits refinement of our knowledge of the local interstellar matter in this direction.

We report on the first space-VLBI observations of the OH34.26+0.15 maser in two main line OH transitions at 1665 and 1667 MHz. The observations involved the space radiotelescope on board the Japanese satellite HALCA and an array of ground radio telescopes. The map of the maser region and images of individual maser spots were produced with an angular resolution of 1 milliarcsec which is several times higher than the angular resolution available on the ground. The maser spots were only partly resolved and a lower limit to the brightness temperature 6x10^{12} K was obtained. The maser seems to be located in the direction of low interstellar scattering, an order of magnitude lower than the scattering of a nearby extragalactic source and pulsar.

A new technique has been developed to remove the ionosphere's distorting effects from low frequency VLBI data. By fitting dispersive and non-dispersive components to the phases of multi-frequency data, the ionosphere can be effectively removed from the data without the use of {\em a priori} calibration information. This technique, along with the new gating capability of the VLBA correlator, was used to perform accurate astrometry on pulsar B0950+08, resulting in a much improved measurement of this pulsar's proper motion ($\mu_{\alpha} = -1.6 \pm 0.4$ mas/yr, $\mu_{\delta} = 29.5 \pm 0.5$ mas/yr) and parallax ($\pi = 3.6 \pm 0.3$ mas). This puts the pulsar at a distance of $280 \pm 25$ parsecs, about twice as far as previous estimates, but in good agreement with models of the electron density in the local bubble.

We report on the identifications of radio sources from our survey of the Hubble Deep Field and the SSA13 fields, both of which comprise the deepest radio surveys to date at 1.4 GHz and 8.5 GHz respectively. About 80% of the microjansky radio sources are associated with moderate redshift starburst galaxies or AGNs within the I magnitude range of 17 to 24 with a median of I = 22 mag. Thirty-one (20%) of the radio sources are: 1) fainter than $I>$25 mag, with two objects in the HDF $I_{AB}>$28.5, 2) often identified with very red objects $I-K>$4, and 3) not significantly different in radio properties than the brighter objects. We suggest that most of these objects are associated with heavily obscured starburst galaxies with redshifts between 1 and 3. However, other mechanisms are discussed and cannot be ruled out with the present observations.

We present Very Long Baseline Array phase-referenced measurements of the parallax and proper motion of two pulsars, B0919+06 and B1857-26. Sub-milliarcsecond positional accuracy was obtained by simultaneously observing a weak calibrator source within the 40' field of view of the VLBA at 1.5 GHz. We discuss the merits of using weak close calibrator sources for VLBI observations at low frequencies, and outline a method of observation and data reduction for these type of measurements. For the pulsar B1919+06 we measure a parallax of 0.31 +/- 0.14 mas. The accuracy of the proper motions is approximately 0.5 mas, an order of magnitude improvement over most previous determinations.

We report on sensitive radio observations made with the VLA at 8.5 GHz, centered on the Hubble Deep Field (HDF). We collected data in the A, CnB, C, DnC, and D configurations, corresponding to angular resolutions ranging from 0.3" to 10". We detected 29 radio sources in a complete sample within 4.6' of the HDF center and above a flux density limit of 9.0 microjy (5 sigma). Seven of these sources are located within the HDF itself, while the remaining 22 sources are covered by the Hubble Flanking Fields (HFFs) or ground based optical images. All of the sources in the HDF are identified with galaxies with a mean magnitude R = 21.7, while the mean magnitude of the identifications outside the HDF is R = 22.1. Three radio sources have no optical counterparts to R = 27. Based on a radio and optical positional coincidence, we detected an additional 19 radio sources in this field (seven of which are contained in the HDF) with 6.3 microjy < S < 9.0 (3.5 sigma < S < 5 sigma) and and R < 25, but which are not included in the complete sample. At least 2/3 of the optical identifications of these microjansky radio sources are with disk systems (i.e., irregulars, mergers and spirals). We briefly discuss the possible physical mechnisms driving the strong evolution in the microjansky radio galaxy population. [truncated]

We report on the detection of a feature with negative radio flux in a sensitive, low resolution observation with the VLA. This morphologically peculiar feature is approximately 25" 65" in size with a peak central amplitude of about -0.25 mK compared to the brightness temperature of the background sky at 3.6 cm. Within about 1' of this microwave decrement, we also found two radio quiet quasars, both at 2.561, with a projected physical separation of about 1 Mpc ($q_o$ = 0.5, h = 0.5), suggestive of a possible galaxy cluster in the region. We are unable to account for this negative source by systematic phase fluctuations within our observations, sidelobe artifacts, or instrumental noise. We discuss possible physical origins of this microwave source, in particular the inverse-Compton scattering of the cosmic microwave background by hot gas in a distant cluster of galaxies.