Source author record

J. S. B. Wyithe

J. S. B. Wyithe appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

39works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

39 published item(s)

preprint2021arXiv

The Prevalence of Galaxy Overdensities Around UV-Luminous Lyman $\mathbfα$ Emitters in the Epoch of Reionization

Before the end of the epoch of reionization, the Hydrogen in the Universe was predominantly neutral. This leads to a strong attenuation of Ly$α$ lines of $z\gtrsim6$ galaxies in the intergalactic medium. Nevertheless, Ly$α$ has been detected up to very high redshifts ($z\sim9$) for several especially UV luminous galaxies. Here, we test to what extent the galaxy's local environment might impact the Ly$α$ transmission of such sources. We present an analysis of dedicated Hubble Space Telescope (HST) imaging in the CANDELS/EGS field to search for fainter neighbours around three of the most UV luminous and most distant spectroscopically confirmed Ly$α$ emitters: EGS-zs8-1, EGS-zs8-2 and EGSY-z8p7 at $z_\mathrm{spec}=7.73$, 7.48, and 8.68, respectively. We combine the multi-wavelength HST imaging with Spitzer data to reliably select $z\sim7-9$ galaxies around the central, UV-luminous sources. In all cases, we find a clear enhancement of neighbouring galaxies compared to the expected number in a blank field (by a factor $\sim 3-9\times$). Our analysis thus reveals ubiquitous overdensities around luminous Ly$α$ emitting sources in the heart of the cosmic reionization epoch. We show that our results are in excellent agreement with expectations from the Dragons simulation, confirming the theoretical prediction that the first ionized bubbles preferentially formed in overdense regions. JWST follow-up observations of the neighbouring galaxies identified here will be needed to confirm their physical association and to map out the ionized regions produced by these sources.

preprint2020arXiv

Deep multi-redshift limits on Epoch of Reionisation 21cm Power Spectra from Four Seasons of Murchison Widefield Array Observations

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$σ$ 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.

preprint2016arXiv

A High Reliability Survey of Discrete Epoch of Reionization Foreground Sources in the MWA EoR0 Field

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.

preprint2016arXiv

CHIPS: The Cosmological HI Power Spectrum Estimator

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$σ$ upper limit on the EoR dimensionless power at $k=0.05$~h.Mpc$^{-1}$ of $Δ_k^2<7.6\times{10^4}$~mK$^2$ in the redshift range $z=[6.2-6.6]$, consistent with previous estimates.

preprint2016arXiv

Delay Spectrum with Phase-Tracking Arrays: Extracting the HI power spectrum from the Epoch of Reionization

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.

preprint2016arXiv

First Limits on the 21 cm Power Spectrum during the Epoch of X-ray heating

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.

preprint2016arXiv

First Season MWA EoR Power Spectrum Results at Redshift 7

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 $Δ^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.

preprint2016arXiv

Low frequency observations of linearly polarized structures in the interstellar medium near the south Galactic pole

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.

preprint2016arXiv

Parametrising Epoch of Reionization foregrounds: A deep survey of low-frequency point-source spectra with the MWA

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.

preprint2016arXiv

The Importance of Wide-field Foreground Removal for 21 cm Cosmology: A Demonstration With Early MWA Epoch of Reionization Observations

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.

preprint2016arXiv

The Murchison Widefield Array 21 cm Power Spectrum Analysis Methodology

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.

preprint2016arXiv

The relation between star formation rate and stellar mass of galaxies at z $\sim$ 1-4

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.

preprint2015arXiv

A spectroscopically confirmed z=1.327 galaxy-scale deflector magnifying a z~8 Lyman-Break galaxy in the Brightest of Reionizing Galaxies survey

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-$α$, H-$β$, [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σ$, 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érsic source magnified by a factor of $μ=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$.

preprint2015arXiv

Confirmation of Wide-Field Signatures in Redshifted 21 cm Power Spectra

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.

preprint2015arXiv

Foregrounds in Wide-Field Redshifted 21 cm Power Spectra

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.

preprint2015arXiv

Gravitational waves from binary supermassive black holes missing in pulsar observations

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.

preprint2015arXiv

Intensity Mapping Cross-Correlations: Connecting the Largest Scales to Galaxy Evolution

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.

preprint2015arXiv

The Murchison Widefield Array Correlator

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.

preprint2014arXiv

Binary super-massive black hole environments diminish the gravitational-wave signal in the pulsar timing band

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.

preprint2014arXiv

Measurement of Galaxy Clustering at z~7.2 and the Evolution of Galaxy Bias from 3.8<z<8 in the XDF, GOODS-S AND GOODS-N

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σ$, 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σ$ 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, $ε_{DC}$). At $z=7.2$ values near unity are preferred, which may be explained by the shortened halo assembly time at high-redshift.

preprint2014arXiv

Prospects for gravitational-wave detection and supermassive black hole astrophysics with pulsar timing arrays

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.

preprint2013arXiv

A 189 MHz, 2400 square degree polarization survey with the Murchison Widefield Array 32-element prototype

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 $(α,δ) = (4^h,-27^\circ.6)$ in terms of three dimensional power spectra

preprint2013arXiv

A survey for transients and variables with the Murchison Widefield Array 32-tile prototype at 154 MHz

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.

preprint2013arXiv

Gravitational-wave Limits from Pulsar Timing Constrain Supermassive Black Hole Evolution

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 $Ω_{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.

preprint2013arXiv

On the detection and tracking of space debris using the Murchison Widefield Array. I. Simulations and test observations demonstrate feasibility

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.

preprint2013arXiv

The giant lobes of Centaurus A observed at 118 MHz with the Murchison Widefield Array

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.

preprint2013arXiv

The Murchison Widefield Array: solar science with the low frequency SKA Precursor

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.

preprint2012arXiv

A new layout optimization technique for interferometric arrays, applied to the MWA

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.

preprint2012arXiv

Does a 'stochastic' background of gravitational waves exist in the pulsar timing band?

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.

preprint2012arXiv

Low Frequency Observations of the Moon with the Murchison Widefield Array

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.

preprint2012arXiv

Near-infrared observations of type Ia supernovae: The best known standard candle for cosmology

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.

preprint2012arXiv

On the role of feedback in shaping the cosmic abundance and clustering of neutral atomic hydrogen in galaxies

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.

preprint2012arXiv

Realisation of a low frequency SKA Precursor: The Murchison Widefield Array

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.

preprint2012arXiv

The EoR Sensitivity of the Murchison Widefield Array

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σ detection of the EoR signal with one full season of observation on two fields (900 and 700 hours).

preprint2012arXiv

The Murchison Widefield Array: the Square Kilometre Array Precursor at low radio frequencies

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.

preprint2011arXiv

A microlensing measurement of dark matter fractions in three lensing galaxies

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.

preprint2010arXiv

Subtraction of point sources from interferometric radio images through an algebraic forward modeling scheme

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.

preprint2000arXiv

The rate of caustic crossing microlensing events for Q2237+0305

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.

preprint1999arXiv

Determination of the transverse velocity of Q2237+0305

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.