Search

Around a million quasars have been catalogued in the Universe by probing deeper and using new methods for discovery. However, the hardest ones to find seem to be the rarest and brightest specimen. In this work, we study the properties of the most luminous of all quasars found so far. It has been overlooked until recently, which demonstrates that modern all-sky surveys have much to reveal. The black hole in this quasar accretes around one solar mass per day onto an existing mass of $\sim$17 billion solar masses. In this process its accretion disc alone releases a radiative energy of $2\times 10^{41}$ Watts. If the quasar is not strongly gravitationally lensed, then its broad line region (BLR) is expected to have the largest physical and angular diameter occurring in the Universe, and will allow the Very Large Telescope Interferometer to image its rotation and measure its black hole mass directly. This will be an important test for BLR size-luminosity relations, whose extrapolation has underpinned common black-hole mass estimates at high redshift., Comment: This version of the article has been accepted for publication after peer review and is subject to Springer Nature's AM terms of use, but is not the Version of Record and does not reflect post-acceptance improvements, or any corrections. The Version of Record is available online at http://dx.doi.org/10.1038/s41550-024-02195-x

Precise instrumental calibration is of crucial importance to 21-cm cosmology experiments. The Murchison Widefield Array's (MWA) Phase II compact configuration offers us opportunities for both redundant calibration and sky-based calibration algorithms; using the two in tandem is a potential approach to mitigate calibration errors caused by inaccurate sky models. The MWA Epoch of Reionization (EoR) experiment targets three patches of the sky (dubbed EoR0, EoR1, and EoR2) with deep observations. Previous work in \cite{Li_2018} and \cite{Wenyang_2019} studied the effect of tandem calibration on the EoR0 field and found that it yielded no significant improvement in the power spectrum over sky-based calibration alone. In this work, we apply similar techniques to the EoR1 field and find a distinct result: the improvements in the power spectrum from tandem calibration are significant. To understand this result, we analyze both the calibration solutions themselves and the effects on the power spectrum over three nights of EoR1 observations. We conclude that the presence of the bright radio galaxy Fornax A in EoR1 degrades the performance of sky-based calibration, which in turn enables redundant calibration to have a larger impact. These results suggest that redundant calibration can indeed mitigate some level of model-incompleteness error., Comment: 20 pages, 11 figures, accepted for publication in PASA

We use the latest multi-redshift ($z=6.5-8.7$) upper limits on the 21-cm signal from the Murchison Widefield Array (MWA) to explore astrophysical models which are inconsistent with the data. These upper limits are achieved using 298 h of carefully excised data over four observing seasons. To explore these upper limits in the context of reionisation astrophysics, we use 21CMMC. We then connect the disfavoured regions of parameter space to existing observational constraints on reionisation such as high-$z$ galaxy ultra-violet (UV) luminosity functions, background UV photoionisation rate, intergalactic medium (IGM) neutral fraction, the electron scattering optical depth and the soft-band X-ray emissivity. We find the vast majority of models disfavoured by the MWA limits are already inconsistent with existing observational constraints. These inconsistent models arise from two classes of models: (i) `cold' reionisation and (ii) pure matter density fluctuations (i.e. no reionisation). However, a small subsample of models are consistent implying the existing MWA limits provide unique information in disfavouring models of reionisation, albeit extremely weakly. We also provide the first limits on the soft-band X-ray emissivity from galaxies at high redshifts, finding $1\sigma$ lower limits of $\epsilon_{{\rm X},0.5-2~{\rm keV}}\gtrsim10^{34.5}$ erg s$^{-1}$ Mpc$^{-3}$. Finally, we recover 95 per cent disfavoured limits on the IGM spin temperature of $\bar{T}_{\rm S}\lesssim$ 1.3, 1.4, 1.5, 1.8, 2.1, 2.4 K at $z=6.5, 6.8, 7.1, 7.8, 8.2, 8.7$. With this we infer the IGM must have undergone, at the very least, a small amount of X-ray heating. Note, the limits on $\epsilon_{{\rm X},0.5-2~{\rm keV}}$ and $\bar{T}_{\rm S}$ are conditional on the IGM neutral fraction., Comment: 14 pages, 7 figures, 1 table. Submitted to MNRAS, comments welcome. arXiv admin note: text overlap with arXiv:2006.03203

Since the discovery of active galactic nuclei (AGN) and their subclasses, a unification scheme of AGN has been long sought. Orientation-based unified models predict that some of the diversity within AGN subclasses can be explained by the different viewing angles of the observer. Several orientation categorisations have been suggested, but a widely applicable measure has yet to be found. Using the properties of the ultraviolet and optical broad emission lines of quasars, in particular the velocity offsets and line widths of high-ionisation CIV and low-ionisation MgII lines, a correlation has been measured. It is postulated that this correlation is due to the viewing angle of the observer. Comparison with other orientation tracers shows consistency with this interpretation. Using a simulation of a wide angle disk-wind model for the broad emission line region, we successfully replicate the observed correlation with inclination. Future more detailed modelling will not only enable improved accuracy in the determination of the viewing angle to individual AGN, but will also substantially increase our understanding of the emitting regions of AGN., Comment: Accepted for publication in MNRAS. 15 pages, 17 figures, 2 table

We present new HST WFPC3 imaging of four gravitationally lensed quasars: MG 0414+0534; RXJ 0911+0551; B 1422+231; WFI J2026-4536. In three of these systems we detect wavelength-dependent microlensing, which we use to place constraints on the sizes and temperature profiles of the accretion discs in each quasar. Accretion disc radius is assumed to vary with wavelength according to the power-law relationship $r\propto \lambda^p$, equivalent to a radial temperature profile of $T\propto r^{-1/p}$. The goal of this work is to search for deviations from standard thin disc theory, which predicts that radius goes as wavelength to the power $p=4/3$. We find a wide range of power-law indices, from $p=1.4^{+0.5}_{-0.4}$ in B 1422+231 to $p=2.3^{+0.5}_{-0.4}$ in WFI J2026-4536. The measured value of $p$ appears to correlate with the strength of the wavelength-dependent microlensing. We explore this issue with mock simulations using a fixed accretion disc with $p=1.5$, and find that cases where wavelength-dependent microlensing is small tend to under-estimate the value of $p$. This casts doubt on previous ensemble single-epoch measurements which have favoured low values using samples of lensed quasars that display only moderate chromatic effects. Using only our systems with strong chromatic microlensing we prefer $p>4/3$, corresponding to shallower temperature profiles than expected from standard thin disc theory., Comment: 20 pages, 11 figures, accepted for publication in MNRAS

A key to understanding quasar unification paradigms is the emission properties of broad absorption line quasars (BALQs). The fact that only a small fraction of quasar spectra exhibit deep absorption troughs blueward of the broad permitted emission lines provides a crucial clue to the structure of quasar emitting regions. To learn whether it is possible to discriminate between the BALQ and non-BALQ populations given the observed spectral properties of a quasar, we employ two approaches: one based on statistical methods and the other supervised machine learning classification, applied to quasar samples from the Sloan Digital Sky Survey. The features explored include continuum and emission line properties, in particular the absolute magnitude, redshift, spectral index, line width, asymmetry, strength, and relative velocity offsets of high-ionisation CIV $\lambda 1549$ and low-ionisation MgII $\lambda 2798$ lines. We consider a complete population of quasars, and assume that the statistical distributions of properties represent all angles where the quasar is viewed without obscuration. The distributions of the BALQ and non-BALQ sample properties show few significant differences. None of the observed continuum and emission line features are capable of differentiating between the two samples. Most published narrow disk-wind models are inconsistent with these observations, and an alternative disk-wind model is proposed. The key feature of the proposed model is a disk-wind filling a wide opening angle with multiple radial streams of dense clumps., Comment: Accepted for publication in MNRAS. 19 pages, 3 figures, 6 table

The structure and kinematics of the broad line region (BLR) in quasars are still not well established. One popular BLR model is the disk-wind model that offers a geometric unification of a quasar based on the angle of viewing. We construct a simple kinematical disk-wind model with a narrow outflowing wind angle. The model is combined with radiative transfer in the Sobolev, or high velocity, limit. We examine how angle of viewing affects the observed characteristics of the emission line, especially the line widths and velocity offsets. The line profiles exhibit distinct properties depending on the orientation, wind opening angle, and region of the wind where the emission arises. At low inclination angle (close to face-on), we find the shape of the emission line is asymmetric with narrow width and significantly blueshifted. As the inclination angle increases (close to edge-on), the line profile becomes more symmetric, broader, and less blueshifted. Additionally, lines that arise close to the base of the disk wind, near the accretion disk, tend to be broad and symmetric. The relative increase in blueshift of the emission line with increasing wind vertical distance is larger for polar winds compared with equatorial winds. By considering the optical thickness of the wind, single-peaked line profiles are recovered for the intermediate and equatorial outflowing wind. The model is able to reproduce a faster response in either the red and blue sides of the line profile found in reverberation mapping studies. A quicker response in the red side is achieved in the model with a polar wind and intermediate wind opening angle at low viewing angle. The blue side response is faster for an equatorial wind seen at high inclination., Comment: Accepted for publication in PASA. 15 pages, 6 figures, 1 table

Low radio frequency solar observations using the Murchison Widefield Array have recently revealed the presence of numerous weak, short-lived and narrow-band emission features, even during moderately quiet solar conditions. These non-thermal features occur at rates of many thousands per hour in the 30.72 MHz observing bandwidth, and hence, necessarily require an automated approach for their detection and characterization. Here, we employ continuous wavelet transform using a mother Ricker wavelet for feature detection from the dynamic spectrum. We establish the efficacy of this approach and present the first statistically robust characterization of the properties of these features. In particular, we examine distributions of their peak flux densities, spectral spans, temporal spans and peak frequencies. We can reliably detect features weaker than 1 SFU, making them, to the best of our knowledge, the weakest bursts reported in literature. The distribution of their peak flux densities follows a power law with an index of -2.23 in the 12-155 SFU range, implying that they can provide an energetically significant contribution to coronal and chromospheric heating. These features typically last for 1-2 seconds and possess bandwidths of about 4-5 MHz. Their occurrence rate remains fairly flat in the 140-210 MHz frequency range. At the time resolution of the data, they appear as stationary bursts, exhibiting no perceptible frequency drift. These features also appear to ride on a broadband background continuum, hinting at the likelihood of them being weak type-I bursts., Comment: 15 pages, 14 figures, Accepted for publication in ApJ

Using the Murchison Widefield Array (MWA), the low-frequency Square Kilometre Array (SKA1 LOW) precursor located in Western Australia, we have completed the GaLactic and Extragalactic All-sky MWA (GLEAM) survey, and present the resulting extragalactic catalogue, utilising the first year of observations. The catalogue covers 24,831 square degrees, over declinations south of $+30^\circ$ and Galactic latitudes outside $10^\circ$ of the Galactic plane, excluding some areas such as the Magellanic Clouds. It contains 307,455 radio sources with 20 separate flux density measurements across 72--231MHz, selected from a time- and frequency- integrated image centred at 200MHz, with a resolution of $\approx 2$'. Over the catalogued region, we estimate that the catalogue is 90% complete at 170mJy, and 50% complete at 55mJy, and large areas are complete at even lower flux density levels. Its reliability is 99.97% above the detection threshold of $5\sigma$, which itself is typically 50mJy. These observations constitute the widest fractional bandwidth and largest sky area survey at radio frequencies to date, and calibrate the low frequency flux density scale of the southern sky to better than 10%. This paper presents details of the flagging, imaging, mosaicking, and source extraction/characterisation, as well as estimates of the completeness and reliability. All source measurements and images are available online (http://www.mwatelescope.org/science/gleam-survey). This is the first in a series of publications describing the GLEAM survey results., Comment: 30 pages, 18 figures, 6 tables, published in Monthly Notices of the Royal Astronomical Society

Black hole mass is a key factor in determining how a black hole interacts with its environment. However, the determination of black hole masses at high redshifts depends on secondary mass estimators, which are based on empirical relationships and broad approximations. A dynamical disk wind broad line region (BLR) model of active galactic nuclei (AGN) is built in order to test the impact on the black hole mass calculation due to different BLR geometries and the inclination of the AGN. Monte Carlo simulations of two disk wind models are constructed to recover the virial scale factor, $f$, at various inclination angles. The resulting $f$ values strongly correlate with inclination angle, with large $f$ values associated with small inclination angles (close to face-on) and small $f$ values with large inclination angles (close to edge-on). The $f$ factors are consistent with previously determined $f$ values, found from empirical relationships. Setting $f$ as a constant may introduce a bias into virial black hole mass estimates for a large sample of AGN. However, the extent of the bias depends on the line width characterisation (e.g. full width at half maximum (FWHM) or line dispersion). Masses estimated using $f_{\text{FWHM}}$ tend to biased towards larger masses, but this can be corrected by calibrating for the width or shape of the emission line., Comment: Accepted for publication in PASA. 9 pages, 5 figures, 2 table

We have investigated, using both a theoretical and an empirical approach, the frequency of low redshift galaxy-galaxy lensing systems in which the signature of weak lensing might be directly detectable. We find good agreement between these two approaches. In order to make a theoretical estimate of the weak lensing shear, $\gamma$, for each galaxy in a catalogue, we have made an estimate of the asymptotic circular velocity from the stellar mass using three different approaches: from a simulation based relation, from an empirically-derived relation, and using the baryonic Tully-Fisher relation. Using data from the Galaxy and Mass Assembly redshift survey we estimate the frequency of detectable weak lensing at low redshift. We find that to a redshift of $z\sim 0.6$, the probability of a galaxy being weakly lensed by at least $\gamma = 0.02$ is $\sim 0.01$. A scatter in the $M_*-M_h$ relation results in a shift towards higher measured shears for a given population of galaxies. Given this, and the good probability of weak lensing at low redshifts, we have investigated the feasibility of measuring the scatter in the $M_*-M_h$ relation using shear statistics. This is a novel measurement, and is made possible because DSM is able to make individual \itshape direct~\upshape shear measurements, in contrast to traditional weak lensing techniques which can only make statistical measurements. We estimate that for a shear measurement error of $\Delta\gamma = 0.02$ (consistent with the sensitivity of DSM), a sample of $\sim$50,000 spatially and spectrally resolved galaxies would allow a measurement of the scatter in the $M_*-M_h$ relation to be made. While there are no currently existing IFU surveys of this size, there are upcoming surveys which will provide this data (e.g The Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), surveys with Hector, and the Square Kilometre Array (SKA))., Comment: Accepted for publication in PASA. 18 pages, 12 figures, 1 table

We have developed a new technique called Direct Shear Mapping (DSM) to measure gravitational lensing shear directly from observations of a single background source. The technique assumes the velocity map of an un-lensed, stably-rotating galaxy will be rotationally symmetric. Lensing distorts the velocity map making it asymmetric. The degree of lensing can be inferred by determining the transformation required to restore axisymmetry. This technique is in contrast to traditional weak lensing methods, which require averaging an ensemble of background galaxy ellipticity measurements, to obtain a single shear measurement. We have tested the efficacy of our fitting algorithm with a suite of systematic tests on simulated data. We demonstrate that we are in principle able to measure shears as small as 0.01. In practice, we have fitted for the shear in very low redshift (and hence un-lensed) velocity maps, and have obtained null result with an error of $\pm 0.01$. This high sensitivity results from analysing spatially resolved spectroscopic images (i.e. 3D data cubes), including not just shape information (as in traditional weak lensing measurements) but velocity information as well. Spirals and rotating ellipticals are ideal targets for this new technique. Data from any large IFU or radio telescope is suitable, or indeed any instrument with spatially resolved spectroscopy such as SAMI, ALMA, HETDEX and SKA., Comment: 14 pages, 12 figures, 2 tables

We present new integral field spectroscopy of the gravitationally lensed broad absorption line (BAL) quasar H1413+117, covering the ultraviolet to visible rest-frame spectral range. We observe strong microlensing signatures in lensed image D, and we use this microlensing to simultaneously constrain both the broad emission and broad absorption line gas. By modeling the lens system over the range of probable lensing galaxy redshifts and using on a new argument based on the wavelength-independence of the broad line lensing magnifications, we determine that there is no significant broad line emission from smaller than ~20 light days. We also perform spectral decomposition to derive the intrinsic broad emission line (BEL) and continuum spectrum, subject to BAL absorption. We also reconstruct the intrinsic BAL absorption profile, whose features allow us to constrain outflow kinematics in the context of a disk-wind model. We find a very sharp, blueshifted onset of absorption of 1,500 km/s in both C IV and N V that may correspond to an inner edge of a disk-wind's radial outflow. The lower ionization Si IV and Al III have higher-velocity absorption onsets, consistent with a decreasing ionization parameter with radius in an accelerating outflow. There is evidence of strong absorption in the BEL component which indicates a high covering factor for absorption over two orders of magnitude in outflow radius., Comment: 29 pages, 8 figures

An FPGA-based digital-receiver has been developed for a low-frequency imaging radio interferometer, the Murchison Widefield Array (MWA). The MWA, located at the Murchison Radio-astronomy Observatory (MRO) in Western Australia, consists of 128 dual-polarized aperture-array elements (tiles) operating between 80 and 300\,MHz, with a total processed bandwidth of 30.72 MHz for each polarization. Radio-frequency signals from the tiles are amplified and band limited using analog signal conditioning units; sampled and channelized by digital-receivers. The signals from eight tiles are processed by a single digital-receiver, thus requiring 16 digital-receivers for the MWA. The main function of the digital-receivers is to digitize the broad-band signals from each tile, channelize them to form the sky-band, and transport it through optical fibers to a centrally located correlator for further processing. The digital-receiver firmware also implements functions to measure the signal power, perform power equalization across the band, detect interference-like events, and invoke diagnostic modes. The digital-receiver is controlled by high-level programs running on a single-board-computer. This paper presents the digital-receiver design, implementation, current status, and plans for future enhancements., Comment: 14 pages, 7 figures

The detection of the Epoch of Reionization (EoR) in the redshifted 21-cm line is a challenging task. Here we formulate the detection of the EoR signal using the drift scan strategy. This method potentially has better instrumental stability as compared to the case where a single patch of sky is tracked. We demonstrate that the correlation time between measured visibilities could extend up to 1-2 hr for an interferometer array such as the Murchison Widefield Array (MWA), which has a wide primary beam. We estimate the EoR power based on cross-correlation of visibilities across time and show that the drift scan strategy is capable of the detection of the EoR signal with comparable/better signal-to-noise as compared to the tracking case. We also estimate the visibility correlation for a set of bright point sources and argue that the statistical inhomogeneity of bright point sources might allow their separation from the EoR signal., Comment: 23 pages, 14 figures, accepted for publication in ApJ

In this paper we explore for the first time the relative magnitudes of three fundamental sources of uncertainty, namely, foreground contamination, thermal noise and sample variance in detecting the HI power spectrum from the Epoch of Reionization (EoR). We derive limits on the sensitivity of a Fourier synthesis telescope to detect EoR based on its array configuration and a statistical representation of images made by the instrument. We use the Murchison Widefield Array (MWA) configuration for our studies. Using a unified framework for estimating signal and noise components in the HI power spectrum, we derive an expression for and estimate the contamination from extragalactic point-like sources in three-dimensional k-space. Sensitivity for EoR HI power spectrum detection is estimated for different observing modes with MWA. With 1000 hours of observing on a single field using the 128-tile MWA, EoR detection is feasible (S/N > 1 for $k\lesssim 0.8$ Mpc$^{-1}$). Bandpass shaping and refinements to the EoR window are found to be effective in containing foreground contamination, which makes the instrument tolerant to imaging errors. We find that for a given observing time, observing many independent fields of view does not offer an advantage over a single field observation when thermal noise dominates over other uncertainties in the derived power spectrum., Comment: 19 pages, 14 figures (19 figures in total including subfigures), 2 appendices, in press for publication in ApJ

We present techniques for bridging the gap between idealized inverse covariance weighted quadratic estimation of 21 cm power spectra and the real-world challenges presented universally by interferometric observation. By carefully evaluating various estimators and adapting our techniques for large but incomplete data sets, we develop a robust power spectrum estimation framework that preserves the so-called "EoR window" and keeps track of estimator errors and covariances. We apply our method to observations from the 32-tile prototype of the Murchinson Widefield Array to demonstrate the importance of a judicious analysis technique. Lastly, we apply our method to investigate the dependence of the clean EoR window on frequency--especially the frequency dependence of the so-called "wedge" feature--and establish upper limits on the power spectrum from z = 6.2 to z = 11.7. Our lowest limit is Delta(k) < 0.3 Kelvin at 95% confidence at a comoving scale k = 0.046 Mpc^-1 and z = 9.5., Comment: 27 pages, 10 figures. Updated to match the accepted Phys. Rev. D version

Significant new opportunities for astrophysics and cosmology have been identified at low radio frequencies. The Murchison Widefield Array is the first telescope in the Southern Hemisphere designed specifically to explore the low-frequency astronomical sky between 80 and 300 MHz with arcminute angular resolution and high survey efficiency. The telescope will enable new advances along four key science themes, including searching for redshifted 21 cm emission from the epoch of reionisation in the early Universe; Galactic and extragalactic all-sky southern hemisphere surveys; time-domain astrophysics; and solar, heliospheric, and ionospheric science and space weather. The Murchison Widefield Array is located in Western Australia at the site of the planned Square Kilometre Array (SKA) low-band telescope and is the only low-frequency SKA precursor facility. In this paper, we review the performance properties of the Murchison Widefield Array and describe its primary scientific objectives., Comment: Updated with revisions. 32 pages including figures and references. Submitted to PASA

We demonstrate a new technique for determining the physical conditions of the broad line emitting gas in quasars, using near-infrared hydrogen emission lines. Unlike higher ionisation species, hydrogen is an efficient line emitter for a very wide range of photoionisation conditions, and the observed line ratios depend strongly on the density and photoionisation state of the gas present. A locally optimally emitting cloud model of the broad emission line region was compared to measured emission lines of four nearby ($z\approx0.2$) quasars that have optical and NIR spectra of sufficient signal-to-noise to measure their Paschen lines. The model provides a good fit to three of the objects, and a fair fit to the fourth object, a ULIRG. We find that low incident ionising fluxes ($\phih<10^{18}$\cmsqs), and high gas densities ($\nh>10^{12}$\cmcu) are required to reproduce the observed hydrogen emission line ratios. This analysis demonstrates that the use of composite spectra in photoionisation modelling is inappropriate; models must be fitted to the individual spectra of quasars., Comment: 12 pages, accepted for publication by ApJ

The Murchison Widefield Array (MWA) is a new low-frequency, wide field-of-view radio interferometer under development at the Murchison Radio-astronomy Observatory (MRO) in Western Australia. We have used a 32-element MWA prototype interferometer (MWA-32T) to observe two 50-degree diameter fields in the southern sky in the 110 MHz to 200 MHz band in order to evaluate the performance of the MWA-32T, to develop techniques for epoch of reionization experiments, and to make measurements of astronomical foregrounds. We developed a calibration and imaging pipeline for the MWA-32T, and used it to produce ~15' angular resolution maps of the two fields. We perform a blind source extraction using these confusion-limited images, and detect 655 sources at high significance with an additional 871 lower significance source candidates. We compare these sources with existing low-frequency radio surveys in order to assess the MWA-32T system performance, wide field analysis algorithms, and catalog quality. Our source catalog is found to agree well with existing low-frequency surveys in these regions of the sky and with statistical distributions of point sources derived from Northern Hemisphere surveys; it represents one of the deepest surveys to date of this sky field in the 110 MHz to 200 MHz band., Comment: 20 pages, 6 tables, 12 figures. 1 online-only machine readable table. Submitted to ApJ

We present the first spectroscopic images of solar radio transients from the prototype for the Murchison Widefield Array (MWA), observed on 2010 March 27. Our observations span the instantaneous frequency band 170.9-201.6 MHz. Though our observing period is characterized as a period of `low' to `medium' activity, one broadband emission feature and numerous short-lived, narrowband, non-thermal emission features are evident. Our data represent a significant advance in low radio frequency solar imaging, enabling us to follow the spatial, spectral, and temporal evolution of events simultaneously and in unprecedented detail. The rich variety of features seen here reaffirms the coronal diagnostic capability of low radio frequency emission and provides an early glimpse of the nature of radio observations that will become available as the next generation of low frequency radio interferometers come on-line over the next few years., Comment: 15 pages, 5 figures, accepted for publication in Astrophysical Journal Letters. Movies for figures 4 and 5 available at http://www.mwatelescope.org/info/mwa_proto.html

The Murchison Widefield Array (MWA) is a dipole-based aperture array synthesis telescope designed to operate in the 80-300 MHz frequency range. It is capable of a wide range of science investigations, but is initially focused on three key science projects. These are detection and characterization of 3-dimensional brightness temperature fluctuations in the 21cm line of neutral hydrogen during the Epoch of Reionization (EoR) at redshifts from 6 to 10, solar imaging and remote sensing of the inner heliosphere via propagation effects on signals from distant background sources,and high-sensitivity exploration of the variable radio sky. The array design features 8192 dual-polarization broad-band active dipoles, arranged into 512 tiles comprising 16 dipoles each. The tiles are quasi-randomly distributed over an aperture 1.5km in diameter, with a small number of outliers extending to 3km. All tile-tile baselines are correlated in custom FPGA-based hardware, yielding a Nyquist-sampled instantaneous monochromatic uv coverage and unprecedented point spread function (PSF) quality. The correlated data are calibrated in real time using novel position-dependent self-calibration algorithms. The array is located in the Murchison region of outback Western Australia. This region is characterized by extremely low population density and a superbly radio-quiet environment,allowing full exploitation of the instrumental capabilities., Comment: 9 pages, 5 figures, 1 table. Accepted for publication in Proceedings of the IEEE

Several gravitationally lensed quasars are observed with anomalous magnifications in pairs of images that straddle a critical curve. Simple theoretical arguments suggest that the magnification of these images should be approximately equivalent, whereas one image is observed to be significantly demagnified. Microlensing provides a possible explanation for this discrepancy. There are two key parameters when modelling this effect. The first, the fraction of smooth matter in the lens at the image positions, has been explored by Schechter and Wambsganss (2002). They have shown that the anomalous flux ratio observed in the lensed quasar MG 0414+0534 is a priori a factor of 5 more likely if the assumed smooth matter content in the lens model is increased from 0% to 93%. The second parameter, the size of the emission region, is explored in this paper, and shown to be more significant. We find that the broadening of the magnification probability distributions due to smooth matter content is washed out for source sizes that are predicted by standard models for quasars. We apply our model to the anomalous lensed quasar MG 0414+0534, and find a 95% upper limit of 2.62 x 10^(16) h^(-1/2) (M/Msun)^(1/2) cm on the radius of the I-band emission region. The smooth matter percentage in the lens is unconstrained., Comment: 6 pages, 6 figures. To be published in MNRAS

We report on the detection of giant pulses from the Crab Nebula pulsar at a frequency of 200 MHz using the field deployment system designed for the Mileura Widefield Array's Low Frequency Demonstrator (MWA-LFD). Our observations are among the first high-quality detections at such low frequencies. The measured pulse shapes are deconvolved for interstellar pulse broadening, yielding a pulse-broadening time of 670$\pm$100 $\mu$s, and the implied strength of scattering (scattering measure) is the lowest that is estimated towards the Crab nebula from observations made so far. The sensitivity of the system is largely dictated by the sky background, and our simple equipment is capable of detecting pulses that are brighter than $\sim$9 kJy in amplitude. The brightest giant pulse detected in our data has a peak amplitude of $\sim$50 kJy, and the implied brightness temperature is $10^{31.6}$ K. We discuss the giant pulse detection prospects with the full MWA-LFD system. With a sensitivity over two orders of magnitude larger than the prototype equipment, the full system will be capable of detecting such bright giant pulses out to a wide range of Galactic distances; from $\sim$8 to $\sim$30 kpc depending on the frequency. The MWA-LFD will thus be a highly promising instrument for the studies of giant pulses and other fast radio transients at low frequencies., Comment: 10 pages, 6 figures, Accepted for publication in the Astrophysical Journal

Experiments were performed with prototype antenna tiles for the Mileura Widefield Array--Low Frequency Demonstrator (MWA-LFD) to better understand the widefield, wideband properties of their design and to characterize the radio frequency interference (RFI) between 80 and 300 MHz at the site in Western Australia. Observations acquired during the six month deployment confirmed the predicted sensitivity of the antennas, sky-noise dominated system temperatures, and phase-coherent interferometric measurements. The radio spectrum is remarkably free of strong terrestrial signals, with the exception of two narrow frequency bands allocated to satellite downlinks and rare bursts due to ground-based transmissions being scattered from aircraft and meteor trails. Results indicate the potential of the MWA-LFD to make significant achievements in its three key science objectives: epoch of reionziation science, heliospheric science, and radio transient detection., Comment: Accepted by AJ. 17 pages with figures

We have developed a new software tool, Lensview, for modelling resolved gravitational lens images. Based on the LensMEM algorithm, the software finds the best fitting lens mass model and source brightness distribution using a maximum entropy constraint. The method can be used with any point spread function or lens model. We review the algorithm and introduce some significant improvements. We also investigate and discuss issues associated with the statistical uncertainties of models and model parameters and the issues of source plane size and source pixel size. We test the software on simulated optical and radio data to evaluate how well lens models can be recovered and with what accuracy. For optical data, lens model parameters can typically be recovered with better than 1% accuracy and the degeneracy between mass ellipticity and power law is reduced. For radio data, we find that systematic errors associated with using processed radio maps, rather than the visibilities, are of similar magnitude to the random errors. Hence analysing radio data in image space is still useful and meaningful. The software is applied to the optical arc HST J15433+5352 and the radio ring MG1549+3047 using a simple elliptical isothermal lens model. For HST J15433+5352, the Einstein radius is 0.525" +/- 0.015 which probably includes a substantial convergence contribution from a neighbouring galaxy. For MG1549+3047, the model has Einstein radius 1.105" +/- 0.005 and core radius 0.16" 0.03. The total mass enclosed in the critical radius is 7.06 x 10^{10} Solar masses for our best model., Comment: 21 pages, 24 figures, appearing in MNRAS. Software available from http://www.cfa.harvard.edu/~rwayth/lensview/Lensview_Home.html

We examine the clustering properties of HI-selected galaxies through an analysis of the HI Parkes All-Sky Survey Catalogue (HICAT) two-point correlation function. Various sub-samples are extracted from this catalogue to study the overall clustering of HI-rich galaxies and its dependence on luminosity, HI gas mass and rotational velocity. These samples cover the entire southern sky Dec < 0 deg, containing up to 4,174 galaxies over the radial velocity range 300-12,700 km/s. A scale length of r_0 = 3.45 +/- 0.25 Mpc/h and slope of gamma = 1.47 +/- 0.08 is obtained for the HI-rich galaxy real-space correlation function, making gas-rich galaxies among the most weakly clustered objects known. HI-selected galaxies also exhibit weaker clustering than optically selected galaxies of comparable luminosities. Good agreement is found between our results and those of synthetic HI-rich galaxy catalogues generated from the Millennium Run CDM simulation. Bisecting HICAT using different parameter cuts, clustering is found to depend most strongly on rotational velocity and luminosity, while the dependency on HI mass is marginal. Splitting the sample around v_rot = 108 km/s, a scale length of r_0 = 2.86 +/- 0.46 Mpc/h is found for galaxies with low rotational velocities compared to r_0 = 3.96 +/- 0.33 Mpc/h for the high rotational velocity sample., Comment: Accepted for publication in the Astrophysical Journal

The high-redshift radio-loud quasar PKS 2126-158 is found to have a large number of red galaxies in close apparent proximity. We use the Gemini Multi-Object Spectrograph (GMOS) on Gemini South to obtain optical spectra for a large fraction of these sources. We show that there is a group of galaxies at $z\sim0.66$, coincident with a metal-line absorption system seen in the quasar's optical spectrum. The multiplexing capabilities of GMOS also allow us to measure redshifts of many foreground galaxies in the field surrounding the quasar. The galaxy group has five confirmed members, and a further four fainter galaxies are possibly associated. All confirmed members exhibit early-type galaxy spectra, a rare situation for a Mg II absorbing system. We discuss the relationship of this group to the absorbing gas, and the possibility of gravitational lensing of the quasar due to the intervening galaxies., Comment: Monthly Notices of the Royal Astronomical Society, in press. 10 pages, 8 figures

We have investigated the effect of spatial resolution on determining pencil-beam like velocity widths and column densities in galaxies. Three 21-cm datasets are used, the HIPASS galaxy catalogue, a subset of HIPASS galaxies with ATCA maps and a high-resolution image of the LMC. Velocity widths measured from 21-cm emission in local galaxies are compared with those measured in intermediate redshift Damped Lyman alpha (DLA) absorbers. We conclude that spatial resolution has a severe effect on measuring pencil-beam like velocity widths in galaxies. Spatial smoothing by a factor of 240 is shown to increase the median velocity width by a factor of two. Thus any difference between velocity widths measured from global profiles or low spatial resolution 21-cm maps at z=0 and DLAs at z>1 cannot unambiguously be attributed to galaxy evolution. The effect on column density measurements is less severe and the values of dN/dz from local low-resolution 21-cm measurements are expected to be overestimated by only ~10 per cent., Comment: 5 pages, 6 figures, accepted for publication in MNRAS letters

We present spatially resolved spectroscopic images of the gravitationally lensed QSO 2237+0305 taken with the GMOS Integral Field Unit (IFU) on the Gemini North telescope. These observations have the best spatial resolution of any IFU observations of this object to date and include the redshifted CIII] and MgII QSO broad lines. Unlike Mediavilla et al. 1998, we find no evidence for an arc of resolved broad line emission in either the CIII] or MgII lines. We calculate the image flux ratios of both the integrated emission lines and the surrounding continua. The flux ratios of the CIII] and MgII emission lines are consistent with each other but differ substantially from their corresponding continuum flux ratios and the radio/mid-IR flux ratios previously published. We argue that the broad emission line region must be microlensed and the CIII] and MgII emission regions must be approximately the same size and co-located along the line-of-sight. Assuming a simple model for the broad emission line region and the continuum region, we show the size of the CIII] / MgII broad line region is ~0.06 h_{70}^{1/2} pc and the continuum region is <= 0.02 h_{70}^{1/2} pc., Comment: 7 pages, 7 figures. Accepted for publication in MNRAS

We investigate whether a baryonic dark galaxy or `galaxy without stars' could persist indefinitely in the local universe, while remaining stable against star formation. To this end, a simple model has been constructed to determine the equilibrium distribution and composition of a gaseous protogalactic disk. Specifically, we determine the amount of gas that will transit to a Toomre unstable cold phase via the H2 cooling channel in the presence of a UV--X-ray cosmic background radiation field. All but one of the models are predicted to become unstable to star formation. Moreover, we find that all our model objects would be detectable via HI line emission, even in the case that star formation is potentially avoided. These results are consistent with the non-detection of isolated extragalactic HI clouds with no optical counterpart (galaxies without stars) by HIPASS. Additionally, where star formation is predicted to occur, we determine the minimum interstellar radiation field required to restore gravothermal stability, which we then relate to a minimum global star formation rate. This leads to the prediction of a previously undocumented relation between HI mass and star formation rate that is observed for a wide variety of dwarf galaxies in the HI mass range 10^8--10^10 M_sun. The existence of such a relation strongly supports the notion that the well observed population of dwarf galaxies represent the minimum rates of self-regulating star formation in the universe. (Barely abridged), Comment: 19 pages, 8 figures, TeX using emulateapj.cls, v2 accepted for publication in ApJ (16/8/5) with one figure deleted and a number of minor clarifying revisions

H1517+656 is an unusual source, even for a BL Lac object. It is one of the most luminous BL Lacs known, with extreme emission properties at radio, optical, and X-ray frequencies. Furthermore, in our WFPC2 snapshot survey we discovered a series of peculiar arcs describing a 2.4 arcsec radius ring surrounding the source. This paper describes follow-up observations with additional WFPC2 bands and the STIS longpass filter, which have revealed this structure to be the remnants of a very recent galaxy merger. Population synthesis modelling has shown that regions of the arcs have stellar populations with age < 20 Myrs. Additionally, the circularity of the arcs indicates that the plane of the collision and hence accretion is very close to the plane of the sky. Given that BL Lac jets are thought to be aligned with the line of sight, this observation may provide a direct link between the transfer of angular momentum in an interaction and the generation of a radio source., Comment: 29 pages, 9 figures, accepted for publication in ApJ

The location of the images in a multiple-image gravitational lens system are strongly dependent on the orientation angle of the mass distribution. As such, we can use the location of the images and the photometric properties of the visible matter to constrain the properties of the dark halo. We apply this to the optical Einstein Ring system 0047-2808 and find that the dark halo is almost spherical and is aligned in the same direction as the stars to within a few degrees., Comment: 4 pages, 2 figures. Conference proceedings of IAU 8th Asian-Pacific regional meeting in Tokyo

We determine the mass distribution of a spiral galaxy, 2237+0305 using both gravitational lensing and dynamical constraints. We find that lensing can break the disc-halo degeneracy. 2237+0305 has a sub-maximal disc, contributing 57+/-3 per cent of the rotational support at the disc maximum. The disc mass-to-light ratio is 1.1+/-0.2 in the I-band and the bulge, 2.9+/-0.5. The dark matter halo, modelled as a softened isothermal sphere, has a large core radius (13.4+/-0.4 kpc, equivalent to 1.4r_d) to high accuracy for the best-fitting solution. The image positions are reasonably well fitted, but require further rotation information to obtain a unique solution., Comment: 11 pages, 6 figures, accepted by MNRAS, in press

We measure the spectral properties of a representative sub-sample of 187 quasars, drawn from the Parkes Half-Jansky, Flat-radio-spectrum Sample (PHFS). Quasars with a wide range of rest-frame optical/UV continuum slopes are included in the analysis: their colours range from 2 < B-K < 7. The median H-beta and [O III] emission-line equivalent widths of the red quasar sub-sample are a factor of ten weaker than those of the blue quasar sub-sample. Both the colours and the emission-line equivalent widths of the red quasars can be explained by the addition of a featureless red synchrotron continuum component to an otherwise normal blue quasar spectrum. The relative strengths of the blue and red components span two orders of magnitude at rest-frame 500nm. The blue component is weaker relative to the red component in low optical luminosity sources. This suggests that the fraction of accretion energy going into optical emission from the jet is greater in low luminosity quasars. This synchrotron model does not, however, fit around 10% of the quasars, which have both red colours and high equivalent width emission-lines. We hypothesise that these red, strong-lined quasars have intrinsically weak Big Blue Bumps. There is no discontinuity in spectral properties between the BL Lac objects in our sample and the other quasars. The synchrotron emission component only dominates the spectrum at longer wavelengths, so existing BL Lac surveys will be biassed against high redshift objects., Comment: 22 pages, 12 figures, accepted for publication in PASA. Data tables and composite spectra from the paper can be found at http://msowww.anu.edu.au/~pfrancis/

We present model fits to spectral energy distributions in the optical and NIR of >100 flat-spectrum radio quasars from the Parkes Half-Jansky Flat-spectrum Sample. We find that approximately 40% of the sources have power-law SEDs, while a similar number show evidence for two primary components: a blue power law and optical synchrotron emission. The blue power law is similar to the dominant component observed in the spectra of optically-selected QSOs. There is strong evidence that the synchrotron component has a turnover in the UV-optical rest frame of the spectrum. In the remaining sources it is likely that the synchrotron peaks at longer wavelengths. This mixture of two components is supported by optical polarisation measurements in a subgroup of the sources. The sources with power law SEDs show evidence for an excess number of red power law slopes compared to optically-selected quasars. There are additional spectral components in some of the sources, such as dust and the underlying galaxy, which have not been considered here., Comment: MNRAS, in press; 16 pages, 13 figures

If the statistics of quasar lensing are to be used to constrain the cosmological model, the various potential systematic uncertainties in the deflector population must be accounted for. One ambiguity in the lensing properties of the dominant deflectors - elliptical galaxies - is their dynamical normalisation. Assuming they can be modelled as isothermal spheres, it is shown that their mass normalisation is sufficiently high that current upper limits on the cosmological constant remain valid, and in fact probably strengthened, even if ellipticals have significant core radii., Comment: New Cosmological Data and the Values of the Fundamental Parameters, A.N. Lasenby & A. Wilkinson, eds.; 2 pages, 1 figure

Gravitationally-lensed quasars should be discovered as a by-product of large galaxy redshift surveys, being discovered spectroscopically when a low-redshift galaxy exhibits high-redshift quasar emission lines. The number of lenses expected is higher than previously estimated, mainly due to the fact that the presence of the quasar images brings faint deflector galaxies above the survey limit. Thus the a posteriori likelihood of the discovery of Q 2237+0305 in the Center for Astrophysics redshift survey is approximately 0.03. In the future, the 2 degree Field survey should yield at least 10 lensed quasars, and the Sloan Digitial Sky Survey up to 100., Comment: Gravitational Lensing: Recent Progress and Future Goals, C.S. Kochanek & T.G. Brainerd, eds., in press; 2 pages, 1 figure

Galaxy redshift surveys can be used to detect gravitationally-lensed quasars if the spectra obtained are searched for the quasars' emission lines. Previous investigations of this possibility have used simple models to show that the 2 degree Field (2dF) redshift survey could yield several tens of new lenses, and that the larger Sloan Digital Sky Survey should contain an order of magnitude more. However the particular selection effects of the samples were not included in these calculations, limiting the robustness of the predictions; thus a more detailed simulation of the 2dF survey was undertaken here. The use of an isophotal magnitude limit reduces both the depth of the sample and the expected number of lenses, but more important is the Automatic Plate Measuring survey's star-galaxy separation algorithm, used to generate the 2dF input catalogue. It is found that most quasar lenses are classed as merged stars, with only the few lenses with low-redshift deflectors likely to be classified as galaxies. Explicit inclusion of these selection effects implies that the 2dF survey should contain 10 lenses on average. The largest remaining uncertainty is the lack of knowledge of the ease with which any underlying quasars can be extracted from the survey spectra., Comment: MNRAS, in press; 14 pages, 19 figures

The fraction of high-redshift sources which are multiply-imaged by intervening galaxies is strongly dependent on the cosmological constant, and so can be a useful probe of the cosmological model. However its power is limited by various systematic (and random) uncertainties in the calculation of lensing probabilities, one of the most important of which is the dynamical normalisation of elliptical galaxies. Assuming ellipticals' mass distributions can be modelled as isothermal spheres, the mass normalisation depends on: the velocity anisotropy; the luminosity density; the core radius; and the area over which the velocity dispersion is measured. The differences in the lensing probability and optical depth produced by using the correct normalisation can be comparable to the differences between even the most extreme cosmological models. The existing data is not sufficient to determine the correct normalisation with enough certainty to allow lensing statistics to be used to their full potential. However, as the correct lensing probability is almost certainly higher than is usually assumed, upper bounds on the cosmological constant are not weakened by these possibilities., Comment: MNRAS, in press; 13 pages, 22 figures

Gravitationally-lensed quasars can be discovered as a by-product of galaxy redshift surveys. Lenses discovered spectroscopically in this way should require less observational effort per event than those found in dedicated lens surveys. Further, the lens galaxies should be relatively nearby, facilitating a number of detailed observations that are impossible for the more common high-redshift lenses. This is epitomised by the wide range of results that have been obtained from Q 2237+0305, which was discovered as part of the Center for Astrophysics redshift survey, and remains the only quasar lens discovered in this way. The likelihood of this survey yielding a lens is calculated to be about 0.03, which is an order of magnitude larger than previous estimates due to two effects. Firstly, the quasar images themselves increase the observed flux of the lens, so that lens galaxies up to a magnitude fainter than the nominal survey limit must be included in the calculation. Secondly, it is possible for lensed quasars with extremely faint deflectors to enter the survey due to the extended morphology of the multiple images. Extrapolating these results to future surveys, the 2 degree Field galaxy redshift survey should contain between 10 and 50 lenses and the Sloan Digital Sky Survey should yield between 50 and 300 lenses, depending on the cosmological model and the observing conditions., Comment: MNRAS, in press; 15 pages, 17 figures

The absence of any wide-separation gravitational lenses in the Large Bright Quasar Survey is used to place limits on the population of cluster-sized halos in the universe, and hence constrain a number of cosmological parameters. The results agree with previous investigations in strongly ruling out the standard cold dark matter model but they are consistent with low-density universes in which the primordial fluctuation spectrum matches both cluster abundances and cosmic microwave background measurements. These conclusions are essentially independent of the cosmological constant, which is in stark contrast to the statistics of galaxy lenses. The constraints presented here are nullified if clusters have core radii of more than 10 kpc, but are free of a number of potential systematic errors, due to the homogeneity of the data., Comment: MNRAS, in press; 7 pages, 8 figures

Quasar pairs are either physically distinct binary quasars or the result of gravitational lensing. The majority of known pairs are in fact lenses, with a few confirmed as binaries, leaving a population of objects that have not yet been successfully classified. Building on the arguments of Kochanek, Falco & Munoz (1999), it is shown that there are no objective reasons to reject the binary interpretation for most of these. In particular, the similarity of the spectra of the quasar pairs appears to be an artifact of the generic nature of quasar spectra. The two ambiguous pairs discovered as part of the Large Bright Quasar Survey (Q 1429-053 and Q 2153-0256) are analysed using principle components analysis, which shows that their spectral similarities are not greater than expected for a randomly chosen pair of quasars from the survey. The assumption of the binary hypothesis allows the dynamics, time-scales and separation distribution of binary quasars to be investigated and constrained. The most plausible model is that the quasars' activity is triggered by tidal interactions in a galatic merger, but that the (re-)activation of the galactic nuclei occurs quite late in the interaction, when the nuclei are within 80+/-30 kpc of each other. A simple dynamical friction model for the decaying orbits reproduces the observed distribution of projected separations, but the decay time inferred is comparable to a Hubble time. Hence it is predicted that binary quasars are only observable as such in the early stages of galactic collisions, after which the quiescent super-massive black holes orbit in the merger remnant for some time., Comment: 12 pages, 12 figures

We present quasi-simultaneous multi-colour optical/near-IR photometry for 157 radio selected quasars, forming an unbiassed sub-sample of the Parkes Flat-Spectrum Sample. Data are also presented for 12 optically selected QSOs, drawn from the Large Bright QSO Survey. The spectral energy distributions of the radio- and optically-selected sources are quite different. The optically selected QSOs are all very similar: they have blue spectral energy distributions curving downwards at shorter wavelengths. Roughly 90% of the radio-selected quasars have roughly power-law spectral energy distributions, with slopes ranging from Fnu=nu^0 to Fnu=nu^-2. The remaining 10% have spectral energy distributions showing sharp peaks: these are radio galaxies and highly reddened quasars. Four radio sources were not detected down to magnitude limits of H ~ 19.6. These are probably high redshift (z >3) galaxies or quasars. We show that the colours of our red quasars lie close to the stellar locus in the optical: they will be hard to identify in surveys such as the Sloan Digital Sky Survey. If near-IR photometry is added, however, the red power-law sources can be clearly separated from the stellar locus: IR surveys such as 2MASS should be capable of finding these sources on the basis of their excess flux in the K-band., Comment: 20 pages text, four pages (separate) of data tables, 14 figures. Accepted for publication in PASA. Tables can be obtained in machine format from http://msowww.anu.edu.au/~pfrancis/publications.html

We describe an algorithm that measures self-consistently the relative galaxy contribution in a sample of radio-quasars from their optical spectra alone. This is based on a spectral fitting method which uses the size of the characteristic 4000\AA~ feature of elliptical galaxy SEDs. We apply this method to the Parkes Half-Jansky Flat Spectrum sample of Drinkwater et al. (1997) to determine whether emission from the host galaxy can significantly contribute to the very red optical-to-near-infrared colours observed. We find that at around $2\sigma$ confidence, most of the reddening in unresolved (mostly quasar-like) sources is unlikely to be due to contamination by a red stellar component., Comment: 11 pages, 11 figures. Accepted for Publication in Monthly Notices of the Royal Astronomical Society

The Parkes Half-Jansky Flat Spectrum Sample contains a large number of sources with unusually red optical-to-near-infrared continua. If this is to be interpreted as extinction by dust in the line-of-sight, then associated material might also give rise to absorption in the soft X-ray regime. This hypothesis is tested using broadband (0.1-2.4 keV) data from the {\it ROSAT} All-Sky Survey provided by Siebert et al. (1998). Significant ($>3\sigma$ confidence level) correlations between optical (and near-infrared)--to--soft X-ray continuum slope and optical extinction are found in the data, consistent with absorption by material with metallicity and a range in gas-to-dust ratio as observed in the local ISM. Under this simple model, the soft X-rays are absorbed at a level consistent with the range of extinctions ($0< A_{V}< 6$ magnitudes) implied by the observed optical reddening. Excess X-ray absorption by warm (ionised) gas, (ie. a `warm absorber') is not required., Comment: 23 pages of text, 3 figures, to appear in Jan 10 (1999) issue of The Astrophysical Journal

In the gravitational lens system Q2237+0305 the cruciform quasar image geometry is twisted by ten degrees by the lens effect of a bar in the lensing galaxy. This effect can be used to measure the mass of the bar. We construct a new lensing model for this system with a power-law elliptical bulge and a Ferrers bar. The observed ellipticity of the optical isophotes of the galaxy leads to a nearly isothermal elliptical profile for the bulge with a total quasar magnification of 16^{+5}_{-4}. We measure a bar mass of 7.5\pm1.5 \times 10^8 h^{-1}_{75} M_{\sun} in the region inside the quasar images., Comment: 9 pages, 5 Postscript figures, uses mn.sty and eps.sty, submitted to MNRAS