Correlations between black holes and their host galaxies provide insight into what drives black hole-host co-evolution. We use the Meraxes semi-analytic model to investigate the growth of black holes and their host galaxies from high redshift to the present day. Our modelling finds no significant evolution in the black hole-bulge and black hole-total stellar mass relations out to a redshift of 8. The black hole-total stellar mass relation has similar but slightly larger scatter than the black hole-bulge relation, with the scatter in both decreasing with increasing redshift. In our modelling the growth of galaxies, bulges and black holes are all tightly related, even at the highest redshifts. We find that black hole growth is dominated by instability-driven or secular quasar-mode growth and not by merger-driven growth at all redshifts. Our model also predicts that disc-dominated galaxies lie on the black hole-total stellar mass relation, but lie offset from the black hole-bulge mass relation, in agreement with recent observations and hydrodynamical simulations.

We compute the spherically-averaged power spectrum from four seasons of data obtained for the Epoch of Reionisation (EoR) project observed with the Murchison Widefield Array (MWA). We measure the EoR power spectrum over $k= 0.07-3.0~h$Mpc$^{-1}$ at redshifts $z=6.5-8.7$. The largest aggregation of 110 hours on EoR0 high-band (3,340 observations), yields a lowest measurement of (43~mK)$^2$ = 1.8$\times$10$^3$ mK$^2$ at $k$=0.14~$h$Mpc$^{-1}$ and $z=6.5$ (2$\sigma$ thermal noise plus sample variance). Using the Real-Time System to calibrate and the CHIPS pipeline to estimate power spectra, we select the best observations from the central five pointings within the 2013--2016 observing seasons, observing three independent fields and in two frequency bands. This yields 13,591 2-minute snapshots (453 hours), based on a quality assurance metric that measures ionospheric activity. We perform another cut to remove poorly-calibrated data, based on power in the foreground-dominated and EoR-dominated regions of the two-dimensional power spectrum, reducing the set to 12,569 observations (419 hours). These data are processed in groups of 20 observations, to retain the capacity to identify poor data, and used to analyse the evolution and structure of the data over field, frequency, and data quality. We subsequently choose the cleanest 8,935 observations (298 hours of data) to form integrated power spectra over the different fields, pointings and redshift ranges.

The compact configuration of Phase II of the Murchison Widefield Array (MWA) consists of both a redundant subarray and pseudo-random baselines, offering unique opportunities to perform sky-model and redundant interferometric calibration. The highly redundant hexagonal cores give improved power spectrum sensitivity. In this paper, we present the analysis of nearly 40 hours of data targeting one of the MWA's EoR fields observed in 2016. We use both improved analysis techniques presented in Barry et al. (2019) as well as several additional techniques developed for this work, including data quality control methods and interferometric calibration approaches. We show the EoR power spectrum limits at redshift 6.5, 6.8 and 7.1 based on our deep analysis on this 40-hour data set. These limits span a range in $k$ space of $0.18$ $h$ $\mathrm{Mpc^{-1}}$ $<k<1.6$ $h$ $\mathrm{Mpc^{-1}}$, with a lowest measurement of $\Delta^2\leqslant2.39\times 10^3$ $\mathrm{mK}^2$ at $k=0.59$ $h$ $\mathrm{Mpc^{-1}}$ and $z=6.5$.

We examine the properties of the host galaxies of $z=7$ quasars using the large volume, cosmological hydrodynamical simulation BlueTides. We find that the most massive black holes and quasars are hosted by massive galaxies with stellar masses $\log(M_\ast/M_\odot)=10.8\pm0.2$, and $10.2\pm0.4$, which have a wide range of star formation rates, of $50\substack{+119 \\ -35}M_\odot/\rm{yr}$ and $19\substack{+28 \\ -12}M_\odot/\rm{yr}$, respectively. The hosts of the most massive black holes and quasars in BlueTides are generally bulge-dominated, with bulge-to-total mass ratio $B/T\simeq0.85\pm0.1$, however their morphologies are not biased relative to the overall $z=7$ galaxy sample. We find that the hosts of the most massive black holes and quasars are significantly more compact, with half-mass radii $R_{0.5}=0.29\substack{+0.15 \\ -0.10}$ kpc and $0.28\substack{+0.08 \\ -0.06}$ kpc respectively, relative to galaxies with similar masses, which have $R_{0.5}=0.5\substack{+0.3 \\ -0.2}$ kpc. We make mock James Webb Space Telescope (JWST) images of these quasars and their host galaxies. We find that distinguishing the host from the quasar emission will be possible but still challenging with JWST, due to the small sizes of quasar hosts. We find that currently observable quasar samples are biased tracers of the intrinsic black hole--stellar mass relations, following a relation that is 0.27 dex higher than that of the full galaxy sample. Finally, we find that black hole hosts are generally quite isolated. However, the most massive black holes are more likely to be found in denser environments than the typical $M_{\textrm{BH}}>10^{6.5}M_\odot$ black hole, indicating that minor mergers at least play some role in growing black holes in the early Universe.

Measurements of 21 cm Epoch of Reionization (EoR) structure are subject to systematics originating from both the analysis and the observation conditions. Using 2013 data from the Murchison Widefield Array (MWA), we show the importance of mitigating both sources of contamination. A direct comparison between results from Beardsley et al. 2016 and our updated analysis demonstrates new precision techniques, lowering analysis systematics by a factor of 2.8 in power. We then further lower systematics by excising observations contaminated by ultra-faint RFI, reducing by an additional factor of 3.8 in power for the zenith pointing. With this enhanced analysis precision and newly developed RFI mitigation, we calculate a noise-dominated upper limit on the EoR structure of $\Delta^2 \leq 3.9 \times 10^3$ mK$^2$ at $k=0.20$ $\textit{h}$ Mpc$^{-1}$ and $z=7$ using 21 hr of data, improving previous MWA limits by almost an order of magnitude.

Intergalactic medium temperature is a powerful probe of the epoch of reionisation, as information is retained long after reionisation itself. However, mean temperatures are highly degenerate with the timing of reionisation, with the amount heat injected during the epoch, and with the subsequent cooling rates. We post-process a suite of semi-analytic galaxy formation models to characterise how different thermal statistics of the intergalactic medium can be used to constrain reionisation. Temperature is highly correlated with redshift of reionisation for a period of time after the gas is heated. However as the gas cools, thermal memory of reionisation is lost, and a power-law temperature-density relation is formed, $T = T_0(1+\delta)^{1-\gamma}$ with $\gamma \approx 1.5$. Constraining our model against observations of electron optical depth and temperature at mean density, we find that reionisation likely finished at $z_{\rm{reion}} = 6.8 ^{+ 0.5} _{-0.8}$ with a soft spectral slope of $\alpha = 2.8 ^{+ 1.2} _{-1.0}$. By restricting spectral slope to the range $[0.5,2.5]$ motivated by population II synthesis models, reionisation timing is further constrained to $z_{\rm{reion}} = 6.9 ^{+ 0.4} _{-0.5}$. We find that, in the future, the degeneracies between reionisation timing and background spectrum can be broken using the scatter in temperatures and integrated thermal history.

We present a new analysis of high-redshift UV observations using a semi-analytic galaxy formation model, and provide self-consistent predictions of the infrared excess (IRX) -- $\beta$ relations and cosmic star formation rate density. We combine the Charlot & Fall dust attenuation model with the Meraxes semi-analytic model, and explore three different parametrisations for the dust optical depths, linked to star formation rate, dust-to-gas ratio and gas column density respectively. A Bayesian approach is employed to statistically calibrate model free parameters including star formation efficiency, mass loading factor, dust optical depths and reddening slope directly against UV luminosity functions and colour-magnitude relations at z ~ 4-7. The best-fit models show excellent agreement with the observations. We calculate IRX using energy balance arguments, and find that the large intrinsic scatter in the IRX -$\beta$ plane is driven by the specific star formation rate. Additionally, the difference among the three dust models suggests a factor of two systematic uncertainty in the dust-corrected star formation rate when using the Meurer IRX - $\beta$ relation at z > 4.

We study the sizes, angular momenta and morphologies of high-redshift galaxies using an update of the Meraxes semi-analytic galaxy evolution model. Our model successfully reproduces a range of observations from redshifts $z=0$-$10$. We find that the effective radius of a galaxy disc scales with UV luminosity as $R_e\propto L_{\textrm{UV}}^{0.33}$ at $z=5$-$10$, and with stellar mass as $R_e\propto M_\ast^{0.24}$ at $z=5$ but with a slope that increases at higher redshifts. Our model predicts that the median galaxy size scales with redshift as $R_e \propto (1+z)^{-m}$, where $m=1.98\pm0.07$ for galaxies with $(0.3$-$1)L^\ast_{z=3}$ and $m=2.15\pm0.05$ for galaxies with $(0.12$-$0.3)L^\ast_{z=3}$. We find that the ratio between stellar and halo specific angular momentum is typically less than one and decreases with halo and stellar mass. This relation shows no redshift dependence, while the relation between specific angular momentum and stellar mass decreases by $\sim0.5$ dex from $z=7$ to $z=2$. Our model reproduces the distribution of local galaxy morphologies, with bulges formed predominantly through galaxy mergers for low-mass galaxies, disc-instabilities for galaxies with $M_\ast\simeq10^{10}$-$10^{11.5}M_\odot$, and major mergers for the most massive galaxies. At high redshifts, we find galaxy morphologies that are predominantly bulge-dominated.

We apply two methods to estimate the 21~cm bispectrum from data taken within the Epoch of Reionisation (EoR) project of the Murchison Widefield Array (MWA). Using data acquired with the Phase II compact array allows a direct bispectrum estimate to be undertaken on the multiple redundantly-spaced triangles of antenna tiles, as well as an estimate based on data gridded to the $uv$-plane. The direct and gridded bispectrum estimators are applied to 21 hours of high-band (167--197~MHz; $z$=6.2--7.5) data from the 2016 and 2017 observing seasons. Analytic predictions for the bispectrum bias and variance for point source foregrounds are derived. We compare the output of these approaches, the foreground contribution to the signal, and future prospects for measuring the bispectra with redundant and non-redundant arrays. We find that some triangle configurations yield bispectrum estimates that are consistent with the expected noise level after 10 hours, while equilateral configurations are strongly foreground-dominated. Careful choice of triangle configurations may be made to reduce foreground bias that hinders power spectrum estimators, and the 21~cm bispectrum may be accessible in less time than the 21~cm power spectrum for some wave modes, with detections in hundreds of hours.

We explore methods for robust estimation of the 21 cm signal from the Epoch of Reionisation (EoR). A Kernel Density Estimator (KDE) is introduced for measuring the spatial temperature fluctuation power spectrum from the EoR. The KDE estimates the underlying probability distribution function of fluctuations as a function of spatial scale, and contains different systematic biases and errors to the typical approach to estimating the fluctuation power spectrum. Extraction of histograms of visibilities allows moments analysis to be used to discriminate foregrounds from 21 cm signal and thermal noise. We use the information available in the histograms, along with the statistical dis-similarity of foregrounds from two independent observing fields, to robustly separate foregrounds from cosmological signal, while making no assumptions about the Gaussianity of the signal. Using two independent observing fields to robustly discriminate signal from foregrounds is crucial for the analysis presented in this paper. We apply the techniques to 13 hours of Murchison Widefield Array (MWA) EoR data over two observing fields. We compare the output to that obtained with a comparative power spectrum estimation method, and demonstrate the reduced foreground contamination using this approach. Using the second moment obtained directly from the KDE distribution functions yields a factor of 2-3 improvement in power for k < 0.3hMpc^{-1} compared with a matched delay space power estimator, while weighting data by additional statistics does not offer significant improvement beyond that available for thermal noise-only weights.

We investigate the dependence of galaxy clustering at $z \sim 4 - 7$ on UV-luminosity and stellar mass. Our sample consists of $\sim$ 10,000 Lyman-break galaxies (LBGs) in the XDF and CANDELS fields. As part of our analysis, the $M_\star - M_{\rm UV}$ relation is estimated for the sample, which is found to have a nearly linear slope of $d\log_{10} M_\star / d M_{\rm UV} \sim 0.44$. We subsequently measure the angular correlation function and bias in different stellar mass and luminosity bins. We focus on comparing the clustering dependence on these two properties. While UV-luminosity is only related to recent starbursts of a galaxy, stellar mass reflects the integrated build-up of the whole star formation history, which should make it more tightly correlated with halo mass. Hence, the clustering segregation with stellar mass is expected to be larger than with luminosity. However, our measurements suggest that the segregation with luminosity is larger with $\simeq 90\%$ confidence (neglecting contributions from systematic errors). We compare this unexpected result with predictions from the \textsc{Meraxes} semi-analytic galaxy formation model. Interestingly, the model reproduces the observed angular correlation functions, and also suggests stronger clustering segregation with luminosity. The comparison between our observations and the model provides evidence of multiple halo occupation in the small scale clustering.

We describe the motivation and design details of the "Phase II" upgrade of the Murchison Widefield Array (MWA) radio telescope. The expansion doubles to 256 the number of antenna tiles deployed in the array. The new antenna tiles enhance the capabilities of the MWA in several key science areas. Seventy-two of the new tiles are deployed in a regular configuration near the existing MWA core. These new tiles enhance the surface brightness sensitivity of the MWA and will improve the ability of the MWA to estimate the slope of the Epoch of Reionisation power spectrum by a factor of ~3.5. The remaining 56 tiles are deployed on long baselines, doubling the maximum baseline of the array and improving the array u,v coverage. The improved imaging capabilities will provide an order of magnitude improvement in the noise floor of MWA continuum images. The upgrade retains all of the features that have underpinned the MWA's success (large field-of-view, snapshot image quality, pointing agility) and boosts the scientific potential with enhanced imaging capabilities and by enabling new calibration strategies.

Structure imprinted in foreground extragalactic point sources by ionospheric refraction has the potential to contaminate Epoch of Reionisation (EoR) power spectra of the 21~cm emission line of neutral hydrogen. The alteration of the spatial and spectral structure of foreground measurements due to total electron content (TEC) gradients in the ionosphere create a departure from the expected sky signal. We present a general framework for understanding the signatures of ionospheric behaviour in the two-dimensional (2D) neutral hydrogen power spectrum measured by a low-frequency radio interferometer. Two primary classes of ionospheric behaviour are considered, corresponding to dominant modes observed in Murchison Widefield Array (MWA) EoR data; namely, anisotropic structured wave behaviour, and isotropic turbulence. Analytic predictions for power spectrum bias due to this contamination are computed, and compared with simulations. We then apply the ionospheric metric described in Jordan et al. (2017) to study the impact of ionospheric structure on MWA data, by dividing MWA EoR datasets into classes with good and poor ionospheric conditions, using sets of matched 30-minute observations from 2014 September. The results are compared with the analytic and simulated predictions, demonstrating the observed bias in the power spectrum when the ionosphere is active (displays coherent structures or isotropic turbulence). The analysis demonstrates that unless ionospheric activity can be quantified and corrected, active data should not be included in EoR analysis in order to avoid systematic biases in cosmological power spectra. When data are corrected with a model formed from the calibration information, bias reduces below the expected 21~cm signal level. Data are considered `quiet' when the median measured source position offsets are less than 10-15~arcseconds.

We directly compare predictions of dwarf galaxy properties in a semi-analytic model (SAM) with those extracted from a high-resolution hydrodynamic simulation. We focus on galaxies with halo masses of 1e9<Mvir/Msol<1e11 at high redshift ($z\ge5$). We find that, with the modifications previously proposed in Qin et al. (2018), including to suppress the halo mass and baryon fraction, as well as to modulate gas cooling and star formation efficiencies, the SAM can reproduce the cosmic evolution of galaxy properties predicted by the hydrodynamic simulation. These include the galaxy stellar mass function, total baryonic mass, star-forming gas mass and star formation rate at $z\sim5-11$. However, this agreement is only possible by reducing the star formation threshold relative to that suggested by local observations. Otherwise, too much star-forming gas is trapped in quenched dwarf galaxies. We further find that dwarf galaxies rapidly build up their star-forming reservoirs in the early universe ($z>10$), with the relevant time-scale becoming significantly longer towards lower redshifts. This indicates efficient accretion in cold mode in these low-mass objects at high redshift. Note that the improved SAM, which has been calibrated against hydrodynamic simulations, can provide more accurate predictions of high-redshift dwarf galaxy properties that are essential for reionization study.

We study dwarf galaxy formation at high redshift ($z\ge5$) using a suite of high- resolution, cosmological hydrodynamic simulations and a semi-analytic model (SAM). We focus on gas accretion, cooling and star formation in this work by isolating the relevant process from reionization and supernova feedback, which will be further discussed in a companion paper. We apply the SAM to halo merger trees constructed from a collisionless N-body simulation sharing identical initial conditions to the hydrodynamic suite, and calibrate the free parameters against the stellar mass function predicted by the hydrodynamic simulations at z = 5. By making comparisons of the star formation history and gas components calculated by the two modelling techniques, we find that semi-analytic prescriptions that are commonly adopted in the literature of low-redshift galaxy formation do not accurately represent dwarf galaxy properties in the hydrodynamic simulation at earlier times. We propose 3 modifications to SAMs that will provide more accurate high-redshift simulations. These include 1) the halo mass and baryon fraction which are overestimated by collisionless N-body simulations; 2) the star formation efficiency which follows a different cosmic evolutionary path from the hydrodynamic simulation; and 3) the cooling rate which is not well defined for dwarf galaxies at high redshift. Accurate semi-analytic modelling of dwarf galaxy formation informed by detailed hydrodynamical modelling will facilitate reliable semi-analytic predictions over the large volumes needed for the study of reionization.

We investigate the physics that drives the gas accretion rates onto galaxies at the centers of dark matter haloes using the EAGLE suite of hydrodynamical cosmological simulations. We find that at redshifts $z{\le}2$ the accretion rate onto the galaxy increases with halo mass in the halo mass range $10^{10}-10^{11.7}M_{\odot}$, flattens between the halo masses $10^{11.7}-10^{12.7}M_{\odot}$, and increases again for higher-mass haloes. However, the galaxy gas accretion does not flatten at intermediate halo masses when AGN feedback is switched off. To better understand these trends, we develop a physically motivated semi-analytic model of galaxy gas accretion. We show that the flattening is produced by the rate of gas cooling from the hot halo. The ratio of the cooling radius and the virial radius does not decrease continuously with increasing halo mass as generally thought. While it decreases up to ${\sim}10^{13}M_{\odot}$ haloes, it increases for higher halo masses, causing an upturn in the galaxy gas accretion rate. This may indicate that in high-mass haloes AGN feedback is not sufficiently efficient. When there is no AGN feedback, the density of the hot halo is higher, the ratio of the cooling and virial radii does not decrease as much and the cooling rate is higher. Changes in the efficiency of stellar feedback can also increase or decrease the accretion rates onto galaxies. The trends can plausibly be explained by the re-accretion of gas ejected by progenitor galaxies and by the suppression of black hole growth, and hence AGN feedback, by stellar feedback.

HI intensity mapping data traces the large-scale structure matter distribution using the integrated emission of neutral hydrogen gas (HI). The cross-correlation of the intensity maps with optical galaxy surveys can mitigate foreground and systematic effects, but has been shown to significantly depend on galaxy evolution parameters of the HI and the optical sample. Previously, we have shown that the shot noise of the cross-correlation scales with the HI content of the optical samples, such that the shot noise estimation infers the average HI masses of these samples. In this article, we present an adaptive framework for the cross-correlation of HI intensity maps with galaxy samples using our implementation of the halo model formalism (Murray et al 2018, in prep) which utilises the halo occupation distribution of galaxies to predict their power spectra. We compare two HI population models, tracing the spatial halo and the galaxy distribution respectively, and present their auto- and cross-power spectra with an associated galaxy sample. We find that the choice of the HI model and the distribution of the HI within the galaxy sample have minor significance for the shape of the auto- and cross-correlations, but highly impact the measured shot noise amplitude of the estimators, a finding we confirm with simulations. We demonstrate parameter estimation of the HI halo occupation models and advocate this framework for the interpretation of future experimental data, with the prospect of determining the HI masses of optical galaxy samples via the cross-correlation shot noise.

We summarize panchromatic Extragalactic Background Light data to place upper limits on the integrated near-infrared surface brightness (SB) that may come from Population III stars and possible accretion disks around their stellar-mass black holes (BHs) in the epoch of First Light, broadly taken from z$\simeq$7-17. Theoretical predictions and recent near-infrared power-spectra provide tighter constraints on their sky-signal. We outline the physical properties of zero metallicity Population III stars from MESA stellar evolution models through helium-depletion and of BH accretion disks at z$\gtrsim$7. We assume that second-generation non-zero metallicity stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions. We use these near-infrared SB constraints to calculate the number of caustic transits behind lensing clusters that the James Webb Space Telescope and the next generation ground-based telescopes may observe for both Population III stars and their BH accretion disks. Typical caustic magnifications can be $\mu$$\simeq$10$^4$-10$^5$, with rise times of hours and decline times of $\lesssim$1 year for cluster transverse velocities of $v_{T}$$\lesssim$1000 km s$^{-1}$. Microlensing by intracluster medium objects can modify transit magnifications, but lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Population III stars. To observe Population III caustic transits directly may require to monitor 3-30 lensing clusters to AB$\lesssim$29 mag over a decade.

We present the results of a VLT MUSE/FORS2 and Spitzer survey of a unique compact lensing cluster CLIO at z = 0.42, discovered through the GAMA survey using spectroscopic redshifts. Compact and massive clusters such as this are understudied, but provide a unique prospective on dark matter distributions and for finding background lensed high-z galaxies. The CLIO cluster was identified for follow up observations due to its almost unique combination of high mass and dark matter halo concentration, as well as having observed lensing arcs from ground based images. Using dual band optical and infra-red imaging from FORS2 and Spitzer, in combination with MUSE optical spectroscopy we identify 89 cluster members and find background sources out to z = 6.49. We describe the physical state of this cluster, finding a strong correlation between environment and galaxy spectral type. Under the assumption of a NFW profile, we measure the total mass of CLIO to be M$_{200} = (4.49 \pm 0.25) \times 10^{14}$ M$_\odot$. We build and present an initial strong-lensing model for this cluster, and measure a relatively low intracluster light (ICL) fraction of 7.21 $\pm$ 1.53% through galaxy profile fitting. Due to its strong potential for lensing background galaxies and its low ICL, the CLIO cluster will be a target for our 110 hour JWST 'Webb Medium-Deep Field' (WMDF) GTO program.

Motivated by recent measurements of the number density of faint AGN at high redshift, we investigate the contribution of quasars to reionization by tracking the growth of central supermassive black holes in an update of the Meraxes semi-analytic model. The model is calibrated against the observed stellar mass function at $z\sim0.6-7$, the black hole mass function at $z\lesssim0.5$, the global ionizing emissivity at $z\sim2-5$ and the Thomson scattering optical depth. The model reproduces a Magorrian relation in agreement with observations at $z<0.5$ and predicts a decreasing black hole mass towards higher redshifts at fixed total stellar mass. With the implementation of an opening angle of 80 deg for quasar radiation, corresponding to an observable fraction of ${\sim}23.4$ per cent due to obscuration by dust, the model is able to reproduce the observed quasar luminosity function at $z\sim0.6-6$. The stellar light from galaxies hosting faint AGN contributes a significant or dominant fraction of the UV flux. At high redshift, the model is consistent with the bright end quasar luminosity function and suggests that the recent faint $z\sim4$ AGN sample compiled by Giallongo et al. (2015) includes a significant fraction of stellar light. Direct application of this luminosity function to the calculation of AGN ionizing emissivity consequently overestimates the number of ionizing photons produced by quasars by a factor of 3 at $z\sim6$. We conclude that quasars are unlikely to make a significant contribution to reionization.

We investigate the clustering properties of Lyman-break galaxies (LBGs) at $z\sim6$ - $8$. Using the semi-analytical model {\scshape Meraxes} constructed as part of the Dark-ages Reionization And Galaxy-formation Observables from Numerical Simulation (DRAGONS) project, we predict the angular correlation function (ACF) of LBGs at $z\sim6$ - $8$. Overall, we find that the predicted ACFs are in good agreement with recent measurements at $z\sim 6$ and $z\sim 7.2$ from observations consisting of the Hubble eXtreme Deep Field (XDF), the Hubble Ultra-Deep Field (HUDF) and Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) field. We confirm the dependence of clustering on luminosity, with more massive dark matter haloes hosting brighter galaxies, remains valid at high redshift. The predicted galaxy bias at fixed luminosity is found to increase with redshift, in agreement with observations. We find that LBGs of magnitude $M_{{\rm AB(1600)}} < -19.4$ at $6\lesssim z \lesssim 8$ reside in dark matter haloes of mean mass $\sim 10^{11.0}$- $10^{11.5} M_{\rm \odot}$, and this dark matter halo mass does not evolve significantly during reionisation.

We use a suite of hydrodynamical cosmological simulations from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) project to investigate the formation of hot hydrostatic haloes and their dependence on feedback mechanisms. We find that the appearance of a strong bimodality in the probability density function (PDF) of the ratio of the radiative cooling and dynamical times for halo gas provides a clear signature of the formation of a hot corona. Haloes of total mass $10^{11.5}-10^{12}{\rm{M}}_{\odot}$ develop a hot corona independent of redshift, at least in the interval $z=0-4$ where the simulation has sufficiently good statistics. We analyse the build up of the hot gas mass in the halo, $M_{\rm{hot}}$, as a function of halo mass and redshift and find that while more energetic galactic winds powered by SNe increases $M_{\rm{hot}}$, AGN feedback reduces it by ejecting gas from the halo. We also study the thermal properties of gas accreting onto haloes and measure the fraction of shock-heated gas as a function of redshift and halo mass. We develop analytic and semianalytic approaches to estimate a `critical halo mass', $M_{\rm{crit}}$, for hot halo formation. We find that the mass for which the heating rate produced by accretion shocks equals the radiative cooling rate, reproduces the mass above which haloes develop a significant hot atmosphere. This yields a mass estimate of $M_{\rm{crit}} \approx 10^{11.7}{\rm{M}}_{\odot}$ at $z=0$, which agrees with the simulation results. The value of $M_{\rm{crit}}$ depends more strongly on the cooling rate than on any of the feedback parameters.

This work presents a study of galactic outflows driven by stellar feedback. We extract main sequence disc galaxies with stellar mass $10^9\le$ M$_{\star}/$M$_{\odot} \le 5.7\times10^{10}$ at redshift $z=0$ from the highest resolution cosmological simulation of the Evolution and Assembly of GaLaxies and their Environments (EAGLE) set. Synthetic gas rotation velocity and velocity dispersion ($\sigma$) maps are created and compared to observations of disc galaxies obtained with the Sydney-AAO Multi-object Integral field spectrograph (SAMI), where $\sigma$-values greater than $150$ km s$^{-1}$ are most naturally explained by bipolar outflows powered by starburst activity. We find that the extension of the simulated edge-on (pixelated) velocity dispersion probability distribution depends on stellar mass and star formation rate surface density ($\Sigma_{\rm SFR}$), with low-M$_{\star}/$low-$\Sigma_{\rm SFR}$ galaxies showing a narrow peak at low $\sigma$ ($\sim30$ km s$^{-1}$) and more active, high-M$_{\star}/$high-$\Sigma_{\rm SFR}$ galaxies reaching $\sigma>150$ km s$^{-1}$. Although supernova-driven galactic winds in the EAGLE simulations may not entrain enough gas with T $<10^5$ K compared to observed galaxies, we find that gas temperature is a good proxy for the presence of outflows. There is a direct correlation between the thermal state of the gas and its state of motion as described by the $\sigma$-distribution. The following equivalence relations hold in EAGLE: $i)$ low-$\sigma$ peak $\,\Leftrightarrow\,$ disc of the galaxy $\,\Leftrightarrow\,$ gas with T $<10^5$ K; $ii)$ high-$\sigma$ tail $\,\Leftrightarrow\,$ galactic winds $\,\Leftrightarrow\,$ gas with T $\ge 10^5$ K.

Massive quiescent galaxies (MQGs) are thought to have formed stars rapidly at early times followed by a long period of quiescence. The recent discovery of a MQG, ZF-COSMOS-20115 at $z\sim4$, only 1.5 Gyr after the big bang, places new constraints on galaxy growth and the role of feedback in early star formation. Spectroscopic follow-up confirmed ZF-COSMOS-20115 as a MQG at $z=3.717$ with an estimated stellar mass of ${\sim}10^{11}\mathrm{M}_\odot$, showing no evidence of recent star formation. We use the Meraxes semi-analytic model to investigate how ZF-COSMOS-20115 analogues build stellar mass, and why they become quiescent. We identify three analogue galaxies with similar properties to ZF-COSMOS-20115. We find that ZF-COSMOS-20115 is likely hosted by a massive halo with virial mass of ${\sim}10^{13}\mathrm{M}_\odot$, having been through significant mergers at early times. These merger events drove intense growth of the nucleus, which later prevented cooling and quenched star formation. ZF-COSMOS-20115 likely remained quiescent at $z<3.7$. We find that the analogues host the most massive black holes in our simulation and were luminous quasars at $z\sim5$, indicating that ZF-COSMOS-20115 and other MQGs may be the descendants of high-redshift quasars. In addition, the model suggests that ZF-COSMOS-20115 formed in a region of intergalactic medium that was reionized early.

Direct detection of regions of ionized hydrogen (HII) has been suggested as a promising probe of cosmic reionization. Observing the redshifted 21-cm signal of hydrogen from the epoch of reionization (EoR) is a key scientific driver behind new-generation, low-frequency radio interferometers. We investigate the feasibility of combining low-frequency observations with the Square Kilometre Array and near infra-red survey data of the Wide-Field Infrared Survey Telescope to detect cosmic reionization by imaging HII bubbles surrounding massive galaxies during the cosmic dawn. While individual bubbles will be too small to be detected, we find that by stacking redshifted 21-cm spectra centred on known galaxies, it will be possible to directly detect the EoR at $z \sim 9-12$, and to place qualitative constraints on the evolution of the spin temperature of the intergalactic medium (IGM) at $z \geq 9$. In particular, given a detection of ionized bubbles using this technique, it is possible to determine if the IGM surrounding them is typically in absorption or emission. Determining the globally-averaged neutral fraction of the IGM using this method will prove more difficult due to degeneracy with the average size of HII regions.

Using a series of high-resolution hydrodynamical simulations we show that during the rapid growth of high-redshift (z > 5) galaxies, reserves of molecular gas are consumed over a time-scale of 300Myr, almost independent of feedback scheme. We find that there exists no such simple relation for the total gas fractions of these galaxies, with little correlation between gas fractions and specific star formation rates. The bottleneck or limiting factor in the growth of early galaxies is in converting infalling gas to cold star-forming gas. Thus, we find that the majority of high redshift dwarf galaxies are effectively in recession, with demand (of star formation) never rising to meet supply (of gas), irrespective of the baryonic feedback physics modelled. We conclude that the basic assumption of self-regulation in galaxies - that they can adjust total gas consumption within a Hubble time - does not apply for the dwarf galaxies thought to be responsible for providing most UV photons to reionize the high redshift Universe. We demonstrate how this rapid molecular time-scale improves agreement between semi-analytic model predictions of the early Universe and observed stellar mass functions.

We present a theoretical study of intergalactic metal absorption lines imprinted in the spectra of distant quasars during and after the Epoch of Reionization (EoR). We use high resolution hydrodynamical simulations at high redshift ($4 <z<8$), assuming a uniform UV background Haardt--Madau 12, post-processing with CLOUDY photoionization models and Voigt profile fitting to accurately calculate column densities of the ions CII, CIV, SiII, SiIV and OI in the intergalactic medium (IGM). In addition, we generate mock observations of neutral Hydrogen (HI) at $z<6$. Our simulations successfully reproduce the evolution of the cosmological mass density ($\Omega$) of CII and CIV, with $\Omega_{CII}$ exceeding $\Omega_{CIV}$ at $z >6$, consistent with the current picture of the tail of the EoR. The simulated CII exhibits a bimodal distribution with large absorptions in and around galaxies, and some traces in the lower density IGM. We find some discrepancies between the observed and simulated column density relationships among different ionic species at $z=6$, probably due to uncertainties in the assumed UV background. Finally, our simulations are in good agreement with observations of the HI column density distribution function at $z = 4$ and the HI cosmological mass density $\Omega_{HI}$ at $4 < z < 6$.

We present a theoretical model to predict the properties of an observed $z =$ 5.72 Lyman $\alpha$ emitter galaxy - CIV absorption pair separated by 1384 comoving kpc/h. We use the separation of the pair and an outflow velocity/time travelling argument to demonstrate that the observed galaxy cannot be the source of metals for the CIV absorber. We find a plausible explanation for the metal enrichment in the context of our simulations: a dwarf galaxy with $M_{\star} =$ 1.87 $\times$ 10$^{9} M_{\odot}$ located 119 comoving kpc/h away with a wind velocity of $\sim$ 100 km/s launched at $z \sim$ 7. Such a dwarf ($M_{\text{UV}} =$ - 20.5) is fainter than the detection limit of the observed example. In a general analysis of galaxy - CIV absorbers, we find galaxies with -20.5 $< M_{\text{UV}} <$ - 18.8 are responsible for the observed metal signatures. In addition, we find no correlation between the mass of the closest galaxy to the absorber and the distance between them, but a weak anti-correlation between the strength of the absorption and the separation of galaxy - absorber pairs.

We propose an innovative method for measuring the neutral hydrogen (HI) content of an optically-selected spectroscopic sample of galaxies through cross-correlation with HI intensity mapping measurements. We show that the HI-galaxy cross-power spectrum contains an additive shot noise term which scales with the average HI brightness temperature of the optically-selected galaxies, allowing constraints to be placed on the average HI mass per galaxy. This approach can estimate the HI content of populations too faint to directly observe through their 21cm emission over a wide range of redshifts. This cross-correlation, as a function of optical luminosity or colour, can be used to derive HI-scaling relations. We demonstrate that this signal will be detectable by cross-correlating upcoming Australian SKA Pathfinder (ASKAP) observations with existing optically-selected samples. We also use semi-analytic simulations to verify that the HI mass can be successfully recovered by our technique in the range M_HI > 10^8 M_solar, in a manner independent of the underlying power spectrum shape. We conclude that this method is a powerful tool to study galaxy evolution, which only requires a single intensity mapping dataset to infer complementary HI gas information from existing optical and infra-red observations.

We investigate how the hydrostatic suppression of baryonic accretion affects the growth rate of dark matter halos during the Epoch of Reionization. By comparing halo properties in a simplistic hydrodynamic simulation in which gas only cools adiabatically, with its collisionless equivalent, we find that halo growth is slowed as hydrostatic forces prevent gas from collapsing. In our simulations, at the high redshifts relevant for reionization (between ${\sim}6$ and ${\sim}11$), halos that host dwarf galaxies ($\lesssim 10^{9} \mathrm{M_\odot}$) can be reduced by up to a factor of 2 in mass due to the hydrostatic pressure of baryons. Consequently, the inclusion of baryonic effects reduces the amplitude of the low mass tail of the halo mass function by factors of 2 to 4. In addition, we find that the fraction of baryons in dark matter halos hosting dwarf galaxies at high redshift never exceeds ${\sim}90\%$ of the cosmic baryon fraction. When implementing baryonic processes, including cooling, star formation, supernova feedback and reionization, the suppression effects become more significant with further reductions of ${\sim}30\%$ to 60\%. Although convergence tests suggest that the suppression may become weaker in higher resolution simulations, this suppressed growth will be important for semi-analytic models of galaxy formation, in which the halo mass inherited from an underlying N-body simulation directly determines galaxy properties. Based on the adiabatic simulation, we provide tables to account for these effects in N-body simulations, and present a modification of the halo mass function along with explanatory analytic calculations.

Using Hubble data, including new grism spectra, Oesch et al. recently identified GN-z11, an $M_\textrm{UV}$=-21.1 galaxy at $z$=11.1 (just 400Myr after the big bang). With an estimated stellar mass of $\sim$10$^9$M$_{\odot}$, this galaxy is surprisingly bright and massive, raising questions as to how such an extreme object could form so early in the Universe. Using \Meraxes{}, a semi-analytic galaxy-formation model developed as part of the Dark-ages Reionization And Galaxy-formation Observables from Numerical Simulations (DRAGONS) programme, we investigate the potential formation mechanisms and eventual fate of GN-z11. The volume of our simulation is comparable to that of the discovery observations and possesses two analogue galaxies of similar luminosity to this remarkably bright system. Existing in the two most massive subhaloes at $z$=11.1 ($M_\textrm{vir}$=1.4$\times 10^{11}$M$_{\odot}$ and 6.7$\times 10^{10}$M$_{\odot}$), our model analogues show excellent agreement with all available observationally derived properties of GN-z11. Although they are relatively rare outliers from the full galaxy population at high-$z$, they are no longer the most massive or brightest systems by $z$=5. Furthermore, we find that both objects possess relatively smooth, but extremely rapid mass growth histories with consistently high star formation rates and UV luminosities at $z{>}11$, indicating that their brightness is not a transient, merger-driven feature. Our model results suggest that future wide-field surveys with the \textit{James Webb Space Telescope} may be able to detect the progenitors of GN-z11 analogues out to $z{\sim}$13--14, pushing the frontiers of galaxy-formation observations to the early phases of cosmic reionization and providing a valuable glimpse of the first galaxies to reionize the Universe on large scales.

We investigate the properties of the stellar populations of model galaxies as a function of galaxy evolutionary history and angular momentum content. We use the new semi-analytic model presented in Tonini et al. (2016). This new model follows the angular momentum evolution of gas and stars, providing the base for a new star formation recipe, and treatment of the effects of mergers that depends on the central galaxy dynamical structure. We find that the new recipes have the effect of boosting the efficiency of the baryonic cycle in producing and recycling metals, as well as preventing minor mergers from diluting the metallicity of bulges and ellipticals. The model reproduces the stellar mass - stellar metallicity relation for galaxies above 1e10 solar masses, including Brightest Cluster Galaxies. Model disks, galaxies dominated by instability-driven components, and merger-driven objects each stem from different evolutionary channels. These model galaxies therefore occupy different loci in the galaxy mass-size relation, which we find to be in accord with the Atlas 3D classification of disk galaxies, fast rotators and slow rotators. We find that the stellar populations' properties depend on the galaxy evolutionary type, with more evolved stellar populations being part of systems that have lost or dissipated more angular momentum during their assembly history.

We investigate high-redshift galaxy sizes using a semi-analytic model constructed for the Dark-ages Reionization And Galaxy-formation Observables from Numerical Simulation project. Our fiducial model, including strong feedback from supernovae and photoionization background, accurately reproduces the evolution of the stellar mass function and UV luminosity function. Using this model, we study the size--luminosity relation of galaxies and find that the effective radius scales with UV luminosity as $R_\mathrm{e}\propto L^{0.25}$ at $z{\sim}5$--$9$. We show that recently discovered very luminous galaxies at $z{\sim}7$ (Bowler et al. 2016) and $z{\sim}11$ (Oesch et al. 2016) lie on our predicted size--luminosity relations. We find that a significant fraction of galaxies at $z>8$ will not be resolved by JWST, but GMT will have the ability to resolve all galaxies in haloes above the atomic cooling limit. We show that our fiducial model successfully reproduces the redshift evolution of average galaxy sizes at $z>5$. We also explore galaxy sizes in models without supernova feedback. The no-supernova feedback models produce galaxy sizes that are smaller than observations. We therefore confirm that supernova feedback plays an important role in determining the size--luminosity relation of galaxies and its redshift evolution during reionization.

The formation and evolution of galaxies with low neutral atomic hydrogen (HI) masses, M$_{\rm HI}$$<$10$^{8}h^{-2}$M$_{\odot}$, are affected by host dark matter halo mass and photoionisation feedback from the UV background after the end of reionization. We study how the physical processes governing the formation of galaxies with low HI mass are imprinted on the distribution of neutral hydrogen in the Universe using the hierarchical galaxy formation model, GALFORM. We calculate the effect on the correlation function of changing the HI mass detection threshold at redshifts $0 \le z \le 0.5$. We parameterize the clustering as $\xi(r)=(r/r_{0})^{-\gamma}$ and we find that including galaxies with M$_{\rm HI}$$<$10$^{8}h^{-2}$M$_{\odot}$ increases the clustering amplitude $r_{0}$ and slope $\gamma$ compared to samples of higher HI masses. This is due to these galaxies with low HI masses typically being hosted by haloes with masses greater than 10$^{12}{h}^{-1}$M$_{\odot}$, and is in contrast to optically selected surveys for which the inclusion of faint, blue galaxies lowers the clustering amplitude. We show the HI mass function for different host dark matter halo masses and galaxy types (central or satellite) to interpret the values of $r_{0}$ and $\gamma$ of the clustering of HI-selected galaxies. We also predict the contribution of low HI mass galaxies to the 21cm intensity mapping signal. We calculate that a dark matter halo mass resolution better than $\sim$10$^{10}{h}^{-1}$M$_{\odot}$ at redshifts higher than 0.5 is required in order to predict converged 21cm brightness temperature fluctuations.

The Detection of redshifted 21 cm emission from the epoch of reionization (EoR) is a challenging task owing to strong foregrounds that dominate the signal. In this paper, we propose a general method, based on the delay spectrum approach, to extract HI power spectra that is applicable to tracking observations using an imaging radio interferometer (Delay Spectrum with Imaging Arrays (DSIA)). Our method is based on modelling the HI signal taking into account the impact of wide field effects such as the $w$-term which are then used as appropriate weights in cross-correlating the measured visibilities. Our method is applicable to any radio interferometer that tracks a phase center and could be utilized for arrays such as MWA, LOFAR, GMRT, PAPER and HERA. In the literature the delay spectrum approach has been implemented for near-redundant baselines using drift scan observations. In this paper we explore the scheme for non-redundant tracking arrays, and this is the first application of delay spectrum methodology to such data to extract the HI signal. We analyze 3 hours of MWA tracking data on the EoR1 field. We present both 2-dimensional ($k_\parallel,k_\perp$) and 1-dimensional (k) power spectra from the analysis. Our results are in agreement with the findings of other pipelines developed to analyse the MWA EoR data.

The Murchison Widefield Array (MWA) has collected hundreds of hours of Epoch of Reionization (EoR) data and now faces the challenge of overcoming foreground and systematic contamination to reduce the data to a cosmological measurement. We introduce several novel analysis techniques such as cable reflection calibration, hyper-resolution gridding kernels, diffuse foreground model subtraction, and quality control methods. Each change to the analysis pipeline is tested against a two dimensional power spectrum figure of merit to demonstrate improvement. We incorporate the new techniques into a deep integration of 32 hours of MWA data. This data set is used to place a systematic-limited upper limit on the cosmological power spectrum of $\Delta^2 \leq 2.7 \times 10^4$ mK$^2$ at $k=0.27$ h~Mpc$^{-1}$ and $z=7.1$, consistent with other published limits, and a modest improvement (factor of 1.4) over previous MWA results. From this deep analysis we have identified a list of improvements to be made to our EoR data analysis strategies. These improvements will be implemented in the future and detailed in upcoming publications.

We introduce Meraxes, a new, purpose-built semi-analytic galaxy formation model designed for studying galaxy growth during reionization. Meraxes is the first model of its type to include a temporally and spatially coupled treatment of reionization and is built upon a custom (100 Mpc)$^3$ N-body simulation with high temporal and mass resolution, allowing us to resolve the galaxy and star formation physics relevant to early galaxy formation. Our fiducial model with supernova feedback reproduces the observed optical depth to electron scattering and evolution of the galaxy stellar mass function between $z$=5 and 7, predicting that a broad range of halo masses contribute to reionization. Using a constant escape fraction and global recombination rate, our model is unable to simultaneously match the observed ionizing emissivity at $z{\lesssim}6$. However, the use of an evolving escape fraction of 0.05-0.1 at $z{\sim}6$, increasing towards higher redshift, is able to satisfy these three constraints. We also demonstrate that photoionization suppression of low mass galaxy formation during reionization has only a small effect on the ionization history of the inter-galactic medium. This lack of `self-regulation' arises due to the already efficient quenching of star formation by supernova feedback. It is only in models with gas supply-limited star formation that reionization feedback is effective at regulating galaxy growth. We similarly find that reionization has only a small effect on the stellar mass function, with no observationally detectable imprint at $M_{\rm *}{>}10^{7.5}\,{\rm M_{{\odot}}}$. However, patchy reionization has significant effects on individual galaxy masses, with variations of factors of 2-3 at $z$=5 that correlate with environment.

We use the Dark-ages, Reionization And Galaxy-formation Observables from Numerical Simulations (DRAGONS) framework to investigate the effect of galaxy-formation physics on the morphology and statistics of ionized hydrogen (HII) regions during the Epoch of Reioinization (EoR). DRAGONS self-consistently couples a semi-analytic galaxy-formation model with the inhomogeneous ionizing UV background, and can therefore be used to study the dependence of morphology and statistics of reionization on feedback phenomena of the ionizing source galaxy population. Changes in galaxy-formation physics modify the sizes of HII regions and the amplitude and shape of 21-cm power spectra. Of the galaxy physics investigated, we find that supernova feedback plays the most important role in reionization, with HII regions up to $\approx 20$ per cent smaller and a fractional difference in the amplitude of power spectra of up to $\approx 17$ per cent at fixed ionized fraction in the absence of this feedback. We compare our galaxy-formation-based reionization models with past calculations that assume constant stellar-to-halo mass ratios and find that with the correct choice of minimum halo mass, such models can mimic the predicted reionization morphology. Reionization morphology at fixed neutral fraction is therefore not uniquely determined by the details of galaxy formation, but is sensitive to the mass of the haloes hosting the bulk of the ionizing sources. Simple EoR parametrizations are therefore accurate predictors of reionization statistics. However, a complete understanding of reionization using future 21-cm observations will require interpretation with realistic galaxy-formation models, in combination with other observations.

We present deep polarimetric observations at 154 MHz with the Murchison Widefield Array (MWA), covering 625 deg^2 centered on RA=0 h, Dec=-27 deg. The sensitivity available in our deep observations allows an in-band, frequency-dependent analysis of polarized structure for the first time at long wavelengths. Our analysis suggests that the polarized structures are dominated by intrinsic emission but may also have a foreground Faraday screen component. At these wavelengths, the compactness of the MWA baseline distribution provides excellent snapshot sensitivity to large-scale structure. The observations are sensitive to diffuse polarized emission at ~54' resolution with a sensitivity of 5.9 mJy beam^-1 and compact polarized sources at ~2.4' resolution with a sensitivity of 2.3 mJy beam^-1 for a subset (400 deg^2) of this field. The sensitivity allows the effect of ionospheric Faraday rotation to be spatially and temporally measured directly from the diffuse polarized background. Our observations reveal large-scale structures (~1 deg - 8 deg in extent) in linear polarization clearly detectable in ~2 minute snapshots, which would remain undetectable by interferometers with minimum baseline lengths >110 m at 154 MHz. The brightness temperature of these structures is on average 4 K in polarized intensity, peaking at 11 K. Rotation measure synthesis reveals that the structures have Faraday depths ranging from -2 rad m^-2 to 10 rad m^-2 with a large fraction peaking at ~+1 rad m^-2. We estimate a distance of 51+/-20 pc to the polarized emission based on measurements of the in-field pulsar J2330-2005. We detect four extragalactic linearly polarized point sources within the field in our compact source survey. Based on the known polarized source population at 1.4 GHz and non-detections at 154 MHz, we estimate an upper limit on the depolarization ratio of 0.08 from 1.4 GHz to 154 MHz.

Detection of the Epoch of Reionization HI signal requires a precise understanding of the intervening galaxies and AGN, both for instrumental calibration and foreground removal. We present a catalogue of 7394 extragalactic sources at 182 MHz detected in the RA=0 field of the Murchison Widefield Array Epoch of Reionization observation programme. Motivated by unprecedented requirements for precision and reliability we develop new methods for source finding and selection. We apply machine learning methods to self-consistently classify the relative reliability of 9490 source candidates. A subset of 7466 are selected based on reliability class and signal-to-noise ratio criteria. These are statistically cross-matched to four other radio surveys using both position and flux density information. We find 7369 sources to have confident matches, including 90 partially resolved sources that split into a total of 192 sub-components. An additional 25 unmatched sources are included as new radio detections. The catalogue sources have a median spectral index of -0.85. Spectral flattening is seen toward lower frequencies with a median of -0.71 predicted at 182 MHz. The astrometric error is 7 arcsec. compared to a 2.3 arcmin. beam FWHM. The resulting catalogue covers approximately 1400 sq. deg. and is complete to approximately 80 mJy within half beam power. This provides the most reliable discrete source sky model available to date in the MWA EoR0 field for precision foreground subtraction.

We investigate the clustering of Lyman-break galaxies (LBGs) at $z\sim4$. Using the hierarchical galaxy formation model GALFORM, we predict, for the first time using a semi-analytical model with feedback from active galactic nuclei (AGN), the angular correlation function (ACF) of LBGs and find agreement within $3\,\sigma$ with new measurements of the ACF from surveys including the Hubble eXtreme Deep Field (XDF) and CANDELS field. Our simulations confirm the conclusion reached using independent models that although the predicted ACFs reproduce the trend of increased clustering with luminosity, the dependence is less strong than observed. We find that for the detection limits of the XDF field central LBGs at $z\sim 4$ predominantly reside in haloes of mass $\sim 10^{11}-10^{12}h^{-1}M_{\rm \odot}$ and that satellites reside in larger haloes of mass $\sim 10^{12}-10^{13}h^{-1}M_{\rm \odot}$. The model predicts fewer bright satellite LBGs at $z\sim4$ than is inferred from measurements of the ACF at small scales. By analysing the halo occupation distribution (HOD) predicted by the model, we find evidence that AGN feedback affects the HOD of central LBGs in massive haloes. This is a new high-redshift test of this important feedback mechanism. We investigate the effect of photometric errors in the observations on the ACF predictions. We find that the observational uncertainty in the galaxy luminosity reduces the clustering amplitude and that this effect increases towards faint galaxies, particularly on small scales. To compare properties of model with observed LBGs this uncertainty must be considered.

We present the 21 cm power spectrum analysis approach of the Murchison Widefield Array Epoch of Reionization project. In this paper, we compare the outputs of multiple pipelines for the purpose of validating statistical limits cosmological hydrogen at redshifts between 6 and 12. Multiple, independent, data calibration and reduction pipelines are used to make power spectrum limits on a fiducial night of data. Comparing the outputs of imaging and power spectrum stages highlights differences in calibration, foreground subtraction and power spectrum calculation. The power spectra found using these different methods span a space defined by the various tradeoffs between speed, accuracy, and systematic control. Lessons learned from comparing the pipelines range from the algorithmic to the prosaically mundane; all demonstrate the many pitfalls of neglecting reproducibility. We briefly discuss the way these different methods attempt to handle the question of evaluating a significant detection in the presence of foregrounds.

We present first results from radio observations with the Murchison Widefield Array seeking to constrain the power spectrum of 21 cm brightness temperature fluctuations between the redshifts of 11.6 and 17.9 (113 and 75 MHz). Three hours of observations were conducted over two nights with significantly different levels of ionospheric activity. We use these data to assess the impact of systematic errors at low frequency, including the ionosphere and radio-frequency interference, on a power spectrum measurement. We find that after the 1-3 hours of integration presented here, our measurements at the Murchison Radio Observatory are not limited by RFI, even within the FM band, and that the ionosphere does not appear to affect the level of power in the modes that we expect to be sensitive to cosmology. Power spectrum detections, inconsistent with noise, due to fine spectral structure imprinted on the foregrounds by reflections in the signal-chain, occupy the spatial Fourier modes where we would otherwise be most sensitive to the cosmological signal. We are able to reduce this contamination using calibration solutions derived from autocorrelations so that we achieve an sensitivity of $10^4$ mK on comoving scales $k\lesssim 0.5 h$Mpc$^{-1}$. This represents the first upper limits on the $21$ cm power spectrum fluctuations at redshifts $12\lesssim z \lesssim 18$ but is still limited by calibration systematics. While calibration improvements may allow us to further remove this contamination, our results emphasize that future experiments should consider carefully the existence of and their ability to calibrate out any spectral structure within the EoR window.

The relation between the Star Formation Rate (SFR) and stellar mass (${\rm M}_{\star}$) of galaxies represents a fundamental constraint on galaxy formation and has been studied extensively both in observations and cosmological simulations. However, the observed amplitude has not been successfully reproduced in simulations, indicating either that the halo accretion history and baryonic physics are poorly modeled or that observations contain biases. We examine the evolution of the SFR$-{\rm M}_{\star}$ relation of $z\sim1-4 $ galaxies and display the inconsistency between observed relations that are obtained using different techniques. We employ cosmological hydrodynamic simulations and compare these with a range of observed SFR$-{\rm M}_{\star}$ relations. We find that numerical results are consistent with observations that use Spectral Energy Distribution (SED) techniques to estimate star formation rates and dust corrections. On the contrary, simulations are not able to reproduce results that were obtained by combining only UV and IR luminosities. These imply SFRs at a fixed stellar mass that are larger almost by a factor of 5 than those of SED measurements for $z \sim1.5-4$. Furthermore, we find remarkable agreement between the numerical results from various authors who have employed different cosmological codes and run simulations with different resolutions. This is interesting for two reasons. A) simulations can produce realistic populations of galaxies within representative cosmological volumes even at relatively modest resolutions. B) It is likely that current numerical codes that rely on similar subgrid multiphase Inter-Stellar Medium (ISM) models and are tuned to reproduce statistical properties of galaxies, produce similar results for the SFR$-{\rm M}_{\star}$ relation by construction, regardless of resolution, box size and, to some extent, the adopted feedback prescriptions.

In this paper we present calculations of the UV luminosity function from the Dark-ages Reionization And Galaxy-formation Observables from Numerical Simulations (DRAGONS) project, which combines N-body, semi-analytic and semi-numerical modelling designed to study galaxy formation during the Epoch of Reionization. Using galaxy formation physics including supernova feedback, the model naturally reproduces the UV LFs for high-redshift star-forming galaxies from $z{\sim}5$ through to $z{\sim}10$. We investigate the luminosity--star formation rate (SFR) relation, finding that variable SFR histories of galaxies result in a scatter around the median relation of $0.1$--$0.3$ dex depending on UV luminosity. We find close agreement between the model and observationally derived SFR functions. We use our calculated luminosities to investigate the luminosity function below current detection limits, and the ionizing photon budget for reionization. We predict that the slope of the UV LF remains steep below current detection limits and becomes flat at $M_\mathrm{UV}{\gtrsim}{-14}$. We find that $48$ ($17$) per cent of the total UV flux at $z{\sim}6$ ($10$) has been detected above an observational limit of $M_\mathrm{UV}{\sim}{-17}$, and that galaxies fainter than $M_\mathrm{UV}{\sim}{-17}$ are the main source of ionizing photons for reionization. We investigate the luminosity--stellar mass relation, and find a correlation for galaxies with $M_\mathrm{UV}{<}{-14}$ that has the form $M_*{\propto}10^{-0.47M_\mathrm{UV}}$, in good agreement with observations, but which flattens for fainter galaxies. We determine the luminosity--halo mass relation to be $M_\mathrm{vir}{\propto}10^{-0.35M_\mathrm{UV}}$, finding that galaxies with $M_\mathrm{UV}{=}{-20}$ reside in host dark matter haloes of $10^{11.0\pm 0.1}\mathrm{M_\odot}$ at $z{\sim}6$, and that this mass decreases towards high redshift.

We present a theoretical model for the evolution of mass, angular momentum and size of galaxy disks and bulges, and we implement it into the semi-analytic galaxy formation code SAGE. The model follows both secular and violent evolutionary channels, including smooth accretion, disk instabilities, minor and major mergers. We find that the combination of our recipe with hierarchical clustering produces two distinct populations of bulges: merger-driven bulges, akin to classical bulges and ellipticals, and instability-driven bulges, akin to secular (or pseudo-)bulges. The model mostly reproduces the mass-size relation of gaseous and stellar disks, the evolution of the mass-size relation of ellipticals, the Faber-Jackson relation, and the magnitude-colour diagram of classical and secular bulges. The model predicts only a small overlap of merger-driven and instability-driven components in the same galaxy, and predicts different bulge types as a function of galaxy mass and disk fraction. Bulge type also affects the star formation rate and colour at a given luminosity. The model predicts a population of merger-driven red ellipticals that dominate both the low-mass and high-mass ends of the galaxy population, and span all dynamical ages; merger-driven bulges in disk galaxies are dynamically old and do not interfere with subsequent evolution of the star-forming component. Instability-driven bulges dominate the population at intermediate galaxy masses, especially thriving in massive disks. The model green valley is exclusively populated by instability-driven bulge hosts. Through the present implementation the mass accretion history is perceivable in the galaxy structure, morphology and colours.

Experiments that pursue detection of signals from the Epoch of Reionization (EoR) are relying on spectral smoothness of source spectra at low frequencies. This article empirically explores the effect of foreground spectra on EoR experiments by measuring high-resolution full-polarization spectra for the 586 brightest unresolved sources in one of the MWA EoR fields using 45 h of observation. A novel peeling scheme is used to subtract 2500 sources from the visibilities with ionospheric and beam corrections, resulting in the deepest, confusion-limited MWA image so far. The resulting spectra are found to be affected by instrumental effects, which limit the constraints that can be set on source-intrinsic spectral structure. The sensitivity and power-spectrum of the spectra are analysed, and it is found that the spectra of residuals are dominated by PSF sidelobes from nearby undeconvolved sources. We release a catalogue describing the spectral parameters for each measured source.

In this paper we present observations, simulations, and analysis demonstrating the direct connection between the location of foreground emission on the sky and its location in cosmological power spectra from interferometric redshifted 21 cm experiments. We begin with a heuristic formalism for understanding the mapping of sky coordinates into the cylindrically averaged power spectra measurements used by 21 cm experiments, with a focus on the effects of the instrument beam response and the associated sidelobes. We then demonstrate this mapping by analyzing power spectra with both simulated and observed data from the Murchison Widefield Array. We find that removing a foreground model which includes sources in both the main field-of-view and the first sidelobes reduces the contamination in high k_parallel modes by several percent relative to a model which only includes sources in the main field-of-view, with the completeness of the foreground model setting the principal limitation on the amount of power removed. While small, a percent-level amount of foreground power is in itself more than enough to prevent recovery of any EoR signal from these modes. This result demonstrates that foreground subtraction for redshifted 21 cm experiments is truly a wide-field problem, and algorithms and simulations must extend beyond the main instrument field-of-view to potentially recover the full 21 cm power spectrum.

Detection of the cosmological neutral hydrogen signal from the Epoch of Reionization, and estimation of its basic physical parameters, is the principal scientific aim of many current low-frequency radio telescopes. Here we describe the Cosmological HI Power Spectrum Estimator (CHIPS), an algorithm developed and implemented with data from the Murchison Widefield Array (MWA), to compute the two-dimensional and spherically-averaged power spectrum of brightness temperature fluctuations. The principal motivations for CHIPS are the application of realistic instrumental and foreground models to form the optimal estimator, thereby maximising the likelihood of unbiased signal estimation, and allowing a full covariant understanding of the outputs. CHIPS employs an inverse-covariance weighting of the data through the maximum likelihood estimator, thereby allowing use of the full parameter space for signal estimation ("foreground suppression"). We describe the motivation for the algorithm, implementation, application to real and simulated data, and early outputs. Upon application to a set of 3 hours of data, we set a 2$\sigma$ upper limit on the EoR dimensionless power at $k=0.05$~h.Mpc$^{-1}$ of $\Delta_k^2<7.6\times{10^4}$~mK$^2$ in the redshift range $z=[6.2-6.6]$, consistent with previous estimates.

Intensity mapping of the neutral hydrogen (HI) is a new observational tool that can be used to efficiently map the large-scale structure of the Universe over wide redshift ranges. The power spectrum of the intensity maps contains cosmological information on the matter distribution and probes galaxy evolution by tracing the HI content of galaxies at different redshifts and the scale-dependence of HI clustering. The cross-correlation of intensity maps with galaxy surveys is a robust measure of the power spectrum which diminishes systematics caused by instrumental effects and foreground removal. We examine the cross-correlation signature at redshift z=0.9 using a variant of the semi-analytical galaxy formation model SAGE (Croton et al. 2016) applied to the Millennium simulation in order to model the HI gas of galaxies as well as their optical magnitudes based on their star-formation history. We determine the clustering of the cross-correlation power for different types of galaxies determined by their colours, acting as a proxy for their star-formation activity. We find that the cross-correlation coefficient for red quiescent galaxies falls off more quickly on smaller scales k>0.2h/Mpc than for blue star-forming galaxies. Additionally, we create a mock catalogue of highly star-forming galaxies using a selection function to mimic the WiggleZ survey, and use this to predict existing and future cross-correlation measurements of the GBT and Parkes telescope. We find that the cross-power of highly star-forming galaxies shows a higher clustering on small scales than any other galaxy type and that this significantly alters the power spectrum shape on scales k>0.2h/Mpc. We show that the cross-correlation coefficient is not negligible when interpreting the cosmological cross-power spectrum. On the other hand, it contains information about the HI content of the optically selected galaxies.

We use high resolution N-Body simulations to study the concentration and spin parameters of dark matter haloes in the mass range $10^8\, {\rm M}_{\odot}\, h^{-1} < {\rm M} < 10^{11}\, {\rm M}_{\odot}\, h^{-1}$ and redshifts $5{<}z{<}10$, corresponding to the haloes of galaxies thought to be responsible for reionization. We build a sub-sample of equilibrium haloes and contrast their properties to the full population that also includes unrelaxed systems. Concentrations are calculated by fitting both NFW and Einasto profiles to the spherically-averaged density profiles of individual haloes. After removing haloes that are out-of-equilibrium, we find a $z{>}5$ concentration$-$mass ($c(M)$) relation that is almost flat and well described by a simple power-law for both NFW and Einasto fits. The intrinsic scatter around the mean relation is $\Delta c_{\rm{vir}}{\sim1}$ (or 20 per cent) at $z=5$. We also find that the analytic model proposed by Ludlow et al. reproduces the mass and redshift-dependence of halo concentrations. Our best-fit Einasto shape parameter, $\alpha$, depends on peak height, $\nu$, in a manner that is accurately described by $\alpha {=}0.0070\nu^2{+}0.1839$. The distribution of the spin parameter, $\lambda$, has a weak dependence on equilibrium state; $\lambda$ peaks at roughly ${\sim}0.033$ for our relaxed sample, and at ${\sim}0.04$ for the full population. The spin--virial mass relation has a mild negative correlation at high redshift.

Redshifted 21cm measurements of the structure of ionised regions that grow during reionization promise to provide a new probe of early galaxy and structure formation. One of the challenges of modelling reionization is to account both for the sub-halo scale physics of galaxy formation and the regions of ionization on scales that are many orders of magnitude larger. To bridge this gap we first calculate the statistical relationship between ionizing luminosity and Mpc-scale overdensity using detailed models of galaxy formation computed using relatively small volume - ($\sim$100Mpc/$h$)$^{3}$, high resolution dark matter simulations. We then use a Monte-Carlo technique to apply this relationship to reionization of the intergalactic medium within large volume dark matter simulations - ($>$1Gpc/$h$)$^{3}$. The resulting simulations can be used to address the contribution of very large scale clustering of galaxies to the structure of reionization, and show that volumes larger than 500Mpc/$h$ are required to probe the largest reionization features mid-way through reionization. As an example application of our technique, we demonstrate that the predicted 21cm power spectrum amplitude and gradient could be used to determine the importance of supernovae feedback for early galaxy formation.

Gravitational waves are expected to be radiated by supermassive black hole binaries formed during galaxy mergers. A stochastic superposition of gravitational waves from all such binary systems will modulate the arrival times of pulses from radio pulsars. Using observations of millisecond pulsars obtained with the Parkes radio telescope, we constrain the characteristic amplitude of this background, $A_{\rm c,yr}$, to be < $1.0\times10^{-15}$ with 95% confidence. This limit excludes predicted ranges for $A_{\rm c,yr}$ from current models with 91-99.7% probability. We conclude that binary evolution is either stalled or dramatically accelerated by galactic-center environments, and that higher-cadence and shorter-wavelength observations would result in an increased sensitivity to gravitational waves.

We present a detailed analysis of an individual case of gravitational lensing of a $z\sim8$ Lyman-Break galaxy (LBG) in a blank field, identified in Hubble Space Telescope imaging obtained as part of the Brightest of Reionizing Galaxies survey. To investigate the close proximity of the bright ($m_{AB}=25.8$) $Y_{098}$-dropout to a small group of foreground galaxies, we obtained deep spectroscopy of the dropout and two foreground galaxies using VLT/X-Shooter. We detect H-$\alpha$, H-$\beta$, [OIII] and [OII] emission in the brightest two foreground galaxies (unresolved at the natural seeing of $0.8$ arcsec), placing the pair at $z=1.327$. We can rule out emission lines contributing all of the observed broadband flux in $H_{160}$ band at $70\sigma$, allowing us to exclude the $z\sim8$ candidate as a low redshift interloper with broadband photometry dominated by strong emission lines. The foreground galaxy pair lies at the peak of the luminosity, redshift and separation distributions for deflectors of strongly lensed $z\sim8$ objects, and we make a marginal detection of a demagnified secondary image in the deepest ($J_{125}$) filter. We show that the configuration can be accurately modelled by a singular isothermal ellipsoidal deflector and a S\'{e}rsic source magnified by a factor of $\mu=4.3\pm0.2$. The reconstructed source in the best-fitting model is consistent with luminosities and morphologies of $z\sim8$ LBGs in the literature. The lens model yields a group mass of $9.62\pm0.31\times10^{11} M_{\odot}$ and a stellar mass-to-light ratio for the brightest deflector galaxy of $M_{\star}/L_{B}=2.3^{+0.8}_{-0.6} M_{\odot}/L_{\odot}$ within its effective radius. The foreground galaxies' redshifts would make this one of the few strong lensing deflectors discovered at $z>1$.

We explore the galaxy formation physics governing the low mass end of the HI mass function in the local Universe. Specifically, we predict the effects on the HI mass function of varying i) the strength of photoionisation feedback and the redshift of the end of the epoch of reionization, ii) the cosmology, iii) the supernovae feedback prescription, and iv) the efficiency of star formation. We find that the shape of the low-mass end of the HI mass function is most affected by the critical halo mass below which galaxy formation is suppressed by photoionisation heating of the intergalactic medium. We model the redshift dependence of this critical dark matter halo mass by requiring a match to the low-mass end of the HI mass function. The best fitting critical dark matter halo mass decreases as redshift increases in this model, corresponding to a circular velocity of $\sim 50 \, {\rm km \,s}^{-1}$ at $z=0$, $\sim 30 \, {\rm km\, s}^{-1}$ at $z \sim 1$ and $\sim 12 \, {\rm km \, s}^{-1}$ at $z=6$. We find that an evolving critical halo mass is required to explain both the shape and abundance of galaxies in the HI mass function below $M_{\rm HI} \sim 10^{8} h^{-2} {\rm M_{\odot}}$. The model makes specific predictions for the clustering strength of HI-selected galaxies with HI masses > $10^{6} h^{-2} {\rm M_{\odot}}$ and $> 10^{7} h^{-2} {\rm M_{\odot}}$ and for the relation between the HI and stellar mass contents of galaxies which will be testable with upcoming surveys with the Square Kilometre Array and its pathfinders. We conclude that measurements of the HI mass function at $z \ge 0$ will lead to an improvement in our understanding of the net effect of photoionisation feedback on galaxy formation and evolution.

We present a semi-analytic, physically motivated model for dark matter halo concentration as a function of halo mass and redshift. The semi-analytic model combines an analytic model for the halo mass accretion history (MAH), based on extended Press Schechter (EPS) theory, with an empirical relation between concentration and formation time obtained through fits to the results of numerical simulations. Because the semi-analytic model is based on EPS theory, it can be applied to wide ranges in mass, redshift and cosmology. The resulting concentration-mass (c-M) relations are found to agree with the simulations, and because the model applies only to relaxed halos, they do not exhibit the upturn at high masses or high redshifts found by some recent works. We predict a change of slope in the z=0 c-M relation at a mass scale of $10^{11}\rm{M}_{\odot}$. We find that this is due to the change in the functional form of the halo MAH, which goes from being dominated by an exponential (for high-mass halos) to a power-law (for low-mass halos). During the latter phase, the core radius remains approximately constant, and the concentration grows due to the drop of the background density. We also analyse how the c-M relation predicted by this work affects the power produced by dark matter annihilation, finding that at z = 0 the power is two orders of magnitude lower than that obtained from extrapolating best-fitting c-M relations. We provide fits to the c-M relations as well as numerical routines to compute concentrations and MAHs.

We explore the relation between the structure and mass accretion histories of dark matter halos using a suite of cosmological simulations. We confirm that the formation time, defined as the time when the virial mass of the main progenitor equals the mass enclosed within the scale radius, correlates strongly with concentration. We provide a semi-analytic model for halo mass history that combines analytic relations with fits to simulations. This model has the functional form, $M(z) = M_{0}(1+z)^{\alpha}e^{\beta z}$, where the parameters $\alpha$ and $\beta$ are directly correlated with concentration. We then combine this model for the halo mass history with the analytic relations between $\alpha$, $\beta$ and the linear power spectrum derived by Correa et al. (2014) to establish the physical link between halo concentration and the initial density perturbation field. Finally, we provide fitting formulas for the halo mass history as well as numerical routines, we derive the accretion rate as a function of halo mass, and we demonstrate how the halo mass history depends on cosmology and the adopted definition of halo mass.

Understanding the universal accretion history of dark matter halos is the first step towards determining the origin of their structure. We use the extended Press-Schechter formalism to derive the halo mass accretion history from the growth rate of initial density perturbations. We show that the halo mass history is well described by an exponential function of redshift in the high-redshift regime. However, in the low-redshift regime the mass history follows a power law because the growth of density perturbations is halted in the dark energy dominated era due to the accelerated expansion of the Universe. We provide an analytic model that follows the expression $M(z) = M_{0}(1+z)^{af(M_{0})}e^{-f(M_{0})z}$, where $M_{0} = M(z=0)$, $a$ depends on cosmology and $f(M_{0})$ depends only on the linear matter power spectrum. The analytic model does not rely on calibration against numerical simulations and is suitable for any cosmology. We compare our model with the latest empirical models for the mass accretion history in the literature and find very good agreement. We provide numerical routines for the model online.

We confirm our recent prediction of the "pitchfork" foreground signature in power spectra of high-redshift 21 cm measurements where the interferometer is sensitive to large-scale structure on all baselines. This is due to the inherent response of a wide-field instrument and is characterized by enhanced power from foreground emission in Fourier modes adjacent to those considered to be the most sensitive to the cosmological H I signal. In our recent paper, many signatures from the simulation that predicted this feature were validated against Murchison Widefield Array (MWA) data, but this key pitchfork signature was close to the noise level. In this paper, we improve the data sensitivity through the coherent averaging of 12 independent snapshots with identical instrument settings and provide the first confirmation of the prediction with a signal-to-noise ratio > 10. This wide-field effect can be mitigated by careful antenna designs that suppress sensitivity near the horizon. Simple models for antenna apertures that have been proposed for future instruments such as the Hydrogen Epoch of Reionization Array and the Square Kilometre Array indicate they should suppress foreground leakage from the pitchfork by ~40 dB relative to the MWA and significantly increase the likelihood of cosmological signal detection in these critical Fourier modes in the three-dimensional power spectrum.

We investigate the impact of energy released from self-annihilating dark matter on heating of gas in the small, high-redshift dark matter halos thought to host the first stars. A SUSY neutralino like particle is implemented as our dark matter candidate. The PYTHIA code is used to model the final, stable particle distributions produced during the annihilation process. We use an analytic treatment in conjunction with the code MEDEA2 to find the energy transfer and subsequent partition into heating, ionizing and Lyman alpha photon components. We consider a number of halo density models, dark matter particle masses and annihilation channels. We find that the injected energy from dark matter exceeds the binding energy of the gas within a $10^5$ - $10^6$ M$_\odot$ halo at redshifts above 20, preventing star formation in early halos in which primordial gas would otherwise cool. Thus we find that DM annihilation could delay the formation of the first galaxies.

Detection of 21~cm emission of HI from the epoch of reionization, at redshifts z>6, is limited primarily by foreground emission. We investigate the signatures of wide-field measurements and an all-sky foreground model using the delay spectrum technique that maps the measurements to foreground object locations through signal delays between antenna pairs. We demonstrate interferometric measurements are inherently sensitive to all scales, including the largest angular scales, owing to the nature of wide-field measurements. These wide-field effects are generic to all observations but antenna shapes impact their amplitudes substantially. A dish-shaped antenna yields the most desirable features from a foreground contamination viewpoint, relative to a dipole or a phased array. Comparing data from recent Murchison Widefield Array observations, we demonstrate that the foreground signatures that have the largest impact on the HI signal arise from power received far away from the primary field of view. We identify diffuse emission near the horizon as a significant contributing factor, even on wide antenna spacings that usually represent structures on small scales. For signals entering through the primary field of view, compact emission dominates the foreground contamination. These two mechanisms imprint a characteristic "pitchfork" signature on the "foreground wedge" in Fourier delay space. Based on these results, we propose that selective down-weighting of data based on antenna spacing and time can mitigate foreground contamination substantially by a factor ~100 with negligible loss of sensitivity.

Detection of Lyman-Break Galaxies (LBGs) at high-redshift can be affected by gravitational lensing induced by foreground deflectors not only in galaxy clusters, but also in blank fields. We quantify the impact of strong magnification in the samples of $B$, $V$, $i$, $z$ $\&$ $Y$ LBGs ($4\lesssim z \lesssim8$) observed in the XDF and GOODS/CANDELS fields, by investigating the proximity of dropouts to foreground objects. We find that $\sim6\%$ of bright LBGs ($m_{H_{160}}<26$) at $z\sim7$ have been strongly lensed ($\mu>2$) by foreground objects. This fraction decreases from $\sim 3.5\%$ at $z\sim6$ to $\sim1.5\%$ at $z\sim4$. Since the observed fraction of strongly lensed galaxies is a function of the shape of the luminosity function (LF), it can be used to derive Schechter parameters, $\alpha$ and $M_{\star}$, independently from galaxy number counts. Our magnification bias analysis yields Schechter-function parameters in close agreement with those determined from galaxy counts albeit with larger uncertainties. Extrapolation of our analysis to $z\gtrsim 8$ suggests that future surveys with JSWT, WFIRST and EUCLID should find excess LBGs at the bright-end, even if there is an intrinsic exponential cutoff of number counts. Finally, we highlight how the magnification bias measurement near the detection limit can be used as probe of the population of galaxies too faint to be detected. Preliminary results using this novel idea suggest that the magnification bias at $M_{UV}\sim -18$ is not as strong as expected if $\alpha\lesssim -1.7$ extends well below the current detection limits in the XDF. At face value this implies a flattening of the LF at $M_{UV}\gtrsim-16.5$. However, selection effects and completeness estimates are difficult to quantify precisely. Thus, we do not rule out a steep LF extending to $M_{UV}\gtrsim -15$.

The Murchison Widefield Array (MWA) is a Square Kilometre Array (SKA) Precursor. The telescope is located at the Murchison Radio--astronomy Observatory (MRO) in Western Australia (WA). The MWA consists of 4096 dipoles arranged into 128 dual polarisation aperture arrays forming a connected element interferometer that cross-correlates signals from all 256 inputs. A hybrid approach to the correlation task is employed, with some processing stages being performed by bespoke hardware, based on Field Programmable Gate Arrays (FPGAs), and others by Graphics Processing Units (GPUs) housed in general purpose rack mounted servers. The correlation capability required is approximately 8 TFLOPS (Tera FLoating point Operations Per Second). The MWA has commenced operations and the correlator is generating 8.3 TB/day of correlation products, that are subsequently transferred 700 km from the MRO to Perth (WA) in real-time for storage and offline processing. In this paper we outline the correlator design, signal path, and processing elements and present the data format for the internal and external interfaces.

We develop a semi-analytic method for assessing the impact of the large-scale IGM temperature fluctuations expected following He${\rm\,{\scriptstyle II}}$ reionization on three-dimensional clustering measurements of the Ly$\alpha$ forest. Our methodology builds upon the existing large volume, mock Ly$\alpha$ forest survey simulations presented by Greig et al. by including a prescription for a spatially inhomogeneous ionizing background, temperature fluctuations induced by patchy He${\rm\,{\scriptstyle II}}$ photoheating and the clustering of quasars. This approach enables us to achieve a dynamic range within our semi-analytic model substantially larger than currently feasible with computationally expensive, fully numerical simulations. The results agree well with existing numerical simulations, with large-scale temperature fluctuations introducing a scale-dependent increase in the spherically averaged 3D Ly$\alpha$ forest power spectrum of up to 20-30 per cent at wavenumbers $k\sim0.02$ Mpc$^{-1}$. Although these large-scale thermal fluctuations will not substantially impact upon the recovery of the baryon acoustic oscillation scale from existing and forthcoming dark energy spectroscopic surveys, any complete forward modelling of the broad-band term in the Ly$\alpha$ correlation function will none the less require their inclusion.

The Murchison Widefield Array (MWA) is a new low-frequency interferometric radio telescope built in Western Australia at one of the locations of the future Square Kilometre Array (SKA). We describe the automated radio-frequency interference (RFI) detection strategy implemented for the MWA, which is based on the AOFlagger platform, and present 72-231-MHz RFI statistics from 10 observing nights. RFI detection removes 1.1% of the data. RFI from digital TV (DTV) is observed 3% of the time due to occasional ionospheric or atmospheric propagation. After RFI detection and excision, almost all data can be calibrated and imaged without further RFI mitigation efforts, including observations within the FM and DTV bands. The results are compared to a previously published Low-Frequency Array (LOFAR) RFI survey. The remote location of the MWA results in a substantially cleaner RFI environment compared to LOFAR's radio environment, but adequate detection of RFI is still required before data can be analysed. We include specific recommendations designed to make the SKA more robust to RFI, including: the availability of sufficient computing power for RFI detection; accounting for RFI in the receiver design; a smooth band-pass response; and the capability of RFI detection at high time and frequency resolution (second and kHz-scale respectively).

Large-area sky surveys show that massive galaxies undergo at least one major merger in a Hubble time. Ongoing pulsar timing array (PTA) experiments are aimed at measuring the gravitational wave (GW) emission from binary supermassive black holes (SMBHs) at the centres of galaxy merger remnants. In this paper, using the latest observational estimates for a range of galaxy properties and scaling relations, we predict the amplitude of the GW background generated by the binary SMBH population. We also predict the numbers of individual binary SMBH GW sources. We predict the characteristic strain amplitude of the GW background to lie in the range $5.1\times10^{-16}<A_{\rm yr}<2.4\times10^{-15}$ at a frequency of $(1\,{\rm yr})^{-1}$, with 95% confidence. Higher values within this range, which correspond to the more commonly preferred choice of galaxy merger timescale, will fall within the expected sensitivity ranges of existing PTA projects in the next few years. In contrast, we find that a PTA consisting of at least 100 pulsars observed with next-generation radio telescopes will be required to detect continuous-wave GWs from binary SMBHs. We further suggest that GW memory bursts from coalescing SMBH pairs are not viable sources for PTAs. Both the GW background and individual GW source counts are dominated by binaries formed in mergers between early-type galaxies of masses $\gtrsim5\times10^{10}M_{\odot}$ at redshifts $\lesssim 1.5$. Uncertainties in the galaxy merger timescale and the SMBH mass - galaxy bulge mass relation dominate the uncertainty in our predictions.

Lyman-Break Galaxy (LBG) samples observed during reionization ($z\gtrsim6$) with Hubble Space Telescope's Wide Field Camera 3 are reaching sizes sufficient to characterize their clustering properties. Using a combined catalog from the Hubble eXtreme Deep Field and CANDELS surveys, containing $N=743$ LBG candidates at z>6.5 at a mean redshift of $z=7.2$, we detect a clear clustering signal in the angular correlation function (ACF) at $\sim4\sigma$, corresponding to a real-space correlation length $r_{0}=6.7^{+0.9}_{-1.0}h^{-1}$cMpc. The derived galaxy bias $b=8.6^{+0.9}_{-1.0}$ is that of dark-matter halos of $M=10^{11.1^{+0.2}_{-0.3}}$M$_{\odot}$ at $z=7.2$, and highlights that galaxies below the current detection limit ($M_{AB}\sim-17.7$) are expected in lower-mass halos ($M\sim10^{8}-10^{10.5}$M$_{\odot}$). We compute the ACF of LBGs at $z\sim3.8-z\sim5.9$ in the same surveys. A trend of increasing bias is found from $z=3.8$ ($b\sim3.0$) to $z=7.2$ ($b\sim8.6$), broadly consistent with galaxies at fixed luminosity being hosted in dark-matter halos of similar mass at $4<z<6$, followed by a slight rise in halo masses at $z\sim7$ ($\sim2\sigma$ confidence). Separating the data at the median luminosity of the $z=7.2$ sample ($M_{UV}=-19.4$) shows higher clustering at $z=5.9$ for bright galaxies ($r_{0}=5.5^{+1.4}_{-1.5}h^{-1}$cMpc, $b=6.2^{+1.2}_{-1.5}$) compared to faint galaxies ($r_{0}=1.9^{+1.1}_{-1.0}h^{-1}$cMpc, $b=2.7\pm1.2$) implying a constant mass-to-light ratio $\frac{dlogM}{dlogL}\sim1.2^{+1.8}_{-0.8}$. A similar trend is present in the $z=7.2$ sample with larger uncertainty. Finally, our bias measurements allow us to investigate the fraction of dark-matter halos hosting UV-bright galaxies (the duty-cycle, $\epsilon_{DC}$). At $z=7.2$ values near unity are preferred, which may be explained by the shortened halo assembly time at high-redshift.

We use high-resolution simulations of cosmological volumes to model galaxy formation at high-redshift, with the goal of studying the photon budget for reionization. We demonstrate that galaxy formation models that include a strong, thermally coupled supernovae scheme reproduce current observations of star formation rates and specific star formation rates, both during and after the reionization era. These models produce enough UV photons to sustain reionization at z<8 (z<6) through a significant population of faint, unobserved, galaxies for an assumed escape fraction of 20% (5%). This predicted population is consistent with extrapolation of the faint end of observed UV luminosity functions. We find that heating from a global UV/X-ray background after reionization causes a dip in the total global star formation rate density in galaxies below the current observational threshold. Finally, while the currently observed specific star formation rates are incapable of differentiating between supernovae feedback models, sufficiently deep observations will be able to use this diagnostic in the future to investigate galaxy formation at high redshift.

We assess the effects of super-massive black hole (SMBH) environments on the gravitational-wave (GW) signal from binary SMBHs. To date, searches with pulsar timing arrays for GWs from binary SMBHs, in the frequency band $\sim1-100$\,nHz, include the assumptions that all binaries are circular and evolve only through GW emission. However, dynamical studies have shown that the only way that binary SMBH orbits can decay to separations where GW emission dominates the evolution is through interactions with their environments. We augment an existing galaxy and SMBH formation and evolution model with calculations of binary SMBH evolution in stellar environments, accounting for non-zero binary eccentricities. We find that coupling between binaries and their environments causes the expected GW spectral energy distribution to be reduced with respect to the standard assumption of circular, GW-driven binaries, for frequencies up to $\sim20$\,nHz. Larger eccentricities at binary formation further reduce the signal in this regime. We also find that GW bursts from individual eccentric binary SMBHs are unlikely to be detectable with current pulsar timing arrays. The uncertainties in these predictions are large, owing to observational uncertainty in SMBH-galaxy scaling relations and the galaxy stellar mass function, uncertainty in the nature of binary-environment coupling, and uncertainty in the numbers of the most massive binary SMBHs. We conclude, however, that low-frequency GWs from binary SMBHs may be more difficult to detect with pulsar timing arrays than currently thought.

We present a search for transient and variable radio sources at 154 MHz with the Murchison Widefield Array 32-tile prototype. Fifty-one images were obtained that cover a field of view of 1430 deg^2 centred on Hydra A. The observations were obtained over three days in 2010 March and three days in 2011 April and May. The mean cadence of the observations was 26 minutes and there was additional temporal information on day and year timescales. We explore the variability of a sample of 105 low frequency radio sources within the field. Four bright (S > 6 Jy) candidate variable radio sources were identified that displayed low levels of short timescale variability (26 minutes). We conclude that this variability is likely caused by simplifications in the calibration strategy or ionospheric effects. On the timescale of one year we find two sources that show significant variability. We attribute this variability to either refractive scintillation or intrinsic variability. No radio transients were identified and we place an upper limit on the surface density of sources rho < 7.5 x 10^-5 deg^-2 with flux densities > 5.5 Jy, and characteristic timescales of both 26 minutes and one year.

The formation and growth processes of supermassive black holes (SMBHs) are not well constrained. SMBH population models, however, provide specific predictions for the properties of the gravitational-wave background (GWB) from binary SMBHs in merging galaxies throughout the Universe. Using observations from the Parkes Pulsar Timing Array, we constrain the fractional GWB energy density with 95% confidence to be ${\Omega}_{GW}(H_0/73 {\rm km} {\rm s}^{-1} {\rm Mpc}^{-1})^2 < 1.3 \times 10^{-9}$ at a frequency of 2.8 nHz, which is approximately a factor of six more stringent than previous limits. We compare our limit to models of the SMBH population and find inconsistencies at confidence levels between 46% and 91%. For example, the standard galaxy formation model implemented in the Millennium simulations is inconsistent with our limit with 50% probability.

The correlation between 21cm fluctuations and galaxies is sensitive to the astrophysical properties of the galaxies that drove reionization. Thus, detailed measurements of the cross-power spectrum and its evolution could provide a powerful measurement both of the properties of early galaxies and the process of reionization. In this paper, we study the evolution of the cross-power spectrum between 21cm emission and galaxies using a model which combines the hierarchical galaxy formation model GALFORM implemented within the Millennium-II dark matter simulation, with a semi-numerical scheme to describe the resulting ionization structure. We find that inclusion of different feedback processes changes the cross-power spectrum shape and amplitude. In particular, the feature in the cross-power spectrum corresponding to the size of ionized regions is significantly affected by supernovae feedback. We calculate predicted observational uncertainties of the cross-correlation coefficient based on specifications of the Murchison Widefield Array (MWA) combined with galaxy surveys of varying area and depth. We find that the cross-power spectrum could be detected over several square degrees of galaxy survey with galaxy redshift errors less than 0.1.

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.

We present new wide-field observations of Centaurus A (Cen A) and the surrounding region at 118 MHz with the Murchison Widefield Array (MWA) 32-tile prototype, with which we investigate the spectral-index distribution of Cen A's giant radio lobes. We compare our images to 1.4 GHz maps of Cen A and compute spectral indices using temperature-temperature plots and spectral tomography. We find that the morphologies at 118 MHz and 1.4 GHz match very closely apart from an extra peak in the southern lobe at 118 MHz, which provides tentative evidence for the existence of a southern counterpart to the northern middle lobe of Cen A. Our spatially-averaged spectral indices for both the northern and southern lobes are consistent with previous analyses, however we find significant spatial variation of the spectra across the extent of each lobe. Both the spectral-index distribution and the morphology at low radio frequencies support a scenario of multiple outbursts of activity from the central engine. Our results are consistent with inverse-Compton modelling of radio and gamma-ray data that supports a value for the lobe age of between 10 and 80 Myr.

The Murchison Widefield Array (MWA) is a new low frequency interferomeric radio telescope. The MWA is the low frequency precursor to the Square Kilometre Array (SKA) and is the first of three SKA precursors to be operational, supporting a varied science mission ranging from the attempted detection of the Epoch of Reionisation to the monitoring of solar flares and space weather. We explore the possibility that the MWA can be used for the purposes of Space Situational Awareness (SSA). In particular we propose that the MWA can be used as an element of a passive radar facility operating in the frequency range 87.5 - 108 MHz (the commercial FM broadcast band). In this scenario the MWA can be considered the receiving element in a bi-static radar configuration, with FM broadcast stations serving as non-cooperative transmitters. The FM broadcasts propagate into space, are reflected off debris in Earth orbit, and are received at the MWA. The imaging capabilities of the MWA can be used to simultaneously detect multiple pieces of space debris, image their positions on the sky as a function of time, and provide tracking data that can be used to determine orbital parameters. Such a capability would be a valuable addition to Australian and global SSA assets, in terms of southern and eastern hemispheric coverage. We provide a feasibility assessment of this proposal, based on simple calculations and electromagnetic simulations that shows the detection of sub-metre size debris should be possible (debris radius of >0.5 m to ~1000 km altitude). We also present a proof-of-concept set of observations that demonstrate the feasibility of the proposal, based on the detection and tracking of the International Space Station via reflected FM broadcast signals originating in south-west Western Australia. These observations broadly validate our calculations and simulations.

We present a Stokes I, Q and U survey at 189 MHz with the Murchison Widefield Array 32-element prototype covering 2400 square degrees. The survey has a 15.6 arcmin angular resolution and achieves a noise level of 15 mJy/beam. We demonstrate a novel interferometric data analysis that involves calibration of drift scan data, integration through the co-addition of warped snapshot images and deconvolution of the point spread function through forward modeling. We present a point source catalogue down to a flux limit of 4 Jy. We detect polarization from only one of the sources, PMN J0351-2744, at a level of 1.8 \pm 0.4%, whereas the remaining sources have a polarization fraction below 2%. Compared to a reported average value of 7% at 1.4 GHz, the polarization fraction of compact sources significantly decreases at low frequencies. We find a wealth of diffuse polarized emission across a large area of the survey with a maximum peak of ~13 K, primarily with positive rotation measure values smaller than +10 rad/m^2. The small values observed indicate that the emission is likely to have a local origin (closer than a few hundred parsecs). There is a large sky area at 2^h30^m where the diffuse polarized emission rms is fainter than 1 K. Within this area of low Galactic polarization we characterize the foreground properties in a cold sky patch at $(\alpha,\delta) = (4^h,-27^\circ.6)$ in terms of three dimensional power spectra

The observed power spectrum of redshifted 21cm fluctuations is known to be sensitive to the astrophysical properties of the galaxies that drove reionization. Thus, detailed measurements of the 21cm power spectrum and its evolution could lead to measurements of the properties of early galaxies that are otherwise inaccessible. In this paper, we study the effect of mass and redshift dependent escape fractions of ionizing radiation on the ability of forthcoming experiments to constrain galaxy formation via the redshifted 21cm power spectrum. We use a model for reionization which combines the hierarchical galaxy formation model GALFORM implemented within the Millennium-II dark matter simulation, with a semi-numerical scheme to describe the resulting ionization structure. Using this model we show that the structure and distribution of ionised regions at fixed neutral fraction, and hence the slope and amplitude of the 21 cm power spectrum, is dependent on the variation of ionising photon escape fraction with galaxy mass and redshift. However, we find that the influence of the unknown escape fraction and its evolution is smaller than the dominant astrophysical effect provided by SNe feedback strength in high redshift galaxies. The unknown escape fraction of ionizing radiation from galaxies is therefore unlikely to prevent measurement of the properties of high redshift star formation using observations of the 21cm power spectrum.

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.

The Murchison Widefield Array is a low frequency (80 - 300 MHz) SKA Precursor, comprising 128 aperture array elements (known as tiles) distributed over an area of 3 km diameter. The MWA is located at the extraordinarily radio quiet Murchison Radioastronomy Observatory in the mid-west of Western Australia, the selected home for the Phase 1 and Phase 2 SKA low frequency arrays. The MWA science goals include: 1) detection of fluctuations in the brightness temperature of the diffuse redshifted 21 cm line of neutral hydrogen from the epoch of reionisation; 2) studies of Galactic and extragalactic processes based on deep, confusion-limited surveys of the full sky visible to the array; 3) time domain astrophysics through exploration of the variable radio sky; and 4) solar imaging and characterisation of the heliosphere and ionosphere via propagation effects on background radio source emission. This paper concentrates on the capabilities of the MWA for solar science and summarises some of the solar science results to date, in advance of the initial operation of the final instrument in 2013.

The number of strong (equivalent width > 1A) MgII absorbers observed towards Gamma-ray bursts (GRBs) has been found to be statistically larger than the number of strong absorbers towards quasi-stellar objects (QSOs). We formalize this "MgII problem" and present a detailed explanation of the statistical tools required to assess the significance of the discrepancy. We find the problem exists at the 4{\sigma} level for GRBs with high-resolution spectra. It has been suggested that the discrepancy can be resolved by the combination of a dust obscuration bias towards QSOs, and a strong gravitational lensing bias towards GRBs. We investigate one of the two most probable lensed GRBs that we presented in our previous work (GRB020405; Rapoport et al.) and find it not to be strongly gravitationally lensed, constraining the percentage of lensed GRBs to be < 35% (2{\sigma}). Dust obscuration of QSOs has been estimated to be a significant effect with dusty MgII systems removing ~20% of absorbed objects from flux-limited QSO samples. We find that if ~30% of the strong MgII systems towards QSOs are missing from the observed samples, then GRBs and QSOs would have comparable numbers of absorbers per unit redshift. Thus, gravitational lensing bias is likely to make only a modest contribution to solving the MgII problem. However, if the dust obscuration bias has been slightly underestimated, the MgII problem would no longer persist.

The Murchison Widefield Array is a low frequency (80 - 300 MHz) SKA Precursor, comprising 128 aperture array elements distributed over an area of 3 km diameter. The MWA is located at the extraordinarily radio quiet Murchison Radioastronomy Observatory in the mid-west of Western Australia, the selected home for the Phase 1 and Phase 2 SKA low frequency arrays. The MWA science goals include: 1) detection of fluctuations in the brightness temperature of the diffuse redshifted 21 cm line of neutral hydrogen from the epoch of reionisation; 2) studies of Galactic and extragalactic processes based on deep, confusion-limited surveys of the full sky visible to the array; 3) time domain astrophysics through exploration of the variable radio sky; and 4) solar imaging and characterisation of the heliosphere and ionosphere via propagation effects on background radio source emission. This paper will focus on a brief discussion of the as-built MWA system, highlighting several novel characteristics of the instrument, and a brief progress report (as of June 2012) on the final construction phase. Practical completion of the MWA is expected in November 2012, with commissioning commencing from approximately August 2012 and operations commencing near mid 2013. A brief description of recent science results from the MWA prototype instrument is given.

Large surveys for Lyman-alpha emitting (LAE) galaxies have been proposed as a new method for measuring clustering of the galaxy population at high redshift with the goal of determining cosmological parameters. However, Lyman-alpha radiative transfer effects may modify the observed clustering of LAE galaxies in a way that mimics gravitational effects, potentially reducing the precision of cosmological constraints. For example, the effect of the linear redshift-space distortion on the power spectrum of LAE galaxies is potentially degenerate with Lyman-alpha radiative transfer effects owing to the dependence of observed flux on intergalactic medium velocity gradients. In this paper, we show that the three-point function (bispectrum) can distinguish between gravitational and non-gravitational effects, and thus breaks these degeneracies, making it possible to recover cosmological parameters from LAE galaxy surveys. Constraints on the angular diameter distance and the Hubble expansion rate can also be improved by combining power spectrum and bispectrum measurements.

A new generation of low frequency radio telescopes is seeking to observe the redshifted 21 cm signal from the Epoch of Reionization (EoR), requiring innovative methods of calibration and imaging to overcome the difficulties of widefield low frequency radio interferometry. Precise calibration will be required to separate the small expected EoR signal from the strong foreground emission at the frequencies of interest between 80 and 300 MHz. The Moon may be useful as a calibration source for detection of the EoR signature, as it should have a smooth and predictable thermal spectrum across the frequency band of interest. Initial observations of the Moon with the Murchison Widefield Array 32 tile prototype show that the Moon does exhibit a similar trend to that expected for a cool thermally emitting body in the observed frequency range, but that the spectrum is corrupted by reflected radio emission from Earth. In particular, there is an abrupt increase in the observed flux density of the Moon within the internationally recognised Frequency Modulated (FM) radio band. The observations have implications for future low frequency surveys and EoR detection experiments that will need to take this reflected emission from the Moon into account. The results also allow us to estimate the equivalent isotropic power emitted by the Earth in the FM band and to determine how bright the Earth might appear at metre wavelengths to an observer beyond our own solar system.

We investigate the impact of feedback - from supernovae (SNe), active galactic nuclei (AGN) and a photo-ionizing background at high redshifts - on the neutral atomic hydrogen (HI) mass function, the $b_{\rm J}$ band luminosity function, and the spatial clustering of these galaxies at $z$=0. We use a version of the semi-analytical galaxy formation model GALFORM that calculates self-consistently the amount of HI in a galaxy as a function of cosmic time and links its star formation rate to its mass of molecular hydrogen (H$_2$). We find that a systematic increase or decrease in the strength of SNe feedback leads to a systematic decrease or increase in the amplitudes of the luminosity and HI mass functions, but has little influence on their overall shapes. Varying the strength of AGN feedback influences only the numbers of the brightest or most HI massive galaxies, while the impact of varying the strength of photo-ionization feedback is restricted to changing the numbers of the faintest or least HI massive galaxies.Our results suggest that the HI mass function is a more sensitive probe of the consequences of cosmological reionization for galaxy formation than the luminosity function. We find that increasing the strength of any of the modes of feedback acts to weaken the clustering strength of galaxies, regardless of their HI-richness. In contrast, weaker AGN feedback has little effect on the clustering strength whereas weaker SNe feedback increases the clustering strength of HI-poor galaxies more strongly than HI-rich galaxies. These results indicate that forthcoming HI surveys on next generation radio telescopes such as the Square Kilometre Array and its pathfinders will be exploited most fruitfully as part of multiwavelength survey campaigns.

Understanding the epoch of reionization and the properties of the first galaxies represents an important goal for modern cosmology. The structure of reionization, and hence the observed power spectrum of redshifted 21cm fluctuations are known to be sensitive to the astrophysical properties of the galaxies that drove reionization. Thus, detailed measurements of the 21cm power spectrum and its evolution could lead to measurements of the properties of early galaxies that are otherwise inaccessible. In this paper, we make predictions for the ionised structure during reionization and the 21cm power spectrum based on detailed models of galaxy formation. We combine the semi-analytic GALFORM model implemented within the Millennium-II dark matter simulation, with a semi-numerical scheme to describe the resulting ionization structure. Using these models we show that the details of SNe and radiative feedback affect the structure and distribution of ionised regions, and hence the slope and amplitude of the 21 cm power spectrum. These results indicate that forthcoming measurements of the 21cm power-spectrum could be used to uncover details of early galaxy formation. We find that the strength of SNe feedback is the dominant effect governing the evolution of structure during reionization. In particular we show SNe feedback to be more important than radiative feedback, the presence of which we find does not influence either the total stellar mass or overall ionising photon budget. Thus, if SNe feedback is effective at suppressing star formation in high redshift galaxies, we find that photoionization feedback does not lead to self-regulation of the reionization process as has been thought.

We investigate the effects of gravitational waves (GWs) from a simulated population of binary super-massive black holes (SMBHs) on pulsar timing array datasets. We construct a distribution describing the binary SMBH population from an existing semi-analytic galaxy formation model. Using realizations of the binary SMBH population generated from this distribution, we simulate pulsar timing datasets with GW-induced variations. We find that the statistics of these variations do not correspond to an isotropic, stochastic GW background. The "Hellings & Downs" correlations between simulated datasets for different pulsars are recovered on average, though the scatter of the correlation estimates is greater than expected for an isotropic, stochastic GW background. These results are attributable to the fact that just a few GW sources dominate the GW-induced variations in every Fourier frequency bin of a 5-year dataset. Current constraints on the amplitude of the GW signal from binary SMBHs will be biased. Individual binary systems are likely to be detectable in 5-year pulsar timing array datasets where the noise is dominated by GW-induced variations. Searches for GWs in pulsar timing array data therefore need to account for the effects of individual sources of GWs.

Using the final 128 antenna locations of the Murchison Widefield Array (MWA), we calculate its sensitivity to the Epoch of Reionization (EoR) power spectrum of red- shifted 21 cm emission for a fiducial model and provide the tools to calculate the sensitivity for any model. Our calculation takes into account synthesis rotation, chro- matic and asymmetrical baseline effects, and excludes modes that will be contaminated by foreground subtraction. For the fiducial model, the MWA will be capable of a 14{\sigma} detection of the EoR signal with one full season of observation on two fields (900 and 700 hours).

We present an analysis of the Hubble diagram for 12 Type Ia supernovae (SNe Ia) observed in the near-infrared J and H bands. We select SNe exclusively from the redshift range 0.03 < z < 0.09 to reduce uncertainties coming from peculiar velocities while remaining in a cosmologically well-understood region. All of the SNe in our sample exhibit no spectral or B-band light-curve peculiarities and lie in the B-band stretch range of 0.8-1.15. Our results suggest that SNe Ia observed in the near-infrared (NIR) are the best known standard candles. We fit previously determined NIR light-curve templates to new high-precision data to derive peak magnitudes and to determine the scatter about the Hubble line. Photometry of the 12 SNe is presented in the natural system. Using a standard cosmology of (H_0, Omega_m, Lambda) = (70,0.27,0.73) we find a median J-band absolute magnitude of M_J = -18.39 with a scatter of 0.116 and a median H-band absolute magnitude of M_H = -18.36 with a scatter of 0.085. The scatter in the H band is the smallest yet measured. We search for correlations between residuals in the J- and H-band Hubble diagrams and SN properties, such as SN colour, B-band stretch and the projected distance from host-galaxy centre. The only significant correlation is between the J-band Hubble residual and the J-H pseudo-colour. We also examine how the scatter changes when fewer points in the near-infrared are used to constrain the light curve. With a single point in the H band taken anywhere from 10 days before to 15 days after B-band maximum light and a prior on the date of H-band maximum set from the date of B-band maximum, we find that we can measure distances to an accuracy of 6%. The precision of SNe Ia in the NIR provides new opportunities for precision measurements of both the expansion history of the universe and peculiar velocities of nearby galaxies.

The Murchison Widefield Array (MWA) is one of three Square Kilometre Array Precursor telescopes and is located at the Murchison Radio-astronomy Observatory in the Murchison Shire of the mid-west of Western Australia, a location chosen for its extremely low levels of radio frequency interference. The MWA operates at low radio frequencies, 80-300 MHz, with a processed bandwidth of 30.72 MHz for both linear polarisations, and consists of 128 aperture arrays (known as tiles) distributed over a ~3 km diameter area. Novel hybrid hardware/software correlation and a real-time imaging and calibration systems comprise the MWA signal processing backend. In this paper the as-built MWA is described both at a system and sub-system level, the expected performance of the array is presented, and the science goals of the instrument are summarised.

Sixty percent of gamma-ray bursts (GRBs) reveal strong MgII absorbing systems, which is a factor of ~2 times the rate seen along lines-of-sight to quasars. The discrepancy in the covering factor is most likely to be the result of either quasars being obscured due to dust, or the consequence of many GRBs being strongly gravitationally lensed. We analyze observations of GRBs that show strong foreground MgII absorption. We Monte-Carlo the distances between foreground galaxies in the HUDF and lines of sight distributed randomly within a radius derived from the covering factor, and find that galaxies are located systematically closer to the position of the observed GRBs than expected for random lines of sight. This over-density at small impact parameters is statistically more robust than the well known excess of MgII absorbers among GRB afterglow spectra, and presents a new puzzle for MgII absorption studies. The over-density cannot be explained by obscuration in the GRB sample, but is a natural consequence of gravitational lensing. Upon examining the particular configurations of galaxies near a sample of GRBs with strong MgII absorption, we find several intriguing lensing candidates. Our results suggest that lensing provides a viable contribution to the observed enhancement of strong MgII absorption along lines of sight to GRBs, and we outline the future observations required to test this hypothesis conclusively.

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.

Antenna layout is an important design consideration for radio interferometers because it determines the quality of the snapshot point spread function (PSF, or array beam). This is particularly true for experiments targeting the 21 cm Epoch of Reionization signal as the quality of the foreground subtraction depends directly on the spatial dynamic range and thus the smoothness of the baseline distribution. Nearly all sites have constraints on where antennas can be placed---even at the remote Australian location of the MWA (Murchison Widefield Array) there are rock outcrops, flood zones, heritages areas, emergency runways and trees. These exclusion areas can introduce spatial structure into the baseline distribution that enhance the PSF sidelobes and reduce the angular dynamic range. In this paper we present a new method of constrained antenna placement that reduces the spatial structure in the baseline distribution. This method not only outperforms random placement algorithms that avoid exclusion zones, but surprisingly outperforms random placement algorithms without constraints to provide what we believe are the smoothest constrained baseline distributions developed to date. We use our new algorithm to determine antenna placements for the originally planned MWA, and present the antenna locations, baseline distribution, and snapshot PSF for this array choice.

Measurements of the intergalactic medium (IGM) temperature provide a potentially powerful constraint on the reionisation history due to the thermal imprint left by the photo-ionisation of neutral hydrogen. However, until recently IGM temperature measurements were limited to redshifts 2 < z < 4.8, restricting the ability of these data to probe the reionisation history at z > 6. In this work, we use recent measurements of the IGM temperature in the near-zones of seven quasars at z ~ 5.8 - 6.4, combined with a semi-numerical model for inhomogeneous reionisation, to establish new constraints on the redshift at which hydrogen reionisation completed. We calibrate the model to reproduce observational constraints on the electron scattering optical depth and the HI photo-ionisation rate, and compute the resulting spatially inhomogeneous temperature distribution at z ~ 6 for a variety of reionisation scenarios. Under standard assumptions for the ionising spectra of population-II sources, the near-zone temperature measurements constrain the redshift by which hydrogen reionisation was complete to be z > 7.9 (6.5) at 68 (95) per cent confidence. We conclude that future temperature measurements around other high redshift quasars will significantly increase the power of this technique, enabling these results to be tightened and generalised.

The 21cm forest -- HI absorption features in the spectra of high-redshift radio sources -- can potentially provide a unique probe of the largely neutral intergalactic medium (IGM) during the epoch of reionization. We present simulations of the 21cm forest due to the large scale structure of the reionization-era IGM, including a prescription for x-ray heating and the percolation of photoionization bubbles. We show that, if detected with future instruments such as the Square Kilometer Array (SKA), the 21cm forest can provide a significant constraint on the thermal history of the IGM. Detection will be aided by consideration of the sudden increase in signal variance at the onset of 21cm absorption. If radio foregrounds and the intrinsic source spectra are well understood, the flux decrement over wide bandwidths can also improve detection prospects. Our analysis accounts for the possibility of narrow absorption lines from intervening dense regions, but, unlike previous studies, our results do not depend on their properties. Assuming x-ray heating corresponding to a local stellar population, we estimate that a statistically significant detection of 21cm absorption could be made by SKA in less than a year of observing against a Cygnus A-type source at $z \sim 9$, as opposed to nearly a decade for a significant detection of the detailed forest features. We discuss observational challenges due to uncertainties regarding the abundance of background sources and the strength of the 21cm absorption signal.

We present measurements of the intergalactic medium (IGM) temperature within ~5 proper Mpc of seven luminous quasars at z~6. The constraints are obtained from the Doppler widths of Lyman-alpha absorption lines in the quasar near-zones and build upon our previous measurement for the z=6.02 quasar SDSS J0818+1722. The expanded data set, combined with an improved treatment of systematic uncertainties, yields an average temperature at the mean density of log (T_0/K) = 4.21 \pm 0.03 (\pm^0.06_0.07) at 68 (95) per cent confidence for a flat prior distribution over 3.2 < log (T_0/K) < 4.8. In comparison, temperatures measured from the general IGM at z~5 are ~0.3 dex cooler, implying an additional source of heating around these quasars which is not yet present in the general IGM at slightly lower redshift. This heating is most likely due to the recent reionisation of He-II in vicinity of these quasars, which have hard, non-thermal ionising spectra. The elevated temperatures may therefore represent evidence for the earliest stages of He-II reionisation in the most biased regions of the high-redshift Universe. The temperature as a function of distance from the quasars is consistent with being constant, log(T_0/K)~4.2, with no evidence for a line-of-sight thermal proximity effect. However, the limited extent of the quasar near-zones prevents the detection of He-III regions larger than ~5 proper Mpc. Under the assumption the quasars have reionised the He-II in their vicinity, we infer the data are consistent with an average optically bright phase of duration in excess of 10^6.5 yr. These measurements represent the highest-redshift IGM temperature constraints to date, and thus provide a valuable data set for confronting models of H-I reionisation.

High redshift measurements of the baryonic acoustic oscillation scale (BAO) from large Ly-alpha forest surveys represent the next frontier of dark energy studies. As part of this effort, efficient simulations of the BAO signature from the Ly-alpha forest will be required. We construct a model for producing fast, large volume simulations of the Ly-alpha forest for this purpose. Utilising a calibrated semi-analytic approach, we are able to run very large simulations in 1 Gpc^3 volumes which fully resolve the Jeans scale in less than a day on a desktop PC using a GPU enabled version of our code. The Ly-alpha forest spectra extracted from our semi-analytical simulations are in excellent agreement with those obtained from a fully hydrodynamical reference simulation. Furthermore, we find our simulated data are in broad agreement with observational measurements of the flux probability distribution and 1D flux power spectrum. We are able to correctly recover the input BAO scale from the 3D Ly-alpha flux power spectrum measured from our simulated data, and estimate that a BOSS-like 10^4 deg^2 survey with ~15 background sources per square degree and a signal-to-noise of ~5 per pixel should achieve a measurement of the BAO scale to within ~1.4 per cent. We also use our simulations to provide simple power-law expressions for estimating the fractional error on the BAO scale on varying the signal-to-noise and the number density of background sources. The speed and flexibility of our approach is well suited for exploring parameter space and the impact of observational and astrophysical systematics on the recovery of the BAO signature from forthcoming large scale spectroscopic surveys.

Measurement of redshifted 21-cm emission from neutral hydrogen promises to be the most effective method for studying the reionisation history of hydrogen and, indirectly, the first galaxies. These studies will be limited not by raw sensitivity to the signal, but rather, by bright foreground radiation from Galactic and extragalactic radio sources and the Galactic continuum. In addition, leakage due to gain errors and non-ideal feeds conspire to further contaminate low-frequency radio obsevations. This leakage leads to a portion of the complex linear polarisation signal finding its way into Stokes I, and inhibits the detection of the non-polarised cosmological signal from the epoch of reionisation. In this work, we show that rotation measure synthesis can be used to recover the signature of cosmic hydrogen reionisation in the presence of contamination by polarised foregrounds. To achieve this, we apply the rotation measure synthesis technique to the Stokes I component of a synthetic data cube containing Galactic foreground emission, the effect of instrumental polarisation leakage, and redshifted 21-cm emission by neutral hydrogen from the epoch of reionisation. This produces an effective Stokes I Faraday dispersion function for each line of sight, from which instrumental polarisation leakage can be fitted and subtracted. Our results show that it is possible to recover the signature of reionisation in its late stages (z ~ 7) by way of the 21-cm power spectrum, as well as through tomographic imaging of ionised cavities in the intergalactic medium.

Narrowband observations have detected no Ly-alpha emission within a 70 pMpc^3 volume centered on the z = 2.168 quasar PKS 0424-131. This is in contrast to surveys of Ly-alpha emitters in the field at similar redshifts and flux limits, which indicate that tens of sources should be visible within the same volume. The observed difference indicates that the quasar environment has a significant influence on the observed density of Ly-alpha emitters. To quantify this effect we have constructed a semi-analytic model to simulate the effect of a luminous quasar on nearby Ly-alpha emitters. We find the null detection around PKS 0424-131 implies that the minimum isothermal temperature of Ly-alpha emitter host halos is greater than 3.4 x 10^6 K (68% level), corresponding to a virial mass of ~1.2 x 10^12 M_solar. This indicates that the intense UV emission of the quasar may be suppressing the star formation in nearby galaxies. Our study illustrates that low redshift quasar environments may serve as a surrogate for studying the radiative suppression of galaxy formation during the epoch of reionization.

The measurement of the light scattered from extrasolar planets informs atmospheric and formation models. With the discovery of many hot Jupiter planets orbiting nearby stars, this motivates the development of robust methods of characterisation from follow up observations. In this paper we discuss two methods for determining the planetary albedo in transiting systems. First, the most widely used method for measuring the light scattered by hot Jupiters (Collier Cameron et al.) is investigated for application for typical echelle spectra of a transiting planet system, showing that detection requires high signal-to-noise ratio data of bright planets. Secondly a new Fourier analysis method is also presented, which is model-independent and utilises the benefits of the reduced number of unknown parameters in transiting systems. This approach involves solving for the planet and stellar spectra in Fourier space by least-squares. The sensitivities of the methods are determined via Monte Carlo simulations for a range of planet-to-star fluxes. We find the Fourier analysis method to be better suited to the ideal case of typical observations of a well constrained transiting system than the Collier Cameron et al. method. We apply the Fourier analysis method for extracting the light scattered by transiting hot Jupiters from high resolution spectra to echelle spectra of HD 209458 and HD 189733. Unfortunately we are unable to improve on the previous upper limit of the planet-to-star flux for HD 209458b set by space-based observations. A 1{\sigma}upper limit on the planet-to-star flux of HD 189733b is measured in the wavelength range of 558.83-599.56 nm yielding {\epsilon} < 4.5 \times 10-4. Improvement in the measurement of the upper limit of the planet-to-star flux of this system, with ground-based capabilities, requires data with a higher signal-to-noise ratio, and increased stability of the telescope.

Direct measurements of dark matter distributions in galaxies are currently only possible through the use of gravitational lensing observations. Combinations of lens modelling and stellar velocity dispersion measurements provide the best constraints on dark matter distributions in individual galaxies, however they can be quite complex. In this paper, we use observations and simulations of gravitational microlensing to measure the smooth (dark) matter mass fraction at the position of lensed images in three lens galaxies: MG 0414+0534, SDSS J0924+0219 and Q2237+0305. The first two systems consist of early-type lens galaxies, and both display a flux ratio anomaly in their close image pair. Anomalies such as these suggest a high smooth matter percentage is likely, and indeed we prefer ~50 per cent smooth matter in MG 0414+0534, and ~80 per cent in SDSS J0924+0219 at the projected locations of the lensed images. Q2237+0305 differs somewhat in that its lensed images lie in the central kiloparsec of the barred spiral lens galaxy, where we expect stars to dominate the mass distribution. In this system, we find a smooth matter percentage that is consistent with zero.

The observed number counts of high-redshift galaxy candidates have been used to build up a statistical description of star-forming activity at redshift z >~ 7, when galaxies reionized the Universe. Standard models predict that a high incidence of gravitational lensing will probably distort measurements of flux and number of these earliest galaxies. The raw probability of this happening has been estimated to be ~ 0.5 percent, but can be larger owing to observational biases. Here we report that gravitational lensing is likely to dominate the observed properties of galaxies with redshifts of z >~ 12, when the instrumental limiting magnitude is expected to be brighter than the characteristic magnitude of the galaxy sample. The number counts could be modified by an order of magnitude, with most galaxies being part of multiply imaged systems, located less than 1 arcsec from brighter foreground galaxies at z ~ 2. This lens-induced association of high-redshift and foreground galaxies has perhaps already been observed among a sample of galaxy candidates identified at z ~ 10.6. Future surveys will need to be designed to account for a significant gravitational lensing bias in high-redshift galaxy samples.

We present a method for subtracting point sources from interferometric radio images via forward modeling of the instrument response and involving an algebraic nonlinear minimization. The method is applied to simulated maps of the Murchison Wide-field Array but is generally useful in cases where only image data are available. After source subtraction, the residual maps have no statistical difference to the expected thermal noise distribution at all angular scales, indicating high effectiveness in the subtraction. Simulations indicate that the errors in recovering the source parameters decrease with increasing signal-to-noise ratio, which is consistent with the theoretical measurement errors. In applying the technique to simulated snapshot observations with the Murchison Wide-field Array, we found that all 101 sources present in the simulation were recovered with an average position error of 10 arcsec and an average flux density error of 0.15%. This led to a dynamic range increase of approximately 3 orders of magnitude. Since all the sources were deconvolved jointly, the subtraction was not limited by source sidelobes but by thermal noise. This technique is a promising deconvolution method for upcoming radio arrays with a huge number of elements, and a candidate for the difficult task of subtracting foreground sources from observations of the 21 cm neutral Hydrogen signal from the epoch of reionization.

Bright point sources associated with extragalactic AGN and radio galaxies are an important foreground for low frequency radio experiments aimed at detecting the redshifted 21cm emission from neutral hydrogen during the epoch of reionization. The frequency dependence of the synthesized beam implies that the sidelobes of these sources will move across the field of view as a function of observing frequency, hence frustrating line-of-sight foreground subtraction techniques. We describe a method for subtracting these point sources from dirty maps produced by an instrument such as the MWA. This technique combines matched filters with an iterative centroiding scheme to locate and characterize point sources in the presence of a diffuse background. Simulations show that this technique can improve the dynamic range of EOR maps by 2-3 orders of magnitude.

The thermal state of the intergalactic medium (IGM) provides an indirect probe of both the HI and HeII reionisation epochs. Current constraints on the IGM temperature from the Lya forest are restricted to the redshift range 2<z<4.5, limiting the ability to probe the thermal memory of HI reionisation toward higher redshift. In this work, we present the first direct measurement of the IGM temperature around a z=6 quasar by analysing the Doppler widths of Lya absorption lines in the proximity zone of SDSS J0818+1722. We use a high resolution (R= 40000) Keck/HIRES spectrum in combination with detailed numerical modelling to obtain the temperature at mean density, T_0=23600\pm^5000_6900K (\pm^9200_9300K) at 68 (95) per cent confidence assuming a prior probability 13500K<T_0<38500 K following HI and HeII reionisation. This enables us to place an upper limit on the redshift of HI reionisation, z_H, within 33 comoving Mpc of SDSS J0818+1722. If the quasar reionises the HeII in its vicinity, then in the limit of instantaneous reionisation we infer z_H<9.0 (11.0) at 68 (95) per cent confidence assuming photoheating is the dominant heat source and that HI reionisation is driven by ionising sources with soft spectra, typical of population II stars. If the HI and HeII in the IGM around SDSS J0818+1722 are instead reionised simultaneously by a population of massive metal-free stars, characterised by very hard ionising spectra, we obtain a tighter upper limit of z_H<8.4 (9.4). Initiating reionisation at higher redshifts produces temperatures which are too low with respect to our constraint unless the HI ionising sources or the quasar itself have spectra significantly harder than typically assumed.

Measurement of the spatial distribution of neutral hydrogen via the redshifted 21 cm line promises to revolutionize our knowledge of the epoch of reionization and the first galaxies, and may provide a powerful new tool for observational cosmology from redshifts 1<z<4 . In this review we discuss recent advances in our theoretical understanding of the epoch of reionization (EoR), the application of 21 cm tomography to cosmology and measurements of the dark energy equation of state after reionization, and the instrumentation and observational techniques shared by 21 cm EoR and post reionization cosmology machines. We place particular emphasis on the expected signal and observational capabilities of first generation 21 cm fluctuation instruments.

High-redshift gamma-ray bursts (GRBs) offer an extraordinary opportunity to study aspects of the early Universe, including the cosmic star formation rate (SFR). Motivated by the two recent highest-z GRBs, GRB 080913 at z = 6.7 and GRB 090423 at z = 8.1, and more than four years of Swift observations, we first confirm that the GRB rate does not trace the SFR in an unbiased way. Correcting for this, we find that the implied SFR to beyond z = 8 is consistent with LBG-based measurements after accounting for unseen galaxies at the faint end of the UV luminosity function. We show that this provides support for the integrated star formation in the range 6 < z < 8 to have been alone sufficient to reionize the Universe.

In this paper, we explore the constraints on the reionization history that are provided by current observations of the Lyman-alpha forest and the CMB. Rather than using a particular semi-analytic model, we take the novel approach of parametrizing the ionizing sources with arbitrary functions, and perform likelihood analyses to constrain possible reionization histories. We find model independent conclusions that reionization is likely to be mostly complete by z=8 and that the IGM was 50% ionized at z=9-10. Upcoming low-frequency observations of the redshifted 21 cm line of neutral hydrogen are expected to place significantly better constraints on the hydrogen neutral fraction at 6<z<12. We use our constraints on the reionization history to predict the likely amplitude of the 21 cm power spectrum and show that observations with the highest signal-to-noise ratio will most likely be made at frequencies corresponding to z=9-10. This result provides an important guide to the upcoming 21 cm observations. Finally, we assess the impact that measurement of the neutral fraction will have on our knowledge of reionization and the early source population. Our results show that a single measurement of the neutral fraction mid-way through the reionization era will significantly enhance our knowledge of the entire reionization history.

We assess the effect of a population of high-redshift quasars on the 21-cm power spectrum during the epoch of reionisation. Our approach is to implement a semi-numerical scheme to calculate the three-dimensional structure of ionised regions surrounding massive halos at high redshift. We include the ionising influence of luminous quasars by populating a simulated overdensity field with quasars using a Monte Carlo Markov Chain algorithm. We find that quasars modify both the amplitude and shape of the power spectrum at a level which is of the same order as the fractional contribution to reionisation. The modification is found both at constant redshift and at constant global neutral fraction, and arises because ionising photons produced by quasars are biased relative to the density field at a level that is higher than steller ionising photons. Our results imply that quasar ionisation will need to be included in detailed modelling of observed 21-cm power spectra.

The identification of an extrasolar planet as Earth-like will depend on the detection of atmospheric signatures or surface non-uniformities. In this paper we present spatially unresolved flux light curves of Earth for the purpose of studying a prototype extrasolar terrestrial planet. Our monitoring of the photometric variability of earthshine revealed changes of up to 23 % per hour in the brightness of Earth's scattered light at around 600 nm, due to the removal of specular reflection from the view of the Moon. This variability is accompanied by reddening of the spectrum, and results from a change in surface properties across the continental boundary between the Indian Ocean and Africa's east coast. Our results based on earthshine monitoring indicate that specular reflection should provide a useful tool in determining the presence of liquid water on extrasolar planets via photometric observations.

Observations of distances to Type-Ia supernovae can be explained by cosmological models that include either a gigaparsec-scale void, or a cosmic flow, without the need for Dark Energy. Instead of invoking dark energy, these inhomogeneous models instead violate the Copernican Principle. we show that current cosmological observations (Supernovae, Baryon Acoustic Oscillations and estimates of the Hubble parameters based on the age of the oldest stars) are not able to rule out inhomogeneous anti-Copernican models. The next generation of surveys for baryonic acoustic oscillations will be sufficiently precise to either validate the Copernican Principle or determine the existence of a local Gpc scale inhomogeneity.

It has recently been suggested that the power spectrum of redshifted 21cm fluctuations could be used to measure the scale of baryonic acoustic oscillations (BAOs) during the reionisation era. The resulting measurements are potentially as precise as those offered by the next generation of galaxy redshift surveys at lower redshift. However unlike galaxy redshift surveys, which in the linear regime are subject to a scale independent galaxy bias, the growth of ionised regions during reionisation is thought to introduce a strongly scale dependent relationship between the 21cm and mass power spectra. We use a semi-numerical model for reionisation to assess the impact of ionised regions on the precision and accuracy with which the BAO scale could be measured using redshifted 21cm observations. For a model in which reionisation is completed at z~6, we find that the constraints on the BAO scale are not systematically biased at z > 6.5. In this scenario, and assuming the sensitivity attainable with a low-frequency array comprising 10 times the collecting area of the Murchison Widefield Array, the BAO scale could be measured to within 1.5 per cent in the range 6.5 < z < 7.5.

Observations of gravitational microlensing in multiply imaged quasars currently provide the only direct probe of quasar emission region structure on sub-microarcsecond scales. Analyses of microlensing variability are observationally expensive, requiring long-term monitoring of lensed systems. Here we demonstrate a technique for constraining the size of the quasar continuum emission region as a function of wavelength using single-epoch multi-band imaging. We have obtained images of the lensed quasar MG 0414+0534 in five wavelength bands using the Magellan 6.5-metre Baade telescope at Las Campanas Observatory, Chile. These data, in combination with two existing epochs of Hubble Space Telescope data, are used to model the size of the continuum emission region $\sigma$ as a power-law in wavelength, $\sigma\propto\lambda^\nu$. We place an upper limit on the Gaussian width of the $r^\prime$-band emission region of $1.80 \times 10^{16} h_{70}^{-1/2} (<M>/\rmn{M}_{\odot})^{1/2}$cm, and constrain the power-law index to $0.77\leq\nu\leq2.67$ (95 per cent confidence range). These results can be used to constrain models of quasar accretion discs. As a example, we find that the accretion disc in MG 0414+0534 is statistically consistent with a Shakura-Sunyaev thin disc model.

We test statistically the hypothesis that radio pulsar glitches result from an avalanche process, in which angular momentum is transferred erratically from the flywheel-like superfluid in the star to the slowly decelerating, solid crust via spatially connected chains of local, impulsive, threshold-activated events, so that the system fluctuates around a self-organised critical state. Analysis of the glitch population (currently 285 events from 101 pulsars) demonstrates that the size distribution in individual pulsars is consistent with being scale invariant, as expected for an avalanche process. The waiting-time distribution is consistent with being exponential in seven out of nine pulsars where it can be measured reliably, after adjusting for observational limits on the minimum waiting time, as for a constant-rate Poisson process. PSR J0537$-$6910 and PSR J0835$-$4510 are the exceptions; their waiting-time distributions show evidence of quasiperiodicity. In each object, stationarity requires that the rate $\lambda$ equals $- \epsilon \dot{\nu} / <\Delta\nu>$, where $\dot{\nu}$ is the angular acceleration of the crust, $<\Delta\nu>$ is the mean glitch size, and $\epsilon\dot{\nu}$ is the relative angular acceleration of the crust and superfluid. There is no evidence that $\lambda$ changes monotonically with spin-down age. The rate distribution itself is fitted reasonably well by an exponential for $\lambda \geq 0.25 {\rm yr^{-1}}$. For $\lambda < 0.25 {\rm yr^{-1}}$, its exact form is unknown; the exponential overestimates the number of glitching pulsars observed at low $\lambda$, where the limited total observation time exercises a selection bias.

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.

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.

The gravitational waves generated during supermassive black hole (SMBH) coalescence are prime candidates for detection by the satellite LISA. We use the extended Press-Schechter formalism combined with empirically motivated estimates for the SMBH--dark matter halo mass relation and SMBH occupation fraction to estimate the maximum coalescence rate for major SMBH mergers. Assuming efficient binary coalescence, and guided by the lowest nuclear black hole mass inferred in local galactic bulges and nearby low-luminosity active galactic nuclei (10^5 Msun) we predict approximately 15 detections per year at a signal to noise greater than five, in each of the inspiral and ringdown phases. Rare coalescences between SMBHs having masses in excess of 10^7 Msun will be more readily detected via gravitational waves from the ringdown phase.

We present seven epochs of spectroscopy on the quadruply imaged quasar SDSS J1004+4112, spanning observed-frame time delays from 1 to 322 days. The spectra reveal differences in the emission lines between the lensed images. Specifically, component A showed a strong enhancement in the blue wings of several high-ionization lines relative to component B, which lasted at least 28 days (observed frame) then faded. Since the predicted time delay between A and B is <30 days, our time coverage suggests that the event was not intrinsic to the quasar. We attribute these variations to microlensing of part of the broad emission line region of the quasar, apparently resolving structure in the source plane on a scale of ~10^{16} cm at z=1.734. In addition, we observed smaller differences in the emission line profiles between components A and B that persisted throughout the time span, which may also be due to microlensing or millilensing. Further spectroscopic monitoring of this system holds considerable promise for resolving the structure of the broad emission line region in quasars.

The coalescence of massive black hole (BH) binaries due to galaxy mergers provides a primary source of low-frequency gravitational radiation detectable by pulsar timing measurements and by the proposed LISA (Laser Interferometry Space Antenna) observatory. We compute the expected gravitational radiation signal from sources at all redshifts by combining the predicted merger rate of galactic halos with recent measurements of the relation between BH mass, M_bh, and the velocity dispersion of its host galaxy, sigma. Our main findings are as follows: (i) the nHz frequency background is dominated by BH binaries at redshifts z<2, and existing limits from pulsar timing data place tight constraints on the allowed normalization and power-law slope of the M_bh-sigma relation or on the fraction of BH binaries that proceed to coalescence; (ii) more than half of all discrete mHz massive BH sources detectable by LISA are likely to originate at redshifts z>7; (iii) the number of LISA sources per unit redshift per year should drop substantially after reionization as long as BH formation is triggered by gas cooling in galaxies. Studies of the highest redshift sources among the few hundred detectable events per year, will provide unique information about the physics and history of black hole growth in galaxies.

We present a generalization of the concept of magnification bias for gravitationally-lensed quasars, in which the quasars are selected by flux in more than one wavelength band. To illustrate the principle, we consider the case of two-band selection, in which the fluxes in the two bands are uncorrelated, perfectly correlated, or correlated with scatter. For uncorrelated fluxes, we show that the previously-held result - that the bias is the product of the single-band biases - is generally false. We demonstrate some important properties of the multi-band magnification bias using model luminosity functions inspired by observed correlations among X-ray, optical, infrared and radio fluxes of quasars. In particular, the bias need not be an increasing function of each flux, and the bias can be extremely large for non-linear correlations. The latter fact may account for the high lensing rates found in some X-ray/optical and infrared/radio selected samples.

The fraction of quasar's and gamma-ray burst (GRB) afterglows that vary due to microlensing by the stellar populations of intervening elliptical/S0 galaxies is computed by combining the joint distribution of effective microlensing convergence (kappa) and shear (gamma) with microlensing magnification patterns. Microlensing is common in multiply imaged sources. We find that 1 in 3 multiply imaged quasars should vary by more than 0.5 magnitudes per decade due to microlensing, while 10% of macrolensed GRB afterglows should show a departure of more than 0.5 magnitudes from their intrinsic light-curve during the first 30 days. However microlensing by stars is rare in general, with only 1 source in ~500 varying by more than 0.5 magnitudes during the same periods. We find that most microlensing by stars will be observed in a regime where gamma>0.1. Thus point-mass lenses do not provide an adequate description for most microlensing events. If dark matter halos contain a large fraction of mass in compact objects, the fraction of microlensed (by 0.5 magnitudes) images rises significantly to ~1 in 10 for quasars and ~1 in 5 for GRB afterglows. Comparison of variability between macrolensed and normal quasar images, and a moderate number of well sampled GRB afterglow light-curves should therefore discover or refute the existence of stellar mass compact objects in galaxy halos. While microlensing results in departures of the distribution of magnifications from that of a smooth model, the effect on the macrolensing magnification bias for the discovery of lenses in quasar surveys is small. On the other hand, microlensing significantly broadens the distribution of macrolensed image flux ratios.

Published measurements of distances to detached eclipsing binaries in the Large Magellanic Cloud combine stellar surface areas (in absolute units) determined from photometric light and radial velocity curves with surface brightnesses from model atmospheres and observed spectra. The method does not require the stars to be normal or undistorted, and is not limited in its applicability to the well detached systems that have traditionally been considered. With this in mind we carry out semi-detached light curve solutions for SMC binaries discovered by the OGLE collaboration, identify candidates for distance estimation, and tabulate OGLE light curve solutions. We point out important advantages for semi-detached binaries as standard candles; (1) light curve solutions can be strengthened by exploiting lobe-filling configurations, (2) only single-lined spectra may be needed because the mass ratio can often be determined from photometry, and (3) nearly all semi-detached binaries have sensibly circular orbits. We carry out simulations with synthetic data to show that semi-detached binaries can often be reliably identified and to quantify the accuracy of solutions. The simulations demonstrate two additional advantages for semi-detached distance determination candidates; (4) the well-known difficulty in distinguishing solutions with interchanged radii (aliasing) is much less severe for semi-detached than for detached binaries, and (5) the desirable condition of complete eclipse is identified with improved reliability. In many cases we find that detached and semi-detached systems can be distinguished, and we select 36 candidate semi-detached systems from the OGLE SMC eclipsing binary catalog. We also note the near absence of Algol-like light curves in the OGLE sample and discuss possible explanations.

Published parametric models of the Einstein Cross gravitational lens demonstrate that the image geometry can be reproduced by families of models. In particular, the slope of the mass-profile for the lens galaxy is unconstrained. However, recent models predict a dependence of image flux ratios on the slope of the mass profile. We use this dependence to constrain the mass profile by calculating the likelihood of the slope using published mid-IR flux ratios (including microlensing variability). We find that the galaxy is likely to be flatter than isothermal, and therefore that the mass-to-light ratio is decreasing in the inner kpc.

We use published mid-IR and V-band flux ratios for images A and B of Q2237+0305 to demonstrate that the size of the mid-IR emission region has a scale comparable to or larger than the microlens Einstein Radius (ER, ~10^17 cm for solar mass stars). Q2237+0305 has been monitored extensively in the R and V-bands for ~15 years. The variability record shows significant microlensing variability of the optical emission region, and has been used by several studies to demonstrate that the optical emission region is much smaller than the ER for solar-mass objects. For the majority of the monitoring history, the optical flux ratios have differed significantly from those predicted by macro-models. In contrast, recent observations in mid-IR show flux ratios similar to those measured in the radio, and to predictions of some lens models, implying that the mid-IR flux is emitted from a region that is at least 2 orders of magnitude larger than the optical emission region. We have calculated the likeli-hood of the observed mid-IR flux ratio as a function of mid-IR source size given the observed V-band flux ratio. The expected flux ratio for a source having dimensions of ~1 ER is a sensitive function of the macro model adopted. However we find that the probability of source size given the observed flux ratios is primarily sensitive to the ratio of the macro-model magnifications. The majority of published macro models for Q2237+0305 yield a flux ratio for images B and A of 0.8 - 1.1. By combining probabilities from the ratios A/B and C/D we infer that the diameter of a circular IR emission region is >1ER with >95% confidence. For microlensing by low-mass stars, this source size limit rules out non-thermal processes such as synchrotron as mechanisms for mid-IR emission.

It is conventional to calculate the probability of microlensing for a cosmologically distant source based on the Press-Gunn approximation that the lensing objects are uniformly and randomly distributed in the intervening space with a constant comoving density. Here we investigate more realistic cosmological microlensing statistics by considering the strong spatial clustering of likely lensing objects with each other in galaxies and their association with the clumps of dark matter that make up the massive halos of galaxies. The distribution of microlensing optical depth (kappa) along randomly chosen sight lines is calculated as is the conditional distribution of kappa along sight lines near one which is strongly microlensed. Our overall result is that the Press-Gunn approximation is a useful order-of-magnitude approximation if the massive halos of galaxies are made of dark compact objects but that it fails badly and can be qualitatively misleading in the more likely case in which only the ordinary stellar populations of galaxies are the dominant source of cosmological microlensing events. In particular, we find that microlensing by stars is limited to of order 1 percent of high redshift sources at any one time. Furthermore, even though only a small fraction of high redshift sources are multiply-imaged (by galaxies), it is these sources that are most likely to be microlensed by stars. Consequently, microlensing by stars is usually observed at kappa's near 1 where the simple isolated point mass lens approximation is not appropriate. However, if CDM halos are composed of condensed objects, then more than 10 percent of high redshift sources are microlensed at any given time. The vast majority of these sources are not multiply-imaged, and have kappa's smaller than 0.01.

Strong lensing is a powerful probe of the distribution of matter in the cores of clusters of galaxies. Recent studies suggest that the cold dark matter model predicts cores that are denser than those observed in galaxies, groups and clusters. One possible resolution of the discrepancy is that the dark matter has strong interactions (SIDM), which leads to lower central densities. A generalized form of the Navarro, Frenk and White profile (Zhao profile) may be used to describe these halos. In this paper we examine gravitational lensing statistics for this class of model. The optical depth to multiple imaging is a very sensitive function of the profile parameters in the range of interest for SIDM halos around clusters of galaxies. Less concentrated profiles, which result from larger self-interaction cross-sections, can produce many fewer lensed pairs. Lensing statistics provide a powerful test for SIDM. More realistic and observationally oriented calculations remain to be done, however larger self-interaction cross-sections may well be ruled out by the very existence of strong lenses on galaxy cluster scales. The inclusion of centrally dominant cluster galaxies should boost the cross-section to multiple imaging. However our preliminary calculations suggest that the additional multiple imaging rate is small with respect to the differences in multiple imaging rate for different halo profiles. In future statistical studies, it will be important to properly account for the scatter among halo profiles since the optical depth to multiple imaging is dominated by the most concentrated members of a cluster population.

Detached eclipsing binary stars provide a robust one-step distance determination to nearby galaxies. As a by-product of Galactic microlensing searches, catalogs of thousands of variable stars including eclipsing binaries have been produced by the OGLE, MACHO and EROS collaborations. We present photometric solutions for detached eclipsing binaries in the Small Magellanic Cloud (SMC) discovered by the OGLE collaboration. The solutions were obtained with an automated version of the Wilson-Devinney program. By fitting mock catalogs of eclipsing binaries we find that the normalized stellar radii (particularly their sum) and the surface brightness ratio are accurately described by the fitted parameters and estimated standard errors, despite various systematic uncertainties. In many cases these parameters are well constrained. In addition we find that systems exhibiting complete eclipses can be reliably identified where the fractional standard errors in the radii are small. We present two quantitatively selected sub-samples of eclipsing binaries that will be excellent distance indicators. These can be used both for computation of the distance to the SMC and to probe its structure. One particularly interesting binary has a very well determined solution, exhibits complete eclipses, and is comprised of well detached G-type, class $II$ giants.

If cluster dark matter is in the form of compact objects it will introduce fluctuations into the light curves of distant sources. Current searches for MACHOs in clusters of galaxies focus on monitoring quasars behind nearby systems. This paper considers the effect of such a compact population on the surface brightness distribution of giant gravitationally lensed arcs. As the microlensing optical depth is significant in these clusters, the expected fluctuations are substantial and are observable. Focusing on the giant arc seen in Abell 370, we demonstrate that several `extreme' events would be visible in a comparison of HST observations at two epochs. Utilizing NGST, long term monitoring should reveal a ubiquitous twinkling of brightness over the surface of the arcs.

The first investigations of the response of the microlensing magnification pattern (at an optical depth of order unity) to the mass function of the microlenses found that the resulting statistics depend only on the mean microlens mass <m>. In particular the mean microlensing caustic crossing rate was found to be proportional to the squareroot of <m>. We show that while this is true in the limit of mass functions with a narrow range of mass, in general the magnification pattern shows structure that reflects the contribution to the optical depth of microlenses with different masses. We present a better approximation, relating the microlens mass function to light-curve statistics. We show that the variability statistics of quasar microlensing light-curves can (in principle) be inverted to obtain the mass function of the microlenses in the mass range over which the mass density remains comparable, ie. p(m)dm~C/m. A preliminary analysis of the structure function for Q2237+0305 suggests that there is not a significant contribution to the optical depth from very low mass objects (0.001 solar mass). However observations of multiple microlensed quasars for a period of ~20 years may in the future yield a detailed p(m)dm. In the mass range where the number density is comparable, ie. p(m)dm~const., the distribution of flux factors could be inverted to find the microlens mass function. This may be used as a probe of the abundance of planets with orbital radii >100 AU.

Gravitational lensing of a gamma ray burst (GRB) by a single point mass will produce a second, delayed signal. Several authors have discussed using microlensed GRBs to probe a possible cosmological population of compact objects. We analyse a closely related phenomenon; the effect of microlensing by low to medium optical depth in compact objects on the averaged observed light-curve of a sample of GRBs. We discuss the cumulative measured flux as a function of time resulting from delays due to microlensing by cosmological compact objects. The time-scale and curvature of this function describe unique values for the compact object mass and optical depth. For GRBs with durations larger than the detector resolution, limits could be placed on the mass and optical depth of cosmological compact objects. The method does not rely on the separation of lensed bursts from those which are spatially coincident.

The central regions of the gravitationally lensed quasar Q2237+0305 can be indirectly resolved on nano-arcsecond scales if viewed spectrophotometricly during a microlensing high magnification event (HME). Q2237+0305 is currently being monitored from the ground (eg. OGLE collaboration, Apache Point Observatory), with the goal, among others, of triggering ground and spacecraft based target of opportunity (TOO) observations of an HME. In this work we investigate the rate of change (trigger) in image brightness that signals an imminent HME and importantly, the separation between the trigger and the event peak. In addition, we produce colour dependent model light-curves by combining high-resolution microlensing simulations with a realistic model for a thermal accretion disc source. We make hypothetical target of opportunity spectroscopic observations using our determination of the appropriate trigger as a guide. We find that if the source spectrum varies with source radius, a 3 observation TOO program should be able to observe a microlensing change in the continuum slope following a light-curve trigger with a success rate of >80%.

Microlensing in Q2237+0305 between 1985 and 1995 (eg. Irwin et al. 1989; Corrigan et al. 1991; Ostensen et al. 1996) has been interpreted in two different ways; as microlensing by stellar mass objects of a continuum source having dimensions significantly smaller than the microlens Einstein radius (ER) (eg. Wambsganss, Paczynski & Schneider 1990; Rauch & Blandford 1991), and as microlensing by very low mass objects of a source as large as 5 ER (Refsdal & Stabell 1993; Haugan 1996). In this paper we present evidence in favour of a small source. Limits on the source size (in units of ER) are obtained from the combination of limits on the number of microlens Einstein radii crossed by the source during the monitoring period with two separate light-curve features. Firstly, recently published monitoring data (Wozniak et al. 2000; OGLE web page) show large variations (~0.8-1.5 magnitudes) between image brightnesses over a period of 700 days or ~15% of the monitoring period. Secondly, the 1988 peak in the image A light-curve had a duration that is a small fraction (<0.02) of the monitoring period. Such rapid microlensing rises and short microlensing peaks only occur for small sources. We find that the observed large-rapid variation limits the source size to be <0.2 ER (95% confidence). The width of the light-curve peak provides a stronger constraint of <0.02 ER (99% confidence). The Einstein radius (projected into the source plane) of the average microlens mass (m) in Q2237+0305 is ER ~ 10^{17}\sqrt{m} cm. The interpretation that stars are responsible for microlensing in Q2237+0305 therefore results in limits on the continuum source size that are consistent with current accretion disc theory.

Spectrophotometric observation of the gravitationally microlensed quasar Q2237+0305 during a High Magnification Event (HME) is potentially a very powerful tool for probing the structure of the quasars accretion disc on scales of less than 10^-8 arc seconds. How often we can expect to observe a HME is dependent on the lens system parameters of galactic transverse velocity, mean microlens mass and the size of the magnified continuum source. We have previously used published microlensed light-curves to obtain expressions for the likely-hood of the values for these parameters (Wyithe, Webster & Turner 1999b,2000; Wyithe, Webster, Turner & Mortlock 2000). Here we use this information to investigate the expected rate of SHMEs. We find the average rate of SHMEs as well as the number that we can expect to observe over periods of a decade and of a single observing season. We find that the average SHME rate summed over all images in Q2237+0305 is 1.5\pm0.6 - 6.3\pm1.3 events per decade. During the period following a caustic crossing we find that the event rate in the corresponding image is enhanced by 50-100%, and therefore that the overall event rate may be higher during these periods. From the distribution of events expected during a 6 month period we find that there is 1 chance in 4-10 of observing a SHME per observing season. The systematic dependence in these values arises from the different assumptions for smooth matter content, orientation of the galactic transverse velocity and the size of photometric error in the monitoring data. The results support continued monitoring of Q2237+0305 with the aim of obtaining detailed spectroscopic and photometric observations of a SHME.

The four bright images of the gravitationally lensed quasar Q2237+0305 are being monitored from the ground (eg. OGLE collaboration, Apache Point Observatory) in the hope of observing a high magnification event (HME). Over the past three seasons (1997-1999) the OGLE collaboration has produced microlensing light-curves with unprecedented coverage. These demonstrate smooth, independent (therefore microlensing) variability between the images (Wozniak et al. 2000a,b; OGLE web page). We have retrospectively compared probability functions for high-magnification event parameters with several observed light-curve features. We conclude that the 1999 image C peak was due to the source having passed outside of a cusp rather than to a caustic crossing. In addition, we find that the image C light-curve shows evidence for a caustic crossing between the 1997 and 1998 observing seasons involving the appearance of new critical images. Our models predict that the next image C event is most likely to arrive 500 days following the 1999 peak, but with a large uncertainty (100-2000 days). Finally, given the image A light-curve derivative at the end of the 1999 observing season, our modelling suggests that a caustic crossing will occur between the 1999 and 2000 observing seasons, implying a minimum for the image A light-curve ~1-1.5 magnitudes fainter than the November 1999 level.

Gravitational microlensing at cosmological distances is potentially a powerful tool for probing the mass functions of stars and compact objects in other galaxies. In the case of multiply-imaged quasars, microlensing data has been used to determine the average microlens mass. However the measurements have relied on an assumed transverse velocity for the lensing galaxy. Since the measured mass scales with the square of the transverse velocity, published mass limits are quite uncertain. In the case of Q2237+0305 we have properly constrained this uncertainty. The distribution of light curve derivatives allows quantitative treatment of the relative rates of microlensing due to proper motions of microlenses, the orbital stream motion of microlenses and the bulk galactic transverse velocity. By demanding that the microlensing rate due to the motions of microlenses is the minimum that should be observed we determine lower limits for the average mass of stars and compact objects in the bulge of Q2237+0305. If microlenses are assumed to move in an orbital stream the lower limit ranges between 0.005 and 0.023 solar masses where the the systematic dependence is due to the fraction of smooth matter and the size of photometric error assumed for published monitoring data. However, if the microlenses are assumed to move according to an isotropic velocity dispersion then a larger lower limit of 0.019-0.11 solar masses is obtained. A significant contribution of Jupiter mass compact objects to the mass distribution of the galactic bulge of Q2237+0305 is therefore unambiguously ruled out.

Following the detection of a gravitational microlensing high magnification event (HME) in Q2237+0305A attempts have been made to place limits on the dimensions of the quasar continuum source. The analyses have studied either the observed event magnitude or the event duration. The latter approach has been hampered by lack of knowledge about the transverse velocity of the lensing galaxy. We obtain both upper and lower statistical limits on the size of the continuum source from the observed HME using determinations of transverse velocity obtained from the published monitoring data. Our calculations take account of the caustic orientation as well as the component of the caustic velocity that results from stellar proper motions. Our determination of source size relies on an estimated duration of 52 days for the HME, and so will be refined when more HMEs are observed. We find that the upper and lower limits on the magnified region of the R-band continuum source are 6 times 10^15 and 2 times 10^13 cm respectively (99% confidence). Through consideration of the joint probability for source size and mean microlens mass we find that the mean mass lies between ~0.01 solar masses and ~1 solar mass (95% confidence).

We present a method for computing the probability distribution of microlensed light curve derivatives both in the case of a static lens with a transverse velocity, and in the case of microlensing that is produced through stellar proper motions. The distributions are closely related in form, and can be considered equivalent after appropriate scaling of the input transverse velocity. The comparison of the distributions in this manner provides a consistent way to consider the relative contribution to microlensing (both large and small fluctuations) of the two classes of motion, a problem that is otherwise an extremely expensive computational exercise. We find that the relative contribution of stellar proper motions to the microlensing rate is independent of the mass function assumed for the microlenses, but is a function of optical depth and shear. We find that stellar proper motions produce a higher overall microlensing rate than a transverse velocity of the same magnitude. This effect becomes more pronounced at higher optical depth. With the introduction of shear, the relative rates of microlensing become dependent on the direction of the transverse velocity. This may have important consequences in the case of quadruply lensed quasars such as Q2237+0305, where the alignment of the shear vector with the source trajectory varies between images.

The largest systematic uncertainty present in the analysis of gravitationally microlensed quasar light curves is that of the galactic transverse velocity. We describe a method for determining the transverse velocity as well as its application to published monitoring data of Q2237+0305 (Irwin et al. 1989; Corrigan et al. 1991; Ostensen et al 1995). We find that while this data displays strong evidence for microlensing, it limits the transverse velocity to be < 500 km/sec (for microlens masses of 0.5 solar masses). In addition, by combining these results with a new method for modelling the microlensing contribution of stellar proper motions, we find that the masses of microlenses in Q2237+0305 are consistent with those of low mass stars.

Determination of microlensing parameters in the gravitationally lensed quasar Q2237+0305 from the statistics of high magnification events will require monitoring for more than 100 years (Wambsganss, Paczynski & Schneider 1990). However we show that the effective transverse velocity of the lensing galaxy can be determined on a more realistic time-scale through consideration of the distribution of light-curve derivatives. The 10 years of existing monitoring data for Q2237+0305 are analysed. These data display strong evidence for microlensing that is not associated with a high magnification event. An upper limit of v < 500 km/sec is obtained for the galactic transverse velocity which is smaller than previously assumed values. The analysis suggests that the observed microlensing variation may be predominantly due to stellar proper motions. The statistical significance of the results obtained from our method will be increased by the addition of data points from current and future monitoring campaigns. However reduced photometric errors will be more valuable than an increased sampling rate.

The method devised by Lewis et al. (1993) for calculating the light curve of a microlensed point source is expanded to two dimensions to enable the calculation of light curves of extended sources. This method is significantly faster than the ray shooting method that has been used in the past. The increased efficiency is used to obtain much higher resolution light curves over increased timescales. We investigate the signatures arising from different source geometries in a realistic microlensing model. We show that a large fraction of high magnification events (HMEs) in image A of Q2237+0305 involve only one caustic, and could therefore yield information on the structure of the quasar continuum through the recognition of a characteristic event shape. In addition, the cataloguing of HMEs into morphological type will, in theory, enable the direction of the transverse motion, as well as the source size to be obtained from long term monitoring.