The recent demonstration of a real-time direct imaging radio interferometry correlator represents a new capability in radio astronomy. However wide field imaging with this method is challenging since wide-field effects and array non-coplanarity degrade image quality if not compensated for. Here we present an alternative direct imaging correlation strategy using a Direct Fourier Transform (DFT), modelled as a linear operator facilitating a matrix multiplication between the DFT matrix and a vector of the electric fields from each antenna. This offers perfect correction for wide field and non-coplanarity effects. When implemented with data from the Long Wavelength Array (LWA), it offers comparable computational performance to previously demonstrated direct imaging techniques, despite having a theoretically higher floating point cost. It also has additional benefits, such as imaging sparse arrays and control over which sky co-ordinates are imaged, allowing variable pixel placement across an image. It is in practice a highly flexible and efficient method of direct radio imaging when implemented on suitable arrays. A functioning Electric Field Direct imaging architecture using the DFT is presented, alongside an exploration of techniques for wide-field imaging similar to those in visibility based imaging, and an explanation of why they do not fit well to imaging directly with the digitized electric field data. The DFT imaging method is demonstrated on real data from the LWA telescope, alongside a detailed performance analysis, as well as an exploration of its applicability to other arrays.

In radio interferometry imaging, the gridding procedure of convolving visibilities with a chosen gridding function is necessary to transform visibility values into uniformly sampled grid points. We propose here a parameterised family of "least-misfit gridding functions" which minimise an upper bound on the difference between the DFT and FFT dirty images for a given gridding support width and image cropping ratio. When compared with the widely used spheroidal function with similar parameters, these provide more than 100 times better alias suppression and RMS misfit reduction over the usable dirty map. We discuss how appropriate parameter selection and tabulation of these functions allow for a balance between accuracy, computational cost and storage size. Although it is possible to reduce the errors introduced in the gridding or degridding process to the level of machine precision, accuracy comparable to that achieved by CASA requires only a lookup table with 300 entries and a support width of 3, allowing for a greatly reduced computation cost for a given performance.

Methods currently in use for locating and characterising sources in radio interferometry maps are designed for processing images, and require interferometric maps to be preprocessed so as to resemble conventional images. We demonstrate a Bayesian code - BaSC - that takes into account the interferometric visibility data despite working with more computationally manageable image domain data products. This method is better able to discriminate nearby sources than the commonly used SExtractor, and has potential even in more complicated cases. We also demonstrate the correctness of the Bayesian resolving formula for simulated data, and its implications for source discrimination at distances below the full width half maximum of the restoring beam.

A main science goal for the Large Synoptic Survey Telescope (LSST) is to measure the cosmic shear signal from weak lensing to extreme accuracy. One difficulty, however, is that with the short exposure time ($\simeq$15 seconds) proposed, the spatial variation of the Point Spread Function (PSF) shapes may be dominated by the atmosphere, in addition to optics errors. While optics errors mainly cause the PSF to vary on angular scales similar or larger than a single CCD sensor, the atmosphere generates stochastic structures on a wide range of angular scales. It thus becomes a challenge to infer the multi-scale, complex atmospheric PSF patterns by interpolating the sparsely sampled stars in the field. In this paper we present a new method, PSFent, for interpolating the PSF shape parameters, based on reconstructing underlying shape parameter maps with a multi-scale maximum entropy algorithm. We demonstrate, using images from the LSST Photon Simulator, the performance of our approach relative to a 5th-order polynomial fit (representing the current standard) and a simple boxcar smoothing technique. Quantitatively, PSFent predicts more accurate PSF models in all scenarios and the residual PSF errors are spatially less correlated. This improvement in PSF interpolation leads to a factor of 3.5 lower systematic errors in the shear power spectrum on scales smaller than $\sim13'$, compared to polynomial fitting. We estimate that with PSFent and for stellar densities greater than $\simeq1/{\rm arcmin}^{2}$, the spurious shear correlation from PSF interpolation, after combining a complete 10-year dataset from LSST, is lower than the corresponding statistical uncertainties on the cosmic shear power spectrum, even under a conservative scenario.

We present a method for accelerating the calculation of CMB power spectra, matter power spectra and likelihood functions for use in cosmological parameter estimation. The algorithm, called CosmoNet, is based on training a multilayer perceptron neural network and shares all the advantages of the recently released Pico algorithm of Fendt & Wandelt, but has several additional benefits in terms of simplicity, computational speed, memory requirements and ease of training. We demonstrate the capabilities of CosmoNet by computing CMB power spectra over a box in the parameter space of flat \Lambda CDM models containing the 3\sigma WMAP1 confidence region. We also use CosmoNet to compute the WMAP3 likelihood for flat \Lambda CDM models and show that marginalised posteriors on parameters derived are very similar to those obtained using CAMB and the WMAP3 code. We find that the average error in the power spectra is typically 2-3% of cosmic variance, and that CosmoNet is \sim 7 \times 10^4 faster than CAMB (for flat models) and \sim 6 \times 10^6 times faster than the official WMAP3 likelihood code. CosmoNet and an interface to CosmoMC are publically available at www.mrao.cam.ac.uk/software/cosmonet.

A new gauge theory of gravity is presented. The theory is constructed in a flat background spacetime and employs gauge fields to ensure that all relations between physical quantities are independent of the positions and orientations of the matter fields. In this manner all properties of the background spacetime are removed from physics, and what remains are a set of `intrinsic' relations between physical fields. The properties of the gravitational gauge fields are derived from both classical and quantum viewpoints. Field equations are then derived from an action principle, and consistency with the minimal coupling procedure selects an action that is unique up to the possible inclusion of a cosmological constant. This in turn singles out a unique form of spin-torsion interaction. A new method for solving the field equations is outlined and applied to the case of a time-dependent, spherically-symmetric perfect fluid. A gauge is found which reduces the physics to a set of essentially Newtonian equations. These equations are then applied to the study of cosmology, and to the formation and properties of black holes. The existence of global solutions enables one to discuss the properties of field lines inside the horizon due to a point charge held outside it. The Dirac equation is studied in a black hole background and provides a quick derivation of the Hawking temperature.

We develop the maximum-entropy weak shear mass reconstruction method presented in earlier papers by taking each background galaxy image shape as an independent estimator of the reduced shear field and incorporating an intrinsic smoothness into the reconstruction. The characteristic length scale of this smoothing is determined by Bayesian methods. Within this algorithm the uncertainties due to both the intrinsic distribution of galaxy shapes and galaxy shape estimation are carried through to the final mass reconstruction, and the mass within arbitrarily shaped apertures can be calculated with corresponding uncertainties. We apply this method to two clusters taken from n-body simulations using mock observations corresponding to Keck LRIS and mosaiced HST WFPC2 fields. We demonstrate that the Bayesian choice of smoothing length is sensible and that masses within apertures (including one on a filamentary structure) are reliable, provided the field of view is not too small. We apply the method to data taken on the cluster MS1054-03 using the Keck LRIS (Clowe et al. 2000) and HST (Hoekstra et al. 2000), finding results in agreement with this previous work; we also present reconstructions with optimal smoothing lengths, and mass estimates which do not rely on any assumptions of circular symmetry. The code used in this work (LensEnt2) is available from the web.