We explore the cosmological multi-tracer synergies between a DESI-like Emission Line Galaxy (ELG) distribution and a Tianlai-like 21cm intensity map. We use simulated maps generated from a particle simulation in the light-cone volume (Horizon Run 4), sky-trimmed and including the effects of foreground contamination, its removal, and instrument noise. We first validate how the foreground residual affects the recovered 21cm signal by putting different levels of foreground contamination into the 21cm maps. We find that the contamination cannot be ignored in the angular auto-power spectra of HI even when it is small, but has no influence on the accuracy of the angular cross-power spectra between HI and galaxies. In the foreground-cleaned map case, as information is lost in the cleaning procedure, there is also a bias in the cross-power spectrum. However, we found that the bias from the cross-power spectrum is scale-independent, which is easily parameterised as part of the model, while the offset in the HI auto-power spectrum is non-linear. In particular, we tested that the cross power also benefits from the cancellation of the bias in the power spectrum measurement that is induced by the instrument noise, which changes the shape of the auto-power spectra but leaves the cross-power unaffected. We then modelled the angular cross-power spectra to fit the BAO feature in broadband shape of the angular cross power spectrum, including contamination from the residual foreground and the effect of instrument noise. We forecast a constraint on the angular diameter distance $D_\mathrm{A}$ for the Tianlai Pathfinder redshift $0.775<z<1.03$, giving a distance measurement with a precision of 3.3% at that redshift.

We present a description of the Australian Dark Energy Survey (OzDES) and summarise the results from its six years of operations. Using the 2dF fibre positioner and AAOmega spectrograph on the 3.9-metre Anglo-Australian Telescope, OzDES has monitored 771 AGN, classified hundreds of supernovae, and obtained redshifts for thousands of galaxies that hosted a transient within the 10 deep fields of the Dark Energy Survey. We also present the second OzDES data release, containing the redshifts of almost 30,000 sources, some as faint as $r_{\mathrm AB}=24$ mag, and 375,000 individual spectra. These data, in combination with the time-series photometry from the Dark Energy Survey, will be used to measure the expansion history of the Universe out to $z\sim1.2$ and the masses of hundreds of black holes out to $z\sim4$. OzDES is a template for future surveys that combine simultaneous monitoring of targets with wide-field imaging cameras and wide-field multi-object spectrographs.

The physics of gravity on cosmological scales affects both the rate of assembly of large-scale structure, and the gravitational lensing of background light through this cosmic web. By comparing the amplitude of these different observational signatures, we can construct tests that can distinguish General Relativity from its potential modifications. We use the latest weak gravitational lensing dataset from the Kilo-Degree Survey, KiDS-1000, in conjunction with overlapping galaxy spectroscopic redshift surveys BOSS and 2dFLenS, to perform the most precise existing amplitude-ratio test. We measure the associated E_G statistic with 15-20% errors, in five dz = 0.1 tomographic redshift bins in the range 0.2 < z < 0.7, on projected scales up to 100 Mpc/h. The scale-independence and redshift-dependence of these measurements are consistent with the theoretical expectation of General Relativity in a Universe with matter density Omega_m = 0.27 +/- 0.04. We demonstrate that our results are robust against different analysis choices, including schemes for correcting the effects of source photometric redshift errors, and compare the performance of angular and projected galaxy-galaxy lensing statistics.

We study the impact of different bias and redshift-space models on the halo power spectrum, quantifying their effect by comparing the fit to a subset of realizations taken from the WizCOLA suite. These provide simulated power spectrum measurements between $k_{\rm min}$ = 0.03 h/Mpc and $k_{\rm max}$ = 0.29 h/Mpc, constructed using the comoving Lagrangian acceleration method. For the bias prescription we include (i) simple linear bias; (ii) the McDonald & Roy model and (iii) its coevolution variant introduced by Saito et al.; and (iv) a very general model including all terms up to one-loop and corrections from advection. For the redshift-space modelling we include the Kaiser formula with exponential damping and the power spectrum provided by (i) tree-level perturbation theory and (ii) the Halofit prescription; (iii) one-loop perturbation theory, also with exponential damping; and (iv) an effective field theory description, also at one-loop, with damping represented by the EFT subtractions. We quantify the improvement from each layer of modelling by measuring the typical improvement in chi-square when fitting to a member of the simulation suite. We attempt to detect overfitting by testing for compatibility between the best-fit power spectrum per realization and the best-fit over the entire WizCOLA suite. For both bias and the redshift-space map we find that increasingly permissive models yield improvements in chi-square but with diminishing returns. The most permissive models show modest evidence for overfitting. Accounting for model complexity using the Bayesian Information Criterion, we argue that standard perturbation theory up to one-loop, or a related model such as that of Taruya, Nishimichi & Saito, coupled to the coevolution bias model, is likely to provide a good compromise for near-future galaxy surveys operating with comparable $k_{\rm max}$.

A violation of the distance-duality relation is directly linked with a temporal variation of the electromagnetic fine-structure constant. We consider a number of well-studied $f(T)$ gravity models and we revise the theoretical prediction of their corresponding induced violation of the distance-duality relationship. We further extract constraints on the involved model parameters through fine-structure constant variation data, alongside with supernovae data, and Hubble parameter measurements. Moreover, we constrain the evolution of the effective $f(T)$ gravitational constant. Finally, we compare with revised constraints on the phenomenological parametrisations of the violation of the equivalence principle in the electromagnetic sector.

The distribution of cosmological neutral hydrogen will provide a new window into the large-scale structure of the Universe with the next generation of radio telescopes and surveys. The observation of this material, through 21cm line emission, will be confused by foreground emission in the same frequencies. Even after these foregrounds are removed, the reconstructed map may not exactly match the original cosmological signal, which will introduce systematic errors and offset into the measured correlations. In this paper, we simulate future surveys of neutral hydrogen using the Horizon Run 4 (HR4) cosmological N-body simulation. We generate HI intensity maps from the HR4 halo catalogue, and combine with foreground radio emission maps from the Global Sky Model, to create accurate simulations over the entire sky. We simulate the HI sky for the frequency range 700-800 MHz, matching the sensitivity of the Tianlai pathfinder. We test the accuracy of the fastICA, PCA and log-polynomial fitting foreground removal methods to recover the input cosmological angular power spectrum and measure the parameters. We show the effect of survey noise levels and beam sizes on the recovered the cosmological constraints. We find that while the reconstruction removes power from the cosmological 21cm distribution on large-scales, we can correct for this and recover the input parameters in the noise-free case. However, the effect of noise and beam size of the Tianlai pathfinder prevents accurate recovery of the cosmological parameters when using only intensity mapping information.

We describe the observations and resultant galaxy cluster catalog from the 2770 deg$^2$ SPTpol Extended Cluster Survey (SPT-ECS). Clusters are identified via the Sunyaev-Zel'dovich (SZ) effect, and confirmed with a combination of archival and targeted follow-up data, making particular use of data from the Dark Energy Survey (DES). With incomplete followup we have confirmed as clusters 244 of 266 candidates at a detection significance $\xi \ge 5$ and an additional 204 systems at $4<\xi<5$. The confirmed sample has a median mass of $M_{500c} \sim {4.4 \times 10^{14} M_\odot h_{70}^{-1}}$, a median redshift of $z=0.49$, and we have identified 44 strong gravitational lenses in the sample thus far. Radio data are used to characterize contamination to the SZ signal; the median contamination for confirmed clusters is predicted to be $\sim$1% of the SZ signal at the $\xi>4$ threshold, and $<4\%$ of clusters have a predicted contamination $>10\% $ of their measured SZ flux. We associate SZ-selected clusters, from both SPT-ECS and the SPT-SZ survey, with clusters from the DES redMaPPer sample, and find an offset distribution between the SZ center and central galaxy in general agreement with previous work, though with a larger fraction of clusters with significant offsets. Adopting a fixed Planck-like cosmology, we measure the optical richness-to-SZ-mass ($\lambda-M$) relation and find it to be 28% shallower than that from a weak-lensing analysis of the DES data---a difference significant at the 4 $\sigma$ level---with the relations intersecting at $\lambda=60$ . The SPT-ECS cluster sample will be particularly useful for studying the evolution of massive clusters and, in combination with DES lensing observations and the SPT-SZ cluster sample, will be an important component of future cosmological analyses.

The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy and spacetime.

We present DASH (Deep Automated Supernova and Host classifier), a novel software package that automates the classification of the type, age, redshift, and host galaxy of supernova spectra. DASH makes use of a new approach that does not rely on iterative template matching techniques like all previous software, but instead classifies based on the learned features of each supernova's type and age. It has achieved this by employing a deep convolutional neural network to train a matching algorithm. This approach has enabled DASH to be orders of magnitude faster than previous tools, being able to accurately classify hundreds or thousands of objects within seconds. We have tested its performance on four years of data from the Australian Dark Energy Survey (OzDES). The deep learning models were developed using TensorFlow, and were trained using over 4000 supernova spectra taken from the CfA Supernova Program and the Berkeley SN Ia Program as used in SNID (Supernova Identification software, Blondin & Tonry 2007). Unlike template matching methods, the trained models are independent of the number of spectra in the training data, which allows for DASH's unprecedented speed. We have developed both a graphical interface for easy visual classification and analysis of supernovae, and a Python library for the autonomous and quick classification of several supernova spectra. The speed, accuracy, user-friendliness, and versatility of DASH presents an advancement to existing spectral classification tools. We have made the code publicly available on GitHub and PyPI (pip install astrodash) to allow for further contributions and development. The package documentation is available at https://astrodash.readthedocs.io.

Following on our purpose of developing a unified pipeline for large-scale structure data analysis with angular (i.e. harmonic-space) power spectra, we now include the weak lensing effect of magnification bias on galaxy clustering in a publicly available, modular parameter estimation code. We thus forecast constraints on the parameters of the concordance cosmological model, dark energy, and modified gravity theories from galaxy clustering tomographic angular power spectra. We find that a correct modelling of magnification is crucial in order not to bias the estimation of cosmological parameters, especially in the case of deep galaxy surveys. Our case study adopts specifications of the Evolutionary Map of the Universe (EMU), which is a full-sky, deep radio-continuum survey, and is expected to probe the Universe up to redshift $z\sim6$. We assume the Limber approximation, and include magnification bias on top of density fluctuations and redshift-space distortions. By restricting our analysis to the regime where the Limber approximation holds true, we significantly minimise the computational time needed, compared to that of the exact calculation. We also show that there is a trend for more biased parameter estimates from neglecting magnification when the redshift bins are very wide. We conclude that this result implies a strong dependence on the lensing contribution, which is an integrated effect and becomes dominant when wide redshift bins are considered. Finally, we note that instead of being considered a contaminant, magnification bias encodes important cosmological information, and its inclusion leads leading to an alleviation of its degeneracy between the galaxy bias and the amplitude normalisation of the matter fluctuations.

The volume of data that will be produced by the next generation of astrophysical instruments represents a significant opportunity for making unplanned and unexpected discoveries. Conversely, finding unexpected objects or phenomena within such large volumes of data presents a challenge that may best be solved using computational and statistical approaches. We present the application of a coarse-grained complexity measure for identifying interesting observations in large astronomical data sets. This measure, which has been termed apparent complexity, has been shown to model human intuition and perceptions of complexity. Apparent complexity is computationally efficient to derive and can be used to segment and identify interesting observations in very large data sets based on their morphological complexity. We show, using data from the Australia Telescope Large Area Survey, that apparent complexity can be combined with clustering methods to provide an automated process for distinguishing between images of galaxies which have been classified as having simple and complex morphologies. The approach generalizes well when applied to new data after being calibrated on a smaller data set, where it performs better than tested classification methods using pixel data. This generalizability positions apparent complexity as a suitable machine-learning feature for identifying complex observations with unanticipated features.

Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and analyze such data. On the positive side, it has become clear that the outstanding issues are purely technical in nature and can be solved with sufficient development activity. Such activity will enable science across redshifts, from early galaxy evolution in the pre-reionization era to dark energy evolution at low redshift.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

The Evolutionary Map of the Universe (EMU) is an all-sky survey in radio-continuum which uses the Australian SKA Pathfinder (ASKAP). Using galaxy angular power spectrum and the integrated Sachs-Wolfe effect, we study the potential of EMU to constrain models beyond $\Lambda$CDM (i.e., local primordial non-Gaussianity, dynamical dark energy, spatial curvature and deviations from general relativity), for different design sensitivities. We also include a multi-tracer analysis, distinguishing between star-forming galaxies and galaxies with an active galactic nucleus, to further improve EMU's potential. We find that EMU could measure the dark energy equation of state parameters around 35\% more precisely than existing constraints, and that the constraints on $f_{\rm NL}$ and modified gravity parameters will improve up to a factor $\sim2$ with respect to Planck and redshift space distortions measurements. With this work we demonstrate the promising potential of EMU to contribute to our understanding of the Universe.

We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.

We present a new measurement of $E_{\rm G}$, which combines measurements of weak gravitational lensing, real-space galaxy clustering and redshift space distortions. This statistic was proposed as a consistency test of General Relativity (GR) that is insensitive to linear, deterministic galaxy bias and the matter clustering amplitude. We combine deep imaging data from KiDS with overlapping spectroscopy from 2dFLenS, BOSS DR12 and GAMA and find $E_{\rm G}(\overline{z}=0.267)=0.43 \pm 0.13$ (GAMA), $E_{\rm G}(\overline{z}=0.305)=0.27 \pm 0.08$ (LOWZ+2dFLOZ) and $E_{\rm G}(\overline{z}=0.554)=0.26 \pm 0.07$ (CMASS+2dFHIZ). We demonstrate that the existing tension in the value of the matter density parameter hinders the robustness of this statistic as solely a test of GR. We find that our $E_{\rm G}$ measurements, as well as existing ones in the literature, favour a lower matter density cosmology than the Cosmic Microwave Background. For a flat $\Lambda$CDM Universe and assuming GR, we find $\Omega_{\rm m}(z=0)=0.25\pm0.03$. With this paper we publicly release the 2dFLenS dataset at: \url{http://2dflens.swin.edu.au}.

We present a weak gravitational lensing analysis of 815 square degree of $i$-band imaging from the Kilo-Degree Survey (KiDS-$i$-800). In contrast to the deep $r$-band observations, which take priority during excellent seeing conditions and form the primary KiDS dataset (KiDS-$r$-450), the complementary yet shallower KiDS-$i$-800 spans a wide range of observing conditions. The overlapping KiDS-$i$-800 and KiDS-$r$-450 imaging therefore provides a unique opportunity to assess the robustness of weak lensing measurements. In our analysis, we introduce two new `null' tests. The `nulled' two-point shear correlation function uses a matched catalogue to show that the calibrated KiDS-$i$-800 and KiDS-$r$-450 shear measurements agree at the level of $1 \pm 4$\%. We use five galaxy lens samples to determine a `nulled' galaxy-galaxy lensing signal from the full KiDS-$i$-800 and KiDS-$r$-450 surveys and find that the measurements agree to $7 \pm 5$\% when the KiDS-$i$-800 source redshift distribution is calibrated using either spectroscopic redshifts, or the 30-band photometric redshifts from the COSMOS survey.

We present a machine-learning photometric redshift analysis of the Kilo-Degree Survey Data Release 3, using two neural-network based techniques: ANNz2 and MLPQNA. Despite limited coverage of spectroscopic training sets, these ML codes provide photo-zs of quality comparable to, if not better than, those from the BPZ code, at least up to zphot<0.9 and r<23.5. At the bright end of r<20, where very complete spectroscopic data overlapping with KiDS are available, the performance of the ML photo-zs clearly surpasses that of BPZ, currently the primary photo-z method for KiDS. Using the Galaxy And Mass Assembly (GAMA) spectroscopic survey as calibration, we furthermore study how photo-zs improve for bright sources when photometric parameters additional to magnitudes are included in the photo-z derivation, as well as when VIKING and WISE infrared bands are added. While the fiducial four-band ugri setup gives a photo-z bias $\delta z=-2e-4$ and scatter $\sigma_z<0.022$ at mean z = 0.23, combining magnitudes, colours, and galaxy sizes reduces the scatter by ~7% and the bias by an order of magnitude. Once the ugri and IR magnitudes are joined into 12-band photometry spanning up to 12 $\mu$, the scatter decreases by more than 10% over the fiducial case. Finally, using the 12 bands together with optical colours and linear sizes gives $\delta z<4e-5$ and $\sigma_z<0.019$. This paper also serves as a reference for two public photo-z catalogues accompanying KiDS DR3, both obtained using the ANNz2 code. The first one, of general purpose, includes all the 39 million KiDS sources with four-band ugri measurements in DR3. The second dataset, optimized for low-redshift studies such as galaxy-galaxy lensing, is limited to r<20, and provides photo-zs of much better quality than in the full-depth case thanks to incorporating optical magnitudes, colours, and sizes in the GAMA-calibrated photo-z derivation.

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.

We present results for the first three years of OzDES, a six-year programme to obtain redshifts for objects in the Dark Energy Survey (DES) supernova fields using the 2dF fibre positioner and AAOmega spectrograph on the Anglo-Australian Telescope. OzDES is a multi-object spectroscopic survey targeting multiple types of targets at multiple epochs over a multi-year baseline, and is one of the first multi-object spectroscopic surveys to dynamically include transients into the target list soon after their discovery. At the end of three years, OzDES has spectroscopically confirmed almost 100 supernovae, and has measured redshifts for 17,000 objects, including the redshifts of 2,566 supernova hosts. We examine how our ability to measure redshifts for targets of various types depends on signal-to-noise, magnitude, and exposure time, finding that our redshift success rate increases significantly at a signal-to-noise of 2 to 3 per 1-Angstrom bin. We also find that the change in signal-to-noise with exposure time closely matches the Poisson limit for stacked exposures as long as 10 hours. We use these results to predict the redshift yield of the full OzDES survey, as well as the potential yields of future surveys on other facilities such as the 4m Multi-Object Spectroscopic Telescope (4MOST), the Subaru Prime Focus Spectrograph (PFS), and the Maunakea Spectroscopic Explorer (MSE). This work marks the first OzDES data release, comprising 14,693 redshifts. OzDES is on target to obtain over a yield of approximately 5,700 supernova host-galaxy redshifts.

In 'modified' gravity the observed acceleration of the universe is explained by changing the gravitational force law or the number of degrees of freedom in the gravitational sector. Both possibilities can be tested by measurements of cosmological structure formation. In this paper we elaborate the details of such tests using the Galileon model as a case study. We pay attention to the possibility that each new degree of freedom may have stochastically independent initial conditions, generating different types of potential well in the early universe and breaking complete correlation between density and velocity power spectra. This 'stochastic bias' can confuse schemes to parametrize the predictions of modified gravity models, such as the use of the growth parameter f alone. Using data from the WiggleZ Dark Energy Survey we show that it will be possible to obtain constraints using information about the cosmological-scale force law embedded in the multipole power spectra of redshift-space distortions. As an example, we obtain an upper limit on the strength of the conformal coupling to matter in the cubic Galileon model, giving |1/M| < 200 / Mp. This allows the fifth-force to be stronger than gravity, but is consistent with zero coupling.

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 the impact of astrophysical systematics on cosmic shear cosmological parameter constraints from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), and the concordance with cosmic microwave background measurements by Planck. We present updated CFHTLenS cosmic shear tomography measurements extended to degree scales using a covariance calibrated by a new suite of N-body simulations. We analyze these measurements with a new model fitting pipeline, accounting for key systematic uncertainties arising from intrinsic galaxy alignments, baryonic effects in the nonlinear matter power spectrum, and photometric redshift uncertainties. We examine the impact of the systematic degrees of freedom on the cosmological parameter constraints, both independently and jointly. When the systematic uncertainties are considered independently, the intrinsic alignment amplitude is the only degree of freedom that is substantially preferred by the data. When the systematic uncertainties are considered jointly, there is no consistently strong preference in favor of the more complex models. We quantify the level of concordance between the CFHTLenS and Planck datasets by employing two distinct data concordance tests, grounded in Bayesian evidence and information theory. We find that the two data concordance tests largely agree with one another, and that the level of concordance between the CFHTLenS and Planck datasets is sensitive to the exact details of the systematic uncertainties included in our analysis, ranging from decisive discordance to substantial concordance as the treatment of the systematic uncertainties becomes more conservative. The least conservative scenario is the one most favored by the cosmic shear data, but it is also the one that shows the greatest degree of discordance with Planck. The data and analysis code are public at https://github.com/sjoudaki/cfhtlens_revisited

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

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.

We explore the distance-redshift relation using a cosmographic methodology, and show how the cosmographic parameters can be used to determine the redshift of transition from deceleration to acceleration. Such a transition at a low redshift occupies only a small region of the available parameter space, and the prior assumption of an early period of deceleration can significantly change the posterior constraints. We use available type Ia Supernovae (SN-Ia) and Baryon Acoustic Oscillation (BAO) data sets to determine the cosmographic deceleration $q_0$, jerk $j_0$, snap $s_0$ and lerk $l_0$ parameters. The parameters are consistent with the $\Lambda$CDM model for a flat universe within 2-sigma. We derive constraints on the redshift of transition from deceleration to acceleration for the different expansions, and find $z_{\rm acc} > 0.14$ at 95% confidence in the most conservative case.

As an alternative explanation for the cosmic acceleration, $f(R)$ theories of gravity can predict an almost identical expansion history to standard $\Lambda$CDM, yet make very different predictions for the growth of cosmological structures. Measurements of the cosmic bulk flow provides a method for determining the strength of gravity over the history of structure formation. We use the modified gravity N-body code ECOSMOG to simulate dark matter particles and make predictions for the bulk flow magnitude in both $\Lambda$CDM and $f(R)$ gravity. With the peculiar velocities output by ECOSMOG we determine the bulk flow at depths ranging from $20h^{-1}$Mpc to $50h^{-1}$Mpc, following the redshift and sky distribution of the 2MASS Tully-Fisher survey (2MTF). At each depth, we find that the $\Lambda$CDM and $f_{R0} = 10^{-5}$ simulations produce bulk flow measurements that are consistent with $\Lambda$CDM predictions and the 2MTF survey at a $1\sigma$ level. We also find that adopting an $f(R)$ strength of $f_{R0} = 10^{-3}$ predict a much larger value for the bulk flow, which disagree with $\Lambda$CDM predictions at all depths considered. We conclude that $f_{R0}$ must be constrained to a level no greater than $10^{-4}$ to agree with bulk flow measurements.

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

In this work we model the quintessence potential in a Taylor series expansion, up to second order, around the present-day value of the scalar field. The field is evolved in a thawing regime assuming zero initial velocity. We use the latest data from the Planck satellite, baryonic acoustic oscillations observations from the Sloan Digital Sky Survey, and Supernovae luminosity distance information from Union2.1 to constrain our models parameters, and also include perturbation growth data from the WiggleZ, BOSS and the 6dF surveys. The supernova data provide the strongest individual constraint on the potential parameters. We show that the growth data performance is competitive with the other datasets in constraining the dark energy parameters we introduce. We also conclude that the combined constraints we obtain for our model parameters, when compared to previous works of nearly a decade ago, have shown only modest improvement, even with new growth of structure data added to previously-existent types of data.

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

Dark energy may be the first sign of new fundamental physics in the Universe, taking either a physical form or revealing a correction to Einsteinian gravity. Weak gravitational lensing and galaxy peculiar velocities provide complementary probes of General Relativity, and in combination allow us to test modified theories of gravity in a unique way. We perform such an analysis by combining measurements of cosmic shear tomography from the Canada-France Hawaii Telescope Lensing Survey (CFHTLenS) with the growth of structure from the WiggleZ Dark Energy Survey and the Six-degree-Field Galaxy Survey (6dFGS), producing the strongest existing joint constraints on the metric potentials that describe general theories of gravity. For scale-independent modifications to the metric potentials which evolve linearly with the effective dark energy density, we find present-day cosmological deviations in the Newtonian potential and curvature potential from the prediction of General Relativity to be (Delta Psi)/Psi = 0.05 \pm 0.25 and (Delta Phi)/Phi = -0.05 \pm 0.3 respectively (68 per cent CL).

We investigate the prior dependence on the inferred spectrum of primordial tensor perturbations, in light of recent results from BICEP2 and taking into account a possible dust contribution to polarized anisotropies. We highlight an optimized parameterization of the tensor power spectrum, and adoption of a logarithmic prior on its amplitude $A_T$, leading to results that transform more evenly under change of pivot scale. In the absence of foregrounds the tension between the results of BICEP2 and Planck drives the tensor spectral index $n_T$ to be blue-tilted in a joint analysis, which would be in contradiction to the standard inflation prediction ($n_T<0$). When foregrounds are accounted for, the BICEP2 results no longer require non-standard inflationary parameter regions. We present limits on primordial $A_T$ and $n_T$, adopting foreground scenarios put forward by Mortonson & Seljak and motivated by Planck 353 GHz observations, and assess what dust contribution leaves a detectable cosmological signal. We find that if there is sufficient dust for the signal to be compatible with standard inflation, then the primordial signal is too weak to be robustly detected by BICEP2 if Planck+WMAP upper limits from temperature and E-mode polarization are correct.

We investigate how well future large-scale radio surveys could measure different shapes of primordial non-Gaussianity; in particular we focus on angle-dependent non-Gaussianity arising from primordial anisotropic sources, whose bispectrum has an angle dependence between the three wavevectors that is characterized by Legendre polynomials $\mathcal{P}_L$ and expansion coefficients $c_L$. We provide forecasts for measurements of galaxy power spectrum, finding that Large-Scale Structure (LSS) data could allow measurements of primordial non-Gaussianity competitive or improving upon current constraints set by CMB experiments, for all the shapes considered. We argue that the best constraints will come from the possibility to assign redshift information to radio galaxy surveys, and investigate a few possible scenarios for the EMU and SKA surveys. A realistic (futuristic) modeling could provide constraints of $f_{\rm NL}^{\rm loc} \approx 1 (0.5)$ for the local shape, $f_{\rm NL}$ of $\mathcal{O}(10) (\mathcal{O}(1))$ for the orthogonal, equilateral and folded shapes, and $c_{L=1} \approx 80 (2)$, $c_{L=2} \approx 400 (10)$ for angle-dependent non-Gaussianity. The more futuristic forecasts show the potential of LSS analyses to considerably improve current constraints on non-Gaussianity, and so on models of the primordial Universe. Finally, we find the minimum requirements that would be needed to reach $\sigma(c_{L=1})=10$, which can be considered as a typical (lower) value predicted by some (inflationary) models.

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

We present measurements of both scale- and time-dependent deviations from the standard gravitational field equations. These late-time modifications are introduced separately for relativistic and non-relativistic particles, by way of the parameters $G_{\rm matter}(k,z)$ and $G_{\rm light}(k,z)$ using two bins in both scale and time, with transition wavenumber $0.01$ Mpc$^{-1}$ and redshift 1. We emphasize the use of two dynamical probes to constrain this set of parameters, galaxy power spectrum multipoles and the direct peculiar velocity power spectrum, which probe fluctuations on different scales. The multipole measurements are derived from the WiggleZ and BOSS Data Release 11 CMASS galaxy redshift surveys and the velocity power spectrum is measured from the velocity sub-sample of the 6-degree Field Galaxy Survey. We combine with additional cosmological probes including baryon acoustic oscillations, Type Ia SNe, the cosmic microwave background (CMB), lensing of the CMB, and the temperature--galaxy cross-correlation. Using a Markov Chain Monte Carlo likelihood analysis, we find the inferred best-fit parameter values of $G_{\rm matter}(k,z)$ and $G_{\rm light}(k,z)$ to be consistent with the standard model at the $95\%$ confidence level. Furthermore, accounting for the Alcock-Paczynski effect, we perform joint fits for the expansion history and growth index gamma; we measure $\gamma = 0.665 \pm 0.0669$ ($68\%$ C.L) for a fixed expansion history, and $\gamma = 0.73^{+0.08}_{-0.10}$ ($68\%$ C.L) when the expansion history is allowed to deviate from $\Lambda$CDM. With a fixed expansion history the inferred value is consistent with GR at the $95\%$ C.L; alternatively, a $2\sigma$ tension is observed when the expansion history is not fixed, this tension is worsened by the combination of growth and SNe data.

OzDES is a five-year, 100-night, spectroscopic survey on the Anglo-Australian Telescope, whose primary aim is to measure redshifts of approximately 2,500 Type Ia supernovae host galaxies over the redshift range 0.1 < z < 1.2, and derive reverberation-mapped black hole masses for approximately 500 active galactic nuclei and quasars over 0.3 < z < 4.5. This treasure trove of data forms a major part of the spectroscopic follow-up for the Dark Energy Survey for which we are also targeting cluster galaxies, radio galaxies, strong lenses, and unidentified transients, as well as measuring luminous red galaxies and emission line galaxies to help calibrate photometric redshifts. Here we present an overview of the OzDES program and our first-year results. Between Dec 2012 and Dec 2013, we observed over 10,000 objects and measured more than 6,000 redshifts. Our strategy of retargeting faint objects across many observing runs has allowed us to measure redshifts for galaxies as faint as m_r=25 mag. We outline our target selection and observing strategy, quantify the redshift success rate for different types of targets, and discuss the implications for our main science goals. Finally, we highlight a few interesting objects as examples of the fortuitous yet not totally unexpected discoveries that can come from such a large spectroscopic survey.

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

In order to explain cosmic acceleration without invoking "dark" physics, we consider $f(R)$ modified gravity models, which replace the standard Einstein-Hilbert action in General Relativity with a higher derivative theory. We use data from the WiggleZ Dark Energy survey to probe the formation of structure on large scales which can place tight constraints on these models. We combine the large-scale structure data with measurements of the cosmic microwave background from the Planck surveyor. After parameterising the modification of the action using the Compton wavelength parameter $B_0$, we constrain this parameter using ISiTGR, assuming an initial non-informative log prior probability distribution of this cross-over scale. We find that the addition of the WiggleZ power spectrum provides the tightest constraints to date on $B_0$ by an order of magnitude, giving ${\rm log}_{10}(B_0) < -4.07$ at 95% confidence limit. Finally, we test whether the effect of adding the lensing amplitude $A_{\rm Lens}$ and the sum of the neutrino mass $\sum m_\nu$ is able to reconcile current tensions present in these parameters, but find $f(R)$ gravity an inadequate explanation.

We present the strongest current cosmological upper limit on the sum of neutrino masses of < 0.18 (95% confidence). It is obtained by adding observations of the large-scale matter power spectrum from the WiggleZ Dark Energy Survey to observations of the cosmic microwave background data from the Planck surveyor, and measurements of the baryon acoustic oscillation scale. The limit is highly sensitive to the priors and assumptions about the neutrino scenario. We explore scenarios with neutrino masses close to the upper limit (degenerate masses), neutrino masses close to the lower limit where the hierarchy plays a role, and addition of massive or massless sterile species.

Combining galaxy cluster data from the ROSAT All-Sky Survey and the Chandra X-ray Observatory, cosmic microwave background data from the Wilkinson Microwave Anisotropy Probe, and galaxy clustering data from the WiggleZ Dark Energy Survey, the 6-degree Field Galaxy Survey and the Sloan Digital Sky Survey III, we test for consistency the cosmic growth of structure predicted by General Relativity (GR) and the cosmic expansion history predicted by the cosmological constant plus cold dark matter paradigm (LCDM). The combination of these three independent, well studied measurements of the evolution of the mean energy density and its fluctuations is able to break strong degeneracies between model parameters. We model the key properties of cosmic growth with the normalization of the matter power spectrum, sigma_8, and the cosmic growth index, gamma, and those of cosmic expansion with the mean matter density, Omega_m, the Hubble constant, H_0, and a kinematical parameter equivalent to that for the dark energy equation of state, w. For a spatially flat geometry, w=-1, and allowing for systematic uncertainties, we obtain sigma_8=0.785+-0.019 and gamma=0.570+0.064-0.063 (at the 68.3 per cent confidence level). Allowing both w and gamma to vary we find w=-0.950+0.069-0.070 and gamma=0.533+-0.080. To further tighten the constraints on the expansion parameters, we also include supernova, Cepheid variable and baryon acoustic oscillation data. For w=-1, we have gamma=0.616+-0.061. For our most general model with a free w, we measure Omega_m=0.278+0.012-0.011, H_0=70.0+-1.3 km s^-1 Mpc^-1 and w=-0.987+0.054-0.053 for the expansion parameters, and sigma_8=0.789+-0.019 and gamma=0.604+-0.078 for the growth parameters. These results are in excellent agreement with GR+LCDM (gamma~0.55; w=-1) and represent the tightest and most robust simultaneous constraint on cosmic growth and expansion to date.

We review the use of Bayesian Model Averaging in astrophysics. We first introduce the statistical basis of Bayesian Model Selection and Model Averaging. We discuss methods to calculate the model-averaged posteriors, including Markov Chain Monte Carlo (MCMC), nested sampling, Population Monte Carlo, and Reversible Jump MCMC. We then review some applications of Bayesian Model Averaging in astrophysics, including measurements of the dark energy and primordial power spectrum parameters in cosmology, cluster weak lensing and Sunyaev-Zel'dovich effect data, estimating distances to Cepheids, and classifying variable stars.

Neutrinos are one of the major puzzles in modern physics. Despite measurements of mass differences, the Standard Model of particle physics describes them as exactly massless. Additionally, recent measurements from both particle physics experiments and cosmology indicate the existence of more than the three Standard Model species. Here we review the cosmological evidence and its possible interpretations.

Recent indications from both particle physics and cosmology suggest the existence of more than three neutrino species. In cosmological analyses the effects of neutrino mass and number of species can in principle be disentangled for fixed cosmological parameters. However, since we do not have perfect measurements of the standard Lambda Cold Dark Matter model parameters some correlation remains between the neutrino mass and number of species, and both parameters should be included in the analysis. Combining the newest observations of several cosmological probes (cosmic microwave background, large scale structure, expansion rate) we obtain N_eff=3.58(+0.15/-0.16 at 68% CL) (+0.55/-0.53 at 95% CL) and a sum of neutrino masses of less than 0.60 eV (95 CL), which are currently the strongest constraints on N_eff and M_nu from an analysis including both parameters. The preference for N_eff >3 is now at a 2sigma level.

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 present an analysis of supernova light curves simulated for the upcoming Dark Energy Survey (DES) supernova search. The simulations employ a code suite that generates and fits realistic light curves in order to obtain distance modulus/redshift pairs that are passed to a cosmology fitter. We investigated several different survey strategies including field selection, supernova selection biases, and photometric redshift measurements. Using the results of this study, we chose a 30 square degree search area in the griz filter set. We forecast 1) that this survey will provide a homogeneous sample of up to 4000 Type Ia supernovae in the redshift range 0.05<z<1.2, and 2) that the increased red efficiency of the DES camera will significantly improve high-redshift color measurements. The redshift of each supernova with an identified host galaxy will be obtained from spectroscopic observations of the host. A supernova spectrum will be obtained for a subset of the sample, which will be utilized for control studies. In addition, we have investigated the use of combined photometric redshifts taking into account data from both the host and supernova. We have investigated and estimated the likely contamination from core-collapse supernovae based on photometric identification, and have found that a Type Ia supernova sample purity of up to 98% is obtainable given specific assumptions. Furthermore, we present systematic uncertainties due to sample purity, photometric calibration, dust extinction priors, filter-centroid shifts, and inter-calibration. We conclude by estimating the uncertainty on the cosmological parameters that will be measured from the DES supernova data.

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 demonstrate a methodology for optimizing the ability of future dark energy surveys to answer model selection questions, such as `Is acceleration due to a cosmological constant or a dynamical dark energy model?'. Model selection Figures of Merit are defined, exploiting the Bayes factor, and surveys optimized over their design parameter space via a Monte Carlo method. As a specific example we apply our methods to generic multi-fibre baryon acoustic oscillation spectroscopic surveys, comparable to that proposed for SuMIRe PFS, and present implementations based on the Savage-Dickey Density Ratio that are both accurate and practical for use in optimization. It is shown that whilst the optimal surveys using model selection agree with those found using the Dark Energy Task Force (DETF) Figure of Merit, they provide better informed flexibility of survey configuration and an absolute scale for performance; for example, we find survey configurations with close to optimal model selection performance despite their corresponding DETF Figure of Merit being at only 50% of its maximum. This Bayes factor approach allows us to interpret the survey configurations that will be good enough for the task at hand, vital especially when wanting to add extra science goals and in dealing with time restrictions or multiple probes within the same project.

Our goals are (i) to search for BAO and large-scale structure in current QSO survey data and (ii) to use these and simulation/forecast results to assess the science case for a new, >10x larger, QSO survey. We first combine the SDSS, 2QZ and 2SLAQ surveys to form a survey of ~60000 QSOs. We find a hint of a peak in the QSO 2-point correlation function, xi(s), at the same scale (~105h^-1 Mpc) as detected by Eisenstein et al (2005) in their sample of DR5 LRGs but only at low statistical significance. We then compare these data with QSO mock catalogues from the Hubble Volume simulation used by Hoyle et al (2002) and find that both routes give statistical error estimates that are consistent at ~100h^-1 Mpc scales. Mock catalogues are then used to estimate the nominal survey size needed for a 3-4 sigma detection of the BAO peak. We find that a redshift survey of ~250000 z<2.2 QSOs is required over ~3000 deg^2. This is further confirmed by static log-normal simulations where the BAO are clearly detectable in the QSO power spectrum and correlation function. The nominal survey would on its own produce the first detection of, for example, discontinuous dark energy evolution in the so far uncharted 1<z<2.2 redshift range. A survey with ~50% higher QSO sky densities and 50% bigger area will give an ~6sigma BAO detection, leading to an error ~60% of the size of the BOSS error on the dark energy evolution parameter, w_a. Another important aim for a QSO survey is to place new limits on primordial non-Gaussianity at large scales, testing tentative evidence we have found for the evolution of the linear form of the combined QSO xi(s) at z~1.6. Such a QSO survey will also determine the gravitational growth rate at z~1.6 via z-space distortions, allow lensing tomography via QSO magnification bias while also measuring the exact luminosity dependence of small-scale QSO clustering.

We investigate the prior dependence of constraints on cosmic tensor perturbations. Commonly imposed is the strong prior of the single-field inflationary consistency equation, relating the tensor spectral index nT to the tensor-to-scalar ratio r. Dropping it leads to significantly different constraints on nT, with both positive and negative values allowed with comparable likelihood, and substantially increases the upper limit on r on scales k = 0.01 Mpc^-1 to 0.05 Mpc^-1, by a factor of ten or more. Even if the consistency equation is adopted, a uniform prior on r on one scale does not correspond to a uniform one on another; constraints therefore depend on the pivot scale chosen. We assess the size of this effect and determine the optimal scale for constraining the tensor amplitude, both with and without the consistency relation.

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

We constrain the evolution of Newton's constant using the growth rate of large-scale structure measured by the WiggleZ Dark Energy Survey in the redshift range $0.1 < z < 0.9$. We use this data in two ways. Firstly we constrain the matter density of the Universe, $\omms$ (assuming General Relativity), and use this to construct a diagnostic to detect the presence of an evolving Newton's constant. Secondly we directly measure the evolution of Newton's constant, $\Geff$, that appears in Modified Gravity theories, without assuming General Relativity to be true. The novelty of these approaches are that, contrary to other methods, they do not require knowledge of the expansion history of the Universe, $H(z)$, making them model independent tests. Our constraints for the second derivative of Newton's constant at the present day, assuming it is slowly evolving as suggested by Big Bang Nucleosynthesis constraints, using the WiggleZ data is $\ddotGeff(t_0)=-1.19\pm 0.95\cdot 10^{-20}h^2 \textrm{yr}^{-2}$, where $h$ is defined via $H_0=100 h km s^{-1}Mpc^{-1}$, while using both the WiggleZ and the Sloan Digital Sky Survey Luminous Red Galaxy (SDSS LRG) data is $\ddotGeff(t_0)=-3.6\pm 6.8\cdot 10^{-21}h^2 \rm{yr}^{-2}$, both being consistent with General Relativity. Finally, our constraint for the rms mass fluctuation $\sigma_8$ using the WiggleZ data is $\sigma_8=0.75 \pm 0.08$, while using both the WiggleZ and the SDSS LRG data $\sigma_8=0.77 \pm 0.07$, both in good agreement with the latest measurements from the Cosmic Microwave Background radiation.

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.

Cosmological parameter uncertainties are often stated assuming a particular model, neglecting the model uncertainty, even when Bayesian model selection is unable to identify a conclusive best model. Bayesian model averaging is a method for assessing parameter uncertainties in situations where there is also uncertainty in the underlying model. We apply model averaging to the estimation of the parameters associated with the primordial power spectra of curvature and tensor perturbations. We use CosmoNest and MultiNest to compute the model Evidences and posteriors, using cosmic microwave data from WMAP, ACBAR, BOOMERanG and CBI, plus large-scale structure data from the SDSS DR7. We find that the model-averaged 95% credible interval for the spectral index using all of the data is 0.940 < n_s < 1.000, where n_s is specified at a pivot scale 0.015 Mpc^{-1}. For the tensors model averaging can tighten the credible upper limit, depending on prior assumptions.

We report results from the Supernova Photometric Classification Challenge (SNPCC), a publicly released mix of simulated supernovae (SNe), with types (Ia, Ibc, and II) selected in proportion to their expected rate. The simulation was realized in the griz filters of the Dark Energy Survey (DES) with realistic observing conditions (sky noise, point-spread function and atmospheric transparency) based on years of recorded conditions at the DES site. Simulations of non-Ia type SNe are based on spectroscopically confirmed light curves that include unpublished non-Ia samples donated from the Carnegie Supernova Project (CSP), the Supernova Legacy Survey (SNLS), and the Sloan Digital Sky Survey-II (SDSS-II). A spectroscopically confirmed subset was provided for training. We challenged scientists to run their classification algorithms and report a type and photo-z for each SN. Participants from 10 groups contributed 13 entries for the sample that included a host-galaxy photo-z for each SN, and 9 entries for the sample that had no redshift information. Several different classification strategies resulted in similar performance, and for all entries the performance was significantly better for the training subset than for the unconfirmed sample. For the spectroscopically unconfirmed subset, the entry with the highest average figure of merit for classifying SNe~Ia has an efficiency of 0.96 and an SN~Ia purity of 0.79. As a public resource for the future development of photometric SN classification and photo-z estimators, we have released updated simulations with improvements based on our experience from the SNPCC, added samples corresponding to the Large Synoptic Survey Telescope (LSST) and the SDSS, and provided the answer keys so that developers can evaluate their own analysis.

Future photometric supernova surveys will produce vastly more candidates than can be followed up spectroscopically, highlighting the need for effective classification methods based on lightcurves alone. Here we introduce boosting and kernel density estimation techniques which have minimal astrophysical input, and compare their performance on 20,000 simulated Dark Energy Survey lightcurves. We demonstrate that these methods are comparable to the best template fitting methods currently used, and in particular do not require the redshift of the host galaxy or candidate. However both methods require a training sample that is representative of the full population, so typical spectroscopic supernova subsamples will lead to poor performance. To enable the full potential of such blind methods, we recommend that representative training samples should be used and so specific attention should be given to their creation in the design phase of future photometric surveys.

We explore unification of dark matter and dark energy in a theory containing a scalar field of non-Lagrangian type, obtained by direct insertion of a kinetic term into the energy-momentum tensor. This scalar is different from quintessence, having an equation of state between -1 and 0 and a zero sound speed in its rest frame. We solve the equations of motion for an exponential potential via a rewriting as an autonomous system, and demonstrate the observational viability of the scenario, for sufficiently small exponential potential parameter \lambda, by comparison to a compilation of kinematical cosmological data.

Cosmological observations are normally fit under the assumption that the dark sector can be decomposed into dark matter and dark energy components. However, as long as the probes remain purely gravitational, there is no unique decomposition and observations can only constrain a single dark fluid; this is known as the dark degeneracy. We use observations to directly constrain this dark fluid in a model-independent way, demonstrating in particular that the data cannot be fit by a dark fluid with a single constant equation of state. Parameterizing the dark fluid equation of state by a variety of polynomials in the scale factor $a$, we use current kinematical data to constrain the parameters. While the simplest interpretation of the dark fluid remains that it is comprised of separate dark matter and cosmological constant contributions, our results cover other model types including unified dark energy/matter scenarios.

We extend our study of the optimization of large baryon acoustic oscillation (BAO) surveys to return the best constraints on the dark energy, building on Paper I of this series (Parkinson et al. 2007). The survey galaxies are assumed to be pre-selected active, star-forming galaxies observed by their line emission with a constant number density across the redshift bin. Star-forming galaxies have a redshift desert in the region 1.6 < z < 2, and so this redshift range was excluded from the analysis. We use the Seo & Eisenstein (2007) fitting formula for the accuracies of the BAO measurements, using only the information for the oscillatory part of the power spectrum as distance and expansion rate rulers. We go beyond our earlier analysis by examining the effect of including curvature on the optimal survey configuration and updating the expected `prior' constraints from Planck and SDSS. We once again find that the optimal survey strategy involves minimizing the exposure time and maximizing the survey area (within the instrumental constraints), and that all time should be spent observing in the low-redshift range (z<1.6) rather than beyond the redshift desert, z>2. We find that when assuming a flat universe the optimal survey makes measurements in the redshift range 0.1 < z <0.7, but that including curvature as a nuisance parameter requires us to push the maximum redshift to 1.35, to remove the degeneracy between curvature and evolving dark energy. The inclusion of expected other data sets (such as WiggleZ, BOSS and a stage III SN-Ia survey) removes the necessity of measurements below redshift 0.9, and pushes the maximum redshift up to 1.5. We discuss considerations in determining the best survey strategy in light of uncertainty in the true underlying cosmological model.

One of the biggest problems faced by current and next-generation astronomical surveys is trying to produce large numbers of accurate cross identifications across a range of wavelength regimes with varying data quality and positional uncertainty. Until recently simple spatial "nearest neighbour" associations have been sufficient for most applications. However as advances in instrumentation allow more sensitive images to be made the rapid increase in the source density has meant that source confusion across multiple wavelengths is a serious problem. The field of far-IR and sub-mm astronomy has been particularly hampered by such problems. The poor angular resolution of current sub-mm and far-IR instruments is such that in a lot of cases there are multiple plausible counterparts for each source at other wavelengths. Here we present a new automated method of producing associations between sources at different wavelengths using a combination of spatial and SED information set in the Bayesian framework presented by Budavari & Szalay (2008). Testing of the technique is performed on both simulated catalogues of sources from GaLICS and real data from multi-wavelength observations of the SXDF. It is found that a single figure of merit, the Bayes factor, can be effectively used to describe the confidence in the match. Further applications of this technique to future Herschel datasets are discussed.

Although cosmic microwave background (CMB) anisotropy data alone cannot constrain simultaneously the spatial curvature and the equation of state of dark energy, CMB data provide a valuable addition to other experimental results. However computing a full CMB power spectrum with a Boltzmann code is quite slow; for instance if we want to work with many dark energy and/or modified gravity models, or would like to optimize experiments where many different configurations need to be tested, it is possible to adopt a quicker and more efficient approach. In this paper we consider the compression of the projected Planck CMB data into four parameters, R (scaled distance to last scattering surface), l_a (angular scale of sound horizon at last scattering), Omega_b h^2 (baryon density fraction) and n_s (powerlaw index of primordial matter power spectrum), all of which can be computed quickly. We show that, although this compression loses information compared to the full likelihood, such information loss becomes negligible when more data is added. We also demonstrate that the method can be used for scalar field dark energy independently of the parametrisation of the equation of state, and discuss how this method should be used for other kinds of dark energy models.

We use model selection forecasting to assess the ability of the Planck satellite to make a positive detection of spectral index running. We simulate Planck data for a range of assumed cosmological parameter values, and carry out a three-way Bayesian model comparison of a Harrison-Zel'dovich model, a power-law model, and a model including running. We find that Planck will be able to strongly support running only if its true value satisfies |dn/d ln k| > 0.02.

We consider the extent to which future imaging surveys of galaxies can distinguish between dark energy and modified gravity models for the origin of the cosmic acceleration. Dynamical dark energy models may have similar expansion rates as models of modified gravity, yet predict different growth of structure histories. We parameterize the cosmic expansion by the two parameters, $w_0$ and $w_a$, and the linear growth rate of density fluctuations by Linder's $\gamma$, independently. Dark energy models generically predict $\gamma \approx 0.55$, while the DGP model $\gamma \approx 0.68$. To determine if future imaging surveys can constrain $\gamma$ within 20 percent (or $\Delta\gamma<0.1$), we perform the Fisher matrix analysis for a weak lensing survey such as the on-going Hyper Suprime-Cam (HSC) project. Under the condition that the total observation time is fixed, we compute the Figure of Merit (FoM) as a function of the exposure time $\texp$. We find that the tomography technique effectively improves the FoM, which has a broad peak around $\texp\simeq {\rm several}\sim 10$ minutes; a shallow and wide survey is preferred to constrain the $\gamma$ parameter. While $\Delta\gamma < 0.1$ cannot be achieved by the HSC weak-lensing survey alone, one can improve the constraints by combining with a follow-up spectroscopic survey like WFMOS and/or future CMB observations.

We optimize the design of future spectroscopic redshift surveys for constraining the dark energy via precision measurements of the baryon acoustic oscillations (BAO), with particular emphasis on the design of the Wide-Field Multi-Object Spectrograph (WFMOS). We develop a model that predicts the number density of possible target galaxies as a function of exposure time and redshift. We use this number counts model together with fitting formulae for the accuracy of the BAO measurements to determine the effectiveness of different surveys and instrument designs. We search through the available survey parameter space to find the optimal survey with respect to the dark energy equation-of-state parameters according to the Dark Energy Task Force Figure-of-Merit, including predictions of future measurements from the Planck satellite. We optimize the survey to test the LambdaCDM model, assuming that galaxies are pre-selected using photometric redshifts to have a constant number density with redshift, and using a non-linear cut-off for the matter power spectrum that evolves with redshift. We find that line-emission galaxies are strongly preferred as targets over continuum emission galaxies. The optimal survey covers a redshift range 0.8 < z < 1.4, over the widest possible area (6000 sq. degs from 1500 hours observing time). The most efficient number of fibres for the spectrograph is 2,000, and the survey performance continues to improve with the addition of extra fibres until a plateau is reached at 10,000 fibres. The optimal point in the survey parameter space is not highly peaked and is not significantly affected by including constraints from upcoming supernovae surveys and other BAO experiments.

In astro-ph/0702542v2, Linder and Miquel seek to criticize the use of Bayesian model selection for data analysis and for survey forecasting and design. Their discussion is based on three serious misunderstandings of the conceptual underpinnings and application of model-level Bayesian inference, which invalidate all their main conclusions. Their paper includes numerous further inaccuracies, including an erroneous calculation of the Bayesian Information Criterion. Here we seek to set the record straight.

We update and extend our previous work reconstructing the potential of a quintessence field from current observational data. We extend the cosmological dataset to include new supernova data, plus information from the cosmic microwave background and from baryon acoustic oscillations. We extend the modelling by considering Pad\'e approximant expansions as well as Taylor series, and by using observations to assess the viability of the tracker hypothesis. We find that parameter constraints have improved by a factor of two, with a strengthening of the preference of the cosmological constant over evolving quintessence models. Present data show some signs, though inconclusive, of favouring tracker models over non-tracker models under our assumptions.

We compute the Bayesian evidences for one- and two-parameter models of evolving dark energy, and compare them to the evidence for a cosmological constant, using current data from Type Ia supernova, baryon acoustic oscillations, and the cosmic microwave background. We use only distance information, ignoring dark energy perturbations. We find that, under various priors on the dark energy parameters, LambdaCDM is currently favoured as compared to the dark energy models. We consider the parameter constraints that arise under Bayesian model averaging, and discuss the implication of our results for future dark energy projects seeking to detect dark energy evolution. The model selection approach complements and extends the figure-of-merit approach of the Dark Energy Task Force in assessing future experiments, and suggests a significantly-modified interpretation of that statistic.

Model selection aims to determine which theoretical models are most plausible given some data, without necessarily asking about the preferred values of the model parameters. A common model selection question is to ask when new data require introduction of an additional parameter, describing a newly-discovered physical effect. We review several model selection statistics, and then focus on use of the Bayesian evidence, which implements the usual Bayesian analysis framework at the level of models rather than parameters. We describe our CosmoNest code, which is the first computationally-efficient implementation of Bayesian model selection in a cosmological context. We apply it to recent WMAP satellite data, examining the need for a perturbation spectral index differing from the scale-invariant (Harrison-Zel'dovich) case.

We use the three-year WMAP observations to determine the normalization of the matter power spectrum in inflationary cosmologies. In this context, the quantity of interest is not the normalization marginalized over all parameters, but rather the normalization as a function of the inflationary parameters n and r with marginalization over the remaining cosmological parameters. We compute this normalization and provide an accurate fitting function. The statistical uncertainty in the normalization is 3 percent, roughly half that achieved by COBE. We use the k-l relation for the standard cosmological model to identify the pivot scale for the WMAP normalization. We also quote the inflationary energy scale corresponding to the WMAP normalization.

The recent WMAP3 results have placed measurements of the spectral index n_S in an interesting position. While parameter estimation techniques indicate that the Harrison-Zel'dovich spectrum n_S=1 is strongly excluded (in the absence of tensor perturbations), Bayesian model selection techniques reveal that the case against n_S=1 is not yet conclusive. In this paper, we forecast the ability of the Planck satellite mission to use Bayesian model selection to convincingly exclude (or favour) the Harrison-Zel'dovich model.

We present a Bayesian model selection analysis of WMAP3 data using our code CosmoNest. We focus on the density perturbation spectral index $n_S$ and the tensor-to-scalar ratio $r$, which define the plane of slow-roll inflationary models. We find that while the Bayesian evidence supports the conclusion that $n_S \neq 1$, the data are not yet powerful enough to do so at a strong or decisive level. If tensors are assumed absent, the current odds are approximately 8 to 1 in favour of $n_S \neq 1$ under our assumptions, when WMAP3 data is used together with external data sets. WMAP3 data on its own is unable to distinguish between the two models. Further, inclusion of $r$ as a parameter weakens the conclusion against the Harrison-Zel'dovich case (n_S = 1, r=0), albeit in a prior-dependent way. In appendices we describe the CosmoNest code in detail, noting its ability to supply posterior samples as well as to accurately compute the Bayesian evidence. We make a first public release of CosmoNest, now available at http://www.cosmonest.org.

We introduce a statistical measure of the effective model complexity, called the Bayesian complexity. We demonstrate that the Bayesian complexity can be used to assess how many effective parameters a set of data can support and that it is a useful complement to the model likelihood (the evidence) in model selection questions. We apply this approach to recent measurements of cosmic microwave background anisotropies combined with the Hubble Space Telescope measurement of the Hubble parameter. Using mildly non-informative priors, we show how the 3-year WMAP data improves on the first-year data by being able to measure both the spectral index and the reionization epoch at the same time. We also find that a non-zero curvature is strongly disfavored. We conclude that although current data could constrain at least seven effective parameters, only six of them are required in a scheme based on the Lambda-CDM concordance cosmology.

A key science goal of upcoming dark energy surveys is to seek time evolution of the dark energy. This problem is one of {\em model selection}, where the aim is to differentiate between cosmological models with different numbers of parameters. However, the power of these surveys is traditionally assessed by estimating their ability to constrain parameters, which is a different statistical problem. In this paper we use Bayesian model selection techniques, specifically forecasting of the Bayes factors, to compare the abilities of different proposed surveys in discovering dark energy evolution. We consider six experiments -- supernova luminosity measurements by the Supernova Legacy Survey, SNAP, JEDI, and ALPACA, and baryon acoustic oscillation measurements by WFMOS and JEDI -- and use Bayes factor plots to compare their statistical constraining power. The concept of Bayes factor forecasting has much broader applicability than dark energy surveys.

The abundance of new cosmological data becoming available means that a wider range of cosmological models are testable than ever before. However, an important distinction must be made between parameter fitting and model selection. While parameter fitting simply determines how well a model fits the data, model selection statistics, such as the Bayesian Evidence, are now necessary to choose between these different models, and in particular to assess the need for new parameters. We implement a new evidence algorithm known as nested sampling, which combines accuracy, generality of application and computational feasibility, and apply it to some cosmological datasets and models. We find that a five-parameter model with Harrison-Zel'dovich initial spectrum is currently preferred.

The next generation of galaxy surveys will attempt to measure the baryon oscillations in the clustering power spectrum with high accuracy. These oscillations encode a preferred scale which may be used as a standard ruler to constrain cosmological parameters and dark energy models. In this paper we present simple analytical fitting formulae for the accuracy with which the preferred scale may be determined in the tangential and radial directions by future spectroscopic and photometric galaxy redshift surveys. We express these accuracies as a function of survey parameters such as the central redshift, volume, galaxy number density and (where applicable) photometric redshift error. These fitting formulae should greatly increase the efficiency of optimizing future surveys, which requires analysis of a potentially vast number of survey configurations and cosmological models. The formulae are calibrated using a grid of Monte Carlo simulations, which are analyzed by dividing out the overall shape of the power spectrum before fitting a simple decaying sinusoid to the oscillations. The fitting formulae reproduce the simulation results with a fractional scatter of 7% (10%) in the tangential (radial) directions over a wide range of input parameters. We also indicate how sparse-sampling strategies may enhance the effective survey area if the sampling scale is much smaller than the projected baryon oscillation scale.

We highlight the flexibility of the IPSO experiment-design framework by contrasting its application to CMB, weak lensing and redshift surveys. We illustrate the latter with a 10 parameter MCMC D-optimisation of a dark energy redshift survey. When averaged over a standard dark energy model space the resulting optimal survey typically has only one or two redshift bins, located at z<2. By exploiting optimisation we show how the statistical power of such surveys is significantly enhanced. Experiment design is aided by the richness of the figure of merit landscape which means one can impose secondary optimisation criteria at little cost. For example, one may choose either to maximally test a single model (such as \Lambda CDM) or to get the most general model-independent constraints possible (e.g. on a whole space of dark energy models). Such freedom points to a future where cosmological experiments become increasingly specialised and optimisation increasingly important.

We describe an algorithm which directly determines the quintessence potential from observational data, without using an equation of state parametrisation. The strategy is to numerically determine observational quantities as a function of the expansion coefficients of the quintessence potential, which are then constrained using a likelihood approach. We further impose a model selection criterion, the Bayesian Information Criterion, to determine the appropriate level of the potential expansion. In addition to the potential parameters, the present-day quintessence field velocity is kept as a free parameter. Our investigation contains unusual model types, including a scalar field moving on a flat potential, or in an uphill direction, and is general enough to permit oscillating quintessence field models. We apply our method to the `gold' Type Ia supernovae sample of Riess et al. (2004), confirming the pure cosmological constant model as the best description of current supernovae luminosity-redshift data. Our method is optimal for extracting quintessence parameters from future data.

Detecting dark energy dynamics is the main quest of current dark energy research. Addressing the issue demands a fully consistent analysis of CMB, large scale structure and SN-Ia data with multi-parameter freedom valid for all redshifts. Here we undertake a ten parameter analysis of general dark energy confronted with the first year WMAP, 2dF galaxy survey and latest SN-Ia data. Despite the huge freedom in dark energy dynamics there are no new degeneracies with standard cosmic parameters apart from a mild degeneracy between reionisation and the redshift of acceleration, both of which effectively suppress small scale power. Breaking this degeneracy will help significantly in detecting dynamics, if it exists. Our best-fit model to the data has significant late-time evolution at z<1.5. Phantom models are also considered and we find that the best-fit crosses w=-1 which, if confirmed, would be a clear signal for radically new physics. Treatment of such rapidly varying models requires careful integration of the dark energy density usually not implemented in standard codes, leading to crucial errors of up to 5%. Nevertheless cosmic variance means that standard LCDM models are still a very good fit to the data and evidence for dynamics is currently very weak. Independent tests of reionisation or the epoch of acceleration (e.g. ISW-LSS correlations) or reduction of cosmic variance at large scales (e.g. cluster polarisation at high redshift) may prove key in the hunt for dynamics.

With the WMAP data we can now begin to test realistic models of inflation involving multiple scalar fields. These naturally lead to correlated adiabatic and isocurvature (entropy) perturbations with a running spectral index. We present the first full (9 parameter) likelihood analysis of double inflation with WMAP data and find that despite the extra freedom, supersymmetric hybrid potentials are strongly constrained with less than 7% correlated isocurvature component allowed when standard priors are imposed on the cosomological parameters. As a result we also find that Akaike & Bayesian model selection criteria rather strongly prefer single-field inflation, just as equivalent analysis prefers a cosmological constant over dynamical dark energy in the late universe. It appears that simplicity is the best guide to our universe.

By combining the recent WMAP measurements of the cosmic microwave background anisotropies and the results of the recent luminosity distance measurements to type-Ia supernovae, we find that the normalization of the matter power spectrum on cluster scales, sigma_8, can be used to discriminate between dynamical models of dark energy (quintessence models) and a conventional cosmological constant model (LCDM).

Dark energy affects the CMB through its perturbations and affects both CMB and Sn-Ia through its background evolution. Using recent CMB and Sn-Ia data sets, together with the most general parameterization of the dark energy equation of state available, we find that today w < -0.8 (2 sigma). We also find that the value of the normalization of the power spectrum on cluster scales, sigma8, can be used to discriminate between dynamical models of dark energy (Quintessence models) and a cosmological constant model (Lambda CDM).

Cosmological dark energy is a natural source of variation of the fine structure constant. Using a model-independent approach we show that once general assumptions about the alpha-varying interactions are made, astronomical probes of its variation constrain the dark energy equation of state today to satisfy -1 < w_f < -0.96 at 3-sigma and significantly disfavour late-time changes in the equation of state. We show how dark-energy-induced spatial perturbations of alpha are linked to violations of the Equivalence Principle and are thus negligible at low-redshift, in stark contrast to the BSBM theories. This provides a new test of dark energy as the source of alpha variation.

Imagine a scenario in which the dark energy forms via the condensation of dark matter at some low redshift. The Compton wavelength therefore changes from small to very large at the transition, unlike quintessence or metamorphosis. We study CMB, large scale structure, supernova and radio galaxy constraints on condensation by performing a 4 parameter likelihood analysis over the Hubble constant and the three parameters associated with Q, the condensate field: Omega_Q, w_f and z_t (energy density and equation of state today, and redshift of transition). Condensation roughly interpolates between Lambda CDM (for large z_t) and sCDM (low z_t) and provides a slightly better fit to the data than Lambda CDM. We confirm that there is no degeneracy in the CMB between H and z_t and discuss the implications of late-time transitions for the Lyman-alpha forest. Finally we discuss the nonlinear phase of both condensation and metamorphosis, which is much more interesting than in standard quintessence models.

We show explicitly some exciting features of double-inflation: (i) it can often lead to strongly correlated adiabatic and entropy (isocurvature) power spectra. (ii) The two-field slow-roll consistency relations can be violated when the correlation is large at Hubble crossing. (iii) The spectra of adiabatic and entropy perturbations can be strongly scale-dependent and tilted toward either the red or blue. These effects are typically due to a light or time-dependent entropy mass and a non-negligible angular velocity in field space during inflation. They are illustrated via a multi-parameter numerical search for correlations in two concrete models. The correlation is found to be particularly strong in a supersymmetric scenario due to rapid growth of entropy perturbations in the tachyonic region separating the two inflationary stages. Our analysis suggests that realistic double-inflation models will provide a rich and fruitful arena for the application of future cosmic data sets and new approximation schemes which go beyond slow-roll.

Until now all well-studied dark energy models proposed have very large Compton wavelengths, $\lambda_c > 100 Mpc$ at all time and hence show no clustering on small scales. Evidence against this would disfavour both $\Lambda$CDM and quintessence. We show however, that large $\lambda_c$ is actually favoured by current CMB, LSS and SN1a data by comparing quintessence with a condensate scenario with the same background dynamics. The condensate scenario has a very small $\lambda_c$ before making a transition to negative pressure and large $\lambda_c$ and therefore exhibits a rich CMB phenomenology which roughly interpolates between $\Lambda$CDM and flat CDM models.