Issue 
A&A
Volume 523, NovemberDecember 2010



Article Number  A1  
Number of page(s)  21  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/200913657  
Published online  10 November 2010 
Simultaneous measurement of cosmology and intrinsic alignments using joint cosmic shear and galaxy number density correlations
^{1}
Department of Physics and Astronomy, University College
London, London
WC1E 6BT,
UK
email: joachimi@astro.unibonn.de
^{2}
ArgelanderInstitut für Astronomie (AIfA), Universität
Bonn, Auf dem Hügel
71, 53121
Bonn,
Germany
Received:
12
November
2009
Accepted:
9
August
2010
Aims. Cosmic shear is a powerful method to constrain cosmology, provided that any systematic effects are under control. The intrinsic alignment of galaxies is expected to severely bias parameter estimates if not taken into account. We explore the potential of a joint analysis of tomographic galaxy ellipticity, galaxy number density, and ellipticitynumber density crosscorrelations to simultaneously constrain cosmology and selfcalibrate unknown intrinsic alignment and galaxy bias contributions.
Methods. We treat intrinsic alignments and galaxy biasing as free functions of scale and redshift and marginalise over the resulting parameter sets. Constraints on cosmology are calculated by combining the likelihoods from all twopoint correlations between galaxy ellipticity and galaxy number density. The information required for these calculations is already available in a standard cosmic shear data set. We include contributions to these functions from cosmic shear, intrinsic alignments, galaxy clustering and magnification effects.
Results. In a Fisher matrix analysis we compare our constraints with those from cosmic shear alone in the absence of intrinsic alignments. For a potential future large area survey, such as Euclid, the extra information from the additional correlation functions can make up for the additional free parameters in the intrinsic alignment and galaxy bias terms, depending on the flexibility in the models. For example, the dark energy task force figure of merit is recovered even when more than 100 free parameters are marginalised over. We find that the redshift quality requirements are similar to those calculated in the absence of intrinsic alignments.
Key words: cosmology: theory / gravitational lensing: weak / largescale structure of the Universe / cosmological parameters / methods: data analysis
© ESO, 2010
1. Introduction
On its way to Earth the light from distant galaxies is continuously deflected by the matter density inhomogeneities which it passes by. This induces distortions of the shapes of the projected galaxy images on the sky, causing modifications of these shapes of order 1%, known as weak gravitational lensing, or cosmic shear. Hence, to measure this distortion effect on distant galaxies, one requires statistical methods; see e.g. Bartelmann & Schneider (2001) and Schneider (2006) for detailed reviews.
Detected in 2000 (Bacon et al. 2000; Kaiser et al. 2000; van Waerbeke et al. 2000; Wittman et al. 2000), cosmic shear has since rapidly evolved into a mature technique that produces increasingly stringent constraints on cosmological parameters (see e.g. Jarvis et al. 2006; Hoekstra et al. 2006; Semboloni et al. 2006; Hetterscheidt et al. 2007; Benjamin et al. 2007; Massey et al. 2007; Schrabback et al. 2007; Fu et al. 2008; Schrabback et al. 2010). Probing both the evolution of structure and the geometry of the Universe at low redshifts, it is complementary to observations of the cosmic microwave background (e.g. Tereno et al. 2005; Das & Spergel 2009) and considered as potentially the most powerful method to pin down the properties of dark energy (Hu 2002; Albrecht et al. 2006; Peacock et al. 2006). Upcoming and future surveys like PanSTARRS^{1}, DES^{2}, LSST^{3}, JDEM^{4}, and Euclid^{5} will produce deep imaging over large fractions of the sky and thereby yield unprecedented insight into the properties of dark matter, dark energy and gravitation (e.g. Takada & Jain 2004; Réfrégier et al. 2004, 2006; Kitching et al. 2008a; Thomas et al. 2009).
The great statistical power of cosmic shear demands a careful assessment of possible systematic errors that might bias the results. A serious limitation may arise from a physical systematic caused by the intrinsic alignment of galaxies. The matter structure around galaxies can modify their intrinsic shape and their orientation. Firstly, this can result in correlations between the intrinsic shapes of galaxies which are close both on the sky and in redshift (intrinsic ellipticity correlations, or II correlations). Moreover, a dark matter halo can intrinsically align a physically close galaxy in the foreground and at the same time contribute to the lensing signal of a background object, which induces gravitational shearintrinsic ellipticity correlation (GI, Hirata & Seljak 2004).
Detailed investigations have been performed on the alignment between haloes (Croft & Metzler 2000; Heavens et al. 2000; Lee & Pen 2000; Catelan et al. 2001; Crittenden et al. 2001; Jing 2002; Mackey et al. 2002; Hirata & Seljak 2004; Bridle & Abdalla 2007; Schneider & Bridle 2010), as well as the alignment of the spin or the shape of a galaxy with its own dark matter halo (e.g. Pen et al. 2000; van den Bosch et al. 2002; Okumura et al. 2009; Okumura & Jing 2009; Brainerd et al. 2009; see also Schäfer 2009). Intrinsic alignments have also been investigated observationally, where nonvanishing II and GI signals have been detected in several surveys (Brown et al. 2002; Heymans et al. 2004; Mandelbaum et al. 2006; Hirata et al. 2007; Brainerd et al. 2009).
The results of both theoretical studies and observations show large variations, but most are consistent with a contamination of the order 10% by both II and GI correlations for future surveys that further divide the galaxy sample into redshift slices (cosmic shear tomography). Hence, the control of intrinsic alignments in cosmic shear studies is crucial to obtain unbiased results on cosmological parameters. Accurate models would solve the problem, but progress is hampered due to the dependence of intrinsic alignments on the intricacies of galaxy formation and evolution within their dark matter environment. Currently, the level of models is crude, and partly only motivated phenomenologically (see e.g. Schneider & Bridle 2010, for recent progress).
The II contamination can be controlled relatively easily by excluding close pairs of galaxies from the analysis (King & Schneider 2002, 2003; Heymans & Heavens 2003; Takada & White 2004). Joachimi & Schneider (2008, 2009) introduced a nulling technique which transforms the cosmic shear data vector and discards all entries of the transformed data set that are potentially contaminated by the GI signal. While this approach only relies on the wellknown redshift dependence of gravitational lensing, King (2005) projects out the GI term by making use of template functions. Furthermore the work by Mandelbaum et al. (2006) and Hirata et al. (2007) suggests that the intrinsic alignment may be dominated by luminous red galaxies which could be eliminated from the cosmic shear catalogues. All these removal techniques require excellent redshift information, and still they can cause a significant reduction in the constraints on cosmology.
Deep imaging surveys not only provide information about the shape of galaxies, but allow in addition for a measurement of galaxy number densities, as well as crosscorrelations between shape and number density information. This substantial extension of the set of observables increases the cosmological information to be extracted and, more importantly, enables one to internally calibrate systematic effects (Hu & Jain 2004; Bernstein 2009). By adding galaxy number density information one adds signals that are capable of pinning down the functional form of intrinsic alignments, but one also introduces as another systematic, the galaxy bias, which quantifies the lack of knowledge about how galaxies, i.e. the visible baryonic matter, follow the underlying dark matter distribution.
It is the scope of this work to elucidate the performance of a joint analysis of galaxy shape and number density information as regards the ability to constrain cosmological parameter in presence of general and flexible models of intrinsic alignments and galaxy bias. In doing so we incorporate several cosmological signals which have been considered before as promising probes of cosmology themselves, including galaxy clustering from photometric redshift surveys (Blake & Bridle 2005; Dolney et al. 2006; Zhan 2006; Blake et al. 2007; Padmanabhan et al. 2007) galaxygalaxy lensing (e.g. Schneider & Rix 1997; Guzik & Seljak 2001, 2002; Seljak 2002; Seljak et al. 2005; Yoo et al. 2006; Johnston et al. 2007; Cacciato et al. 2009) and lensing magnification (Broadhurst et al. 1995; Zhang & Pen 2005, 2006; van Waerbeke 2010). We follow the ansatz outlined in Bernstein (2009) and extend the investigation by Bridle & King (2007) who considered the residual information content in galaxy shape correlations after marginalising over the parameters of two loglinear grid models representing the II and GI terms. We quantify the crosscalibration properties of the joint set of observables and determine the requirements on cosmological surveys to efficiently apply this joint approach.
This paper is organised as follows: in Sect.2 we give an overview on the twopoint correlations that form part of the galaxy shape and number density observables, and we derive their explicit form. Two appendices provide further details. Section 3 demonstrates how we model the different signals and their dependence on cosmology. We introduce a general grid parametrisation for the intrinsic alignments and the galaxy bias. Furthermore we summarise our Fisher matrix formalism and the figures of merit we employ. In Sect.4 we present our results on the dependence of the parameter constraints on the freedom in the model of intrinsic alignments and galaxy bias, the characteristics of the redshift distributions, and the priors on the different sets of nuisance parameters. Finally, in Sect.5 we summarise our findings and conclude.
2. Twopoint correlations from cosmological surveys
Cosmological imaging surveys observe the angular positions and the projected shapes of huge numbers of galaxies over increasingly large areas on the sky. In addition, by means of multicolour photometry, it is possible to perform a tomographic analysis, i.e. obtain coarse information about the lineofsight dimension in terms of photometric redshifts (photoz). From the galaxy shapes in a given region of space, one can infer the ellipticity (1)where the superscript in parentheses assigns a photoz bin i. The observed ellipticity ϵ has contributions from the gravitational shear γ_{G} and an intrinsic shear γ_{I}, which is caused by the alignment of a galaxy in its surrounding gravitational field. Moreover, ϵ is assumed to have an uncorrelated component ϵ_{rnd}, which accounts for the purely random part of the intrinsic orientations and shapes of galaxies. Note that (1) is only valid if the gravitational shear is weak, see e.g. Seitz & Schneider (1997); Bartelmann & Schneider (2001) and for certain definitions of ellipticity.
Likewise, the positions of galaxies can be used to construct an estimate of the number density contrast (2)which is determined by the intrinsic number density contrast of galaxies n_{g} and the alteration of galaxy counts due to lensing magnification n_{m}. An uncorrelated shot noise contribution is added via n_{rnd}. In contrast to ϵ^{(i)}(θ) the number density contrast n^{(i)}(θ) can obviously not be estimated from individual galaxies. One can understand n^{(i)}(θ) as the ensemble average over a hypothetical, Poissondistributed random field of which the observed galaxy distribution is one particular representation. The formal relation between the projected number density contrast as used in (2) and the threedimensional galaxy number density fluctuations will be provided below, see (12).
As was already noted in Bernstein (2009), (1) and (2) are symmetric in the sense that they both contain an intrinsic contribution and a term caused by gravitational lensing effects. Under usual circumstances the correlated part of the ellipticity is dominated by the gravitational shear, whereas the largest term in (2) is due to the intrinsic number density contrast.
Both ellipticity and number overdensity vanish if averaged over sufficiently large scales. Thus, one considers to lowest order twopoint statistics of these quantities. Since all realspace twopoint measures are related to the power spectrum (see e.g. Kaiser 1992), we can work in terms of power spectra without loss of generality, which is desirable in particular due to a simpler structure of the signal covariances in Fourier space. Denoting the Fourier transform by a tilde, the power spectrum between redshift bins i and j can then be defined by (3)where is the twodimensional Dirac deltadistribution, and where ℓ denotes the angular frequency, the Fourier variable on the sky. The measures x_{a} and x_{b} can correspond to any of the set . The random contributions in (1) and (2) are not correlated with any of the other measures and only yield a contribution to the noise, see Sect.3.4.
Inserting (1) and (2) into (3), one obtains the complete set of tomographic twopoint observables which are available from shape and number density information see Bernstein (2009). We name signals stemming from galaxy shape information by capital letters (“G” for gravitational shear, “I” for intrinsic shear) and signals related to galaxy number densities by small letters (“g” for intrinsic number density fluctuations, “m” for lensing magnification). An overview of the nomenclature of the correlations in (4) to (6) is provided in Table 1. Note that (4) and (5) are symmetric with respect to their photoz bin arguments. Hence, if N_{zbin} denotes the number of available photoz bins, one has N_{zbin}(N_{zbin} + 1)/2 observables for every considered angular frequency. In contrast, one can exploit ellipticitynumber density crosscorrelation power spectra (6) per ℓ.
Overview on the twopoint correlations considered in this work.
Intrinsic galaxy clustering (gg) add a strong signal to the correlations of galaxy number densities (5), but its use to obtain cosmological parameter estimates is limited due by poor knowledge of the galaxy bias (e.g. Lahav & Suto 2004). Gravitational lensing modifies the flux of objects and thus reduces or increases number counts of galaxies above a certain limiting magnitude, depending on the form of the galaxy luminosity function close to the limiting magnitude. This produces magnification correlations (mm) and intrinsic number densitymagnification crosscorrelations (gm). The gm correlations occur when a foreground mass overdensity (underdensity) contains an overdensity (underdensity) of galaxies and (de)magnifies background objects along the same line of sight causing an apparent over or underdensity of galaxies at higher redshift.
Crosscorrelations between galaxy number densities and ellipticities (6) contain contributions from cross terms between intrinsic clustering and shear (gG), intrinsic clustering and intrinsic shear (gI), magnification and shear (mG), and magnification and intrinsic shear (mI). For instance, one expects to find gI and gG signals when a mass structure leads to an overdensity in the local galaxy distribution and influences the intrinsic shape of galaxies at the same redshift or contributes to the shear of background objects. The latter is the usual galaxygalaxy lensing signal. Because a foreground overdensity can in addition enhance galaxy counts due to lensing magnification, the mG and mI signals will also be nonvanishing. The form of all these correlations will be further discussed in Sect.3.
All nonrandom terms in (1) and (2), given for a photometric redshift bin i, can be related to a source term 𝒮 , which is a function of spatial coordinates, i.e. (7)where we defined a weight function w that depends on the photoz bin i (for a similar approach see Hu & Jain 2004). Here, χ denotes comoving distance, and f_{K}(χ) is the comoving angular diameter distance. If (7) holds for two quantities and , their projected power spectrum is given by the lineofsight integral of the threedimensional source power spectrum via Limber’s equation in Fourier space (Kaiser 1992), (8)By identifying weights and source terms for gravitational and intrinsic shear, as well as intrinsic clustering and magnification, we can derive Limber equations for all power spectra entering (4)–(6).
To compute the equivalent of (7) for the cosmic shear case, we first note that in Fourier space the shear and the convergence are related by the simple equation , where ϕ_{ℓ} is the polar angle of ℓ. As a consequence, the power spectra of shear and convergence are identical. Therefore, we can equivalently use the convergence κ^{(i)}(θ) as the cosmic shear observable. It is related to the threedimensional matter density contrast δ via (9)where the weight is given by (10)see Bartelmann & Schneider (2001); Schneider (2006) for details. Here a denotes the scale factor and p^{(i)}(χ) the comoving distance probability distribution of those galaxies in bin i for which shape information is available.
Analogously to the lensing case, one can define a convergence of the intrinsic shear field , which is directly related to the intrinsic shear via . This intrinsic convergence is a projection of the threedimensional intrinsic shear field , which can be written as (11)see e.g. Hirata & Seljak (2004) for the analogous expression in terms of intrinsic shear. Here we have assumed that the intrinsic shear field is – like the gravitational shear field – curlfree to good approximation. This holds for instance for the linear alignment model developed in Hirata & Seljak (2004). Then corresponds to the Fourier transform of as defined in Schneider & Bridle (2010).
Likewise, angular galaxy number density fluctuations are given by the lineofsight projection of threedimensional number density fluctuations δ_{g} as (e.g. Hu & Jain 2004) (12)It is important to note that p^{(i)}(χ) is the same as in (10) and (11), i.e. the number counts are restricted to those galaxies with shape measurements, which require a higher signaltonoise than the position determination. In principle, number density information could be obtained for a larger number of galaxies, in particular fainter ones. But, to determine the contribution to number density correlations by magnification, it is necessary to measure the slope of the luminosity function α^{(i)} at the faint end of the used galaxy distribution. We will detail the exact definition and the determination of α^{(i)} in Sect.3.2. Since it is desirable to extract the values of the α^{(i)} internally from the survey, one needs to be able to measure fluxes down to values slightly below the magnitude limit of the galaxies included in p^{(i)}(χ). Moreover galaxy number density measurements may require photometric redshifts which are of the same or better quality than for cosmic shear studies, limiting the number of faint usable galaxies. Hence, we argue that the choice of identical distance probability distributions for both shape and number density signals is a fair assumption. We add the warning that one may have to account for selection biases, for instance if one investigates cosmic shear and magnification effects with the same galaxy sample (e.g. Schmidt et al. 2009; Krause & Hirata 2009).
We write the contribution of magnification effects to the number density measurement as (13)As several assumptions enter this equation, in particular concerning the treatment of the slope α^{(i)}, we provide a detailed derivation of (13) in Appendix A. Comparing the projection Eqs. (9)–(13) to the general form (7), one can derive all possible cross and autopower spectra in the form of the general Limber Eq. (8). For completeness we have given these Limber equations in Appendix B.
Good models of the threedimensional source power spectra in the Limber equations (see also the right hand column of Table 1) are unknown except for the nonlinear theory matter power spectrum P_{δδ}. The distribution of galaxies is expected to follow the distribution of dark matter, so that the galaxy clustering power spectra should be related to P_{δδ}. However, to date it is unknown how much the galaxy clustering deviates from dark matter clustering, in particular on small scales. This is usually expressed in terms of the galaxy bias b_{g}, which is a function of both angular scale k and redshift or lineofsight distance χ. Hence, one can write (14)where to describe the crosscorrelation between matter and galaxy clustering, we introduced a correlation coefficient r_{g} in the second equality.
The intrinsic alignment power spectra depend on the intricacies of galaxy formation and evolution within their dark matter environment. Again, precise models of the intrinsic alignment have to rely on baryonic physics and are currently not available. For symmetry reasons we parametrise our lack of knowledge about the intrinsic alignment power spectra similarly to the galaxy bias as (15)with the intrinsic alignment bias b_{I} and correlation coefficient r_{I} (following Bernstein 2009). Although the power spectrum P_{gI} could in principle contain a third, independent correlation coefficient, we assume that it is sufficient to write (16)i.e. we hypothesise that correlations between intrinsic number density fluctuations and intrinsic alignments can entirely be traced back to the effects of the intrinsic alignment bias and the galaxy bias. This is a strong assumption since instead of introducing a fifth completely unconstrained bias term, (16) establishes a link between the galaxy bias and intrinsic alignment biases.
Our equation (16) is effectively included within the last term in curly brackets of Bernstein (2009), Eq. (19), and we have effectively set s^{gκ} = 0 in Bernstein (2009), Eq. (35). Bernstein (2009) fixes his s^{gκ} to be a single unknown scalar across the survey, stating that we expect this type of crosscorrelation to have a minimal effect on cosmological constraints. It would be interesting to check this by comparing results in which s^{gκ} is allowed to vary with those in which s^{gκ} = 0. However this is beyond the scope of this paper.
We note that the unknown quantity in question is the crosscorrelation between the intrinsic alignment field and the galaxy position field (r^{gκ} in the Bernstein 2009, notation), which is precisely the correlation measured in the observational constraints papers (e.g. Mandelbaum et al. 2006; Hirata et al. 2007; Mandelbaum et al. 2010). Our ansatz states that this is simply related to the crosscorrelation between the intrinsic alignment field and the mass, and the crosscorrelation between the galaxy position field and the mass (r^{gκ} = r^{g}r^{κ} in the Bernstein 2009, notation). These latter two quantities are much harder to measure because we do not know the mass field well. Indeed gravitational lensing and the observations we discuss in this paper are likely the best probe of these quantities. Ultimately these data could be used to constrain this additional freedom.
To shed further light on this question, we now discuss a simple toy model to illustrate the nature of our assumption. In the simplest case the galaxy distribution may trace the mass distribution, and the intrinsic shear field is the gradient of the mass distribution. In a toy universe we could displace the position mass field on the sky relative to the other two fields which would render the correlation coefficient r_{I} (or r^{κ} in the notation of Bernstein 2009) zero and also the correlation r_{g} = 0 (for large displacements), leaving a strong correlation between the intrinsic shear field and the galaxy position field (unity r^{gκ} in the Bernstein 2009, notation, thus s^{gκ} = 1).
Similar more physical arguments could be made by considering stochasticity instead of displacements. However, we might expect the physical origin of both the galaxy position field and the intrinsic shear field to both lie in the mass field, since most of the interactions are mediated by gravity. Therefore it seems reasonable to expect there to be no additional interplay between the intrinsic shear field and galaxy position fields that would produce further crosscorrelations (nonzero s^{gκ} in the Bernstein 2009, notation).
We will not limit the values of the correlation coefficients to the interval [−1; + 1]. It is formally possible that r > 1 if our assumption about the statistics of the galaxy distribution, usually taken to be Poissonian, is incorrect (Bernstein 2009). Treating the correlation coefficients as completely free parameters, our choice of parametrisation in (14) and (15) is equivalent to modelling P_{gg} and P_{gδ}, or likewise P_{II} and P_{δI}, independently.
We insert the parametrisations (14), (15), and (16) into the set of Limber equations and can this way relate all power spectra entering (4)–(6) to the threedimensional matter power spectrum: The matter power spectrum and the distances f_{K}(χ), which also enter q^{(i)}(χ), depend on cosmology and can therefore be exploited to constrain cosmological parameters. While distances and P_{δδ} are well known from theory, the probability distribution of galaxies p^{(i)}(χ) has to be measured by using additional spectroscopic redshift information (e.g. Huterer et al. 2006; Ma et al. 2006; Abdalla et al. 2007; Bridle & King 2007; Bernstein & Huterer 2010) with a certain level of uncertainty. It may also be possible to infer some additional information from the cosmic shear data itself (Newman 2008; Schneider et al. 2006; Zhang et al. 2010). The same holds for the slopes of the luminosity function α^{(i)}, which can be determined from the survey by studying the flux of galaxies close to the magnitude limit. The least known quantities in the equations above are the bias terms , for which we will thus introduce a very general parametrisation in Sect.3.3.
3. Modelling
In this section we detail the modelling of the terms entering (17) to (26). We specify how we parametrise the uncertainty in the galaxy redshift distributions, the slope of the luminosity function, and the bias terms. Moreover we describe our Fisher matrix approach and the way we infer the resulting errors on cosmological parameters.
3.1. Matter power spectrum and survey characteristics
As the basis for our analysis we compute matter power spectra for a spatially flat CDM universe with fiducial parameters Ω_{m} = 0.25, Ω_{DE} = 0.75, and H_{0} = 100h_{100}kms^{1}Mpc^{1} with h_{100} = 0.7. We incorporate a variable dark energy model by parametrising its equation of state, relating pressure p_{DE} to density ρ_{DE}, as (Chevallier & Polarski 2001; Linder 2003) (27)where the ΛCDM Universe is chosen as the fiducial model, i.e. w_{0} = − 1 and w_{a} = 0. The dark energy density parameter is then given by integrating Eq.(3) of Linder (2003), (28)The primordial power spectrum of matter density fluctuations is assumed to be a power law with fiducial slope n_{s} = 1. We employ the fiducial normalisation σ_{8} = 0.8. The transfer function of Eisenstein & Hu (1998) is used without baryonic wiggles, computing the shape parameter with a fiducial value of Ω_{b} = 0.05.
The nonlinear corrections to the power spectrum are computed by means of the fit formula by Smith et al. (2003). We account for the influence of dark energy on structure growth by modifying the halo model fitting routine of Smith et al. (2003) following the approach of Réfrégier et al. (2008). We provide a summary of this modification in Appendix D.
The survey characteristics follow the rough specifications of a Stage IV experiment (Albrecht et al. 2006) such as the ESA Euclid satellite mission. To compute the noise properties, we assume the maximum extragalactic sky coverage of A = 20000deg^{2} and a total number density of galaxies n = 35 arcmin^{2}. Shape noise is characterised by a total dispersion of intrinsic ellipticities of σ_{ϵ} = 0.35. We refer to this survey as Euclidlike in the remainder of this paper.
According to Smail et al. (1994) we assume an overall number of galaxies per unit redshift, per square arcminute (29)with the galaxy surface density Σ_{0} and β = 1.5. The probability distribution over all the galaxies p_{tot}(z) is proportional to the number density n_{tot}(z). We set , which produces a distribution with median redshift z_{med} = 0.9. The distribution is cut at z_{max} = 3 and then normalised to unity. For the tomography we define photometric bins by dividing the distribution (29) such that every bin contains the same number of galaxies. This choice is merely for computational convenience and to allow for an easy comparison between results with a different number of bins. As default we will use N_{zbin} = 10 bins.
To account for photometric redshift errors, we assume that the fraction of catastrophic failures in the assignment of photometric redshifts is negligible, but include the spread of the true redshifts in the binwise distributions by writing the conditional probability of obtaining a photometric redshift z_{ph} given the true redshift z as (30)where σ_{ph} denotes the photometric redshift dispersion. The redshift distribution of an individual photoz bin p^{(i)}(z) is then obtained by integrating (30) over the bin width and by weighting the result by the overall redshift distribution (29), see Joachimi & Schneider (2009) for details. We use σ_{ph} = 0.05 as our default value. This also follows fiducial Model 1 of Ma et al. (2006).
Since the underlying redshift distributions p^{(i)}(z) are determined by measurement, they are not perfectly known, but introduce further uncertainty into the analysis. A detailed analysis of the dependence of the joint analysis of galaxy shape and number density information on redshift parameters, and also the potential of calibrating these errors internally, will be investigated elsewhere (e.g. Zhang et al. 2010). Bridle & King (2007) have undertaken a more detailed study of the effect of redshift errors in the case of ellipticity correlations only. For the purpose of this work we assume that the value of σ_{ph} is unknown, i.e. we use it as a single, global parameter to account for the uncertainty in the redshift distributions. We employ a wide Gaussian prior on σ_{ph} of 10 for reasons of numerical stability.
3.2. Galaxy luminosity function
In order to calculate power spectra which include the lensing magnification signal, we need to model the slope of the cumulative galaxy luminosity function at the magnitude limit of the galaxy number density catalogue. In Appendix C we extend observational results for the normalisation and redshift scaling of the galaxy redshift distribution (29) by Blake & Bridle (2005) to provide a fitting formula for the luminosity function slope as a function of redshift and survey magnitude limit.
We use the fit given by (C.2) and (C.3) with the parameters listed in Table C.1 to compute the slope of the luminosity function at r_{lim} = 24. The discussion in this work applies to groundbased surveys because the COMBO17 luminosity functions are calculated for the SDSS r filter as observed from the ground. A space mission to a depth of r_{lim} = 24 will have a different luminosity function slope, corresponding more closely to a deeper ground based survey, depending on the resolution of the spacebased survey. We use results for r_{lim} = 24 throughout this paper for both ground and space surveys. We note that from Fig.C.1, top panel, the slope of the luminosity function is changed little on increasing the survey depth beyond r_{lim} = 24.
The fiducial slope in a photoz bin i is defined as , where is the median redshift of bin i, see Appendix A. We assume α^{(i)} is also measured from the survey itself, and therefore adds another source of uncertainty to the analysis which we account for by setting α^{(i)} to be free parameters for all i = 1,...,N_{zbin}. Again we apply a wide Gaussian prior of 10 on every slope parameter, which does not have a measurable influence on our results, but merely ensures numerical stability.
3.3. Galaxy and intrinsic alignment bias
As already outlined in Sect.2, the bias terms encoding the galaxy bias and intrinsic alignments are the least accurately known contributions to (17) to (26). We parametrise each of these terms on a grid in k and z, following Bridle & King (2007) whose ansatz is in turn similar to the recommendations by the Dark Energy Task Force (Albrecht et al. 2006) and Bernstein (2009). Every bias term is assumed to vary around a fiducial functional form as (31)where Q_{X}(k,z(χ)) is an unknown twodimensional function which comprises the aforementioned grid, and where A_{X} denotes an additional free overall amplitude. We use N_{K} bins in k and N_{Z} bins in redshift for each bias term and linearly interpolate in the logarithms of Q_{X}, k, and 1 + z, so that Q_{X} is given by (32)for k_{i} < k ≤ k_{i + 1} and z_{j} < z ≤ z_{j + 1}, where we defined (33)The free parameters are the grid nodes . Since Q_{X} is a multiplicative function, one reproduces the base model if A_{X} = 1 and if all vanish. The effect of this parametrisation on the observable projected power spectra is illustrated in Bridle & King (2007) for the case of intrinsic alignments.
The indices in (32) run from i = 0,...,N_{K} + 1 and 0 = 1,...,N_{Z} + 1. We fix all parameters at the edge of the grid by setting the parameters with indices i,j = 0; i = N_{K} + 1 or j = N_{Z} + 1 to , so that we have N_{K} × N_{Z} free grid parameters per bias term. We place the lowest and highest grid nodes at the limits of our integration ranges, so k_{0} = 3.3 × 10^{7} h Mpc^{1} and k_{NK + 1} = 3.3 × 10^{4} h Mpc^{1} in k, and in the redshift dimension z_{0} = 0 and z_{NZ + 1} = 19. The grid nodes, which are free to vary, are loglinearly spaced in a smaller range, respectively. We use k_{1} = 10^{3} h Mpc^{1} and k_{NK} = 2 h Mpc^{1}, and for the redshift range z_{1} = z_{0} and z_{NZ} = 3. In the special case of N_{K} = 1 we position the only free parameter in the k dimension at the centre between k_{1} and k_{NK}, and proceed likewise for redshifts.
Overview on the total number of nuisance parameters used for different setups.
It is important to note that while Bridle & King (2007) limit the flexible grid parametrisation to the nonlinear regime of the power spectra, we attempt to cover all k ranges which substantially contribute to the observable power spectra. As we fix the grid values on the edges, the overall scaling of the bias terms is not free, so that we use the amplitude A_{X} as a further varying parameter throughout. To all bias term parameters we add a very wide Gaussian prior of standard deviation 50 to ensure numerical stability. Together with the global uncertainty on the redshift distributions, expressed in terms of σ_{ph}, and the values of the slope of the galaxy luminosity function per photoz bin, we obtain a large number of nuisance parameters that we determine simultaneously with the cosmological parameters of interest. For later reference, we have summarised the total number of nuisance parameters for different setups in Table 2. While our parametrisation is fairly general and should capture most of the variability, it is of course possible that the bias terms depend on more parameters than k and z. For instance, it is well known that both intrinsic alignments and galaxy bias are a function of galaxy colour and luminosity which could be incorporated into our approach in the future. For observational constraints on this effect using intrinsic alignments see Mandelbaum et al. (2006), Hirata et al. (2007) and in galaxy biasing see McCracken et al. (2008), Swanson et al. (2008), Simon et al. (2009), Cresswell & Percival (2009), Wang et al. (2007) for recent examples.
Fig. 1 Fiducial power spectra for all considered correlations. The upper right panels depict the contributions to ϵϵ (in black) and nn (in magenta) correlations. The lower left panels show the contributions to correlations between number density fluctuations and ellipticity. Since we only show correlations with i ≤ j, we make in this plot a distinction between nϵ (in red; number density contribution in the foreground, e.g. gG) and ϵn (in blue; number density contribution in the background, e.g. Gg) correlations. In each subpanel a different tomographic redshift bin correlation is shown. For clarity only odd bins are displayed. In the upper right panels the usual cosmic shear signal (GG) is shown as a black solid lines; the intrinsic alignment GI term is shown by the black dashed lines; the intrinsic alignment II term is shown by the dotted black line; the usual galaxy clustering signal (gg) is shown by the magenta solid line; the cross correlation between galaxy clustering and lensing magnification (gm) is shown by the magenta dashed line; the lensing magnification correlation functions (mm) are shown by the magenta dotted line. In the lower left panels the solid blue line shows the correlation between lensing shear and galaxy clustering (Gg); the blue dashed line shows the correlation between lensing shear and lensing magnification (gm); the blue dotdashed line shows the correlation between intrinsic alignment and galaxy clustering (Ig or equivalently gI); the red solid line shows the correlation between galaxy clustering and lensing shear (gG), which is equivalent to the blue solid line with redshift bin indices i and j reversed; similarly the red dashed line shows the correlation between lensing magnification and lensing shear (mG), for cases where the magnification occurs at lower redshift than the shear (i < j); finally the dotted line shows the correlation between lensing magnification and intrinsic alignment (mI). 
To compute the fiducial models for all the power spectra entering the observables (4) to (6), we set A_{X} = 1 and all , i.e. they are fully determined by the base models. We set (34)With the choice for and we reproduce the nonlinear modification of the linear alignment model by Bridle & King (2007). Lacking solid physical motivation, it is yet in agreement with current observational evidence (Mandelbaum et al. 2006; see Bridle & King 2007, for a comparison) and the halo model studies by Schneider & Bridle (2010). It is based on the linear alignment model (Hirata & Seljak 2004) which is expected to provide a good description of intrinsic alignments on the largest scales. We assume the galaxy bias to be of order unity for our fiducial model, and set . Note that most investigations of galaxy clustering consider much less flexibility in the galaxy bias.
In Fig.1 we plot the fiducial angular power spectra of all considered signals for different combinations of photoz bins. Note that the ellipticitynumber density crosscorrelations are not symmetric under exchanging the photometric redshift bins. Hence, in this figure we treat nϵ and ϵn correlations, as well as all signals contributing to them, separately, keeping i ≤ j for all .
The GG signal shows the usual behaviour of moderate increase with increasing redshift of the contributing photozbins. The redshift scaling of the IG term is similar, but peaks when the source galaxies are at high redshift in the background (large j), while the galaxies that are intrinsically aligned are at low redshift (small i). For the model used here the IG contribution can even surpass the shear signal in this case. Due to the narrow kernel containing p^{(i)}(χ)p^{(j)}(χ), see (19), the II signal is strong in the autocorrelations i = j, but drops off quickly as soon as the overlap of the redshift distributions decreases.
Due to the similar kernel, the scaling of the galaxy clustering contribution (gg) resembles the II term, but gg constitutes a much stronger signal. Lensing magnification (mm) adds the largest fraction of the galaxy number correlations at the highest redshifts, showing a slightly stronger redshift scaling than GG. However, the mm term always remains subdominant with respect to signals with a contribution from galaxy clustering; even for widely separated galaxy redshift distributions, say i = 1 and j = 9, and the gm crossterm is considerably stronger than mm. Such contributions might be a serious obstacle for probing cosmology with the lensing magnification signal as proposed by Broadhurst et al. (1995), Zhang & Pen (2005, 2006), van Waerbeke (2010). Yet in our approach, where the galaxy bias is taken into account and parametrised, the magnification signal yields a valuable contribution to the galaxy number correlations, which helps constraining the cosmological model.
The signals within the number densityellipticity crosscorrelations are not symmetric when swapping the photoz bins. When the contribution by number density fluctuations stems from the foreground, the gG signal is strong, in particular if the photoz bins are far apart in redshift, whereas the “Gg” (in the notation of Fig.1) drops off fast if i < j because the shear signal of foreground galaxies is not correlated with the clustering of galaxies at much higher redshift. The mG, Gm, and GG signals differ only by the term including the slope of the luminosity function and thus have similar amplitudes. Correlations between intrinsic alignment and magnification (mI and Im) are subdominant throughout, obtaining their largest amplitudes if intrinsic alignments at low redshifts are combined with the magnification signal from galaxies far in the background, i.e. in the figure for Im at i = 1 and j = 9. Finally, the symmetric Ig term is the largest contribution for autocorrelations of number densityellipticity observables, but decreases quickly in the cross terms, again due to the kernel p^{(i)}(χ)p^{(j)}(χ) since we have assumed the photometric redshift errors are reasonably well behaved, without catastrophic outliers. Note that we have plotted the absolute values of the power spectra in Fig.1 and that the correlations IG, gm, Gm, mG, and gI are negative.
3.4. Parameter constraints
We determine constraints on our parameters using a Fisher matrix analysis. To account for the errors and correlations of our observables, we compute covariances of the power spectra (4) to (6) in the Gaussian approximation, extending the results of Joachimi et al. (2008), see also Hu & Jain (2004). If we denote the difference between estimator and its ensemble average by , one can write for the covariance (35)where A is the survey size and Δℓ the width of the corresponding angular frequency bin. As the Kronecker symbol δ_{ℓℓ′} indicates, the covariance is diagonal in ℓ in the Gaussian limit, which keeps the computation and inversion of (35) tractable. The subscripts {α,β,γ,δ} can be either ϵ or n, where holds. To account for the shot and shape noise contributions induced by the random terms in (1) and (2), we have defined (36)the second term given by (37)Here denotes the total intrinsic ellipticity dispersion, and is the average galaxy number density per steradian in photoz bin i.
Fig. 2 Left panels: figures of merit as a function of the number of free parameters as a function of wave vector N_{K} in the bias terms For each line type, the upper curve is obtained for a number of free bias parameters as a function of redshift N_{Z} = 2, the lower is for N_{Z} = 4. Right panels: same as on the left, but as a function of N_{Z}, i.e. the number of redshift parameters in the bias terms. The upper curves for each set correspond now to N_{K} = 2 and the lower ones to N_{K} = 4, respectively. Upper panels: figure of merit taking into account the full cosmological parameter space, FoM_{TOT}, see (45). Lower panels: dark energy figure of merit from the Dark Energy Task Force FoM_{DETF}, see (44). Dashed curves correspond to results using galaxy ellipticity correlations (ϵϵ) only, dotted black curves to galaxy number density correlations (nn) only, and solid black curves to results using all correlations (ϵϵ, nn and ϵn). The grey dotted lines show results for nn correlations without imposing cuts in angular frequency. The constant grey line marks the FoM computed for the pure lensing, i.e. GG, signal, assuming intrinsic alignments do not exist. In addition we show the resulting figures of merit when using our most flexible parametrisation with N_{K} = N_{Z} = 7 as filled symbols. Circles correspond to ϵϵ, triangles to nn, and diamonds to all correlations. 
Combining the observable power spectra, we compose the total data vector (38)for every angular frequency considered. The corresponding covariance, again for every ℓ, reads (39)with the block matrices given by (35). The number of galaxy ellipticity (ϵϵ) and number density (nn) observables entering is N_{zbin}(N_{zbin} + 1)/2, respectively, while there are ellipticitynumber density cross terms (nϵ), which are not symmetric. In the analysis that follows we will also consider ϵϵ and nn correlations only. In these cases the covariance (39) is reduced accordingly to its upper left or lower right block.
For reasons of computational time the total number of parameters that we can consider is limited to a few hundred. As a consequence the k dependence of the galaxy bias can not be parametrised by more than about ten parameters per redshift grid node. This number might not provide enough freedom in b_{g} (and r_{g}) to represent a sufficiently general set of functional forms, which inadvertently may cause strong constraints on cosmological parameters due to the strong signal of galaxy clustering. Hence, we follow existing studies of galaxy clustering by discarding the clustering contribution in the nonlinear regime where the signal is largest and the form of the galaxy bias most uncertain.
Rassat et al. (2008) calculated wave vectors as a function of redshift at which the threedimensional power spectrum has to be cut off to avoid number density fluctuations above a certain threshold, used as an indicator for nonlinearity. Since we do still have a fairly general parametrisation of the bias terms associated with galaxy bias, we can afford to include the mildly nonlinear regime into our analysis. Consulting Fig.2 of Rassat et al. (2008), we choose a simple linear parametrisation of the form (40)This relation roughly coincides with the fiducial curve in the figure, producing slightly more conservative cuts at low redshifts.
We do not cut the threedimensional power spectrum in kspace, but instead exclude projected power spectra above a threshold angular frequency from the Fisher matrix analysis. This maximum angular frequency is computed via (41)where we choose as a characteristic redshift of bin i the median redshift . Hence, we obtain a cutoff ℓ for every photoz bin. We choose that ϵϵ correlations are not at all affected by this cutoff because they are not dominated by terms involving galaxy bias. We impose on nϵ correlations, where i is the photoz bin from which the number density signal stems. For observables we use the cutoff calculated for bin j. Note that, due to the fast dropoff of the galaxy clustering signal with increasingly different median redshifts of bins i and j, the more optimistic choice of bin j over i in the latter case should not influence our results decisively.
Overview on the cuts in angular frequency.
To compute the Fisher matrix, we use N_{ℓ} = 50 angular frequency bins, spaced logarithmically between ℓ_{min} = 10 and ℓ_{max} = 3000, the latter value being a conservative maximum for future surveys. We assume that the covariance of the power spectra is independent of the cosmological parameters, so that it does not contribute to the constraints. Then the Fisher matrix reads (e.g. Tegmark et al. 1997) (42)where N_{d} is the dimension of . The first summation in (42) runs over all usable angular frequency bins. The number of usable bins depends on the type of correlation and is determined by the cutoff angular frequency as described above. For the default setup we have summarised , , and the number of usable bins for every photoz bin in Table 3. The derivatives in (42) are taken with respect to the elements of the parameter vector (43)The first seven entries of p correspond to the cosmological parameters that we are interested in, while the remaining nuisance parameters account for the uncertainty in the galaxy redshift distribution, the slope of the galaxy luminosity function, the intrinsic alignments, and the galaxy bias, see Table 2. In summary, we use a maximum of observables (actually significantly less due to the ℓcuts of the galaxy number density signals) to measure a total of 4N_{Z}N_{K} + N_{zbin} + 8 parameters. Note that since we have referred all signals contributing to the observables (4) to (6) to the matter power spectrum, they all constrain the set of cosmological parameters; none of them is fixed when calculating (42).
The minimum variance bound of the error on a parameter p_{μ}, if determined simultaneously with all other parameters, is given by . This error provides us with a lower bound on the marginalised 1σerror on p_{μ}. To assess the statistical power of the survey by means of a single number, we use the figure of merit (FoM) suggested by the Dark Energy Task Force (DETF) Report (Albrecht et al. 2006), (44)where the subscript “w_{0}w_{a}” denotes the 2 × 2 submatrix of the inverse Fisher matrix that corresponds to the entries belonging to the two dark energy parameters. Note that different prefactors for (44) are used in the literature. To allow for direct comparison with Bridle & King (2007), divide our findings for the FoM_{DETF} by four.
While (44) is restricted to the quality of constraints on dark energy, we also seek to consider the errors on all cosmological parameters of interest. We are interested in the total volume of the error ellipsoid in parameter space, which is measured by the determinant of the Fisher matrix. Hence, we define (45)where only the submatrix of the inverse Fisher matrix that corresponds to the seven cosmological parameters under investigation is used in the determinant, as indicated by the subscript. The determinant of the inverse is computed in order to take the effect of marginalising over nuisance parameters into account.
4. Results
Based on the Fisher matrix formalism described in the foregoing section, we will now analyse the performance of a cosmological galaxy survey with combined number density and shear information. We are going to investigate the residual information content in the data after marginalising over models of the intrinsic alignments and the galaxy bias with varying degrees of freedom. Furthermore, we will study the dependence of our FoM on the number of photoz bins and the width of the binwise redshift distributions as well as on the priors imposed on the nuisance parameters. The information contained in the individual signals and their susceptibility to the nuisance parameters is also assessed. Throughout this section we use the default survey characteristics and parameter values unless specified otherwise.
4.1. Dependence on intrinsic alignments and galaxy bias
The four bias terms each comprise N_{K} × N_{Z} + 1 nuisance parameters. The galaxy ellipticity (ϵϵ) power spectra only contain the two intrinsic alignment bias terms, whereas the number density correlations (nn) are only affected by the galaxy bias terms. The crosscorrelations between ellipticity and number density link those signals which depend on both galaxy bias and intrinsic alignments and thus allow for their internal crosscalibration. An example is the study by Zhang (2010) which investigates the interrelations between the IG, gI, and gg terms.
In Fig.2 both figures of merit are plotted as a function of N_{K} and N_{Z}, respectively. If we restrict our analysis to ϵϵ correlations only, our setup is similar to the most realistic setup considered in Bridle & King (2007). We reproduce their result that the FoM_{DETF} drops as a function of N_{K}, dropping most sharply at small N_{K} and then levelling off. It falls significantly below the reference value of the FoM_{DETF} computed for a pure GG signal. However, while Bridle & King (2007) find a decrease by about a factor of 2 for N_{K} = 10 and N_{Z} = 2 compared to the lensing only case, our FoM_{DETF} decreases by more than a factor of 4. This discrepancy can be traced back to the fact that Bridle & King (2007) limit the nuisance parameter gridding to the nonlinear regime in kspace. We additionally plot the figures of merit as a function of N_{Z}, finding that the ϵϵ results do in fact flatten as N_{Z} reaches high values around 10. We find a very similar behaviour in terms of the FoM_{TOT}. The dark energy parameters suffer more than other parameters from the uncertainty of the bias terms in the redshift direction, as the FoM_{DETF} deteriorates faster than the FoM_{TOT} as a function of N_{Z}.
With our default settings the pure nn correlations constrain cosmology much more weakly than ellipticity correlations. Recall that our galaxy clustering analysis uses a much more flexible bias parametrisation than most other work. Using a small number of nuisance parameters we get a FoM_{DETF} which is marginally above unity. This result is of the same order of magnitude as the findings of Rassat et al. (2008) who determine FoM_{DETF} = 4.2 for a spectroscopic, spacebased survey in the spherical harmonics approach. Contrary to us, Rassat et al. (2008) use a less flexible bias parametrisation and include BAO features, but neglect magnification effects in their study. The decrease in both figures of merit for a larger number of nuisance parameters is more pronounced than for ϵϵ correlations, in particular as a function of N_{Z}.
We have also shown the results obtained without the additional angular frequency cuts in Fig.2 (i.e. the maximum angular frequency is ℓ_{max} = 3000 for all angular power spectra). Because the full galaxy clustering signal is strong and a comparatively direct probe of the dark matter power spectrum if the galaxy bias is well known, the values of the figures of merit are much higher in this case and even surpass the lensingonly level if few nuisance parameters are used. This corresponds to the case where we have a reasonably good understanding of galaxy biasing over a large range of scales, into the deeply nonlinear regime. For larger N_{K} and N_{Z} the decrease in FoM is considerable, but weaker than with ℓcuts.
In the limit of a large number of nuisance parameters we expect that the galaxy bias is parametrised in a sufficiently flexible way, such that the curves with and without ℓcuts should coincide or at least be of the same order of magnitude. Since this is not the case, and since there is no obvious sign of a lower boundary that the figures of merit are approaching, we hypothesise that this limit is only achieved for a very large, numerically and computationally prohibitive, number of nuisance parameters. Thus, the use of the additional angular frequency cuts in our default analysis is the most practical way to take into account our lack of understanding of galaxy biasing on small scales.
The simultaneous use of all available observables significantly boosts the parameter constraints, due to the addition of extra information from the survey, and due to breaking degeneracies between cosmological and nuisance parameters. The FoM_{DETF} is up to a factor of about 50 higher than the lensingonly value. Both figures of merit decrease for larger N_{K} and N_{Z}, attaining mostly shallow slopes at N_{K},N_{Z} ≈ 10. Considering the most flexible model in each panel of Fig.2, the FoM_{TOT} remains above, but close to the value for a pure GG signal while the FoM_{DETF} degrades slightly below the lensingonly value. In contrast to the setup with ϵϵ or nn correlations only, the decrease is less pronounced with N_{Z} than N_{K} here, i.e. the crosscalibration between all correlations can partially compensate the loss in redshift information that affected the dark energy parameters in the former cases.
Marginalised parameter errors and the two figures of Merit for a survey with ten photometric redshift bins used for tomography N_{zbin} = 10 and a photometric redshift uncertainty parameter σ_{ph} = 0.05.
The most flexible configuration we consider has N_{K} = N_{Z} = 7 which corresponds to 200 nuisance parameters within the four bias terms. This limit is not inherent to our method, but is merely set for computational practicality. Both the intrinsic alignments and the galaxy bias have a physical origin and therefore are expected to produce smooth signals, which should not oscillate strongly or feature sharp peaks. Hence on the default angular scales used this model with 100 free parameters for each of intrinsic alignments and galaxy bias should yield a fairly general representation of the signals if one can rely on coarse prior information on the fiducial base model.
Fig. 3 1σcontours for all pairs of cosmological parameters considered, marginalised over all other parameters. We have used a photometric redshift uncertainty parameter value σ_{ph} = 0.05, ten photometric redshift bins for tomography N_{zbin} = 10, and the most flexible intrinsic alignment and bias model considered in this paper, with over two hundred free parameters (N_{K} = 7, and N_{Z} = 7). Orange (light hatched) confidence regions result from using galaxy number density correlations (nn) (excluding the nonlinear regime) only, red (dark hatched) regions use ellipticity correlations (ϵϵ) alone, and blue (filled) regions correspond to using all available information including densityellipticity crosscorrelations. For reference, the contours obtained from a pure lensing signal are shown as black lines. Flat priors on cosmological parameters have been applied. 
The figures of merit and the individual parameter errors for the most flexible model are given in Table 4, for the pure GG signal, ϵϵ, nn, and all correlations. Compared to lensing alone, the FoM_{DETF} decreases by about a factor of 20 for ϵϵ correlations, which means that cosmic shear is severely affected if one assumes very little prior knowledge about intrinsic alignments. There are hardly any dark energy constraints for nn correlations only, given this freedom in the galaxy bias. However, using all available correlation simultaneously, we can recover just over half of the pure cosmic shear DETF figure of merit. We find that the total error volume of the cosmological parameters, i.e. the FoM_{TOT}, is the same using all the correlation information as for the lensingonly case. Looking at the marginalised parameter errors, w_{0} and w_{a} suffer particularly strongly in the ϵϵ only and nn only case. Less information about the dark energy parameters is lost if all correlations are used, and the constraints on h_{100} and Ω_{b} improve over the lensingonly case. In Fig.3 we show in addition the marginalised 1σerror ellipses of all possible pairs of cosmological parameters for the same setup. The degeneracy directions are fairly similar for the usual lensing only case and the results with all correlations presented in this paper.
We repeat this analysis for a nearer term StageIIIlike survey such as the Dark Energy Survey (DES). To this end, we use a survey size of A = 5000 deg^{2}, a median redshift of the overall redshift distribution of z_{med} = 0.8, a total number density of galaxies of n = 10 arcmin^{2}, and a total ellipticity dispersion of σ_{ϵ} = 0.23, keeping all other parameters at their default values. We refer to this survey as DESlike for the rest of the paper.
Note that different groups use different conventions for quoting the number density of galaxies and the ellipticity dispersion. The Euclid team define σ_{ϵ} = 0.35, which is motivated by considering it to be the typical uncertainty on the shear taking into account both shape noise and measurement errors. The effective number density has then been estimated relative to this value and found to be between 30 and 40 arcmin^{2} (Laureijs et al. 2009). Because of the choice of σ_{ϵ}, this number density value is approximately equal to the actual number of usable galaxies in the image. The DES team start by defining σ_{ϵ} = 0.23, motivated by considering the shape noise alone. The effective number density is then found to be n = 10 arcmin^{2} (Annis et al. 2005). Due to this choice of σ_{ϵ}, the number density value is smaller than the number density of usable galaxies in the image, and represents the effective number density of galaxies with negligible measurement errors. The different definitions of σ_{ϵ} also explain the low number density quoted for DES in comparison with shallower surveys such as the CHFTLS.
In this work we choose to take the number density and ellipticity dispersion values as provided by the collaborations and do not adjust them, which would be beyond the scope of this paper. In any case the physically meaningful quantity is the ratio σ_{ϵ}/ which determines the uncertainties on the power spectra, see (37). Thus, as long as we use consistent definitions of n and σ_{ϵ}, we obtain the correct shape noise contributions to the cosmic shear covariance.
Our findings are also shown in Table 4. For the DESlike survey the FoM_{TOT} for all correlations is again about the same as in the lensingonly case, and the ratios of FoM_{DETF} values are slightly smaller than for the StageIVlike survey. The FoM_{DETF} values for ϵϵ, nn, and all possible correlations are about a factor of ten larger from the StageIVlike cosmic shear survey than the StageIIIlike survey, offering significant benefit beyond the minimal requirement of Albrecht et al. (2006).
4.2. Dependence on characteristics of the redshift distribution
It is well known that dividing the galaxy sample into several redshift photoz bins greatly improves constraints from cosmic shear, but due to the broad lensing kernel (10) there is little benefit in having more distributions than three to five (e.g. Hu 1999; Simon et al. 2004; Ma et al. 2006). This result does not hold true anymore if one aims at controlling the intrinsic alignment contamination in the cosmic shear signal, mainly manifest via its characteristic redshift dependence. Using both marginalisation (Bridle & King 2007) and parameterfree approaches (Joachimi & Schneider 2008), one finds that the figures of merit only start to stabilise when using ten redshift distributions or more.
Fig. 4 Upper panel: the figure of Merit for all cosmological parameter space FoM_{TOT} as a function of the number of photoz bins used for tomography N_{zbin}, shown for ϵϵ (dashed line), nn (dotted line), and all (solid) correlations. The grey line corresponds to results for lensing only (GG). Throughout, N_{K} = N_{Z} = 5 nuisance parameters for the bias terms are used. These results are obtained for the standard set of parameters and σ_{ph} = 0.05. Lower panel: same as above, but in terms of the dark energy figure of merit FoM_{DETF}. 
In Fig.4 our figures of merit are plotted as a function of the number of photoz bins, using N_{K} = N_{Z} = 5 nuisance parameters per bias term. In agreement with the findings mentioned above both the FoM_{TOT} and the FoM_{DETF} in the case of a pure cosmic shear signal become approximately constant for N_{zbin} ≳ 5, and increasing N_{zbin} beyond three has little effect. For ϵϵ correlations the gain in FoM is considerable up to N_{zbin} ~ 7; for larger N_{zbin} the curves rise only slowly.
Since nn correlations are more localised because of the compact kernel of the dominating galaxy clustering (gg) signal, increasing N_{zbin} proves beneficial for these observables up to the maximum number of photoz bins we have considered, although the bins feature an increasing overlap of their corresponding redshift distributions as we keep σ_{ph} = 0.05 fixed. As one would expect, we obtain an intermediate scaling with N_{zbin} for the complete set of available correlations. Our fiducial choice of N_{zbin} = 10 is beyond the regime of strongly varying figures of merit at small N_{zbin}, but the further increase in FoM is more pronounced than for a pure lensing signal or ϵϵ correlations only with FoM_{DETF} rising by an additional 80% on increasing the number of photometric redshift bins from 10 to 20.
Figure 5 shows the figures of merit as a function of the photoz dispersion, normalised to the value at σ_{ph} = 0.02. Since the FoM_{TOT} is a logarithmic quantity, we compute differences rather than ratios of the figure of merit, i.e. Hence the difference d_{FoM} is directly related to the change in volume of the error ellipsoid spanned by the set of cosmological parameters. We have returned to our default value of N_{zbin} = 10.
The pure lensing signal needs merely coarse redshift information to attain its full statistical power, and hence its FoM hardly suffers from the increasing spread in the redshift distributions. However, as redshift information is vital to account for intrinsic alignments, the figures of merit for the cosmic shear signal after marginalisation over intrinsic alignments decrease considerably by more than 80% in terms of the FoM_{DETF} on changing from σ_{ph} = 0.02 to σ_{ph} = 0.1, which is in line with Bridle & King (2007). For the same change in σ_{ph} the FoM_{TOT} is reduced by about 12. It is interesting to note that the dependence of the FoM_{TOT} on σ_{ph} is close to linear for both lensing and ϵϵ signal.
Fig. 5 Upper panel: the difference d_{FoM}, defined in (47), as a function of the photoz dispersion σ_{ph}, shown for ϵϵ correlations (dashed line), all correlations (black solid line), and the lensing only signal (grey solid line). Throughout, nuisance parameters N_{K} = N_{Z} = 5 are used. These results are obtained for the standard set of parameters and N_{zbin} = 10. Lower panel: Same as above, but in terms of the ratio r_{FoM}, given in (46). 
Adding galaxy number density information largely alleviates the information loss. The degradation in FoM is only slightly stronger than for the pure lensing signals in the expected regime of highquality photoz information with σ_{ph} ≲ 0.06. For larger σ_{ph} the increasing overlap of the redshift distributions causes both figures of merit to decrease further. We have to add the caveat that in this investigation we consider only a single parameter that accounts for the uncertainty in the shape of the redshift distributions, σ_{ph}.
It remains to be seen how sensitive these conclusions are to greater levels of uncertainty in the photometric redshift distribution calibration. For example Zhang et al. (2010) consider the use of nn and ϵn terms for selfcalibrating these uncertainties but do not include intrinsic alignments in their main calculations. If we assumed e.g. uncertain median redshifts of each individual distribution, which is beyond the scope of this paper and under investigation elsewhere, we would have obtained significantly lower figures of merit; see Kitching et al. (2008b) and also Bridle & King (2007) who consider this for the case of ϵϵ correlations only. The curve for all correlations in Fig.5 may hence approach the one for ϵϵ correlations, especially for σ_{ph} ≲ 0.06. Consequently, redshift distributions with a small spread may turn out to be even more desirable when taking galaxy number density information into account.
4.3. Dependence on nuisance parameter priors
On all our nuisance parameters we can expect to have prior information to a certain extent, at least by the time large spacebased cosmological surveys will be undertaken. Here we investigate the dependence of the resulting FoM_{TOT} on tightening the priors on the different sets of nuisance parameters, employing our most general configuration with N_{K} = N_{Z} = 7. We show , i.e. the FoM_{TOT}, referred to its value for the widest priors we apply as default, as a function of the Gaussian prior width σ_{prior} in Fig.6.
Fig. 6 The difference d_{FOM} between the FoM_{TOT} for a given prior value σ_{prior} and the fiducial FoM_{TOT}, obtained for the least stringent, default values of the different priors we apply. Upper panel: effect of tightening the priors on the nuisance parameters in the bias terms for ϵϵ (dashed line), nn (dotted line), and all (solid) correlations. Lower panel: effect of tightening the priors on the slopes of the luminosity function α^{(i)}, and for the prior on the uncertainty of σ_{ph}. Grey curves correspond to priors on α^{(i)}, black curves to those on σ_{ph}. Except for the ϵϵ case the latter remain very close to zero. As above, ϵϵ correlations are shown as dashed lines, nn correlations as dotted lines, and all correlations as solid lines. Note that ϵϵ correlations do not depend on the α^{(i)}. For all curves, the remaining priors are each set to their default values of 50 for the bias term parameters, and 10 for priors on both α^{(i)} and σ_{ph}. 
For reasons of numerical stability, we have imposed wide Gaussian priors on all nuisance parameters, i.e. 10 on both the α^{(i)} and σ_{ph}, and 50 on the parameters within the intrinsic alignment and galaxy bias terms. As the curves clearly indicate, these priors are noninformative since all of the curves have flattened off by this value. The strongest effect is achieved by tightening the priors on the bias terms (note the different scaling of the two panels). Prior knowledge on the galaxy bias tremendously improves constraints by nn correlations, raising the FoM_{TOT} close to the values of the lensingonly case if the uncertainty of the bias term nuisance parameters is of the order 10^{3}. Since in a situation with such excellent prior knowledge about galaxy bias, the cuts in angular frequency could be much less stringent, the increase in FoM would be even more pronounced.
The FoM_{TOT} for ϵϵ correlations changes from its value at the default configuration to the lensingonly value^{6} (compare to Table 4) in the limited range between σ_{prior} ~1 and σ_{prior} ~10^{2}. Thus coarse prior information about intrinsic alignments has little effect on the FoM, but, compared to the galaxy bias, a model with comparatively moderate precision can already bring the FoM_{TOT} back close to its optimum. If all available correlations are used, the internal calibration works well to constrain the nuisance parameters, so that the priors only mildly improve the FoM_{TOT}. However, bear in mind that in spite of our flexible parametrisation, the choice of a fiducial base model for both intrinsic alignments and galaxy bias influences the results. Hence the base models should be as realistic as possible. It is questionable whether our knowledge about the bias terms will ever suffice to impose priors on the bias term nuisance parameters beyond the approximately correct form of the base models.
As already discussed in the foregoing section, the global parameter σ_{ph} is excellently constrained by the nn correlations. Therefore even tight priors do not have any effect on the FoM_{TOT} if nn correlations form part of the data vector. Priors below σ_{prior} = 10^{2} increase the information content in ϵϵ correlations marginally. As far as the uncertainty in the redshift distributions is concerned, our investigation is still to be regarded as idealistic. We defer the joint analysis of shape and number density correlations in presence of unknown parameters in each individual redshift distribution to future work.
The lack of knowledge in the slopes of the luminosity function α^{(i)} entering the lensing magnification signal is negligible in comparison with the effect of the galaxy bias nuisance parameters. If all correlations are considered, the α^{(i)} are well constrained since the improvement in FoM due to priors is below the 1% level.
Fig. 7 The FoM_{TOT} for different subsets of correlations that form the nn signal. Throughout, we used N_{zbin} = 10 and σ_{ph} = 0.05. The marginalisation has been performed for N_{K} = N_{Z} = 7 nuisance parameters (black dashed bars) or N_{K} = N_{Z} = 3 nuisance parameters (grey solid bars). The two leftmost columns show the full set of ϵϵ and nn correlations for reference. The two centre columns stand for gg and mm correlations only, employing the cuts in angular frequency. The rightmost columns show again gg and mm, but here the full information up to ℓ_{max} = 3000 is considered. Note that mm correlations depend neither on intrinsic alignments nor on galaxy bias and hence are independent of N_{K} and N_{Z}. 
4.4. Information content in the individual signals
The question arises of which signals contributing to the observable power spectra contain most of the information about cosmology or suffer most strongly from the uncertainty in intrinsic alignments and galaxy bias. Bridle & King (2007) have studied the effect of the two intrinsic alignment terms on ϵϵ correlations with a parametrisation very similar to ours. In the following we consider nn correlations by repeating the Fisher matrix analysis twice: once assuming that there is no lensing magnification signal, i.e. only gg contributes to nn; and once assuming that there is no intrinsic galaxy clustering, i.e. only mm contributes to nn.
For this calculation we consider constraints from galaxy clustering information alone (nn), and assume that galaxy shape information (ϵϵ and nϵ) is not used. We show the resulting FoM_{TOT} for the most flexible model (N_{K} = N_{Z} = 7) and a model with somewhat less flexibility (N_{K} = N_{Z} = 3) in Fig.7. On comparing the second, third and fourth columns from the left we see that the model with less freedom (solid grey lines) has approximately the same FoM_{TOT} for the complete nn signal as when the nn signal is made up of either one of the pure gg and mm terms alone. We now examine how this conclusion is changed when increased flexibility is allowed in the galaxy bias model (dashed lines). Lensing magnification does not depend on any of our bias terms, so that its FoM_{TOT} value is the same for the flexible model (the dashed and solid lines are on top of each other in the fourth column from the left). In contrast, constraints on cosmology from the intrinsic galaxy clustering information alone (gg) weaken significantly as greater flexibility is included in the galaxy bias model. Specifically, the constraint from the gg signal drops by more than 20 in the logarithmic FoM_{TOT} when increasing the number of nuisance parameters per bias term by a factor of five (compare the solid and dashed lines in the third column from the left).
The nn correlations, containing both clustering and magnification, unsurprisingly have a FoM_{TOT} in between (second column from the left). Hence, if the galaxy bias is well known, the galaxy clustering signal dominates the information from the nn correlations, but for very flexible bias models lensing magnification does allow more information on cosmology to be extracted, as compared to what might be expected if magnification did not take place (compare the dashed lines in the second and third columns).
The two right hand columns of Fig.7 show the FoM_{TOT} obtained when the full range of angular scales is used in the power spectra, thus including information from nonlinear scales. As already seen in Fig.2, the clustering constraints improve strongly when adding the signal from the nonlinear regime. The magnification contribution remains trustworthy far into the nonlinear regime, but is of course also affected by the ℓcuts. When using the full range of angular scales the FoM_{TOT} of the mm signal is increased by about 7. The cosmological constraints from intrinsic galaxy clustering (gg) alone are now much tighter than from the magnification effect alone (mm), when the smaller number of bias parameters are used (compare solid lines in the two right hand columns). However, the cosmological constraints from the maximally flexible model are still tighter from magnification alone than from galaxy clustering (dashed lines in the right hand two columns).
Note that the constraints from the isolated mm and gg signals are hypothetical – the full nn correlations including both terms are the only true observables. It is only possible to separate the contributions in an approximate fashion when making use of their characteristic scaling with redshift, see Fig.1. But all of this available information is already included in our nn results when a large enough number of tomographic redshift bins are used.
We now assess the cosmological information available in the different subsets of observables when performing the Fisher matrix analysis for all correlations simultaneously. To this end, we split up the summation of the Fisher matrix (42) into three parts, corresponding to ϵϵ, nn, and nϵ correlations. We also consider pairs of observables e.g. ϵϵ, and nn. The Fisher matrices of each part are then inverted separately to yield individual parameter errors and figures of merit. Thereby we split the total information into subsets of the data vector, which could in principle be observed independently.
However, the figures of merit we compute do not correspond to those which one would obtain for an independent analysis of ϵϵ, nn, or nϵ correlations because we extract the relevant rows and columns of the inverse of the full covariance matrix (39) to insert into the Fisher matrix calculation (42). Due to the inversion, the covariance terms of the different subsets mix, thus accounting for the crosscorrelations between the observables in the different subsets. This is desirable for our purposes because these terms add together in the full calculation. Consequently, it is possible that the FoM of a subset is larger than the one obtained for the complete data vector if there are anticorrelations with observables of other subsets which are not taken into account due to the splitting. Formally speaking, this means that due to anticorrelations of different observables, offdiagonal terms of the Fisher matrix can have negative entries, which produce negative terms in the sum in (42). We indeed observe this behaviour for one of the subsets.
Fig. 8 Upper panel: the FoM_{TOT} for different subsets of correlations that are observable from a galaxy survey with galaxy shape and number density measurements. The Fisher matrix has been computed for the total data set containing all correlations, i.e. these results represent the information content of the subsets as a part of the total signal. Black dashed bars correspond to a marginalisation over N_{K} = N_{Z} = 7 nuisance parameters, grey solid bars to N_{K} = N_{Z} = 3 nuisance parameters. The labels on the abscissa indicate the different combinations of correlations used, where ϵn stands for the crosscorrelations between number density and ellipticity. Lower panel: difference Δ FoM, see (48), for the same subsets as above. This difference can be understood as a measure of the depletion of information in the subsets due to the marginalisation over nuisance parameters. 
The FoM_{TOT} for the subsets of ϵϵ, nn, and nϵ correlations, as well as all possible combinations thereof, are given in the upper panel of Fig.8, again for the two parametrisations N_{K} = N_{Z} = {3;7}. As part of the full set of observables, the ϵϵ correlations, governed by the cosmic shear signal, yield the highest FoM_{TOT} and contribute most to the cosmological information in the full set (second column). They are followed by crosscorrelations between number density and ellipticity, which have a moderate FoM_{TOT} (fourth column). This occurs in spite of the smaller range in angular frequency used for this observable. In addition, the least amount of information is lost when switching to the more flexible bias model. The subset of nn correlations has by far the lowest FoM_{TOT}, which becomes even more pronounced when the bias model has more nuisance parameters. This hierarchy in FoM is also mirrored in the results for the different combinations of two subsets.
To get a more explicit measure of the effect of the marginalisation over the galaxy bias and intrinsic alignments nuisance parameters, we compute the quantity (48)for each of the subsets. The first term in (48) is of a similar form as the FoM_{TOT}, but the cosmological elements are extracted before the inverse or determinant is taken i.e. it does not include the marginalisation over nuisance parameters. Hence, ΔFoM quantifies the depletion in FoM due to the marginalisation. As is evident from the lower panel in Fig.8, the susceptibility of the subsets to the nuisance parameters is closely related to their contribution to the FoM_{TOT} of the full set. Accordingly, ΔFoM is largest for the nn subset while the ellipticitynumber density crosscorrelations have substantially smaller ΔFoM although they contain all four bias terms. The effect of intrinsic alignments on the ϵϵ subset is relatively small when compared to the effect of galaxy bias on the nn correlations.
Since the ϵϵ correlations only depend on intrinsic alignments (b_{I} and r_{I}) whereas nn correlations only feature galaxy bias (b_{g} and b_{g}), combining them helps little in reducing ΔFoM, see the column labelled “ϵϵ + nn”. Both lines in this column remain higher than for the other combinations “ϵϵ + ϵn” and “nn + ϵn”. The ellipticitynumber density crosscorrelations instead have great potential in breaking degeneracies between cosmological and nuisance parameters. In particular, adding their subsets to the ϵϵ correlations further decreases ΔFoM below the ϵϵonly value for the most flexible bias term model. This synergy is presumably related to the internal calibration between the IG and gI signals, as investigated by Zhang (2010). Finally, the full set of shape and number density observables clearly calibrates the nuisance parameters best as it produces the smallest ΔFoM for both the most flexible and the more rigid bias term models.
5. Conclusions
In this work we studied the joint analysis of galaxy number density and shape correlations to constrain cosmological parameters in presence of contaminations by the intrinsic alignment of galaxies and the galaxy bias. We considered the signals due to gravitational shear, intrinsic shear, intrinsic galaxy clustering, and lensing magnification, explicitly computing all possible twopoint correlations thereof. We introduced a twodimensional grid parametrisation to account for the unknown scale and redshift dependence of both intrinsic alignments and galaxy bias. Further nuisance parameters were used to describe the uncertainty in the width of the photoz binwise redshift distributions and in the slope of the galaxy luminosity functions within each bin.
Our Fisher matrix analysis demonstrates that the simultaneous use of ellipticity correlations, numberdensity correlations, and in particular ellipticitynumberdensity crosscorrelations allows for a substantial amount of internal calibration of the bias terms. With flexible models that contain in total more than 200 nuisance parameters we can recover the volume of the error ellipsoid in parameter space when compared to assuming a pure gravitational lensing signal and using just the ellipticity correlation information. The dark energy parameters w_{0} and w_{a} suffer more than other parameters on marginalisation over nuisance parameters, so that 56% of the FoM_{DETF} are retained for a Euclidlike survey in this most flexible setup considered. The FoM_{DETF} for the combined set of shape and number density correlations is close to the pure lensing FoM_{DETF} if we choose a model which uses about 100 nuisance parameters to describe intrinsic alignments and galaxy bias. Our approach also proves beneficial for upcoming groundbased surveys with DESlike survey characteristics. The benefit is greatest for the more ambitious survey. In addition, we assumed the slopes of the galaxy luminosity functions in each photoz bin to be unknown nuisance parameters and found that they are well calibrated internally from the data.
The information which we added on top of the standard cosmic shear analysis comes without any extra cost since galaxy number density measurements are directly available from imaging data. Given our encouraging findings, we hence suggest that the joint consideration of galaxy shape and number density information could become the standard technique whenever intrinsic alignments are suspected to make a significant contribution to the cosmic shear signal.
When interpreting our results, one has to keep in mind that our grid parametrisation of the bias terms has limited flexibility. On increasing the number of nuisance parameters in the grids, the curves for the figures of merit flatten off. However we are at present not able to say whether they approach a lower limit or continue to decrease for a large number of grid nodes. This issue, which is of considerable theoretical and practical interest, is currently under investigation, as well as a comparison with the performance of removal techniques such as nulling (Joachimi & Schneider 2008, 2009).
Due to the finite number of grid points the resulting angular power spectra will not span perfectly the range of possible physical models. Moreover, we found that of order 10 nuisance parameters in k per grid node in redshift do not suffice to represent all relevant functional forms of the galaxy bias. Therefore we removed all observables from the parameter estimation which have a significant contribution from the galaxy clustering signal at large k where the signal is strongest and the bias least known. However, since in the near future we expect to have at least coarse, but reliable knowledge about the functional forms of the galaxy bias on linear scales and the intrinsic alignment signals, the models used in our approach should still yield realistic results.
While the loglinear grid parametrisation is fairly general and intuitive, it may not be the most efficient way to represent freedom in intrinsic alignments and the galaxy bias. As all bias terms originate from physical processes, they are smooth and do not feature strong oscillations or isolated peaks. Thus, we presume that the bias terms can efficiently be parametrised in terms of complete sets of smooth functions such as the Fourier, Chebyshev, or Legendre series. Truncating the higher orders of these series will only limit the model to represent highly oscillatory or smallscale features. These parametrisations might therefore be more efficient in comprising the full set of realistic bias terms (see also Kitching et al. 2009) for a given number of nuisance parameters, this number in turn being dictated by the available computational power. The latter will play a critical role when the approach suggested here is performed in a full likelihood analysis with several hundreds of nuisance parameters. We add that it might be necessary to furthermore consider the bias terms individually for different galaxy types and luminosities because it is known that both the intrinsic alignments and the galaxy bias vary considerably with galaxy type and luminosity (for recent examples see Mandelbaum et al. 2006; Hirata et al. 2007; McCracken et al. 2008; Swanson et al. 2008; Simon et al. 2009; Cresswell & Percival 2009; Wang et al. 2007).
Redshift information is crucial for discriminating the different signals that contribute to the observables. We investigated the dependence of the parameter constraints from the joint set of correlations on characteristics of the redshift distributions. We confirm the observation by foregoing works that the number of photoz bins needed to retrieve the bulk of information about cosmology increases by at least a factor of two when using only ellipticity correlations and marginalising over the intrinsic alignment signals. Using the complete set of correlations, the figures of merit do not level off so quickly, but continue to increase as the number of tomographic bins is increased. This can be explained by the narrower kernel in the redshift integrations that link the projected and the threedimensional power spectra involving galaxy numberdensity signals.
Moreover the figures of merit decrease substantially more slowly as a function of the photometric redshift dispersion σ_{ph} when adding galaxy number density information to the data. However, in our marginalisation we only used a single parameter accounting for the uncertainty in the spread of all binwise redshift distributions. This way the redshift crossterms of the numberdensity correlations efficiently calibrate the shape of the redshift distributions. Therefore this result alone cannot be interpreted as a potential relaxation of the requirements for photometric redshift accuracy in cosmological surveys featuring cosmic shear. Conversely, the continuing increase in the figure of merit as a function of N_{zbin} makes an even larger number of wellseparated and compact redshift distributions desirable in the case of the joint data set.
The recent work by Zhang et al. (2010) uses densityellipticity correlations to selfcalibrate photometric redshift parameters but a full approach may need to simultaneously deal with intrinsic alignments and photometric redshift properties if sufficient spectra cannot be obtained to calibrate the photometric redshifts independently. In a forthcoming paper we will investigate into a more realistic approach including uncertainty and outliers in the redshift distributions of each individual photoz bin, as well as the benefits of spectroscopic redshift information for a subsample of the galaxy catalogue.
We emphasise that in the approach suggested here all considered correlations help in constraining cosmological parameters. Hence, none of the contributions is regarded as a systematic signal, and the different levels of uncertainty concerning the exact form of the signals are represented by nuisance parameters. By means of the joint analysis we increase the statistical power via the crosscalibration abilities of the signals, and reduce the risk of undetected systematic effects one faces when considering signals such as cosmic shear, galaxy clustering, or lensing magnification individually.
This integrative ansatz is not limited to twopoint correlations, but can be generalised to the threepoint level in a straightforward way. Future surveys will provide excellent data for studying threepoint correlations whose exploitation can break parameter degeneracies and improve constraints considerably (e.g. Takada & Jain 2004). Our knowledge about intrinsic alignments and galaxy bias at the threepoint level is currently even more limited than in the twopoint case (however see Semboloni et al. 2008), so that a joint investigation of shape and number density observables, including a general parametrisation of the various bias terms, may be an appropriate way forward. However, limitations due to computational power will most likely play a dominant role in this case.
One can also think of incorporating further sets of observables into the analysis. For instance, planned surveys like Euclid will also include a spectroscopic survey of a subset of the galaxies in order to determine the underlying redshift distributions of the photoz bins. The spectroscopic data can be used to measure baryonic acoustic oscillations (note that the acoustic peaks are not included in our galaxy numberdensity correlations) and galaxy peculiar velocities. The latter allow for the measurement of redshift space distortions. Similar to the case considered here, the joint analysis of galaxy number density, shape, and velocity information (Guzik et al. 2010) will efficiently crosscalibrate nuisance terms such as the galaxy bias and tighten constraints on the cosmological model.
Panoramic Survey Telescope & Rapid Response System, http://panstarrs.ifa.hawaii.edu
Dark Energy Survey, https://www.darkenergysurvey.org
Large Synoptic Survey Telescope, http://www.lsst.org
Joint Dark Energy Mission, http://jdem.gdfc.nasa.gov
In fact we expect the for to be slightly above the lensingonly case since the intrinsic alignment signals also constrain cosmology if their form is perfectly known, see also Bridle & King (2007).
Acknowledgments
We would like to thank Peter Schneider, Tom Kitching, Adam Amara, Gary Bernstein, Bhuvnesh Jain and Ofer Lahav for helpful discussions. B.J. acknowledges support by the Deutsche Telekom Stiftung and the BonnCologne Graduate School of Physics and Astronomy. S.L.B. thanks the Royal Society for support in the form of a University Research Fellowship.
References
 Abdalla, F. B., Amara, A., Capak, P., et al. 2007, MNRAS, 387, 969 [Google Scholar]
 Albrecht, A., Bernstein, G., Cahn, R., et al. 2006, unpublished [arXiv:astroph/0609591] [Google Scholar]
 Annis, J., Bridle, S. L., Castander, F. J., et al. 2005, unpublished [arXiv:astroph/0510195] [Google Scholar]
 Bacon, D. J., Réfrégier, A. R., & Ellis, R. S. 2000, MNRAS, 318, 625 [NASA ADS] [CrossRef] [Google Scholar]
 Bartelmann, M., & Schneider, P. 2001, Phys. Reports, 340, 291 [NASA ADS] [CrossRef] [Google Scholar]
 Benjamin, J., Heymans, C., Semboloni, E., et al. 2007, MNRAS, 381, 702 [NASA ADS] [CrossRef] [Google Scholar]
 Bernstein, G. M. 2009, ApJ, 695, 652 [NASA ADS] [CrossRef] [Google Scholar]
 Bernstein, G. M., & Huterer, D. 2010, MNRAS, 401, 1399 [NASA ADS] [CrossRef] [Google Scholar]
 Blake, C., & Bridle, S. 2005, MNRAS, 363, 1329 [NASA ADS] [CrossRef] [Google Scholar]
 Blake, C., Collister, A., Bridle, S., & Lahav, O. 2007, MNRAS, 374, 1527 [NASA ADS] [CrossRef] [Google Scholar]
 Brainerd, T., Agustsson, I., Madsen, C. A., & Edmonds, J. A. 2009, ApJ, submitted [arXiv:astroph/0904.3095] [Google Scholar]
 Bridle, S., & Abdalla, F. B. 2007, ApJ, 655, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Bridle, S., & King, L. 2007, NJPh, 9, 444 [NASA ADS] [CrossRef] [Google Scholar]
 Broadhurst, T. J., Taylor, A. N., & Peacock, J. A. 1995, ApJ, 438, 49 [NASA ADS] [CrossRef] [Google Scholar]
 Brown, M. L., Taylor, A. N., Hambly, N. C., & Dye, S. 2002, MNRAS, 333, 501 [NASA ADS] [CrossRef] [Google Scholar]
 Cacciato, M., van den Bosch, F. C., More, S., Li, R., et al. 2009, MNRAS, 394, 929 [NASA ADS] [CrossRef] [Google Scholar]
 Catelan, P., Kamionkowski, M., & Blandford, R. D. 2001, MNRAS, 320, 7 [Google Scholar]
 Chevallier, M., & Polarski, D. 2001, Int. J. Mod. Phys., D10 [Google Scholar]
 Cresswell, J. G., & Percival, W. J. 2009, MNRAS, 392, 682 [NASA ADS] [CrossRef] [Google Scholar]
 Crittenden, R. G., Natarajan, P., Pen, U., & Theuns, T. 2001, ApJ, 559, 552 [NASA ADS] [CrossRef] [Google Scholar]
 Croft, R. A. C., & Metzler, C. A. 2000, ApJ, 545, 561 [NASA ADS] [CrossRef] [Google Scholar]
 Das, S., & Spergel, D. N. 2009, Phys. Rev. D, 79, 043509 [NASA ADS] [CrossRef] [Google Scholar]
 Dolney, D., Jain, B., & Takada, M. 2006, MNRAS, 366, 884 [NASA ADS] [Google Scholar]
 Eisenstein, D. J., & Hu, W. 1998, ApJ, 496, 605 [NASA ADS] [CrossRef] [Google Scholar]
 Fu, L., Semboloni, E., Hoekstra, H., et al. 2008, A&A, 479, 9 [Google Scholar]
 Guzik, J., & Seljak, U. 2001, MNRAS, 321, 439 [NASA ADS] [CrossRef] [Google Scholar]
 Guzik, J., & Seljak, U. 2002, MNRAS, 335, 311 [NASA ADS] [CrossRef] [Google Scholar]
 Guzik, J., Jain, B., & Takada, M. 2010, Phys. Rev. D, 81, 023503 [NASA ADS] [CrossRef] [Google Scholar]
 Heavens, A., Réfrégier, A., & Heymans, C. 2000, MNRAS, 319 [Google Scholar]
 Hetterscheidt, M., Simon, P., Schirmer, M., et al. 2007, A&A, 468, 859 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Heymans, C., & Heavens, A. 2003, MNRAS, 339, 711 [NASA ADS] [CrossRef] [Google Scholar]
 Heymans, C., Brown, M., Heavens, A., et al. 2004, MNRAS, 347, 895 [NASA ADS] [CrossRef] [Google Scholar]
 Hirata, C. M., & Seljak, U. 2004, Phys. Rev. D, 70, 063526 [NASA ADS] [CrossRef] [Google Scholar]
 Hirata, C. M., Mandelbaum, R., Ishak, M., et al. 2007, MNRAS, 381, 1197 [NASA ADS] [CrossRef] [Google Scholar]
 Hoekstra, H., Mellier, Y., van Waerbeke, L., et al. 2006, ApJ, 647, 116 [NASA ADS] [CrossRef] [Google Scholar]
 Hu, W. 1999, ApJ, 522, 21 [NASA ADS] [CrossRef] [Google Scholar]
 Hu, W. 2002, Phys. Rev. D, 65, 023003 [NASA ADS] [CrossRef] [Google Scholar]
 Hu, W., & Jain, B. 2004, Phys. Rev. D, 70, 043009 [NASA ADS] [CrossRef] [Google Scholar]
 Huterer, D., Takada, M., Bernstein, G., & Jain, B. 2006, MNRAS, 366, 101 [NASA ADS] [CrossRef] [Google Scholar]
 Jarvis, M., Jain, B., Bernstein, G., & Dolney, D. 2006, ApJ, 644, 71 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Jing, Y. P. 2002, MNRAS, 335, 89 [NASA ADS] [CrossRef] [Google Scholar]
 Joachimi, B., & Schneider, P. 2008, A&A, 488, 829 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Joachimi, B., & Schneider, P. 2009, A&A, 507, 105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Joachimi, B., Schneider, P., & Eifler, T. 2008, A&A, 477, 43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Johnston, D. E., Sheldon, E. S., Tasitsiomi, A., et al. 2007, ApJ, 656, 27 [NASA ADS] [CrossRef] [Google Scholar]
 Kaiser, N. 1992, ApJ, 388, 272 [NASA ADS] [CrossRef] [Google Scholar]
 Kaiser, N., Wilson, G., & Luppino, G. 2000, unpublished [arXiv:astroph/0003338] [Google Scholar]
 King, L. J. 2005, A&A, 441, 47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 King, L. J., & Schneider, P. 2002, A&A, 396, 411 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 King, L. J., & Schneider, P. 2003, A&A, 398, 23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kitching, T. D., Heavens, A. F., Verde, L., Serra, P., & Melchiorri, A. 2008a, Phys. Rev. D, 77, 103008 [NASA ADS] [CrossRef] [Google Scholar]
 Kitching, T. D., Taylor, A. N., & Heavens, A. F. 2008b, MNRAS, 389, 173 [NASA ADS] [CrossRef] [Google Scholar]
 Kitching, T. D., Amara, A., Abdalla, F. B., Joachimi, B., & Réfrégier, A. 2009, MNRAS, 399, 2107 [NASA ADS] [CrossRef] [Google Scholar]
 Krause, E., & Hirata, C. 2010, A&A, [arXiv:astroph/0910.3786] [Google Scholar]
 Lahav, O., & Suto, Y. 2004, Liv. Rev. Rel., 7, 8 [Google Scholar]
 Laureijs, R., et al. 2009, Euclid Assessment Study Report for the ESA Cosmic Visions, ESA/SRE(2009)2 [arXiv:astroph/0912.0914] [Google Scholar]
 Lee, J., & Pen, U.L. 2000, ApJ, 532, L5 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Linder, E. V. 2003, Phys. Rev. Lett., 90, 091301 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Liu, C. T., Capak, P., Mobasher, B., et al. 2008, ApJ, 672, 198 [NASA ADS] [CrossRef] [Google Scholar]
 Ma, Z., Hu, W., & Huterer, D. 2006, ApJ, 636, 21 [NASA ADS] [CrossRef] [Google Scholar]
 Mackey, J., White, M., & Kamionkowski, M. 2002, MNRAS, 332, 788 [NASA ADS] [CrossRef] [Google Scholar]
 Mandelbaum, R., Hirata, C. M., Ishak, M., Seljak, U., & Brinkmann, J. 2006, MNRAS, 367, 611 [NASA ADS] [CrossRef] [Google Scholar]
 Mandelbaum, R., Blake, C., Bridle, S., Abdalla, F. B., et al. 2010, MNRAS, accepted [Google Scholar]
 Massey, R., Rhodes, J., Leauthaud, A., et al. 2007, ApJS, 172, 239 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 McCracken, H. J., Ilbert, O., Mellier, Y., Bertin, E., et al. 2008, A&A, 479, 321 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 McDonald, P., Trac, H., & Contaldi, C. 2006, MNRAS, 366, 547 [NASA ADS] [CrossRef] [Google Scholar]
 Newman, J. A. 2008, ApJ, 684, 88 [NASA ADS] [CrossRef] [Google Scholar]
 Okumura, T., & Jing, Y. P. 2009, ApJ, 694, L83 [NASA ADS] [CrossRef] [Google Scholar]
 Okumura, T., Jing, Y. P., & Li, C. 2009, ApJ, 694, 214 [NASA ADS] [CrossRef] [Google Scholar]
 Padmanabhan, N., Schlegel, D. J., Seljak, U., Makarov, A., et al. 2007, MNRAS, 378, 852 [NASA ADS] [CrossRef] [Google Scholar]
 Peacock, J. A., Schneider, P., Efstathiou, G., et al. 2006, in ESAESO Working Group on Fundamental Cosmology, ed. E. J. A. Peacock et al. [Google Scholar]
 Pen, U.L., Lee, J., & Seljak, U. 2000, ApJ, 543, L107 [NASA ADS] [CrossRef] [Google Scholar]
 Rassat, A., Amara, A., Amendola, L., et al. 2008, MNRAS, submitted [arXiv:astroph/0810.0003] [Google Scholar]
 Réfrégier, A., Massey, R., Rhodes, J., et al. 2004, AJ, 127, 3102 [NASA ADS] [CrossRef] [Google Scholar]
 Réfrégier, A., Boulade, O., Mellier, Y., et al. 2006, SPIE, 6265, 62651Y [CrossRef] [Google Scholar]
 Réfrégier, A., Amara, A., Kitching, T., & Rassat, A. 2008, A&A, submitted [arXiv:astroph/0810.1285] [Google Scholar]
 Schäfer, B. M. 2009, IJMPD, 18, 173 [Google Scholar]
 Schmidt, F., Rozo, E., Dodelson, S., Hui, L., & Sheldon, E. 2009, Phys. Rev. Lett., 103, 051301 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Schneider, P. 2006, in SaasFee Advanced Course 33: Gravitational Lensing: Strong, Weak and Micro, ed. G. Meylan, P. Jetzer, P. North, P. Schneider, C. S. Kochanek, & J. Wambsganss, 269 [Google Scholar]
 Schneider, P., & Rix, H.W. 1997, ApJ, 474, 25 [NASA ADS] [CrossRef] [Google Scholar]
 Schneider, M. D., & Bridle, S. 2010, MNRAS, 402, 2127 [NASA ADS] [CrossRef] [Google Scholar]
 Schneider, M., Knox, L., Zhan, H., & Connolly, A. 2006, ApJ, 651, 14 [NASA ADS] [CrossRef] [Google Scholar]
 Schrabback, T., Erben, T., Simon, P., et al. 2007, A&A, 468, 823 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Schrabback, T., Hartlap, J., Joachimi, B., et al. 2010, A&A, 516, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Seitz, C., & Schneider, P. 1997, A&A, 318, 687 [NASA ADS] [Google Scholar]
 Seljak, U. 2002, MNRAS, 337, 769 [NASA ADS] [CrossRef] [Google Scholar]
 Seljak, U., Makarov, A., Mandelbaum, R., Hirata, C. M., et al. 2005, Phys. Rev. D, 71, 043511 [NASA ADS] [CrossRef] [Google Scholar]
 Semboloni, E., Mellier, Y., van Waerbeke, L., et al. 2006, A&A, 452, 51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Semboloni, E., Heymans, C., van Waerbeke, L., & Schneider, P. 2008, MNRAS, 388, 991 [NASA ADS] [CrossRef] [Google Scholar]
 Simon, P., King, L. J., & Schneider, P. 2004, A&A, 417, 873 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Simon, P., Hetterscheidt, M., Wolf, C., et al. 2009, MNRAS, 398, 807 [NASA ADS] [CrossRef] [Google Scholar]
 Smail, I., Ellis, R. S., & Fitchett, M. 1994, MNRAS, 270, 245 [NASA ADS] [CrossRef] [Google Scholar]
 Smith, R. E., Peacock, J. A., Jenkins, A., et al. 2003, MNRAS, 341, 1311 [NASA ADS] [CrossRef] [Google Scholar]
 Swanson, M. E. C., Tegmark, M., Blanton, M., & Zehavi, I. 2008, MNRAS, 385, 1635 [NASA ADS] [CrossRef] [Google Scholar]
 Takada, M., & Jain, B. 2004, MNRAS, 348, 897 [NASA ADS] [CrossRef] [Google Scholar]
 Takada, M., & White, M. 2004, ApJ, 601, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Tegmark, M., Taylor, A. N., & Heavens, A. F. 1997, ApJ, 480, 22 [NASA ADS] [CrossRef] [Google Scholar]
 Tereno, I., Doré, O., van Waerbeke, L., & Mellier, Y. 2005, A&A, 429, 383 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Thomas, S. A., Abdalla, F. B., & Weller, J. 2009, MNRAS, 395, 197 [NASA ADS] [CrossRef] [Google Scholar]
 van den Bosch, F. C., Abel, T., Croft, R. A. C., Hernquist, L., & White, S. D. M. 2002, ApJ, 576, 21 [NASA ADS] [CrossRef] [Google Scholar]
 van Waerbeke, L. 2010, MNRAS, 401, 2093 [NASA ADS] [CrossRef] [Google Scholar]
 van Waerbeke, L., Mellier, Y., Erben, T., et al. 2000, A&A, 358, 30 [NASA ADS] [Google Scholar]
 Wang, Y., Yang, X.H., Mo, H. J., & van den Bosch, F. C. 2007, ApJ, 664, 608 [NASA ADS] [CrossRef] [Google Scholar]
 Wittman, D. M., Tyson, J. A., Kirkman, D., Antonio, I. D., & Bernstein, G. 2000, Nature, 405, 143 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Wolf, C., Meisenheimer, K., Rix, H.W., et al. 2003, A&A, 401, 73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Yoo, J., Tinker, J. L., Weinberg, D. H., et al. 2006, ApJ, 652, 26 [NASA ADS] [CrossRef] [Google Scholar]
 Zhan, H. 2006, JCAP, 0608, 008 [Google Scholar]
 Zhang, P. 2010, ApJ, 720, 1090 [NASA ADS] [CrossRef] [Google Scholar]
 Zhang, P., & Pen, U.L. 2005, Phys. Rev. Lett., 95, 241302 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Zhang, P., & Pen, U.L. 2006, MNRAS, 367, 169 [NASA ADS] [CrossRef] [Google Scholar]
 Zhang, P., Pen, U.L., & Bernstein, G. 2010, MNRAS, 405, 359 [NASA ADS] [Google Scholar]
Appendix A: Projected lensing magnification signal
This appendix details the derivation of the enhancement or depletion of projected galaxy counts due to lensing magnification. The number density of galaxies n, counted above a flux threshold S at angular position θ and comoving distance χ is altered by gravitational lensing according to (A.1)where n_{0} is the original galaxy number density, and where μ denotes the magnification (Bartelmann & Schneider 2001). One assumes that the galaxy luminosity function close to the flux limit of the survey can locally be written as a power law, n( > S,f_{K}(χ)θ,χ) ∝ S^{ − α(χ)}. The slope α depends on the lineofsight distance, or equivalently, redshift, but should not depend on angular dimensions due to isotropy. However, it is a function of the magnitude limit in the observed filter, in this work denoted by r_{lim}. This dependence is dealt with in Sect.3.2, but for ease of notation we drop r_{lim} as an argument of α for the remainder of this section. Plugging the powerlaw form of the luminosity function into (A.1) yields (A.2)Again following Bartelmann & Schneider (2001), one can approximate the magnification in the weak lensing regime as μ ≈ 1 + 2 κ_{G}. Since κ_{G} ≪ 1, we can in addition do a Taylor approximation to arrive at (A.3)Defining the excess galaxy density contrast due to magnification effects as (A.4)one obtains for the corresponding projected density contrast (A.5)In exact analogy to the standard derivation of (9) one can now insert the relation between the convergence and the threedimensional matter density contrast, (A.6)which, after swapping the order of integration and the names of the integration variables, yields (A.7)Here, we have defined the weight (A.8)Given that the slope of the luminosity function should be a smooth function of comoving distance, α(χ) varies only weakly over the range of the integration in (A.8), being determined by the distribution p^{(i)}(χ), which has relatively compact support. Hence, the mean value theorem constitutes a good approximation, so that we can write (A.9)where we define α^{(i)} to be the slope of the luminosity function, evaluated at the median redshift of the photometric bin i. Inserting (A.9) into (A.7) results in (13), which we employ throughout this work.
Appendix B: Limber equations of cosmological signals
For reference, we have collected here the Fourierspace Limber equations of all projected power spectra that contribute to the cosmological signals considered in this work: The terminology of both projected and threedimensional power spectra is summarised in Table 1. The weights that enter the foregoing equations are the probability distribution of galaxies with comoving distance p^{(i)}(χ), and the ones defined in (10), and (A.8).
Appendix C: The slope of the galaxy luminosity function
To compute the lensing magnification signal, we need to model the powerlaw slope of the cumulative galaxy luminosity function at the magnitude limit of the galaxy number density samples under consideration. We base our modelling on Blake & Bridle (2005) who have determined galaxy redshift distributions for a given magnitude limit, using COMBO17 luminosity functions for the SDSS r filter (Wolf et al. 2003).
They fitted these distributions with functions of the form (29) with β = 1.5, using two free parameters, the redshift scaling and the normalisation given by the galaxy surface density Σ_{0}. We set Σ_{0} and as a function of survey depth making use of Table 1 of Blake & Bridle (2005) and fit a power law to each quantity as a function of the limiting magnitude r_{lim}, (C.1)where we find good fits using Σ_{0,c} = 9.83, η_{Σ} = 19, and . This allows us to extrapolate beyond the range of their Table, which stops at r_{lim} = 24.
Note that our definition of α is not the same as the exponent of the Schechter function (sometimes also denoted α, see e.g. Wolf et al. 2003, Eq. (4)). Our α is the negative of the slope of the cumulative luminosity function. Therefore for comparison to the Schechter function, in the faint galaxy (powerlaw) regime, one must take the negative of our alpha and subtract unity. Note that our use of the cumulative luminosity function is consistent with the literature on lensing magnification (e.g. Schmidt et al. 2009). Our typical α values of around 0.5 are also therefore consistent with the luminosity function literature (e.g. Liu et al. 2008 find Schechter function exponent values between −1 and −2 depending on the spectral type from COSMOS).
Fig. C.1 Top panel: slope of the galaxy luminosity function α as a function of redshift, shown for different magnitude limits r_{lim} as indicated in the legend. Bottom panel: fit parameters a_{1} (solid), a_{2} (dotted), and a_{3} (dashed) as a function of r_{lim}. These parameters are obtained from the polynomial fit in (C.2). 
We are interested in the slope of the luminosity function α(z,r_{lim}) at the cosmic shear survey magnitude limit. This slope is a function of redshift and magnitude limit. From (29) and (C.1) we have the number of galaxies as a function of redshift and magnitude limit. We convert each magnitude limit into a flux limit S and set the number of galaxies above the flux limit equal to S^{ − α} (see also Appendix A). The resulting curves for α(z,r_{lim}) are shown in Fig.C.1, upper panel.
For convenience, and to extrapolate the slope α(z,r_{lim}) to values z > 1, we now provide a fitting formula . First we expand the slope using a polynomial in redshift, with coefficients that depend on the limiting magnitude. Then we find an approximate equation for these coefficients as a function of limiting magnitude. This results in equations for the slope as a function of redshift and magnitude limit in terms of 15 numbers given in Table C.1.
We fit the slope of the luminosity function as a function of redshift with a secondorder polynomial (C.2)The polynomial coefficients a_{i} are functions of the limiting magnitude, and are shown in Fig. C.1, lower panel. We find that these coefficients are in turn well fit by a function of the form (C.3)with parameters b_{ij} given in Table C.1. We chose not to use b_{i2} as a free parameter for the fit, but set it as b_{i2} = ±1, to determine the sign of the term in parentheses. By means of (C.2) and (C.3) we have condensed the dependence of α on redshift and r_{lim} into the 15 parameters summarised in Table C.1.
In Fig. C.2 we plot the relative accuracy of this set of fit formulae with respect to α(z,r_{lim}) as given in Fig. C.1, lower panel. Over the dominant part of the considered parameter space the fit formulae provide an excellent approximation, which deviates less than 1% from the original fits (C.1). Thus one can expect that within the framework of this approach (C.2) and (C.3) extrapolate α(z,r_{lim}) reasonably well to z > 1. Significantly larger deviations can only be found for the brightest limiting magnitudes at redshifts z ≲ 0.1, a region of the parameter plane which is irrelevant for a competitive cosmological survey.
Fig. C.2 Percentage deviation of the simplified fit as defined in (C.2) and (C.3) from the original fits described by (C.1) and shown in the upper panel of Fig. C.1. The contour levels correspond to the percentages given on the curves. Negative deviations are indicated by dashed contours. Note that across most of the parameter plane the modulus of the deviation is less than 1%. 
We have used and extrapolated the COMBO17 luminosity functions despite the incompleteness beyond R = 23 and relative unreliability in the redshift range 1.2 < z < 2 due to the lack of spectral features in the observing bands used. This could be improved by using deeper data which have infrared observations, for example COSMOS ground and space data and CFHTLSWide. We emphasise that this is simply used for our choice of fiducial model, and within the modelling it is assumed that the luminosity function slopes are an unknown function of redshift and are marginalised over. We do not expect the choice of fiducial model to have a large effect on our analysis. Note that the number of galaxies as a function of redshift is inserted into the analysis using the assumptions made by the project teams and is not linked to our luminosity function calculations.
It would be more elegant to simultaneously derive the redshift distribution and luminosity function slope from luminosity functions derived from data as a function of redshift. However this would make it difficult to compare our results with those in the literature from the imaging survey project teams, who have already chosen a galaxy redshift distribution. We recommend that the project teams also make available luminosity function slopes to be used in analyses such as these.
Appendix D: Modelling the effect of dark energy on nonlinear growth
In this appendix we summarise the approach made in the publicly available icosmo code^{7} (Réfrégier et al. 2008) to model the nonlinear evolution of structure in presence of a dark energy equation of state w(z) ≠ −1, which we have adopted for this work. It is based upon a modification of the halofit routine of Smith et al. (2003) in which the interpolation between open and flat cosmological models is determined by the parameter f = Ω_{Λ}/(1 − Ω_{m}), where f = 0 corresponds to open universes without dark energy (we assume Ω_{m} < 1 here), and f = 1 to flat ΛCDM models.
Now the fact is exploited that certain variable dark energy models mimic the expansion history of open CDM universes. Using w_{0} = −1/3 and w_{a} = 0 in (28) for a flat ΛCDM model, it is readily seen that the same Hubble parameter is obtained as for an open CDM universe with identical Ω_{m} but without dark energy, where (28) plays the role of the curvature term. Motivated by this coincidence, the interpolation parameter f in the halofit routine is replaced by (D.1)where w(z) = w_{0} + w_{a}z/(1 + z), and f = Ω_{Λ}/(1 − Ω_{m}) as before. If a model does not feature dark energy, f′ = f = 0. For a flat ΛCDM model the interpolation now takes place between w(z) = −1/3, mimicking the case of an open CDM in the original halofit (f′ = 0), and the cosmological constant w(z) ≡ −1 (f′ = 1).
The performance of this simplistic ansatz has been tested in Fig.10 of Schrabback et al. (2010), finding fair agreement with the fit formula of McDonald et al. (2006). Note that the fits to the simulations by McDonald et al. (2006), which include a dependence on w_{0}, are not suitable for direct use in this article due to the limited range in cosmological parameters, most notably σ_{8}.
All Tables
Marginalised parameter errors and the two figures of Merit for a survey with ten photometric redshift bins used for tomography N_{zbin} = 10 and a photometric redshift uncertainty parameter σ_{ph} = 0.05.
All Figures
Fig. 1 Fiducial power spectra for all considered correlations. The upper right panels depict the contributions to ϵϵ (in black) and nn (in magenta) correlations. The lower left panels show the contributions to correlations between number density fluctuations and ellipticity. Since we only show correlations with i ≤ j, we make in this plot a distinction between nϵ (in red; number density contribution in the foreground, e.g. gG) and ϵn (in blue; number density contribution in the background, e.g. Gg) correlations. In each subpanel a different tomographic redshift bin correlation is shown. For clarity only odd bins are displayed. In the upper right panels the usual cosmic shear signal (GG) is shown as a black solid lines; the intrinsic alignment GI term is shown by the black dashed lines; the intrinsic alignment II term is shown by the dotted black line; the usual galaxy clustering signal (gg) is shown by the magenta solid line; the cross correlation between galaxy clustering and lensing magnification (gm) is shown by the magenta dashed line; the lensing magnification correlation functions (mm) are shown by the magenta dotted line. In the lower left panels the solid blue line shows the correlation between lensing shear and galaxy clustering (Gg); the blue dashed line shows the correlation between lensing shear and lensing magnification (gm); the blue dotdashed line shows the correlation between intrinsic alignment and galaxy clustering (Ig or equivalently gI); the red solid line shows the correlation between galaxy clustering and lensing shear (gG), which is equivalent to the blue solid line with redshift bin indices i and j reversed; similarly the red dashed line shows the correlation between lensing magnification and lensing shear (mG), for cases where the magnification occurs at lower redshift than the shear (i < j); finally the dotted line shows the correlation between lensing magnification and intrinsic alignment (mI). 

In the text 
Fig. 2 Left panels: figures of merit as a function of the number of free parameters as a function of wave vector N_{K} in the bias terms For each line type, the upper curve is obtained for a number of free bias parameters as a function of redshift N_{Z} = 2, the lower is for N_{Z} = 4. Right panels: same as on the left, but as a function of N_{Z}, i.e. the number of redshift parameters in the bias terms. The upper curves for each set correspond now to N_{K} = 2 and the lower ones to N_{K} = 4, respectively. Upper panels: figure of merit taking into account the full cosmological parameter space, FoM_{TOT}, see (45). Lower panels: dark energy figure of merit from the Dark Energy Task Force FoM_{DETF}, see (44). Dashed curves correspond to results using galaxy ellipticity correlations (ϵϵ) only, dotted black curves to galaxy number density correlations (nn) only, and solid black curves to results using all correlations (ϵϵ, nn and ϵn). The grey dotted lines show results for nn correlations without imposing cuts in angular frequency. The constant grey line marks the FoM computed for the pure lensing, i.e. GG, signal, assuming intrinsic alignments do not exist. In addition we show the resulting figures of merit when using our most flexible parametrisation with N_{K} = N_{Z} = 7 as filled symbols. Circles correspond to ϵϵ, triangles to nn, and diamonds to all correlations. 

In the text 
Fig. 3 1σcontours for all pairs of cosmological parameters considered, marginalised over all other parameters. We have used a photometric redshift uncertainty parameter value σ_{ph} = 0.05, ten photometric redshift bins for tomography N_{zbin} = 10, and the most flexible intrinsic alignment and bias model considered in this paper, with over two hundred free parameters (N_{K} = 7, and N_{Z} = 7). Orange (light hatched) confidence regions result from using galaxy number density correlations (nn) (excluding the nonlinear regime) only, red (dark hatched) regions use ellipticity correlations (ϵϵ) alone, and blue (filled) regions correspond to using all available information including densityellipticity crosscorrelations. For reference, the contours obtained from a pure lensing signal are shown as black lines. Flat priors on cosmological parameters have been applied. 

In the text 
Fig. 4 Upper panel: the figure of Merit for all cosmological parameter space FoM_{TOT} as a function of the number of photoz bins used for tomography N_{zbin}, shown for ϵϵ (dashed line), nn (dotted line), and all (solid) correlations. The grey line corresponds to results for lensing only (GG). Throughout, N_{K} = N_{Z} = 5 nuisance parameters for the bias terms are used. These results are obtained for the standard set of parameters and σ_{ph} = 0.05. Lower panel: same as above, but in terms of the dark energy figure of merit FoM_{DETF}. 

In the text 
Fig. 5 Upper panel: the difference d_{FoM}, defined in (47), as a function of the photoz dispersion σ_{ph}, shown for ϵϵ correlations (dashed line), all correlations (black solid line), and the lensing only signal (grey solid line). Throughout, nuisance parameters N_{K} = N_{Z} = 5 are used. These results are obtained for the standard set of parameters and N_{zbin} = 10. Lower panel: Same as above, but in terms of the ratio r_{FoM}, given in (46). 

In the text 
Fig. 6 The difference d_{FOM} between the FoM_{TOT} for a given prior value σ_{prior} and the fiducial FoM_{TOT}, obtained for the least stringent, default values of the different priors we apply. Upper panel: effect of tightening the priors on the nuisance parameters in the bias terms for ϵϵ (dashed line), nn (dotted line), and all (solid) correlations. Lower panel: effect of tightening the priors on the slopes of the luminosity function α^{(i)}, and for the prior on the uncertainty of σ_{ph}. Grey curves correspond to priors on α^{(i)}, black curves to those on σ_{ph}. Except for the ϵϵ case the latter remain very close to zero. As above, ϵϵ correlations are shown as dashed lines, nn correlations as dotted lines, and all correlations as solid lines. Note that ϵϵ correlations do not depend on the α^{(i)}. For all curves, the remaining priors are each set to their default values of 50 for the bias term parameters, and 10 for priors on both α^{(i)} and σ_{ph}. 

In the text 
Fig. 7 The FoM_{TOT} for different subsets of correlations that form the nn signal. Throughout, we used N_{zbin} = 10 and σ_{ph} = 0.05. The marginalisation has been performed for N_{K} = N_{Z} = 7 nuisance parameters (black dashed bars) or N_{K} = N_{Z} = 3 nuisance parameters (grey solid bars). The two leftmost columns show the full set of ϵϵ and nn correlations for reference. The two centre columns stand for gg and mm correlations only, employing the cuts in angular frequency. The rightmost columns show again gg and mm, but here the full information up to ℓ_{max} = 3000 is considered. Note that mm correlations depend neither on intrinsic alignments nor on galaxy bias and hence are independent of N_{K} and N_{Z}. 

In the text 
Fig. 8 Upper panel: the FoM_{TOT} for different subsets of correlations that are observable from a galaxy survey with galaxy shape and number density measurements. The Fisher matrix has been computed for the total data set containing all correlations, i.e. these results represent the information content of the subsets as a part of the total signal. Black dashed bars correspond to a marginalisation over N_{K} = N_{Z} = 7 nuisance parameters, grey solid bars to N_{K} = N_{Z} = 3 nuisance parameters. The labels on the abscissa indicate the different combinations of correlations used, where ϵn stands for the crosscorrelations between number density and ellipticity. Lower panel: difference Δ FoM, see (48), for the same subsets as above. This difference can be understood as a measure of the depletion of information in the subsets due to the marginalisation over nuisance parameters. 

In the text 
Fig. C.1 Top panel: slope of the galaxy luminosity function α as a function of redshift, shown for different magnitude limits r_{lim} as indicated in the legend. Bottom panel: fit parameters a_{1} (solid), a_{2} (dotted), and a_{3} (dashed) as a function of r_{lim}. These parameters are obtained from the polynomial fit in (C.2). 

In the text 
Fig. C.2 Percentage deviation of the simplified fit as defined in (C.2) and (C.3) from the original fits described by (C.1) and shown in the upper panel of Fig. C.1. The contour levels correspond to the percentages given on the curves. Negative deviations are indicated by dashed contours. Note that across most of the parameter plane the modulus of the deviation is less than 1%. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.