We used the mark weighted correlation functions (MCFs), $W(s)$, to study the large scale structure of the Universe. We studied five types of MCFs with the weighting scheme $\rho^\alpha$, where $\rho$ is the local density, and $\alpha$ is taken as $-1,\ -0.5,\ 0,\ 0.5$, and 1. We found that different MCFs have very different amplitudes and scale-dependence. Some of the MCFs exhibit distinctive peaks and valleys that do not exist in the standard correlation functions. Their locations are robust against the redshifts and the background geometry, however it is unlikely that they can be used as ``standard rulers'' to probe the cosmic expansion history. Nonetheless we find that these features may be used to probe parameters related with the structure formation history, such as the values of $\sigma_8$ and the galaxy bias. Finally, after conducting a comprehensive analysis using the full shapes of the $W(s)$s and $W_{\Delta s}(\mu)$s, we found that, combining different types of MCFs can significantly improve the cosmological parameter constraints. Compared with using only the standard correlation function, the combinations of MCFs with $\alpha=0,\ 0.5,\ 1$ and $\alpha=0,\ -1,\ -0.5,\ 0.5,\ 1$ can improve the constraints on $\Omega_m$ and $w$ by $\approx30\%$ and $50\%$, respectively. We find highly significant evidence that MCFs can improve cosmological parameter constraints.

We propose a light-weight deep convolutional neural network (CNN) to estimate the cosmological parameters from simulated 3-dimensional dark matter distributions with high accuracy. The training set is based on 465 realizations of a cubic box with a side length of $256\ h^{-1}\ \rm Mpc$, sampled with $128^3$ particles interpolated over a cubic grid of $128^3$ voxels. These volumes have cosmological parameters varying within the flat $\Lambda$CDM parameter space of $0.16 \leq \Omega_m \leq 0.46$ and $2.0 \leq 10^9 A_s \leq 2.3$. The neural network takes as an input cubes with $32^3$ voxels and has three convolution layers, three dense layers, together with some batch normalization and pooling layers. In the final predictions from the network we find a $2.5\%$ bias on the primordial amplitude $\sigma_8$ that can not easily be resolved by continued training. We correct this bias to obtain unprecedented accuracy in the cosmological parameter estimation with statistical uncertainties of $\delta \Omega_m$=0.0015 and $\delta \sigma_8$=0.0029, which are several times better than the results of previous CNN works. Compared with a 2-point analysis method using clustering region of 0-130 and 10-130 $h^{-1}$ Mpc, the CNN constraints are several times and an order of magnitude more precise, respectively. Finally, we conduct preliminary checks of the error-tolerance abilities of the neural network, and find that it exhibits robustness against smoothing, masking, random noise, global variation, rotation, reflection, and simulation resolution. Those effects are well understood in typical clustering analysis, but had not been tested before for the CNN approach. Our work shows that CNN can be more promising than people expected in deriving tight cosmological constraints from the cosmic large scale structure.

The tomographic AP method is so far the best method in separating the Alcock-Paczynski (AP) signal from the redshift space distortion (RSD) effects and deriving powerful constraints on cosmological parameters using the $\lesssim40h^{-1}\ \rm Mpc$ clustering region. To guarantee that the method can be easily applied to the future large scale structure (LSS) surveys, we study the possibility of estimating the systematics of the method using fast simulation method. The major contribution of the systematics comes from the non-zero redshift evolution of the RSD effects, which is quantified by $\hat\xi_{\Delta s}(\mu,z)$ in our analysis, and estimated using the BigMultidark exact N-body simulation and approximate COLA simulation samples. We find about 5\%/10\% evolution when comparing the $\hat\xi_{\Delta s}(\mu,z)$ measured as $z=0.5$/$z=1$ to the measurements at $z=0$. We checked the inaccuracy in the 2pCFs computed using COLA, and find it 5-10 times smaller than the intrinsic systematics of the tomographic AP method, indicating that using COLA to estimate the systematics is good enough. Finally, we test the effect of halo bias, and find $\lesssim$1.5\% change in $\hat\xi_{\Delta s}$ when varying the halo mass within the range of $2\times 10^{12}$ to $10^{14}$ $M_{\odot}$. We will perform more studies to achieve an accurate and efficient estimation of the systematics in redshift range of $z=0-1.5$.

The tomographic Alcock-Paczynski (AP) method utilizes the redshift evolution of the AP distortion to place constraints on cosmological parameters. It has proved to be a robust method that can separate the AP signature from the redshift space distortion (RSD) effect, and deliver powerful cosmological constraints using the $\lesssim 40h^{-1}\ \rm Mpc$ clustering region. In previous works, the tomographic AP method was performed via the anisotropic 2-point correlation function statistic. In this work we consider the feasibility of conducting the analysis in the Fourier domain and examine the pros and cons of this approach. We use the integrated galaxy power spectrum (PS) as a function of direction, $\hat P_{\Delta k}(\mu)$, to quantify the magnitude of anisotropy in the large-scale structure clustering, and use its redshift variation to do the AP test. The method is tested on the large, high resolution Big-MultiDark Planck (BigMD) simulation at redshifts $z=0-1$, using the underlying true cosmology $\Omega_m=0.3071,\ w=-1$. Testing the redshift evolution of $\hat P_{\Delta k}(\mu)$ in the true cosmology and cosmologies deviating from the truth with $\delta \Omega_m=0.1,\ \delta w=0.3$, we find that the redshift evolution of the AP distortion overwhelms the effects created by the RSD by a factor of $\sim1.7-3.6$. We test the method in the range of $k\in(0.2,1.8)\ h\ \rm Mpc^{-1}$, and find that it works well throughout the entire regime. We tune the halo mass within the range $2\times 10^{13}$ to $10^{14}\ M_{\odot}$, and find that the change of halo bias results in $\lesssim 5 \%$ change in $\hat P_{\Delta k}(\mu)$, which is less significant compared with the cosmological effect. Our work shows that it is feasible to conduct the tomographic AP analysis in the Fourier space.

We develop an improved Alcock-Paczynski (AP) test method that uses the redshift-space two-point correlation function (2pCF) of galaxies. Cosmological constraints can be obtained by examining the redshift dependence of the normalized 2pCF, which should not change apart from the expected small non-linear evolution. An incorrect choice of cosmology used to convert redshift to comoving distance will manifest itself as redshift-dependent 2pCF. Our method decomposes the redshift difference of the two-dimensional correlation function into the Legendre polynomials whose amplitudes are modeled by radial fitting functions. Our likelihood analysis with this 2-D fitting scheme tightens the constraints on $\Omega_m$ and ${w}$ by $\sim 40\%$ compared to the method of Li et al. (2016, 2017, 2018) that uses one dimensional angular dependence only. We also find that the correction for the non-linear evolution in the 2pCF has a non-negligible cosmology dependence, which has been neglected in previous similar studies by Li et al.. With an accurate accounting for the non-linear systematics and use of full two-dimensional shape information of the 2pCF down to scales as small as $5~h^{-1}{\rm Mpc}$ it is expected that the AP test with redshift-space galaxy clustering anisotropy can be a powerful method to constrain the expansion history of the universe.

The tomographic Alcock-Paczynski (AP) method can result in tight cosmological constraints by using small and intermediate clustering scales of the large scale structure (LSS) of the galaxy distribution. By focusing on the redshift dependence, the AP distortion can be distinguished from the distortions produced by the redshift space distortions (RSD). In this work, we combine the tomographic AP method with other recent observational datasets of SNIa+BAO+CMB+$H_0$ to reconstruct the dark energy equation-of-state $w$ in a non-parametric form. The result favors a dynamical DE at $z\lesssim1$, and shows a mild deviation ($\lesssim2\sigma$) from $w=-1$ at $z=0.5-0.7$. We find the addition of the AP method improves the low redshift ($z\lesssim0.7$) constraint by $\sim50\%$.

We present an algorithm for the fast computation of the general $N$-point spatial correlation functions of any discrete point set embedded within an Euclidean space of $\mathbb{R}^n$. Utilizing the concepts of kd-trees and graph databases, we describe how to count all possible $N$-tuples in binned configurations within a given length scale, e.g. all pairs of points or all triplets of points with side lengths $<r_{max}$. Through bench-marking we show the computational advantage of our new graph based algorithm over more traditional methods. We show that all 3-point configurations up to and beyond the Baryon Acoustic Oscillation scale ($\sim$200 Mpc in physical units) can be performed on current SDSS data in reasonable time. Finally we present the first measurements of the 4-point correlation function of $\sim$0.5 million SDSS galaxies over the redshift range $0.43<z<0.7$.

We perform an anisotropic clustering analysis of 1,133,326 galaxies from the Sloan Digital Sky Survey (SDSS-III) Baryon Oscillation Spectroscopic Survey (BOSS) Data Release (DR) 12 covering the redshift range $0.15<z<0.69$. The geometrical distortions of the galaxy positions, caused by incorrect cosmological model assumptions, are captured in the anisotropic two-point correlation function on scales 6 -- 40 $h^{-1}\rm Mpc$. The redshift evolution of this anisotropic clustering is used to place constraints on the cosmological parameters. We improve the methodology of Li et al. 2016, to enable efficient exploration of high dimensional cosmological parameter spaces, and apply it to the Chevallier-Polarski-Linder parametrization of dark energy, $w=w_0+w_a{z}/({1+z})$. In combination with the CMB, BAO, SNIa and $H_0$ from Cepheid data, we obtain $\Omega_m = 0.301 \pm 0.008,\ w_0 = -1.042 \pm 0.067,\ $ and $w_a = -0.07 \pm 0.29$ (68.3\% CL). Adding our new AP measurements to the aforementioned results reduces the error bars by $\sim$30 -- 40\% and improves the dark energy figure of merit by a factor of $\sim$2. We check the robustness of the results using realistic mock galaxy catalogues.

Baryon Acoustic Oscillations (BAOs) in the early Universe are predicted to leave an as yet undetected signature on the relative clustering of total mass versus luminous matter. This signature, a modulation of the relative large-scale clustering of baryons and dark matter, offers a new angle to compare the large scale distribution of light versus mass. A detection of this effect would provide an important confirmation of the standard cosmological paradigm and constrain alternatives to dark matter as well as non-standard fluctuations such as Compensated Isocurvature Perturbations (CIPs). The first attempt to measure this effect in the SDSS-III BOSS Data Release 10 CMASS sample remained inconclusive but allowed to develop a method, which we detail here and use to conduct the second observational search. When using the same model as in our previous study and including CIPs in the model, the DR12 data are consistent with a null-detection, a result in tension with the strong evidence previously measured with the DR10 data. This tension remains when we use a more realistic model taking into account our knowledge of the survey flux limit, as the data then privilege a zero effect. In the absence of CIPs, we obtain a null detection consistent with both the absence of the effect and the amplitude predicted in previous theoretical studies. This shows the necessity of more accurate data in order to prove or disprove the theoretical predictions.

We develop a methodology to use the redshift dependence of the galaxy 2-point correlation function (2pCF) across the line-of-sight, $\xi(r_{\bot})$, as a probe of cosmological parameters. The positions of galaxies in comoving Cartesian space varies under different cosmological parameter choices, inducing a {\it redshift-dependent scaling} in the galaxy distribution. This geometrical distortion can be observed as a redshift-dependent rescaling in the measured $\xi(r_{\bot})$. We test this methodology using a sample of 1.75 billion mock galaxies at redshifts 0, 0.5, 1, 1.5, 2, drawn from the Horizon Run 4 N-body simulation. The shape of $\xi(r_{\bot})$ can exhibit a significant redshift evolution when the galaxy sample is analyzed under a cosmology differing from the true, simulated one. Other contributions, including the gravitational growth of structure, galaxy bias, and the redshift space distortions, do not produce large redshift evolution in the shape. We show that one can make use of this geometrical distortion to constrain the values of cosmological parameters governing the expansion history of the universe. This method could be applicable to future large scale structure surveys, especially photometric surveys such as DES, LSST, to derive tight cosmological constraints. This work is a continuation of our previous works as a strategy to constrain cosmological parameters using redshift-invariant physical quantities.

We present the first constraints on cosmology from the Dark Energy Survey (DES), using weak lensing measurements from the preliminary Science Verification (SV) data. We use 139 square degrees of SV data, which is less than 3\% of the full DES survey area. Using cosmic shear 2-point measurements over three redshift bins we find $\sigma_8 (\Omega_{\rm m}/0.3)^{0.5} = 0.81 \pm 0.06$ (68\% confidence), after marginalising over 7 systematics parameters and 3 other cosmological parameters. We examine the robustness of our results to the choice of data vector and systematics assumed, and find them to be stable. About $20$\% of our error bar comes from marginalising over shear and photometric redshift calibration uncertainties. The current state-of-the-art cosmic shear measurements from CFHTLenS are mildly discrepant with the cosmological constraints from Planck CMB data; our results are consistent with both datasets. Our uncertainties are $\sim$30\% larger than those from CFHTLenS when we carry out a comparable analysis of the two datasets, which we attribute largely to the lower number density of our shear catalogue. We investigate constraints on dark energy and find that, with this small fraction of the full survey, the DES SV constraints make negligible impact on the Planck constraints. The moderate disagreement between the CFHTLenS and Planck values of $\sigma_8 (\Omega_{\rm m}/0.3)^{0.5}$ is present regardless of the value of $w$.

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

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

We apply the methodology developed in \cite{Li2014,Li2015} to BOSS DR12 galaxies and derive cosmological constraints from the redshift dependence of the Alcock-Paczynski (AP) effect. The apparent anisotropy in the distribution of observed galaxies arise from two main sources, the redshift-space distortion (RSD) effect due to the galaxy peculiar velocities, and the geometric distortion when incorrect cosmological models are assumed for transforming redshift to comoving distance, known as the AP effect. Anisotropies produced by the RSD effect are, although large, maintaining a nearly uniform magnitude over a large range of redshift, while the degree of anisotropies from the AP effect varies with redshift by much larger magnitude. We split the DR12 galaxies into six redshift bins, measure the 2-point correlation function in each bin, and assess the redshift evolution of anisotropies. We obtain constraints of $\Omega_m=0.290 \pm 0.053,\ \ w = -1.07 \pm 0.15$, which are comparable with the current constraints from other cosmological probes such as type Ia supernovae, cosmic microwave background, and baryon acoustic oscillation (BAO). Combining these cosmological probes with our method yield tight constraints of $ \Omega_m = 0.301 \pm 0.006,\ w=-1.054 \pm 0.025$. Our method is complementary to the other large scale structure probes like BAO and topology. We expect this technique will play an important role in deriving cosmological constraints from large scale structure surveys.

We present measurements of weak gravitational lensing cosmic shear two-point statistics using Dark Energy Survey Science Verification data. We demonstrate that our results are robust to the choice of shear measurement pipeline, either ngmix or im3shape, and robust to the choice of two-point statistic, including both real and Fourier-space statistics. Our results pass a suite of null tests including tests for B-mode contamination and direct tests for any dependence of the two-point functions on a set of 16 observing conditions and galaxy properties, such as seeing, airmass, galaxy color, galaxy magnitude, etc. We furthermore use a large suite of simulations to compute the covariance matrix of the cosmic shear measurements and assign statistical significance to our null tests. We find that our covariance matrix is consistent with the halo model prediction, indicating that it has the appropriate level of halo sample variance. We compare the same jackknife procedure applied to the data and the simulations in order to search for additional sources of noise not captured by the simulations. We find no statistically significant extra sources of noise in the data. The overall detection significance with tomography for our highest source density catalog is 9.7sigma. Cosmological constraints from the measurements in this work are presented in a companion paper (DES et al. 2015).

We present measurements of the spatial clustering statistics in redshift space of various scalar field modified gravity simulations. We utilise the two-point and the three-point correlation functions to quantify the spatial distribution of dark matter halos within these simulations and thus discern between the models. We compare $\Lambda$CDM simulations to various modified gravity scenarios and find consistency with previous work in terms of 2-point statistics in real and redshift-space. However using higher order statistics such as the three-point correlation function in redshift space we find significant deviations from $\Lambda$CDM hinting that higher order statistics may prove to be a useful tool in the hunt for deviations from General Relativity.

The cosmic distance can be precisely determined using a `standard ruler' imprinted by primordial baryon acoustic oscillation (hereafter BAO) in the early Universe. The BAO at the targeted epoch is observed by analysing galaxy clustering in redshift space (hereafter RSD) for which a theoretical formulation is not yet fully understood, and thus makes this methodology unsatisfactory. The BAO analysis following a full RSD modelling is contaminated by systematic uncertainties due to a non-linear smearing effects such as non-linear corrections and by random viral velocity of galaxies. However, the BAO can be probed independently of RSD contamination using the BAO peak positions located in the 2D anisotropic correlation function. A new methodology is presented to measure peak positions, to test whether it is also contaminated by the same systematics in RSD, and to provide the radial and transverse cosmic distances determined by the 2D BAO peak positions. We find that in our model independent anisotropic clustering analysis we can obtain about $2\%$ and $5\%$ constraints on $D_A$ and $H^{-1}$ respectively with current BOSS data, which is competitive with other analysis.

Baryon Acoustic Oscillations (BAOs) in the early Universe are predicted to leave an as yet undetected signature on the relative clustering of total mass versus luminous matter. A detection of this effect would provide an important confirmation of the standard cosmological paradigm and constrain alternatives to dark matter as well as non-standard fluctuations such as Compensated Isocurvature Perturbations (CIPs). We conduct the first observational search for this effect, by comparing the number-weighted and luminosity-weighted correlation functions, using the SDSS-III BOSS Data Release 10 CMASS sample. When including CIPs in our model, we formally obtain evidence at $3.2\sigma$ of the relative clustering signature and a limit that matches the existing upper limits on the amplitude of CIPs. However, various tests suggest that these results are not yet robust, perhaps due to systematic biases in the data. The method developed in this Letter, used with more accurate future data such as that from DESI, is likely to confirm or disprove our preliminary evidence.

There are two redshifts in cosmology: $z_{obs}$, the observed redshift computed via spectral lines, and the model redshift, $z$, defined by the effective FLRW scale factor. In general these do not coincide. We place observational constraints on the allowed distortions of $z$ away from $z_{obs}$ - a possibility we dub redshift remapping. Remapping is degenerate with cosmic dynamics for either $d_L(z)$ or $H(z)$ observations alone: for example, the simple remapping $z = \alpha_1 z_{obs} +\alpha_2 z_{obs}^2$ allows a decelerating Einstein de Sitter universe to fit the observed supernova Hubble diagram as successfully as $\Lambda$CDM, highlighting that supernova data alone cannot prove that the universe is accelerating. We show however, that redshift remapping leads to apparent violations of cosmic distance duality that can be used to detect its presence even when neither a specific theory of gravity nor the Copernican Principle are assumed. Combining current data sets favours acceleration but does not yet rule out redshift remapping as an alternative to dark energy. Future surveys, however, will provide exquisite constraints on remapping and any models -- such as backreaction -- that predict it.

We analyse the clustering of cosmic large scale structure using a consistent modified gravity perturbation theory, accounting for anisotropic effects along and transverse to the line of sight. The growth factor has a particular scale dependence in f(R) gravity and we fit for the shape parameter f_{R0} simultaneously with the distance and the large scale (general relativity) limit of the growth function. Using more than 690,000 galaxies in the Baryon Oscillation Spectroscopy Survey Data Release 11, we find no evidence for extra scale dependence, with the 95\% confidence upper limit |f_{R0}| <8 \times 10^{-4}. Future clustering data, such as from the Dark Energy Spectroscopic Instrument, can use this consistent methodology to impose tighter constraints.

We propose a method using the redshift dependence of the Alcock-Paczynski (AP) test and volume effect to measure the cosmic expansion history. The galaxy two-point correlation function as a function of angle, $\xi(\mu)$, is measured at different redshifts. Assuming an incorrect cosmological model to convert galaxy redshifts to distances, the shape of $\xi(\mu)$ appears anisotropic due to the AP effect, and the amplitude shifted by the change in comoving volume. Due to the redshift dependence of the AP and volume effect, both the shape and amplitude of $\xi(\mu)$ exhibit redshift dependence. Similar to Li et.al (2014), we find the redshift-space distortions (RSD) caused by galaxy peculiar velocities, although significantly distorts $\xi(\mu)$, exhibit much less redshift evolution compared to the AP and volume effects. By focusing on the redshift dependence of $\xi(\mu)$, we can correctly recover the cosmological parameters despite the contamination of RSD. The method is tested by using the Horizon Run 3 N-body simulation, from which we made a series of $1/8$-sky mock surveys having 8 million physically self-bound halos and sampled to have roughly a uniform number density in $z=0-1.5$. We find the AP effect results in tight, unbiased constraints on the density parameter and dark energy equation of state, with 68.3% CL intervals $\delta \Omega_m\sim0.03$ and $\delta w\sim0.1$, and the volume effect leads to much tighter constraints of $\delta \Omega_m\sim0.007$ and $\delta w\sim0.035$.

We address the problem of separating stars from galaxies in future large photometric surveys. We focus our analysis on simulations of the Dark Energy Survey (DES). In the first part of the paper, we derive the science requirements on star/galaxy separation, for measurement of the cosmological parameters with the Gravitational Weak Lensing and Large Scale Structure probes. These requirements are dictated by the need to control both the statistical and systematic errors on the cosmological parameters, and by Point Spread Function calibration. We formulate the requirements in terms of the completeness and purity provided by a given star/galaxy classifier. In order to achieve these requirements at faint magnitudes, we propose a new method for star/galaxy separation in the second part of the paper. We first use Principal Component Analysis to outline the correlations between the objects parameters and extract from it the most relevant information. We then use the reduced set of parameters as input to an Artificial Neural Network. This multi-parameter approach improves upon purely morphometric classifiers (such as the classifier implemented in SExtractor), especially at faint magnitudes: it increases the purity by up to 20% for stars and by up to 12% for galaxies, at i-magnitude fainter than 23.

We present the measurements of the redshift-space three-point correlation function (3PCF) for z~0.5 luminous red galaxies of the CMASS sample in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Data Release 11. The 3PCF measurements are interpreted within the halo occupation distribution framework using high-resolution N-body simulations, and the model successfully reproduces the 3PCF on scales larger than 1Mpc/h. As with the case for the redshift-space two-point correlation functions, we find that the redshift-space 3PCF measurements also favour the inclusion of galaxy velocity bias in the model. In particular, the central galaxy in a halo is on average in motion with respect to the core of the halo. We discuss the potential of the small-scale 3PCF to tighten the constraints on the relation between galaxies and dark matter haloes and on the phase-space distribution of galaxies.

We analyze the clustering of large scale structure in the Universe in a model independent method, accounting for anisotropic effects along and transverse to the line of sight. The Baryon Oscillation Spectroscopy Survey Data Release 11 provides a large sample of 690,000 galaxies, allowing determination of the Hubble expansion H, angular distance D_A, and growth rate G_T at an effective redshift of z=0.57. After careful bias and convergence studies of the effects from small scale clustering, we find that cutting transverse separations below 40 Mpc/h delivers robust results while smaller scale data leads to a bias due to unmodelled nonlinear and velocity effects. The converged results are in agreement with concordance LCDM cosmology, general relativity, and minimal neutrino mass, all within the 68% confidence level. We also present results separately for the northern and southern hemisphere sky, finding a slight tension in the growth rate -- potentially a signature of anisotropic stress, or just covariance with small scale velocities -- but within 68% CL.

We present a one per cent measurement of the cosmic distance scale from the detections of the baryon acoustic oscillations in the clustering of galaxies from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). Our results come from the Data Release 11 (DR11) sample, containing nearly one million galaxies and covering approximately $8\,500$ square degrees and the redshift range $0.2<z<0.7$. We also compare these results with those from the publicly released DR9 and DR10 samples. Assuming a concordance $\Lambda$CDM cosmological model, the DR11 sample covers a volume of 13\,Gpc${}^3$ and is the largest region of the Universe ever surveyed at this density. We measure the correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of over $7\,\sigma$ in both the correlation function and power spectrum. Fitting for the position of the acoustic features measures the distance relative to the sound horizon at the drag epoch, $r_d$, which has a value of $r_{d,{\rm fid}}=149.28\,$Mpc in our fiducial cosmology. We find $D_V=(1264\pm25\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.32$ and $D_V=(2056\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.57$. At 1.0 per cent, this latter measure is the most precise distance constraint ever obtained from a galaxy survey. Separating the clustering along and transverse to the line-of-sight yields measurements at $z=0.57$ of $D_A=(1421\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ and $H=(96.8\pm3.4\,{\rm km/s/Mpc})(r_{d,{\rm fid}}/r_d)$. Our measurements of the distance scale are in good agreement with previous BAO measurements and with the predictions from cosmic microwave background data for a spatially flat cold dark matter model with a cosmological constant.

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

Our observations of the Universe are fundamentally anisotropic, with data from galaxies separated transverse to the line of sight coming from the same epoch while that from galaxies separated parallel to the line of sight coming from different times. Moreover, galaxy velocities along the line of sight change their redshift, giving redshift space distortions. We perform a full two-dimensional anisotropy analysis of galaxy clustering data, fitting in a substantially model independent manner the angular diameter distance D_A, Hubble parameter H, and growth rate ddelta/dln a without assuming a dark energy model. The results demonstrate consistency with LCDM expansion and growth, hence also testing general relativity. We also point out the interpretation dependence of the effective redshift z_eff, and its cosmological impact for next generation surveys.

We present non-linear weak lensing predictions for coupled dark energy models using the CoDECS simulations. We calculate the shear correlation function and error covariance expected for these models, for forthcoming ground-based (such as DES) and space-based (Euclid) weak lensing surveys. We obtain predictions for the discriminatory power of a ground-based survey similar to DES and a space-based survey such as Euclid in distinguishing between $\Lambda$CDM and coupled dark energy models; we show that using the non-linear lensing signal we could discriminate between $\Lambda$CDM and exponential constant coupling models with $\beta_0\geq0.1$ at $4\sigma$ confidence level with a DES-like survey, and $\beta_0\geq0.05$ at $5\sigma$ confidence level with Euclid. We also demonstrate that estimating the coupled dark energy models' non-linear power spectrum, using the $\Lambda$CDM Halofit fitting formula, results in biases in the shear correlation function that exceed the survey errors.

We present a catalogue of galaxy groups and clusters selected using a friends-of-friends algorithm with a dynamic linking length from the 2dF-SDSS and QSO (2SLAQ) luminous red galaxy survey. The linking parameters for the code are chosen through an analysis of simulated 2SLAQ haloes. The resulting catalogue includes 313 clusters containing 1,152 galaxies. The galaxy groups and clusters have an average velocity dispersion of sigma_v = 467.97 km/s and an average size of R_clt = 0.78 Mpc/h. Galaxies from regions of one square degree and centred on the galaxy clusters were downloaded from the Sloan Digital Sky Survey Data Release 6 (SDSS DR6). Investigating the photometric redshifts and cluster red-sequence of these galaxies shows that the galaxy clusters detected with the FoF algorithm are reliable out to z~0.6. We estimate masses for the clusters using their velocity dispersions. These mass estimates are shown to be consistent with 2SLAQ mock halo masses. Further analysis of the simulation haloes shows that clipping out low richness groups with large radii improves the purity of catalogue from 52% to 88%, while retaining a completeness of 94%. Finally, we test the two-point correlation function of our cluster catalogue. We find a best-fitting power law model with parameters r0 = 24\pm4 Mpc/h and gamma = -2.1\pm 0.2, which are in agreement with other low redshift cluster samples and consistent with a {\Lambda}CDM universe.

We present new measurements of the coherent motion of galaxies based on observations of the large-scale redshift-space distortions seen in the two-dimensional two-point correlation function of Luminous Red Galaxies in Data Release Seven of the Sloan Digital Sky Survey. We have developed a new methodology for estimating these coherent motions, which is less dependent on the details of galaxy bias and of the cosmological model to explain the late-time acceleration of the expansion of the Universe. We measure a one-dimensional velocity dispersion of galaxies on large-scales of \sigma_v=3.01^{+0.45}_{-0.46} Mpc/h and \sigma_v=3.69^{+0.47}_{-0.47} \mpcoh$ at a mean redshift of z=0.25 and 0.38 respectively. These values are fully consistent with predictions for a WMAP7-normalised LCDM Universe and inconsistent at confidence of 3.8\sigma with a Dvali-Gabadadze-Porrati (DGP) model for the Universe. We can convert the units of these $\sigma_v$ measurements to 270^{+40}_{-41} km/s and 320^{+41}_{-41} km/s respectively (assuming a $\Lambda$CDM universe), which are lower that expected based on recent low redshift (z<0.2) measurements of the peculiar velocity field (or "bulk flows"). It is difficult to directly compare these measurements as they cover different redshift ranges and different areas of the sky. However, one possible cosmological explanation for this discrepancy is that our Galaxy is located in unusually over, or under, dense region of the Universe.

We study the halo bispectrum from non-Gaussian initial conditions. Based on a set of large $N$-body simulations starting from initial density fields with local type non-Gaussianity, we find that the halo bispectrum exhibits a strong dependence on the shape and scale of Fourier space triangles near squeezed configurations at large scales. The amplitude of the halo bispectrum roughly scales as $f_nl^2$. The resultant scaling on the triangular shape is consistent with that predicted by Jeong & Komatsu based on perturbation theory. We systematically investigate this dependence with varying redshifts and halo mass thresholds. It is shown that the $f_nl$ dependence of the halo bispectrum is stronger for more massive haloes at higher redshifts. This feature can be a useful discriminator of inflation scenarios in future deep and wide galaxy redshift surveys.

We present here a new parameterization for the bulk motions of galaxies and clusters (in the linear regime) that can be measured statistically from the shape and amplitude of the two-dimensional two-point correlation function. We further propose the one-dimensional velocity dispersion (v_p) of the bulk flow as a complementary measure of redshift-space distortions, which is model-independent and not dependent on the normalisation method. As a demonstration, we have applied our new methodology to the C4 cluster catalogue constructed from Data Release Three (DR3) of the Sloan Digital Sky Survey. We find v_p=270^{+433}km/s (also consistent with v_p=0) for this cluster sample (at z=0.1), which is in agreement with that predicted for a WMAP5-normalised LCDM model (i.e., v_p(LCDM=203km/s). This measurement does not lend support to recent claims of excessive bulk motions (\simeq1000 km/s) which appear in conflict with LCDM, although our large statistical error cannot rule them out. From the measured coherent evolution of v_p, we develop a technique to re-construct the perturbed potential, as well as estimating the unbiased matter density fluctuations and scale--independent bias.

We present a new method for dealing with geometrical selection effects in galaxy surveys while using a multifractal framework. The power of multifractal analysis lies in its connection to higher order moments, in that it not only probes clustering on different scales but also different densities. Therefore any incompleteness issues must be correctly addressed before blindly applying this technique to real survey data.

Modern cosmological observations clearly reveal that the universe contains a hierarchy of clustering. However, recent surveys show a transition to homogeneity on large scales. The exact scale at which this transition occurs is still a topic of much debate. There has been much work done in trying to characterise the galaxy distribution using multifractals. However, for a number of years the size, depth and accuracy of galaxy surveys was regarded as insufficient to give a definitive answer. One of the main problems which arises in a multifractal analysis is how to deal with observational selection effects: ie `masks' in the survey region and a geometric boundary to the survey itself. In this thesis I will introduce a volume boundary correction which is rather similar to the approach developed by Pan and Coles in 2001, but which improves on their angular boundary correction in two important respects: firstly, our volume correction `throws away' fewer galaxies close the boundary of a given data set and secondly it is computationally more efficient. After application of our volume correction, I will then show how the underlying generalised dimensions of a given point set can be computed. I will apply this procedure to calculate the generalised fractal dimensions of both simulated fractal point sets and mock galaxy surveys which mimic the properties of the recent IRAS PSCz catalogue.