This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, demonstrate the EVN's potential to generate new and exciting results that will transform our view of the cosmos. Together with e-MERLIN, the EVN provides a range of baseline lengths that permit unique studies of faint radio sources to be made over a wide range of spatial scales. The science cases are reviewed in six chapters that cover the following broad areas: cosmology, galaxy formation and evolution, innermost regions of active galactic nuclei, explosive phenomena and transients, stars and stellar masers in the Milky Way, celestial reference frames and space applications. The document concludes with identifying the synergies with other radio, as well as multi-band/multi-messenger instruments, and provide the recommendations for future improvements. The appendices briefly describe other radio VLBI arrays, the technological framework for EVN developments, and a selection of spectral lines of astrophysical interest below 100 GHz. The document includes a glossary for non-specialists, and a list of acronyms at the end.

We present new high-resolution and high-sensitivity studies of the jets in the WAT source 3C 465, using deep transverse-resolved radio observations from e-MERLIN, and with complementary observations from the VLA. We derive a lower limit $\beta_{\rm j}$ = ($\nu_{\rm j}$/$c$) $\gtrsim$ 0.5 for the jet speed, and an upper limit $\theta_{\rm j}$ $\lesssim$ 61$^{\circ}$ for the jet angle to the line of sight. The jet spectral index ($\alpha$, defined in the sense $S \propto \nu^{\alpha}$) is fairly constant (<$\alpha_{\rm jet}$> = $-$0.7), and spectral flattening within 4.4 kpc of the core coincides with bright knots and is consistent with the site of X-ray particle acceleration at the base of the radio jet found in previous studies. There is little difference between the spectra of the two hotspot components, plausibly indicating that electron populations of the same properties are injected there. The NW and SE plumes are approximately homologous structures, with variations in mass injection and propagation in external pressure and density gradients in the two regions plausibly accounting for the slightly steeper spectrum in the NW plume, <$\alpha_{\rm NWp}$> = $-$1.43 compared with the SE plume, <$\alpha_{\rm SEp}$> = $-$1.38. Our synchrotron lifetime model supports plausible reacceleration of particles within the plume materials. Overall, our results show that the first-order Fermi process at mildly relativistic and non-relativistic shocks is the most likely acceleration mechanism at play in 3C 465 and distinguish differences between the acceleration at $\beta_{\rm j}$ $>$ 0.5 and $\beta_{\rm j}$ $<$ 0.5. The former case can accelerate electrons to higher Lorentz factors.

We discuss the science drivers for ALMA Band 2 which spans the frequency range from 67 to 90 GHz. The key science in this frequency range are the study of the deuterated molecules in cold, dense, quiescent gas and the study of redshifted emission from galaxies in CO and other species. However, Band 2 has a range of other applications which are also presented. The science enabled by a single receiver system which would combine ALMA Bands 2 and 3 covering the frequency range 67 to 116 GHz, as well as the possible doubling of the IF bandwidth of ALMA to 16 GHz, are also considered.

This is the second paper of a series exploring the multi-component (stars, warm and cold gas and radio jets) properties of a sample of eleven nearby low excitation radio galaxies (LERGs), with the aim of better understanding the AGN fuelling/feedback cycle in these objects. Here we present a study of the molecular gas kinematics of six sample galaxies detected in $^{12}$CO(2-1) with ALMA. In all cases, our modelling suggests that the bulk of the gas in the observed (sub-)kpc CO discs is in ordered rotation. Nevertheless, low-level distortions are ubiquitous, indicating that the molecular gas is not fully relaxed into the host galaxy potential. The majority of the discs, however, are only marginally resolved, preventing us from drawing strong conclusions. NGC 3557 and NGC 3100 are special cases. The features observed in the CO velocity curve of NGC 3557 allow us to estimate a super-massive black hole (SMBH) mass of $(7.10\pm0.02)\times10^{8}$ M$_{\odot}$, in agreement with expectations from the M$_{\rm SMBH}- \sigma_{*}$ relation. The rotation pattern of NGC 3100 shows distortions that appear to be consistent with the presence of both a position angle and inclination warp. Non-negligible radial motions are also found in the plane of the CO disc, likely consistent with streaming motions associated with the spiral pattern found in the inner regions of the disc. The dominant radial motions are likely to be inflows, supporting a scenario in which the cold gas is contributing to the fuelling of the AGN.

One of the main considerations in the ALMA Development Roadmap for the future of operations beyond 2030 is to at least double its on-sky instantaneous bandwidth capabilities. Thanks to the technological innovations of the past two decades, we can now produce wider bandwidth receivers than were foreseen at the time of the original ALMA specifications. In several cases, the band edges set by technology at that time are also no longer relevant. In this memo, we look into the scientific advantages of beginning with Band 2 when implementing such wideband technologies. The Band 2 receiver system will be the last of the original ALMA bands, completing ALMA's coverage of the atmospheric windows from 35-950 GHz, and is not yet covered by any other ALMA receiver. New receiver designs covering and significantly extending the original ALMA Band 2 frequency range (67-90 GHz) can now implement these technologies. We explore the scientific and operational advantages of a receiver covering the full 67-116 GHz atmospheric window. In addition to technological goals, the ALMA Development Roadmap provides 3 new key science drivers for ALMA, to probe: 1) the Origins of Galaxies, 2) the Origins of Chemical Complexity, and 3) the Origins of Planets. In this memo, we describe how the wide RF Band 2 system can help achieve these goals, enabling several high-profile science programmes to be executed uniquely or more effectively than with separate systems, requiring an overall much lower array time and achieving more consistent calibration accuracy: contiguous broad-band spectral surveys, measurements of deuterated line ratios, and more generally fractionation studies, improved continuum measurements (also necessary for reliable line flux measurements), simultaneous broad-band observations of transient phenomena, and improved bandwidth for 3 mm very long baseline interferometry (VLBI).

This is the first paper of a series exploring the multi-frequency properties of a sample of eleven nearby low excitation radio galaxies (LERGs) in the southern sky. We are conducting an extensive study of different galaxy components (stars, warm and cold gas, radio jets) with the aim of improving our understanding of the AGN fuelling/feedback cycle in LERGs. We present ALMA Band 6 $^{12}$CO(2-1) and continuum observations of nine sources. Continuum emission from the radio cores was detected in all objects. Six sources also show mm emission from jets on kpc/sub-kpc scales. The jet structures are very similar at mm and cm wavelengths. We conclude that synchrotron emission associated with the radio jets dominates the continuum spectra up to 230 GHz. The $^{12}$CO(2-1) line was detected in emission in six out of nine objects, with molecular gas masses ranging from $2 \times 10^{7}$ to $2 \times 10^{10}$ M$_{\rm \odot}$. The CO detections show disc-like structures on scales from $\approx$0.2 to $\approx$10 kpc. In one case (NGC 3100) the CO disc presents some asymmetries and is disrupted in the direction of the northern radio jet, indicating a possible jet/disc interaction. In IC 4296, CO is detected in absorption against the radio core, as well as in emission. In four of the six galaxies with CO detections, the gas rotation axes are roughly parallel to the radio jets in projection; the remaining two cases show large misalignments. In those objects where optical imaging is available, dust and CO appear to be co-spatial.

We present a comparison of Doppler-shifted H-alpha line emission observed by the Global Jet Watch from freshly-launched jet ejecta at the nucleus of the Galactic microquasar SS433 with subsequent ALMA imaging at mm-wavelengths of the same jet ejecta. There is a remarkable similarity between the transversely-resolved synchrotron emission and the prediction of the jet trace from optical spectroscopy: this is an a priori prediction not an a posteriori fit, confirming the ballistic nature of the jet propagation. The mm-wavelength of the ALMA polarimetry is sufficiently short that the Faraday rotation is negligible and therefore that the observed E-vector directions are accurately orthogonal to the projected local magnetic field. Close to the nucleus the B-field vectors are perpendicular to the direction of propagation. Further out from the nucleus, the B-field vectors that are coincident with the jet instead become parallel to the ridge line; this occurs at a distance where the jet bolides are expected to expand into one another. X-ray variability has also been observed at this location; this has a natural explanation if shocks from the expanding and colliding bolides cause particle acceleration. In regions distinctly separate from the jet ridge line, the fractional polarisation approaches the theoretical maximum for synchrotron emission.

We use 2.0 Msec of Chandra observations to investigate the cocoon shocks of Cygnus A and some implications for its lobes and jet. Measured shock Mach numbers vary in the range 1.18-1.66 around the cocoon. We estimate a total outburst energy of $\simeq 4.7\times10^{60}\rm\ erg$, with an age of $\simeq 2 \times 10^{7}\rm\ yr$. The average postshock pressure is found to be $8.6 \pm 0.3 \times 10^{-10}\rm\ erg\ cm^{-3}$, which agrees with the average pressure of the thin rim of compressed gas between the radio lobes and shocks, as determined from X-ray spectra. However, average rim pressures are found to be lower in the western lobe than in the eastern lobe by $\simeq 20\%$. Pressure estimates for hotspots A and D from synchrotron self-Compton models imply that each jet exerts a ram pressure $\gtrsim$ 3 times its static pressure, consistent with the positions of the hotspots moving about on the cocoon shock over time. A steady, one-dimensional flow model is used to estimate jet properties, finding mildly relativistic flow speeds within the allowed parameter range. Models in which the jet carries a negligible flux of rest mass are consistent with with the observed properties of the jets and hotspots. This favors the jets being light, implying that the kinetic power and momentum flux are carried primarily by the internal energy of the jet plasma rather than by its rest mass.

We use new and archival Chandra observations of Cygnus A, totalling $\sim$1.9 Ms, to investigate the distribution and temperature structure of gas lying within the projected extent of the cocoon shock and exhibiting a rib-like structure. We confirm that the X-rays are dominated by thermal emission with an average temperature of around 4 keV, and have discovered an asymmetry in the temperature gradient, with the southwestern part of the gas cooler than the rest by up to 2 keV. Pressure estimates suggest that the gas is a coherent structure of single origin located inside the cocoon, with a mass of roughly $2\times10^{10} M_{\odot}$. We conclude that the gas is debris resulting from disintegration of the cool core of the Cygnus A cluster after the passage of the jet during the early stages of the current epoch of activity. The 4 keV gas now lies on the central inside surface of the hotter cocoon rim. The temperature gradient could result from an offset between the centre of the cluster core and the Cygnus A host galaxy at the switch-on of current radio activity.

Einstein's General Theory of Relativity (GR) successfully describes gravity. The most fundamental predictions of GR are black holes (BHs), but in spite of many convincing BH candidates in the Universe, there is no conclusive experimental proof of their existence using astronomical observations in the electromagnetic spectrum. Are BHs real astrophysical objects? Does GR hold in its most extreme limit or are alternatives needed? The prime target to address these fundamental questions is in the center of our own Galaxy, which hosts the closest and best-constrained supermassive BH candidate in the Universe, Sagittarius A* (Sgr A*). Three different types of experiments hold the promise to test GR in a strong-field regime using observations of Sgr A* with new-generation instruments. The first experiment aims to image the relativistic plasma emission which surrounds the event horizon and forms a "shadow" cast against the background, whose predicted size (~50 microarcseconds) can now be resolved by upcoming VLBI experiments at mm-waves such as the Event Horizon Telescope (EHT). The second experiment aims to monitor stars orbiting Sgr A* with the upcoming near-infrared interferometer GRAVITY at the Very Large Telescope (VLT). The third experiment aims to time a radio pulsar in tight orbit about Sgr A* using radio telescopes (including the Atacama Large Millimeter Array or ALMA). The BlackHoleCam project exploits the synergy between these three different techniques and aims to measure the main BH parameters with sufficient precision to provide fundamental tests of GR and probe the spacetime around a BH in any metric theory of gravity. Here, we review our current knowledge of the physical properties of Sgr A* as well as the current status of such experimental efforts towards imaging the event horizon, measuring stellar orbits, and timing pulsars around Sgr A*.

The Chandra X-ray observatory has discovered several dozen anomalously X-ray-bright jets associated with powerful quasars. A popular explanation for the X-ray flux from the knots in these jets is that relativistic synchrotron-emitting electrons inverse-Compton scatter Cosmic Microwave Background (CMB) photons to X-ray energies (the IC/CMB model). This model predicts a high gamma-ray flux which should be detectable by the Fermi Large Area Telescope (LAT) for many sources. GeV-band upper limits from Fermi/LAT for the well-known anomalous X-ray jet in PKS 0637-752 were previously shown in Meyer et al., (2015) to violate the predictions of the IC/CMB model. Previously, measurements of the jet synchrotron spectrum, important for accurately predicting the gamma-ray flux level, were lacking between radio and infrared wavelengths. Here we present new Atacama Large Millimeter/submillimeter Array (ALMA) observations of the large-scale jet at 100, 233, and 319 GHz which further constrain the synchrotron spectrum, supporting the previously published empirical model. We also present updated limits from the Fermi/LAT using the new `Pass 8' calibration and approximately 30% more time on source. With these deeper limits we rule out the IC/CMB model at the 8.7 sigma level. Finally, we demonstrate that complete knowledge of the synchrotron SED is critical in evaluating the IC/CMB model.

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

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

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

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

We explore the use of SKA to deduce the physical parameters of kiloparsec-scale jet flows in radio galaxies. Jets in Active Galactic Nuclei are relativistic where they are first formed, but their speeds and compositions change as they propagate. It has long been known that kiloparsec-scale jets in radio galaxies can be divided into two flavours: strong (found in powerful sources, narrow and terminating in compact hot-spots) and weak (found in low-luminosity sources, geometrically flaring, unable to form hot-spots and terminating in diffuse lobes or tails). We have developed methods to model AGN jets as intrinsically symmetrical, relativistic flows by fitting to deep, well-resolved radio images in Stokes I, Q and U. This has yielded a wealth of information about the brightest few weak-flavour jets. Our first key objective is to observe large samples of weak and transition jets at 0.1 - 0.5 arcsec resolution with SKA1-MID. This would allow us to see how jet propagation depends on power and environment and to quantify the energy and momentum input into the IGM. We will require typical noise levels of 1 microJy/beam, and may be able to exploit survey imaging in some cases. Our second, more challenging, application is to determine the velocity fields in strong-flavour jets. Do they have very fast spines with bulk Lorentz factors of 5 - 10? Is there evidence for magnetic confinement by a toroidal field? What are their energy fluxes? This is a major imaging challenge for SKA2: we need resolution better than 0.05 arcsec, ideally in the 1 - 10 GHz frequency range, with rms noise levels of roughly 10 nJy/beam and extremely high dynamic range, imaging fidelity and polarization purity.

Relativistic jets in active galactic nuclei (AGN) are among the most powerful astrophysical objects discovered to date. Indeed, jetted AGN studies have been considered a prominent science case for SKA, and were included in several different chapters of the previous SKA Science Book (Carilli & Rawlings 2004). Most of the fundamental questions about the physics of relativistic jets still remain unanswered, and await high-sensitivity radio instruments such as SKA to solve them. These questions will be addressed specially through analysis of the massive data sets arising from the deep, all-sky surveys (both total and polarimetric flux) from SKA1. Wide-field very-long-baseline-interferometric survey observations involving SKA1 will serve as a unique tool for distinguishing between extragalactic relativistic jets and star forming galaxies via brightness temperature measurements. Subsequent SKA1 studies of relativistic jets at different resolutions will allow for unprecedented cosmological studies of AGN jets up to the epoch of re-ionization, enabling detailed characterization of the jet composition, magnetic field, particle populations, and plasma properties on all scales. SKA will enable us to study the dependence of jet power and star formation on other properties of the AGN system. SKA1 will enable such studies for large samples of jets, while VLBI observations involving SKA1 will provide the sensitivity for pc-scale imaging, and SKA2 (with its extraordinary sensitivity and dynamic range) will allow us for the first time to resolve and model the weakest radio structures in the most powerful radio-loud AGN.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

Jets in low-luminosity radio galaxies are known to decelerate from relativistic speeds on parsec scales to mildly or sub-relativistic speeds on kiloparsec scales. Several mechanisms have been proposed to explain this effect, including strong reconfinement shocks and the growth of instabilities (both leading to boundary-layer entrainment) and mass loading from stellar winds or molecular clouds. We have performed a series of axisymmetric simulations of the early evolution of jets in a realistic ambient medium to probe the effects of mass loading from stellar winds using the code Ratpenat. We study the evolution of Fanaroff-Riley Class I (FRI) jets, with kinetic powers L_j \sim 1.e41-1.e44 erg/s, within the first 2 kpc of their evolution, where deceleration by stellar mass loading should be most effective. Mass entrainment rates consistent with present models of stellar mass loss in elliptical galaxies produce deceleration and effective decollimation of weak FRI jets within the first kiloparsec. However, powerful FRI jets are not decelerated significantly. In those cases where the mass loading is important, the jets show larger opening angles and decollimate at smaller distances, but the overall structure and dynamics of the bow-shock are similar to those of unloaded jets with the same power and thrust. According to our results, the flaring observed on kpc scales is initiated by mass loading in the weaker FRI jets and by reconfinement shocks or the growth of instabilities in the more powerful jets. The final mechanism of decollimation and deceleration is always the development of disruptive pinching modes.

We model the kiloparsec-scale synchrotron emission from jets in 10 Fanaroff-Riley Class I radio galaxies for which we have sensitive, high-resolution imaging and polarimetry from the Very Large Array. We assume that the jets are intrinsically symmetrical, axisymmetric, decelerating, relativistic outflows and we infer their inclination angles and the spatial variations of their flow velocities, magnetic field structures and emissivities using a common set of fitting functions. The inferred inclinations agree well with independent indicators. The spreading rates increase rapidly, then decrease, in a flaring region. The jets then recollimate to form conical outer regions at distance r0 from the active galactic nucleus (AGN). The flaring regions are homologous when scaled by r0. At ~0.1r0, the jets brighten abruptly at the onset of a high-emissivity region and we find an outflow speed of ~0.8c, with a uniform transverse profile. Jet deceleration first becomes detectable at ~0.2r0 and the outflow often becomes slower at its edges than it is on-axis. Deceleration continues until ~0.6r0, after which the outflow speed is usually constant. The dominant magnetic-field component is longitudinal close to the AGN and toroidal after recollimation, but the field evolution is initially much slower than predicted by flux-freezing. In the flaring region, acceleration of ultrarelativistic particles is required to counterbalance the effects of adiabatic losses and account for observed X-ray synchrotron emission, but the brightness evolution of the outer jets is consistent with adiabatic losses alone. We interpret our results as effects of the interaction between the jets and their surroundings. (Slightly abridged).

An international consortium is presently constructing a beamformer for the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile that will be available as a facility instrument. The beamformer will aggregate the entire collecting area of the array into a single, very large aperture. The extraordinary sensitivity of phased ALMA, combined with the extremely fine angular resolution available on baselines to the Northern Hemisphere, will enable transformational new very long baseline interferometry (VLBI) observations in Bands 6 and 7 (1.3 and 0.8 mm) and provide substantial improvements to existing VLBI arrays in Bands 1 and 3 (7 and 3 mm). The ALMA beamformer will have impact on a variety of scientific topics, including accretion and outflow processes around black holes in active galactic nuclei (AGN), tests of general relativity near black holes, jet launch and collimation from AGN and microquasars, pulsar and magnetar emission processes, the chemical history of the universe and the evolution of fundamental constants across cosmic time, maser science, and astrometry.

We present accurate, spatially resolved imaging of radio spectra at the bases of jets in eleven low-luminosity (Fanaroff-Riley I) radio galaxies, derived from Very Large Array (VLA) observations. We pay careful attention to calibration and to the effects of random and systematic errors, and we base the flux-density scale on recent measurements of VLA primary amplitude calibrators by Perley & Butler (2013). We show images and profiles of spectral index over the frequency range 1.4 - 8.5 GHz, together with values integrated over fiducial regions defined by our relativistic models of the jets. We find that the spectral index decreases (the spectrum flattens) with distance from the nucleus in all of the jets. The mean spectral indices are 0.66 +/- 0.01 where the jets first brighten abruptly and 0.59 +/- 0.01 after they recollimate. The mean change in spectral index between these locations, which is independent of calibration and flux-density scale errors and is therefore more accurately and robustly measured, is -0.067 +/- 0.006. Our jet models associate the decrease in spectral index with a bulk deceleration of the flow from 0.8c to <0.5c. We suggest that the decrease is the result of a change in the characteristics of ongoing particle acceleration. One possible acceleration mechanism is the first-order Fermi process in mildly relativistic shocks: in the Bohm limit, the index of the electron energy spectrum, p, is slightly larger than 2 and decreases with velocity upstream of the shock. This possibility is consistent with our measurements, but requires shocks throughout the jet volume rather than at a few discrete locations. A second possibility is that two acceleration mechanisms operate in these jets: one (with p = 2.32) dominant close to the galactic nucleus and associated with high flow speeds, another (with p = 2.18) taking over at larger distances and slower flow speeds.

We examine the relative orientation of radio jets and dusty tori surrounding the AGN in powerful radio galaxies at z > 1. The radio core dominance R = P(20 GHz) /P(500 MHz) serves as an orientation indicator, measuring the ratio between the anisotropic Doppler-beamed core extended core emission and the isotropic lobe emission. Assuming a fixed cylindrical geometry for the hot, dusty torus, we derive its inclination i by fitting optically-thick radiative transfer models to spectral energy distributions obtained with the Spitzer Space Telescope. We find a highly significant anti-correlation (p < 0.0001) between R and i in our sample of 35 type 2 AGN combined with a sample of 18 z = 1 3CR sources containing both type 1 and 2 AGN. This analysis provides observational evidence both for the Unified scheme of AGN and for the common assumption that radio jets are in general perpendicular to the plane of the torus. The use of inclinations derived from mid-infrared photometry breaks several degeneracies which have been problematic in earlier analyses. We illustrate this by deriving the core Lorentz factor Gamma from the R-i anti-correlation, finding Gamma > 1.3.

We show that asymmetries in total intensity and linear polarization between the radio jets and counter-jets in two lobed Fanaroff-Riley Class I (FR I) radio galaxies, B2 0206+35 (UGC 1651) and B2 0755+37 (NGC 2484), can be accounted for if these jets are intrinsically symmetrical, with decelerating relativistic outflows surrounded by mildly relativistic backflows. Our interpretation is motivated by sensitive, well-resolved Very Large Array imaging which shows that both jets in both sources have a two-component structure transverse to their axes. Close to the jet axis, a centrally-darkened counter-jet lies opposite a centrally-brightened jet, but both are surrounded by broader collimated emission that is brighter on the counter-jet side. We have adapted our previous models of FR I jets as relativistic outflows to include an added component of symmetric backflow. We find that the observed radio emission, after subtracting contributions from the extended lobes, is well described by models in which decelerating outflows with parameters similar to those derived for jets in plumed FR I sources are surrounded by backflows containing predominantly toroidal magnetic fields. These return to within a few kpc of the galaxies with velocities of roughly 0.25c and radiate with a synchrotron spectral index close to 0.55. We discuss whether such backflow is to be expected in lobed FR I sources and suggest ways in which our hypothesis can be tested by further observations.

We explore the magneto-ionic environment of the isolated radio galaxy B2 0755+37 using detailed imaging of the distributions of Faraday rotation and depolarization over the radio source from Very Large Array observations at 1385,1465 and 4860 MHz and new X-ray data from XMM-Newton. The Rotation Measure (RM) distribution is complex, with evidence for anisotropic fluctuations in two regions. The approaching lobe shows low and uniform RM in an unusual `stripe' along an extension of the jet axis and a linear gradient transverse to this axis over its Northern half. The leading edge of the receding lobe shows arc-like RM structures with sign reversals. Elsewhere, the RM structures are reasonably isotropic. The RM power spectra are well described by cut-off power laws with slopes ranging from 2.1 to 3.2 in different sub-regions. The corresponding magnetic-field autocorrelation lengths, where well-determined, range from 0.25 to 1.4 kpc. It is likely that the fluctuations are mostly produced by compressed gas and field around the leading edges of the lobes. We identify areas of high depolarization around the jets and inner lobes. These could be produced by dense gas immediately surrounding the radio emission containing a magnetic field which is tangled on small scales. We also identify four ways in which the well known depolarization (Faraday depth) asymmetry between jetted and counter-jetted lobes of extended radio sources can be modified by interactions with the surrounding medium.

We present deep, high-resolution imaging of the nearby Fanaroff-Riley Class I (FR I) radio galaxies NGC 193, B2 0206+35, B2 0755+37 and M 84 at frequencies of 4.9 and 1.4 GHz using new and archival multi-configuration observations from the Very Large Array. In addition, we describe lower-resolution observations of B2 0326+39 and a reanalysis of our published images of 3C 296. All of these radio galaxies show twin jets and well-defined lobes or bridges of emission, and we examine the common properties of this class of source. We show detailed images of total intensity, brightness gradient, spectral index, degree of polarization and projected magnetic-field direction. The jet bases are very similar to those in tailed twin-jet sources and show the characteristics of decelerating, relativistic flows. Except on one side of M 84, we find that the jets can be traced at least as far as the ends of the lobes, where they often form structures which we call "caps" with sharp outer brightness gradients. Continuing, but less well collimated flows back into the lobes from the caps can often be identified by their relatively flat spectral indices. The lobes in these radio galaxies are similar in morphology, spectral-index distribution and magnetic-field structure to those in more powerful (FR II) sources, but lack hot-spots or other evidence for strong shocks at the ends of the jets. M 84 may be an intermediate case between lobed and tailed sources, in which one jet does not reach the end of its lobe, but disrupts to form a "bubble".

We present detailed imaging of Faraday rotation and depolarization for the radio galaxies 0206+35, 3C 270, 3C 353 and M 84, based on Very Large Array observations at multiple frequencies in the range 1365 to 8440 MHz. This work suggests a more complex picture of the magneto-ionic environments of radio galaxies than was apparent from earlier work. All of the sources show spectacular banded rotation measure (RM) structures with contours of constant RM perpendicular to the major axes of their radio lobes. We give a comprehensive description of the banded RM phenomenon and present an initial attempt to interpret it as a consequence of interactions between the sources and their surroundings. We show that the material responsible for the Faraday rotation is in front of the radio emission and that the bands are likely to be caused by magnetized plasma which has been compressed by the expanding radio lobes. A two-dimensional magnetic structure in which the field lines are a family of ellipses draped around the leading edge of the lobe can produce RM bands in the correct orientation for any source orientation. We also report the first detections of rims of high depolarization at the edges of the inner radio lobes of M 84 and 3C 270. These are spatially coincident with shells of enhanced X-ray surface brightness, in which both the field strength and the thermal gas density are likely to be increased by compression.

Context. It has recently been proposed that the jets of low-luminosity radio galaxies are powered by direct accretion of the hot phase of the IGM onto the central black hole. Cold gas remains a plausible alternative fuel supply, however. The most compelling evidence that cold gas plays a role in fueling radio galaxies is that dust is detected more commonly and/or in larger quantities in (elliptical) radio galaxies compared with radio-quiet elliptical galaxies. On the other hand, only small numbers of radio galaxies have yet been detected in CO (and even fewer imaged), and whether or not all radio galaxies have enough cold gas to fuel their jets remains an open question. If so, then the dynamics of the cold gas in the nuclei of radio galaxies may provide important clues to the fuelling mechanism. Aims. The only instrument capable of imaging the molecular component on scales relevant to the accretion process is ALMA, but very little is yet known about CO in southern radio galaxies. Our aim is to measure the CO content in a complete volume-limited sample of southern radio galaxies, in order to create a well-defined list of nearby targets to be imaged in the near future with ALMA. Methods. APEX has recently been equipped with a receiver (APEX-1) able to observe the 230 GHz waveband. This allows us to search for CO(2-1) line emission in our target galaxies. Results. Here we present the results for our first three southern targets, proposed for APEX-1 spectroscopy during science verification: NGC3557, IC4296 and NGC1399. The experiment was successful with two targets detected, and possible indications for a double-horned CO line profile, consistent with ordered rotation. These early results are encouraging, demonstrating that APEX can efficiently detect CO in nearby radio galaxies. We therefore plan to use APEX to obtain CO spectroscopy for all our southern targets.

We describe a survey in the ELAIS N2 region with the VLA at 43.4 GHz, carried out with 1627 independent snapshot observations in D-configuration and covering about 0.5 square degrees. One certain source is detected, a previously-catalogued flat-spectrum QSO at z=2.2. A few (<5) other sources may be present at about the 3sigma level, as determined from positions of source-like deflections coinciding with blue stellar objects, or with sources from lower-frequency surveys. Independently we show how all the source-like detections identified in the data can be used with a maximum-likelihood technique to constrain the 43-GHz source counts at a level of ~7 mJy. Previous estimates of the counts at 43 GHz, based on lower-frequency counts and spectral measurements, are consistent with these constraints, although the present results are suggestive of somewhat higher surface densities at the 7 mJy level. They do not provide direct evidence of intrusion of a previously unknown source population, although the several candidate sources need examination before such a population can be ruled out.

The goal of this work is to constrain the strength and structure of the magnetic field associated with the environment of the radio source 3C 449, using observations of Faraday rotation, which we model with a structure function technique and by comparison with numerical simulations. We assume that the magnetic field is a Gaussian, isotropic random variable and that it is embedded in the hot intra-group plasma surrounding the radio source. For this purpose, we present detailed rotation measure images for the polarized radio source 3C 449, previously observed with the Very Large Array at seven frequencies between 1.365 and 8.385 GHz. We quantify the statistics of the magnetic-field fluctuations by deriving rotation measure structure functions, which we fit using models derived from theoretical power spectra. We quantify the errors due to sampling by making multiple two-dimensional realizations of the best-fitting power spectrum.We also use depolarization measurements to estimate the minimum scale of the field variations. We then make three-dimensional models with a gas density distribution derived from X-ray observations and a random magnetic field with this power spectrum. Under these assumptions we find that both rotation measure and depolarization data are consistent with a broken power-law magnetic-field power spectrum, with a break at about 11 kpc and slopes of 2.98 and 2.07 at smaller and larger scales respectively. The maximum and minimum scales of the fluctuations are around 65 and 0.2 kpc, respectively. The average magnetic field strength at the cluster centre is 3.5 +/-1.2 micro-G, decreasing linearly with the gas density within about 16 kpc of the nucleus.

We present an analytical model for jets in Fanaroff & Riley Class I (FRI) radio galaxies, in which an initially laminar, relativistic flow is surrounded by a shear layer. We apply the appropriate conservation laws to constrain the jet parameters, starting the model where the radio emission is observed to brighten abruptly. We assume that the laminar flow fills the jet there and that pressure balance with the surroundings is maintained from that point outwards. Entrainment continuously injects new material into the jet and forms a shear layer, which contains material from both the environment and the laminar core. The shear layer expands rapidly with distance until finally the core disappears, and all of the material is mixed into the shear layer. Beyond this point, the shear layer expands in a cone and decelerates smoothly. We apply our model to the well-observed FRI source 3C31 and show that there is a self-consistent solution. We derive the jet power, together with the variations of mass flux and and entrainment rate with distance from the nucleus. The predicted variation of bulk velocity with distance in the outer parts of the jets is in good agreement with model fits to VLA observations. Our prediction for the shape of the laminar core can be tested with higher-resolution imaging.

We use high-quality VLA images of the Fanaroff & Riley Class I radio galaxy 3C 31 at six frequencies in the range 1365 to 8440MHz to explore the spatial scale and origin of the rotation measure (RM) fluctuations on the line of sight to the radio source. We analyse the distribution of the degree of polarization to show that the large depolarization asymmetry between the North and South sides of the source seen in earlier work largely disappears as the resolution is increased. We show that the depolarization seen at low resolution results primarily from unresolved gradients in a Faraday screen in front of the synchrotron-emitting plasma. We establish that the residual degree of polarization in the short-wavelength limit should follow a Burn law and we fit such a law to our data to estimate the residual depolarization at high resolution. We show that the observed RM variations over selected areas of 3C 31 are consistent with a power spectrum of magnetic fluctuations in front of 3C 31 whose power-law slope changes significantly on the scales sampled by our data. The power spectrum can only have the form expected for Kolmogorov turbulence on scales <5 kpc. On larger scales we find a flatter slope. We also compare the global variations of RM across 3C 31 with the results of three-dimensional simulations of the magnetic-field fluctuations in the surrounding magnetoionic medium. We show that our data are consistent with a field distribution that favours the plane perpendicular to the jet axis - probably because the radio source has evacuated a large cavity in the surrounding medium. We also apply our analysis techniques to the case of Hydra A, where the shape and the size of the cavities produced by the source in the surrounding medium are known from X-ray data. (Abridged)

We present high-quality VLA images of the FR I radio galaxy 3C 31 in the frequency range 1365 to 8440 MHz with angular resolutions from 0.25 to 40 arcsec. Our new images reveal complex, well resolved filamentary substructure in the radio jets and tails. We also use these images to explore the spectral structure of 3C 31 on large and small scales. We infer the apparent magnetic field structure by correcting for Faraday rotation. Some of the intensity substructure in the jets is clearly related to structure in their apparent magnetic field: there are arcs of emission where the degree of linear polarization increases, with the apparent magnetic field parallel to the ridges of the arcs. The spectral indices are significantly steeper (0.62) within 7 arcsec of the nucleus than between 7 and 50 arcsec (0.52 - 0.57). The spectra of the jet edges are also slightly flatter than the average for their surroundings. At larger distances, the jets are clearly delimited from surrounding larger-scale emission both by their flatter radio spectra and by sharp brightness gradients. The spectral index of 0.62 in the first 7 arcsec of 3C 31's jets is very close to that found in other FR I galaxies where their jets first brighten in the radio and where X-ray synchrotron emission is most prominent. Farther from the nucleus, where the spectra flatten, X-ray emission is fainter relative to the radio. The brightest X-ray emission from FR I jets is therefore not associated with the flattest radio spectra, but with a particle-acceleration process whose characteristic energy index is 2.24. The spectral flattening with distance from the nucleus occurs where our relativistic jet models require deceleration, and the flatter-spectra at the jet edges may be associated with transverse velocity shear. (Slightly abridged)

We present a detailed study of the environments of a sample of nine low-power radio galaxies, based on new and archival XMM-Newton observations. We report new detections of group-scale environments around three radio galaxies, 3C 296, NGC 1044 and 3C 76.1. As with previous studies, we find that FR-I radio galaxies inhabit group environments ranging over nearly two orders of magnitude in bolometric X-ray luminosity, but find no evidence for a tight relationship between large-scale X-ray environment and radio-source properties such as size, radio luminosity, and axial ratio. We confirm earlier work showing that equipartition internal pressures are typically lower than the external pressures acting on the radio lobes, so that additional non-radiating particles must be present or the lobes must be magnetically dominated. We present the first direct observational evidence that entrainment may provide this missing pressure, in the form of a relationship between radio-source structure and apparent pressure imbalance. Finally, our study provides further support for the presence of an apparent temperature excess in radio-loud groups compared to the group population as a whole. For five of eight temperature excesses, the energy required to inflate the radio lobes is comparable to the energy required to produce this excess by heating of the group gas; however, in three cases the current radio source appears too weak to produce the temperature excess. It remains unclear whether the temperature excess can be directly associated with the current phase of AGN activity, or whether it is instead either a signature of previous AGN activity or simply an indicator of the particular set of group properties most conducive to the growth of an FR-I radio galaxy.

There is now unequivocal evidence that the jets in FR I radio galaxies are initially relativistic, decelerating flows. On the assumption that they are axisymmetric and intrinsically symmetrical (a good approximation close to the nucleus), we can make models of their geometry, velocity, emissivity and field structure whose parameters can be determined by fitting to deep VLA observations. Mass entrainment - either from stellar mass loss within the jet volume or via a boundary layer at the jet surface - is the most likely cause for deceleration. This idea is quantitatively consistent with the velocity field and geometry inferred from kinematic modelling and the external gas density and pressure profiles derived from X-ray observations. The jets must initially be very light, perhaps with an electron-positron composition.

We describe very accurate imaging of radio spectral index for the inner jets in three FR I radio galaxies. Where the jets first brighten, there is a remarkably small dispersion around a spectral index of 0.62. This is also the region where bright X-ray emission is detected. Further from the nucleus, the spectral index flattens slightly to 0.50 - 0.55 and X-ray emission, although still detectable, is fainter relative to the radio. The brightest X-ray emission from the jets is therefore not associated with the flattest radio spectra, but rather with some particle-acceleration process whose characteristic energy index is 2.24. The change in spectral index occurs roughly where our relativistic jet models require rapid deceleration. Flatter-spectrum edges can be seen where the jets are isolated from significant surrounding diffuse emission and we suggest that these are associated with shear.

We present an analysis of the magnetic-field fluctuations in the magnetoionic medium in front of the radio galaxy 3C 31 derived from rotation-measure (RM) fits to high-resolution polarization images. We first show that the Faraday rotation must be due primarily to a foreground medium. We determine the RM structure functions for different parts of the source and infer that the simplest form for the power spectrum is a power law with a high-frequency cutoff. We also present three-dimensional simulations of RM produced by a tangled magnetic field in the hot plasma surrounding 3C 31, and show that the observed RM distribution is consistent with a spherical plasma distribution in which the radio source has produced a cavity.

A 10-arcmin field around the HDF(N) contains 92 radio sources >40 uJy, resolved by MERLIN+VLA at 0".2-2".0 resolution. 55 have Chandra X-ray counterparts including 18 with a hard X-ray photon index and high luminosity characteristic of a type-II (obscured) AGN. >70% of the radio sources have been classified as starbursts or AGN using radio morphologies, spectral indices and comparisons with optical appearance and MIR emission. Starbursts outnumber radio AGN 3:1. This study extends the VO methods previously used to identify X-ray-selected obscured type-II AGN to investigate whether very luminous radio and X-ray emission originates from different phenomena in the same galaxy. The high-redshift starbursts have typical sizes of 5--10 kpc and star formation rates of ~1000 Msun/yr. There is no correlation between radio and X-ray luminosities nor spectral indices at z>~1.3. ~70% of both the radio-selected AGN and the starburst samples were detected by Chandra. The X-ray luminosity indicates the presence of an AGN in at least half of the 45 cross-matched radio starbursts, of which 11 are type-II AGN including 7 at z>1.5. This distribution overlaps closely with the X-ray detected radio sources which were also detected by SCUBA. Stacked 1.4-GHz emission at the positions of radio-faint X-ray sources is correlated with X-ray hardness. Most extended radio starbursts at z>1.3 host X-ray selected obscured AGN. Radio emission from most of these ultra-luminous objects is dominated by star formation but it contributes less than 1/3 of their X-ray luminosity. Our results support the inferences from SCUBA and IR data, that at z>1.5, star formation is an order of magnitude more extended and more copious, it is closely linked to AGN activity and it is triggered differently, compared with star formation at lower redshifts.

We present Chandra X-ray results for the jet, nucleus, and gaseous atmosphere of NGC 315, a nearby radio galaxy whose jet kinematics are known through deep radio mapping. Diffuse X-ray synchrotron emission is detected from the jet out to 30 arcsec from the nucleus, through regions both of fast bulk flow and deceleration. The X-ray to radio flux ratio drops considerably where the flow decelerates, but the X-ray and radio emissions show similar transverse extents throughout, requiring distributed particle acceleration to maintain the supply of X-ray-emitting electrons. A remarkable knotty filament within the jet is seen in both the radio and X-ray, contributing roughly 10 per cent of the diffuse emission along its extent at both wavelengths. No completely satisfactory explanation for the filament is found, though its oscillatory appearance, roughly aligned magnetic field, and requirements for particle acceleration, suggest that it is a magnetic strand within a shear layer between fast inner and slower outer flow.

We discuss CO spectral line data of a volume-limited sample of 23 nearby (z<0.03) low luminosity radio galaxies, selected from the B2 catalogue. We investigate whether the CO properties of our sample are correlated with the properties of the host galaxy, and in particular with the dust component. We find strong evidences for a physical link between the dust disks probed by HST in the galaxy cores and the molecular gas probed by the CO spectral lines, which in two cases display a double-horn shape, consistent with ordered rotation. On the other hand, from a preliminary comparison with other samples of radio sources we find no significant differences in molecular gas properties between FRI and FRII radio sources. In order to confirm the suggestion that the CO is dynamically associated with the core dust disks, the most suitable sources of our sample will be proposed for interferometric imaging at PdBI.

We discuss CO spectral line data of a volume-limited sample of 23 nearby (z<0.03) low luminosity radi galaxies, selected from the B2 catalogue. Most of such objects (16/23) have HST imaging. Our aim is to establish the distribution of molecular gas masses in low luminosity radio galaxies, in comparison with other radio source samples, confirm the suggestion that the CO is in ordered rotation, determine its relation to the dust disks observed in these objects with HST and establish targets for future interferometric imaging.

We present new, deep 8.5-GHz VLA observations of the nearby, low-luminosity radio galaxy 3C 296 at resolutions from 0.25 to 5.5 arcsec. These show the intensity and polarization structures of the twin radio jets in detail. We also derive the spectral-index and rotation-measure distributions. We apply our intrinsically symmetrical, decelerating relativistic jet model to the new observations. An optimized model accurately fits the data in both total intensity and linear polarization. We infer that the jets are inclined by 58 deg to the line of sight. On-axis, they decelerate from a (poorly-constrained) initial velocity of 0.8c to approximately 0.4c around 5 kpc from the nucleus, the velocity thereafter remaining constant. The speed at the edge of the jet is low everywhere. The longitudinal profile of proper emissivity has three principal power-law sections and the emission is centre-brightened. Our observations rule out a globally-ordered, helical magnetic-field configuration. Instead, we model the field as random on small scales but anisotropic, with toroidal and longitudinal components only. The ratio of longitudinal to toroidal field falls with distance along the jet, qualitatively but not quantitatively as expected from flux freezing, so that the field is predominantly toroidal far from the nucleus. The toroidal component is relatively stronger at the edges of the jet. A simple adiabatic model fits the emissivity evolution only in the outer region after the jets have decelerated and recollimated; closer to the nucleus, it predicts far too steep an emissivity decline with distance. We also interpret the morphological differences between brightness enhancements ("arcs") in the main and counter-jets as an effect of relativistic aberration. (Abridged)

Our modelling of FR I radio jets as decelerating, relativistic flows allows us to derive their orientations accurately. We present images of Faraday rotation for two of these these objects (3C 31 and NGC 315) and show that the fluctuations of rotation measure (RM) are larger in the fainter (receding) jets, as expected if the rotation occurs in the hot galaxy/group halos. The gas density is much lower in NGC 315 and the RM fluctuations are only just detectable.

We present images of the jets in the nearby radio galaxy NGC 315 made with the VLA at five frequencies between 1.365 and 5 GHz with resolutions between 1.5 and 45 arcsec FWHM. Within 15 arcsec of the nucleus, the spectral index of the jets is 0.61. Further from the nucleus, the spectrum is flatter, with significant transverse structure. Between 15 and 70 arcsec from the nucleus, the spectral index varies from 0.55 on-axis to 0.44 at the edge. This spectral structure suggests a change of dominant particle acceleration mechanism with distance from the nucleus and the transverse gradient may be associated with shear in the jet velocity field. Further from the nucleus, the spectral index has a constant value of 0.47. We derive the distribution of Faraday rotation over the inner +/-400 arcsec of the radio source and show that it has three components: a constant term, a linear gradient (both probably due to our Galaxy) and residual fluctuations at the level of 1 - 2 rad/m^2. These residual fluctuations are smaller in the brighter (approaching) jet, consistent with the idea that they are produced by magnetic fields in a halo of hot plasma that surrounds the radio source. We model this halo, deriving a core radius of approximately 225 arcsec and constraining its central density and magnetic-field strength. We also image the apparent magnetic-field structure over the first +/-200 arcsec from the nucleus.

We report the discovery of a new optical-IR synchrotron jet in the radio galaxy 3C133 from our HST/NICMOS snapshot survey. The jet and eastern hotspot are well resolved, and visible at both optical and IR wavelengths. The IR jet follows the morphology of the inner part of the radio jet, with three distinct knots identified with features in the radio. The radio-IR SED's of the knots are examined, along with those of two more distant hotspots at the eastern extreme of the radio feature. The detected emission appears to be synchrotron, with peaks in the NIR for all except one case, which exhibits a power-law spectrum throughout.

The question of the degree of order in the magnetic fields of relativistic jets is important to any understanding of their production. Both vector-ordered (e.g. helical) and disordered, but anisotropic fields can produce the high observed degrees of polarization. We outline our models of jets in FR I radio galaxies as decelerating relativistic flows. We then present theoretical calculations of the synchrotron emission from different field configurations and compare them with observed emission from FR I jets. We show that large-scale helical fields (with significant poloidal and toroidal components) are inconsistent with observations. The combination of an ordered toroidal and disordered poloidal component is consistent with our data, as is an entirely disordered field. Jets must also contain small, but significant amounts of radial field.

We apply our intrinsically symmetrical, decelerating relativistic jet model to deep VLA imaging of the inner 140 arcsec of the giant low-luminosity radio galaxy NGC 315. An optimized model accurately fits the data in both total intensity and linear polarization. We infer that the velocity, emissivity and field structure in NGC 315 are very similar to those of the other low-luminosity sources we have modelled, but that all of the physical scales are larger by a factor of about 5. We derive an inclination to the line of sight of 38 degrees for the jets. Where they first brighten, their on-axis velocity is approximately v/c = 0.9. They decelerate to v/c = 0.4 between 8 and 18 kpc from the nucleus and the velocity thereafter remains constant. The speed at the edge of the jet is roughly 0.6 of the on-axis value where it is best constrained, but the transverse velocity profile may deviate systematically from the Gaussian form we assume. The proper emissivity profile is split into three power-law regions separated by shorter transition zones. In the first of these, at 3 kpc (the flaring point) the jets expand rapidly at constant emissivity, leading to a large increase in the observed brightness on the approaching side. At 10 kpc, the emissivity drops abruptly by a factor of 2. Where the jets are well resolved their rest-frame emission is centre-brightened. The magnetic field is modelled as random on small scales but anisotropic and we rule out a globally ordered helical configuration. To a first approximation, the field evolves from a mixture of longitudinal and toroidal components to predominantly toroidal, but it also shows variations in structure along and across the jets, with a significant radial component in places. Simple adiabatic models fail to fit the emissivity variations.

Modelling of deep VLA images of the jets in FR I radio galaxies has allowed us to derive their three-dimensional distributions of velocity, emissivity and magnetic-field structure on kiloparsec scales. By combining our models of jet kinematics with measurements of the external pressure and density derived from Chandra observations, we can also determine the jet dynamics via a conservation-law approach. The result is a detailed and quantitative picture of jet deceleration. We discuss the potential application of these techniques on VLBI scales. Our fundamental assumption is that the jets are intrinsically symmetrical, axisymmetric, relativistic flows. Although this is likely to be a good approximation on average, the effects of non-stationarity in the flow (e.g. shocks/knots) may limit its applicability on pc scales. We also stress the need for both the main and counter-jet to be detected with good transverse resolution in linear polarization as well as total intensity. Two foreground effects (free-free absorption and Faraday rotation) must also be corrected. We comment on the implications of observed VLBI polarization and Faraday rotation for the field structure (and, by implication, collimation) of pc-scale jets. Where our VLA observations start to resolve the jets, about 1 kpc from the nucleus, they are travelling at about 0.8c. We briefly discuss the deceleration of jets from a much faster initial speed.

We have observed the twin-jet radio galaxy 3C296 with Chandra. X-ray emission is detected from the nucleus, from the inner parts of the radio jet, and from a small-scale thermal environment around the jet deceleration region. As we have found in previous observations of other twin-jet radio galaxies, the X-ray jet and a steep pressure gradient in the external thermal environment are associated with the region where strong bulk deceleration of the jet material is suggested by radio observations. Our observations provide additional evidence that the inner jets of twin-jet objects are always associated with a relatively cool, dense central X-ray emitting component with a short cooling time.

We summarise the recent progress of a project to determine physical conditions in the jets of low-luminosity (FR I) radio galaxies. We model the jets as intrinsically symmetrical, relativistic, decelerating flows. By fitting to deep, high-resolution VLA images in total intensity and linear polarization, we have been able to derive for the first time the three-dimensional distributions of velocity, rest-frame emissivity and magnetic-field structure. We describe our results for four sources and outline the implications for theories of jet deceleration. We also summarise a conservation-law analysis of 3C 31, which gives profiles of internal pressure, density, Mach number and entrainment rate along the jets, again for the first time. Finally, we discuss the applicability of adiabatic models to FR I jets and identify the regions where distributed particle acceleration is necessary

We apply the intrinsically symmetrical, decelerating relativistic jet model developed by Laing & Bridle for 3C 31 to deep, full-synthesis 8.4-GHz VLA imaging of the two low-luminosity radio galaxies B2 0326+39 and B2 1553+24. After some modifications to the functional forms used to describe the geometry, velocity, emissivity and magnetic-field structure, these models can accurately fit our data in both total intensity and linear polarization. We conclude that the jets in B2 0326+39 and B2 1553+24 are at angles of 64 +/- 5 deg and 7.7 +/- 1.3 deg to the line of sight, respectively. In both objects, we find that the jets decelerate from 0.7 - 0.8c to <0.2c over a distance of approximately 10 kpc, although in B2 1553+24 this transition occurs much further from the nucleus than in B2 0326+39 or 3C 31. The longitudinal emissivity profiles can be divided into sections, each fit accurately by a power law; the indices of these power laws decrease with distance from the nucleus. The magnetic fields in both objects are dominated by the longitudinal component in the high-velocity regions close to the nucleus and by the toroidal component further out, but B2 0326+39 also has a significant radial component at large distances, whereas B2 1553+24 does not. Simple adiabatic models fail to fit the emissivity variations in the regions of high velocity but provide good descriptions of the emissivity after the jets have decelerated. Given the small angle to the line of sight inferred for B2 1553+24, there should be a significant population of similar sources at less extreme orientations. Such objects should have long (>200 kpc), straight, faint jets and we show that their true sizes are likely to have been underestimated in existing images. (Slightly abridged.)

We present a general approach to the modelling of the brightness and polarization structures of adiabatic, decelerating relativistic jets, based on the formalism of Matthews & Scheuer (1990). We compare the predictions of adiabatic jet models with deep, high-resolution observations of the radio jets in the FR I radio galaxy 3C 31. Adiabatic models require coupling between the variations of velocity, magnetic field and particle density. They are therefore more tightly constrained than the models previously presented for 3C 31 by Laing & Bridle (2002). We show that adiabatic models provide a poorer description of the data in two crucial respects: they cannot reproduce the observed magnetic-field structures in detail, and they also predict too steep a brightness decline along the jets for plausible variations of the jet velocity. We find that the innermost regions of the jets show the strongest evidence for non-adiabatic behaviour, and that the adiabatic models provide progressively better descriptions of the jet emission at larger distances from the galactic nucleus. We briefly discuss physical processes which might contribute to this non-adiabatic behaviour. In particular, we develop a parameterized description of distributed particle injection, which we fit to the observed total intensities. We show that particles are preferentially injected where bright X-ray emission is observed, and where we infer that the jets are over-pressured.

We model jets in low-luminosity (FR I) radio galaxies as intrinsically symmetrical, axisymmetric, decelerating relativistic flows with transverse velocity gradients. This allows us to derive velocity fields and the three-dimensional distributions of magnetic-field ordering and rest-frame emissivity. A conservation-law analysis, combining the kinematic model with X-ray observations of the surrounding IGM, gives the profiles of internal density, pressure, Mach number and entrainment rate along the jets. We summarize our recently-published results on 3C 31 and outline new work on other sources and adiabatic jet models.

We present HST observations of previously undetected optical jets in the low-luminosity radio galaxies B2 0755+37 and B2 1553+24. We show that there is accurate spatial coincidence between optical and radio emission, implying that the former is likely to be synchrotron radiation. The physical properties of the jets are similar to those known previously: their radio-optical spectral indices are ~0.7 and in B2 0755+37, the spectrum steepens between optical and X-ray wavelengths. Our results support the hypothesis that optical emission is detectable from jets orientated within ~20 degrees of the line of sight for the B2 sample.

During the course of an investigation of the interaction of the radio galaxy M84 and its ambient cluster gas, we found excess X-ray emission aligned with the northern radio jet. The emission extends from the X-ray core of the host galaxy as a weak bridge and then brightens to a local peak coincident with the first detectable radio knot at approx. 2.5" from the core. The second radio knot at 3.3" is brighter in both radio and X-rays. The X-ray jet terminates 3.9" from the core. Although all the evidence suggests that Doppler favoritism augments the emission of the northern jet, it is unlikely that the excess X-ray emission is produced by inverse Compton emission. We find many similarities between the M84 X-ray jet and recent jet detections from Chandra data of low luminosity radio galaxies. For most of these current detections synchrotron emission is the favored explanation for the observed X-rays.

We present a dynamical analysis of the flow in the jets of the low-luminosity radio galaxy 3C 31 based on our earlier geometrical and kinematic model and on estimates of the external pressure and density distributions from Chandra observations. We apply conservation of particles, energy and momentum to derive the variations of pressure and density along the jets and show that there are self-consistent solutions for deceleration by injection of thermal matter. We initially take the jets to be in pressure equilibrium with the external medium at large distances from the nucleus and the momentum flux to be energy flux/c; we then progressively relax these constraints. We infer that the jets are over-pressured compared with the external medium where they start to expand rapidly. Local minima in the density and pressure and maxima in the mass injection rate and Mach number occur at approximately 3 kpc. Further out, the jets decelerate smoothly with a Mach number close to 1. The mass injection rate we infer is comparable with that expected from stellar mass loss at 1 - 2 kpc, but significantly exceeds it at large distances. We conclude that entrainment from the galactic atmosphere is the dominant mass input process far from the nucleus, but that stellar mass loss may also contribute closer in. The occurrence of a significant over-pressure where expansion begins leads us to suggest that there is a stationary shock system, perhaps caused by reconfinement of an initially free jet. Our results are compatible with a jet consisting of electron-positron plasma on parsec scales which picks up thermal matter from stellar mass loss to reach the inferred density and mass flux at 1 kpc, but we cannot rule out an electron-proton composition with a low-energy cut-off. [Abridged]

We show that the principal differences in appearance of the main and counter-jets within 30 arcsec of the nucleus of the low-luminosity radio galaxy 3C 31 can result entirely from the effects of relativistic aberration in two symmetrical, antiparallel, axisymmetric, time-stationary relativistic flows. We develop empirical parameterized models of the jet geometry and the three-dimensional distributions of the velocity, emissivity and magnetic-field structure and optimize their parameters by fitting the predicted synchrotron intensity and polarization to deep 8.4-GHz VLA observations. We conclude that the jets are at roughly 52 degrees to the line of sight, that they decelerate and that they have transverse velocity gradients. They have three regions with distinct kinematics: a narrow inner region; a flaring region of rapid expansion followed by recollimation and a conical outer region. There is a discontinuity in the flow between inner and flaring regions. The on-axis velocity is close to 0.8c until the end of the flaring region, where it drops abruptly to 0.55c, thereafter falling more slowly to 0.25c at the end of the modelled region. Throughout the flaring and outer regions, the velocity at the edge of the jet is 0.7 of its on-axis value. The magnetic field has primarily toroidal and longitudinal components except in the flaring region, where there is a significant radial component at the jet edge. Simple adiabatic models fail in the inner and flaring regions. The inferred transverse velocity profiles and field structure in the flaring region support the idea that the jets decelerate by entraining the external medium. We demonstrate the appearance of our model at other angles to the line of sight and argue that other low-luminosity radio galaxies resemble 3C 31 seen at different orientations. [Abridged]

We have used a deep Chandra observation of the central regions of the twin-jet FRI radio galaxy 3C31 to resolve the thermal X-ray emission in the central few kiloparsecs of the host galaxy, NGC 383, where the jets are thought to be decelerating rapidly. This allows us to make high-precision measurements of the density, temperature and pressure distributions in this region, and to show that the X-ray emitting gas in the centre of the galaxy has a cooling time of only 5 times 10^7 years. In a companion paper these measurements are used to place constraints on models of the jet dynamics. A previously unknown one-sided X-ray jet in 3C31, extending up to 8 arcsec from the nucleus, is detected and resolved. Its structure and steep X-ray spectrum are similar to those of X-ray jets known in other FRI sources, and we attribute the radiation to synchrotron emission from a high-energy population of electrons. In situ particle acceleration is required in the region of the jet where bulk deceleration is taking place. We also present X-ray spectra and luminosities of the galaxies in the Arp 331 chain of which NGC 383 is a member. The spectrum and spatial properties of the nearby bright X-ray source 1E 0104+3153 are used to argue that the soft X-ray emission is mostly due to a foreground group of galaxies rather than to the background broad absorption-line quasar.

Winged, or X-shaped, radio sources form a small class of morphologically peculiar extragalactic sources. We present multi-frequency radio observations of two such sources. We derive maximum ages since any re-injection of fresh particles of 34 and 17Myr for the wings of 3C223.1 and 3C403, respectively, based on the lack of synchrotron and inverse Compton losses. On morphological grounds we favour an explanation in terms of a fast realignment of the jet axis which occurred within a few Myr. There is no evidence for merger activity, and the host galaxies are found to reside in no more than poor cluster environments. A number of puzzling questions remain about those sources: in particular, although the black hole could realign on sufficiently short timescales, the origin of the realignment is unknown.

We derive expressions for the total and linearly-polarized synchrotron emissivity of an element of plasma containing relativistic particles and disordered magnetic field which has been sheared or compressed along three independent directions. Our treatment follows that given by Matthews & Scheuer (1990) in the special case of a power-law electron energy spectrum. We show that the emissivity integrals depend on a single parameter, making it straightforward to generate one-dimensional look-up tables. We also demonstrate that our formulae give identical results to those in the literature in special cases.

We investigate the correlations between spectral index, jet side and extent of the radio lobes for a sample of nearby FRII radio galaxies. In Dennett-Thorpe et al. (1997) we studied a sample of quasars and found that the high surface brightness regions had flatter spectra on the jet side (explicable as a result of Doppler beaming) whilst the extended regions had spectral asymmetries dependent on lobe length. Unified schemes predict that asymmetries due to beaming will be much smaller in narrow-line radio galaxies than in quasars: we therefore investigate in a similar manner, a sample of radio galaxies with detected jets. We find that spectral asymmetries in these objects are uncorrelated with jet sidedness at all brightness levels, but depend on relative lobe volume. Our results are not in conflict with unified schemes, but suggest that the differences between the two samples are due primarily to power or redshift, rather than to orientation. We also show directly that hotspot spectra steepen as a function of radio power or redshift. Whilst a shift in observed frequency due to the redshift may account for some of the steepening, it cannot account for all of it, and a dependence on radio power is required.

The less depolarized lobe of a radio source is generally the lobe containing the jet (Laing-Garrington correlation) but the less depolarized lobe is also generally that with the flatter radio spectrum (Liu-Pooley correlation). Both effects are strong; taken together they would imply a correlation between jet side and lobe spectral index, i.e. between an orientation-dependent feature and one which is intrinsic. We test this prediction using detailed spectral imaging of a sample of quasars with well-defined jets and investigate whether the result can be reconciled with the standard interpretation of one-sided jets in terms of relativistic aberration. Our central finding is that the spectrum of high surface brightness regions is indeed flatter on the jet side, but that the spectrum of low surface brightness regions is flatter on the side with the longer lobe. We discuss possible causes for these correlations and favour explanations in terms of relativistic bulk motion in the high surface brightness regions and differential synchrotron ageing in the extended lobe material.