Abstract: Accurate power spectrum (or correlation function) covariance matrices are a crucial requirement for cosmological parameter estimation from large scale structure surveys. In order to minimize reliance on computationally expensive mock catalogs, it is important to have a solid analytic understanding of the different components that make up a covariance matrix. Considering the matter power spectrum covariance matrix, it has recently been found that there is a potentially dominant effect on mildly non-linear scales due to power in modes of size equal to and larger than the survey volume. This beat coupling effect has been derived analytically in perturbation theory and while it has been tested with simulations, some questions remain unanswered. Moreover, there is an additional effect of these large modes, which has so far not been included in analytic studies, namely the effect on the estimated average density which enters the power spectrum estimate. In this article, we work out analytic, perturbation theory based expressions including both the beat coupling and this local average effect and we show that while, when isolated, beat coupling indeed causes large excess covariance in agreement with the literature, in a realistic scenario this is compensated almost entirely by the local average effect, leaving only similar to 10% of the excess. We test our analytic expressions by comparison to a suite of large N-body simulations, using both full simulation boxes and subboxes thereof to study cases without beat coupling, with beat coupling and with both beat coupling and the local average effect. For the variances, we find excellent agreement with the analytic expressions for k < 0.2 hMpc(-1) at z = 0.5, while the correlation coefficients agree to beyond k = 0.4 hMpc(-1). As expected, the range of agreement increases towards higher redshift and decreases slightly towards z = 0. We finish by including the large-mode effects in a full covariance matrix description for arbitrary survey geometry and confirming its validity using simulations. This may be useful as a stepping stone towards building an actual galaxy (or other tracer's) power spectrum covariance matrix.

Abstract: We analyse the anisotropic clustering of massive galaxies from the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) sample, which consists of 264-283 galaxies in the redshift range 0.43 < z < 0.7 spanning 3275 deg(2). Both peculiar velocities and errors in the assumed redshiftdistance relation (AlcockPaczynski effect) generate correlations between clustering amplitude and orientation with respect to the line of sight. Together with the sharp baryon acoustic oscillation (BAO) standard ruler, our measurements of the broad-band shape of the monopole and quadrupole correlation functions simultaneously constrain the comoving angular diameter distance (2190 +/- 61 Mpc) to z = 0.57, the Hubble expansion rate at z = 0.57 (92.4 +/- 4.5 km s(-1) Mpc(-1)) and the growth rate of structure at that same redshift (d(sigma 8)/d ln a = 0.43 +/- 0.069). Our analysis provides the best current direct determination of both DA and H in galaxy clustering data using this technique. If we further assume a cold dark matter expansion history, our growth constraint tightens to d(sigma 8)/d ln a = 0.415 +/- 0.034. In combination with the cosmic microwave background, our measurements of D-A,H and d(sigma 8)/d ln a all separately require dark energy at z > 0.57, and when combined imply Omega(A) = 0.74 +/- 0.016, independent of the Universe's evolution at z < 0.57. All of these constraints assume scale-independent linear growth, and assume general relativity to compute both O(10 per cent) non-linear model corrections and our errors. In our companion paper, Samushia et al., we explore further cosmological implications of these observations.