
Recent interest in primordial black holes as a possible dark matter candidate
has motivated the reanalysis of previous methods for constraining massive
astrophysical compact objects in the Milky Way halo and beyond. In order to
derive these constraints, a model for the dark matter distribution around the
Milky Way must be used. Previous microlensing searches have assumed a
semiisothermal density sphere for this task. We show that this model is no
longer consistent with data from the Milky Way rotation curve, and recommend a
replacement. With a new model, we rederive some of the previous microlensing
constraints in the literature. We propagate some of the uncertainties
associated with the size and shape of the Milky Way halo. However, we find that
none of the models tested can sufficiently encapsulate the uncertainties in the
shape of the Milky Way halo. Our analysis reveals that the microlensing
constraints from the Large Magellanic Cloud weaken somewhat for MACHO masses
around $10\, M_\odot$ when this uncertainty is taken into account, but the
constraints tighten at lower masses. Exploring some of the simplifying
assumptions of previous constraints we also study the effect of wide mass
distributions of compact halo objects, as well as the effect of spatial
clustering on microlensing constraints. We find that both effects induce a
shift in the constraints towards smaller masses, and can effectively remove the
microlensing constraints from $M \sim 110 M_\odot$ for certain MACHO
populations.

Taipan is a multiobject spectroscopic galaxy survey starting in 2017 that
will cover 2pi steradians over the southern sky, and obtain optical spectra for
about two million galaxies out to z<0.4. Taipan will use the newlyrefurbished
1.2m UK Schmidt Telescope at Siding Spring Observatory with the new TAIPAN
instrument, which includes an innovative 'Starbugs' positioning system capable
of rapidly and simultaneously deploying up to 150 spectroscopic fibres (and up
to 300 with a proposed upgrade) over the 6deg diameter focal plane, and a
purposebuilt spectrograph operating from 370 to 870nm with resolving power
R>2000. The main scientific goals of Taipan are: (i) to measure the distance
scale of the Universe (primarily governed by the local expansion rate, H_0) to
1% precision, and the structure growth rate of structure to 5%; (ii) to make
the most extensive map yet constructed of the mass distribution and motions in
the local Universe, using peculiar velocities based on improved Fundamental
Plane distances, which will enable sensitive tests of gravitational physics;
and (iii) to deliver a legacy sample of lowredshift galaxies as a unique
laboratory for studying galaxy evolution as a function of mass and environment.
The final survey, which will be completed within 5 years, will consist of a
complete magnitudelimited sample (i<17) of about 1.2x10^6 galaxies,
supplemented by an extension to higher redshifts and fainter magnitudes
(i<18.1) of a luminous red galaxy sample of about 0.8x10^6 galaxies.
Observations and data processing will be carried out remotely and in a
fullyautomated way, using a purposebuilt automated 'virtual observer'
software and an automated data reduction pipeline. The Taipan survey is
deliberately designed to maximise its legacy value, by complementing and
enhancing current and planned surveys of the southern sky at wavelengths from
the optical to the radio.

Presently a ${>}3\sigma$ tension exists between values of the Hubble constant
$H_0$ derived from analysis of fluctuations in the Cosmic Microwave Background
by Planck, and local measurements of the expansion using calibrators of type Ia
supernovae (SNe Ia). We perform a blinded reanalysis of Riess et al. 2011 to
measure $H_0$ from lowredshift SNe Ia, calibrated by Cepheid variables and
geometric distances including to NGC 4258. This paper is a demonstration of
techniques to be applied to the Riess et at. 2016 data. Our endtoend analysis
starts from available CfA3 and LOSS photometry, providing an independent
validation of Riess et al. 2011. We obscure the value of $H_0$ throughout our
analysis and the first stage of the referee process, because calibration of SNe
Ia requires a series of often subtle choices, and the potential for results to
be affected by human bias is significant. Our analysis departs from that of
Riess et al. 2011 by incorporating the covariance matrix method adopted in SNLS
and JLA to quantify SN Ia systematics, and by including a simultaneous fit of
all SN Ia and Cepheid data. We find $H_0 = 72.5 \pm 3.1$ (stat) $\pm 0.77$
(sys) km s$^{1}$ Mpc$^{1}$ with a threegalaxy (NGC 4258+LMC+MW) anchor. The
relative uncertainties are 4.3% statistical, 1.1% systematic, and 4.4% total,
larger than in Riess et al. 2011 (3.3% total) and the Efstathiou 2014
reanalysis (3.4% total). Our error budget for $H_0$ is dominated by statistical
errors due to the small size of the supernova sample, whilst the systematic
contribution is dominated by variation in the Cepheid fits, and for the SNe Ia,
uncertainties in the host galaxy mass dependence and Malmquist bias.

Astronomers are often confronted with funky populations and distributions of
objects: brighter objects are more likely to be detected; targets are selected
based on colour cuts; imperfect classification yields impure samples. Failing
to account for these effects leads to biased analyses. In this paper we present
a simple overview of a Bayesian consideration of sample selection, giving
solutions to both analytically tractable and intractable models. This is
accomplished via a combination of analytic approximations and Monte Carlo
integration, in which dataset simulation is efficiently used to correct for
issues in the observed dataset. This methodology is also applicable for data
truncation, such as requiring densities to be strictly positive. Toy models are
included for demonstration, along with discussions of numerical considerations
and how to optimise for implementation. We provide sample code to demonstrate
the techniques. The methods in this paper should be widely applicable in fields
beyond astronomy, wherever sample selection effects occur.

We consider the dispersion on the supernova distanceredshift relation due to
peculiar velocities and gravitational lensing, and the sensitivity of these
effects to the amplitude of the matter power spectrum. We use the MeMo lensing
likelihood developed by Quartin, Marra & Amendola (2014), which accounts for
the characteristic nonGaussian distribution caused by lensing magnification
with measurements of the first four central moments of the distribution of
magnitudes. We build on the MeMo likelihood by including the effects of
peculiar velocities directly into the model for the moments. In order to
measure the moments from sparse numbers of supernovae, we take a new approach
using Kernel Density Estimation to estimate the underlying probability density
function of the magnitude residuals. We also describe a bootstrap resampling
approach to estimate the data covariance matrix. We then apply the method to
the Joint Lightcurve Analysis (JLA) supernova catalogue. When we impose only
that the intrinsic dispersion in magnitudes is independent of redshift, we find
$\sigma_8=0.44^{+0.63}_{0.44}$ at the one standard deviation level, although
we note that in tests on simulations, this model tends to overestimate the
magnitude of the intrinsic dispersion, and underestimate $\sigma_8$. We note
that the degeneracy between intrinsic dispersion and the effects of $\sigma_8$
is more pronounced when lensing and velocity effects are considered
simultaneously, due to a cancellation of redshift dependence when both effects
are included. Keeping the model of the intrinsic dispersion fixed as a Gaussian
distribution of width 0.14 mag, we find $\sigma_8 = 1.07^{+0.50}_{0.76}$.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV groundbased dark
energy experiment that will study baryon acoustic oscillations (BAO) and the
growth of structure through redshiftspace distortions with a widearea galaxy
and quasar redshift survey. To trace the underlying dark matter distribution,
spectroscopic targets will be selected in four classes from imaging data. We
will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to
even higher redshift, DESI will target bright [O II] emission line galaxies up
to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying
dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the
Ly$\alpha$ forest absorption features in their spectra, which will be used to
trace the distribution of neutral hydrogen. When moonlight prevents efficient
observations of the faint targets of the baseline survey, DESI will conduct a
magnitudelimited Bright Galaxy Survey comprising approximately 10 million
galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy
and quasar redshifts will be obtained to measure the BAO feature and determine
the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV groundbased dark
energy experiment that will study baryon acoustic oscillations and the growth
of structure through redshiftspace distortions with a widearea galaxy and
quasar redshift survey. The DESI instrument is a roboticallyactuated,
fiberfed spectrograph capable of taking up to 5,000 simultaneous spectra over
a wavelength range from 360 nm to 980 nm. The fibers feed ten threearm
spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500,
depending on wavelength. The DESI instrument will be used to conduct a
fiveyear survey designed to cover 14,000 deg$^2$. This powerful instrument
will be installed at prime focus on the 4m Mayall telescope in Kitt Peak,
Arizona, along with a new optical corrector, which will provide a threedegree
diameter field of view. The DESI collaboration will also deliver a
spectroscopic pipeline and data management system to reduce and archive all
data for eventual public use.

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 prereconstruction 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).

Type Ia supernovae are a powerful cosmological probe, that gave the first
strong evidence that the expansion of the universe is accelerating. Here we
provide an overview of how supernovae can go further to reveal information
about what is causing the acceleration, be it dark energy or some modification
to our laws of gravity. We first summarise the many different approaches used
to explain or test the acceleration, including parametric models (like the
standard model, LambdaCDM), nonparametric models, dark fluid models such as
quintessence, and extensions to standard gravity. We also show how supernova
data can be used beyond the Hubble diagram, to give information on
gravitational lensing and peculiar velocities that can be used to distinguish
between models that predict the same expansion history. Finally, we review the
methods of statistical inference that are commonly used, making a point of
separating parameter estimation from model selection.

In this paper we investigate how observational effects could possibly bias
cosmological inferences from peculiar velocity measurements. Specifically, we
look at how bulk flow measurements are compared with theoretical predictions.
Usually bulk flow calculations try to approximate the flow that would occur in
a sphere around the observer. Using the Horizon Run 2 simulation we show that
the traditional methods for bulk flow estimation can overestimate the magnitude
of the bulk flow for two reasons: when the survey geometry is not spherical
(the data do not cover the whole sky), and when the observations undersample
the velocity distributions. Our results may explain why several bulk flow
measurements found bulk flow velocities that seem larger than those expected in
standard {\Lambda}CDM cosmologies. We recommend a different approach when
comparing bulk flows to cosmological models, in which the theoretical
prediction for each bulk flow measurement is calculated specifically for the
geometry and sampling rate of that survey. This means that bulk flow values
will not be comparable between surveys, but instead they are comparable with
cosmological models, which is the more important measure.

The Australian Dark Energy Survey (OzDES) is a 100night spectroscopic survey
underway on the AngloAustralian Telescope using the fibrefed 2degreefield
(2dF) spectrograph. We have developed a new redshifting application Marz with
greater usability, flexibility, and the capacity to analyse a wider range of
object types than the Runz software package previously used for redshifting
spectra from 2dF. Marz is an opensource, clientbased, Javascript
webapplication which provides an intuitive interface and powerful automatic
matching capabilities on spectra generated from the AAOmega spectrograph to
produce high quality spectroscopic redshift measurements. The software can be
run interactively or via the command line, and is easily adaptable to other
instruments and pipelines if conforming to the current FITS file standard is
not possible. Behind the scenes, a modified version of the Autoz
crosscorrelation algorithm is used to match input spectra against a variety of
stellar and galaxy templates, and automatic matching performance for OzDES
spectra has increased from 54% (Runz) to 91% (Marz). Spectra not matched
correctly by the automatic algorithm can be easily redshifted manually by
cycling automatic results, manual template comparison, or marking spectral
features.

The ESSENCE survey discovered 213 Type Ia supernovae at redshifts 0.1 < z <
0.81 between 2002 and 2008. We present their R and Iband photometry, measured
from images obtained using the MOSAIC II camera at the CTIO 4 m Blanco
telescope, along with rapidresponse spectroscopy for each object. We use our
spectroscopic followup observations to determine an accurate, quantitative
classification and a precise redshift. Through an extensive calibration program
we have improved the precision of the CTIO Blanco natural photometric system.
We use several empirical metrics to measure our internal photometric
consistency and our absolute calibration of the survey. We assess the effect of
various potential sources of systematic bias on our measured fluxes, and we
estimate that the dominant term in the systematic error budget from the
photometric calibration on our absolute fluxes is ~1%.

We measure the bulk flow of the local Universe using the 6dF Galaxy Survey
peculiar velocity sample (6dFGSv), the largest and most homogeneous peculiar
velocity sample to date. 6dFGSv is a Fundamental Plane sample of $\sim10^4$
peculiar velocities covering the whole southern hemisphere for galactic
latitude $b > 10^\circ$, out to redshift ${z=0.0537}$. We apply the `Minimum
Variance' bulk flow weighting method, which allows us to make a robust
measurement of the bulk flow on scales of $50$ and $70\,h^{1}{\rm Mpc}$. We
investigate and correct for potential bias due to the lognormal velocity
uncertainties, and verify our method by constructing $\Lambda{\rm CDM}$ 6dFGSv
mock catalogues incorporating the survey selection function. For a hemisphere
of radius $50\,h^{1}{\rm Mpc}$ we find a bulk flow amplitude of
$U=248\pm58\,{\rm km}\,{\rm s}^{1}$ in the direction $(l,b) =
(318^\circ\pm20^\circ,40^\circ\pm13^\circ)$, and for $70\,h^{1}{\rm Mpc}$ we
find $U=243\pm58\,{\rm km}\,{\rm s}^{1}$, in the same direction. Our
measurement gives us a constraint on $\sigma_8$ of $1.01^{+1.07}_{0.58}$. Our
results are in agreement with other recent measurements of the direction of the
bulk flow, and our measured amplitude is consistent with a $\Lambda{\rm CDM}$
prediction.

As part of the OzDES spectroscopic survey we are carrying out a large scale
reverberation mapping study of $\sim$500 quasars over five years in the 30
deg$^2$ area of the Dark Energy Survey (DES) supernova fields. These quasars
have redshifts ranging up to 4 and have apparent AB magnitudes between
$16.8<r<22.5$ mag. The aim of the survey is to measure time lags between
fluctuations in the quasar continuum and broad emission line fluxes of
individual objects in order to measure black hole masses for a broad range of
AGN and constrain the radiusluminosity ($RL$) relationship. Here we
investigate the expected efficiency of the OzDES reverberation mapping campaign
and its possible extensions. We expect to recover lags for $\sim$3545\% of the
quasars. AGN with shorter lags and greater variability are more likely to yield
a lag, and objects with lags $\lesssim$6 months or $\sim$1 year are expected be
recovered the most accurately. The baseline OzDES reverberation mapping
campaign is predicted to produce an unbiased measurement of the $RL$
relationship parameters for H$\beta$, Mg II $\lambda$2798, and C IV
$\lambda$1549. However, extending the baseline survey by either increasing the
spectroscopic cadence, extending the survey season, or improving the emission
line flux measurement accuracy will significantly improve the $RL$ parameter
constraints for all broad emission lines.

Measurements of cosmological parameters via the distanceredshift relation
usually rely on models that assume a homogenous universe. It is commonly
presumed that the largescale structure evident in our Universe has a
negligible impact on the measurement if distances probed in observations are
sufficiently large (compared to the scale of inhomogeneities) and are averaged
over different directions on the sky. This presumption does not hold when
considering the effect of the gravitational redshift caused by our local
gravitational potential, which alters light coming from all distances and
directions in the same way. Despite its small magnitude, this local
gravitational redshift gives rise to noticeable effects in cosmological
inference using SN Ia data. Assuming conservative prior knowledge of the local
potential given by sampling a range of gravitational potentials at locations of
MilkyWaylike galaxies identified in cosmological simulations, we show that
ignoring the gravitational redshift effect in a standard data analysis leads to
an additional systematic error of ~1 per cent in the determination of density
parameters and the dark energy equation of state. We conclude that our local
gravitational field affects our cosmological inference at a level that is
important in future observations aiming to achieve percentlevel accuracy.

We present constraints on testing general relativity (GR) at cosmological
scales using recent data sets and assess the impact of galaxy intrinsic
alignment (IA) in the CFHTLenS lensing data on those constraints. We consider
CMB temperature data from Planck, the galaxy power spectrum from WiggleZ, weak
lensing tomography from the CFHTLenS, ISWgalaxy cross correlations, and BAO
data from 6dF, SDSS DR7, and BOSS DR9. We use a parameterization of the
modified gravity (MG) that is binned in redshift and scale, a parameterization
that evolves monotonically in scale but is binned in redshift, and a functional
parameterization that evolves only in redshift. We present the results in terms
of the MG parameters $Q$ and $\Sigma$. We employ an IA model with an amplitude
$A_{CFHTLenS}$ that is included in the parameter analysis. We find an
improvement in the constraints on the MG parameters corresponding to $4053\%$
increase on the figure of merit compared to previous studies, and GR is found
consistent with the data at the $95\%$ CL. The bounds found on $A_{CFHTLenS}$
are sensitive to the MG parameterization used, and the correlations between
$A_{CFHTLenS}$ and MG parameters are found to be weak to moderate. For all 3 MG
parameterizations $A_{\rm CFHTLenS}$ is found consistent with zero when the
whole lensing sample is used, however, when using the optimized earlytype
galaxy sample a significantly nonzero $A_{\rm CFHTLenS}$ is found for GR and
the third MG parameterization. We find that the tensions observed in previous
studies persist, and there is an indication that CMB data and lensing data
prefer different values for MG parameters, particularly for the parameter
$\Sigma$. The analysis of the confidence contours and probability distributions
suggest that the bimodality found follows that of the known tension in the
$\sigma_8$ parameter. (Abridged)

The Planck CMB experiment has delivered the best constraints so far on
primordial nonGaussianity, ruling out earlyUniverse models of inflation that
generate large nonGaussianity. Although small improvements in the CMB
constraints are expected, the next frontier of precision will come from future
largescale surveys of the galaxy distribution. The advantage of such surveys
is that they can measure many more modes than the CMB  in particular,
forthcoming radio surveys with the SKA will cover huge volumes. Radio continuum
surveys deliver the largest volumes, but with the disadvantage of no redshift
information. In order to mitigate this, we use two additional observables.
First, the integrated SachsWolfe effect  the crosscorrelation of the radio
number counts with the CMB temperature anisotropies  helps to reduce
systematics on the large scales that are sensitive to nonGaussianity. Second,
optical data allows for crossidentification in order to gain some redshift
information. We show that, while the single redshift bin case can provide a
sigma(fNL) ~ 20, and is therefore not competitive with current and future
constraints on nonGaussianity, a tomographic analysis could improve the
constraints by an order of magnitude, even with only two redshift bins. A huge
improvement is provided by the addition of highredshift sources, so having
crossID for highz galaxies and an even higherz radio tail is key to enabling
very precise measurements of fNL. Our results show that SKA continuum surveys
could provide constraints competitive with CMB and forthcoming optical surveys,
potentially allowing a measurement of sigma(fNL) ~ 1 to be made. Moreover,
these measurements would act as a useful check of results obtained with other
probes at other redshift ranges with other methods.

We present scaledependent measurements of the normalised growth rate of
structure $f\sigma_{8}(k, z=0)$ using only the peculiar motions of galaxies. We
use data from the 6degree Field Galaxy Survey velocity sample (6dFGSv)
together with a newlycompiled sample of lowredshift $(z < 0.07)$ type Ia
supernovae. We constrain the growth rate in a series of $\Delta k \sim 0.03
h{\rm Mpc^{1}}$ bins to $\sim35\%$ precision, including a measurement on
scales $>300 h^{1}{\rm Mpc}$, which represents one of the largestscale growth
rate measurement to date. We find no evidence for a scale dependence in the
growth rate, or any statistically significant variation from the growth rate as
predicted by the {\it Planck} cosmology. Bringing all the scales together, we
determine the normalised growth rate at $z=0$ to $\sim15\%$ in a manner {\it
independent} of galaxy bias and in excellent agreement with the constraint from
the measurements of redshiftspace distortions from 6dFGS. We pay particular
attention to systematic errors. We point out that the intrinsic scatter present
in FundamentalPlane and TullyFisher relations is only Gaussian in logarithmic
distance units; wrongly assuming it is Gaussian in linear (velocity) units can
bias cosmological constraints. We also analytically marginalise over zeropoint
errors in distance indicators, validate the accuracy of all our constraints
using numerical simulations, and demonstrate how to combine different
(correlated) velocity surveys using a matrix `hyperparameter' analysis.
Current and forthcoming peculiar velocity surveys will allow us to understand
in detail the growth of structure in the lowredshift universe, providing
strong constraints on the nature of dark energy.

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 shotnoise. 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
volumeaveragedistance, 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.

The way that peculiar velocities are often inferred from measurements of
distances and redshifts makes an approximation, v_p = czH_0 D, that gives
significant errors even at relatively low redshifts (overestimates peculiar
velocity by ~ 100 km/s at z~0.04). Here we demonstrate where the approximation
breaks down, the systematic offset it introduces, and how the exact calculation
should be implemented.

It has been only ~15 years since the discovery of dark energy (although some
may argue there were strong indications even earlier). In the short time since
measurements of type Ia supernovae indicated an accelerating universe, many
other techniques have now confirmed the acceleration is real. The variety of
ways in which dark energy has been confirmed is one of the reasons we are so
confident in the statement that most of the energy in the universe is in a form
we can not see except through its gravitational influence.
This review aims to summarise briefly the many varied ways we now have
measured dark energy. The fact that these different techniques all indicate
that the simplest model remains the best  that dark energy contributes a
constant background acceleration  is remarkable, since each of these
different types of measurements represented opportunities for this simplest
model to fail. Although we currently lack a compelling theoretical explanation
for this acceleration, any explanation will have to explain the wide variety of
complementary observations that we review here.
This is an informal presentation, following the lines of the talk I presented
at the General Relativity and Gravitation (GR20) conference in Warsaw in July
2013.

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 opticallyselected (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.

Higherorder statistics are a useful and complementary tool for measuring the
clustering of galaxies, containing information on the nongaussian evolution
and morphology of largescale structure in the Universe. In this work we
present measurements of the threepoint 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 Nbody simulations to predict the clustering of dark matter, we constrain
the linear and nonlinear 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 1020% accuracies, are
consistent with the predictions of the LCDM concordance cosmology and test this
model in a new way.

The growth history of largescale 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 2point
galaxy correlation function in parallel and transverse directions to the
lineofsight. 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 Nbody simulations for calibration. Our measurements
are consistent with an independent powerspectrum 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 nonparametric
manner, showing that it systematically increases with decreasing redshift, and
investigate the AlcockPaczynski 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 2dimensional
correlation function.

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 widefield infrared surveys such as the near
infrared UKIDSS LAS and the mid infrared WISE AllSky Survey. We find that ~30%
of the UV luminous WiggleZ galaxies are detected at >5sigma in the UKIDSSLAS
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
IRundetected, 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 bestfit 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 bestfit ages. The WiggleZ galaxies have star formation rates of
310 M0/yr and mostly lie at the upper end of the main sequence of starforming
galaxies at these redshifts. Their restframe UV luminosities and stellar
masses are comparable to both local compact UVluminous galaxies as well as
Lyman break galaxies at z~23.(abridged)