Characterizing the epoch of reionization (EoR) at $z\gtrsim 6$ via the redshifted 21 cm line of neutral Hydrogen (HI) is critical to modern astrophysics and cosmology, and thus a key science goal of many current and planned low-frequency radio telescopes. The primary challenge to detecting this signal is the overwhelmingly bright foreground emission at these frequencies, placing stringent requirements on the knowledge of the instruments and inaccuracies in analyses. Results from these experiments have largely been limited not by thermal sensitivity but by systematics, particularly caused by the inability to calibrate the instrument to high accuracy. The interferometric bispectrum phase is immune to antenna-based calibration and errors therein, and presents an independent alternative to detect the EoR HI fluctuations while largely avoiding calibration systematics. Here, we provide a demonstration of this technique on a subset of data from the Hydrogen Epoch of Reionization Array (HERA) to place approximate constraints on the brightness temperature of the intergalactic medium (IGM). From this limited data, at $z=7.7$ we infer "$1\sigma$" upper limits on the IGM brightness temperature to be $\le 316$ "pseudo" mK at $\kappa_\parallel=0.33$ "pseudo" $h$ Mpc$^{-1}$ (data-limited) and $\le 1000$ "pseudo" mK at $\kappa_\parallel=0.875$ "pseudo" $h$ Mpc$^{-1}$ (noise-limited). The "pseudo" units denote only an approximate and not an exact correspondence to the actual distance scales and brightness temperatures. By propagating models in parallel to the data analysis, we confirm that the dynamic range required to separate the cosmic HI signal from the foregrounds is similar to that in standard approaches, and the power spectrum of the bispectrum phase is still data-limited (at $\gtrsim 10^6$ dynamic range) indicating scope for further improvement in sensitivity as the array build-out continues.

The central challenge in 21~cm cosmology is isolating the cosmological signal from bright foregrounds. Many separation techniques rely on the accurate knowledge of the sky and the instrumental response, including the antenna primary beam. For drift-scan telescopes such as the Hydrogen Epoch of Reionization Array \citep[HERA, ][]{DeBoer2017} that do not move, primary beam characterization is particularly challenging because standard beam-calibration routines do not apply \citep{Cornwell2005} and current techniques require accurate source catalogs at the telescope resolution. We present an extension of the method from \citet{Pober2012} where they use beam symmetries to create a network of overlapping source tracks that break the degeneracy between source flux density and beam response and allow their simultaneous estimation. We fit the beam response of our instrument using early HERA observations and find that our results agree well with electromagnetic simulations down to a -20~dB level in power relative to peak gain for sources with high signal-to-noise ratio. In addition, we construct a source catalog with 90 sources down to a flux density of 1.4~Jy at 151~MHz.

The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz.

In 21 cm cosmology, precision calibration is key to the separation of the neutral hydrogen signal from bright astrophysical foregrounds. If not properly calibrated, each antenna element's complex gain can impart spectral structure on otherwise spectrally-smooth foregrounds, complicating that statistical separation. The Hydrogen Epoch of Reionization Array (HERA), an interferometer specialized for 21 cm cosmology and now under construction in South Africa, was designed to be largely calibrated using the self-consistency of repeated measurements of the same interferometric modes. This technique, known as redundant-baseline calibration resolves most of the internal degrees of freedom in the calibration problem. It assumes, however, on antenna elements with identical primary beams placed precisely on a redundant grid. In this work, we review the detailed implementation of the algorithms enabling redundant-baseline calibration and report results with HERA data. We quantify the effects of real-world non-redundancy and how they compare to the idealized scenario in which redundant measurements differ only in their noise realizations. Finally, we study how non-redundancy can produce spurious temporal structure in our calibration solutions---both in data and in simulations---and present strategies for mitigating that structure.

We discuss absolute calibration strategies for Phase I of the Hydrogen Epoch of Reionization Array (HERA), which aims to measure the cosmological 21 cm signal from the Epoch of Reionization (EoR). HERA is a drift-scan array with a 10 degree wide field of view, meaning bright, well-characterized point source transits are scarce. This, combined with HERA's redundant sampling of the uv plane and the modest angular resolution of the Phase I instrument, make traditional sky-based and self-calibration techniques difficult to implement with high dynamic range. Nonetheless, in this work we demonstrate calibration for HERA using point source catalogues and electromagnetic simulations of its primary beam. We show that unmodeled diffuse flux and instrumental contaminants can corrupt the gain solutions, and present a gain smoothing approach for mitigating their impact on the 21 cm power spectrum. We also demonstrate a hybrid sky and redundant calibration scheme and compare it to pure sky-based calibration, showing only a marginal improvement to the gain solutions at intermediate delay scales. Our work suggests that the HERA Phase I system can be well-calibrated for a foreground-avoidance power spectrum estimator by applying direction-independent gains with a small set of degrees of freedom across the frequency and time axes.

We present a study of internal reflection and cross coupling systematics in Phase I of the Hydrogen Epoch of Reionization Array (HERA). In a companion paper, we outlined the mathematical formalism for such systematics and presented algorithms for modeling and removing them from the data. In this work, we apply these techniques to data from HERA's first observing season as a method demonstration. The data show evidence for systematics that, without removal, would hinder a detection of the 21 cm power spectrum for the targeted EoR line-of-sight modes in the range 0.2 < k_parallel < 0.5\ h^-1 Mpc. After systematic removal, we find we can recover these modes in the power spectrum down to the integrated noise-floor of a nightly observation, achieving a dynamic range in the EoR window of 10^-6 in power (mK^2 units) with respect to the bright galactic foreground signal. In the absence of other systematics and assuming the systematic suppression demonstrated here continues to lower noise levels, our results suggest that fully-integrated HERA Phase I may have the capacity to set competitive upper limits on the 21 cm power spectrum. For future observing seasons, HERA will have upgraded analog and digital hardware to better control these systematics in the field.

The detection of the Epoch of Reionization (EoR) delay power spectrum using a "foreground avoidance method" highly depends on the instrument chromaticity. The systematic effects induced by the radio-telescope spread the foreground signal in the delay domain, which contaminates the EoR window theoretically observable. Therefore, it is essential to understand and limit these chromatic effects. This paper describes a method to simulate the frequency and time responses of an antenna, by simultaneously taking into account the analogue RF receiver, the transmission cable, and the mutual coupling caused by adjacent antennas. Applied to the Hydrogen Epoch of Reionization Array (HERA), this study reveals the presence of significant reflections at high delays caused by the 150-m cable which links the antenna to the back-end. Besides, it shows that waves can propagate from one dish to another one through large sections of the array because of mutual coupling. In this more realistic approach, the simulated system time response is attenuated by a factor $10^{4}$ after a characteristic delay which depends on the size of the array and on the antenna position. Ultimately, the system response is attenuated by a factor $10^{5}$ after 1400 ns because of the reflections in the cable, which corresponds to characterizable ${k_\parallel}$-modes above 0.7 $h \rm{Mpc}^{-1}$ at 150 MHz. Thus, this new study shows that the detection of the EoR signal with HERA Phase I will be more challenging than expected. On the other hand, it improves our understanding of the telescope, which is essential to mitigate the instrument chromaticity.

Foreground power dominates the measurements of interferometers that seek a statistical detection of highly-redshifted HI emission from the Epoch of Reionization (EoR). The chromaticity of the instrument creates a boundary in the Fourier transform of frequency (proportional to $k_\parallel$) between spectrally smooth emission, characteristic of the strong synchrotron foreground (the "wedge"), and the spectrally structured emission from HI in the EoR (the "EoR window"). Faraday rotation can inject spectral structure into otherwise smooth polarized foreground emission, which through instrument effects or miscalibration could possibly pollute the EoR window. Using data from the HERA 19-element commissioning array, we investigate the polarization response of this new instrument in the power spectrum domain. We perform a simple image-based calibration based on the unpolarized diffuse emission of the Global Sky Model, and show that it achieves qualitative redundancy between the nominally-redundant baselines of the array and reasonable amplitude accuracy. We construct power spectra of all fully polarized coherencies in all pseudo-Stokes parameters. We compare to simulations based on an unpolarized diffuse sky model and detailed electromagnetic simulations of the dish and feed, confirming that in Stokes I, the calibration does not add significant spectral structure beyond the expected level. Further, this calibration is stable over the 8 days of observations considered. Excess power is seen in the power spectra of the linear polarization Stokes parameters which is not easily attributable to leakage via the primary beam, and results from some combination of residual calibration errors and actual polarized emission. Stokes V is found to be highly discrepant from the expectation of zero power, strongly pointing to the need for more accurate polarized calibration.

Radio Frequency Interference (RFI) is an ever-present limiting factor among radio telescopes even in the most remote observing locations. When looking to retain the maximum amount of sensitivity and reduce contamination for Epoch of Reionization studies, the identification and removal of RFI is especially important. In addition to improved RFI identification, we must also take into account computational efficiency of the RFI-Identification algorithm as radio interferometer arrays such as the Hydrogen Epoch of Reionization Array grow larger in number of receivers. To address this, we present a Deep Fully Convolutional Neural Network (DFCN) that is comprehensive in its use of interferometric data, where both amplitude and phase information are used jointly for identifying RFI. We train the network using simulated HERA visibilities containing mock RFI, yielding a known "ground truth" dataset for evaluating the accuracy of various RFI algorithms. Evaluation of the DFCN model is performed on observations from the 67 dish build-out, HERA-67, and achieves a data throughput of 1.6$\times 10^{5}$ HERA time-ordered 1024 channeled visibilities per hour per GPU. We determine that relative to an amplitude only network including visibility phase adds important adjacent time-frequency context which increases discrimination between RFI and Non-RFI. The inclusion of phase when predicting achieves a Recall of 0.81, Precision of 0.58, and $F_{2}$ score of 0.75 as applied to our HERA-67 observations.

We introduce a new method for performing robust Bayesian estimation of the three-dimensional spatial power spectrum at the Epoch of Reionization (EoR), from interferometric observations. The versatility of this technique allows us to present two approaches. First, when the observations span only a small number of independent spatial frequencies ($k$-modes) we sample directly from the spherical power spectrum coefficients that describe the EoR signal realisation. Second, when the number of $k$-modes to be included in the model becomes large, we sample from the joint probability density of the spherical power spectrum and the signal coefficients, using Hamiltonian Monte Carlo methods to explore this high dimensional ($\sim$ 20000) space efficiently. This approach has been successfully applied to simulated observations that include astrophysically realistic foregrounds in a companion publication (Sims et al. 2016). Here we focus on explaining the methodology in detail, and use simple foreground models to both demonstrate its efficacy, and highlight salient features. In particular, we show that including an arbitrary flat spectrum continuum foreground that is $10^8$ times greater in power than the EoR signal has no detectable impact on our parameter estimates of the EoR power spectrum recovered from the data.

New massively redundant low frequency arrays allow for a novel investigation of closure relations in interferometry. We employ commissioning data from the Hydrogen Epoch of Reionization Array to investigate closure quantities in this densely packed grid array of 14m antennas operating at 100 MHz to 200 MHz. We investigate techniques that utilize closure phase spectra for redundant triads to estimate departures from redundancy for redundant baseline visibilities. We find a median absolute deviation from redundancy in closure phase across the observed frequency range of about 4.5deg. This value translates into a non-redundancy per visibility phase of about 2.6deg, using prototype electronics. The median absolute deviations from redundancy decrease with longer baselines. We show that closure phase spectra can be used to identify ill-behaved antennas in the array, independent of calibration. We investigate the temporal behavior of closure spectra. The Allan variance increases after a one minute stride time, due to passage of the sky through the primary beam of the transit telescope. However, the closure spectra repeat to well within the noise per measurement at corresponding local sidereal times (LST) from day to day. In future papers in this series we will develop the technique of using closure phase spectra in the search for the HI 21cm signal from cosmic reionization.

The Arcminute Microkelvin Imager (AMI) telescopes located at the Mullard Radio Astronomy Observatory near Cambridge have been significantly enhanced by the implementation of a new digital correlator with 1.2 MHz spectral resolution. This system has replaced a 750-MHz resolution analogue lag-based correlator, and was designed to mitigate the effects of radio frequency interference, particularly from geostationary satellites that contaminate observations at low declinations. The upgraded instrument consists of 18 ROACH2 Field Programmable Gate Array platforms used to implement a pair of real-time FX correlators -- one for each of AMI's two arrays. The new system separates the down-converted RF baseband signal from each AMI receiver into two 2.3 GHz-wide sub-bands which are each digitized at 5-Gsps with 8 bits of precision. These digital data streams are filtered into 2048 frequency channels and cross-correlated using FPGA hardware, with a commercial 10 Gb Ethernet switch providing high-speed data interconnect. Images formed using data from the new digital correlator show over an order of magnitude improvement in dynamic range over the previous system. The ability to observe at low declinations has also been significantly improved.

The experimental efforts to detect the redshifted 21 cm signal from the Epoch of Reionization (EoR) are limited predominantly by the chromatic instrumental systematic effect. The delay spectrum methodology for 21 cm power spectrum measurements brought new attention to the critical impact of an antenna's chromaticity on the viability of making this measurement. This methodology established a straightforward relationship between time-domain response of an instrument and the power spectrum modes accessible to a 21 cm EoR experiment. We examine the performance of a prototype of the Hydrogen Epoch of Reionization Array (HERA) array element that is currently observing in Karoo desert, South Africa. We present a mathematical framework to derive the beam integrated frequency response of a HERA prototype element in reception from the return loss measurements between 100-200 MHz and determined the extent of additional foreground contamination in the delay space. The measurement reveals excess spectral structures in comparison to the simulation studies of the HERA element. Combined with the HERA data analysis pipeline that incorporates inverse covariance weighting in optimal quadratic estimation of power spectrum, we find that in spite of its departure from the simulated response, HERA prototype element satisfies the necessary criteria posed by the foreground attenuation limits and potentially can measure the power spectrum at spatial modes as low as $k_{\parallel} > 0.1h$~Mpc$^{-1}$. The work highlights a straightforward method for directly measuring an instrument response and assessing its impact on 21 cm EoR power spectrum measurements for future experiments that will use reflector-type antenna.

The Hydrogen Epoch of Reionization Array (HERA) is a staged experiment to measure 21 cm emission from the primordial intergalactic medium (IGM) throughout cosmic reionization ($z=6-12$), and to explore earlier epochs of our Cosmic Dawn ($z\sim30$). During these epochs, early stars and black holes heated and ionized the IGM, introducing fluctuations in 21 cm emission. HERA is designed to characterize the evolution of the 21 cm power spectrum to constrain the timing and morphology of reionization, the properties of the first galaxies, the evolution of large-scale structure, and the early sources of heating. The full HERA instrument will be a 350-element interferometer in South Africa consisting of 14-m parabolic dishes observing from 50 to 250 MHz. Currently, 19 dishes have been deployed on site and the next 18 are under construction. HERA has been designated as an SKA Precursor instrument. In this paper, we summarize HERA's scientific context and provide forecasts for its key science results. After reviewing the current state of the art in foreground mitigation, we use the delay-spectrum technique to motivate high-level performance requirements for the HERA instrument. Next, we present the HERA instrument design, along with the subsystem specifications that ensure that HERA meets its performance requirements. Finally, we summarize the schedule and status of the project. We conclude by suggesting that, given the realities of foreground contamination, current-generation 21 cm instruments are approaching their sensitivity limits. HERA is designed to bring both the sensitivity and the precision to deliver its primary science on the basis of proven foreground filtering techniques, while developing new subtraction techniques to unlock new capabilities. The result will be a major step toward realizing the widely recognized scientific potential of 21 cm cosmology.

We construct foreground simulations comprising spatially correlated extragalactic and diffuse Galactic emission components and calculate the `intrinsic' (instrument-free) two-dimensional spatial power spectrum and the cylindrically and spherically averaged three-dimensional k-space power spectra of the Epoch of Reionization (EoR) and our foreground simulations using a Bayesian power spectral estimation framework. This leads us to identify a model dependent region of optimal signal estimation for our foreground and EoR models, within which the spatial power in the EoR signal relative to foregrounds is maximised. We identify a target field dependent region, in k-space, of intrinsic foreground power spectral contamination at low k_perp and k_parallel and a transition to a relatively foreground-free intrinsic EoR window in the complement to this region. The contaminated region of k-space demonstrates that simultaneous estimation of the EoR and foregrounds is important for obtaining statistically robust estimates of the EoR power spectrum; biased results will be obtained from methodologies that ignore their covariance. Using simulated observations with frequency dependent uv-coverage and primary beam, with the former derived for HERA in 37-antenna and 331-antenna configuration, we recover instrumental power spectra consistent with their intrinsic counterparts. We discuss the implications of these results for optimal strategies for unbiased estimation of the EoR power spectrum.

Scaling relations for globular clusters (GC) differ from scaling relations for pressure supported (elliptical) galaxies. We show that two-body relaxation is the dominant mechanism in shaping the bivariate dependence of density on mass and Galactocentric distance for Milky Way GCs with masses <10^6 Msun, and it is possible, but not required, that GCs formed with similar scaling relations as ultra-compact dwarf galaxies. We use a fast cluster evolution model to fit a parameterised model for the initial properties of Milky Way GCs to the observed present-day properties. The best-fit cluster initial mass function is substantially flatter (power-law index alpha =- 0.6+/-0.2) than what is observed for young massive clusters (YMCs) forming in the nearby Universe (alpha =~-2). A slightly steeper CIMF is allowed when considering the metal-rich GCs separately (alpha =~-1.2+/-0.4$). If stellar mass loss and two-body relaxation in the Milky Way tidal field are the dominant disruption mechanisms, then GCs formed differently from YMCs.

We have coupled a fast, parametrized star cluster evolution code to a Markov Chain Monte Carlo code to determine the distribution of probable initial conditions of observed star clusters, which may serve as a starting point for future $N$-body calculations. In this paper we validate our method by applying it to a set of star clusters which have been studied in detail numerically with $N$-body simulations and Monte Carlo methods: the Galactic globular clusters M4, 47 Tucanae, NGC 6397, M22, $\omega$ Centauri, Palomar 14 and Palomar 4, the Galactic open cluster M67, and the M31 globular cluster G1. For each cluster we derive a distribution of initial conditions that, after evolution up to the cluster's current age, evolves to the currently observed conditions. We find that there is a connection between the morphology of the distribution of initial conditions and the dynamical age of a cluster and that a degeneracy in the initial half-mass radius towards small radii is present for clusters which have undergone a core collapse during their evolution. We find that the results of our method are in agreement with $N$-body and Monte Carlo studies for the majority of clusters. We conclude that our method is able to find reliable posteriors for the determined initial mass and half-mass radius for observed star clusters, and thus forms an suitable starting point for modeling an observed cluster\rq{}s evolution.

The extremely high sensitivity and resolution of the Square Kilometre Array (SKA) will be useful for addressing a wide set of themes relevant for cosmology, in synergy with current and future cosmic microwave background (CMB) projects. Many of these themes also have a link with future optical-IR and X-ray observations. We discuss the scientific perspectives for these goals, the instrumental requirements and the observational and data analysis approaches, and identify several topics that are important for cosmology and astrophysics at different cosmic epochs.

A new Bayesian method for the analysis of folded pulsar timing data is presented that allows for the simultaneous evaluation of evolution in the pulse profile in either frequency or time, along with the timing model and additional stochastic processes such as red spin noise, or dispersion measure variations. We model the pulse profiles using `shapelets' - a complete ortho-normal set of basis functions that allow us to recreate any physical profile shape. Any evolution in the profiles can then be described as either an arbitrary number of independent profiles, or using some functional form. We perform simulations to compare this approach with established methods for pulsar timing analysis, and to demonstrate model selection between different evolutionary scenarios using the Bayesian evidence. %s The simplicity of our method allows for many possible extensions, such as including models for correlated noise in the pulse profile, or broadening of the pulse profiles due to scattering. As such, while it is a marked departure from standard pulsar timing analysis methods, it has clear applications for both new and current datasets, such as those from the European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA).

We introduce a method for performing a robust Bayesian analysis of non-Gaussianity present in pulsar timing data, simultaneously with the pulsar timing model, and additional stochastic parameters such as those describing red spin noise and dispersion measure variations. The parameters used to define the presence of non-Gaussianity are zero for Gaussian processes, giving a simple method of defining the strength of non-Gaussian behaviour. We use simulations to show that assuming Gaussian statistics when the noise in the data is drawn from a non-Gaussian distribution can significantly increase the uncertainties associated with the pulsar timing model parameters. We then apply the method to the publicly available 15 year Parkes Pulsar Timing Array data release 1 dataset for the binary pulsar J0437$-$4715. In this analysis we present a significant detection of non-Gaussianity in the uncorrelated non-thermal noise, but we find that it does not yet impact the timing model or stochastic parameter estimates significantly compared to analysis performed assuming Gaussian statistics. The methods presented are, however, shown to be of immediate practical use for current European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA) datasets.

A new Bayesian method for performing an image domain search for line-emitting galaxies is presented. The method uses both spatial and spectral information to robustly determine the source properties, employing either simple Gaussian, or other physically motivated models whilst using the evidence to determine the probability that the source is real. In this paper, we describe the method, and its application to both a simulated data set, and a blind survey for cold molecular gas using observations of the Hubble Deep Field North taken with the Plateau de Bure Interferometer. We make a total of 6 robust detections in the survey, 5 of which have counterparts in other observing bands. We identify the most secure detections found in a previous investigation, while finding one new probable line source with an optical ID not seen in the previous analysis. This study acts as a pilot application of Bayesian statistics to future searches to be carried out both for low-$J$ CO transitions of high redshift galaxies using the JVLA, and at millimeter wavelengths with ALMA, enabling the inference of robust scientific conclusions about the history of the molecular gas properties of star-forming galaxies in the Universe through cosmic time.

There is mounting evidence that mechanical radio source feedback is important in galaxy evolution and in order to quantify this feedback, detailed models of radio source evolution are required. We present an extension to current analytic models that encompasses young radio sources with physical sizes on sub-kiloparsec scales. This work builds on an existing young source dynamical model to include radiative losses in a flat environment, and as such, is the best physically-motivated Compact Symmetric Object model to date. Results predict that young radio sources experience significant radiative loss on length scales and spectral scales consistent with observed Compact Steep-Spectrum sources. We include full expressions for the transition to self-similar expansion and present this complete model of radio source evolution from first cocoon formation to end of source lifetime around 10^8 years within the context of a simplified King profile external atmosphere.

We present a new version of the fast star cluster evolution code Evolve Me A Cluster of StarS (EMACSS). While previous versions of EMACSS reproduced clusters of single-mass stars, this version models clusters with an evolving stellar content. Stellar evolution dominates early evolution, and leads to: (1) reduction of the mean mass of stars due to the mass loss of high-mass stars; (2) expansion of the half-mass radius; (3) for (nearly) Roche Volume filling clusters, the induced escape of stars. Once sufficient relaxation has occurred (~ 10 relaxation times-scales), clusters reach a second, 'balanced' state whereby the core releases energy as required by the cluster as a whole. In this state: (1) stars escape due to tidal effects faster than before balanced evolution; (2) the half-mass radius expands or contracts depending on the Roche volume filling factor; and (3) the mean mass of stars increases due to the preferential ejection of low-mass stars. We compare the EMACSS results of several cluster properties against N-body simulations of clusters spanning a range of initial number of stars, mass, half-mass radius, and tidal environments, and show that our prescription accurately predicts cluster evolution for this database. Finally, we consider applications for EMACSS, such as studies of galactic globular cluster populations in cosmological simulations.

We introduce version two of the fast star cluster evolution code Evolve Me A Cluster of StarS (EMACSS). The first version (Alexander & Gieles) assumed that cluster evolution is balanced for the majority of the life-cycle, meaning that the rate of energy generation in the core of the cluster equals the diffusion rate of energy by two-body relaxation, which makes the code suitable for modelling clusters in weak tidal fields. In this new version we extend the model to include an unbalanced phase of evolution to describe the pre-collapse evolution and the accompanying escape rate such that clusters in strong tidal fields can also be modelled. We also add a prescription for the evolution of the core radius and density and a related cluster concentration parameter. The model simultaneously solves a series of first-order ordinary differential equations for the rate of change of the core radius, half-mass radius and the number of member stars N. About two thousand integration steps in time are required to solve for the entire evolution of a star cluster and this number is approximately independent of N. We compare the model to the variation of these parameters following from a series of direct N-body calculations of single-mass clusters and find good agreement in the evolution of all parameters. Relevant time-scales, such as the total lifetimes and core collapse times, are reproduced with an accuracy of about 10% for clusters with various initial half-mass radii (relative to their Jacobi radii) and a range of different initial N up to N = 65536. We intend to extend this framework to include more realistic initial conditions, such as a stellar mass spectrum and mass loss from stars. The EMACSS code can be used in star cluster population studies and in models that consider the co-evolution of (globular) star clusters and large scale structures.

Here we present a Bayesian method of including discrete measurements of dispersion measure due to the interstellar medium in the direction of a pulsar as prior information in the analysis of that pulsar. We use a simple simulation to show the efficacy of this method, where the inclusion of the additional measurements results in both a significant increase in the precision with which the timing model parameters can be obtained, and an improved upper limit on the amplitude of any red noise in the dataset. We show that this method can be applied where no multi-frequency data exists across much of the dataset, and where there is no simultaneous multi-frequency data for any given observing epoch. Including such information in the analysis of upcoming International Pulsar Timing Array (IPTA) and European Pulsar Timing Array (EPTA) data releases could therefore prove invaluable in obtaining the most constraining limits on gravitational wave signals within those datasets.

A new Bayesian software package for the analysis of pulsar timing data is presented in the form of TempoNest which allows for the robust determination of the non-linear pulsar timing solution simultaneously with a range of additional stochastic parameters. This includes both red spin noise and dispersion measure variations using either power law descriptions of the noise, or through a model-independent method that parameterises the power at individual frequencies in the signal. We use TempoNest to show that at noise levels representative of current datasets in the European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA) the linear timing model can underestimate the uncertainties of the timing solution by up to an order of magnitude. We also show how to perform Bayesian model selection between different sets of timing model and stochastic parameters, for example, by demonstrating that in the pulsar B1937+21 both the dispersion measure variations and spin noise in the data are optimally modelled by simple power laws. Finally we show that not including the stochastic parameters simultaneously with the timing model can lead to unpredictable variation in the estimated uncertainties, compromising the robustness of the scientific results extracted from such analysis.

A new model independent method is presented for the analysis of pulsar timing data and the estimation of the spectral properties of an isotropic gravitational wave background (GWB). We show that by rephrasing the likelihood we are able to eliminate the most costly aspects of computation normally associated with this type of data analysis. When applied to the International Pulsar Timing Array Mock Data Challenge data sets this results in speedups of approximately 2 to 3 orders of magnitude compared to established methods. We present three applications of the new likelihood. In the low signal to noise regime we sample directly from the power spectrum coefficients of the GWB signal realization. In the high signal to noise regime, where the data can support a large number of coefficients, we sample from the joint probability density of the power spectrum coefficients for the individual pulsars and the GWB signal realization. Critically in both these cases we need make no assumptions about the form of the power spectrum of the GWB, or the individual pulsars. Finally we present a method for characterizing the spatial correlation between pulsars on the sky, making no assumptions about the form of that correlation, and therefore providing the only truly general Bayesian method of confirming a GWB detection from pulsar timing data.

Studies of extra-galactic globular clusters have shown that the peak size of the globular cluster (GC) radius distribution (RD) depends only weakly on galactic environment, and can be used as a standard ruler. We model RDs of GC populations using a simple prescription for a Hubble time of relaxation driven evolution of cluster mass and radius, and explore the conditions under which the RD can be used as a standard ruler. We consider a power-law cluster initial mass function (CIMF) with and without an exponential truncation, and focus in particular on a flat and a steep CIMF (power-law indices of 0 and -2, respectively). For the initial half-mass radii at birth we adopt either Roche-lobe filling conditions ('filling',meaning that the ratio of half-mass to Jacobi radius is approximately rh/rJ ~ 0.15) or strongly Roche-lobe under-filling conditions ('under-filling', implying that initially rh/rJ << 0.15). Assuming a constant orbital velocity about the galaxy centre we find for a steep CIMF that the typical half-light radius scales with galactocentric radius RG as RG^1/3. This weak scaling is consistent with observations, but this scenario has the (well known) problem that too many low-mass clusters survive. A flat CIMF with 'filling' initial conditions results in the correct mass function at old ages, but with too many large (massive) clusters at large RG. An 'underfilling' GC population with a flat CIMF also results in the correct mass function, and can also successfully reproduce the shape of the RD, with a peak size that is (almost) independent of RG. In this case, the peak size depends (almost) only on the peak mass of the GC mass function. The (near) universality of the GC RD is therefore because of the (near) universality of the CIMF. There are some extended GCs in the outer halo of the Milky Way that cannot be explained by this model.

We use sensitive observations of three high redshift sources; [CII] fine structure and CO(2-1) rotational transitions for the z=6.4 Quasar host galaxy (QSO) J1148+5251, and [CII] and CO(5-4) transitions from the QSO BR1202-0725 and its sub-millimeter companion (SMG) galaxy at z=4.7. We use these observations to place constraints on the quantity Dz = z(CO) - z(CII) for each source where z(CO) and z(CII) are the observed redshifts of the CO rotational transition and [CII] fine structure transition respectively, using a combination of approaches; 1) By modelling the emission line profiles using `shapelets' to compare both the emission redshifts and the line profiles themselves, in order to make inferences about the intrinsic velocity differences between the molecular and atomic gas, and 2) By performing a marginalisation over all model parameters in order to calculate a non-parametric estimate of Dz. We derive 99% confidence intervals for the marginalised posterior of Dz of (-1.9 pm 1.3) x10^-3, (-3 pm 8) x10^-4 and (-2 pm 4) x10^-3 for J1148+5251, and the BR1202-0725 QSO and SMG respectively. We show the [CII] and CO(2-1) line profiles for J1148+5251 are consistent with each other within the limits of the data, whilst the [CII] and CO(5-4) line profiles from the BR1202-0725 QSO and SMG respectively have 65 and >99.9% probabilities of being inconsistent, with the CO(5-4) lines ~ 30% wider than the [CII] lines. Therefore whilst the observed values of Dz can correspond to variations in the quantity Delta F/F with cosmic time, where F=alpha^2/mu, with alpha the fine structure constant, and mu the proton-to-electron mass ratio, of both (-3.3 pm 2.3) x10^-4 for a look back time of 12.9 Gyr and of (-5 pm 15) x10^-5 for a look back time of 12.4 Gyr we propose that they are the result of the two species of gas being spatially separated as indicated by the inconsistencies in their line profiles.

The Arcminute Microkelvin Imager (AMI) is a telescope specifically designed for high sensitivity measurements of low-surface-brightness features at cm-wavelength and has unique, important capabilities. It consists of two interferometer arrays operating over 13.5-18 GHz that image structures on scales of 0.5-10 arcmin with very low systematics. The Small Array (AMI-SA; ten 3.7-m antennas) couples very well to Sunyaev-Zel'dovich features from galaxy clusters and to many Galactic features. The Large Array (AMI-LA; eight 13-m antennas) has a collecting area ten times that of the AMI-SA and longer baselines, crucially allowing the removal of the effects of confusing radio point sources from regions of low surface-brightness, extended emission. Moreover AMI provides fast, deep object surveying and allows monitoring of large numbers of objects. In this White Paper we review the new science - both Galactic and extragalactic - already achieved with AMI and outline the prospects for much more.

In the case of an initially conical jet, we study the relation between jet collimation by the external pressure and large-scale morphology. We first consider the important length-scales in the problem, and then carry out axisymmetric hydrodynamic simulations that include, for certain parameters, all these length-scales. We find three important scales related to the collimation region: (i) where the sideways ram-pressure equals the external pressure, (ii) where the jet density equals the ambient density, and (iii) where the forward ram-pressure falls below the ambient pressure. These scales are set by the external Mach-number and opening angle of the jet. We demonstrate that the relative magnitudes of these scales determine the collimation, Mach-number, density and morphology of the large scale jet. Based on analysis of the shock structure, we reproduce successfully the morphology of Fanaroff-Riley (FR) class I and II radio sources. Within the framework of the model, an FR I radio source must have a large intrinsic opening angle. Entrainment of ambient gas might also be important. We also show that all FR I sources with radio lobes or similar features must have had an earlier FR II phase.

The James Clerk Maxwell Telescope Nearby Galaxies Legacy Survey (NGLS) comprises an HI-selected sample of 155 galaxies spanning all morphological types with distances less than 25 Mpc. We describe the scientific goals of the survey, the sample selection, and the observing strategy. We also present an atlas and analysis of the CO J=3-2 maps for the 47 galaxies in the NGLS which are also part of the Spitzer Infrared Nearby Galaxies Survey. We find a wide range of molecular gas mass fractions in the galaxies in this sample and explore the correlation of the far-infrared luminosity, which traces star formation, with the CO luminosity, which traces the molecular gas mass. By comparing the NGLS data with merging galaxies at low and high redshift which have also been observed in the CO J=3-2 line, we show that the correlation of far-infrared and CO luminosity shows a significant trend with luminosity. This trend is consistent with a molecular gas depletion time which is more than an order of magnitude faster in the merger galaxies than in nearby normal galaxies. We also find a strong correlation of the L(FIR)/L(CO3-2) ratio with the atomic to molecular gas mass ratio. This correlation suggests that some of the far-infrared emission originates from dust associated with atomic gas and that its contribution is particularly important in galaxies where most of the gas is in the atomic phase.

We introduce the star cluster evolution code Evolve Me A Cluster of StarS (EMACSS), a simple yet physically motivated computational model that describes the evolution of some fundamental properties of star clusters in static tidal fields. We base our prescription upon the flow of energy within the cluster, which is a constant fraction of the total energy per half-mass relaxation time. According to Henon's predictions, this flow is independent of the precise mechanisms for energy production within the core, and therefore does not require a complete description of the many-body interactions therein. For a cluster of equal-mass stars, we thence use dynamical theory and analytic descriptions of escape mechanisms to construct a series of coupled differential equations expressing the time evolution of cluster mass and radius. These equations are numerically solved using a fourth-order Runge-Kutta integration kernel, and the results benchmarked against a data base of direct N-body simulations. We use simulations containing a modest initial number of stars (1024 <= N <= 65 536) and point-mass tidal fields of various strengths. Our prescription is publicly available and reproduces the N-body results to within ~10 per cent accuracy for the entire post-collapse evolution of star clusters.

We present 3-D MHD and synthetic numerical simulations to follow the evolution of randomly magnetized intra-cluster medium plasma under the effects of powerful, light, hypersonic and bipolar jets. We prescribe the cluster magnetic field (CMF) as a Gaussian random field with power law energy spectrum tuned to the expectation for Kolmogorov turbulence. We investigate the power of jets and the viewing angle used for the synthetic Rotation Measure (RM) observations. We find the model radio sources introduce and amplify fluctuations on the RM statistical properties; the average RM and the RM standard deviation are increased by the action of the jets. This may lead to overestimations of the CMFs' strength up to 70%. The effect correlates with the jet power. Jets distort and amplify CMFs especially near the edges of the lobes and the jets' heads. Thus the RM structure functions are flattened at scales comparable to the source size. Jet-produced RM enhancements depend on the orientation of the jet axis to the line of sight and are thus more apparent in quasars than in radio galaxies.

A combination of three-dimensional (3D) magnetohydrodynamics (MHD) and synthetic numerical simulations are presented to follow the evolution of a randomly magnetised plasma that models the intra-cluster medium (ICM), under the isolated effects of powerful, light, hypersonic and bipolar Fanaroff-Riley class II (FR II) jets. We prescribe the cluster magnetic field (CMF) as a Gaussian random field with a Kolmogorov-like energy spectrum. Both the power of the jets and the viewing angle that is used for the synthetic Rotation Measure (RM) observations are investigated. We find the model radio sources introduce and amplify fluctuations on the RM statistical properties which we analyse as a function of time as well as the viewing angle. The average RM and the RM standard deviation are increased by the action of the jets. Energetics, RM statistics and magnetic power spectral analysis consistently show that the effects also correlate with the jets' power, and that the lightest, fastest jets produce the strongest changes in their environment. We see jets distort and amplify the CMFs especially near the edges of the lobes and the jets' heads. This process leads to a flattening of the RM structure functions at scales comparable to the source size. The edge features we find are similar to ones observed in Hydra A. The results show that jet-produced RM enhancements are more apparent in quasars than in radio galaxies. Globally, jets tend to enhance the RM standard deviation which may lead to overestimations of the CMFs' strength by about 70%. This study means to serve as a pathfinder for the SKA, EVLA and LOFAR to follow the evolution of cosmic magnetic fields.

Radio observations of Fanaroff-Riley class II sources often show correlations between the synchrotron emission and the linear-polarimetric distributions. Magnetic position vectors seem to align with the projected emission of both the radio jets and the sources' edges. Using statistics we study such relation as well as its unknown time evolution via synthetic polarisation maps of model FR II sources formed in 3D-MHD numerical simulations of bipolar, hypersonic and weakly magnetised jets. The magnetic field is initially random with a Kolmogorov power spectrum, everywhere. We investigate the structure and evolution of magnetic fields in the sources as a function of the power of jets and the observational viewing angle. Our synthetic polarisation maps agree with observations, showing B-field vectors which are predominantly aligned with the jet axis, and show that magnetic fields inside sources are shaped by the jets' backflow. Polarimetry is found to correlate with time, the viewing angle and the jet-to-ambient density contrast. The magnetic structure inside thin elongated sources is more uniform than inside more spherical ones. We see jets increase the magnetic energy in cocoons in proportion to the jet velocity and the cocoon width. Filaments in the synthetic emission maps suggest turbulence develops in evolved sources.

3D-MHD numerical simulations of bipolar, hypersonic, weakly magnetized jets and synthetic synchrotron observations are presented to study the structure and evolution of magnetic fields in FR II radio sources. The magnetic field setup in the jet is initially random. The power of the jets as well as the observational viewing angle are investigated. We find that synthetic polarization maps agree with observations and show that magnetic fields inside the sources are shaped by the jets' backflow. Polarimetry statistics correlates with time, the viewing angle and the jet-to-ambient density contrast. The magnetic structure inside thin elongated sources is more uniform than for ones with fatter cocoons. Jets increase the magnetic energy in cocoons, in proportion to the jet velocity. Both, filaments in synthetic emission maps and 3D magnetic power spectra suggest that turbulence develops in evolved sources.

We present further observations of the Lockman Hole field, made with the Giant Metrewave Radio Telescope at 610 MHz with a resolution of 6 x 5 arcsec^2. These complement our earlier observations of the central approx 5 deg^2 by covering a further approx 8 deg^2, with an r.m.s. noise down to ~80 microJy beam^-1. A catalogue of 4934 radio sources is presented.

We develop an analytical model to follow the cosmological evolution of magnetic fields in disk galaxies. Our assumption is that fields are amplified from a small seed field via magnetohydrodynamical (MHD) turbulence. We further assume that this process is fast compared to other relevant timescales, and occurs principally in the cold disk gas. We follow the turbulent energy density using the Shabala & Alexander (2009) galaxy formation and evolution model. Three processes are important to the turbulent energy budget: infall of cool gas onto the disk and supernova feedback increase the turbulence; while star formation removes gas and hence turbulent energy from the cold gas. Finally, we assume that field energy is continuously transferred from the incoherent random field into an ordered field by differential galactic rotation. Model predictions are compared with observations of local late type galaxies by Fitt & Alexander (1993) and Shabala et al. (2008). The model reproduces observed magnetic field strengths and luminosities in low and intermediate-mass galaxies. These quantities are overpredicted in the most massive hosts, suggesting that inclusion of gas ejection by powerful AGNs is necessary in order to quench gas cooling and reconcile the predicted and observed magnetic field strengths.

We studied Faraday rotation measure (RM) in turbulent media with the rms Mach number of unity, using isothermal, magnetohydrodynamic turbulence simulations. Four cases with different values of initial plasma beta were considered. Our main findings are as follows. (1) There is no strong correlation between the fluctuations of magnetic field strength and gas density. So the magnetic field strength estimated with RM/DM (DM is the dispersion measure) correctly represents the true mean strength of the magnetic field along the line of sight. (2) The frequency distribution of RMs is well fitted to the Gaussian. In addition, there is a good relation between the width of the distribution of RM/$\bar{\rm RM}$ ($\bar{\rm RM}$ is the average value of RMs) and the strength of the regular field along the line of sight; the width is narrower, if the field strength is stronger. We discussed the implications of our findings in the warm ionized medium where the Mach number of turbulent motions is around unity.

We investigate the potential of the Square Kilometre Array (SKA) for measuring the magnetic fields in clusters of galaxies via Faraday rotation of background polarised sources. [...] We find that about 10 per cent of the sky is covered by a significant extragalactic Faraday screen. Most of it has rotation measures between 10 and 100 rad/m/m. We argue that the cluster centres should have up to about 5000 rad/m/m. We show that the proposed mid frequency aperture array of the SKA as well as the lowest band of the SKA dish array are well suited to make measurements for most of these rotation measure values, typically requiring a signal-to-noise of ten. We calculate the spacing of sources forming a grid for the purpose of measuring foreground rotation measures: it reaches a spacing of 36 arcsec for a 100 hour SKA observation per field. We also calculate the statistics for background RM measurements in clusters of galaxies. We find that a first phase of the SKA would allow us to take stacking experiments out to high redshifts (>1), and provide improved magnetic field structure measurements for individual nearby clusters. The full SKA aperture array would be able to make very detailed magnetic field structure measurements of clusters with more than 100 background sources per cluster up to a redshift of 0.5 and more than 1000 background sources per cluster for nearby clusters, and could for reasonable assumptions about future measurements of electron densities in high redshift clusters constrain the power law index for the magnetic field evolution to better than dm=0.4, if the magnetic field in clusters should follow B ~ (1+z)^m.

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 present a galaxy evolution model which incorporates a physically motivated implementation of AGN feedback. Intermittent jets inflate cocoons of radio plasma which then expand supersonically, shock heating the ambient gas. The model reproduces observed star formation histories to the highest redshifts for which reliable data exists, as well as the observed galaxy colour bimodality. Intermittent radio source feedback also naturally provides a way of keeping the black hole and spheroid growth in step. We find possible evidence for a top-heavy Initial Mass Function (IMF) for $z>2$, consistent with observations of element abundances, and sub-mm and Lyman break galaxy counts.

We probe the relationship between star formation rate (SFR) and radio synchrotron luminosity in galaxies at 0 < z < 2 within the northern Spitzer Wide-area Infrared Extragalactic survey (SWIRE) fields, in order to investigate some of the assumptions that go into calculating the star formation history of the Universe from deep radio observations. We present new 610-MHz Giant Metrewave Radio Telescope (GMRT) observations of the European Large-Area ISO Survey (ELAIS)-N2 field, and using this data, along with previous GMRT surveys carried out in the ELAIS-N1 and Lockman Hole regions, we construct a sample of galaxies which have redshift and SFR information available from the SWIRE survey. We test whether the local relationship between SFR and radio luminosity is applicable to z = 2 galaxies, and look for evolution in this relationship with both redshift and SFR in order to examine whether the physical processes which lead to synchrotron radiation have remained the same since the peak of star formation in the Universe. We find that the local calibration between radio luminosity and star formation can be successfully applied to radio-selected high-redshift, high-SFR galaxies, although we identify a small number of sources where this may not be the case; these sources show evidence for inaccurate estimations of their SFR, but there may also be some contribution from physical effects such as the recent onset of starburst activity, or suppression of the radio luminosity within these galaxies.

For over four decades, synchrotron-radiating sources have played a series of pathfinding roles in the study of galaxy clusters and large scale structure. Such sources are uniquely sensitive to the turbulence and shock structures of large-scale environments, and their cosmic rays and magnetic fields often play important dynamic and thermodynamic roles. They provide essential complements to studies at other wavebands. Over the next decade, they will fill essential gaps in both cluster astrophysics and the cosmological growth of structure in the universe, especially where the signatures of shocks and turbulence, or even the underlying thermal plasma itself, are otherwise undetectable. Simultaneously, synchrotron studies offer a unique tool for exploring the fundamental question of the origins of cosmic magnetic fields. This work will be based on the new generation of m/cm-wave radio telescopes now in construction, as well as major advances in the sophistication of 3-D MHD simulations.

We investigate the infrared / radio correlation using the technique of source stacking, in order to probe the average properties of radio sources that are too faint to be detected individually. We compare the two methods used in the literature to stack sources, and demonstrate that the creation of stacked images leads to a loss of information. We stack infrared sources in the Spitzer extragalactic First Look Survey (xFLS) field, and the three northern Spitzer Wide-area Infrared Extragalactic survey (SWIRE) fields, using radio surveys created at 610 MHz and 1.4 GHz, and find a variation in the absolute strength of the correlation between the xFLS and SWIRE regions, but no evidence for significant evolution in the correlation over the 24-um flux density range 150 uJy - 2 mJy. We carry out the first radio source stacking experiment using 70-um-selected galaxies, and find no evidence for significant evolution over the 70-um flux density range 10 mJy - 100 mJy.

We investigate the effects of intermittent strong jets from an Active Galactic Nuclei (AGN) of a massive galaxy in the core of a cool core galaxy cluster, on the dynamics and metal distribution of the intra-cluster medium (ICM). We use a simple model for the metal distribution within the host galaxy which includes metal injection via star formation. We carry out 2D axisymmetric hydrodynamic simulations of these systems. After having established a cooling flow, two light jets are injected in opposite directions with a range of (intermittent) active phases. We follow the time evolution of the system from the jets' active phases up to 3 Gyr. The general metallicity evolution for all our simulations is very similar on large-scales. The convective flows driven by the jets advect gas and metals from the central galaxy to distances beyond 1.5 Mpc within the cluster. Intermittent jets are able to distribute metals to grater radii. The metal injection has effects on the general metal abundances present in the ICM, the faster the metal replenishment time the higher the metal abundance. Although metallicity gradients in the very central regions of cool core clusters are likely to be shaped by less-energetic phenomena, we find evidence in our simulations for metallicity gradients similar to those observed out to ~400 kpc radius. The metal distribution details of the central galaxy have no effects on the dynamical evolution of the ICM metals.

Infrared-Faint Radio Sources (IFRSs) are a class of source which are bright at radio frequencies, but do not appear in deep infrared images. We report the detection of 14 IFRSs within the Spitzer extragalactic First Look Survey field, eight of which are detected near to the limiting magnitude of a deep R-band image of the region, at R ~ 24.5. Sensitive Spitzer Space Telescope images are stacked in order to place upper limits on their mid-infrared flux densities, and using recent 610-MHz and 1.4-GHz observations we find that they have spectral indices which vary between alpha = 0.05 and 1.38, where we define alpha such that S = S_0 nu^(- alpha), and should not be thought of as a single source population. We place constraints on the luminosity and linear size of these sources, and through comparison with well-studied local objects in the 3CRR catalogue demonstrate that they can be modelled as being compact (< 20 kpc) Fanaroff-Riley Class II (FRII) radio galaxies located at high redshift (z > 4).

As part of an ongoing series of deep GMRT surveys we have observed the Spitzer extragalactic First Look Survey field, producing the deepest wide-field 610-MHz survey published to date. We reach an rms noise of 30 microJy/beam before primary beam correction, with a resolution of ~6 arcsec over an area of ~4 square degrees. By combining these observations with the existing 1.4-GHz VLA survey produced by Condon et al. (2003), along with infrared data in up to seven wavebands from the Spitzer Space Telescope, optical photometry from SDSS and a range of spectroscopic redshift surveys, we are able to study the relationship between radio luminosity and star formation rate in star-forming galaxies up to z ~ 1. The large amount of multi-wavelength data available allows accurate k-corrections to be performed in the radio due to the knowledge of the radio spectral index, and in the infrared through the use of a semi-empirical radiative transfer model which models star-forming regions, warm dust surrounding these regions, and diffuse interstellar dust, taking into account the star formation rate, star formation history and hydrogen column density within each galaxy. A strong correlation is seen between radio luminosity and infrared-derived star formation rates, which is best fit by a slightly non-linear power-law. We look for cosmic evolution in the comparative radio brightness of star-forming galaxies by searching for deviations away from the global relationship. Any such deviation would indicate a systematic variation in one or more of the properties controlling synchrotron radiation, in particular an increase in the magnetic field strengths of star-forming galaxies. The data shows no evidence for such an effect, suggesting that there has been little evolution in the magnetic fields of galaxies since z ~ 1.

We use a volume- and flux-limited sample of local ($0.03 \leq z \leq 0.1$) radio galaxies with optical counterparts to address the question of how long a typical galaxy spends in radio-active and quiescent states. The length of the active phase has a strong dependence on the stellar mass of the host galaxy. Radio sources in the most massive hosts are also retriggered more frequently. The time spent in the active phase has the same dependence on stellar mass as does the gas cooling rate, suggesting the onset of the quiescent phase is due to fuel depletion. We find radio and emission line AGN activity to be independent, consistent with these corresponding to different accretion states.

We present a semi-empirical model for the infrared emission of dust around star-forming sites in galaxies. Our approach combines a simple model of radiative transfer in dust clouds with a state-of-the-art model of the microscopic optical properties of dust grains pioneered by Draine & Li. In combination with the Starburst 99 stellar spectral synthesis package, this framework is able to produce synthetic spectra for galaxies which extend from the Lyman limit through to the far-infrared. We use it to probe how model galaxy spectra depend upon the physical characteristics of their dust grain populations, and upon the energy sources which heat that dust. We compare the predictions of our model with the 8- and 24-micron luminosities of sources in the Spitzer First Look Survey, and conclude by using the models to analyse the relative merits of various colour diagnostics in distinguishing systems out to a redshift of 2 with ongoing star formation from those with only old stellar populations.

We present observations of the Lockman Hole with the Giant Metrewave Radio Telescope (GMRT). Twelve pointings were observed, covering a total area of ~ five square degrees with a resolution of 6 x 5 arcsec, position angle +45deg. The majority of the pointings have an rms noise of 60 microJy/beam before correction for the attenuation of the GMRT primary beam. Techniques used for data reduction and production of a mosaicked image of the region are described, and the final mosaic is presented, along with a catalogue of 2845 sources detected above 6 sigma. Radio source counts are calculated at 610 MHz and combined with existing 1.4-GHz source counts, in order to show that pure luminosity evolution of the local radio luminosity functions for active galactic nuclei and starburst galaxies is sufficient to account for the two source counts simultaneously.

Observations of the ELAIS-N1 field taken at 610 MHz with the Giant Metrewave Radio Telescope are presented. Nineteen pointings were observed, covering a total area of 9 square degrees with a resolution of 6" x 5", PA +45 deg. Four of the pointings were deep observations with an rms of 40 microJy before primary beam correction, with the remaining fifteen pointings having an rms of 70 microJy. The techniques used for data reduction and production of a mosaicked image of the region are described, and the final mosaic is presented, along with a catalogue of 2500 sources detected above 6 sigma. This work complements the large amount of optical and infrared data already available on the region. We calculate 610-MHz source counts down to 270 microJy, and find further evidence for the turnover in differential number counts below 1 mJy, previously seen at both 610 MHz and 1.4 GHz.

Observations of the Spitzer extragalactic First Look Survey field taken at 610 MHz with the Giant Metrewave Radio Telescope are presented. Seven individual pointings were observed, covering an area of 4 square degrees with a resolution of 5.8'' x 4.7'', PA 60 deg. The r.m.s. noise at the centre of the pointings is between 27 and 30 microJy before correction for the GMRT primary beam. The techniques used for data reduction and production of a mosaicked image of the region are described, and the final mosaic, along with a catalogue of 3944 sources detected above 5 sigma, are presented. The survey complements existing radio and infrared data available for this region.

We present data probing the spatial and kinematical distribution of both the atomic (HI) and molecular (CO) gas in NGC 5218, the late-type barred spiral galaxy in the spiral-elliptical interacting pair, Arp 104. We consider these data in conjunction with far-infrared and radio continuum data, and N-body simulations, to study the galaxies interactions, and the star formation properties of NGC 5218. We use these data to assess the importance of the bar and tidal interaction on the evolution of NGC 5218, and the extent to which the tidal interaction may have been important in triggering the bar. The molecular gas distribution of NGC 5218 appears to have been strongly affected by the bar; the distribution is centrally condensed with a very large surface density in the central region. The N-body simulations indicate a timescale since perigalacticon of approximately 3 x 10**8 yr, which is consistent with the interaction having triggered or enhanced the bar potential in NGC 5218, leading to inflow and the large central molecular gas density observed. Whilst NGC 5218 appears to be undergoing active star formation, its star formation efficiency is comparable to a `normal' SBb galaxy. We propose that this system may be on the brink of a more active phase of star formation.

We present 2D axisymmetric hydrodynamic simulations investigating the long term effect of FR II radio galaxies on the metal distribution of the surrounding intra-cluster medium (ICM). A light jet is injected into a cooling flow atmosphere for 10-30 Myr. We then follow the subsequent evolution for 3 Gyr on a spherical grid spanning 3 Mpc in radius. A series of passive tracer particles were placed in an annulus about the cluster core to simulate metal carrying clouds in order to calculate the metallicity (Z) as a function of time and radial distance from the cluster centre. The jet has a significant effect on the ICM over the entire 3 Gyr period. By the end of the simulations, the jets produced metallicities of ~10% of the initial metallicity of the cluster core throughout much of the cluster. The jets transport the metals not only in mixing regions, but also through upwelling ICM behind the jet, enriching the cluster over both long and short distances.

The interaction of optically emitting clouds with warm X-ray gas and hot, tenuous radio plasma in radio jet cocoons is modelled by 2D compressible hydrodynamic simulations. The initial setup is the Kelvin-Helmholtz instability at a contact surface of density contrast 10,000. The denser medium contains clouds of higher density. Optically thin radiation is realised via a cooling source term. The cool phase effectively extracts energy from the other gas which is both, radiated away and used for acceleration of the cold phase. This increases the system's cooling rate substantially and leads to a massively amplified cold mass dropout. We show that it is feasible, given small seed clouds of order 100 solar masses, that all of the optically emitting gas in a radio jet cocoon may be produced by this mechanism on the propagation timescale of the jet. The mass is generally distributed as T^-1/2 with temperature, with a prominent peak at 14,000 K. This peak is likely to be related to the counteracting effects of shock heating and a strong rise in the cooling function. The volume filling factor of cold gas in this peak is of the order 10^-5 to 10^-3 and generally increases during the simulation time. The simulations tend towards an isotropic scale free Kolmogorov-type energy spectrum over the simulation timescale. We find the same Mach-number density relation as Kritsuk & Norman (2004) and show that this relation may explain the velocity widths of emission lines associated with high redshift radio galaxies, if the environmental temperature is lower, or the jet-ambient density ratio is less extreme than in their low redshift counterparts.

We present Giant Metrewave Radio Telescope 21-cm HI observations towards a sample of compact radio sources behind galaxy groups, to search for cool HI. The results -- from high dynamic range spectra for 8 lines-of-sight through 7 galaxy groups -- do not show any evidence for absorption by cool HI. At a resolution of 20 km/s, the optical depth upper limits obtained were between 0.0075 and 0.035 (3sigma); these correspond to upper limits of a few times 10**23 m**-2 for the column density of any cool HI along these lines of sight (assuming a spin temperature of 100 K).

We investigate the atomic and molecular interstellar medium and star formation of NGC 275, the late-type spiral galaxy in Arp 140, which is interacting with NGC 274, an early-type system. The atomic gas (HI) observations reveal a tidal tail from NGC 275 which extends many optical radii beyond the interacting pair. The HI morphology implies a prograde encounter between the galaxy pair approximately 1.5 x 10**8 years ago. The Halpha emission from NGC 275 indicates clumpy irregular star-formation, clumpiness which is mirrored by the underlying mass distribution as traced by the Ks-band emission. The molecular gas distribution is striking in its anti-correlation with the {HII regions. Despite the evolved nature of NGC 275's interaction and its barred potential, neither the molecular gas nor the star formation are centrally concentrated. We suggest that this structure results from stochastic star formation leading to preferential consumption of the gas in certain regions of the galaxy. In contrast to the often assumed picture of interacting galaxies, NGC 275, which appears to be close to merger, does not display enhanced or centrally concentrated star formation. If the eventual merger is to lead to a significant burst of star formation it must be preceded by a significant conversion of atomic to molecular gas as at the current rate of star formation all the molecular gas will be exhausted by the time the merger is complete.

We present multi-frequency radio observations of a sample of z~2 obscured (type-2) quasars in the Spitzer extragalactic First Look Survey area. We combine the public data at 1.4 GHz, used in the selection of these sources, with new observations at 610 MHz (GMRT) and at 4.9 GHz (VLA). We find the sample includes sources with steep, flat and gigahertz-peaked spectra. There are no strong correlations between the presence or absence of emission lines in the optical spectra and the radio spectral properties of the sample. However, there are no secure flat-spectrum type-2 quasars with narrow emission lines which would be problematic for unified schemes. Most of the population have straight radio spectra with spectral index alpha~1 as is expected for developed, potentially FRI-like, jets in which continous injection of relativistic electrons is accompanied by inverse-Compton losses against the cosmic microwave background.

Using the third data release of the Sloan Digital Sky Survey (SDSS) we have rigorously defined a volume limited sample of early-type galaxies in the redshift range z < 0.1. We have defined the density of the local environment for each galaxy using a method which takes account of the redshift bias introduced by survey boundaries if traditional methods are used. At luminosities greater than our absolute r-band magnitude cutoff of -20.45 the mean density of environment shows no trend with redshift. We calculate the Lick indices for the entire sample and correct for aperture effects and velocity dispersion in a model independent way. Although we find no dependence of redshift or luminosity with environment we do find that the mean velocity dispersion, sigma, of early-type galaxies in dense environments tends to be higher than in low density environments. Taking account of this effect we find that several indices show small but very significant trends with environment that are not the result of the correlation between indices and velocity dispersion. The statistical significance of the data is sufficiently high to reveal that models accounting only for alpha-enhancement struggle to produce a consistent picture of age and metallicity of the sample galaxies, whereas a model that also includes carbon enhancement fares much better. We find that early-type galaxies in the field are younger than those in environments typical of clusters but that neither metallicity, alpha-enhancement nor carbon enhancement are influenced by the environment. The youngest early-type galaxies in both field and cluster environments are those with the lowest sigma. However, there is some evidence that the objects with the largest sigma are slightly younger, especially in denser environments.

We have combined up-to-date stellar population synthesis models, a simple radiative transfer approach, and a fully comprehensive dust model with the aim of developing a simple but quantitative way of interpreting the mid-infrared spectra of galaxies. We apply these models to the observed correlations of mid-infrared luminosities (at 8 and 24 micron) with the star-formation rate of normal galaxies and find that the observations are naturally reproduced by our models. We further find that the observed 24 micron correlation places a weak constraint on relative distribution of dust and stars.

The effect of galaxy interactions on star formation has been investigated using Data Release 1 of the Sloan Digital Sky Survey (SDSS). Both the imaging and spectroscopy data products have been used to construct a catalogue of nearest companions to a volume limited (0.03 < z < 0.1) sample of galaxies drawn from the main galaxy sample of SDSS. Of the 13973 galaxies in the volume limited sample, we have identified 12492 systems with companions at projected separations less than 300 kpc. Star-formation rates for the volume-limited sample have been calculated from extinction and aperture corrected H-alpha luminosities and, where available, IRAS data. Specific star formation rates were calculated by estimating galaxy masses from z-band luminosities, and r-band concentration indices were used as an indicator of morphological class. The mean specific star-formation rate is found to be strongly enhanced for projected separations of less than 30 kpc. For late-type galaxies the correlation extends out to projected separations of 300 kpc and is most pronounced in actively star-forming systems. The specific star-formation rate is observed to decrease with increasing recessional velocity difference, but the magnitude of this effect is small compared to that associated with the projected separation. We also observe a tight relationship between the concentration index and pair separation; the mean concentration index is largest for pairs with separations of approximately 75 kpc and declines rapidly for separations smaller than this. This is interpreted as being due to the presence of tidally-triggered nuclear starbursts in close pairs. Further, we find no dependence of star formation enhancement on the morphological type or mass of the companion galaxy

We present high-sensitivity radio continuum observations with the VLA and Ryle Telescope at 1.5, 4.9, 8.4 and 15.4 GHz of the dwarf irregular galaxy NGC 1569. The radio data show an extended, irregularly-shaped halo with filamentary structure around the galaxy. The spectral index maps reveal an unusually patchy distribution with regions of flat spectral index extending into the halo. The data allow us to perform a spatially-resolved spectral fitting analysis of the continuum emission from which we derive maps of the thermal and synchrotron emission. The thermal radio emission is concentrated towards the brightest \hii region west of the super star clusters A and B whereas the distribution of the synchrotron emission peaks in a bar-like structure in the disk extending between the two clusters. The total flux density of the thermal radio emission allows us to derive the integrated synchrotron spectrum and we confirm the break in the spectrum that was found by \citet{isr-deb88}. We discuss various possibilities that could produce such a break and conclude that the only mechanism able to fit the radio data and being consistent with data at other wavelengths is a convective wind allowing cosmic ray electrons to escape from the halo.

Both Arp 140 and Arp 104 reveal extended tidal tails in the HI. The Halpha and FIR fluxes of Arp 140 yield similar SFR of ~0.8 solar masses per year. In contrast the Halpha flux of Arp 104 yields a SFR of ~0.05 solar masses per year, ~ 20 times smaller than that obtained from the FIR flux. Spectra were used to examine the changing velocity of atomic and molecular gas in NGC 5218 (Arp 104). The atomic and molecular gas were found to be dynamically similar with comparable velocities and velocity widths across the galaxy; consistent with the two phases responding similarly to the interaction, or enhanced HI to CO conversion in the centre of the galaxy.

We have performed 3D hydrodynamical simulations of FR-II radio sources in beta-profile cooling-flow clusters. The effects of cooling of the cluster gas were incorporated into a modified version of the ZEUS-MP code. The simulations followed not only the active phase of the radio source, but also the long term behaviour for up to 2 Gyr after the jets of the radio source were switched off. We find as expected that the radio source has a significant effect on the cooling flow while it is active, however we also find that the effects of the radio source on the cluster are long-lived. A buoyancy driven convective flow is established as the remnants of the radio source rise through the cluster dragging material from the cluster core. Although the central Mpc of the cluster reverts to having a cooling flow, this asymmetric convective flow is able to remove the cool gas accumulating at the cluster core and indeed there is a net outflow persisting for timescales of about an order of magnitude longer than the time for which the source is active or longer. The convective flow may also provide a mechanism to enhance the metallicity of the cluster gas at large cluster radii.

The effect of the expansion of powerful FR-II radio sources into a cluster environment is discussed. The analysis considers both the thermal and temporal evolution of the ICM which has passed through the bow shock of the radio source and the effect of this swept-up gas on the dynamics of the radio source itself. The final state of the swept-up ICM is critcally dependent on the thermal conductivity of the gas. If the gas behind the bow shock expands adiabatically, and the source is expanding into a steeply falling atmosphere, then a narrow dense layer will form as the radio source lifts gas out of the cluster potential. This layer has a cooling time very much less than that of the gas just ahead of the radio source. This effect does not occur if the thermal conductivity of the gas is high, or if the cluster atmosphere is shallow. The swept-up gas also affects the dyamics of the radio source especially as it slows towards sub-sonic expansion. The preferential accumulation of the swept-up gas to the sides of the cocoon leads to the aspect ratio of the source increasing. Eventually the contact surface must become Rayleigh-Taylor unstable leading both to inflow of the swept-up ICM into the cavity created by the cocoon, but also substantial mixing of the cooler denser swept-up gas with the ambient ICM thereby creating a multi-phase ICM. The radio source is likely to have a marked effect on the cluster on timescales long compared to the age of the source.

We have used SCUBA to observe a complete sample of 104 galaxies selected at 60 microns from the IRAS BGS and we present here the 850 micron measurements. Fitting the 60,100 and 850 micron fluxes with a single temperature dust model gives the sample mean temperature T=36 K and beta = 1.3. We do not rule out the possibility of dust which is colder than this, if a 20 K component was present then our dust masses would increase by factor 1.5-3. We present the first measurements of the luminosity and dust mass functions, which were well fitted by Schechter functions (unlike those 60 microns). We have correlated many global galaxy properties with the submillimetre and find that there is a tendancy for less optically luminous galaxies to contain warmer dust and have greater star formation efficiencies (cf. Young 1999). The average gas-to-dust ratio for the sample is 581 +/- 43 (using both atomic and molecular hydrogen), significantly higher than the Galactic value of 160. We believe this discrepancy is due to a cold dust component at T < 20 K. There is a suprisingly tight correlation between dust mass and the mass of molecular hydrogen as estimated from CO measurements, with an intrinsic scatter of ~50%.

We present a model for the compression and heating of the ICM by powerful radio galaxies and quasars. Based on a self-similar model of the dynamical evolution of FRII-type objects we numerically integrate the hydrodynamic equations governing the flow of the shocked ICM in between the bow shock and the radio lobes of these sources. The resulting gas properties are presented and discussed. The X-ray emission of the shocked gas is calculated and is found to be in agreement with observations. The enhancement of the X-ray emission of cluster gas due to the presence of powerful radio galaxies may play an important role in the direct detection of cluster gas at high redshifts.

We present results of a numerical integration of the hydrodynamical equations governing the self-similar, two-dimensional gas flow behind the bow shock of an FRII radio source embedded in an IGM with a power law density profile. The model predicts pressure gradients within the cocoons consistent with modest backflow. For very steep external density profiles sources may well not expand in a self-similar fashion and in this case the model is not self-consistent. The assumption of ram pressure confinement of the cocoons perpendicular to the jet axis is found to overestimate the ratio of the pressure in front of the radio hot spots and that in the cocoons. Based on the properties of the gas between bow shock and cocoon we calculate the X-ray surface brightness of the flow. This emission is found to be a good tracer of the density distribution within the flow and varies significantly with the properties of the unshocked IGM. The cooling-time of the shocked IGM is found to be comparable to, or greater than, the Hubble time. The influence of a radio source on the evolution of its gaseous surroundings therefore extends well beyond the limited life time of the source itself. We compare our results with the X-ray map of Cygnus A and find some evidence for cold, dense gas clumps in the surroundings of this object. The extended X-ray emission observed around 3C 356 may also be caused by the bow shock of this radio source. We also present an empirical model for the X-ray emission of the shocked IGM due to thermal bremsstrahlung.

In a previous paper we presented measurements of the properties of jets and cores in a large sample of FRII radio galaxies with z<0.3. Here we test, by means of Monte Carlo simulations, the consistency of those data with models in which the prominences of cores and jets are determined by relativistic beaming. We conclude that relativistic beaming is needed to explain the relationships between core and jet prominences, and that speeds between 0.5 and 0.7c on kpc scales provide the best fits to the data.

We present an analytical model for the cosmological evolution of the FRII source population. Based on an earlier model for the intrinsic radio luminosity - linear size evolution of these objects, we construct theoretical source samples. The source distributions in the radio power - linear size plane of these samples are then compared with that of an observed flux-limited sample. We find that the source parameters determining the radio luminosity of FRII objects can not be independent of each other. The best-fitting models predict the jet power to be correlated either with the life time of the source or with the shape of the density distribution of the source environment. The latter case is consistent with the observed tendency of the most luminous radio sources at high redshift to be located in richer and more extended environments than their low redshift counterparts. We also find evidence for a class of FRII sources distinctly different from the main population. These sources are extremely old and/or are located in very underdense environments. The luminosity function of FRII sources resulting from the model is in good agreement with previous results for high luminosity sources. The apparent luminosity evolution of the radio luminosity function is not reproduced because of the high flux limit of the used comparison sample. The cosmological evolution of the median linear size of FRII sources is found to be mild.

Using a self-similar model for the expansion of cocoons surrounding the jets in powerful extragalactic radio sources (type FRII), we investigate the influence of the properties of the gas surrounding these objects on their evolution. The variation of their radio luminosity as a function of linear size for individual sources is determined in dependence of the density distribution of the IGM the source is embedded in. Based on these results, the cosmological evolution of the FRII population and hence that of the environments of these sources can be constrained. The environments are found to be denser and/or more extended at high redshifts. The bow shock surrounding the cocoon of FRII sources compresses and heats the IGM. A numerical integration of the hydrodynamic equations governing the gas flow between the bow shock and the cocoon boundary is presented. From this we determine the appearance of the large scale structure of radio galaxies in X-rays and the cooling times for the IGM heated by the passage of the bow shock. This has important implications for the cosmological evolution of the IGM.

We present new radio continuum data at 4 frequencies on the supermassive, peculiar galaxy NGC 1961. These observations allow us to separate the thermal and the nonthermal radio emission and to determine the nonthermal spectral index distribution. This spectral index distribution in the galactic disk is unusual: at the maxima of the radio emission the synchrotron spectrum is very steep, indicating aged cosmic ray electrons. Away from the maxima the spectrum is much flatter. The steep spectrum of the synchrotron emission at the maxima indicates that a strong decline of the star formation rate has taken place at these sites. The extended radio emission is a sign of recent cosmic ray acceleration, probably by recent star formation. We suggest that a violent event in the past, most likely a merger or a collision with an intergalactic gas cloud, has caused the various unusual features of the galaxy.

In previous papers we have discussed high-resolution observations of a large sample of powerful radio galaxies with z<0.3. Jets are detected in up to 80 per cent of the sample, and radio cores in nearly all the objects; in addition, we are able to resolve the hot spots in most sources. In this paper we present measurements of the radio properties of these components. The prominences of the jets detected do not appear to be a function of radio luminosity, providing the clearest evidence yet that the reported low detection rate of jets in radio galaxies has been an artefact of low-sensitivity observations. We find a positive correlation between the total source length and core prominence in the narrow-line radio galaxies. We have found evidence for a relationship between hot spot size and total source size, but few other significant relationships between hot spot properties and those of the jets or lobes. We compare our measurements to those of Bridle et al (1994), based on observations of a sample of quasars, and argue that the results are consistent with a modification of the unified model in which the broad-line radio galaxies are the low-luminosity counterparts of quasars, although the situation is complicated by contamination with low-excitation radio galaxies which appear to have radio properties different from the high-excitation objects. We discuss the classes of empirical model that can be fitted to the dataset.

This is the second in a series of papers presenting an analytical model for the evolution of FRII radio sources. In this paper we evaluate the expected radio emission from a radio source incorporating energy loss processes for the relativistic electrons. By combining these results with our earlier dynamical model we calculate evolutionary tracks through the Power-Linear size diagram. These tracks are in good agreement with the observed distribution of sources in this diagram. The effects of different forms for the evolution of the magnetic field in the cocoon, the redshift of the source, the environment of the source and the defining parameters of the jet are investigated. The evolutionary tracks are found to be insensitive to the assumed form of the magnetic field evolution. Some evidence against protons as a major constituent of the jet material is also found.

We present observations made with the VLA at 1.5 and 8.4 GHz of the nearby FRI radio galaxy 3C296. The most recent models of FRI radio galaxies suggest that substantial deceleration must take place in their jets, with strongly relativistic velocities on parsec scales giving place to at most mildly relativistic velocities on scales of tens of kiloparsecs. The region over which this deceleration takes place is therefore of considerable interest. By considering the side-to-side asymmetries of the jets of 3C296, we constrain the region of strong deceleration in the source. Our observations show evidence that the jets have slow edges surrounding faster central spines. We discuss the implications of our observations for models of the magnetic field structure in these objects.

We present high-sensitivity multi-frequency radio continuum observations of the starburst galaxy NGC~2146. We have fitted these data with a three-dimensional diffusion model. The model can describe the radio emission from the inner disk of NGC~2146 well, indicating that diffusion is the dominant mode of propagation in this region. Our results are indicating that NGC~2146 has recently undergone a starburst, the star forming activity being located in a central bar. The spatial variation of the radio emission and of the spectral index yield tight constraints on the diffusion coefficient $D_0$ and the energy dependence of the diffusion. Away from the central bar of the galaxy the radio emission becomes filamentary and the diffusion model was found to be a poor fit to the data in these regions; we attribute this to different transport processes being important in the halo of the galaxy.