Search

Dark Matter (DM) halos with masses below $\sim10^{8}$ $M_{\odot}$, which would help to discriminate between DM models, may be detected through their gravitational effect on distant sources. The same applies to primordial black holes, considered as an alternative scenario to DM particle models. However, there is still no evidence for the existence of such objects. With the aim of finding compact objects in the mass range $\sim$ 10$^{6}$ -- 10$^{9}$$M_{\odot}$, we search for strong gravitational lenses on milli (mas)-arcseconds scales (< 150 mas). For our search, we used the Astrogeo VLBI FITS image database -- the largest publicly available database, containing multi-frequency VLBI data of 13828 individual sources. We used the citizen science approach to visually inspect all sources in all available frequencies in search for images with multiple compact components on mas-scales. At the final stage, sources were excluded based on the surface brightness preservation criterion. We obtained a sample of 40 sources that passed all steps and therefore are judged to be milli-arcsecond lens candidates. These sources are currently followed-up with on-going European VLBI Network (EVN) observations at 5 and 22 GHz. Based on spectral index measurements, we suggest that two of our candidates have a higher probability to be associated with gravitational lenses., Comment: 6 pages, 1 figure; accepted for publication

We find evidence for large-scale clustering amongst Fermi-selected BL Lac objects but not amongst Fermi-selected FSRQs. Using two-point correlation functions we have investigated the clustering properties of different classes of objects from the Fermi LAT 4FGL catalogue. We wanted to test the idea based on optical polarization observations that there might be large volumes of space in which AGN axes are aligned. To do this we needed a clean sample of blazars as these are objects with their jet axes pointing towards the observer and Fermi sources provide such a sample. We find that high latitude Fermi sources taken as a whole show a significant clustering signal on scales up to 30 degrees. To investigate if all blazars behave in the same way we used the machine learning classifications of Kovacevic, et al. (2020), which are based only on gamma-ray information, to separate BL Lac-like objects from FSRQ-like objects. A possible explanation for the clustering signal we find amongst the BL Lac-like objects is that there are indeed large volumes of space in which AGN axes are aligned. This signal might be washed out in FSRQs since they occupy a much larger volume of space. Thus our results support the idea that large scale polarization alignments could originate from coherent alignments of AGN axes. We speculate that these axis alignments may be related to the well-known intrinsic alignments of galaxy optical position angles.

PKS 1413+135 is one of the most peculiar blazars known. Its strange properties led to the hypothesis almost four decades ago that it is gravitationally lensed by a mass concentration associated with an intervening galaxy. It exhibits symmetric achromatic variability, a rare form of variability that has been attributed to gravitational milli-lensing. It has been classified as a BL Lac object, and is one of the rare objects in this class with a visible counterjet. BL Lac objects have jet axes aligned close to the line of sight. It has also been classified as a compact symmetric object, which have jet axes not aligned close to the line of sight. Intensive efforts to understand this blazar have hitherto failed to resolve even the questions of the orientation of the relativistic jet, and the host galaxy. Answering these two questions is important as they challenge our understanding of jets in active galactic nuclei and the classification schemes we use to describe them. We show that the jet axis is aligned close to the line of sight and PKS 1413+135 is almost certainly not located in the apparent host galaxy, but is a background object in the redshift range $0.247 < z < 0.5$. The intervening spiral galaxy at $z = 0.247$ provides a natural host for the putative lens responsible for symmetric achromatic variability and is shown to be a Seyfert 2 galaxy. We also show that, as for the radio emission, a "multizone" model is needed to account for the high-energy emission., Comment: Accepted for publication in The Astrophysical Journal. 28 pages, 15 figures

The radio source at the center of the galaxy cluster Abell 2626, also known as the Kite, stands out for its unique morphology composed of four, symmetric arcs. Previous studies have probed the properties of this source at different frequencies and its interplay with the surrounding thermal plasma, but the puzzle of its origin is still unsolved. We use new LOw Frequency ARray (LOFAR) observation from the LOFAR Two-meter Sky Survey at 144 MHz to investigate the origin of the Kite.} We present a detailed analysis of the new radio data which we combined with archival radio and X-ray observations. We have produced a new, resolved spectral index map of the source with a resolution of 7$''$ and we studied the spatial correlation of radio and X-ray emission to investigate the interplay between thermal and non-thermal plasma. The new LOFAR data have changed our view of the Kite by discovering two steep-spectrum ($\alpha<-1.5$) plumes of emission connected to the arcs. The spectral analysis shows, for the first time, a spatial trend of the spectrum along the arcs with evidence of curved synchrotron spectra and a spatial correlation with the X-ray surface brightness. On the basis of our results, we propose that the Kite was originally an X-shaped radio galaxy whose fossil radio plasma, after the end of the activity of the central active galactic nucleus, has been compressed due to motions of the thermal plasma in which it is encompassed. The interplay between the compression and advection of the fossil plasma, with the restarting of the nuclear activity of the central galaxy, could have enhanced the radio emission of the fossil plasma producing the arcs of the Kite. We present also the first, low-frequency observation of a jellyfish galaxy in the same field, in which we detect extended, low-frequency emission without a counterpart at higher frequencies., Comment: 10 pages, 7 figures. Accepted for publication on Astronomy & Astrophysics on 23/09/2020

The Baryon acoustic oscillations from Integrated Neutral Gas Observations (BINGO) telescope is a 40-m~class radio telescope under construction that has been designed to measure the large-angular-scale intensity of HI emission at 980--1260 MHz and hence to constrain dark energy parameters. A large focal plane array comprising of 1.7-metre diameter, 4.3-metre length corrugated feed horns is required in order to optimally illuminate the telescope. Additionally, very clean beams with low sidelobes across a broad frequency range are required, in order to facilitate the separation of the faint HI emission from bright Galactic foreground emission. Using novel construction methods, a full-sized prototype horn has been assembled. It has an average insertion loss of around 0.15 dB across the band, with a return loss around -25 dB. The main beam is Gaussian with the first sidelobe at around $-25 dB. A septum polariser to separate the signal into the two hands of circular polarization has also been designed, built and tested., Comment: This is a pre-print of an article published in Experimental Astronomy. The final authenticated version is available online at: https://doi.org/10.1007/s10686-020-09666-9

The Baryon acoustic oscillations from Integrated Neutral Gas Observations (BINGO) telescope is a new 40-m class radio telescope to measure the large-angular-scale intensity of Hi emission at 980-1260 MHz to constrain dark energy parameters. As it needs to measure faint cosmological signals at the milliKelvin level, it requires a site that has very low radio frequency interference (RFI) at frequencies around 1 GHz. We report on measurement campaigns across Uruguay and Brazil to find a suitable site, which looked at the strength of the mobile phone signals and other radio transmissions, the location of wind turbines, and also included mapping airplane flight paths. The site chosen for the BINGO telescope is a valley at Serra do Urubu, a remote part of Paraiba in North-East Brazil, which has sheltering terrain. During our measurements with a portable receiver we did not detect any RFI in or near the BINGO band, given the sensitivity of the equipment. A radio quiet zone around the selected site has been requested to the Brazilian authorities ahead of the telescope construction., Comment: Preprint of an article accepted in the Journal of Astronomical Instrumentation, copyright 2018 World Scientific Publishing Company https://www.worldscientific.com/worldscinet/jai

We have reanalysed the 1996/1997 VLA monitoring data of the gravitational lens system JVAS B0218+357 to produce improved total flux density and polarization variability curves at 15, 8.4 and 5 GHz. This has been done using improved calibration techniques, accurate subtraction of the emission from the Einstein ring and careful correction of various systematic effects, especially an offset in polarization position angle that is hour-angle dependent. The variations in total and polarized flux density give the best constraints and we determine a combined delay estimate of $11.3 \pm 0.2$ d (1$\sigma$). This is consistent with the $\gamma$-ray value recently derived using the Fermi Gamma-ray Space Telescope and thus we find no evidence for a positional shift between the radio and $\gamma$-ray emitting regions. Combined with the previously published lens model found using LensClean, the new delay gives a value for the Hubble constant of $H_0 = 72.9 \pm 2.6$ km s$^{-1}$ Mpc$^{-1}$ (1$\sigma$)., Comment: 17 pages and 14 figures. Accepted for publication in MNRAS

HI intensity mapping is an emerging tool to probe dark energy. Observations of the redshifted HI signal will be contaminated by instrumental noise, atmospheric and Galactic foregrounds. The latter is expected to be four orders of magnitude brighter than the HI emission we wish to detect. We present a simulation of single-dish observations including an instrumental noise model with 1/f and white noise, and sky emission with a diffuse Galactic foreground and HI emission. We consider two foreground cleaning methods: spectral parametric fitting and principal component analysis. For a smooth frequency spectrum of the foreground and instrumental effects, we find that the parametric fitting method provides residuals that are still contaminated by foreground and 1/f noise, but the principal component analysis can remove this contamination down to the thermal noise level. This method is robust for a range of different models of foreground and noise, and so constitutes a promising way to recover the HI signal from the data. However, it induces a leakage of the cosmological signal into the subtracted foreground of around 5%. The efficiency of the component separation methods depends heavily on the smoothness of the frequency spectrum of the foreground and the 1/f noise. We find that as, long as the spectral variations over the band are slow compared to the channel width, the foreground cleaning method still works., Comment: 14 pages, 12 figures. Submitted to MNRAS

We outline the prospects for performing pioneering radio weak gravitational lensing analyses using observations from a potential forthcoming JVLA Sky Survey program. A large-scale survey with the JVLA can offer interesting and unique opportunities for performing weak lensing studies in the radio band, a field which has until now been the preserve of optical telescopes. In particular, the JVLA has the capacity for large, deep radio surveys with relatively high angular resolution, which are the key characteristics required for a successful weak lensing study. We highlight the potential advantages and unique aspects of performing weak lensing in the radio band. In particular, the inclusion of continuum polarisation information can greatly reduce noise in weak lensing reconstructions and can also remove the effects of intrinsic galaxy alignments, the key astrophysical systematic effect that limits weak lensing at all wavelengths. We identify a VLASS "deep fields" program (total area ~10-20 square degs), to be conducted at L-band and with high-resolution (A-array configuration), as the optimal survey strategy from the point of view of weak lensing science. Such a survey will build on the unique strengths of the JVLA and will remain unsurpassed in terms of its combination of resolution and sensitivity until the advent of the Square Kilometre Array. We identify the best fields on the JVLA-accessible sky from the point of view of overlapping with existing deep optical and near infra-red data which will provide crucial redshift information and facilitate a host of additional compelling multi-wavelength science., Comment: Submitted in response to NRAO's recent call for community white papers on the VLA Sky Survey (VLASS)

Unified schemes of radio sources, which account for different types of radio AGN in terms of anisotropic radio and optical emission, together with different orientations of the ejection axis to the line of sight, have been invoked for many years. Recently, large samples of optical quasars, mainly from the Sloan Digital Sky Survey, together with large radio samples, such as FIRST, have become available. These hold the promise of providing more stringent tests of unified schemes but, compared to previous samples, lack high resolution radio maps. Nevertheless they have been used to investigate unified schemes, in some cases yielding results which appear inconsistent with such theories. Here we investigate using simulations how the selection effects to which such investigations are subject can influence the conclusions drawn. In particular, we find that the effects of limited resolution do not allow core-dominated radio sources to be fully represented in the samples, that the effects of limited sensitivity systematically exclude some classes of sources and the lack of deep radio data make it difficult to decide to what extent closely separated radio sources are associated. Nevertheless, we conclude that relativistic unified schemes are entirely compatible with the current observational data. For a sample selected from SDSS and FIRST which includes weak-cored triples we find that the equivalent width of the [OIII] emission line decreases as core-dominance increases, as expected, and also that core-dominated quasars are optically brighter than weak-cored quasars., Comment: 12 pages, 9 figures: accepted by MNRAS

BINGO is a concept for performing a 21cm intensity mapping survey using a single dish telescope. We briefly discuss the idea of intensity mapping and go on to define our single dish concept. This involves a \sim 40 m dish with an array of \sim 50 feed horns placed \sim 90 m above the dish using a pseudo-correlation detection system based on room temperature LNAs and one of the celestial poles as references. We discuss how such an array operating between 960 and 1260 MHz could be used to measure the acoustic scale to 2.4% over the redshift range 0.13

We discuss the detection of large scale HI intensity fluctuations using a single dish approach with the ultimate objective of measuring the Baryonic Acoustic Oscillations and constraining the properties of dark energy. We present 3D power spectra, 2D angular power spectra for individual redshift slices, and also individual line-of-sight spectra, computed using the S^3 simulated HI catalogue which is based on the Millennium Simulation. We consider optimal instrument design and survey strategies for a single dish observation at low and high redshift for a fixed sensitivity. For a survey corresponding to an instrument with T_sys=50 K, 50 feed horns and 1 year of observations, we find that at low redshift (z \approx 0.3), a resolution of 40 arc min and a survey of 5000 deg^2 is close to optimal, whereas at higher redshift (z \approx 0.9) a resolution of 10 arcmin and 500 deg^2 would be necessary. Continuum foreground emission from the Galaxy and extragalactic radio sources are potentially a problem. We suggest that it could be that the dominant extragalactic foreground comes from the clustering of very weak sources. We assess its amplitude and discuss ways by which it might be mitigated. We then introduce our concept for a single dish telescope designed to detect BAO at low redshifts. It involves an under-illumintated static 40 m dish and a 60 element receiver array held 90 m above the under-illuminated dish. Correlation receivers will be used with each main science beam referenced against an antenna pointing at one of the Celestial Poles for stability and control of systematics. We make sensitivity estimates for our proposed system and projections for the uncertainties on the power spectrum after 1 year of observations. We find that it is possible to measure the acoustic scale at z\approx 0.3 with an accuracy 2.4% and that w can be measured to an accuracy of 16%., Comment: 20 pages, 20 figures, submitted to MNRAS

The Fermi-LAT survey provides a large sample of blazars selected on the strength of their inverse Compton emission. We cross-correlate the first Fermi-LAT catalogue with the CRATES radio catalogue and use this sample to investigate whether blazar gamma-ray luminosities are influenced by the availability of external photons to be up-scattered. Using the 8.4 GHz flux densities of their compact radio cores as a proxy for their jet power, we calculate their Compton Efficiency parameters, which measure the ability of jets to convert power in the form of ultra-relativistic electrons into Compton gamma-rays. We find no clear differences in Compton efficiencies between BL Lac objects and FSRQs and no anti-correlation between Compton efficiency and synchrotron peak frequency. This suggests that the scattering of external photons is energetically unimportant compared to the synchrotron self-Compton process. These results contradict the predictions of the blazar sequence., Comment: 6 pages, 5 figures. Submitted to MNRAS

The Hubble constant value is currently known to 10% accuracy unless assumptions are made for the cosmology (Sandage et al. 2006). Gravitational lens systems provide another probe of the Hubble constant using time delay measurements. However, current investigations of ~20 time delay lenses, albeit of varying levels of sophistication, have resulted in different values of the Hubble constant ranging from 50-80 km/s/Mpc. In order to reduce uncertainties, more time delay measurements are essential together with better determined mass models (Oguri 2007, Saha et al. 2006). We propose a more efficient technique for measuring time delays which does not require regular monitoring with a high-resolution interferometer array. The method uses double image and long-axis quadruple lens systems in which the brighter component varies first and dominates the total flux density. Monitoring the total flux density with low-resolution but high sensitivity radio telescopes provides the variation of the brighter image and is used to trigger high-resolution observations which can then be used to see the variation in the fainter image. We present simulations of this method together with a pilot project using the WSRT (Westerbork Radio Synthesis Telescope) to trigger VLA (Very Large Array) observations. This new method is promising for measuring time delays because it uses relatively small amounts of time on high-resolution telescopes. This will be important because many SKA pathfinder telescopes, such as MeerKAT (Karoo Array Telescope) and ASKAP (Australian Square Kilometre Array Pathfinder), have high sensitivity but limited resolution., Comment: 4 pages, 6 figures, Contribution to the 10th European VLBI Network Symposium and EVN Users Meeting: VLBI and the new generation of radio arrays

Knowledge of the population of radio sources in the range ~2-200 GHz is important for understanding their effects on measurements of the Cosmic Microwave Background power spectrum. We report measurements of the 30 GHz flux densities of 605 radio sources from the Combined Radio All-sky Targeted Eight-GHz Survey (CRATES), which have been made with the One Centimetre Receiver Array prototype (OCRA-p) on the Torun 32-m telescope. The flux densities of sources that were also observed by WMAP and previous OCRA surveys are in broad agreement with those reported here, however a number of sources display intrinsic variability. We find a good correlation between the 30 GHz and Fermi gamma-ray flux densities for common sources. We examine the radio spectra of all observed sources and report a number of Gigahertz-peaked and inverted spectrum sources. These measurements will be useful for comparison to those from the Low Frequency Instrument of the Planck satellite, which will make some of its most sensitive observations in the region covered here., Comment: 21 pages (9 pages of text, 12 pages of table), 7 figures. Erratum appended to end (page 20). Accepted by MNRAS. The definitive version is available at www.blackwell-synergy.com

We have studied the implications of high sensitivity polarization measurements of objects from the WMAP point source catalogue made using the VLA at 8.4, 22 and 43 GHz. The fractional polarization of sources is almost independent of frequency with a median of ~2 per cent and an average, for detected sources, of ~3.5 per cent. These values are also independent of the total intensity over the narrow range of intensity we sample. Using a contemporaneous sample of 105 sources detected at all 3 VLA frequencies, we have investigated the spectral behaviour as a function of frequency by means of a 2-colour diagram. Most sources have power-law spectra in total intensity, as expected. On the other hand they appear to be almost randomly distributed in the polarized intensity 2-colour diagram. This is compatible with the polarized spectra being much less smooth than those in intensity and we speculate on the physical origins of this. We have performed an analysis of the correlations between the fractional polarization and spectral indices including computation of the principal components. We find that there is little correlation between the fractional polarization and the intensity spectral indices. This is also the case when we include polarization measurements at 1.4 GHz from the NVSS. In addition we compute 45 rotation measures from polarization position angles which are compatible with a \lambda^2 law. We use our results to predict the level of point source confusion noise that contaminates CMB polarization measurements aimed at detecting primordial gravitational waves from inflation. We conclude that some level of source subtraction will be necessary to detect r~0.1 below 100 GHz and at all frequencies to detect r~0.01. We present estimates of the level of contamination expected and the number of sources which need to be subtracted as a function of the imposed cut flux density and frequency., Comment: 21 pages, 17 figures. MNRAS (accepted)

We have investigated the correlations between optical and radio isophotal position angles for 14302 SDSS galaxies with r magnitudes brighter than 18. All the galaxies are identified with extended FIRST radio sources. For passive early-type galaxies, which we distinguish from the others by using the colour, concentration and their principal components, we find a strong statistical alignment of the radio axes with the optical minor axes. Since the radio emission is driven by accretion on to a nuclear black hole we argue that the observed correlation gives new insight into the structure of elliptical galaxies, for example, whether or not the nuclear kinematics are decoupled from the rest of the galaxy. Our results imply that a significant fraction of the galaxies are oblate spheroids, perhaps rotationally supported, with their radio emission aligned with the stellar minor axis. Remarkably, the strength of the correlation of the radio major axis with the optical minor axis depends on radio loudness. Those objects with a low ratio of FIRST radio flux density to total stellar light show a strong minor axis correlation while the stronger radio sources do not. This split may reflect different formation histories and we suggest this may be a new manifestation of the better known dichotomy between slow rotating and fast rotating ellipticals., Comment: To be published in the proceedings of the conference "Accretion and Ejection in AGN: A global view, June 22-26 2009 - Como, Italy"

We present polarization measurements at 8.4, 22, and 43 GHz made with the VLA of a complete sample of extragalactic sources stronger than 1 Jy in the 5-year WMAP catalogue and with declinations north of -34 degrees. The observations were motivated by the need to know the polarization properties of radio sources at frequencies of tens of GHz in order to subtract polarized foregrounds for future sensitive Cosmic Microwave Background (CMB) experiments. The total intensity and polarization measurements are generally consistent with comparable VLA calibration measurements for less-variable sources, and within a similar range to WMAP fluxes for unresolved sources. A further paper will present correlations between measured parameters and derive implications for CMB measurements., Comment: 16 pages, 5 figures. Accepted by MNRAS. Two large figures are available at http://www.jb.man.ac.uk/~njj/pol_maps and http://www.jb.man.ac.uk/~njj/pol_spec

We present the first results from the new Survey of Extragalactic Nuclear Spectral Energies (SENSE) sample of "blazars". The sample has been chosen with minimal selection effects and is therefore ideal to probe the intrinsic properties of the blazar population. We report evidence for negative cosmological evolution in this radio selected sample and give an outline of future work related to the SENSE sample., Comment: 2 pages, to appear in the proceedings of "Accretion and ejection in AGN: a global view" (Como, 22-26 June 2009)

Small angular scale (high l) studies of cosmic microwave background anisotropies require accurate knowledge of the statistical properties of extragalactic sources at cm-mm wavelengths. We have used a 30 GHz dual-beam receiver (OCRA-p) on the Torun 32-m telescope to measure the flux densities of 121 sources in VSA fields selected at 15 GHz with the Ryle Telescope. We have detected 57 sources above a limiting flux density of 5mJy, of which 31 sources have a flux density greater than 10mJy, which is our effective completeness limit. From these measurements we derive a surface density of sources above 10mJy at 30 GHz of 2.0+/-0.4 per square degree. This is consistent with the surface density obtained by Mason et al. (2009) who observed a large sample of sources selected at a much lower frequency (1.4 GHz). We have also investigated the dependence of the spectral index distribution on flux density by comparing our results with those for sources above 1 Jy selected from the WMAP 22 GHz catalogue. We conclude that the proportion of steep spectrum sources increases with decreasing flux density, qualitatively consistent with the predictions of deZotti et al. (2005). We find no evidence for an unexpected population of sources whose spectra rise towards high frequencies, which would affect our ability to interpret current high resolution CMB observations at 30 GHz and above., Comment: 12 pages, 12 figures, submitted to MNRAS

Aims. We investigate the properties of the host galaxy of the blazar J0324+3410 (B2 0321+33) by the analysis of B and R images obtained with the NOT under good photometric conditions. Methods: The galaxy was studied using different methods: Sersic model fitting, unsharp-masked images, B-R image and B-R profile analysis. Results: The images show that the host galaxy has a ring-like morphology. The B-R colour image reveals two bluish zones: one that coincides with the nuclear region, interpreted as the signature of emission related to the active nucleus, the other zone is extended and is located in the host ring-structure. We discuss the hypothesis that the later is thermal emission from a burst of star formation triggered by an interacting/merging process.

We report the results of a survey of 442 planetary nebulae at 30 GHz. The purpose of the survey is to develop a list of planetary nebulae as calibration sources which could be used for high frequency calibration in future. For 41 PNe with sufficient data, we test the emission mechanisms in order to evaluate whether or not spinning dust plays an important role in their spectra at 30 GHz. The 30-GHz data were obtained with a twin-beam differencing radiometer, OCRA-p, which is in operation on the Torun 32-m telescope. Sources were scanned both in right ascension and declination. We estimated flux densities at 30 GHz using a free-free emission model and compared it with our data. The primary result is a catalogue containing the flux densities of 93 planetary nebulae at 30 GHz. Sources with sufficient data were compared with a spectral model of free-free emission. The model shows that free-free emission can generally explain the observed flux densities at 30 GHz thus no other emission mechanism is needed to account for the high frequency spectra., Comment: 10 pages, 7 Postscript figures, to be published in A&A

We investigate the correlations between optical and radio isophotal position angles for 14302 SDSS galaxies with $r$ magnitudes brighter than 18 and which have been associated with extended FIRST radio sources. We identify two separate populations of galaxies using the colour, concentration and their principal components. Surprisingly strong statistical alignments are found: late-type galaxies are overwhelmingly biased towards a position angle differences of $0^{\circ}$ and early-type galaxies to $90^{\circ}$. The late-type alignment can be easily understood in terms of the standard picture in which the radio emission is intimately related to areas of recent star-formation. In early-type galaxies the radio emission is expected to be driven by accretion on to a nuclear black hole. We argue that the observed correlation of the radio axis with the minor axis of the large-scale stellar distribution gives a fundamental insight into the structure of elliptical galaxies, for example, whether or not the nuclear kinematics are decoupled form the rest of the galaxy. Our results imply that the galaxies are oblate spheroids with their radio emission aligned with the minor axis. Remarkably the strength of the correlation of the radio major axis with the optical minor axis depends on radio loudness. Those objects with a low ratio of FIRST radio flux density to total stellar light show a strong minor axis correlation while the stronger radio sources do not. This may reflect different formation histories for the different objects and we suggest we may be seeing the different behaviour of rationally supported and non-rotationally supported ellipticals., Comment: Version to appear in MNRAS

We have examined the statistics of the polarization position angles determined for point sources in the NRAO-VLA sky survey (NVSS) and find that there is a statistically significant bias toward angles which are multiples of 45 degrees. The formal probability that the polarization angles are drawn from a uniform distribution is exponentially small. When the sample of those NVSS sources with polarizations detected with a signal to noise $\geq$3 is split either around the median polarized flux density or the median fractional polarization, the effect appears to be stronger for the more highly polarized sources. Regions containing strong sources and regions at low galactic latitudes are not responsible for the non-uniform distribution of position angles. We identify CLEAN bias as the probable cause of the dominant effect, coupled with small multiplicative and additive offsets on each of the Stokes parameters. Our findings have implications for the extraction of science, such as information concerning galactic magnetic fields, from large scale polarization surveys.

Context: We present the soft X-ray properties obtained in the ROSAT All-Sky survey and from pointed PSPC observations for the AGN in the complete flux-density limited Caltech-Jodrell Bank flat spectrum sample (hereafter CJF). CJF is a VLBI survey (VLBA observations at 5 GHz) of 293 AGN with detailed information on jet component motion. Aims: We investigate and discuss the soft X-ray properties of this AGN sample and examine the correlations between X-ray and VLBI properties, test beaming scenarios, and search for the discriminating properties between the sub-samples detected and not detected by ROSAT. Methods: Comparing the observed and the predicted X-ray fluxes by assuming an Inverse Compton (IC) origin for the observed X-rays, we compute the beaming or Doppler factor and contrast the Doppler factors with other beaming indicators derived from the VLBI observations, such as the value of the expansion velocity, and the observed and intrinsic brightness temperature. In addition, we investigate the large-scale radio structure of the AGN and the difference between the pc- and kpc-scale structure (misalignment) with regard to the X-ray observations.Results: We find a nearly linear relation between X-ray and radio luminosities, and a similar but less stringent behaviour for the relation between optical and X-ray luminosities. The quasars detected by ROSAT have a different $\beta_{\rm app}$-redshift relationship compared to the non-detected ones. ROSAT-detected sources tend to reveal extended large-scale radio structures more often. Conclusions: We conclude that beaming alone cannot explain the observed dichotomy of ROSAT detection or non-detection and assume that the large-scale jet structure plays a decisive role.

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., Comment: Figure 1 and Table 2 are only available in electronic form at the CDS and soon at http://www.mpifr-bonn.mpg.de/staff/sbritzen/cjf.html

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., Comment: recommended for publication in A&A; 27 pages; Table C.1. and Fig. C.1. will be available in electronic form at the CDS and soon at http://www.mpifr-bonn.mpg.de/staff/sbritzen/cjf.html

To measure the 30-GHz flux densities of the 293 sources in the Caltech-Jodrell Bank flat-spectrum (CJF) sample. The measurements are part of an ongoing programme to measure the spectral energy distributions of flat spectrum radio sources and to correlate them with the milliarcsecond structures from VLBI and other measured astrophysical properties. The 30-GHz data were obtained with a twin-beam differencing radiometer system mounted on the Torun 32-m telescope. The system has an angular resolution of 1.2 arcmin. Together with radio spectral data obtained from the literature, the 30-GHz data have enabled us to identify 42 of the CJF sources as Giga-hertz Peaked Spectrum (GPS) sources. Seventeen percent of the sources have rising spectra (alpha > 0) between 5 and 30 GHz., Comment: 9 pages, 2 figures, 2 Tables, submitted to Astronomy & Astrophysics

We report a probable gravitational lens J0316+4328, one of 19 candidate asymmetric double lenses (2 images at a high flux density ratio) from CLASS. Observations with the Very Large Array (VLA), MERLIN and the Very Long Baseline Array (VLBA) imply that J0316+4328 is a lens with high confidence. It has 2 images separated by 0.40", with 6 GHz flux densities of 62 mJy and 3.2 mJy. The flux density ratio of ~19 (constant over the frequency range 6-22 GHz) is the largest for any 2 image gravitational lens. High resolution optical imaging and deeper VLBI maps should confirm the lensing interpretation and provide inputs to detailed lens models. The unique configuration will give strong constraints on the lens galaxy's mass profile., Comment: Accepted to MNRAS Letters. 5 pages, 6 figures, 3 tables

We have used the very large JVAS/CLASS 8.4-GHz surveys of flat-spectrum radio sources to test the hypothesis that there is a systematic alignment of polarization position angle vectors on cosmological scales of the type claimed by Hutsemekers et al (2005). The polarization position angles of 4290 sources with polarized flux density >=1 mJy have been examined. They do not reveal large-scale alignments either as a whole or when split in half into high-redshift (z >= 1.24) and low-redshift sub-samples. Nor do the radio sources which lie in the specific areas covered by Hutsemekers et al (2005). show any significant effect. We have also looked at the position angles of parsec-scale jets derived from VLBI observations and again find no evidence for systematic alignments. Finally, we have investigated the correlation between the polarization position angle and those of the parsec-scale jets. As expected, we find that there is a tendency for the polarization angles to be perpendicular to the jet angles. However, the difference in jet and polarization position angles does not show any systematic trend in different parts of the sky., Comment: 13 pages, 11 figures, 1 table, accepted for publication by MNRAS

We have used the very large JVAS/CLASS 8.4-GHz surveys of flat-spectrum radio sources to obtain a large, uniformly observed and calibrated, sample of radio source polarizations. These are useful for many investigations of the properties of radio sources and the interstellar medium. We discuss comparisons with polarization measurements from this survey and from other large-scale surveys of polarization in flat-spectrum sources., Comment: Accepted by MNRAS. 8 pages, 5 figures. Full version of Table 2 available at http://www.jb.man.ac.uk/~njj/classqu_pol

Central gravitational image detection is very important for the study of the mass distribution of the inner parts ($\sim 100$ pc) of lens galaxies. However, the detection of such images is extremely rare and difficult. We present a 1.7-GHz High Sensitivity Array (HSA) observation of the double-image radio lens system B1030+074. The data are combined with archive VLBA and global-VLBI observations, and careful consideration is given to the effects of noise, {\sc clean}ing and self-calibration. An upper limit is derived for the strength of the central image of 180 $\mu$Jy (90% confidence level), considerably greater than would have been expected on the basis of a simple analysis. This gives a lower limit of $\sim 10^3$ for the ratio of the brightest image to the central image. For cusped models of lens mass distributions, we have made use of this non-detection to constrain the relation between inner power-law slope $\beta$ of the lensing galaxy mass profile, and its break radius $r_b$. For $r_b>130$ pc the power-law slope is required to be close to isothermal ($\beta>1.8$). A flatter inner slope is allowed if a massive black hole is present at the centre of the lensing galaxy, but the effect of the black hole is small unless it is $\sim 10$ times more massive than that implied by the relation between black hole mass and stellar velocity dispersion. By comparing four epochs of VLBI observations, we also detected possible superluminal motion in the jet in the brighter A image. The B jet remains unresolved, as expected from a simple lens model of the system., Comment: Accepted by MNRAS; 13 pages, 11 figures, 2 tables; some figures bitmapped

We present the selection of the Jodrell Bank Flat-spectrum (JBF) radio source sample, which is designed to reduce the uncertainties in the Cosmic Lens All-Sky Survey (CLASS) gravitational lensing statistics arising from the lack of knowledge about the parent population luminosity function. From observations at 4.86 GHz with the Very Large Array, we have selected a sample of 117 flat-spectrum radio sources with flux densities greater than 5 mJy. These sources were selected in a similar manner to the CLASS complete sample and are therefore representative of the parent population at low flux densities. The vast majority (~90 per cent) of the JBF sample are found to be compact on the arcsecond scales probed here and show little evidence of any extended radio jet emission. Using the JBF and CLASS complete samples we find the differential number counts slope of the parent population above and below the CLASS 30 mJy flux density limit to be -2.07+/-0.02 and -1.96+/-0.12, respectively., Comment: 10 pages, 4 figures, accepted for publication in MNRAS

Many lens surveys have hitherto used observations of large samples of background sources to select the small minority which are multiply imaged by lensing galaxies along the line of sight. Recently surveys such as SLACS and OLS have improved the efficiency of surveys by pre-selecting double-redshift systems from SDSS. We explore other ways to improve survey efficiency by optimum use of astrometric and morphological information in existing large-scale optical and radio surveys. The method exploits the small position differences between FIRST radio positions of lensed images and the SDSS lens galaxy positions, together with the marginal resolution of some larger gravitational lens systems by the FIRST beam. We present results of a small pilot study with the VLA and MERLIN, and discuss the desirable criteria for future surveys., Comment: Accepted by MNRAS. 9 pages, 5 figures

We address the question of whether or not the properties of all low-luminosity flat spectrum radio sources, not just the obvious BL Lac objects, are consistent with them being the relativistically beamed counterparts of the low radio luminosity radio galaxies (the Fanaroff-Riley type 1 - FRI). We compare the emission-line luminosities of a well studied sample of core-dominated radio sources with a matched sample of FR I radio galaxies. The emission lines in the core-dominated objects are on average significantly more luminous than those in the comparison sample, inconsistent with the simplest unified models in which there is no orientation dependence of the line emission. However, the results are equally consistent with a scenario in which, for the majority of objects, beaming has minimal effect on the observed core luminosities of a large fraction of the FRI population and that intrinsically stronger cores simply give rise to stronger emission lines. We conclude that FRI unification is much more complex than usually portrayed, and models combining beaming with an intrinsic relationship between core and emission line strengths need to be explored., Comment: 9 pages, accepted for publication in MNRAS

We report the discovery of a new gravitational lens system from the CLASS survey, CLASS B0631+519. VLA, MERLIN and VLBA observations show a doubly-imaged radio core, a doubly-imaged lobe and a second lobe that is probably quadruply-imaged. The maximum image separation is 1.16 arcseconds. The VLBA resolves the most magnified image of the flat-spectrum radio core into a number of sub-components spread across approximately 20 milli-arcseconds. Optical and near-infrared imaging with the ACS and NICMOS cameras on the HST show that there are two galaxies along the line of sight to the lensed source, as was previously discovered by optical spectroscopy. The nearer galaxy at z=0.0896 is a small blue irregular, while the more distant galaxy at z=0.6196 is an elliptical type and appears to contribute most of the lensing effect. The host galaxy of the lensed source is visible in NICMOS imaging as a set of arcs that form an almost complete Einstein ring. Mass modelling using non-parametric techniques can reproduce the ring and indicates that the irregular galaxy has a (localised) effect on the flux density distribution in the Einstein ring at the 5-10% level., Comment: 15 pages, 13 figures, accepted for publication in MNRAS

We present observations of CLASS B2108+213, the widest separation gravitational lens system discovered by the Cosmic Lens All-Sky Survey. Radio imaging using the VLA at 8.46 GHz and MERLIN at 5 GHz shows two compact components separated by 4.56 arcsec with a faint third component in between which we believe is emission from a lensing galaxy. 5-GHz VLBA observations reveal milliarcsecond-scale structure in the two lensed images that is consistent with gravitational lensing. Optical emission from the two lensed images and two lensing galaxies within the Einstein radius is detected in Hubble Space Telescope imaging. Furthermore, an optical gravitational arc, associated with the strongest lensed component, has been detected. Surrounding the system are a number of faint galaxies which may help explain the wide image separation. A plausible mass distribution model for CLASS B2108+213 is also presented., Comment: 9 pages, 8 figures, accepted for publication in MNRAS

We analyse a group of radio sources, a subset of the 200 mJy sample, all of which have core-jet radio structures measured with VLBI and have flat spectra stretching from the radio to the millimetre/sub-millimetre band. Thus the objects have most of the properties expected of blazars. However, they display varied optical properties ranging from ``Seyfert-like'' objects, through BL Lac objects, to ``normal'' elliptical galaxies. We investigate the distribution of synchrotron peak frequencies in their Spectral Energy Distributions (SEDs) and find a broad distribution between 10^12 and 10^16 Hz. Our conclusion is that we should consider virtually all objects in the sample as blazars since much of the diversity in their classification based on traditional optical criteria arises from differences in the frequency at which the non-thermal emission begins to decline. Specifically, an object is only classified as BL Lac when its peak frequency falls in the near IR/optical range. We determine peak frequencies using the same method for objects from other blazar samples. An important result is that our objects do not follow the blazar spectral sequence proposed by Fossati et al. and Ghisellini et al. in which, on average, peak frequencies increase as the radio luminosity decreases. Most of our low radio-luminosity sources have peaks in their SEDs at low frequencies, not at the expected high frequencies. We suggest that at least part of the systematic trend seen by Fossati et al. and Ghisellini et al. results from selection effects., Comment: Accepted for publication in MNRAS

The advent of new observational facilities in the last two decades has allowed the rapid discovery and high-resolution optical imaging of many strong lens systems from galaxy to cluster scales, as well as their spectroscopic follow-up. Radio telescopes have played the dominant role in the systematic detection of dozens of new arcsec-scale lens systems. For the future, we expect nothing less! The next major ground- and space-based facilities, especially the Square Kilometer Array can discover tens of thousands of new lens systems in large sky surveys. For optical imaging and spectroscopic follow-up a strong synergy with planned optical facilities is needed. Here, we discuss the field where strong gravitational lensing is expected to play the dominant role and where SKA can have a major impact: The study of the internal mass structure and evolution of galaxies and clusters to z~1. In addition, studies of more exotic phenomena are contemplated. For example, milli- and microlensing can provide a way to measure the mass-functions of stars and CDM substructure at cosmological distances. All-sky radio monitoring will also rapidly develop the field of time-domain lensing., Comment: 9 pages, 3 figures; to appear in "Science with the Square Kilometer Array," eds. C. Carilli and S. Rawlings, New Astronomy Reviews (Elsevier: Amsterdam)

We present deep optical observations of the gravitational lens system CLASS B0218+357, from which we derive an estimate of the Hubble constant. Extensive radio observations have reduced the degeneracies between Hubble's constant and the mass model in this lens to one involving only the position of the radio-quiet lensing galaxy relative to the lensed images. B0218+357 has an image separation of only 334 mas, so optical observations have previously been unable to resolve the lens galaxy from the bright lensed images. Using the new Advanced Camera for Surveys installed on the Hubble Space Telescope, we have measured the separation between the lens galaxy centre and the brightest image. The position found, and hence our estimate of Hubble's constant, depends on our approach to the spiral arms in B0218+357. If the most prominent arms are left unmasked, we find a value for Hubble's constant of 70+/-5 km/s /Mpc (95% confidence). If the spiral arms are masked out, we find a value of 61+/-7 km/s /Mpc (95% confidence)., Comment: Published in MNRAS. 15 pages, 7 figures. Astro-ph version now consistent with published version. Major changes include an improved noise model and combination of optical galaxy position with VLBI constraints

We address the question of why low-luminosity radio sources with similar flat radio spectra show a range of optical activity. The investigation is based on the spectral energy distributions (SEDs) of objects from the 200 mJy sample. We gathered new data from the VLA at 43 GHz, from SCUBA in the JCMT at 2000, 1350 and 850 micron and from the ISOPHOT instrument on ISO at 170, 90, 60 and 25 micron. There is considerable diversity amongst the SEDs of the objects: there are objects with steep broadband spectra between centimetre and millimetre bands (14% of the sample), there are those with flat broadband spectra over most of the spectral range (48% of the sample) and there are those which show pronounced sub-mm/infrared excesses (27% of the sample). Some objects of the first group have two sided radio morphology indicating that their pc-scale emission is not dominated by beamed jet emission. Amongst the objects that have smooth broadband spectra from the radio to the infrared there are passive elliptical galaxies as well as the expected BL Lacs objects. The most pronounced sub-mm/infrared excesses are shown by the broad emission line objects., Comment: Accepted for publication in MNRAS

We present new observations of the gravitational lens system CLASS B0128+437. HST observations detect a very faint, extended object in I-band with no emission from the lensed images visible; no detection at all is made in V-band. The lens system is detected with much higher signal to noise with UKIRT in K-band, but the resolution is not sufficient to allow the lensed images and the lens galaxy to be separated. A careful astrometric calibration, however, suggests that the peak of the infrared emission corresponds to the two merging images A and B and therefore that the lensed images dominate at infrared wavelengths. The new radio data consist of VLBI radio images at three frequencies, 2.3, 5 and 8.4GHz, made with the VLBA and the 100-m Effelsberg telescope. The lensed source consists of three well-defined sub-components embedded in a more extended jet. Due to the fact that the sub-components have different spectral indices it is possible to determine which part of each image corresponds to the same source sub-component. Our main finding is that one of the images, B, looks very different to the others, there being no obvious division into separate sub-components and the image being apparently both broader and smoother. This is a consequence we believe of scatter-broadening in the ISM of the lensing galaxy. The large number of multiply-imaged source sub-components also provide an abundance of modelling constraints and we have attempted to fit an SIE+external shear model to the data, as well as utilising the novel method of Evans & Witt. It proves difficult in both cases, however, to obtain a satisfactory fit which strongly suggests the presence of sub-structure in the mass distribution of the lensing galaxy, perhaps of the kind that is predicted by CDM theories of structure formation., Comment: Accepted for publication in MNRAS. 13 pages, 8 figures. Full resolution versions of Figures 1, 2, 3, 4 and 5 are obtainable from ftp://ftp.jive.nl/pub/biggs/0128_highres.tar.gz

We present the optical spectra of four newly discovered gravitational lenses from the Cosmic Lens All-Sky Survey (CLASS). These observations were carried out using the Low Resolution Imaging Spectrograph on the W. M. Keck-I Telescope as part of a program to study galaxy-scale gravitational lenses. From our spectra we found the redshift of the background source in CLASS B0128+437 (z_s=3.1240+-0.0042) and the lensing galaxy redshifts in CLASS B0445+123 (z_l=0.5583+-0.0003) and CLASS B0850+054 (z_l=0.5883+-0.0006). Intriguingly, we also discovered that CLASS B0631+519 may have two lensing galaxies (z_l,1=0.0896+-0.0001, z_l,2=0.6196+-0.0004). We also found a single unidentified emission line from the lensing galaxy in CLASS B0128+437 and the lensed source in CLASS B0850+054. We find the lensing galaxies in CLASS B0445+123 and CLASS B0631+519 (l,2) to be early-type galaxies with Einstein Radii of 2.8-3.0 h^{-1} kpc. The deflector in CLASS B0850+054 is a late-type galaxy with an Einstein Radius of 1.6 h^{-1} kpc., Comment: 9 pages, 5 figures, accepted for publication in MNRAS

B0218+357 is one of the most promising systems to determine the Hubble constant from gravitational lenses. Consisting of two bright resolved images plus an Einstein ring, it provides better constraints for the mass model than other systems. The main problem left until now was the poorly determined position of the lensing galaxy. After presenting results from classical lens modelling, we apply our improved version of LensClean which utilizes the Einstein ring for lens modelling purposes. The primary result using isothermal models is a well defined lens position which allows the first reliable measurement of the Hubble constant from this system. The result of H0=(78+-6) km/s/Mpc (2 sigma) is high compared with other lenses but compatible with the HST key project and WMAP results. We furthermore discuss effects of different radial mass profiles. The power-law exponent of the potential is constrained by VLBI data to be beta=1.04+-0.02, very close to isothermal. The effect on H0 is expected to be very small. We also present a composite map (lensed and unlensed) which shows the rich structure of B0218+357 on scales from milli-arcseconds to arcseconds. Finally we use a comparison of observations at different frequencies to investigate the question of possible weakening of one of the images by propagation effects and/or source shifts with frequency. The data clearly favour the model of significant extinction without noticeable source position shifts. The technical details of our variant of the LensClean method are presented in the accompanying Paper I., Comment: To appear in MNRAS together with Paper I (astro-ph/0312262). 18 Pages. Full resolution version available at http://www.astro.physik.uni-potsdam.de/~olaf or on request

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., Comment: To be published in "Radio Astronomy at the Fringe", ASP Conf. Ser. Vol. 300, J. A. Zensus, M. H. Cohen, & E. Ros (eds.). 6 pages, 1 figure, needs rafringe.sty

We present flux-ratio curves of the fold and cusp (i.e. close multiple) images of six JVAS/CLASS gravitational lens systems. The data were obtained over a period of 8.5 months in 2001 with the Multi-Element Radio-Linked Interferometer Network (MERLIN) at 5-GHz with 50 mas resolution, as part of a MERLIN Key-Project. Even though the time delays between the fold and cusp images are small (<~1 day) compared to the time-scale of intrinsic source variability, all six lens systems show evidence that suggests the presence of extrinsic variability. In particular, the cusp images of B2045+265 -- regarded as the strongest case of the violation of the cusp relation (i.e. the sum of the magnifications of the three cusp images add to zero) -- show extrinsic variations in their flux-ratios up to ~40 percent peak-to-peak on time scales of several months. Its low Galactic latitude of b=-10 degree and a line-of-sight toward the Cygnus superbubble region suggest that Galactic scintillation is the most likely cause. The cusp images of B1422+231 at b=+69 degree do not show strong extrinsic variability. Galactic scintillation can therefore cause significant scatter in the cusp and fold relations of some radio lens systems (up to 10 percent rms), even though these relations remain violated when averaged over a <~1 year time baseline., Comment: 17pag, 2 figs, 2 tables; version as accepted for publication in ApJ