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

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.

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

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

Bayesian inference is the workhorse of gravitational-wave astronomy, for example, determining the mass and spins of merging black holes, revealing the neutron star equation of state, and unveiling the population properties of compact binaries. The science enabled by these inferences comes with a computational cost that can limit the questions we are able to answer. This cost is expected to grow. As detectors improve, the detection rate will go up, allowing less time to analyze each event. Improvement in low-frequency sensitivity will yield longer signals, increasing the number of computations per event. The growing number of entries in the transient catalog will drive up the cost of population studies. While Bayesian inference calculations are not entirely parallelizable, key components are embarrassingly parallel: calculating the gravitational waveform and evaluating the likelihood function. Graphical processor units (GPUs) are adept at such parallel calculations. We report on progress porting gravitational-wave inference calculations to GPUs. Using a single code - which takes advantage of GPU architecture if it is available - we compare computation times using modern GPUs (NVIDIA P100) and CPUs (Intel Gold 6140). We demonstrate speed-ups of $\sim 50 \times$ for compact binary coalescence gravitational waveform generation and likelihood evaluation and more than $100\times$ for population inference within the lifetime of current detectors. Further improvement is likely with continued development. Our python-based code is publicly available and can be used without familiarity with the parallel computing platform, CUDA.

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

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

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

We perform a combined analysis of cosmic shear tomography, galaxy-galaxy lensing tomography, and redshift-space multipole power spectra (monopole and quadrupole) using 450 deg$^2$ of imaging data by the Kilo Degree Survey (KiDS) overlapping with two spectroscopic surveys: the 2-degree Field Lensing Survey (2dFLenS) and the Baryon Oscillation Spectroscopic Survey (BOSS). We restrict the galaxy-galaxy lensing and multipole power spectrum measurements to the overlapping regions with KiDS, and self-consistently compute the full covariance between the different observables using a large suite of $N$-body simulations. We methodically analyze different combinations of the observables, finding that galaxy-galaxy lensing measurements are particularly useful in improving the constraint on the intrinsic alignment amplitude (by 30%, positive at $3.5\sigma$ in the fiducial data analysis), while the multipole power spectra are useful in tightening the constraints along the lensing degeneracy direction (e.g. factor of two stronger matter density constraint in the fiducial analysis). The fully combined constraint on $S_8 \equiv \sigma_8 \sqrt{\Omega_{\rm m}/0.3} = 0.742 \pm 0.035$, which is an improvement by 20% compared to KiDS alone, corresponds to a $2.6\sigma$ discordance with Planck, and is not significantly affected by fitting to a more conservative set of scales. Given the tightening of the parameter space, we are unable to resolve the discordance with an extended cosmology that is simultaneously favored in a model selection sense, including the sum of neutrino masses, curvature, evolving dark energy, and modified gravity. The complementarity of our observables allows for constraints on modified gravity degrees of freedom that are not simultaneously bounded with either probe alone, and up to a factor of three improvement in the $S_8$ constraint in the extended cosmology compared to KiDS alone.

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

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

We introduce gbpTrees: an algorithm for constructing merger trees from cosmological simulations, designed to identify and correct for pathological cases introduced by errors or ambiguities in the halo finding process. gbpTrees is built upon a halo matching method utilising pseudo-radial moments constructed from radially-sorted particle ID lists (no other information is required) and a scheme for classifying merger tree pathologies from networks of matches made to-and-from haloes across snapshots ranging forward-and-backward in time. Focusing on Subfind catalogs for this work, a sweep of parameters influencing our merger tree construction yields the optimal snapshot cadence and scanning range required for converged results. Pathologies proliferate when snapshots are spaced by $\lesssim{0.128}$ dynamical times; conveniently similar to that needed for convergence of semi-analytical modelling, as established by Benson etal Total merger counts are converged at the level of $\sim{5}$% for friends-of-friends (FoF) haloes of size $n_{\rm p}\gtrsim{75}$ across a factor of 512 in mass resolution, but substructure rates converge more slowly with mass resolution, reaching convergence of $\sim{10}$% for $n_{\rm p}\gtrsim{100}$ and particle mass $m_{\rm p}{\lesssim}10^{9}M_\odot$. We present analytic fits to FoF and substructure merger rates across nearly all observed galactic history ($z{\le}8.5$). While we find good agreement with the results presented by Fakhouri etal for FoF haloes, a slightly flatter dependance on merger ratio and increased major merger rates are found, reducing previously reported discrepancies with extended Press-Schechter estimates. When appropriately defined, substructure merger rates show a similar mass ratio dependance as FoF rates, but with stronger mass and redshift dependencies for their normalisation.

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

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

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

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

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

We present a new training set for estimating empirical photometric redshifts of galaxies, which was created as part of the 2dFLenS project. This training set is located in a 700 sq deg area of the KiDS South field and is randomly selected and nearly complete at r<19.5. We investigate the photometric redshift performance obtained with ugriz photometry from VST-ATLAS and W1/W2 from WISE, based on several empirical and template methods. The best redshift errors are obtained with kernel-density estimation, as are the lowest biases, which are consistent with zero within statistical noise. The 68th percentiles of the redshift scatter for magnitude-limited samples at r<(15.5, 17.5, 19.5) are (0.014, 0.017, 0.028). In this magnitude range, there are no known ambiguities in the colour-redshift map, consistent with a small rate of redshift outliers. In the fainter regime, the KDE method produces p(z) estimates per galaxy that represent unbiased and accurate redshift frequency expectations. The p(z) sum over any subsample is consistent with the true redshift frequency plus Poisson noise. Further improvements in redshift precision at r<20 would mostly be expected from filter sets with narrower passbands to increase the sensitivity of colours to small changes in redshift.

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 develop a statistical estimator to infer the redshift probability distribution of a photometric sample of galaxies from its angular cross-correlation in redshift bins with an overlapping spectroscopic sample. This estimator is a minimum variance weighted quadratic function of the data: a quadratic estimator. This extends and modifies the methodology presented by McQuinn & White (2013). The derived source redshift distribution is degenerate with the source galaxy bias, which must be constrained via additional assumptions. We apply this estimator to constrain source galaxy redshift distributions in the Kilo-Degree imaging survey through cross-correlation with the spectroscopic 2-degree Field Lensing Survey, presenting results first as a binned step-wise distribution in the range z < 0.8, and then building a continuous distribution using a Gaussian process model. We demonstrate the robustness of our methodology using mock catalogues constructed from N-body simulations, and comparisons with other techniques for inferring the redshift distribution.

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

Using a suite of cosmology simulations of a sample of $> 120$ galaxy clusters with $\log(M_{DM, vir}) \le 14.5$. We compare clusters that form in purely dark matter run and their counterparts in hydro runs and investigate 4 independent parameters, that are normally used to classify dynamical state. We find that the virial ratio $\eta$ in hydro-dynamical runs is $\sim 10$ per cent lower than in the DM run, and there is no clear separation between the relaxed and unrelaxed clusters for any parameter. Further, using the velocity dispersion deviation parameter $\zeta$, which is defined as the ratio between cluster velocity dispersion $\sigma$ and the theoretical prediction $\sigma_t = \sqrt{G M_{total}/R}$, we find that there is a linear correlation between the virial ratio $\eta$ and this $\zeta$ parameter. We propose to use this $\zeta$ parameter, which can be easily derived from observed galaxy clusters, as a substitute of the $\eta$ parameter to quantify the cluster dynamical state.

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

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

We present the 2-degree Field Lensing Survey (2dFLenS), a new galaxy redshift survey performed at the Anglo-Australian Telescope. 2dFLenS is the first wide-area spectroscopic survey specifically targeting the area mapped by deep-imaging gravitational lensing fields, in this case the Kilo-Degree Survey. 2dFLenS obtained 70,079 redshifts in the range z < 0.9 over an area of 731 sq deg, and is designed to extend the datasets available for testing gravitational physics and promote the development of relevant algorithms for joint imaging and spectroscopic analysis. The redshift sample consists first of 40,531 Luminous Red Galaxies (LRGs), which enable analyses of galaxy-galaxy lensing, redshift-space distortion, and the overlapping source redshift distribution by cross-correlation. An additional 28,269 redshifts form a magnitude-limited (r < 19.5) nearly-complete sub-sample, allowing direct source classification and photometric-redshift calibration. In this paper, we describe the motivation, target selection, spectroscopic observations, and clustering analysis of 2dFLenS. We use power spectrum multipole measurements to fit the redshift-space distortion parameter of the LRG sample in two redshift ranges 0.15 < z < 0.43 and 0.43 < z < 0.7 as beta = 0.49 +/- 0.15 and beta = 0.26 +/- 0.09, respectively. These values are consistent with those obtained from LRGs in the Baryon Oscillation Spectroscopic Survey. 2dFLenS data products will be released via our website http://2dflens.swin.edu.au.

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

We present the Dark-ages Reionization and Galaxy-formation Observables from Numerical Simulations (DRAGONS) program and Tiamat, the collisionless N-body simulation program upon which DRAGONS is built. The primary trait distinguishing Tiamat from other large simulation programs is its density of outputs at high redshift (100 from z=35 to z=5; roughly one every 10 Myr) enabling the construction of very accurate merger trees at an epoch when galaxy formation is rapid and mergers extremely frequent. We find that the friends-of-friends halo mass function agrees well with the prediction of Watson et al. at high masses, but deviates at low masses, perhaps due to our use of a different halo finder or perhaps indicating a break from "universal" behaviour. We then analyse the dynamical evolution of galaxies during the Epoch of Reionization finding that only a small fraction (~20%) of galactic halos are relaxed. We illustrate this using standard relaxation metrics to establish two dynamical recovery time-scales: i) halos need ~1.5 dynamical times following formation, and ii) ~2 dynamical times following a major (3:1) or minor (10:1) merger to be relaxed. This is remarkably consistent across a wide mass range. Lastly, we use a phase-space halo finder to illustrate that major mergers drive long-lived massive phase-space structures which take many dynamical times to dissipate. This can yield significant differences in the inferred mass build-up of galactic halos and we suggest that care must be taken to ensure a physically meaningful match between the galaxy-formation physics of semi-analytic models and the halo finders supplying their input.

This paper describes a new publicly available codebase for modelling galaxy formation in a cosmological context, the "Semi-Analytic Galaxy Evolution" model, or SAGE for short. SAGE is a significant update to that used in Croton et al. (2006) and has been rebuilt to be modular and customisable. The model will run on any N-body simulation whose trees are organised in a supported format and contain a minimum set of basic halo properties. In this work we present the baryonic prescriptions implemented in SAGE to describe the formation and evolution of galaxies, and their calibration for three N-body simulations: Millennium, Bolshoi, and GiggleZ. Updated physics include: gas accretion, ejection due to feedback, and reincorporation via the galactic fountain; a new gas cooling--radio mode active galactic nucleus (AGN) heating cycle; AGN feedback in the quasar mode; a new treatment of gas in satellite galaxies; and galaxy mergers, disruption, and the build-up of intra-cluster stars. Throughout, we show the results of a common default parameterization on each simulation, with a focus on the local galaxy population.

We introduce the Theoretical Astrophysical Observatory (TAO), an online virtual laboratory that houses mock observations of galaxy survey data. Such mocks have become an integral part of the modern analysis pipeline. However, building them requires an expert knowledge of galaxy modelling and simulation techniques, significant investment in software development, and access to high performance computing. These requirements make it difficult for a small research team or individual to quickly build a mock catalogue suited to their needs. To address this TAO offers access to multiple cosmological simulations and semi-analytic galaxy formation models from an intuitive and clean web interface. Results can be funnelled through science modules and sent to a dedicated supercomputer for further processing and manipulation. These modules include the ability to (1) construct custom observer light-cones from the simulation data cubes; (2) generate the stellar emission from star formation histories, apply dust extinction, and compute absolute and/or apparent magnitudes; and (3) produce mock images of the sky. All of TAO's features can be accessed without any programming requirements. The modular nature of TAO opens it up for further expansion in the future.

Galaxy clusters are an established and powerful test-bed for theories of both galaxy evolution and cosmology. Accurate interpretation of cluster observations often requires robust identification of the location of the centre. Using a statistical sample of clusters drawn from a suite of cosmological simulations in which we have explored a range of galaxy formation models, we investigate how the location of this centre is affected by the choice of observable - stars, hot gas, or the full mass distribution as can be probed by the gravitational potential. We explore several measures of cluster centre: the minimum of the gravitational potential, which would expect to define the centre if the cluster is in dynamical equilibrium; the peak of the density; the centre of BCG; and the peak and centroid of X-ray luminosity. We find that the centre of BCG correlates more strongly with the minimum of the gravitational potential than the X-ray defined centres, while AGN feedback acts to significantly enhance the offset between the peak X-ray luminosity and minimum gravitational potential. These results highlight the importance of centre identification when interpreting clusters observations, in particular when comparing theoretical predictions and observational data.

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

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.

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

We use the Gigaparsec WiggleZ (GiggleZ) simulations to characterise galaxy bias and its scale dependence for a range of redshifts and halo masses in a standard $\Lambda$LCDM cosmology. Assuming bias converges to a scale independent form at large scales, we develop a phenomenological model which fully expresses the mass and redshift dependence of bias and its scale dependence in real or redshift space. We then use this to illustrate how scale-dependent bias can systematically skew measurements of the growth-rate of cosmic structure obtained from redshift-space distortion measurements. When data is fit only to scales $k_{\rm max}{\le}0.1$ $[h^{-1} \rm{Mpc}]^{-1}$, we find that these effects are significant only for large biases ($b{\gtrsim}3$) at large redshifts ($z{\gtrsim}1$). However, when smaller scales are incorporated ($k_{\rm max}{\le}0.2$ $[h^{-1} \rm{Mpc}]^{-1}$) to increase measurement precision, the combination of reduced statistical uncertainty and increased scale dependent bias can result in highly significant systematics for most large halos across all redshifts. We identify several new interesting aspects of bias, including a significant large-scale bias boost for small halos at low-redshifts due to substructure effects ($\sim$20\% for Milky Way-like systems) and a nearly redshift-independent halo mass (corresponding to a redshift-space bias of ${\sim}1.5$) for which halo bias has little-or-no scale dependence on scales greater than $3$ $[h^{-1} {\rm Mpc}]$. This suggests an optimal strategy of targeting bias ${\sim}{1.5}$ systems for clustering studies which are dominated by systematic uncertainties in how observed halo (or galaxy) distributions map to their underlying mass distribution, such as cosmological measurements of neutrino masses. Code for generating our fitting formula is publicly available at http://gbpoole.github.io/Poole_2014a_code/ (Abridged)

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.

We have undertaken an ambitious program to visually classify all galaxies in the five CANDELS fields down to H<24.5 involving the dedicated efforts of 65 individual classifiers. Once completed, we expect to have detailed morphological classifications for over 50,000 galaxies up to z<4 over all the fields. Here, we present our detailed visual classification scheme, which was designed to cover a wide range of CANDELS science goals. This scheme includes the basic Hubble sequence types, but also includes a detailed look at mergers and interactions, the clumpiness of galaxies, $k$-corrections, and a variety of other structural properties. In this paper, we focus on the first field to be completed -- GOODS-S. The wide area coverage spanning the full field includes 7634 galaxies that have been classified by at least three different people. In the deep area of the field, 2534 galaxies have been classified by at least five different people at three different depths. With this paper, we release to the public all of the visual classifications in GOODS-S along with the GUI that we developed to classify galaxies. We find that the level of agreement among classifiers is good and depends on both the galaxy magnitude and the galaxy type, with disks showing the highest level of agreement and irregulars the lowest. A comparison of our classifications with the Sersic index and rest-frame colors shows a clear separation between disk and spheroid populations. Finally, we explore morphological k-corrections between the V-band and H-band observations and find that a small fraction (84 galaxies in total) are classified as being very different between these two bands. These galaxies typically have very clumpy and extended morphology or are very faint in the V-band.

We present the spectroscopic confirmation of a galaxy cluster at $z=2.095$ in the COSMOS field. This galaxy cluster was first reported in the ZFOURGE survey as harboring evolved massive galaxies using photometric redshifts derived with deep near-infrared (NIR) medium-band filters. We obtain medium resolution ($R \sim$ 3600) NIR spectroscopy with MOSFIRE on the Keck 1 telescope and secure 180 redshifts in a $12'\times12'$ region. We find a prominent spike of 57 galaxies at $z=2.095$ corresponding to the galaxy cluster. The cluster velocity dispersion is measured to be $\sigma_{\rm v1D}$ = 552 $\pm$ 52 km/s. This is the first study of a galaxy cluster in this redshift range ($z \gt 2.0$) with the combination of spectral resolution ($\sim$26 km/s) and the number of confirmed members (${>}50$) needed to impose a meaningful constraint on the cluster velocity dispersion and map its members over a large field of view. Our $\Lambda$CDM cosmological simulation suggests that this cluster will most likely evolve into a Virgo-like cluster with ${\rm M_{vir}}{=}10^{14.4\pm0.3} {\rm M_\odot}$ ($68\%$ confidence) at $z\sim$ 0. The theoretical expectation of finding such a cluster is $\sim$ $4\%$. Our results demonstrate the feasibility of studying galaxy clusters at $z > 2$ in the same detailed manner using multi-object NIR spectrographs as has been done in the optical in lower redshift clusters.

We study the relationship between the structure and star-formation rate (SFR) of X-ray selected low and moderate luminosity active galactic nuclei (AGNs) in the two Chandra Deep Fields, using Hubble Space Telescope imaging from the Cosmic Assembly Near Infrared Extragalactic Legacy Survey (CANDELS) and deep far-infrared maps from the PEP+GOODS-Herschel survey. We derive detailed distributions of structural parameters and FIR luminosities from carefully constructed control samples of galaxies, which we then compare to those of the AGNs. At z~1, AGNs show slightly diskier light profiles than massive inactive (non-AGN) galaxies, as well as modestly higher levels of gross galaxy disturbance (as measured by visual signatures of interactions and clumpy structure). In contrast, at z~2, AGNs show similar levels of galaxy disturbance as inactive galaxies, but display a red central light enhancement, which may arise due to a more pronounced bulge in AGN hosts or due to extinguished nuclear light. We undertake a number of tests of these alternatives, but our results do not strongly favour one interpretation over the other. The mean SFR and its distribution among AGNs and inactive galaxies are similar at z>1.5. At z<1, however, clear and significant enhancements are seen in the SFRs of AGNs with bulge-dominated light profiles. These trends suggest an evolution in the relation between nuclear activity and host properties with redshift, towards a minor role for mergers and interactions at z>1.5.

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.

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

We explore several correlations between various large-scale galaxy properties, particularly total globular cluster population (N_GCS), the central black hole mass (M_BH), velocity dispersion (nominally sigma_e), and bulge mass (M_dyn). Our data sample of 49 galaxies, for which both N_GC and M_BH are known, is larger than used in previous discussions of these two parameters and we employ the same sample to explore all pairs of correlations. Further, within this galaxy sample we investigate the scatter in each quantity, with emphasis on the range of published values for sigma_e and effective radius (R_e). We find that these two quantities in particular are difficult to measure consistently and caution that precise intercomparison of galaxy properties involving R_e and sigma_e is particularly difficult. Using both chi^2 and Monte Carlo Markov Chain (MCMC) fitting techniques, we show that quoted observational uncertainties for all parameters are too small to represent the true scatter in the data. We find that the correlation between M_dyn and N_GC is stronger than either the M_BH-sigma_e or M_BH-N_GC relations. We suggest that this is because both the galaxy bulge population ans N_GC were fundamentally established at an early epoch during the same series of star-forming events. By contrast, although the seed for M_BH was likely formed at a similar epoch, its growth over time is less similar from galaxy to galaxy and thus less predictable.

Peculiar velocity surveys, which measure galaxy peculiar velocities directly from standard candles in addition to redshifts, can provide strong constraints on the linear growth rate of cosmological large-scale structure at low redshift. The improvement originates from the physical relationship between galaxy density and peculiar velocity, which substantially reduces cosmic variance. We present the results of Fisher matrix forecasts of correlated fields of galaxy density and velocity. Peculiar velocity can improve the growth rate constraints by about a factor of two compared to density alone, if we can use all the information for wavenumber k < 0.2 h/Mpc. Future peculiar velocity surveys, TAIPAN, and the all-sky HI surveys, WALLABY and WNSHS, can measure the growth rate, f*sigma8 to 3 per cent at z ~ 0.025. Although the velocity subsample is about an order of magnitude smaller than the redshift sample from the same survey, it improves the constraint by 40 per cent compared to the same survey without velocity measurements. Peculiar velocity surveys can also measure the growth rate as a function of wavenumber with 15-30 per cent uncertainties in bins with widths 0.01 h/Mpc in the range 0.02 h/Mpc < k < 0.1 h/Mpc, which is a large improvement over galaxy density only. Such measurements on very large scales can detect signatures of modified gravity or non-Gaussianity through scale-dependent growth rate or galaxy bias. We use N-body simulations to improve the modelling of auto- and cross-power spectra of galaxy density and peculiar velocity by introducing a new redshift-space distortion term to the velocity, which has been neglected in previous studies. The velocity power spectrum has a damping in redshift space, which is larger than that naively expected from the similar effect in the galaxy power spectrum.

We present new and accurate measurements of the cosmic distance-redshift relation, spanning 0.2 < z < 1, using the topology of large-scale structure as a cosmological standard ruler. Our results derive from an analysis of the Minkowski functionals of the density field traced by the WiggleZ Dark Energy Survey. The Minkowski functionals are a set of statistics which completely describe the topological nature of each isodensity surface within the field, as a function of the density value. Given the shape of the underlying matter power spectrum, measured by fluctuations in the Cosmic Microwave Background radiation, the expected amplitudes of the Minkowski functionals are specified as an excursion set of a Gaussian random field, with minimal non-Gaussian corrections for the smoothing scales > 10 Mpc/h considered in this analysis. The measured amplitudes then determine the cosmic distance D_V(z), which we obtain with 3-7% accuracies in six independent redshift slices, with the standard ruler originating in the known curvature of the model power spectrum at the smoothing scale. We introduce a new method for correcting the topological statistics for the sparse-sampling of the density field by the galaxy tracers, and validate our overall methodology using mock catalogues from N-body simulations. Our distance measurements are consistent with standard models which describe the cosmic expansion history, and with previous analyses of baryon acoustic oscillations (BAOs) detected by the WiggleZ Survey, with the topological results yielding a higher distance precision by a factor of 2. However, the full redshift-space power-spectrum shape is required to recover the topological distances, in contrast to the preferred length scale imprinted by BAOs, which is determined by simpler physics.

(abridged) DYNAMO is a multi-wavelength, spatially-resolved survey of local ($z \sim 0.1$) star-forming galaxies designed to study evolution through comparison with samples at z~2. Half of the sample has integrated H-alpha luminosities of >$10^{42}$ erg/s, the typical lower limit for resolved spectroscopy at z~2. The sample covers a range in stellar mass ($10^9$-$10^{11}$ Msun) and star-formation rate (0.2-100 Msun/yr). In this first paper of a series, we present integral-field spectroscopy of H-alpha emission for the sample of 67 galaxies. We infer gas fractions in our sample as high as ~0.8, higher than typical for local galaxies. Gas fraction correlates with stellar mass in galaxies with star-formation rates below 10 Msun/yr, as found by COLDGASS, but galaxies with higher star-formation rates have higher than expected gas fractions. There is only a weak correlation, if any, between gas fraction and gas velocity dispersion. Galaxies in the sample visually classified as disc-like are offset from the local stellar-mass Tully-Fisher relation to higher circular velocities, but this offset vanishes when both gas and stars are included in the baryonic Tully-Fisher relation. The mean gas velocity dispersion of the sample is ~50 km/s, and V/sigma ranges from 2 to 10 for most of the discs, similar to 'turbulent' galaxies at high redshift. Half of our sample show disc-like rotation, while ~20 percent show no signs of rotation. The division between rotating and non-rotating is approximately...

We present the first application of a "multiple-tracer" redshift-space distortion (RSD) analysis to an observational galaxy sample, using data from the Galaxy and Mass Assembly survey (GAMA). Our dataset is an r < 19.8 magnitude-limited sample of 178,579 galaxies covering redshift interval z < 0.5 and area 180 deg^2. We obtain improvements of 10-20% in measurements of the gravitational growth rate compared to a single-tracer analysis, deriving from the correlated sample variance imprinted in the distributions of the overlapping galaxy populations. We present new expressions for the covariances between the auto-power and cross-power spectra of galaxy samples that are valid for a general survey selection function and weighting scheme. We find no evidence for a systematic dependence of the measured growth rate on the galaxy tracer used, justifying the RSD modelling assumptions, and validate our results using mock catalogues from N-body simulations. For multiple tracers selected by galaxy colour, we measure normalized growth rates in two independent redshift bins f*sigma_8(z=0.18) = 0.36 +/- 0.09 and f*sigma_8(z=0.38) = 0.44 +/- 0.06, in agreement with standard GR gravity and other galaxy surveys at similar redshifts.

We introduce a simple model to self-consistently connect the growth of galaxies to the formation history of their host dark matter haloes. Our model is defined by two simple functions: the "baryonic growth function" which controls the rate at which new baryonic material is made available for star formation, and the "physics function" which controls the efficiency with which this material is converted into stars. Using simple, phenomenologically motivated forms for both functions that depend only on a single halo property, we demonstrate the model's ability to reproduce the z=0 red and blue stellar mass functions. Furthermore, by adding redshift as a second input variable to the physics function we show that the reproduction of the global stellar mass function out to z=3 is improved. We conclude by discussing the general utility of our new model, highlighting its usefulness for creating mock galaxy samples which have a number of key advantages over those generated by other techniques.

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.

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.

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)

We use high-resolution N-body simulations to develop a new, flexible, empirical approach for measuring the growth rate from redshift-space distortions (RSD) in the 2-point galaxy correlation function. We quantify the systematic error in measuring the growth rate in a $1 \, h^{-3}$ Gpc$^3$ volume over a range of redshifts, from the dark matter particle distribution and a range of halo-mass catalogues with a number density comparable to the latest large-volume galaxy surveys such as the WiggleZ Dark Energy Survey and the Baryon Oscillation Spectroscopic Survey (BOSS). Our simulations allow us to span halo masses with bias factors ranging from unity (probed by emission-line galaxies) to more massive haloes hosting Luminous Red Galaxies. We show that the measured growth rate is sensitive to the model adopted for the small-scale real-space correlation function, and in particular that the "standard" assumption of a power-law correlation function can result in a significant systematic error in the growth rate determination. We introduce a new, empirical fitting function that produces results with a lower (5-10%) amplitude of systematic error. We also introduce a new technique which permits the galaxy pairwise velocity distribution, the quantity which drives the non-linear growth of structure, to be measured as a non-parametric stepwise function. Our (model-independent) results agree well with an exponential pairwise velocity distribution, expected from theoretical considerations, and are consistent with direct measurements of halo velocity differences from the parent catalogues. In a companion paper we present the application of our new methodology to the WiggleZ Survey dataset.

We investigate the ability of the Croton et al. (2006) semi-analytic model to reproduce the evolution of observed galaxies across the final 7 billion years of cosmic history. Using Monte-Carlo Markov Chain techniques we explore the available parameter space to produce a model which attempts to achieve a statistically accurate fit to the observed stellar mass function at z=0 and z~0.8, as well as the local black hole-bulge relation. We find that in order to be successful we are required to push supernova feedback efficiencies to extreme limits which are, in some cases, unjustified by current observations. This leads us to the conclusion that the current model may be incomplete. Using the posterior probability distributions provided by our fitting, as well as the qualitative details of our produced stellar mass functions, we suggest that any future model improvements must act to preferentially bolster star formation efficiency in the most massive halos at high redshift.

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.

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.

We present a detailed analysis of redshift-space distortions in the two-point correlation function of the 6dF Galaxy Survey (6dFGS). The K-band selected sub-sample which we employ in this study contains 81971 galaxies distributed over 17000deg^2 with an effective redshift z = 0.067. By modelling the 2D galaxy correlation function, xi(r_p,pi), we measure the parameter combination f(z)sigma_8(z) = 0.423 +/- 0.055. Alternatively, by assuming standard gravity we can break the degeneracy between sigma_8 and the galaxy bias parameter, b. Combining our data with the Hubble constant prior from Riess et al (2011), we measure sigma_8 = 0.76 +/- 0.11 and Omega_m = 0.250 +/- 0.022, consistent with constraints from other galaxy surveys and the Cosmic Microwave Background data from WMAP7. Combining our measurement of fsigma_8 with WMAP7 allows us to test the relationship between matter and gravity on cosmic scales by constraining the growth index of density fluctuations, gamma. Using only 6dFGS and WMAP7 data we find gamma = 0.547 +/- 0.088, consistent with the prediction of General Relativity. We note that because of the low effective redshift of 6dFGS our measurement of the growth rate is independent of the fiducial cosmological model (Alcock-Paczynski effect). We also show that our conclusions are not sensitive to the model adopted for non-linear redshift-space distortions. Using a Fisher matrix analysis we report predictions for constraints on fsigma_8 for the WALLABY survey and the proposed TAIPAN survey. The WALLABY survey will be able to measure fsigma_8 with a precision of 4-10%, depending on the modelling of non-linear structure formation. This is comparable to the predicted precision for the best redshift bins of the Baryon Oscillation Spectroscopic Survey (BOSS), demonstrating that low-redshift surveys have a significant role to play in future tests of dark energy and modified gravity.

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.

We report the first results from the Z-FOURGE survey: the discovery of a candidate galaxy cluster at z = 2.2 consisting of two compact overdensities with red galaxies detected at >20{\sigma} above the mean surface density. The discovery was made possible by a new deep (Ks < 24.8 AB 5{\sigma}) Magellan/FourStar near-IR imaging survey with 5 custom medium-bandwidth filters. The filters pinpoint the location of the Balmer/4000\AA\ break in evolved stellar populations at 1.5 < z < 3.5, yielding significantly more accurate photometric redshifts than possible with broadband imaging alone. The overdensities are within 1' of each other in the COSMOS field and appear to be embedded in a larger structure that contains at least one additional overdensity (\sim10{\sigma}). Considering the global properties of the overdensities, the z = 2.2 system appears to be the most distant example of a galaxy cluster with a population of red galaxies. A comparison to a large {\Lambda}CDM simulation suggests that the system may consist of merging subclusters, with properties in between those of z > 2 protoclusters with more diffuse distributions of blue galaxies and the lower-redshift galaxy clusters with prominent red sequences. The structure is completely absent in public optical catalogs in COSMOS and only weakly visible in a shallower near-IR survey. The discovery showcases the potential of deep near-IR surveys with medium-band filters to advance the understanding of environment and galaxy evolution at z > 1.5.

We present the properties of 8 star-forming regions, or 'clumps,' in 3 galaxies at z~1.3 from the WiggleZ Dark Energy Survey, which are resolved with the OSIRIS integral field spectrograph. Within turbulent discs, \sigma~90 km/s, clumps are measured with average sizes of 1.5 kpc and average Jeans masses of 4.2 x 10^9 \Msolar, in total accounting for 20-30 per cent of the stellar mass of the discs. These findings lend observational support to models that predict larger clumps will form as a result of higher disc velocity dispersions driven-up by cosmological gas accretion. As a consequence of the changes in global environment, it may be predicted that star-forming regions at high redshift should not resemble star-forming regions locally. Yet despite the increased sizes and dispersions, clumps and HII regions are found to follow tight scaling relations over the range z=0-2 for size, velocity dispersion, luminosity, and mass when comparing >2000 HII regions locally and 30 clumps at z>1 (\sigma \propto r^{0.42+/-0.03}, L(H\alpha) \propto r^{2.72+/-0.04}, L(H\alpha) \propto \sigma^{4.18+/-0.21}, and L(H\alpha) \propto M_{Jeans}^{1.24+/-0.05}). We discuss these results in the context of the existing simulations of clump formation and evolution, with an emphasis on the processes that drive-up the turbulent motions in the interstellar medium. Our results indicate that while the turbulence of discs may have important implications for the size and luminosity of regions which form within them, the same processes govern their formation from high redshift to the current epoch.

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

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

Upcoming next generation galactic surveys, such as GAIA and HERMES, will deliver unprecedented detail about the structure and make-up of our Galaxy, the Milky Way, and promise to radically improve our understanding of it. However, to benefit our broader knowledge of galaxy formation and evolution we first need to quantify how typical the Galaxy is with respect to other galaxies of its type. Through modeling and comparison with a large sample of galaxies drawn from the Sloan Digital Sky Survey and Galaxy Zoo, we provide tentative yet tantalizing evidence to show that both the Milky Way and nearby M31 are undergoing a critical transformation of their global properties. Both appear to possess attributes that are consistent with galaxies midway between the distinct blue and red bimodal color populations. In extragalactic surveys, such `green valley' galaxies are transition objects whose star formation typically will have all but extinguished in less than ~5 Gyrs. This finding reveals the possible future of our own galactic home, and opens a new window of opportunity to study such galactic transformations up close.

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.

High spatial and spectral resolution observations of star formation and kinematics in early galaxies have shown that two-thirds are massive rotating disk galaxies with the remainder being less massive non-rotating objects. The line of sight averaged velocity dispersions are typically five times higher than in today's disk galaxies. This has suggested that gravitationally-unstable, gas-rich disks in the early Universe are fuelled by cold, dense accreting gas flowing along cosmic filaments and penetrating hot galactic gas halos. However these accreting flows have not been observed, and cosmic accretion cannot power the observed level of turbulence. Here we report on a new sample of rare high-velocity-dispersion disk galaxies we have discovered in the nearby Universe where cold accretion is unlikely to drive their high star-formation rates. We find that the velocity dispersion is most fundamentally correlated with their star-formation rates, and not their mass nor gas fraction, which leads to a new picture where star formation itself is the energetic driver of galaxy disk turbulence at all cosmic epochs.

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 a combined X-ray and optical analysis of the cold front cluster Abell 1201 using archival Chandra, data and multi-object spectroscopy taken with the 3.9m Anglo Australian and 6.5m Multiple Mirror Telescopes. This paper represents the first in a series presenting a study of a sample of cold front clusters selected from the Chandra, archives with the aim of relating cold fronts to merger activity, understanding the dynamics of mergers and their effect on the cluster constituents. The Chandra X-ray imagery of Abell 1201 reveals two conspicuous surface brightness discontinuities, which are shown to be cold fronts, and a remnant core structure. Temperature maps reveal a complex multi-phase temperature structure with regions of hot gas interspersed with fingers of cold gas. Our optical analysis is based on a sample of 321 confirmed members, whose mean redshift is z=0.1673 +/- 0.0002 and velocity dispersion is 778 +/- 36 km/s. We search for dynamical substructure and find clear evidence for multiple localized velocity substructures coincident with over-densities in the galaxy surface density. Most notably, we find structure coincident with the remnant X-ray core. Despite the clear evidence for dynamical activity, we find the peculiar velocity distribution does not deviate significantly from Gaussian. We apply two-body dynamical analyses in order to assess which of the substructures are bound, and thus dynamically important in terms of the cluster merger history. We propose that the cold fronts in Abell 1201 are a consequence of its merger with a smaller subunit, which has induced gas motions that gave rise to `sloshing' cold fronts. Abell 1201 illustrates the value of combining multi-wavelength data and multiple substructure detection techniques when attempting to ascertain the dynamical state of a cluster.

(Abridged) We use the simulations presented in Poole et al. 2006 to examine the effects of mergers on compact cool cores in X-ray clusters. We propose a scheme for classifying the morphology of clusters based on their surface brightness and entropy profiles. Three dominant morphologies emerge: two hosting compact cores and central temperatures which are cool (CCC systems) or warm (CWC systems) and one hosting extended cores which are warm (EWC systems). We find that CCC states are disrupted only after direct collisions between cluster cores in head-on collisions or during second pericentric passage in off-axis mergers. By the time they relax, our remnant cores have generally been heated to warm core (CWC or EWC) states but subsequently recover CCC states. The only case resulting in a long-lived EWC state is a slightly off-axis 3:1 merger for which the majority of shock heating occurs during the accretion of a low-entropy stream formed from the disruption of the secondary's core. Compression prevents core temperatures from falling until after relaxation thus explaining the observed population of relaxed CWC systems with no need to invoke AGN feedback. The morphological segregation observed in the L_x-T_x and beta-r_c scaling relations is reflected in our simulations as well. However, none of the cases we have studied produce sufficiently high remnant central entropies to account for the most under-luminous EWC systems observed. Lastly, systems which initially host central metallicity gradients do not yield merger remnants with flat metallicity profiles. Taken together, these results suggest that once formed, compact core systems are remarkably stable against disruption from mergers. It remains to be demonstrated exactly how the sizable observed population of extended core systems was formed.

We use the suite of simulations presented in Poole et al (2006) to examine global X-ray and Sunyaev-Zel'dovich (SZ) observables for systems of merging relaxed X-ray clusters. In all of the scaling relations generated from the properties we have studied, we observe a generic evolution (in good qualitative agreement with previous authors): a rapid transient roughly along the mass scaling relations, a subsequent slow drift across the scatter until virialization, followed by a slow evolution along and up the mass scaling relations as cooling recovers in the cluster cores. However, this drift is not sufficient to account for the observed scatter in the scaling relations. We also study the effects of mergers on several theoretical temperature measures of the intracluster medium: emission weighted measures (T_ew), the spectroscopic-like measure proposed by Mazzotta et al (2004; T_sl) and plasma model fits to the integrated spectrum of the system (T_spec). We find that T_sl tracks T_spec for the entire duration of our mergers, illustrating that it remains a good tool for observational comparison even for highly disturbed systems. Furthermore, the transient temperature increases produced during first and second pericentric passage are 15-40% larger for T_ew than for T_sl or T_spec. This suggests that the effects of transient temperature increases on sigma_8 and Omega_M derived by Randall et al (2002) are over estimated. Lastly, we examine the X-ray SZ proxy proposed by Kravtsov et al (2006) and find that the tight mass scaling relation they predict remains secure through the entire duration of a merger event, independent of projection effects.

We present optical measurements of the faint end of the luminosity function in the core of the Coma cluster. Dwarf galaxies are detected down to a limiting magnitude of R approx. 25.75 in images taken with the Hubble Space Telescope. This represents the faintest determination of the Coma luminosity function to date. With the assumption that errors due to cosmic variance are small, evidence is found for a steep faint end slope with alpha < -2. Such a value is expected in theories in which reionization and other feedback processes are dependent on density.

We report on the analysis of a suite of SPH simulations (incorporating cooling and star formation) of mergers involving idealised X-ray clusters whose initial conditions resemble relaxed clusters with cool compact cores observed by Chandra and XMM. The simulations sample the most interesting, theoretically plausible, range of impact parameters and progenitor mass ratios. We find that all mergers evolve via a common progression. We illustrate this progression in the projected gas density, X-ray surface brightness, SZ, temperature, and gas entropy maps. Several different classes of transient ``cold front''-like features can arise over the course of a merger. We find that all of these classes are present in Chandra and XMM observations of merging systems and propose a naming scheme for these features: ``comet-like'' tails, bridges, plumes, streams and edges. In none of the cases considered do the initial cool compact cores of the primary and the secondary get destroyed during the course of the mergers. We quantify the evolving morphology of our mergers using centroid variance, power ratios and offset between the X-ray and the projected mass maps. We find that the centroid variance best captures the dynamical state of the cluster. Placing the system at z=0.1, we find that all easily identified observable traces of the secondary disappear from a simulated 50 ks Chandra image following the second pericentric passage. The system, however, takes approximately 2 additional Gyrs to relax and virialize. Temperature fluctuations at the level of 20% can persist in the final systems well past the point of virialization, suggesting that that the existence of temperature fluctuations, in and of themselves, do not necessarily indicate a disturbed or unrelaxed system.

(Abridged) Theoretical models that include only gravitationally-driven processes fail to match the observed mean X-ray properties of clusters. As a result, there has recently been increased interest in models in which either radiative cooling or entropy injection play a central role in mediating the properties of the intracluster medium. Both sets of models give reasonable fits to the mean properties of clusters, but cooling only models result in fractions of cold baryons in excess of observationally established limits and the simplest entropy injection models do not treat the "cooling core" structure present in many clusters and cannot account for entropy profiles revealed by recent X-ray observations. We consider models that marry radiative cooling with entropy injection, and confront model predictions for the global and structural properties of massive clusters with the latest X-ray data. The models successfully and simultaneously reproduce the observed L-T and L-M relations, yield detailed entropy, surface brightness, and temperature profiles in excellent agreement with observations, and predict a cooled gas fraction that is consistent with observational constraints. The model also provides a possible explanation for the significant intrinsic scatter present in the L-T and L-M relations and provides a natural way of distinguishing between clusters classically identified as "cooling flow" clusters and dynamically relaxed "non-cooling flow" clusters. The former correspond to systems that had only mild levels (< 300 keV cm^2) of entropy injection, while the latter are identified as systems that had much higher entropy injection. This is borne out by the entropy profiles derived from Chandra and XMM-Newton.

Recent X-ray observations have highlighted clusters that lack entropy cores. At first glance, these results appear to invalidate the preheated ICM models. We show that a self-consistent preheating model, which factors in the effects of radiative cooling, is in excellent agreement with the observations. Moreover, the model naturally explains the intrinsic scatter in the L-T relation, with ``cooling flow'' and ``non-cooling flow'' systems corresponding to mildly and strongly preheated systems, respectively. We discuss why preheating ought to be favoured over merging as a mechanism for the origin of ``non-cooling flow'' clusters.

Within the standard framework of structure formation, where clusters and groups of galaxies are built up from the merging of smaller systems, the physical properties of the intracluster medium, such as the gas temperature and the total X-ray luminosity, are predicted to possess well defined self-similar scaling relations. Observed clusters and groups, however, show strong deviations from these predicted relations. We argue that these deviations are unlikely to be entirely due to observational biasses; we assume they are physically based, due to the presence of excess entropy in the intracluster medium in addition to that generated by accretion shocks during the formation of the cluster. Several mechanisms have been suggested as a means of generating this entropy. Focussing on those mechanisms that preheat the gas before it becomes a constituent of the virialized cluster environment, we present a simple, intuitive, physically motivated, analytic model that successfully captures the important physics associated with the accretion of high entropy gas onto group and cluster-scale systems. We use the model to derive the new relationships between the observable properties of clusters and groups of galaxies, as well as the evolution of these relations. These include the luminosity-temperature and luminosity-$\sigma$ relations, as well as the temperature distribution function and X-ray luminosity function. These properties are found to be a more accurate description of the observations than those predicted from the standard framework. Future observations will that further test the efficacy of the preheated gas scenario are also discussed.

The Hubble Space Telescope has been used to obtain WFPC2 (V,I) photometry for a sample of 10,000 red giant stars in the outer halo of the giant elliptical NGC 5128. Comparison with the fiducial RGBs of Milky Way globular clusters and model isochrones demonstrates that this outer-halo population is completely dominated by old stars, with an extremely broad metallicity range extending from the most metal-poor Galactic globulars at [Fe/H] = -2 up to Solar abundance. By interpolation within the fiducial lines we derive the metallicity distribution function (MDF) for the halo, which turns out to have a distinct two-component structure: about 2/3 of the stars belong to a sharply peaked metal-richer component ([Fe/H] = -0.3, sigma = 0.22 dex), while the remaining one-third belong to a more extended metal-poor distribution sloping gradually down to low metallicity. The MDF can be remarkably well matched by a simple closed-box model of chemical enrichment, where the first component starts with an initial abundance Z_0 =0 and the second component with Z_0 = 0.18 Z_Sun. A basic "in situ" formation picture with two distinct epochs of star formation fits the observations better than other models involving major contributions from accretions or mergers.

The Hubble Space Telescope has been used to obtain WFPC2 (V,I) photometry for a large sample of stars in the outer halo of the giant elliptical NGC 5128 (d = 4 Mpc). The globular cluster N5128-C44, at the center of the Planetary Camera field, is well enough resolved to permit the construction of a color-magnitude diagram (CMD) for it which covers the brightest two magnitudes of the giant branch. The CMD is consistent with that of a normal old, moderately low-metallicity ([Fe/H] = -1.30 globular cluster, distinctly more metal-poor than most of the field halo stars at the same projected location (which average [Fe/H] ~ -0.5). This is the most distant globular cluster in which direct color-magnitude photometry has been achieved to date, and the first one belonging to a giant E galaxy.