The Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) is a next-generation survey of neutral hydrogen (HI) in the Local Universe. It uses the widefield, high-resolution capability of the Australian Square Kilometer Array Pathfinder (ASKAP), a radio interferometer consisting of 36 x 12-m dishes equipped with Phased-Array Feeds (PAFs), located in an extremely radio-quiet zone in Western Australia. WALLABY aims to survey three-quarters of the sky (-90 degr < Dec < +30 degr) to a redshift of z < 0.26, and generate spectral line image cubes at ~30 arcsec resolution and ~1.6 mJy/beam per 4 km/s channel sensitivity. ASKAP's instantaneous field of view at 1.4 GHz, delivered by the PAF's 36 beams, is about 30 sq deg. At an integrated signal-to-noise ratio of five, WALLABY is expected to detect over half a million galaxies with a mean redshift of z ~ 0.05 (~200 Mpc). The scientific goals of WALLABY include: (a) a census of gas-rich galaxies in the vicinity of the Local Group; (b) a study of the HI properties of galaxies, groups and clusters, in particular the influence of the environment on galaxy evolution; and (c) the refinement of cosmological parameters using the spatial and redshift distribution of low-bias gas-rich galaxies. For context we provide an overview of previous large-scale HI surveys. Combined with existing and new multi-wavelength sky surveys, WALLABY will enable an exciting new generation of panchromatic studies of the Local Universe. - First results from the WALLABY pilot survey are revealed, with initial data products publicly available in the CSIRO ASKAP Science Data Archive (CASDA).

The WiggleZ Dark Energy Survey measured the redshifts of over 200,000 UV-selected (NUV<22.8 mag) galaxies on the Anglo-Australian Telescope. The survey detected the baryon acoustic oscillation signal in the large scale distribution of galaxies over the redshift range 0.2<z<1.0, confirming the acceleration of the expansion of the Universe and measuring the rate of structure growth within it. Here we present the final data release of the survey: a catalogue of 225415 galaxies and individual files of the galaxy spectra. We analyse the emission-line properties of these UV-luminous Lyman-break galaxies by stacking the spectra in bins of luminosity, redshift, and stellar mass. The most luminous (-25 mag < MFUV <-22 mag) galaxies have very broad H-beta emission from active nuclei, as well as a broad second component to the [OIII] (495.9 nm, 500.7 nm) doublet lines that is blue shifted by 100 km/s, indicating the presence of gas outflows in these galaxies. The composite spectra allow us to detect and measure the temperature-sensitive [OIII] (436.3 nm) line and obtain metallicities using the direct method. The metallicities of intermediate stellar mass (8.8<log(M*/Msun)<10) WiggleZ galaxies are consistent with normal emission-line galaxies at the same masses. In contrast, the metallicities of high stellar mass (10<log(M*/Msun)<12) WiggleZ galaxies are significantly lower than for normal emission-line galaxies at the same masses. This is not an effect of evolution as the metallicities do not vary with redshift; it is most likely a property specific to the extremely UV-luminous WiggleZ galaxies.

We present the results of a deep survey of the nearby Sculptor group and the associated Sculptor filament taken with the Parkes 64-m radio telescope in the 21-cm emission line of neutral hydrogen. We detect 31 HI sources in the Sculptor group/filament, eight of which are new HI detections. We derive a slope of the HI mass function along the Sculptor filament of $\alpha = -1.10^{+0.20}_{-0.11}$, which is significantly flatter than the global mass function and consistent with the flat slopes previously found in other low-density group environments. Some physical process, such as star formation, photoionisation or ram-pressure stripping, must therefore be responsible for removing neutral gas predominantly from low-mass galaxies. All of our HI detections have a confirmed or tentative optical counterpart and are likely associated with luminous rather than 'dark' galaxies. Despite a column density sensitivity of about $4 \times 10^{17}~\mathrm{cm}^{-2}$, we do not find any traces of extragalactic gas or tidal streams, suggesting that the Sculptor filament is, at the current time, a relatively quiescent environment that has not seen any recent major interactions or mergers.

We present results from the 2D anisotropic Baryon Acoustic Oscillation (BAO) signal present in the final dataset from the WiggleZ Dark Energy Survey. We analyse the WiggleZ data in two ways: firstly using the full shape of the 2D correlation function and secondly focussing only on the position of the BAO peak in the reconstructed data set. When fitting for the full shape of the 2D correlation function we use a multipole expansion to compare with theory. When we use the reconstructed data we marginalise over the shape and just measure the position of the BAO peak, analysing the data in wedges separating the signal along the line of sight from that parallel to the line of sight. We verify our method with mock data and find the results to be free of bias or systematic offsets. We also redo the pre-reconstruction angle averaged (1D) WiggleZ BAO analysis with an improved covariance and present an updated result. The final results are presented in the form of $\Omega_c h^2$, $H(z)$, and $D_A(z)$ for three redshift bins with effective redshifts $z = 0.44$, $0.60$, and $0.73$. Within these bins and methodologies, we recover constraints between 5% and 22% error. Our cosmological constraints are consistent with Flat $\Lambda$CDM cosmology and agree with results from the Baryon Oscillation Spectroscopic Survey (BOSS).

We present the Large Area Radio Galaxy Evolution Spectroscopic Survey (LARGESS), a spectroscopic catalogue of radio sources designed to include the full range of radio AGN populations out to redshift z = 0.8. The catalogue covers roughly 800 square degrees of sky, and provides optical identifications for 19,179 radio sources from the 1.4 GHz Faint Images of the Radio Sky at Twenty-cm (FIRST) survey down to an optical magnitude limit of i_mod < 20.5 in Sloan Digital Sky Survey (SDSS) images. Both galaxies and point-like objects are included, and no colour cuts are applied. In collaboration with the WiggleZ and Galaxy And Mass Assembly (GAMA) spectroscopic survey teams, we have obtained new spectra for over 5,000 objects in the LARGESS sample. Combining these new spectra with data from earlier surveys provides spectroscopic data for 12,329 radio sources in the survey area, of which 10,856 have reliable redshifts. 85% of the LARGESS spectroscopic sample are radio AGN (median redshift z = 0.44), and 15% are nearby star-forming galaxies (median z = 0.08). Low-excitation radio galaxies (LERGs) comprise the majority (83%) of LARGESS radio AGN at z < 0.8, with 12% being high-excitation radio galaxies (HERGs) and 5% radio-loud QSOs. Unlike the more homogeneous LERG and QSO sub-populations, HERGs are a heterogeneous class of objects with relatively blue optical colours and a wide dispersion in mid-infrared colours. This is consistent with a picture in which most HERGs are hosted by galaxies with recent or ongoing star formation as well as a classical accretion disk.

We present radio Active Galactic Nuclei (AGN) luminosity functions over the redshift range 0.005 < z < 0.75. The sample from which the luminosity functions are constructed is an optical spectroscopic survey of radio galaxies, identified from matched Faint Images of the Radio Sky at Twenty-cm survey (FIRST) sources and Sloan Digital Sky Survey (SDSS) images.The radio AGN are separated into Low Excitation Radio Galaxies (LERGs) and High Excitation Radio Galaxies (HERGs) using the optical spectra. We derive radio luminosity functions for LERGs and HERGs separately in the three redshift bins (0.005 < z < 0.3, 0.3 < z < 0.5 and 0.5 < z <0.75). The radio luminosity functions can be well described by a double power-law. Assuming this double power-law shape the LERG population displays little or no evolution over this redshift range evolving as ~$(1+z)^{0.06}$ assuming pure density evolution or ~ $(1+z)^{0.46}$ assuming pure luminosity evolution. In contrast, the HERG population evolves more rapidly, best fitted by ~$(1+z)^{2.93}$ assuming a double power-law shape and pure density evolution. If a pure luminosity model is assumed the best fitting HERG evolution is parameterised by ~$(1+z)^{7.41}$. The characteristic break in the radio luminosity function occurs at a significantly higher power (~1 dex) for the HERG population in comparison to the LERGs. This is consistent with the two populations representing fundamentally different accretion modes.

We present a new redshift survey, the 2dF Quasar Dark Energy Survey pilot (2QDESp), which consists of ${\approx}10000$ quasars from ${\approx}150$ deg$^2$ of the southern sky, based on VST-ATLAS imaging and 2dF/AAOmega spectroscopy. Combining our optical photometry with the WISE (W1,W2) bands we can select essentially contamination free quasar samples with $0.8{<}z{<}2.5$ and $g{<}20.5$. At fainter magnitudes, optical UVX selection is still required to reach our $g{\approx}22.5$ limit. Using both these techniques we observed quasar redshifts at sky densities up to $90$ deg$^{-2}$. By comparing 2QDESp with other surveys (SDSS, 2QZ and 2SLAQ) we find that quasar clustering is approximately luminosity independent, with results for all four surveys consistent with a correlation scale of $r_{0}{=}6.1{\pm}0.1 \: h^{-1}$Mpc, despite their decade range in luminosity. We find a significant redshift dependence of clustering, particularly when BOSS data with $r_{0}{=}7.3{\pm}0.1 \: h^{-1}$Mpc are included at $z{\approx}2.4$. All quasars remain consistent with having a single host halo mass of ${\approx}2{\pm}1{\times}10^{12} \: h^{-1}M_\odot$. This result implies that either quasars do not radiate at a fixed fraction of the Eddington luminosity or AGN black hole and dark matter halo masses are weakly correlated. No significant evidence is found to support fainter, X-ray selected quasars at low redshift having larger halo masses as predicted by the `hot halo' mode AGN model of Fanidakis et al. 2013. Finally, although the combined quasar sample reaches an effective volume as large as that of the original SDSS LRG sample, we do not detect the BAO feature in these data.

We place the most robust constraint to date on the scale of the turnover in the cosmological matter power spectrum using data from the WiggleZ Dark Energy Survey. We find this feature to lie at a scale of $k_0=0.0160^{+0.0041}_{-0.0035}$ [h/Mpc] (68% confidence) for an effective redshift of 0.62 and obtain from this the first-ever turnover-derived distance and cosmology constraints: a measure of the cosmic distance-redshift relation in units of the horizon scale at the redshift of radiation-matter equality (r_H) of D_V(z=0.62)/r_H=18.3 (+6.3/-3.3) and, assuming a prior on the number of extra relativistic degrees of freedom $N_{eff}=3$, constraints on the matter density parameter $\Omega_Mh^2=0.136^{+0.026}_{-0.052}$ and on the redshift of matter-radiation equality $z_{eq}=3274^{+631}_{-1260}$. All results are in excellent agreement with the predictions of standard LCDM models. Our constraints on the logarithmic slope of the power spectrum on scales larger than the turnover is bounded in the lower limit with values only as low as -1 allowed, with the prediction of standard LCDM models easily accommodated by our results. Lastly, we generate forecasts for the achievable precision of future surveys at constraining $k_0$, $\Omega_Mh^2$, $z_{eq}$ and $N_{eff}$. We find that BOSS should substantially improve upon the WiggleZ turnover constraint, reaching a precision on $k_0$ of $\pm$9% (68% confidence), translating to precisions on $\Omega_Mh^2$ and $z_{eq}$ of $\pm$10% (assuming a prior $N_{eff}=3$) and on $N_{eff}$ of (+78/-56)% (assuming a prior $\Omega_Mh^2=0.135$). This is sufficient precision to sharpen the constraints on $N_{eff}$ from WMAP, particularly in its upper limit. For Euclid, we find corresponding attainable precisions on $(k_0, \Omega_Mh^2, N_eff)$ of (3,4,+17/-21)%. This represents a precision approaching our forecasts for the Planck Surveyor.

The Evolutionary Map of the Universe (EMU) is a proposed radio continuum survey of the Southern Hemisphere up to declination +30 deg., with the Australian Square Kilometre Array Pathfinder (ASKAP). EMU will use an automated source identification and measurement approach that is demonstrably optimal, to maximise the reliability, utility and robustness of the resulting radio source catalogues. As part of the process of achieving this aim, a "Data Challenge" has been conducted, providing international teams the opportunity to test a variety of source finders on a set of simulated images. The aim is to quantify the accuracy of existing automated source finding and measurement approaches, and to identify potential limitations. The Challenge attracted nine independent teams, who tested eleven different source finding tools. In addition, the Challenge initiators also tested the current ASKAPsoft source-finding tool to establish how it could benefit from incorporating successful features of the other tools. Here we present the results of the Data Challenge, identifying the successes and limitations for this broad variety of the current generation of radio source finding tools. As expected, most finders demonstrate completeness levels close to 100% at 10sigma dropping to levels around 10% by 5sigma. The reliability is typically close to 100% at 10sigma, with performance to lower sensitivities varying greatly between finders. All finders demonstrate the usual trade-off between completeness and reliability, whereby maintaining a high completeness at low signal-to-noise comes at the expense of reduced reliability, and vice-versa. We conclude with a series of recommendations for improving the performance of the ASKAPsoft source-finding tool.

We use a new, improved version of the HI Parkes All-Sky Survey to search for HI emission from nine new, ultra-faint Milky Way satellite galaxy candidates recently discovered in data from the Dark Energy Survey. None of the candidates is detected in HI, implying upper limits for their HI masses of typically several hundred to a few thousand solar masses. The resulting upper limits on M_HI / L_V and M_HI / M_star suggest that at least some of the new galaxy candidates are HI deficient. This finding is consistent with the general HI deficiency of satellite galaxies located within the Milky Way's virial radius and supports the hypothesis that gas is being removed from satellites by tidal and ram-pressure forces during perigalactic passages. In addition, some of the objects may be embedded in, and interacting with, the extended neutral and ionised gas filaments of the Magellanic Stream.

We introduce SoFiA, a flexible software application for the detection and parameterization of sources in 3D spectral-line datasets. SoFiA combines for the first time in a single piece of software a set of new source-finding and parameterization algorithms developed on the way to future HI surveys with ASKAP (WALLABY, DINGO) and APERTIF. It is designed to enable the general use of these new algorithms by the community on a broad range of datasets. The key advantages of SoFiA are the ability to: search for line emission on multiple scales to detect 3D sources in a complete and reliable way, taking into account noise level variations and the presence of artefacts in a data cube; estimate the reliability of individual detections; look for signal in arbitrarily large data cubes using a catalogue of 3D coordinates as a prior; provide a wide range of source parameters and output products which facilitate further analysis by the user. We highlight the modularity of SoFiA, which makes it a flexible package allowing users to select and apply only the algorithms useful for their data and science questions. This modularity makes it also possible to easily expand SoFiA in order to include additional methods as they become available. The full SoFiA distribution, including a dedicated graphical user interface, is publicly available for download.

The lenticular galaxy ESO 243-49 hosts the ultraluminous X-ray source HLX-1, the best candidate intermediate mass black hole (IMBH) currently known. The environments of IMBHs remain unknown, however the proposed candidates include the nuclei of dwarf galaxies or globular clusters. Evidence at optical wavelengths points at HLX-1 being the remnant of an accreted dwarf galaxy. Here we report the Australia Telescope Compact Array radio observations of HI emission in and around ESO 243-49 searching for signatures of a recent merger event. No HI line emission is detected in ESO 243-49 with a 5$\sigma$ upper limit on the HI gas mass of a few $10^8 M_{\odot}$. A likely reason for this non-detection is the cluster environment depleting ESO 243-49's HI gas reservoir. The upper limit is consistent with an interpretation of HLX-1 as a dwarf satellite of ESO 243-49, however more sensitive observations are required for a detection. We detect ~$5 \times 10^8 M_{\odot}$ of HI gas in the peculiar spiral galaxy AM 0108-462, located at a projected distance of ~170 kpc from ESO 243-49. This amount of HI gas is ~10 times less than in spiral galaxies with similar optical and near-infrared properties in the field, strengthening the conclusion that the cluster environment indeed depletes the HI gas reservoir of these two galaxies. Here we also report observations of AM 0108-462 in several optical and near-infrared bands using the Magellan 6.5 m telescopes, and archival X-ray and ultraviolet observations with XMM-Newton and Swift. These data combined with the HI line data suggest it is likely that AM 0108-462 is experiencing a merger event.

We present significant improvements in cosmic distance measurements from the WiggleZ Dark Energy Survey, achieved by applying the reconstruction of the baryonic acoustic feature technique. We show using both data and simulations that the reconstruction technique can often be effective despite patchiness of the survey, significant edge effects and shot-noise. We investigate three redshift bins in the redshift range 0.2<$z$<1, and in all three find improvement after reconstruction in the detection of the baryonic acoustic feature and its usage as a standard ruler. We measure model independent distance measures $D_{\mathrm V}(r_{\mathrm s}^\mathrm{fid}/r_{\mathrm s})$ of 1716 $\pm$ 83 Mpc, 2221 $\pm$ 101 Mpc, 2516 $\pm$ 86 Mpc (68% CL) at effective redshifts z = 0.44, 0.6, 0.73, respectively, where $D_{\mathrm V}$ is the volume-average-distance, and $r_{\mathrm s}$ is the sound horizon at the end of the baryon drag epoch. These significantly improved 4.8, 4.5 and 3.4 percent accuracy measurements are equivalent to those expected from surveys with up to 2.5 times the volume of WiggleZ. These measurements are fully consistent with cosmologies allowed by the analyses of the Planck Collaboration and the Sloan Digital Sky Survey.We provide the $D_{\mathrm V}(r_{\mathrm s}^\mathrm{fid}/r_{\mathrm s})$ posterior probability distributions and their covariances. When combining these measurements with temperature fluctuations measurements of Planck, the polarization of WMAP9, and the 6dF Galaxy Survey baryonic acoustic feature, we do not detect deviations from a flat LCDM model. Assuming this model we constrain the current expansion rate to $H_0$ = 67.15 $\pm$ 0.98 kms$^{-1}$Mpc$^{-1}$. Allowing the equation of state of dark energy to vary we obtain $w_\mathrm{DE}$ = -1.080 $\pm$ 0.135. When assuming a curved LCDM model we obtain a curvature value of $\Omega_{\mathrm K}$ = -0.0043 $\pm$ 0.0047.

Accurate parametrization of galaxies detected in the 21-cm HI emission is of fundamental importance to the measurement of commonly used indicators of galaxy evolution, including the Tully-Fisher relation and the HI mass function. Here, we propose a new analytic function, named the 'busy function', that can be used to accurately describe the characteristic double-horn HI profile of many galaxies. The busy function is a continuous, differentiable function that consists of only two basic functions, the error function, erf(x), and a polynomial, |x|^n, of degree n >= 2. We present the basic properties of the busy function and illustrate its great flexibility in fitting a wide range of HI profiles from the Gaussian profiles of dwarf galaxies to the broad, asymmetric double-horn profiles of spiral galaxies. Applications of the busy function include the accurate and efficient parametrization of observed HI spectra of galaxies and the construction of spectral templates for simulations and matched filtering algorithms. We demonstrate the busy function's power by automatically fitting it to the HI spectra of 1000 galaxies from the HIPASS Bright Galaxy Catalog, using our own C/C++ implementation, and comparing the resulting parameters with the catalogued ones. We also present two methods for determining the uncertainties of observational parameters derived from the fit.

We present the ultraviolet (UV) luminosity function of galaxies from the GALEX Medium Imaging Survey with measured spectroscopic redshifts from the first data release of the WiggleZ Dark Energy Survey. This sample selects galaxies with high star formation rates: at 0.6 < z < 0.9 the median star formation rate is at the upper 95th percentile of optically-selected (r<22.5) galaxies and the sample contains about 50 per cent of all NUV < 22.8, 0.6 < z < 0.9 starburst galaxies within the volume sampled. The most luminous galaxies in our sample (-21.0>M_NUV>-22.5) evolve very rapidly with a number density declining as (1+z)^{5\pm 1} from redshift z = 0.9 to z = 0.6. These starburst galaxies (M_NUV<-21 is approximately a star formation rate of 30 \msuny) contribute about 1 per cent of cosmic star formation over the redshift range z=0.6 to z=0.9. The star formation rate density of these very luminous galaxies evolves rapidly, as (1+z)^{4\pm 1}. Such a rapid evolution implies the majority of star formation in these large galaxies must have occurred before z = 0.9. We measure the UV luminosity function in 0.05 redshift intervals spanning 0.1<z<0.9, and provide analytic fits to the results. At all redshifts greater than z=0.55 we find that the bright end of the luminosity function is not well described by a pure Schechter function due to an excess of very luminous (M_NUV<-22) galaxies. These luminosity functions can be used to create a radial selection function for the WiggleZ survey or test models of galaxy formation and evolution. Here we test the AGN feedback model in Scannapieco et al. (2005), and find that this AGN feedback model requires AGN feedback efficiency to vary with one or more of the following: stellar mass, star formation rate and redshift.

Higher-order statistics are a useful and complementary tool for measuring the clustering of galaxies, containing information on the non-gaussian evolution and morphology of large-scale structure in the Universe. In this work we present measurements of the three-point correlation function (3PCF) for 187,000 galaxies in the WiggleZ spectroscopic galaxy survey. We explore the WiggleZ 3PCF scale and shape dependence at three different epochs z=0.35, 0.55 and 0.68, the highest redshifts where these measurements have been made to date. Using N-body simulations to predict the clustering of dark matter, we constrain the linear and non-linear bias parameters of WiggleZ galaxies with respect to dark matter, and marginalise over them to obtain constraints on sigma_8(z), the variance of perturbations on a scale of 8 Mpc/h and its evolution with redshift. These measurements of sigma_8(z), which have 10-20% accuracies, are consistent with the predictions of the LCDM concordance cosmology and test this model in a new way.

We characterise the stellar masses and star formation rates in a sample of almost 40000 spectroscopically confirmed UV luminous galaxies at 0.3<z<1.0 selected from within the WiggleZ Dark Energy Survey. In particular, we match this UV bright population to wide-field infrared surveys such as the near infrared UKIDSS LAS and the mid infrared WISE All-Sky Survey. We find that ~30% of the UV luminous WiggleZ galaxies are detected at >5sigma in the UKIDSS-LAS at all redshifts. An even more luminous subset of 15% are also detected in the WISE 3.4 and 4.6um bands. We compute stellar masses for this very large sample of extremely blue galaxies and quantify the sensitivity of the stellar mass estimates to various assumptions made during the SED fitting. The median stellar masses are log10(M*/M0)=9.6\pm0.7, 10.2\pm0.5 and 10.4\pm0.4 for the IR-undetected, UKIDSS detected and UKIDSS+WISE detected galaxies respectively. We demonstrate that the inclusion of NIR photometry can lead to tighter constraints on the stellar masses. The mass estimates are found to be most sensitive to the inclusion of secondary bursts of star formation as well as changes in the stellar population synthesis models, both of which can lead to median discrepancies of the order of 0.3dex in the stellar masses. We find that the best-fit M/LK is significantly lower (by ~0.4 dex) than that predicted by simple optical colour based estimators, in particular for the bluer galaxies with younger best-fit ages. The WiggleZ galaxies have star formation rates of 3-10 M0/yr and mostly lie at the upper end of the main sequence of star-forming galaxies at these redshifts. Their rest-frame UV luminosities and stellar masses are comparable to both local compact UV-luminous galaxies as well as Lyman break galaxies at z~2-3.(abridged)

The growth history of large-scale structure in the Universe is a powerful probe of the cosmological model, including the nature of dark energy. We study the growth rate of cosmic structure to redshift $z = 0.9$ using more than $162{,}000$ galaxy redshifts from the WiggleZ Dark Energy Survey. We divide the data into four redshift slices with effective redshifts $z = [0.2,0.4,0.6,0.76]$ and in each of the samples measure and model the 2-point galaxy correlation function in parallel and transverse directions to the line-of-sight. After simultaneously fitting for the galaxy bias factor we recover values for the cosmic growth rate which are consistent with our assumed $\Lambda$CDM input cosmological model, with an accuracy of around 20% in each redshift slice. We investigate the sensitivity of our results to the details of the assumed model and the range of physical scales fitted, making close comparison with a set of N-body simulations for calibration. Our measurements are consistent with an independent power-spectrum analysis of a similar dataset, demonstrating that the results are not driven by systematic errors. We determine the pairwise velocity dispersion of the sample in a non-parametric manner, showing that it systematically increases with decreasing redshift, and investigate the Alcock-Paczynski effects of changing the assumed fiducial model on the results. Our techniques should prove useful for current and future galaxy surveys mapping the growth rate of structure using the 2-dimensional correlation function.

This paper presents cosmological results from the final data release of the WiggleZ Dark Energy Survey. We perform full analyses of different cosmological models using the WiggleZ power spectra measured at z=0.22, 0.41, 0.60, and 0.78, combined with other cosmological datasets. The limiting factor in this analysis is the theoretical modelling of the galaxy power spectrum, including non-linearities, galaxy bias, and redshift-space distortions. In this paper we assess several different methods for modelling the theoretical power spectrum, testing them against the Gigaparsec WiggleZ simulations (GiggleZ). We fit for a base set of 6 cosmological parameters, {Omega_b h^2, Omega_CDM h^2, H_0, tau, A_s, n_s}, and 5 supplementary parameters {n_run, r, w, Omega_k, sum m_nu}. In combination with the Cosmic Microwave Background (CMB), our results are consistent with the LambdaCDM concordance cosmology, with a measurement of the matter density of Omega_m =0.29 +/- 0.016 and amplitude of fluctuations sigma_8 = 0.825 +/- 0.017. Using WiggleZ data with CMB and other distance and matter power spectra data, we find no evidence for any of the extension parameters being inconsistent with their LambdaCDM model values. The power spectra data and theoretical modelling tools are available for use as a module for CosmoMC, which we here make publicly available at http://smp.uq.edu.au/wigglez-data . We also release the data and random catalogues used to construct the baryon acoustic oscillation correlation function.

We have made the largest-volume measurement to date of the transition to large-scale homogeneity in the distribution of galaxies. We use the WiggleZ survey, a spectroscopic survey of over 200,000 blue galaxies in a cosmic volume of ~1 (Gpc/h)^3. A new method of defining the 'homogeneity scale' is presented, which is more robust than methods previously used in the literature, and which can be easily compared between different surveys. Due to the large cosmic depth of WiggleZ (up to z=1) we are able to make the first measurement of the transition to homogeneity over a range of cosmic epochs. The mean number of galaxies N(<r) in spheres of comoving radius r is proportional to r^3 within 1%, or equivalently the fractal dimension of the sample is within 1% of D_2=3, at radii larger than 71 \pm 8 Mpc/h at z~0.2, 70 \pm 5 Mpc/h at z~0.4, 81 \pm 5 Mpc/h at z~0.6, and 75 \pm 4 Mpc/h at z~0.8. We demonstrate the robustness of our results against selection function effects, using a LCDM N-body simulation and a suite of inhomogeneous fractal distributions. The results are in excellent agreement with both the LCDM N-body simulation and an analytical LCDM prediction. We can exclude a fractal distribution with fractal dimension below D_2=2.97 on scales from ~80 Mpc/h up to the largest scales probed by our measurement, ~300 Mpc/h, at 99.99% confidence.

The Square Kilometre Array and its pathfinders ASKAP and MeerKAT will produce prodigious amounts of data that necessitate automated source finding. The performance of automated source finders can be improved by pre-processing a dataset. In preparation for the WALLABY and DINGO surveys, we have used a test HI datacube constructed from actual Westerbork Telescope noise and WHISP HI galaxies to test the real world improvement of linear smoothing, the {\sc Duchamp} source finder's wavelet de-noising, iterative median smoothing and mathematical morphology subtraction, on intensity threshold source finding of spectral line datasets. To compare these pre-processing methods we have generated completeness-reliability performance curves for each method and a range of input parameters. We find that iterative median smoothing produces the best source finding results for ASKAP HI spectral line observations, but wavelet de-noising is a safer pre-processing technique. In this paper we also present our implementations of iterative median smoothing and mathematical morphology subtraction.

We perform a joint determination of the distance-redshift relation and cosmic expansion rate at redshifts z = 0.44, 0.6 and 0.73 by combining measurements of the baryon acoustic peak and Alcock-Paczynski distortion from galaxy clustering in the WiggleZ Dark Energy Survey, using a large ensemble of mock catalogues to calculate the covariance between the measurements. We find that D_A(z) = (1205 +/- 114, 1380 +/- 95, 1534 +/- 107) Mpc and H(z) = (82.6 +/- 7.8, 87.9 +/- 6.1, 97.3 +/- 7.0) km/s/Mpc at these three redshifts. Further combining our results with other baryon acoustic oscillation and distant supernovae datasets, we use a Monte Carlo Markov Chain technique to determine the evolution of the Hubble parameter H(z) as a stepwise function in 9 redshift bins of width dz = 0.1, also marginalizing over the spatial curvature. Our measurements of H(z), which have precision better than 7% in most redshift bins, are consistent with the expansion history predicted by a cosmological-constant dark-energy model, in which the expansion rate accelerates at redshift z < 0.7.

The absolute neutrino mass scale is currently unknown, but can be constrained from cosmology. The WiggleZ high redshift star-forming blue galaxy sample is less sensitive to systematics from non-linear structure formation, redshift-space distortions and galaxy bias than previous surveys. We obtain a upper limit on the sum of neutrino masses of 0.60eV (95% confidence) for WiggleZ+Wilkinson Microwave Anisotropy Probe. Combining with priors on the Hubble Parameter and the baryon acoustic oscillation scale gives an upper limit of 0.29eV, which is the strongest neutrino mass constraint derived from spectroscopic galaxy redshift surveys.

The large size of the ASKAP HI surveys DINGO and WALLABY necessitates automated 3D source finding. A performance difference of a few percent corresponds to a significant number of galaxies being detected or undetected. As such, the performance of the automated source finding is of paramount importance to both of these surveys. We have analysed the performance of various source finders to determine which will allow us to meet our survey goals during the DINGO and WALLABY design studies. Here we present a comparison of the performance of five different methods of automated source finding. These source finders are Duchamp, the Gamma-finder, CNHI, a 2D-1D Wavelet Reconstruction and S+C finder, a sigma clipping method. Each source finder was applied on the same three-dimensional data cubes containing (a) point sources with a Gaussian velocity profile and (b) spatially extended model-galaxies with inclinations and rotation profiles. We focus on the completeness and reliability of each algorithm when comparing the performance of the different source finders.

We study the evolution of galaxy populations around the spectroscopic WiggleZ sample of starforming galaxies at 0.25 < z < 0.75 using the photometric catalog from the Second Red-Sequence Cluster Survey (RCS2). We probe the optical photometric properties of the net excess neighbor galaxies. The key concept is that the marker galaxies and their neighbors are located at the same redshift, providing a sample of galaxies representing a complete census of galaxies in the neighborhood of star-forming galaxies. The results are compared with those using the RCS WiggleZ Spare-Fibre (RCS-WSF) sample as markers, representing galaxies in cluster environments at 0.25 < z < 0.45. By analyzing the stacked color-color properties of the WiggleZ neighbor galaxies, we find that their optical colors are not a strong function of indicators of star-forming activities such as EW([OII]) or GALEX NUV luminoisty of the markers. The galaxies around the WiggleZ markers exhibit a bimodal distribution on the color-magnitude diagram, with most of them located in the blue cloud. The optical galaxy luminosity functions (GLF) of the blue neighbor galaxies have a faint-end slope \alpha of \sim -1.3, similar to that for galaxies in cluster environments drawn from the RCS-WSF sample. The faint-end slope of the GLF for the red neighbors, however, is \sim -0.4, significantly shallower than the \sim -0.7 found for those in cluster environments. This suggests that the build-up of the faint-end of the red sequence in cluster environments is in a significantly more advanced stage than that in the star-forming and lower galaxy density WiggleZ neighborhoods. We find that the red galaxy fraction (fred) around the star-forming WiggleZ galaxies has similar values from z \sim 0.3 to z \sim 0.6 with fred \sim 0.28, but drops to fred \sim 0.20 at z > \sim0.7. This change of fred with redshift suggests that (and more...)

We describe a simple method to determine the reliability of source finders based on the detection of sources with both positive and negative total flux. Under the assumption that the noise is symmetric and that real sources have positive total flux, negative detections can be used to assign to each positive detection a probability of being real. We discuss this method in the context of upcoming, interferometric HI surveys.

The spectral line datacubes obtained from the Square Kilometre Array (SKA) and its precursors, such as the Australian SKA Pathfinder (ASKAP), will be sufficiently large to necessitate automated detection and parametrisation of sources. Matched filtering is widely acknowledged as the best possible method for the automated detection of sources. This paper presents the Characterised Noise Hi (CNHI) source finder, which employs a novel implementation of matched filtering. This implementation is optimised for the 3-D nature of the planned Wide-field ASKAP Legacy L-band All- sky Blind surveY's (WALLABY) Hi spectral line observations. The CNHI source finder also employs a novel sparse representation of 3-D objects, with a high compression rate, to implement Lutz one-pass algorithm on datacubes that are too large to process in a single pass. WALLABY will use ASKAP's phenomenal 30 square degree field of view to image \sim 70% of the sky. It is expected that WALLABY will find 500 000 Hi galaxies out to z \sim 0.2.

The process of determining the number and characteristics of sources in astronomical images is so fundamental to a large range of astronomical problems that it is perhaps surprising that no standard procedure has ever been defined that has well understood properties with a high degree of statistical rigour on completeness and reliability. There are now a large number of commonly used software tools for accomplishing this task, typically with different tools being used for images acquired using different technologies. Despite this, there have been relatively few quantitative analyses of the robustness or reliability of individual tools, or the details of the techniques they implement. We have an opportunity to redress this omission in the context of surveys planned with the Australian Square Kilometre Array Pathfinder (ASKAP). The Evolutionary Map of the Universe (EMU) survey with ASKAP, a continuum survey of the Southern Hemisphere up to declination +30 deg, aims to utilise an automated source identification and measurement approach that is demonstrably optimal, to maximise the reliability, utility and robustness of the resulting radio source catalogues. A key stage in source extraction methods is the background estimation (background level and noise level) and the choice of a threshold high enough to reject false sources yet not so high that the catalogues are significantly incomplete. In this analysis we present results from testing such algorithms as implemented in the SExtractor, Selavy (Duchamp), and sfind tools on simulated data. In particular the effects of background estimation, threshold and false-discovery rate settings are explored. For parameters that give similar completeness, the false-discovery rate method employed by sfind results in a more reliable catalogue compared to the peak threshold methods of SExtractor and Selavy.

Astronomical observations suggest that today's Universe is dominated by a dark energy of unknown physical origin. One of the most notable consequences in many models is that dark energy should cause the expansion of the Universe to accelerate: but the expansion rate as a function of time has proven very difficult to measure directly. We present a new determination of the cosmic expansion history by combining distant supernovae observations with a geometrical analysis of large-scale galaxy clustering within the WiggleZ Dark Energy Survey, using the Alcock-Paczynski test to measure the distortion of standard spheres. Our result constitutes a robust and non-parametric measurement of the Hubble expansion rate as a function of time, which we measure with 10-15% precision in four bins within the redshift range 0.1 < z < 0.9. We demonstrate that the cosmic expansion is accelerating, in a manner independent of the parameterization of the cosmological model (although assuming cosmic homogeneity in our data analysis). Furthermore, we find that this expansion history is consistent with a cosmological-constant dark energy.

We present measurements of the baryon acoustic peak at redshifts z = 0.44, 0.6 and 0.73 in the galaxy correlation function of the final dataset of the WiggleZ Dark Energy Survey. We combine our correlation function with lower-redshift measurements from the 6-degree Field Galaxy Survey and Sloan Digital Sky Survey, producing a stacked survey correlation function in which the statistical significance of the detection of the baryon acoustic peak is 4.9-sigma relative to a zero-baryon model with no peak. We fit cosmological models to this combined baryon acoustic oscillation (BAO) dataset comprising six distance-redshift data points, and compare the results to similar fits to the latest compilation of supernovae (SNe) and Cosmic Microwave Background (CMB) data. The BAO and SNe datasets produce consistent measurements of the equation-of-state w of dark energy, when separately combined with the CMB, providing a powerful check for systematic errors in either of these distance probes. Combining all datasets we determine w = -1.03 +/- 0.08 for a flat Universe, consistent with a cosmological constant model. Assuming dark energy is a cosmological constant and varying the spatial curvature, we find Omega_k = -0.004 +/- 0.006.

We report evidence of ordered orbital motion in luminous star-forming galaxies at z~1.3. We present integral field spectroscopy (IFS) observations, performed with the OH Suppressing InfraRed Imaging Spectrograph (OSIRIS) system, assisted by laser guide star adaptive optics on the Keck telescope, of 13 star-forming galaxies selected from the WiggleZ Dark Energy Survey. Selected via ultraviolet and [OII] emission, the large volume of the WiggleZ survey allows the selection of sources which have comparable intrinsic luminosity and stellar mass to IFS samples at z>2. Multiple 1-2 kpc size sub-components of emission, or 'clumps', are detected within the Halpha spatial emission which extends over 6-10 kpc in 4 galaxies, resolved compact emission (r<3 kpc) is detected in 5 galaxies, and extended regions of Halpha emission are observed in the remaining 4 galaxies. We discuss these data in the context of different snapshots in a merger sequence and/or the evolutionary stages of coalescence of star-forming regions in an unstable disk. We find evidence of ordered orbital motion in galaxies as expected from disk models and the highest values of velocity dispersion (\sigma>100 km/s) in the most compact sources. This unique data set reveals that the most luminous star-forming galaxies at z>1 are gaseous unstable disks indicating that a different mode of star formation could be feeding gas to galaxies at z>1, and lending support to theories of cold dense gas flows from the intergalactic medium.

We measure the imprint of baryon acoustic oscillations (BAOs) in the galaxy clustering pattern at the highest redshift achieved to date, z=0.6, using the distribution of N=132,509 emission-line galaxies in the WiggleZ Dark Energy Survey. We quantify BAOs using three statistics: the galaxy correlation function, power spectrum and the band-filtered estimator introduced by Xu et al. (2010). The results are mutually consistent, corresponding to a 4.0% measurement of the cosmic distance-redshift relation at z=0.6 (in terms of the acoustic parameter "A(z)" introduced by Eisenstein et al. (2005) we find A(z=0.6) = 0.452 +/- 0.018). Both BAOs and power spectrum shape information contribute toward these constraints. The statistical significance of the detection of the acoustic peak in the correlation function, relative to a wiggle-free model, is 3.2-sigma. The ratios of our distance measurements to those obtained using BAOs in the distribution of Luminous Red Galaxies at redshifts z=0.2 and z=0.35 are consistent with a flat Lambda Cold Dark Matter model that also provides a good fit to the pattern of observed fluctuations in the Cosmic Microwave Background (CMB) radiation. The addition of the current WiggleZ data results in a ~ 30% improvement in the measurement accuracy of a constant equation-of-state, w, using BAO data alone. Based solely on geometric BAO distance ratios, accelerating expansion (w < -1/3) is required with a probability of 99.8%, providing a consistency check of conclusions based on supernovae observations. Further improvements in cosmological constraints will result when the WiggleZ Survey dataset is complete.

We present precise measurements of the growth rate of cosmic structure for the redshift range 0.1 < z < 0.9, using redshift-space distortions in the galaxy power spectrum of the WiggleZ Dark Energy Survey. Our results, which have a precision of around 10% in four independent redshift bins, are well-fit by a flat LCDM cosmological model with matter density parameter Omega_m = 0.27. Our analysis hence indicates that this model provides a self-consistent description of the growth of cosmic structure through large-scale perturbations and the homogeneous cosmic expansion mapped by supernovae and baryon acoustic oscillations. We achieve robust results by systematically comparing our data with several different models of the quasi-linear growth of structure including empirical models, fitting formulae calibrated to N-body simulations, and perturbation theory techniques. We extract the first measurements of the power spectrum of the velocity divergence field, P_vv(k), as a function of redshift (under the assumption that P_gv(k) = -sqrt[P_gg(k) P_vv(k)] where g is the galaxy overdensity field), and demonstrate that the WiggleZ galaxy-mass cross-correlation is consistent with a deterministic (rather than stochastic) scale-independent bias model for WiggleZ galaxies for scales k < 0.3 h/Mpc. Measurements of the cosmic growth rate from the WiggleZ Survey and other current and future observations offer a powerful test of the physical nature of dark energy that is complementary to distance-redshift measures such as supernovae and baryon acoustic oscillations.

Correlations between the intrinsic shapes of galaxy pairs, and between the intrinsic shapes of galaxies and the large-scale density field, may be induced by tidal fields. These correlations, which have been detected at low redshifts (z<0.35) for bright red galaxies in the Sloan Digital Sky Survey (SDSS), and for which upper limits exist for blue galaxies at z~0.1, provide a window into galaxy formation and evolution, and are also an important contaminant for current and future weak lensing surveys. Measurements of these alignments at intermediate redshifts (z~0.6) that are more relevant for cosmic shear observations are very important for understanding the origin and redshift evolution of these alignments, and for minimising their impact on weak lensing measurements. We present the first such intermediate-redshift measurement for blue galaxies, using galaxy shape measurements from SDSS and spectroscopic redshifts from the WiggleZ Dark Energy Survey. Our null detection allows us to place upper limits on the contamination of weak lensing measurements by blue galaxy intrinsic alignments that, for the first time, do not require significant model-dependent extrapolation from the z~0.1 SDSS observations. Also, combining the SDSS and WiggleZ constraints gives us a long redshift baseline with which to constrain intrinsic alignment models and contamination of the cosmic shear power spectrum. Assuming that the alignments can be explained by linear alignment with the smoothed local density field, we find that a measurement of \sigma_8 in a blue-galaxy dominated, CFHTLS-like survey would be contaminated by at most +/-0.02 (95% confidence level, SDSS and WiggleZ) or +/-0.03 (WiggleZ alone) due to intrinsic alignments. [Abridged]

We report one of the most accurate measurements of the three-dimensional large-scale galaxy power spectrum achieved to date, using 56,159 redshifts of bright emission-line galaxies at effective redshift z=0.6 from the WiggleZ Dark Energy Survey at the Anglo-Australian Telescope. We describe in detail how we construct the survey selection function allowing for the varying target completeness and redshift completeness. We measure the total power with an accuracy of approximately 5% in wavenumber bands of dk=0.01 h/Mpc. A model power spectrum including non-linear corrections, combined with a linear galaxy bias factor and a simple model for redshift-space distortions, provides a good fit to our data for scales k < 0.4 h/Mpc. The large-scale shape of the power spectrum is consistent with the best-fitting matter and baryon densities determined by observations of the Cosmic Microwave Background radiation. By splitting the power spectrum measurement as a function of tangential and radial wavenumbers we delineate the characteristic imprint of peculiar velocities. We use these to determine the growth rate of structure as a function of redshift in the range 0.4 < z < 0.8, including a data point at z=0.78 with an accuracy of 20%. Our growth rate measurements are a close match to the self-consistent prediction of the LCDM model. The WiggleZ Survey data will allow a wide range of investigations into the cosmological model, cosmic expansion and growth history, topology of cosmic structure, and Gaussianity of the initial conditions. Our calculation of the survey selection function will be released at a future date via our website wigglez.swin.edu.au.

The WiggleZ Dark Energy Survey is a survey of 240,000 emission line galaxies in the distant universe, measured with the AAOmega spectrograph on the 3.9-m Anglo-Australian Telescope (AAT). The target galaxies are selected using ultraviolet photometry from the GALEX satellite, with a flux limit of NUV<22.8 mag. The redshift range containing 90% of the galaxies is 0.2<z<1.0. The primary aim of the survey is to precisely measure the scale of baryon acoustic oscillations (BAO) imprinted on the spatial distribution of these galaxies at look-back times of 4-8 Gyrs. Detailed forecasts indicate the survey will measure the BAO scale to better than 2% and the tangential and radial acoustic wave scales to approximately 3% and 5%, respectively. This paper provides a detailed description of the survey and its design, as well as the spectroscopic observations, data reduction, and redshift measurement techniques employed. It also presents an analysis of the properties of the target galaxies, including emission line diagnostics which show that they are mostly extreme starburst galaxies, and Hubble Space Telescope images, which show they contain a high fraction of interacting or distorted systems. In conjunction with this paper, we make a public data release of data for the first 100,000 galaxies measured for the project.

The WiggleZ Dark Energy Survey is a large-scale structure survey of intermediate-redshift UV-selected emission-line galaxies scheduled to cover 1000 sq deg, spanning a broad redshift range 0.2 < z < 1.0. The main scientific goal of the survey is the measurement of baryon acoustic oscillations (BAO) in the galaxy clustering pattern at a significantly higher redshift than previous studies. The BAO may be applied as a standard cosmological ruler to constrain dark energy models. Based on the first 20% of the dataset, we present initial results concerning the small-scale clustering of the WiggleZ targets, together with survey forecasts. The WiggleZ galaxy population possesses a clustering length r_0 = 4.40 +/- 0.12 Mpc/h, which is significantly larger than z=0 UV-selected samples, with a slope gamma = 1.92 +/- 0.08. This clustering length is comparable to z=3 Lyman Break Galaxies with similar UV luminosities. The full survey, scheduled for completion in 2010, will map an effective volume V_eff ~ 1 Gpc^3 (evaluated at a scale k = 0.15 h/Mpc) and will measure the angular-diameter distance and Hubble expansion rates in three redshift bins with accuracies ~ 5%. We will determine the value of a constant dark energy equation-of-state parameter, w, with a higher precision than existing supernovae observations using an entirely independent technique. The WiggleZ and supernovae measurements lie in highly complementary directions in the plane of w and the matter density Omega_m. The forecast using the full combination of WiggleZ, supernovae and CMB datasets is a marginalized error sigma(w) = 0.07, providing a robust and precise measurement of the properties of dark energy including cross-checking of systematic errors.

We present an investigation of quasar colour-redshift parameter space in order to search for radio-quiet red quasars and to test the ability of a variant of the KX quasar selection method to detect quasars over a full range of colour without bias. This is achieved by combining IRIS2 imaging with the complete Fornax Cluster Spectroscopic Survey to probe parameter space unavailable to other surveys. We construct a new sample of 69 quasars with measured bJ - K colours. We show that the colour distribution of these quasars is significantly different from that of the Large Bright Quasar Survey's quasars at a 99.9% confidence level. We find 11 of our sample of 69 quasars have signifcantly red colours (bJ - K >= 3.5) and from this, we estimate the red quasar fraction of the K <= 18.4 quasar population to be 31%, and robustly constrain it to be at least 22%. We show that the KX method variant used here is more effective than the UVX selection method, and has less colour bias than optical colour-colour selection methods.

We describe the `WiggleZ' spectroscopic survey of 280,000 star-forming galaxies selected from a combination of GALEX ultra-violet and SDSS + RCS2 optical imaging. The fundamental goal is a detection of the baryonic acoustic oscillations in galaxy clustering at high-redshift (0.5 < z < 1) and a precise measurement of the equation of state of dark energy from this purely geometric and robust method. The survey has already started on the 3.9m Anglo-Australian Telescope using the AAOmega spectrograph, and planned to complete during 2009. The WWW page for the survey can be found at http://astronomy.swin.edu.au/wigglez

We present the results of a search for ultra-compact dwarf galaxies (UCDs) in six different galaxy groups: Dorado, NGC1400, NGC0681, NGC4038, NGC4697 and NGC5084. We searched in the apparent magnitude range 17.5 < b_j < 20.5 (except NGC5084: 19.2 < b_j < 21.0). We found 1 definite plus 2 possible UCD candidates in the Dorado group and 2 possible UCD candidates in the NGC1400 group. No UCDs were found in the other groups. We compared these results with predicted luminosities of UCDs in the groups according to the hypothesis that UCDs are globular clusters formed in galaxies. The theoretical predictions broadly agree with the observational results, but deeper surveys are needed to fully test the predictions.

We have discovered nine ultra-compact dwarf galaxies (UCDs) in the Virgo Cluster, extending samples of these objects outside the Fornax Cluster. Using the 2dF multi-fiber spectrograph on the Anglo-Australian Telescope, the new Virgo members were found among 1500 color-selected, star-like targets with 16.0 < b_j < 20.2 in a two-degree diameter field centered on M87 (NGC4486). The newly-found UCDs are comparable to the UCDs in the Fornax Cluster, with sizes <~ 100 pc, -12.9 < M_B < -10.7, and exhibiting red, absorption-line spectra, indicative of an older stellar population. The properties of these objects remain consistent with the tidal threshing model for the origin of UCDs from the surviving nuclei of nucleated dwarf ellipticals disrupted in the cluster core, but can also be explained as objects that were formed by mergers of star clusters created in galaxy interactions. The discovery that UCDs exist in Virgo shows that this galaxy type is probably a ubiquitous phenomenon in clusters of galaxies; coupled with their possible origin by tidal threshing, the UCD population is a potential indicator and probe of the formation history of a given cluster. We also describe one additional bright UCD with M_B = -12.0 in the core of the Fornax Cluster. We find no further UCDs in our Fornax Cluster Spectroscopic Survey down to b_j = 19.5 in two additional 2dF fields extending as far as 3 degrees from the center of the cluster. All six Fornax bright UCDs identified with 2dF lie within 0.5 degree (projected distance of 170 kpc) of the central elliptical galaxy NGC1399.