This paper presents the results from one of the first observations of ionospheric scintillation taken using the Low-Frequency Array (LOFAR). The observation was of the strong natural radio source Cas A, taken overnight on 18-19 August 2013, and exhibited moderately strong scattering effects in dynamic spectra of intensity received across an observing bandwidth of 10-80MHz. Delay-Doppler spectra (the 2-D FFT of the dynamic spectrum) from the first hour of observation showed two discrete parabolic arcs, one with a steep curvature and the other shallow, which can be used to provide estimates of the distance to, and velocity of, the scattering plasma. A cross-correlation analysis of data received by the dense array of stations in the LOFAR "core" reveals two different velocities in the scintillation pattern: a primary velocity of ~30m/s with a north-west to south-east direction, associated with the steep parabolic arc and a scattering altitude in the F-region or higher, and a secondary velocity of ~110m/s with a north-east to south-west direction, associated with the shallow arc and a scattering altitude in the D-region. Geomagnetic activity was low in the mid-latitudes at the time, but a weak sub-storm at high latitudes reached its peak at the start of the observation. An analysis of Global Navigation Satellite Systems (GNSS) and ionosonde data from the time reveals a larger-scale travelling ionospheric disturbance (TID), possibly the result of the high-latitude activity, travelling in the north-west to south-east direction, and, simultaneously, a smaller--scale TID travelling in a north-east to south-west direction, which could be associated with atmospheric gravity wave activity. The LOFAR observation shows scattering from both TIDs, at different altitudes and propagating in different directions. To the best of our knowledge this is the first time that such a phenomenon has been reported.

The largest Key Science Program of the RadioAstron space VLBI mission is a survey of active galactic nuclei (AGN). The main goal of the survey is to measure and study the brightness of AGN cores in order to better understand the physics of their emission while taking interstellar scattering into consideration. In this paper we present detection statistics for observations on ground-space baselines of a complete sample of radio-strong AGN at the wavelengths of 18, 6, and 1.3 cm. Two-thirds of them are indeed detected by RadioAstron and are found to contain extremely compact, tens to hundreds of $\mu$as structures within their cores.

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45$^\circ$00$'$00$''$ to 57$^\circ$00$'$00$''$) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325,694 sources are detected with a signal of at least five times the noise, and the source density is a factor of $\sim 10$ higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S$_{\rm 144 MHz} = 71\,\mu$Jy beam$^{-1}$ and the point-source completeness is 90% at an integrated flux density of 0.45mJy. The resolution of the images is 6$''$ and the positional accuracy is within 0.2$''$. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further $\sim$20 articles that highlight the scientific potential of LoTSS.

Type II radio bursts are evidence of shocks in the solar atmosphere and inner heliosphere that emit radio waves ranging from sub-meter to kilometer lengths. These shocks may be associated with CMEs and reach speeds higher than the local magnetosonic speed. Radio imaging of decameter wavelengths (20-90 MHz) is now possible with LOFAR, opening a new radio window in which to study coronal shocks that leave the inner solar corona and enter the interplanetary medium and to understand their association with CMEs. To this end, we study a coronal shock associated with a CME and type II radio burst to determine the locations at which the radio emission is generated, and we investigate the origin of the band-splitting phenomenon.

Context. The Sun is an active star that produces large-scale energetic events such as solar flares and coronal mass ejections and numerous smaller-scale events such as solar jets. These events are often associated with accelerated particles that can cause emission at radio wavelengths. The reconfiguration of the solar magnetic field in the corona is believed to be the cause of the majority of solar energetic events and accelerated particles. Aims. Here, we investigate a bright J-burst that was associated with a solar jet and the possible emission mechanism causing these two phenomena. Methods. We used data from the Solar Dynamics Observatory (SDO) to observe a solar jet, and radio data from the Low Frequency Array (LOFAR) and the Nan\c{c}ay Radioheliograph (NRH) to observe a J-burst over a broad frequency range (33-173 MHz) on 9 July 2013 at ~11:06 UT. Results. The J-burst showed fundamental and harmonic components and it was associated with a solar jet observed at extreme ultraviolet wavelengths with SDO. The solar jet occurred at a time and location coincident with the radio burst, in the northern hemisphere, and not inside a group of complex active regions in the southern hemisphere. The jet occurred in the negative polarity region of an area of bipolar plage. Newly emerged positive flux in this region appeared to be the trigger of the jet. Conclusions. Magnetic reconnection between the overlying coronal field lines and the newly emerged positive field lines is most likely the cause of the solar jet. Radio imaging provides a clear association between the jet and the J-burst which shows the path of the accelerated electrons.

(abridged). We outline LBCS (the LOFAR Long-Baseline Calibrator Survey), whose aim is to identify sources suitable for calibrating the highest-resolution observations made with the International LOFAR Telescope, which include baselines >1000 km. Suitable sources must contain significant correlated flux density (50-100mJy) at frequencies around 110--190~MHz on scales of a few hundred mas. At least for the 200--300-km international baselines, we find around 1 suitable calibrator source per square degree over a large part of the northern sky, in agreement with previous work. This should allow a randomly selected target to be successfully phase calibrated on the international baselines in over 50% of cases. Products of the survey include calibrator source lists and fringe-rate and delay maps of wide areas -- typically a few degrees -- around each source. The density of sources with significant correlated flux declines noticeably with baseline length over the range 200--600~km, with good calibrators on the longest baselines appearing only at the rate of 0.5 per square degree. Coherence times decrease from 1--3 minutes on 200-km baselines to about 1 minute on 600-km baselines, suggesting that ionospheric phase variations contain components with scales of a few hundred kilometres. The longest median coherence time, at just over 3 minutes, is seen on the DE609 baseline, which at 227km is close to being the shortest. We see median coherence times of between 80 and 110 seconds on the four longest baselines (580--600~km), and about 2 minutes for the other baselines. The success of phase transfer from calibrator to target is shown to be influenced by distance, in a manner that suggests a coherence patch at 150-MHz of the order of 1 degree.

Cosmic rays are the highest energy particles found in nature. Measurements of the mass composition of cosmic rays between 10^{17} eV and 10^{18} eV are essential to understand whether this energy range is dominated by Galactic or extragalactic sources. It has also been proposed that the astrophysical neutrino signal comes from accelerators capable of producing cosmic rays of these energies. Cosmic rays initiate cascades of secondary particles (air showers) in the atmosphere and their masses are inferred from measurements of the atmospheric depth of the shower maximum, Xmax, or the composition of shower particles reaching the ground. Current measurements suffer from either low precision, or a low duty cycle and a high energy threshold. Radio detection of cosmic rays is a rapidly developing technique, suitable for determination of Xmax with a duty cycle of in principle nearly 100%. The radiation is generated by the separation of relativistic charged particles in the geomagnetic field and a negative charge excess in the shower front. Here we report radio measurements of Xmax with a mean precision of 16 g/cm^2 between 10^{17}-10^{17.5} eV. Because of the high resolution in $Xmax we can determine the mass spectrum and find a mixed composition, containing a light mass fraction of ~80%. Unless the extragalactic component becomes significant already below 10^{17.5} eV, our measurements indicate an additional Galactic component dominating at this energy range.

The QueSERA Study Group (QSG) have been tasked by the RadioNet Board to produce a White Paper on the future organisation and coordination of radio astronomy in Europe. This White Paper describes the options discussed by the QSG, and our conclusions on how to move forward. We propose, that as a first step, RadioNet-work, be established as an entity that persists between EC contracts, and that takes responsibility for preparing or coordinating responses to EC opportunities specific to the field of radio astronomy research infrastructures. RadioNet-work should provide a safety net that ensures that cooperation and collaboration between the various radio astronomy partners in Europe is maintained with or without EC funding.

Air showers induced by cosmic rays create nanosecond pulses detectable at radio frequencies. These pulses have been measured successfully in the past few years at the LOw Frequency ARray (LOFAR) and are used to study the properties of cosmic rays. For a complete understanding of this phenomenon and the underlying physical processes, an absolute calibration of the detecting antenna system is needed. We present three approaches that were used to check and improve the antenna model of LOFAR and to provide an absolute calibration of the whole system for air shower measurements. Two methods are based on calibrated reference sources and one on a calibration approach using the diffuse radio emission of the Galaxy, optimized for short data-sets. An accuracy of 19% in amplitude is reached. The absolute calibration is also compared to predictions from air shower simulations. These results are used to set an absolute energy scale for air shower measurements and can be used as a basis for an absolute scale for the measurement of astronomical transients with LOFAR.

We present the results of a four-month campaign searching for low-frequency radio transients near the North Celestial Pole with the Low-Frequency Array (LOFAR), as part of the Multifrequency Snapshot Sky Survey (MSSS). The data were recorded between 2011 December and 2012 April and comprised 2149 11-minute snapshots, each covering 175 deg^2. We have found one convincing candidate astrophysical transient, with a duration of a few minutes and a flux density at 60 MHz of 15-25 Jy. The transient does not repeat and has no obvious optical or high-energy counterpart, as a result of which its nature is unclear. The detection of this event implies a transient rate at 60 MHz of 3.9 (+14.7, -3.7) x 10^-4 day^-1 deg^-2, and a transient surface density of 1.5 x 10^-5 deg^-2, at a 7.9-Jy limiting flux density and ~10-minute time-scale. The campaign data were also searched for transients at a range of other time-scales, from 0.5 to 297 min, which allowed us to place a range of limits on transient rates at 60 MHz as a function of observation duration.

Context. Observing Jupiter's synchrotron emission from the Earth remains today the sole method to scrutinize the distribution and dynamical behavior of the ultra energetic electrons magnetically trapped around the planet (because in-situ particle data are limited in the inner magnetosphere). Aims. We perform the first resolved and low-frequency imaging of the synchrotron emission with LOFAR at 127 MHz. The radiation comes from low energy electrons (~1-30 MeV) which map a broad region of Jupiter's inner magnetosphere. Methods (see article for complete abstract) Results. The first resolved images of Jupiter's radiation belts at 127-172 MHz are obtained along with total integrated flux densities. They are compared with previous observations at higher frequencies and show a larger extent of the synchrotron emission source (>=4 $R_J$). The asymmetry and the dynamic of east-west emission peaks are measured and the presence of a hot spot at lambda_III=230 {\deg} $\pm$ 25 {\deg}. Spectral flux density measurements are on the low side of previous (unresolved) ones, suggesting a low-frequency turnover and/or time variations of the emission spectrum. Conclusions. LOFAR is a powerful and flexible planetary imager. The observations at 127 MHz depict an extended emission up to ~4-5 planetary radii. The similarities with high frequency results reinforce the conclusion that: i) the magnetic field morphology primarily shapes the brightness distribution of the emission and ii) the radiating electrons are likely radially and latitudinally distributed inside about 2 $R_J$. Nonetheless, the larger extent of the brightness combined with the overall lower flux density, yields new information on Jupiter's electron distribution, that may shed light on the origin and mode of transport of these particles.

LOFAR offers the unique capability of observing pulsars across the 10-240 MHz frequency range with a fractional bandwidth of roughly 50%. This spectral range is well-suited for studying the frequency evolution of pulse profile morphology caused by both intrinsic and extrinsic effects: such as changing emission altitude in the pulsar magnetosphere or scatter broadening by the interstellar medium, respectively. The magnitude of most of these effects increases rapidly towards low frequencies. LOFAR can thus address a number of open questions about the nature of radio pulsar emission and its propagation through the interstellar medium. We present the average pulse profiles of 100 pulsars observed in the two LOFAR frequency bands: High Band (120-167 MHz, 100 profiles) and Low Band (15-62 MHz, 26 profiles). We compare them with Westerbork Synthesis Radio Telescope (WSRT) and Lovell Telescope observations at higher frequencies (350 and1400 MHz) in order to study the profile evolution. The profiles are aligned in absolute phase by folding with a new set of timing solutions from the Lovell Telescope, which we present along with precise dispersion measures obtained with LOFAR. We find that the profile evolution with decreasing radio frequency does not follow a specific trend but, depending on the geometry of the pulsar, new components can enter into, or be hidden from, view. Nonetheless, in general our observations confirm the widening of pulsar profiles at low frequencies, as expected from radius-to-frequency mapping or birefringence theories. We offer this catalog of low-frequency pulsar profiles in a user friendly way via the EPN Database of Pulsar Profiles (http://www.epta.eu.org/epndb/).

We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.

Context. The Sun is an active source of radio emission that is often associated with energetic phenomena ranging from nanoflares to coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), numerous millisecond duration radio bursts have been reported, such as radio spikes or solar S bursts (where S stands for short). To date, these have neither been studied extensively nor imaged because of the instrumental limitations of previous radio telescopes. Aims. Here, Low Frequency Array (LOFAR) observations were used to study the spectral and spatial characteristics of a multitude of S bursts, as well as their origin and possible emission mechanisms. Methods. We used 170 simultaneous tied-array beams for spectroscopy and imaging of S bursts. Since S bursts have short timescales and fine frequency structures, high cadence (~50 ms) tied-array images were used instead of standard interferometric imaging, that is currently limited to one image per second. Results. On 9 July 2013, over 3000 S bursts were observed over a time period of ~8 hours. S bursts were found to appear as groups of short-lived (<1 s) and narrow-bandwidth (~2.5 MHz) features, the majority drifting at ~3.5 MHz/s and a wide range of circular polarisation degrees (2-8 times more polarised than the accompanying Type III bursts). Extrapolation of the photospheric magnetic field using the potential field source surface (PFSS) model suggests that S bursts are associated with a trans-equatorial loop system that connects an active region in the southern hemisphere to a bipolar region of plage in the northern hemisphere. Conclusions. We have identified polarised, short-lived solar radio bursts that have never been imaged before. They are observed at a height and frequency range where plasma emission is the dominant emission mechanism, however they possess some of the characteristics of electron-cyclotron maser emission.

PSR B0823+26, a 0.53-s radio pulsar, displays a host of emission phenomena over timescales of seconds to (at least) hours, including nulling, subpulse drifting, and mode-changing. Studying pulsars like PSR B0823+26 provides further insight into the relationship between these various emission phenomena and what they might teach us about pulsar magnetospheres. Here we report on the LOFAR discovery that PSR B0823+26 has a weak and sporadically emitting 'quiet' (Q) emission mode that is over 100 times weaker (on average) and has a nulling fraction forty-times greater than that of the more regularly-emitting 'bright' (B) mode. Previously, the pulsar has been undetected in the Q-mode, and was assumed to be nulling continuously. PSR B0823+26 shows a further decrease in average flux just before the transition into the B-mode, and perhaps truly turns off completely at these times. Furthermore, simultaneous observations taken with the LOFAR, Westerbork, Lovell, and Effelsberg telescopes between 110 MHz and 2.7 GHz demonstrate that the transition between the Q-mode and B-mode occurs within one single rotation of the neutron star, and that it is concurrent across the range of frequencies observed.

Using observations obtained with the LOw Fequency ARray (LOFAR), the Westerbork Synthesis Radio Telescope (WSRT) and archival Very Large Array (VLA) data, we have traced the radio emission to large scales in the complex source 4C 35.06 located in the core of the galaxy cluster Abell 407. At higher spatial resolution (~4"), the source was known to have two inner radio lobes spanning 31 kpc and a diffuse, low-brightness extension running parallel to them, offset by about 11 kpc (in projection). At 62 MHz, we detect the radio emission of this structure extending out to 210 kpc. At 1.4 GHz and intermediate spatial resolution (~30"), the structure appears to have a helical morphology. We have derived the characteristics of the radio spectral index across the source. We show that the source morphology is most likely the result of at least two episodes of AGN activity separated by a dormant period of around 35 Myr. The AGN is hosted by one of the galaxies located in the cluster core of Abell 407. We propose that it is intermittently active as it moves in the dense environment in the cluster core. Using LOFAR, we can trace the relic plasma from that episode of activity out to greater distances from the core than ever before. Using the the WSRT, we detect HI in absorption against the center of the radio source. The absorption profile is relatively broad (FWHM of 288 km/s), similar to what is found in other clusters. Understanding the duty cycle of the radio emission as well as the triggering mechanism for starting (or restarting) the radio-loud activity can provide important constraints to quantify the impact of AGN feedback on galaxy evolution. The study of these mechanisms at low frequencies using morphological and spectral information promises to bring new important insights in this field.

We present measurements of radio emission from cosmic ray air showers that took place during thunderstorms. The intensity and polarization patterns of these air showers are radically different from those measured during fair-weather conditions. With the use of a simple two-layer model for the atmospheric electric field, these patterns can be well reproduced by state-of-the-art simulation codes. This in turn provides a novel way to study atmospheric electric fields.

Context. The LOw Frequency ARray (LOFAR) radio telescope is a giant digital phased array interferometer with multiple antennas distributed in Europe. It provides discrete sets of Fourier components of the sky brightness. Recovering the original brightness distribution with aperture synthesis forms an inverse problem that can be solved by various deconvolution and minimization methods Aims. Recent papers have established a clear link between the discrete nature of radio interferometry measurement and the "compressed sensing" (CS) theory, which supports sparse reconstruction methods to form an image from the measured visibilities. Empowered by proximal theory, CS offers a sound framework for efficient global minimization and sparse data representation using fast algorithms. Combined with instrumental direction-dependent effects (DDE) in the scope of a real instrument, we developed and validated a new method based on this framework Methods. We implemented a sparse reconstruction method in the standard LOFAR imaging tool and compared the photometric and resolution performance of this new imager with that of CLEAN-based methods (CLEAN and MS-CLEAN) with simulated and real LOFAR data Results. We show that i) sparse reconstruction performs as well as CLEAN in recovering the flux of point sources; ii) performs much better on extended objects (the root mean square error is reduced by a factor of up to 10); and iii) provides a solution with an effective angular resolution 2-3 times better than the CLEAN images. Conclusions. Sparse recovery gives a correct photometry on high dynamic and wide-field images and improved realistic structures of extended sources (of simulated and real LOFAR datasets). This sparse reconstruction method is compatible with modern interferometric imagers that handle DDE corrections (A- and W-projections) required for current and future instruments such as LOFAR and SKA

Adding VLBI capability to the SKA arrays will greatly broaden the science of the SKA, and is feasible within the current specifications. SKA-VLBI can be initially implemented by providing phased-array outputs for SKA1-MID and SKA1-SUR and using these extremely sensitive stations with other radio telescopes, and in SKA2 by realising a distributed configuration providing baselines up to thousands of km, merging it with existing VLBI networks. The motivation for and the possible realization of SKA-VLBI is described in this paper.

Measuring radio emission from air showers offers a novel way to determine properties of the primary cosmic rays such as their mass and energy. Theory predicts that relativistic time compression effects lead to a ring of amplified emission which starts to dominate the emission pattern for frequencies above ~100 MHz. In this article we present the first detailed measurements of this structure. Ring structures in the radio emission of air showers are measured with the LOFAR radio telescope in the frequency range of 110 - 190 MHz. These data are well described by CoREAS simulations. They clearly confirm the importance of including the index of refraction of air as a function of height. Furthermore, the presence of the Cherenkov ring offers the possibility for a geometrical measurement of the depth of shower maximum, which in turn depends on the mass of the primary particle.

Aims. An efficient means of locating calibrator sources for International LOFAR is developed and used to determine the average density of usable calibrator sources on the sky for subarcsecond observations at 140 MHz. Methods. We used the multi-beaming capability of LOFAR to conduct a fast and computationally inexpensive survey with the full International LOFAR array. Sources were pre-selected on the basis of 325 MHz arcminute-scale flux density using existing catalogues. By observing 30 different sources in each of the 12 sets of pointings per hour, we were able to inspect 630 sources in two hours to determine if they possess a sufficiently bright compact component to be usable as LOFAR delay calibrators. Results. Over 40% of the observed sources are detected on multiple baselines between international stations and 86 are classified as satisfactory calibrators. We show that a flat low-frequency spectrum (from 74 to 325 MHz) is the best predictor of compactness at 140 MHz. We extrapolate from our sample to show that the density of calibrators on the sky that are sufficiently bright to calibrate dispersive and non-dispersive delays for the International LOFAR using existing methods is 1.0 per square degree. Conclusions. The observed density of satisfactory delay calibrator sources means that observations with International LOFAR should be possible at virtually any point in the sky, provided that a fast and efficient search using the methodology described here is conducted prior to the observation to identify the best calibrator.

We present LOFAR Low Band observations of the Bootes and 3C295 fields. Our images made at 34, 46, and 62 MHz reach noise levels of 12, 8, and 5 mJy beam$^{-1}$, making them the deepest images ever obtained in this frequency range. In total, we detect between 300 and 400 sources in each of these images, covering an area of 17 to 52 deg$^{2}$. From the observations we derive Euclidean-normalized differential source counts. The 62 MHz source counts agree with previous GMRT 153 MHz and VLA 74 MHz differential source counts, scaling with a spectral index of $-0.7$. We find that a spectral index scaling of $-0.5$ is required to match up the LOFAR 34 MHz source counts. This result is also in agreement with source counts from the 38 MHz 8C survey, indicating that the average spectral index of radio sources flattens towards lower frequencies. We also find evidence for spectral flattening using the individual flux measurements of sources between 34 and 1400 MHz and by calculating the spectral index averaged over the source population. To select ultra-steep spectrum ($\alpha < -1.1$) radio sources, that could be associated with massive high redshift radio galaxies, we compute spectral indices between 62 MHz, 153 MHz and 1.4 GHz for sources in the Bo\"otes field. We cross-correlate these radio sources with optical and infrared catalogues and fit the spectral energy distribution to obtain photometric redshifts. We find that most of these ultra-steep spectrum sources are located in the $ 0.7 \lesssim z \lesssim 2.5$ range.

We have conducted two pilot surveys for radio pulsars and fast transients with the Low-Frequency Array (LOFAR) around 140 MHz and here report on the first low-frequency fast-radio burst limit and the discovery of two new pulsars. The first survey, the LOFAR Pilot Pulsar Survey (LPPS), observed a large fraction of the northern sky, ~1.4 x 10^4 sq. deg, with 1-hr dwell times. Each observation covered ~75 sq. deg using 7 independent fields formed by incoherently summing the high-band antenna fields. The second pilot survey, the LOFAR Tied-Array Survey (LOTAS), spanned ~600 sq. deg, with roughly a 5-fold increase in sensitivity compared with LPPS. Using a coherent sum of the 6 LOFAR "Superterp" stations, we formed 19 tied-array beams, together covering 4 sq. deg per pointing. From LPPS we derive a limit on the occurrence, at 142 MHz, of dispersed radio bursts of < 150 /day/sky, for bursts brighter than S > 107 Jy for the narrowest searched burst duration of 0.66 ms. In LPPS, we re-detected 65 previously known pulsars. LOTAS discovered two pulsars, the first with LOFAR or any digital aperture array. LOTAS also re-detected 27 previously known pulsars. These pilot studies show that LOFAR can efficiently carry out all-sky surveys for pulsars and fast transients, and they set the stage for further surveying efforts using LOFAR and the planned low-frequency component of the Square Kilometer Array.

The Sun is an active source of radio emission which is often associated with energetic phenomena such as solar flares and coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), the Sun has not been imaged extensively because of the instrumental limitations of previous radio telescopes. Here, the combined high spatial, spectral and temporal resolution of the Low Frequency Array (LOFAR) was used to study solar Type III radio bursts at 30-90 MHz and their association with CMEs. The Sun was imaged with 126 simultaneous tied-array beams within 5 solar radii of the solar centre. This method offers benefits over standard interferometric imaging since each beam produces high temporal (83 ms) and spectral resolution (12.5 kHz) dynamic spectra at an array of spatial locations centred on the Sun. LOFAR's standard interferometric output is currently limited to one image per second. Over a period of 30 minutes, multiple Type III radio bursts were observed, a number of which were found to be located at high altitudes (4 solar radii from the solar center at 30 MHz) and to have non-radial trajectories. These bursts occurred at altitudes in excess of values predicted by 1D radial electron density models. The non-radial high altitude Type III bursts were found to be associated with the expanding flank of a CME. The CME may have compressed neighbouring streamer plasma producing larger electron densities at high altitudes, while the non-radial burst trajectories can be explained by the deflection of radial magnetic fields as the CME expanded in the low corona.

This study aims to characterise the polarized foreground emission in the ELAIS-N1 field and to address its possible implications for the extraction of the cosmological 21-cm signal from the Low-Frequency Array - Epoch of Reionization (LOFAR-EoR) data. We use the high band antennas of LOFAR to image this region and RM-synthesis to unravel structures of polarized emission at high Galactic latitudes. The brightness temperature of the detected Galactic emission is on average 4 K in polarized intensity and covers the range from -10 to +13rad m^-2 in Faraday depth. The total polarized intensity and polarization angle show a wide range of morphological features. We have also used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. The LOFAR and WSRT images show a similar complex morphology, at comparable brightness levels, but their spatial correlation is very low. The fractional polarization at 150 MHz, expressed as a percentage of the total intensity, amounts to 1.5%. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies. The wide frequency range, good angular resolution and good sensitivity make LOFAR an exquisite instrument for studying Galactic polarized emission at a resolution of 1-2 rad m^-2 in Faraday depth. The different polarised patterns observed at 150 MHz and 350 MHz are consistent with different source distributions along the line of sight wring in a variety of Faraday thin regions of emission. The presence of polarised foregrounds is a serious complication for Epoch of Reionization experiments. To avoid the leakage of polarized emission into total intensity, which can depend on frequency, we need to calibrate the instrumental polarization across the field of view to a small fraction of 1%.

Extensive air showers, induced by high energy cosmic rays impinging on the Earth's atmosphere, produce radio emission that is measured with the LOFAR radio telescope. As the emission comes from a finite distance of a few kilometers, the incident wavefront is non-planar. A spherical, conical or hyperbolic shape of the wavefront has been proposed, but measurements of individual air showers have been inconclusive so far. For a selected high-quality sample of 161 measured extensive air showers, we have reconstructed the wavefront by measuring pulse arrival times to sub-nanosecond precision in 200 to 350 individual antennas. For each measured air shower, we have fitted a conical, spherical, and hyperboloid shape to the arrival times. The fit quality and a likelihood analysis show that a hyperboloid is the best parametrization. Using a non-planar wavefront shape gives an improved angular resolution, when reconstructing the shower arrival direction. Furthermore, a dependence of the wavefront shape on the shower geometry can be seen. This suggests that it will be possible to use a wavefront shape analysis to get an additional handle on the atmospheric depth of the shower maximum, which is sensitive to the mass of the primary particle.

We present the first detection of carbon radio recombination line absorption along the line of sight to Cygnus A. The observations were carried out with the LOw Frequency ARray in the 33 to 57 MHz range. These low frequency radio observations provide us with a new line of sight to study the diffuse, neutral gas in our Galaxy. To our knowledge this is the first time that foreground Milky Way recombination line absorption has been observed against a bright extragalactic background source. By stacking 48 carbon $\alpha$ lines in the observed frequency range we detect carbon absorption with a signal-to-noise ratio of about 5. The average carbon absorption has a peak optical depth of 2$\times$10$^{-4}$, a line width of 10 km s$^{-1}$ and a velocity of +4 km s$^{-1}$ with respect to the local standard of rest. The associated gas is found to have an electron temperature $T_{e}\sim$ 110 K and density $n_{e}\sim$ 0.06 cm$^{-3}$. These properties imply that the observed carbon $\alpha$ absorption likely arises in the cold neutral medium of the Orion arm of the Milky Way. Hydrogen and helium lines were not detected to a 3$\sigma$ peak optical depth limit of 1.5$\times$10$^{-4}$ for a 4 km s$^{-1}$ channel width. Radio recombination lines associated with Cygnus A itself were also searched for, but are not detected. We set a 3$\sigma$ upper limit of 1.5$\times$10$^{-4}$ for the peak optical depth of these lines for a 4 km s$^{-1}$ channel width.

The low frequency array (LOFAR), is the first radio telescope designed with the capability to measure radio emission from cosmic-ray induced air showers in parallel with interferometric observations. In the first $\sim 2\,\mathrm{years}$ of observing, 405 cosmic-ray events in the energy range of $10^{16} - 10^{18}\,\mathrm{eV}$ have been detected in the band from $30 - 80\,\mathrm{MHz}$. Each of these air showers is registered with up to $\sim1000$ independent antennas resulting in measurements of the radio emission with unprecedented detail. This article describes the dataset, as well as the analysis pipeline, and serves as a reference for future papers based on these data. All steps necessary to achieve a full reconstruction of the electric field at every antenna position are explained, including removal of radio frequency interference, correcting for the antenna response and identification of the pulsed signal.

The characteristic outer scale of turbulence and the ratio of the random to ordered components of the magnetic field are key parameters to characterise magnetic turbulence in the interstellar gas, which affects the propagation of cosmic rays within the Galaxy. We provide new constraints to those two parameters. We use the LOw Frequency ARray (LOFAR) to image the diffuse continuum emission in the Fan region at (l,b) (137.0,+7.0) at 80"x70" resolution in the range [146,174] MHz. We detect multi-scale fluctuations in the Galactic synchrotron emission and compute their power spectrum. Applying theoretical estimates and derivations from the literature for the first time, we derive the outer scale of turbulence and the ratio of random to ordered magnetic field from the characteristics of these fluctuations . We obtain the deepest image of the Fan region to date and find diffuse continuum emission within the primary beam. The power spectrum of the foreground synchrotron fluctuations displays a power law behaviour for scales between 100 and 8 arcmin with a slope of (-1.84+/-0.19). We find an upper limit of about 20 pc for the outer scale of the magnetic interstellar turbulence toward the Fan region. We also find a variation of the ratio of random to ordered field as a function of Galactic coordinates, supporting different turbulent regimes. We use power spectra fluctuations from LOFAR as well as earlier GMRT and WSRT observations to constrain the outer scale of turbulence of the Galactic synchrotron foreground, finding a range of plausible values of 10-20 pc. Then, we use this information to deduce lower limits of the ratio of ordered to random magnetic field strength. These are found to be 0.3, 0.3, and 0.5 for the LOFAR, WSRT and GMRT fields considered respectively. Both these constraints are in agreement with previous estimates.

Faint undetected sources of radio-frequency interference (RFI) might become visible in long radio observations when they are consistently present over time. Thereby, they might obstruct the detection of the weak astronomical signals of interest. This issue is especially important for Epoch of Reionisation (EoR) projects that try to detect the faint redshifted HI signals from the time of the earliest structures in the Universe. We explore the RFI situation at 30-163 MHz by studying brightness histograms of visibility data observed with LOFAR, similar to radio-source-count analyses that are used in cosmology. An empirical RFI distribution model is derived that allows the simulation of RFI in radio observations. The brightness histograms show an RFI distribution that follows a power-law distribution with an estimated exponent around -1.5. With several assumptions, this can be explained with a uniform distribution of terrestrial radio sources whose radiation follows existing propagation models. Extrapolation of the power law implies that the current LOFAR EoR observations should be severely RFI limited if the strength of RFI sources remains strong after time integration. This is in contrast with actual observations, which almost reach the thermal noise and are thought not to be limited by RFI. Therefore, we conclude that it is unlikely that there are undetected RFI sources that will become visible in long observations. Consequently, there is no indication that RFI will prevent an EoR detection with LOFAR.

LOFAR, the LOw-Frequency ARray, is a new-generation radio interferometer constructed in the north of the Netherlands and across europe. Utilizing a novel phased-array design, LOFAR covers the largely unexplored low-frequency range from 10-240 MHz and provides a number of unique observing capabilities. Spreading out from a core located near the village of Exloo in the northeast of the Netherlands, a total of 40 LOFAR stations are nearing completion. A further five stations have been deployed throughout Germany, and one station has been built in each of France, Sweden, and the UK. Digital beam-forming techniques make the LOFAR system agile and allow for rapid repointing of the telescope as well as the potential for multiple simultaneous observations. With its dense core array and long interferometric baselines, LOFAR achieves unparalleled sensitivity and angular resolution in the low-frequency radio regime. The LOFAR facilities are jointly operated by the International LOFAR Telescope (ILT) foundation, as an observatory open to the global astronomical community. LOFAR is one of the first radio observatories to feature automated processing pipelines to deliver fully calibrated science products to its user community. LOFAR's new capabilities, techniques and modus operandi make it an important pathfinder for the Square Kilometre Array (SKA). We give an overview of the LOFAR instrument, its major hardware and software components, and the core science objectives that have driven its design. In addition, we present a selection of new results from the commissioning phase of this new radio observatory.

Faraday rotation measurements using the current and next generation of low-frequency radio telescopes will provide a powerful probe of astronomical magnetic fields. However, achieving the full potential of these measurements requires accurate removal of the time-variable ionospheric Faraday rotation contribution. We present ionFR, a code that calculates the amount of ionospheric Faraday rotation for a specific epoch, geographic location, and line-of-sight. ionFR uses a number of publicly available, GPS-derived total electron content maps and the most recent release of the International Geomagnetic Reference Field. We describe applications of this code for the calibration of radio polarimetric observations, and demonstrate the high accuracy of its modeled ionospheric Faraday rotations using LOFAR pulsar observations. These show that we can accurately determine some of the highest-precision pulsar rotation measures ever achieved. Precision rotation measures can be used to monitor rotation measure variations - either intrinsic or due to the changing line-of-sight through the interstellar medium. This calibration is particularly important for nearby sources, where the ionosphere can contribute a significant fraction of the observed rotation measure. We also discuss planned improvements to ionFR, as well as the importance of ionospheric Faraday rotation calibration for the emerging generation of low-frequency radio telescopes, such as the SKA and its pathfinders.

Cassiopeia A was observed using the Low-Band Antennas of the LOw Frequency ARray (LOFAR) with high spectral resolution. This allowed a search for radio recombination lines (RRLs) along the line-of-sight to this source. Five carbon-alpha RRLs were detected in absorption between 40 and 50 MHz with a signal-to-noise ratio of > 5 from two independent LOFAR datasets. The derived line velocities (v_LSR ~ -50 km/s) and integrated optical depths (~ 13 s^-1) of the RRLs in our spectra, extracted over the whole supernova remnant, are consistent within each LOFAR dataset and with those previously reported. For the first time, we are able to extract spectra against the brightest hotspot of the remnant at frequencies below 330 MHz. These spectra show significantly higher (15-80 %) integrated optical depths, indicating that there is small-scale angular structure on the order of ~1 pc in the absorbing gas distribution over the face of the remnant. We also place an upper limit of 3 x 10^-4 on the peak optical depths of hydrogen and helium RRLs. These results demonstrate that LOFAR has the desired spectral stability and sensitivity to study faint recombination lines in the decameter band.

Some radio pulsars show clear drifting subpulses, in which subpulses are seen to drift in pulse longitude in a systematic pattern. Here we examine how the drifting subpulses of PSR B0809+74 evolve with time and observing frequency. We show that the subpulse period (P3) is constant on timescales of days, months and years, and between 14-5100 MHz. Despite this, the shapes of the driftbands change radically with frequency. Previous studies have concluded that, while the subpulses appear to move through the pulse window approximately linearly at low frequencies (< 500 MHz), a discrete step of 180 degrees in subpulse phase is observed at higher frequencies (> 820 MHz) near to the peak of the average pulse profile. We use LOFAR, GMRT, GBT, WSRT and Effelsberg 100-m data to explore the frequency-dependence of this phase step. We show that the size of the subpulse phase step increases gradually, and is observable even at low frequencies. We attribute the subpulse phase step to the presence of two separate driftbands, whose relative arrival times vary with frequency - one driftband arriving 30 pulses earlier at 20 MHz than it does at 1380 MHz, whilst the other arrives simultaneously at all frequencies. The drifting pattern which is observed here cannot be explained by either the rotating carousel model or the surface oscillation model, and could provide new insight into the physical processes happening within the pulsar magnetosphere.

The aim of the LOFAR Epoch of Reionization (EoR) project is to detect the spectral fluctuations of the redshifted HI 21cm signal. This signal is weaker by several orders of magnitude than the astrophysical foreground signals and hence, in order to achieve this, very long integrations, accurate calibration for stations and ionosphere and reliable foreground removal are essential. One of the prospective observing windows for the LOFAR EoR project will be centered at the North Celestial Pole (NCP). We present results from observations of the NCP window using the LOFAR highband antenna (HBA) array in the frequency range 115 MHz to 163 MHz. The data were obtained in April 2011 during the commissioning phase of LOFAR. We used baselines up to about 30 km. With about 3 nights, of 6 hours each, effective integration we have achieved a noise level of about 100 microJy/PSF in the NCP window. Close to the NCP, the noise level increases to about 180 microJy/PSF, mainly due to additional contamination from unsubtracted nearby sources. We estimate that in our best night, we have reached a noise level only a factor of 1.4 above the thermal limit set by the noise from our Galaxy and the receivers. Our continuum images are several times deeper than have been achieved previously using the WSRT and GMRT arrays. We derive an analytical explanation for the excess noise that we believe to be mainly due to sources at large angular separation from the NCP.

M87 is a giant elliptical galaxy located in the centre of the Virgo cluster, which harbours a supermassive black hole of mass 6.4x10^9 M_sun, whose activity is responsible for the extended (80 kpc) radio lobes that surround the galaxy. The energy generated by matter falling onto the central black hole is ejected and transferred to the intra-cluster medium via a relativistic jet and morphologically complex systems of buoyant bubbles, which rise towards the edges of the extended halo. Here we present the first observations made with the new Low-Frequency Array (LOFAR) of M87 at frequencies down to 20 MHz. Images of M87 were produced at low radio frequencies never explored before at these high spatial resolution and dynamic range. To disentangle different synchrotron models and place constraints on source magnetic field, age and energetics, we also performed a detailed spectral analysis of M87 extended radio-halo using these observations together with archival data. We do not find any sign of new extended emissions; on the contrary the source appears well confined by the high pressure of the intra-cluster medium. A continuous injection of relativistic electrons is the model that best fits our data, and provides a scenario in which the lobes are still supplied by fresh relativistic particles from the active galactic nuclei. We suggest that the discrepancy between the low-frequency radio-spectral slope in the core and in the halo implies a strong adiabatic expansion of the plasma as soon as it leaves the core area. The extended halo has an equipartition magnetic field strength of ~10 uG, which increases to ~13 uG in the zones where the particle flows are more active. The continuous injection model for synchrotron ageing provides an age for the halo of ~40 Myr, which in turn provides a jet kinetic power of 6-10x10^44 erg/s.

X-ray binaries have been an important key in understanding the jet-disc symbiosis in accreting black holes on all mass scales, from stellar-mass to supermassive black holes. SS433 was the first Galactic XRB that has been extensively studied in the radio regime. The radio properties, including the highest angular resolution data can now be better understood in the framework for accretion disc state transitions that is observed in microquasars (black hole X-ray binary systems). SS433 remains unique in various ways to date, and there is still much to learn about black hole accretion phenomena. In the meantime, the electronic very long baseline (e-VLBI) developments at the European VLBI Network (EVN) has allowed us to study microquasars and other transients at milliarcsecond resolutions more flexibly than was possible before. Even more new opportunities will arise as the SKA pathfinders become operational.

Aims: This paper discusses the spectral occupancy for performing radio astronomy with the Low-Frequency Array (LOFAR), with a focus on imaging observations. Methods: We have analysed the radio-frequency interference (RFI) situation in two 24-h surveys with Dutch LOFAR stations, covering 30-78 MHz with low-band antennas and 115-163 MHz with high-band antennas. This is a subset of the full frequency range of LOFAR. The surveys have been observed with a 0.76 kHz / 1 s resolution. Results: We measured the RFI occupancy in the low and high frequency sets to be 1.8% and 3.2% respectively. These values are found to be representative values for the LOFAR radio environment. Between day and night, there is no significant difference in the radio environment. We find that lowering the current observational time and frequency resolutions of LOFAR results in a slight loss of flagging accuracy. At LOFAR's nominal resolution of 0.76 kHz and 1 s, the false-positives rate is about 0.5%. This rate increases approximately linearly when decreasing the data frequency resolution. Conclusions: Currently, by using an automated RFI detection strategy, the LOFAR radio environment poses no perceivable problems for sensitive observing. It remains to be seen if this is still true for very deep observations that integrate over tens of nights, but the situation looks promising. Reasons for the low impact of RFI are the high spectral and time resolution of LOFAR; accurate detection methods; strong filters and high receiver linearity; and the proximity of the antennas to the ground. We discuss some strategies that can be used once low-level RFI starts to become apparent. It is important that the frequency range of LOFAR remains free of broadband interference, such as DAB stations and windmills.

Dispersion in the interstellar medium is a well known phenomenon that follows a simple relationship, which has been used to predict the time delay of dispersed radio pulses since the late 1960s. We performed wide-band simultaneous observations of four pulsars with LOFAR (at 40-190 MHz), the 76-m Lovell Telescope (at 1400 MHz) and the Effelsberg 100-m Telescope (at 8000 MHz) to test the accuracy of the dispersion law over a broad frequency range. In this paper we present the results of these observations which show that the dispersion law is accurate to better than 1 part in 100000 across our observing band. We use this fact to constrain some of the properties of the ISM along the line-of-sight and use the lack of any aberration or retardation effects to determine upper limits on emission heights in the pulsar magnetosphere. We also discuss the effect of pulse profile evolution on our observations, and the implications that it could have for precision pulsar timing projects such as the detection of gravitational waves with pulsar timing arrays.

Abell 2256 is one of the best known examples of a galaxy cluster hosting large-scale diffuse radio emission that is unrelated to individual galaxies. It contains both a giant radio halo and a relic, as well as a number of head-tail sources and smaller diffuse steep-spectrum radio sources. The origin of radio halos and relics is still being debated, but over the last years it has become clear that the presence of these radio sources is closely related to galaxy cluster merger events. Here we present the results from the first LOFAR Low band antenna (LBA) observations of Abell 2256 between 18 and 67 MHz. To our knowledge, the image presented in this paper at 63 MHz is the deepest ever obtained at frequencies below 100 MHz in general. Both the radio halo and the giant relic are detected in the image at 63 MHz, and the diffuse radio emission remains visible at frequencies as low as 20 MHz. The observations confirm the presence of a previously claimed ultra-steep spectrum source to the west of the cluster center with a spectral index of -2.3 \pm 0.4 between 63 and 153 MHz. The steep spectrum suggests that this source is an old part of a head-tail radio source in the cluster. For the radio relic we find an integrated spectral index of -0.81 \pm 0.03, after removing the flux contribution from the other sources. This is relatively flat which could indicate that the efficiency of particle acceleration at the shock substantially changed in the last \sim 0.1 Gyr due to an increase of the shock Mach number. In an alternative scenario, particles are re-accelerated by some mechanism in the downstream region of the shock, resulting in the relatively flat integrated radio spectrum. In the radio halo region we find indications of low-frequency spectral steepening which may suggest that relativistic particles are accelerated in a rather inhomogeneous turbulent region.

When an ultra-high energy neutrino or cosmic ray strikes the Lunar surface a radio-frequency pulse is emitted. We plan to use the LOFAR radio telescope to detect these pulses. In this work we propose an efficient trigger implementation for LOFAR optimized for the observation of short radio pulses.

Low frequency radio waves, while challenging to observe, are a rich source of information about pulsars. The LOw Frequency ARray (LOFAR) is a new radio interferometer operating in the lowest 4 octaves of the ionospheric "radio window": 10-240MHz, that will greatly facilitate observing pulsars at low radio frequencies. Through the huge collecting area, long baselines, and flexible digital hardware, it is expected that LOFAR will revolutionize radio astronomy at the lowest frequencies visible from Earth. LOFAR is a next-generation radio telescope and a pathfinder to the Square Kilometre Array (SKA), in that it incorporates advanced multi-beaming techniques between thousands of individual elements. We discuss the motivation for low-frequency pulsar observations in general and the potential of LOFAR in addressing these science goals. We present LOFAR as it is designed to perform high-time-resolution observations of pulsars and other fast transients, and outline the various relevant observing modes and data reduction pipelines that are already or will soon be implemented to facilitate these observations. A number of results obtained from commissioning observations are presented to demonstrate the exciting potential of the telescope. This paper outlines the case for low frequency pulsar observations and is also intended to serve as a reference for upcoming pulsar/fast transient science papers with LOFAR.

B2352+495 is a prototypical example of a Compact Symmetric Object (CSO). It has a double radio lobe symmetrically located with respect to a central flat spectrum radio core (the location of the AGN) and has a physical extent of less than 200 pc. In this work we report VLBA observation of 21 cm HI absorption toward B2352+495 to investigate the properties of this remarkable radio source, in particular, to explore whether the radio emission can be confined by circumnuclear material (frustration scenario) or whether the source is likely to be young. We confirmed the two HI absorption features previously detected toward B2352+495 - a broad line nearly centered at the systemic velocity of the galaxy and a narrow redshifted component. The atomic gas from the broad absorption component is likely associated with circumnuclear material, consistent with the current paradigm of clumpy HI distribution in toroidal structures around supermassive black holes.

We report on a multi-frequency, multi-epoch campaign of Very Long Baseline Interferometry observations of the radio galaxy 1946+708 using the VLBA and a Global VLBI array. From these high-resolution observations we deduce the kinematic age of the radio source to be $\sim$4000 years, comparable with the ages of other Compact Symmetric Objects (CSOs). Ejections of pairs of jet components appears to take place on time scales of 10 years and these components in the jet travel outward at intrinsic velocities between 0.6 and 0.9 c. From the constraint that jet components cannot have intrinsic velocities faster than light, we derive H_0 > 57 km s^-1 Mpc^-1 from the fastest pair of components launched from the core. We provide strong evidence for the ejection of a new pair of components in ~1997. From the trajectories of the jet components we deduce that the jet is most likely to be helically confined, rather than purely ballistic in nature.

We present images from a long term program (MOJAVE: Monitoring of Jets in AGN with VLBA Experiments) to survey the structure and evolution of parsec-scale jet phenomena associated with bright radio-loud active galaxies in the northern sky. The observations consist of 2424 15 GHz VLBA images of a complete flux-density limited sample of 135 AGN above declination -20 degrees, spanning the period 1994 August to 2007 September. These data were acquired as part of the MOJAVE and 2 cm Survey programs, and from the VLBA archive. The sample selection criteria are based on multi-epoch parsec-scale (VLBA) flux density, and heavily favor highly variable and compact blazars. The sample includes nearly all the most prominent blazars in the northern sky, and is well-suited for statistical analysis and comparison with studies at other wavelengths. Our multi-epoch and stacked-epoch images show 94% of the sample to have apparent one-sided jet morphologies, most likely due to the effects of relativistic beaming. Of the remaining sources, five have two-sided parsec-scale jets, and three are effectively unresolved by the VLBA at 15 GHz, with essentially all of the flux density contained within a few tenths of a milliarcsecond.

Context: This is the second of a series of papers presenting VLBI observations of the 293 Caltech-Jodrell Bank Flat-Spectrum (hereafter CJF) sources and their analysis. Aims: To obtain a consistent motion dataset large enough to allow the systematic properties of the population to be studied. Methods: We present the detailed kinematic analysis of the complete flux-density limited CJF survey. We computed 2-D kinematic models based on the optimal model-fitting parameters of multi-epoch VLBA observations. This allows us to calculate not only radial, but also orthogonal motions, and thus to study curvature and acceleration. Statistical tests of the motions measured and their reliability have been performed. A correlation analysis between the derived apparent motions, luminosities, spectral indices, and core dominance and the resulting consequences is described. Results: With at least one velocity in each of 237 sources, this sample is much larger than any available before and allows a meaningful statistical investigation of apparent motions and any possible correlations with other parameters in AGN jets. Conclusions: This AGN survey provides the basis for any statistical analysis of jet and jet-component properties.

Context: This is the first of a series of papers presenting VLBI observations of the 293 Caltech-Jodrell Bank Flat-Spectrum (hereafter CJF) sources and their analysis. Aims: One of the major goals of the CJF is to make a statistical study of the apparent velocities of the sources. Methods: We have conducted global VLBI and VLBA observations at 5 GHz since 1990, accumulating thirteen separate observing campaigns. Results: We present here an overview of the observations, give details of the data reduction and present the source parameters resulting from a model-fitting procedure. For every source at every observing epoch, an image is shown, built up by restoring the model-fitted components, convolved with the clean beam, into the residual image, which was made by Fourier transforming the visibility data after first subtracting the model-fitted components in the uv-plane. Overplotted we show symbols to represent the model components. Conclusions: We have produced VLBI images of all but 5 of the 293 sources in the complete CJF sample at several epochs and investigated the kinematics of 266 AGN.

We discuss results from a decade long program to study the fine-scale structure and the kinematics of relativistic AGN jets with the aim of better understanding the acceleration and collimation of the relativistic plasma forming AGN jets. From the observed distribution of brightness temperature, apparent velocity, flux density, time variability, and apparent luminosity, the intrinsic properties of the jets including Lorentz factor, luminosity, orientation, and brightness temperature are discussed. Special attention is given to the jet in M87, which has been studied over a wide range of wavelengths and which, due to its proximity, is observed with excellent spatial resolution. Most radio jets appear quite linear, but we also observe curved non-linear jets and non-radial motions. Sometimes, different features in a given jet appear to follow the same curved path but there is evidence for ballistic trajectories as well. The data are best fit with a distribution of Lorentz factors extending up to gamma ~30 and intrinsic luminosity up to ~10^26 W/Hz. In general, gamma-ray quasars may have somewhat larger Lorentz factors than non gamma-ray quasars. Initially the observed brightness temperature near the base of the jet extend up to ~5x10^13 K which is well in excess of the inverse Compton limit and corresponds to a large excess of particle energy over magnetic energy. However, more typically, the observed brightness temperatures are ~2x10^11 K, i.e., closer to equipartition.

We have studied the interrelation of young AGN with their hosts. The objects of study are the young and powerful GPS and CSS radio sources. Due to their small size, GPS and CSS sources are excellent probes of this relation. Furhthermore, their young age allows us to compare them to the larger, old radio sources and establish a time-line evolution of this relation. Combining imaging and spectroscopy at UV, optical and radio wavelengths we find evidence of strong interaction between the host and the radio source. The presence and expansion of the radio source clearly affects the properties and evolution of the host. Furthermore, the radio source and host significantly affect each other's evolution. We describe our results and how these interactions take place.

Relations between the observed quantities for a beamed radio jet, apparent transverse speed and apparent luminosity (beta_app,L), and the intrinsic quantities, Lorentz factor and intrinsic luminosity (gamma,L_o), are investigated. The inversion from measured to intrinsic values is not unique, but approximate limits to gamma and L_o can be found using probability arguments. Roughly half the sources in a flux density--limited, beamed sample have a value of gamma close to the measured beta_app. The methods are applied to observations of 119 AGN jets made with the VLBA at 15 GHz during 1994-2002. The results strongly support the common relativistic beam model for an extragalactic radio jet. The (beta_app,L) data are closely bounded by a theoretical envelope, an aspect curve for gamma=32, L_o= 10^25 W/Hz. This gives limits to the maximum values of gamma and L_o in the sample: gamma_max about 32, and L_o,max ~ 10^26 W/Hz. No sources with both high beta_app and low L are observed. This is not the result of selection effects due to the observing limits, which are flux density S>0.5 Jy, and angular velocity mu<4 mas/yr. Many of the fastest quasars have a pattern Lorentz factor gamma_p close to that of the beam, gamma_b, but some of the slow quasars must have gamma_p<<gamma_b. Three of the 10 galaxies in the sample have a superluminal feature, with speeds up to beta_app about 6. The others are at most mildly relativistic. The galaxies are not off-axis versions of the powerful quasars, but Cygnus A might be an exception.

We present a new method for studying the intrinsic brightness temperatures of the parsec-scale jet cores of Active Galactic Nuclei (AGN). Our method uses observed superluminal motions and observed brightness temperatures for a large sample of AGN to constrain the characteristic intrinsic brightness temperature of the sample as a whole. To study changes in intrinsic brightness temperature, we assume that the Doppler factors of individual jets are constant in time as justified by their relatively small changes in observed flux density. We find that in their median-low brightness temperature state, the sources in our sample have a narrow range of intrinsic brightness temperatures centered on a characteristic temperature, T_int = 3 x 10^10 K, which is close to the value expected for equipartition, when the energy in the radiating particles equals the energy stored in the magnetic fields. However, in their maximum brightness state, we find that sources in our sample have a characteristic intrinsic brightness temperature greater than 2 x 10^11 K, which is well in excess of the equipartition temperature. In this state, we estimate the energy in radiating particles exceeds the energy in the magnetic field by a factor of ~ 10^5. We suggest that the excess of particle energy when sources are in their maximum brightness state is due to injection or acceleration of particles at the base of the jet. Our results suggest that the common method of estimating jet Doppler factors by using a single measurement of observed brightness temperature and/or the assumption of equipartition may lead to large scatter or systematic errors in the derived values.

We present and discuss European VLBI Network UHF band spectral line observations, made to localise the redshifted 21cm HI absorption known to occur in the subgalactic sized compact steep spectrum galaxies 3C 49 and 3C 268.3. We have detected HI absorption towards the western radio lobe of 3C 49 and the northern lobe of 3C 268.3. However, we cannot rule out the presence of similar amounts of HI towards the opposite and much fainter lobes. The radio lobes with detected HI absorption (1) are brighter and closer to the core than the opposite lobes; (2) are depolarized; and (3) are associated with optical emission line gas. The association between the HI absorption and the emission line gas, supports the hypothesis that the HI absorption is produced in the atomic cores of the emission line clouds. Our results are consistent with a picture in which compact steep spectrum sources interact with clouds of dense gas as they propagate through their host galaxy. We suggest that the asymmetries in the radio and optical emission can be due to interaction of a two sided radio source with an asymmetric distribution of dense clouds in their environment.

European VLBI Network spectral imaging of the "compact double" radio source 2050+364 in the UHF band at 1049 MHz has resolved the HI absorbing region, and has shown a faint continuum component to the North (N), in addition to the well-known East-West double (E, W). Re-examination of VLBI continuum images at multiple frequencies suggests that 2050+364 may well be a one-sided core-jet source, which appears as a double over a limited frequency range. One of the dominant features, W, would then be the innermost visible portion of the jet, and could be at or adjacent to the canonical radio core. The other, E, is probably related to shocks at a sudden bend of the jet, towards extended steep-spectrum region N. A remarkably deep and narrow HI absorption line component extends over the entire projected extent of 2050+364. It coincides in velocity with the [OIII] optical doublet lines to within 10 km/s. This HI absorption could arise in the atomic cores of NLR clouds, and the motion in the NLR is then remarkably coherent both along the line-of-sight and across a projected distance of > 300 pc on the plane of the sky. Broader, shallower HI absorption at lower velocities covers only the plausible core area W. This absorption could be due to gas which is either being entrained by the inner jet or is flowing out from the accretion region; it could be related to the BLR.

We have used VLBA fringe visibility data obtained at 15 GHz to examine the compact structure in 250 extragalactic radio sources. For 171 sources in our sample, more than half of the total flux density seen by the VLBA remains unresolved on the longest baselines. There are 163 sources in our list with a median correlated flux density at 15 GHz in excess of 0.5 Jy on the longest baselines. For about 60% of the sources, we have at least one observation in which the core component appears unresolved (generally smaller than 0.05 mas) in one direction, usually transverse to the direction into which the jet extends. BL Lacs are on average more compact than quasars, while active galaxies are on average less compact. Also, in an active galaxy the sub-milliarcsecond core component tends to be less dominant. IDV sources typically have a more compact, more core-dominated structure on sub-milliarcsecond scales than non-IDV sources, and sources with a greater amplitude of intra-day variations tend to have a greater unresolved VLBA flux density. The objects known to be GeV gamma-ray loud appear to have a more compact VLBA structure than the other sources in our sample. This suggests that the mechanisms for the production of gamma-ray emission and for the generation of compact radio synchrotron emitting features are related. The brightness temperature estimates and lower limits for the cores in our sample typically range between 10^11 and 10^13 K, but they extend up to 5x10^13 K, apparently in excess of the equipartition brightness temperature, or the inverse Compton limit for stationary synchrotron sources. The largest component speeds are observed in radio sources with high observed brightness temperatures, as would be expected from relativistic beaming (abridged).

We present a status report of dedicated circular polarization (CP) global VLBI observations of SS433. The total intensity image of the source is presented. We do not detect linear polarization in the core region of the source. We estimate the errors in D-terms for the EVN, and its effect on the final right-left circular polarization gain calibration. The data from this and earlier observations are analyzed with the "zero-V" self-calibration technique, in order to obtain indications for detectable CP in different regions of the source. Our data show no evidence for strong gyro-synchrotron (GS) emission on the size scale probed by our beam.

The radio source 4C12.50 has often been suggested to be a prime candidate for the link between ultraluminous infrared galaxies and young radio galaxies. A VLBI study of the neutral hydrogen in the nuclear regions of this object shows that most of the gas detected close to the systemic velocity is associated with an off-nuclear cloud (~50 to 100 pc from the radio core) with a column density of ~10^22 T_spin/100 K) cm^(-2) and an HI mass of a few times 10^5 to 10^6 M_sun. We consider a number of possibilities to explain the results. In particular, we discus the possibility that this cloud indicates the presence of a rich and clumpy interstellar medium in the centre, likely left over from the merger that triggered the activity and that this medium influences the growth of the radio source. The location of the cloud -- at the edge of the northern radio jet/lobe -- suggests that the radio jet might be interacting with a gas cloud. This interaction could be responsible for bending the young radio jet. The velocity profile of the gas is relatively broad (~150$ km/s) and we interpret this as kinematical evidence for interaction of the radio plasma with the cloud. We also consider the model where the cloud is part of a broader circumnuclear structure. Only a limited region of this structure would have sufficient background radio brightness and large enough column depth in neutral gas to obtain detectable HI absorption against the counterjet. The VLBI study of the neutral hydrogen in 4C12.50 suggests that HI detected near the systemic velocity (as it is often the case in radio galaxies) may not necessarily be connected with a circumnuclear disk or torus (as is very often assumed) but instead could be a tracer of the large-scale medium that surrounds the active nucleus and that may influence the growth of the young radio source.

We report the results of a 15 GHz (2 cm) multi-epoch VLBA program, begun in 1994 to study the outflow in radio jets ejected from quasars and active galaxies. The observed flow of 208 distinct features measured in 110 quasars, active galaxies, and BL Lac objects shows highly collimated relativistic motion with apparent transverse velocities typically between zero and about 15c, with a tail extending up to about 34c. Within individual jets, different features appear to move with a similar characteristic velocity which may represent an underlying continuous jet flow, but we also see some stationary and even apparently inward moving features which co-exist with the main features. Comparison of our observations with published data at other wavelengths suggests that there is a systematic decrease in apparent velocity with increasing wavelength, probably because the observations at different wavelengths sample different parts of the jet structure. The observed distribution of linear velocities is not consistent with any simple ballistic model. Either there is a rather broad range of Lorentz factors, a significant difference between the velocity of the bulk relativistic flow and the pattern speed of underlying shocks, or a combination of these options. Assuming a ballistic flow, comparison of observed apparent velocities and Doppler factors computed from the time scale of flux density variations is consistent with a steep power law distribution of intrinsic Lorentz factors, an isotropic distribution of orientations of the parent population, and intrinsic brightness temperatures about an order of magnitude below the canonical inverse Compton limit. It appears that the parent population of radio jets is not dominated by highly relativistic flows (abridged).

The idea that the radio jets in AGN contain material in relativistic motion is supported by many lines of observational evidence, including morphology, brightness temperature estimated with interferometers and with intrinsic variations, interstellar scintillations, X-rays, and superluminal motion. These are largely independent, and taken together make an irrefutable case for relativistic motion.

We investigate the relation between the Doppler factor determined from variations in total flux at 22 and 37 GHz, and the apparent transverse velocity determined from VLBA observations at 2 cm. The data are consistent with the relativistic beaming theory for compact radio sources, in that the distribution of beta_{app}/delta_{var}, for 30 quasars, is roughly consistent with a Monte Carlo simulation. The intrinsic temperature appears to be ~2x10^{10} K, close to the "equipartition value" calculated by Readhead (1994). We deduce the distribution of Lorentz factors for a group of 48 sources; the values range up to about gamma=40.

We report a change in the trajectory of a well studied jet component of the quasar 3C279. The component changes in apparent projected speed and direction, and we find it to be moving with a Lorentz factor gamma >~ 15 at an initial angle of <~ 1 degree to the line of sight. The new trajectory of the component has approximately the same speed and direction as an earlier superluminal feature, originally seen in the early 1970s. The new direction for the component is also much better aligned with larger scale VLBA and VLA structure out to 0.1 arcseconds. We suggest that the trajectory change is a collimation event occurring at >~ 1 kiloparsec (deprojected) along the jet. While the change in trajectory on the sky appears to be 26 degrees the intrinsic change is <~ 1 degree. We estimate the Doppler factor prior to the change in direction to be D >~ 28 and after the change to be D >~ 23. Comparison to independent constraints on the Doppler factor suggest that the energy in the radiating particles cannot greatly exceed the energy in the magnetic field unless the volume filling factor is very much less than one.

We consider the incidence of HI absorption in intrinsically small sub-galactic sized extragalactic sources selected from sources classified as Gigahertz Peaked Spectrum (GPS) and Compact Steep Spectrum (CSS) sources. We find that the smaller sources (<0.5 kpc) have larger HI column densities than the larger sources (>0.5 kpc). Both a spherical and an axi-symmetric gas distribution, with a radial power law density profile, can be used to explain this anti-correlation between projected linear size and HI column density. Since most detections occur in objects classified as galaxies, we argue that if the unified schemes apply to GPS/CSSs a disk distribution for the HI is more likely. The most favoured explanation for the compact sizes of the GPS/CSSs is that they are young sources evolving in a power law density medium. For the GPSs with measured expansion velocities, our derived densities are within an order of magnitude of those estimated from ram-pressure confinement of the lobes assuming equipartition. Our results therefore support the youth model.

We present an overview of the occurrence and properties of atomic gas associated with compact radio sources at redshifts up to z=0.85. Searches for HI 21cm absorption were made with the Westerbork Synthesis Radio Telescope at UHF-high frequencies (725-1200 MHz). Detections were obtained for 19 of the 57 sources with usable spectra (33%). We have found a large range in line depths, from tau=0.16 to tau<=0.001. There is a substantial variety of line profiles, including Gaussians of less than 10km/s, to more typically 150km/s, as well as irregular and multi-peaked absorption profiles, sometimes spanning several hundred km/s. Assuming uniform coverage of the entire radio source, we obtain column depths of atomic gas between 1e19 and 3.3e21(Tsp/100K)(1/f)cm^(-2). There is evidence for significant gas motions, but in contrast to earlier results at low redshift, there are many sources in which the HI velocity is substantially negative (up to v=-1420km/s) with respect to the optical redshift, suggesting that in these sources the atomic gas, rather than falling into the centre, may be be flowing out, interacting with the jets, or rotating around the nucleus.

The bright radio structure of the LINER elliptical galaxy NGC1052 is dominated by bi-symmetric jets on parsec scales. Features move outward on both sides of the core at v_{app}~0.26c. We have established the occurrence of free-free absorption, and suggest the presence of a patchy, geometrically region oriented roughly orthogonal to the jets; components have a wide range of spectral shapes and brightness temperatures. We distinguish three velocity systems in HI absorption. The ``high velocity'' system is the most prominent of these; it is redshifted from systemic by about 150 km/s. In HI VLBI it is seen towards both jets, but appears to be restricted to a shell 1 to 2 pc away from the core. The central hole might be largely ionised, and could be connected to the free-free absorption. WSRT spectroscopy shows 1667 and 1665 MHz OH main line absorption over at least the full ~250 km/s velocity range seen in HI. In the ``high velocity'' system, the profiles of the OH main lines and HI are similar, which suggests co-location of molecular and atomic gas. The OH satellite lines are also detected in the ``high velocity'' system: 1612 MHz in absorption and 1720 MHz in emission, with complementary strength. But we have no satisfactory model to explain all properties; the connection to H_2O masing gas at the same velocity but apparently a different location is also unclear.

In the Caltech-Jodrell Survey of bright flat-spectrum radio sources, 3-5 epochs have now been observed for nearly all 293 sources; a uniquely large sample. The derivation of component motions has not yet been completed; it is complicated by, for example, variability and jet curvature. Nevertheless, some basic results are clear. The average apparent velocity in CJF quasars is near 3c, while for radio galaxies and BL Lacs it is around 1c. The distribution of velocities is broad, indicating a broad distribution of jet Lorentz factors, with many low values, and/or a decoupling of the Lorentz factors between the cores and the moving jet components, due to bending, speed changes, or due to pattern motions (shocks). A hint has emerged at this workshop that the average apparent velocities might be lower at lower radio frequencies; this will need careful verification.

Combining our own observations with data from the literature, we consider the incidence of HI absorption in Gigahertz Peaked Spectrum (GPS) and Compact Steep Spectrum (CSS) sources. Here we present our preliminary results, where we find that the smaller GPS sources (<1 kpc) on average have larger HI column densities than the larger CSS sources (>1 kpc). Both a spherical and an axi-symmetric gas distribution, with a radial power law density profile, can be used to explain this anti-correlation between projected linear size and HI column density. Since most detections occur in galaxy classified objects, we argue that if the unified schemes apply to the GPS/CSSs, a disk distribution for the HI is more likely.

(Abridged) We discuss multiple VLBI continuum and spectral line observations and WSRT spectroscopy of NGC 1052. Sub-parsec scale features move outward at approximately 0.26c in bi-symmetric jets, most likely oriented near the plane of the sky. Absorption and emission signatures reveal ionised, atomic, and molecular components of the surrounding medium. Seven-frequency (1.4 to 43 GHz) VLBA observations show free-free absorption in the inner parsec, probably together with synchrotron self-absorption. There is apparently a geometrically thick but patchy structure oriented roughly orthogonal to the jets. The western jet is receding: it is covered more deeply and extensively. HI spectral line VLBI reveals atomic gas in front of both jets. There appear to be three velocity systems. The deepest, at "high velocities" (receding by 125 to 200 km/s), seems restricted to a shell 1 to 2 pc away from the core, within which this gas might be largely ionised. WSRT spectroscopy has revealed 1667 and 1665 MHz OH absorption with their line ratio varying roughly from 1:1 to 2:1 between -35 and 200 km/s. In the high velocity system the OH profiles are similar to HI, suggesting co-location of that atomic and molecular gas, and leaving unclear the connection to the H2O masing gas seen elsewhere. We have also detected both 18cm OH satellite lines in the high velocity system. They have conjugate profiles: 1612 MHz is in absorption, and 1720 MHz in emission.

High resolution radio observations remain the most direct way to study the formation and evolution of radio jets associated with the accretion onto massive black holes. We report preliminary results of our seven year VLBA observational program to understand the nature of relativistic beaming in blazars and the surrounding environment of massive black holes. Most blazars show an apparent outward flow away from an active core. However, in a few sources the motion appears inward, most likely the result of projection of a curved trajectory which bends back toward along the line of sight. The apparent motion of jet features is not always oriented along the direction separating the feature from the core, and in a few cases we have observed a clear change in the direction and velocity of a feature as it flows along the jet. In other sources, the motion appears to follow a simple ballistic trajectory. We find no simple relation between the time scales of flux density changes and apparent component velocities.

We present the results of a multi-epoch VLBA study of the powerful radio source 4C +12.50 (PKS 1345+125) at a wavelength of 2 cm. This compact radio source is associated with a hidden quasar whose host galaxy shows signs of a recent merger. It has been classified as a compact symmetric source (CSO) due to its small overall size (~220 pc) and twin-jet morphology, although it also has faint extended emission that may be a relic of previous activity. We report the detection of exceedingly high linear fractional polarization in isolated features of the southern jet (up to 60%), which is highly unusual for a CSO. Given the large amount of gas present in the host galaxy, we would expect significant Faraday depolarization across the whole source, unless the depolarizing gas is fairly clumpy. The southern jet also contains two features that are moving outward from the core at apparent speeds of v/c = 1.0 +/- 0.3 and 1.2 +/- 0.3. These represent the first positive detections of superluminal motion in a CSO, and taken together with the jet to counter-jet flux ratio, suggest an intrinsic flow speed of v/c = 0.84 +/- 0.12. The apparent ridgeline of the jet and counter-jet are consistent with a conical helix of wavelength 280 pc that is the result of Kelvin-Helmholtz instabilities driven by a slow precession of the jet nozzle. A fit to the data implies that the nozzle is precessing around a cone with half-angle 23 degrees, whose axis lies at an angle of 82 degrees to the line of sight. We suggest that the "S"-shaped jet morphologies commonly seen in recent AGN outflows such as 4C +12.50 may simply reflect the fact that their black hole spin axes are still precessing, and have not had sufficient time to align with their accretion disks.

This is a paper on young jet material in a frustratingly complex environment. NGC 1052 has a compact, flat- or GHz-peaked-spectrum radio nucleus consisting of bi-symmetric jets, oriented close to the plane of the sky. Many features on both sides move away at v_app~0.26c (H_0=65 km/s/Mpc). VLBI at seven frequencies shows a wide range of spectral shapes and brightness temperatures; there is clearly free-free absorption, probably together with synchrotron self-absorption, on both sides of the core. The absorbing structure is likely to be geometrically thick and oriented roughly orthogonal to the jets, but it is patchy. HI VLBI shows atomic gas in front of the approaching as well as the receding jet. There appear to be three velocity systems, at least two of which are local to the AGN environment. The "high velocity system", 125 to 200 km/s redward of systemic, seems restricted to a shell 1 to 2 pc away from the core. Closer to the centre, this gas might be largely ionised; it could cause the free-free absorption. WSRT spectroscopy shows 1667 and 1665 MHz OH absorption over a wide velocity range. OH and HI profile similarity suggests co-location of molecular and atomic ``high velocity'' gas; the connection to H2O masing gas is unclear. Further, at ``high velocity'' we detected the OH 1612 MHz satellite line in absorption and the 1720 MHz line in emission, with complementary strengths.

We summarize the results of VLBA and global VLBI observations of SS433 between 1995 and 2000. With these observations we resolve the inner jet of the source and identify an absorption region (~25 AU), the "central radio gap". The radio gap is caused by free-free absorption of the jet radio emission by a flattened outflow from the binary system. Radio emission is detected on 100 AU scales perpendicular to the normal jets ("equatorial emission region"; "equatorial outflow"). At some epochs the emission is smooth but compact features are frequently detected. We suggest that equatorial outflows may be common in microquasars.

We present VLBI imaging results for SS433 from the period 1995-2000. An equatorial emission region is detected and confirmed in 7 experiments at 1.6 and 5 GHz. The structure of the region changes with time on scales of weeks to months. The nature of emission is unknown, but it is certainly non-thermal, maybe optically thin synchrotron. The spectrum however might be attenuated by a thermal population of electrons. There were previous indications for an equatorial outflow from the system. We do detect outward motion of a radio component, and estimate a speed of 1200/sin(i) km/s. Theoretical work and IR observations suggest that the mass-loss rate in SS433 is much higher in this region than in the jets. We suggest that the high-mass microquasars might all have equatorial outflows; in fact, it is the dominant form of mass-loss in these systems.

We are investigating the kinematics of jets in active galactic nuclei on parsec scales by studying a representative population of sources. This study is being carried out using the Very Long Baseline Array at 15 GHz, with more than 800 images taken since 1994. In this contribution we present an overview of the diversity of kinematics for a complete sample of sources.

We report further results from our imaging survey at 15GHz (lambda=2 cm) with the Very Long Baseline Array. This paper presents single epoch images for 39 sources, bringing the total number of objects in the sample to 171. Our sample is representative of a complete unbiased sample and it will be used for statistical analysis of source properties. We compare the observed brightness temperatures derived from our VLBA observations to those derived from total intensity variations at 22 and 37 GHz. These are consistent with intrinsic brightness temperatures in the range 10^10 to 10^12 K. We also present three new spectroscopic redshift values: z=0.517+/-0.001 for 0026+346, z=1.591+/-0.003 for 0727-115, and z=0.2016+/-0.0004 for 1155+251. Images from this VLBA 2 cm survey are available on the Internet under http://www.cv.nrao.edu/2cmsurvey.

We present the results from our deep (16x12 hour) Westerbork Synthesis Radio Telescope (WSRT) observations of the approximately 7 square degree Bootes Deep Field, centered at 14h32m05.75s, 34d16'47.5" (J2000). Our survey consists of 42 discrete pointings, with enough overlap to ensure a uniform sensitivity across the entire field, with a limiting sensitivity of 28 microJy (1 sigma rms). The catalog contains 3172 distinct sources, of which 316 are resolved by the 13"x27" beam. The Bootes field is part of the optical/near-infrared imaging and spectroscopy survey effort conducted at various institutions. The combination of these data sets, and the deep nature of the radio observations will allow unique studies of a large range of topics including the redshift evolution of the luminosity function of radio sources, the K-z relation and the clustering environment of radio galaxies, the radio / far-infrared correlation for distant starbursts, and the nature of obscured radio loud AGN.

Westerbork Radio Synthesis Telescope observations of the redshifted 21cm line absorber against the z_em = 0.871 double lobed quasar 3C196 show that the intervening absorber is an L approximately equal L* spiral galaxy (3C196-G1) and that the absorbing layer of cold gas extends to radii of at least 30 h_50^-1 kpc. The new data solve several long standing puzzles about this system by (1) discovering a second 21cm absorption feature, corresponding to absorption against the NE lobe of the background radio source and (2) spatially ``resolving'' the two absorption features to isolate the absorption along the two lines of sight to the opposing radio lobes. These findings resolve the disagreement in redshift between the UV metal and 21cm lines, and as well as demonstrating that the neutral layer does absorb both lobes of the background radio source. Simple kinematic models with an inclined, rotating gas disk match the observed 21cm profile and are also compatible with both the redshift and velocity spread of the absorption measured in UV resonance lines along a third, independent line of sight to the quasar nucleus and with the lack of 21cm absorption in as earlier VLBI experiment that was sensitive to opacity against the hot spot in the northern lobe. The inferred rotation speed and luminosity for the galaxy are compatible with the z approximately 0 Tully-Fisher Relation. This system illustrates well how 21cm absorption against extended background radio sources is a powerful tool in determining the nature of the damped Lyman-alpha class of QSO absorption line system.

Using Very Long Baseline Interferometry we have searched a sample of 300 compact radio sources for examples of multiple imaging produced by gravitational lensing; no multiple images were found with separations in the angular range 1.5--50 milliarcsec. This null result allows us to place a limit on the cosmological abundance of intergalactic supermassive compact objects in the mass range $\sim 10^{6}$ to $\sim 10^{8}$M$_{\odot}$; such objects cannot make up more than $\sim 1%$ of the closure density (95% confidence). A uniformly distributed population of supermassive black holes forming soon after the Big Bang do not, therefore, contribute significantly to the dark matter content of the Universe.

We present measurements of absorption by the 21cm hyperfine transition of neutral hydrogen toward radio sources at substantial look-back times. These data are used in combination with observations of rotational transitions of common interstellar molecules, to set limits on the evolution of the fine structure constant: (dot alpha)/alpha < 3.5e-15/yr to a look-back time of 4.8 Gyr. The neutral hydrogen observations employed Very Long Baseline Interferometry in order to mitigate the substantial uncertainty arising from the fact that observations at very different wavelengths may probe different lines-of-site due to frequency dependent structure of the background source. We discuss the implication of these results on theories unifying natural forces based on compact quantum dimensions. In the context of string theory, the limit on the secular evolution of the scale factor of these compact dimensions, R, is (dot R)/R < 1e-15/yr. Including terrestrial and other astronomical measurements places limits (2sigma) on slow oscillations of R from the present to the epoch of cosmic nucleosynthesis, just seconds after the big bang, of (Delta R)/R < 1e-5.

WSRT observations have provided a first inventory of the incidence of HI 21 cm line absorption associated with AGN at redshifts up to z=1.0. There is a large range in line depths, from tau=0.44 to tau less than 0.001, and a substantial variety of line profiles, from Gaussians of less than ten km/s to more typically a few hundred km/s, as well as irregular and multi-peaked absorption, sometimes spanning many hundreds of km/s. The chance to detect appreciable HI absorption is greatest in the most compact radio sources, GPSs and CSOs, where it can occur in circumnuclear ``disks'' or ``tori'', as well as in gas enveloping jets and hot spots; inferred densities range at least between 10 per cubic cm and 10000 per cubic cm. But HI absorption occurs also in some CSSs, perhaps associated with jet-cloud interaction regions, and in quasars with a large optical reddening. VLBI observations at the unusual UHF frequencies of redshifted HI 21 cm can give unique ``sight''lines into the physics and evolution of young radio sources and their inner galactic medium.

Compact Symmetric Objects (CSOs) are young radio galaxies whose jet axes lie close to the plane of the sky, and whose appearance is therefore not dominated by relativistic beaming effects. The small linear sizes of CSOs make them valuable for studies of both the evolution of radio galaxies and testing unified schemes for active galactic nuclei (AGN). A parsec-scale region of gas surrounding the central engine is predicted by both accretion and obscuration scenarios. Working surfaces, or ``hot spots,'' and the radio jets of CSOs are close enough to the central engines that this circumnuclear gas can be seen in absorption. The CSOs Observed in the Northern Sky (COINS) sample is comprised of 52 CSO candidates identified in three VLBI surveys. Of these, 27 have now been confirmed as CSOs. Optical redshifts are available in the literature for 28 of the CSO candidates, and HI absorption has been detected toward four. We present new optical spectroscopic redshifts for three of the candidates and summarize the current status of optical identifications. We further report on the discovery of HI in absorption towards the CSO J1816+3457 and summarize the results of neutral hydrogen absorption studies of the sources in this sample.

Nearly simultaneous VLBA observations of 3C 84, the radio source associated with NGC 1275, have been made at multiple frequencies to study the free-free absorption of the northern, or "counterjet", feature found by Walker, Romney, and Benson (1994, Ap.J. 430, L45) and by Vermeulen, Readhead, and Backer (1994, Ap.J. 430, L41). Our observations confirm that the spectra are consistent with free-free absorption and eliminate the possibility that the earlier result was an effect of variability. The northern feature is well resolved spatially, so images have been made showing the distribution of the absorption over a region of about 1.5 pc on a side, beginning about 1.5 pc from the presumed location of the central object. That distribution is dominated by a strong decrease with radial distance. The magnitude of the absoption near 2.5 pc projected distance from the central object is consistent with a 10^4 K gas with an emission measure of about 5 X 10^8 pc cm^-6. The geometry is consistent with absorption by ionized gas associated with an accretion disk. The data provide firm constraints for models of the outer regions of accretion disks and, perhaps, associated winds.

We present multi-frequency VLBA+VLA observations of SS433 at 1.6, 5 and 15 GHz. These observations provide the highest angular resolution radio spectral index maps ever made for this object. Motion of the components of SS433 during the observation is detected. In addition to the usual VLBI jet structure, we detect two radio components in the system at an anomalous position angle. These newly discovered radio emitting regions might be related to a wind-like equatorial outflow or to an extension of the accretion disk. We show that the radio core component is bifurcated with a clear gap between the eastern and western wings of emission. Modelfitting of the precessing jets and the moving knots of SS433 shows that the kinematic centre -- i.e. the binary -- is in the gap between the western and eastern radio core components. Spectral properties and observed core position shifts suggest that we see a combined effect of synchrotron self-absorption and external free-free absorption in the innermost AU-scale region of the source. The spatial distribution of the ionized matter is probably not spherically symmetric around the binary, but could be disk-like.

We observed Sgr A* and two extragalactic radio sources nearby in angle with the VLBA over a period of two years and measured relative positions with an accuracy approaching 0.1 mas. The apparent proper motion of Sgr A* relative to J1745-283 is 5.90 +/- 0.4 mas/yr, almost entirely in the plane of the Galaxy. The effects of the orbit of the Sun around the Galactic Center can account for this motion, and any residual proper motion of Sgr A*, with respect to extragalactic sources, is less than about 20 km/s. Assuming that Sgr A* is at rest at the center of the Galaxy, we estimate that the circular rotation speed in the Galaxy at the position of the Sun is 219 +/- 20 km/s, scaled by Ro/8.0 kpc. Current observations are consistent with Sgr A* containing all of the nearly 2.6 x 10^6 solar masses, deduced from stellar proper motions, in the form of a massive black hole. While the low luminosity of Sgr A*, for example, might possibly have come from a contact binary containing of order 10 solar masses, the lack of substantial motion rules out a "stellar" origin for Sgr A*. The very slow speed of Sgr A* yields a lower limit to the mass of Sgr A* of about 1,000 solar masses. Even for this mass, Sgr A* appears to be radiating at less than 0.1 percent of its Eddington limit.

We have used the VLBA at 15 GHz to image the structure of 132 strong compact AGN and quasars with a resolution better than one milliarcsecond and a dynamic range typically exceeding 1000 to 1. These observations were made as part of a program to investigate the sub-parsec structure of quasars and AGN and to study the changes in their structure with time. Many of the sources included in our study, particularly those located south of +35 degrees, have not been previously imaged with milliarcsecond resolution. Each of the sources has been observed at multiple epochs. In this paper we show images of each of the 132 sources which we have observed. For each source we present data at the epoch which had the best quality data. The milliarcsecond jets generally appear one-sided but two-sided structure is often found in lower luminosity radio galaxies and in high luminosity quasars with gigahertz peaked spectra. Usually the structure is unresolved along the direction perpendicular to the jet, but a few sources have broad plumes. In some low luminosity radio galaxies, the structure appears more symmetric at 2 cm than at long wavelengths. The apparent long wavelength symmetry in these sources is probably due to absorption by intervening material. A few sources contain only a single component with any secondary feature at least a thousand times weaker. We find no obvious correlation of radio morphology and the detection of gamma-ray emission by EGRET.

The Caltech-Jodrell Bank VLBI surveys of bright extragalactic radio sources north of declination 35 deg were carried out between 1990 and 1995 using the Mark-II system, achieving images with a resolution of about 1 mas at 5 GHz. The CJ1 survey (together with the older `PR' survey) includes 200 objects with 5 GHz flux density greater than 0.7 Jy; the CJ2 survey includes 193 flat-spectrum sources with 5 GHz flux density greater than 0.35 Jy; and we have defined a complete flux-density limited sample, CJF, of 293 flat-spectrum sources stronger than 0.35 Jy. We summarize the definition of the samples and the VLBI, VLA, MERLIN, and optical observations, and present some highlights of the astrophysical results. These include: (1) superluminal motion and cosmology; (2) morphology and evolution of the `compact symmetric objects' (CSOs); (3) two-sided motion in some CSOs; (4) the angular-size - redshift diagram; (5) misalignment of parsec-scale and kiloparsec-scale jets.

We present measurements of bi-directional motions in the jets of the radio galaxy 1946+708 at z=0.101. This is a Compact Symmetric Object with striking S-symmetry. Sensitive 15 GHz observations reveal a compact component at the center of symmetry with a strongly inverted spectrum, that we identify as the core. From five 4.9 GHz observations spread over 4 years we have determined the velocities of four compact jet components. If simple kinematic models can be applied then the inclination of the source and the bulk jet velocity can be directly determined for any assumed value of the Hubble constant. Conversely, the measurements already place constraints on the Hubble constant, and we show how further observations of 1946+708 can yield an increasingly accurate determination of H_0.

Polarimetry of the narrow-line radio galaxy PKS 0116+082 at the W.M. Keck telescope shows that it has high and variable optical polarization, presumably due to synchrotron radiation. It is not a BL Lac object because it has strong narrow lines, and it is not an OVV quasar because it has no broad lines and the extended galaxy is prominent. VLA and VLBA images show that it is a compact steep-spectrum radio source with most of the emission coming from a region less than 100 milli-arcsec in size. Of the 25 compact steep-spectrum or gigahertz peaked-spectrum sources measured polarimetrically, four have high optical polarization. One of these has been observed only once but in the other three the polarization is variable. This gives an intriguing hint that variability may be a general property of these objects.