# Reconstruction of the Primordial Power Spectrum using Temperature and Polarisation Data from Multiple Experiments

###### Abstract

We develop a method to reconstruct the primordial power spectrum, , using both temperature and polarisation data from the joint analysis of a number of Cosmic Microwave Background (CMB) observations. The method is an extension of the Richardson-Lucy algorithm, first applied in this context by Shafieloo & Souradeep Shafieloo and Souradeep (2004). We show how the inclusion of polarisation measurements can decrease the uncertainty in the reconstructed power spectrum. In particular, the polarisation data can constrain oscillations in the spectrum more effectively than total intensity only measurements. We apply the estimator to a compilation of current CMB results. The reconstructed spectrum is consistent with the best–fit power spectrum although we find evidence for a ‘dip’ in the power on scales Mpc. This feature appears to be associated with the WMAP power in the region which is consistently below best–fit models. We also forecast the reconstruction for a simulated, Planck–like pla (2006) survey including sample variance limited polarisation data.

## I Introduction

With increasingly precise measurements being made of the CMB (Cosmic Microwave Background) or LSS (Large Scale Structure) it becomes progressively more important to determine how new observations could yield the greatest insight into processes occurring in the early universe. One such observable is the primordial power spectrum of curvature perturbations ,

(1) |

where is the wavenumber.

A generic prediction of the simplest inflationary models is that the density perturbations should be adiabatic and near scale invariant. In these models the spectrum takes the form of a power law . Current limits on this parametrisation place the spectral index Hinshaw et al. (2008). More complex inflationary models such as those with features on the potential Adams et al. (2001); Wang et al. (2005); Hunt and Sarkar (2004); Joy et al. (2008a); Hunt and Sarkar (2007); Pahud et al. (2008); Lerner and McDonald (2009), a small number of -folds Contaldi et al. (2003); Powell and Kinney (2007); Nicholson and Contaldi (2008), or other exotic inflationary models Lesgourgues (2000); Feng and Zhang (2003); Mathews et al. (2004); Jain et al. (2009); Romano and Sasaki (2008) can modify in a manner not compatible with the simple power law description. There have been many parametric searches for the features produced by these models Bridle et al. (2003); Contaldi et al. (2003); Parkinson et al. (2005); Sinha and Souradeep (2006); Sealfon et al. (2005); Mukherjee and Wang (2005); Bridges et al. (2006a, b); Covi et al. (2006); Joy et al. (2008b); Verde and Peiris (2008), although none have proved conclusive. On the other hand, there have been tantalising hints of anomalous features in the data, for example after the first year WMAP results were released, there were strong indications of a cut-off in on large scales. With subsequent data releases the significance of this feature has been reduced, although future observations of the polarisation of the CMB may provide more conclusive evidence Nicholson and Contaldi (2008).

A more thorough search for features in would be one in which the theoretical model is independent from the reconstruction. Methods such as the Richardson-Lucy deconvolution Shafieloo and Souradeep (2004); Shafieloo et al. (2007); Shafieloo and Souradeep (2008), cosmic inversion Matsumiya et al. (2002); Kogo et al. (2004a, b); Nagata and Yokoyama (2008, 2009), and other non-parametric approaches Hannestad (2001); Wang and Mathews (2002); Bridle et al. (2003); Mukherjee and Wang (2003a, b); Hannestad (2004); Tocchini-Valentini et al. (2005); Leach (2006); Paykari and Jaffe (2009) all attempt to overcome theoretical bias. With current computing power these techniques are generally limited to recovering the spectrum for one set of assumed cosmological parameters. This allows one to use a fiducial CMB photon transfer function to integrate the primordial curvature perturbation into today’s photon distribution perturbation. However this assumption hides a significant degeneracy between features in the primordial power spectrum and the physical parameters which determine the height and position of acoustic peaks in the CMB when using only total intensity data in the reconstruction process Hu and Okamoto (2004). As such it is not clear what the significance of any features found in the reconstructed should be. Adding polarisation information into the reconstruction or inversion will significantly reduce this degeneracy since the response of polarisation CMB spectra is phase-shifted with respect to the total intensity response.

The total intensity ( modes) of the CMB have been measured accurately by several instruments Hinshaw et al. (2008); Reichardt et al. (2008); Pryke et al. (2008); Jones et al. (2006); Pearson et al. (2003) to arcminute scales. Measurements in the gradient (-mode) and curl like (-mode) polarisation components, and their correlation with the total intensity , and , lag behind due to their lower amplitude. However a number of experiments are now measuring -mode polarisation with increasing signal-to-noise, starting with the first detection of spectrum Kovac et al. (2002) and subsequent measurements Pearson et al. (2003); Montroy et al. (2006); Hinshaw et al. (2008); Pryke et al. (2008); Bischoff et al. (2008). -mode polarisation has yet to be measured. Primordial tensor fluctuations may have an impact on the reconstruction of on large scales. However, -modes have not been detected yet and we neglect their contribution in this work.

We show in this paper that when information contained in both total intensity and polarisation radiation transfer functions is used in the reconstruction of the primordial power spectrum tighter constraints can be obtained. The paper is organised as follows; in section II we introduce the extension of the Richardson-Lucy algorithm used to estimate . In Section III.1 we test the algorithm using simple forecasts of CMB data. Our tests include template input models with radically broken scale invariance. We explore current limits on in section III.2 and conclude in section IV.

## Ii An extended Richardson-Lucy Estimator

Direct primordial power spectrum reconstruction requires the inversion of the following relations

(2) |

where and represent , , or -type anisotropies, are the angular power spectra for the combination and the are the photon perturbation transfer functions. The transfer functions are obtained by integrating the full Einstein-Boltzmann system of differential equations. These describe the evolution of perturbations in the photon distribution functions in the presence of gravity and other sources of stress–energy. The functions determine all of the structure in the anisotropy spectra which arises after the initial conditions are set. Most notably the contain distinct peaks due to the acoustic oscillation of the tightly coupled photon-baryon fluid in gravitational potential wells at the time of last scattering. The aim of any inversion method is to distinguish such features from any structure in the initial perturbation spectrum.

For a finite sampling of the wavenumber space Eq. (2) can be recast as an operator acting on the primordial spectrum

(3) |

with operator

(4) |

where are the logarithmic intervals for the chosen sampling.

A solution for cannot be obtained from a direct inversion of the as it is numerically singular. This is due to the high level of degeneracy in the transfer functions relating the power at any wavenumber to angular multipoles although the system can be inverted by binning or smoothing appropriately to reduce the degeneracy Nicholson (2009).

An alternative approach is the iterative inversion employing the Richardson-Lucy (RL) method H. (1972); Lucy (1974) for image reconstruction. The RL method has been widely used in enhancing telescope images Jorissen et al. (2001); Surpi and Blandford (2003); Helder and Vink (2008) and it can be shown that the RL method converges to the maximum likelihood estimator in the case of a Poisson distributed signal Shepp and Vardi (1982). In the following we outline the use of the RL estimator in reconstructing the primordial power spectrum as introduced by Shafieloo et al. Shafieloo and Souradeep (2004) and extend it to include properly weighted contributions from polarisation measurements.

Consider the case where the original source plane is the isotropic, primordial Fourier space spanned by the wavenumber and the convolved image plane is the space of angular multipoles . In this case the convolution filter is which relates modes the source power to the image power .

The RL method provides an iterative solution to (3) for , given an observed with

(5) |

where is the image obtained from the iteration and such that in the limit we have . It is also convenient to recast (5) as change in the Shafieloo and Souradeep (2004)

(6) |

such that the cut-offs in multipole space and are not propagated to the iterated solution through the broadening action of .

As it is (6) does not account for errors in the observed image as it applies uniform weighting to all . Since the RL estimator is not a well defined Maximum Likelihood estimator there is not a single choice of optimal weighting. Observational noise must therefore be included through an empirically chosen weighting function. We chose to add a weighting

(7) |

where is the reported error in the . The weights have the properties that if then and if then . The weighted estimator (6) now becomes

(8) |

Most CMB experiments observe only part of the sky which leads to correlated . The estimator does not account for the correlations as it includes only diagonal estimates of the uncertainty in . Indeed even full-sky experiments such as WMAP have correlated due to galaxy and source cuts. Experiments observing only small fractions of the sky usually report results in as a set of bandpowers , where are the bandpower window functions. The bandpowers are less correlated than individual multipoles and allow us to apply the RL estimator to cut–sky experiments by calculating band filtered transfer function operators

(9) |

Using this we can define a generalised iterative estimator with contributions from any number of bandpowers as

(10) |

where the in the full–sky limit (WMAP) and are the reported bandpower window functions for other experiments. In this work we will be using bandpowers from ACBAR Reichardt et al. (2008), QUaD Pryke et al. (2008), BOOMERanG Jones et al. (2006) and CBI Pearson et al. (2003) which increase the range of scales probed past WMAP’s resolution limits.

A problem faced by the RL estimator when including polarisation measurements is that the polarisation bandpowers are correlated with the total intensity ones and cannot be included as further linear contributions to the sum in (10). Following the standard procedure for extending Maximum Likelihood estimators to polarised bandpowers we treat each multipole measurement as a matrix of bandpowers with

(11) |

and similarly elevated the transfer function operator to a matrix

(12) |

and we have disregarded any -mode contribution since current observational limits do not warrant its inclusion. We can now extend the RL estimator to include the full polarisation information as

where is obtained by modifying (7) to sum over banded quantities and is normalised so that the sum over its diagonal elements is unity.

We have implemented the estimator (II) to reconstruct over a discretely sampled grid in from a combination of CMB measurements listed above and for forecasted future data as described below. The full iteration typically takes minutes on a current desktop CPU. As a convergence criterion we iterate until the fractional change in the solution is less than 0.01. Typically this takes iterations to achieve. We reconstruct over some 3000 points between Mpc and Mpc. The sampling is initially logarithmic and switches to linear at Mpc following the typical setup in Einstein-Boltzmann solvers such as CAMB Lewis et al. (2000) but with higher resolution.

As shown in Shafieloo and Souradeep (2004) the RL method requires a smoothing of the recovered . We have found that applying a smoothing kernel at each iteration improves the convergence when using unbinned such as the WMAP results in the estimate. The smoothing helps to damp down large fluctuations at each iteration driven by the scatter in the observed data, due to either sample or noise variance. The behaviour is due to the estimator not being a well defined maximum likelihood one and thus not taking into account the full correlation structure of the data in its weighting. We have chosen a simple smoothing kernel for the applications reported in this work by taking a nearest neighbour running average over the between each iteration. Given our sampling this defines a smoothing size of Mpc at low and Mpc at high . Other smoothing kernels may lead to more efficient behaviour of the estimator such as faster convergence but would be more computationally demanding and would lead to a further loss of resolution in space.

Since the estimator requires minutes to converge it is not currently possible to employ it as part of a full Markov Chain Monte Carlo parameter search in place of the standard primordial spectrum parameter; amplitude and spectral tilt . This may be possible in future if the algorithm is optimised and/or parallelised but in this work we chose to explore the “source” plane assuming a set of fixed fiducial transfer functions based on the best-fit parameters. The transfer functions are obtained from CAMB with parameters , , , and .

To estimate the confidence limits around the sampled solution we Monte Carlo the iterative estimator by generating 1000 simulated sets of . For simplicity we assume Gaussian distributions for all the since more accurate distributions are not easily reconstructed from some of the published data. We iterate to convergence on each of the simulated data sets and obtain a covariance matrix for the samples by averaging over the ensemble

(14) |

with . We quote 1 and 2- errors by taking multiples of the square root of the diagonal elements of the matrix.

## Iii Application of the extended RL estimator

### iii.1 Tests on forecasted CMB data

We have tested the extended RL estimator on a number of input templates using forecasted data with similar experimental properties as the upcoming Planck satellite mission pla (2006). We assume a total of 12 detectors with NETs of 64 observing 80% of the sky over 12 months with a resolution of 7 arcminutes FWHM. We calculate errors around our fiducial CMB best–fit models in both total intensity and polarisation spectra for this experimental setup and use these together with samples on the fiducial model to test the estimator’s convergence. We consider multipoles of for total intensity spectra and for polarisation. We have not taken into account any residual error from foreground subtraction in our forecasts. Thus our forecast are on the optimistic side of the accuracy achievable, particularly in the polarisation where foreground removal will certainly have a significant impact on errors at .

We use four separate templates to generate the s. The first is a standard power law with but with a lower amplitude than the best–fit, the second is a running model with , the third is a power law with a sharp, compensated feature at Mpc Lesgourgues (2000); Adams et al. (2001); Feng and Zhang (2003); Mathews et al. (2004); Joy et al. (2008a); Hunt and Sarkar (2007); Jain et al. (2009); Lerner and McDonald (2009) and the fourth is a power law with superimposed sinusoidal oscillations Wang et al. (2005); Martin and Ringeval (2004, 2005); Pahud et al. (2008); Romano and Sasaki (2008). The extended estimator is run on all four sets of forecasts and the resulting solution is compared to the input one. In Fig. 1 we show the results for runs including only and runs including , , and .

Each run is started with a first guess (green/long–dashed) which is the current best–fit power law spectrum. In all cases the structure in the input is reproduced to some degree over a large range in . Beyond the range shown in Fig. 1 there are significant departures from the correct solution. On the large scales ( Mpc) this is due to the cut–off in the transfer functions due to the horizon scale. At small scales ( Mpc) the cut–off is due to the resolution limit of the forecasts. We find the method is particularly well suited for reconstructing long wavelength structure such as in the running and oscillating model. Although the sharp feature in the third model is present in the reconstructed spectrum, its amplitude is not recovered accurately. This example highlights the limitations of such methods in reconstructing high frequency features in the primordial spectrum. In the first model we also see high frequency features on small scales, this is due to the difference in amplitude between the initial guess and correct solution. A power law initial guess with closer starting amplitude to the correct value reduces this small scale noise. This could be easily obtained by carrying out a standard power law fit to the data before running the reconstruction estimate.

In Fig. 2 we show the same set of reconstructions but include 1 and 2- confidence regions obtained from the Monte Carlo covariance matrix. We have over sampled the and care should be taken in interpreting the significance of any feature as the samples are highly correlated. The plot includes a set of contours for the reconstruction using only total intensity data (blue/dashed) and including polarisation (red/solid). The range in has been modified to emphasise the regions where the best limits are obtained. The addition of polarisation data gives additional constraints on the smaller scales particularly in the oscillating model case. This is not surprising since total intensity and polarisation measurements are complementary in differentiating any structure in the primordial spectrum from acoustic oscillations imprinted on the CMB at last scattering. The polarisation data has little effect in improving the constraints on the spectrum with a compensated feature although a detection of the feature is evident.

### iii.2 Current Limits

In this section we use the extended RL estimator to reconstruct the primordial power spectrum using currently available data. Our combination of data includes the latest WMAP release Hinshaw et al. (2008) including for , the CBIpol results Pearson et al. (2003) including band powers 2 to 12 and all polarisation band powers, the BOOMERanG 2003 Jones et al. (2006) flight band powers, the QUaD Pryke et al. (2008) results, and the ACBAR 2008 Reichardt et al. (2008) excluding the last two band powers due to the excess power. We have excluded the WMAP polarisation measurements as they are too noisy to be used in the reconstruction would require binning to reduce the scatter.

The overall effect of adding in the sub–orbital experiment is to extend the range of the reconstruction past Mpc where the WMAP noise becomes too large for the data to be included in the estimator. Given WMAP’s resolution in multipole space, it strongly dominates the reconstruction on scales larger than Mpc. In Fig. 3 we show the results of the reconstruction. There is no indication of a departure from a pure power law for the recovered spectrum except for a dip at Mpc.

To determine the origin of the feature at Mpc we re-estimate the power spectrum with a cutout of the WMAP data in the range where the lie significantly below the best–fit models. Fig. 4 shows the result of for the reconstruction. No evidence for the feature remains when the WMAP data is removed from the estimate indicating that the dip is associated with the feature in the .

## Iv Discussion

The RL estimator provides a model independent method for reconstruction of the primordial power spectrum of perturbations from measurements of the CMB angular power spectrum . We have extended the RL estimator to use on multiple data sets including properly weighted polarisation data and have used Monte Carlo realisations of the input measurements to estimate confidence limits around the reconstructed spectra.

We have applied the new estimator to current measurements.These include band powers from the QUaD experiment which provide the highest signal–to–noise measurements of power so far. Including the polarisation information increases the constraints on the power spectrum reconstruction as it carries independent information which is complementary in phase to that of the total intensity. However current and measurements are still noise dominated and do not contribute significantly but future sample variance limited measurements will help to constrain any structure in the primordial spectrum as shown in our examples using forecasted data. An exception to this is the current WMAP polarisation data which may have some impact if it were binned to reduce noise scatter. We have not explored this option in this work as this would have required binning of the data too with a resulting loss of resolution in space.

We have found that the addition of polarisation data is particularly helpful in constraining oscillatory structure in . Spectra with superimposed sinusoidal features have been considered in the literature and have been constrained using model dependent fits Martin and Ringeval (2005); Hamann et al. (2007); Pahud et al. (2008); Liu et al. (2009).

Due to the empirical nature of the weighting used the estimator has a limited acceptance range in the signal–to–noise of the band powers it is run on. Low signal–to–noise or cause the estimator to converge very slowly. Conversely if the data is weighted too strongly there is a loss of resolution and signal–to–noise in the reconstructed . The best approach is to exclude the low signal–to–noise measurements from the data being used. It is possible that other weighting schemes may be more efficient in allowing the whole range of measurements to be used. However, in its present form, the estimator will be most useful when all polarisations will be sample variance limited out to the same . As an example of the improvement this will have over current estimates we add our forecasted and data to the current measurements and run the extended RL estimator. The result is shown in Fig. 5 which shows an overall improvement in the Monte Carlo confidence limits over the entire range in being reconstructed.

In this work we have not explored the degeneracy of the method with the physical parameters determining the structure of the CMB transfer functions. The current run–time to convergence for each reconstruction is too long to allow the tens of thousands of runs required to implement the method as part of a full Markov Chain Monte Carlo exploration of the entire parameter space. However with parallelisation and more efficient convergence it may become possible to do this in the near future.

There are further extensions of the RL estimator that could increase its effectiveness. The addition of other observables with different transfer functions such as galaxy redshift surveys or cosmic shear surveys will provide complementary information in the reconstruction. Our extension of the RL estimator is compatible with these other observables as they can be included as further independent data points with properly formatted and binned transfer functions. Such extensions would contribute key information in the reconstruction on scales where the presence of foregrounds in the CMB reduces its effectiveness.

###### Acknowledgements.

This work was supported by a STFC studentship and used the Imperial College high performance computing service^{1}

^{1}1http://www.imperial.ac.uk/ict/services/teachingandresearchservices/highperformancecomputing.

## References

- Shafieloo and Souradeep (2004) A. Shafieloo and T. Souradeep, Phys. Rev. D70, 043523 (2004), eprint astro-ph/0312174.
- pla (2006) Planck: The scientific programme (2006), eprint astro-ph/0604069.
- Hinshaw et al. (2008) G. Hinshaw et al. (WMAP) (2008), eprint astro-ph/0803.0732.
- Adams et al. (2001) J. A. Adams, B. Cresswell, and R. Easther, Phys. Rev. D64, 123514 (2001), eprint astro-ph/0102236.
- Wang et al. (2005) X. Wang, B. Feng, M. Li, X.-L. Chen, and X. Zhang, Int. J. Mod. Phys. D14, 1347 (2005), eprint astro-ph/0209242.
- Hunt and Sarkar (2004) P. Hunt and S. Sarkar, Phys. Rev. D70, 103518 (2004), eprint astro-ph/0408138.
- Joy et al. (2008a) M. Joy, V. Sahni, and A. A. Starobinsky, Phys. Rev. D77, 023514 (2008a), eprint 0711.1585.
- Hunt and Sarkar (2007) P. Hunt and S. Sarkar, Phys. Rev. D76, 123504 (2007), eprint 0706.2443.
- Pahud et al. (2008) C. Pahud, M. Kamionkowski, and A. R. Liddle (2008), eprint 0807.0322.
- Lerner and McDonald (2009) R. Lerner and J. McDonald, Phys. Rev. D79, 023511 (2009), eprint 0811.1933.
- Contaldi et al. (2003) C. R. Contaldi, M. Peloso, L. Kofman, and A. Linde, JCAP 0307, 002 (2003), eprint astro-ph/0303636.
- Powell and Kinney (2007) B. A. Powell and W. H. Kinney, Phys. Rev. D76, 063512 (2007), eprint astro-ph/0612006.
- Nicholson and Contaldi (2008) G. Nicholson and C. R. Contaldi, JCAP 0801, 002 (2008), eprint astro-ph/0701783.
- Lesgourgues (2000) J. Lesgourgues, Nucl. Phys. B582, 593 (2000), eprint hep-ph/9911447.
- Feng and Zhang (2003) B. Feng and X. Zhang, Phys. Lett. B570, 145 (2003), eprint astro-ph/0305020.
- Mathews et al. (2004) G. J. Mathews, D. J. H. Chung, K. Ichiki, T. Kajino, and M. Orito, Phys. Rev. D70, 083505 (2004), eprint astro-ph/0406046.
- Jain et al. (2009) R. K. Jain, P. Chingangbam, J.-O. Gong, L. Sriramkumar, and T. Souradeep, JCAP 0901, 009 (2009), eprint 0809.3915.
- Romano and Sasaki (2008) A. E. Romano and M. Sasaki, Phys. Rev. D78, 103522 (2008), eprint 0809.5142.
- Bridle et al. (2003) S. L. Bridle, A. M. Lewis, J. Weller, and G. Efstathiou, Mon. Not. Roy. Astron. Soc. 342, L72 (2003), eprint astro-ph/0302306.
- Parkinson et al. (2005) D. Parkinson, S. Tsujikawa, B. A. Bassett, and L. Amendola, Phys. Rev. D71, 063524 (2005), eprint astro-ph/0409071.
- Sinha and Souradeep (2006) R. Sinha and T. Souradeep, Phys. Rev. D74, 043518 (2006), eprint astro-ph/0511808.
- Sealfon et al. (2005) C. Sealfon, L. Verde, and R. Jimenez, Phys. Rev. D72, 103520 (2005), eprint astro-ph/0506707.
- Mukherjee and Wang (2005) P. Mukherjee and Y. Wang, JCAP 0512, 007 (2005), eprint astro-ph/0502136.
- Bridges et al. (2006a) M. Bridges, A. N. Lasenby, and M. P. Hobson, Mon. Not. Roy. Astron. Soc. 369, 1123 (2006a), eprint astro-ph/0511573.
- Bridges et al. (2006b) M. Bridges, A. N. Lasenby, and M. P. Hobson (2006b), eprint astro-ph/0607404.
- Covi et al. (2006) L. Covi, J. Hamann, A. Melchiorri, A. Slosar, and I. Sorbera, Phys. Rev. D74, 083509 (2006), eprint astro-ph/0606452.
- Joy et al. (2008b) M. Joy, A. Shafieloo, V. Sahni, and A. A. Starobinsky (2008b), eprint 0807.3334.
- Verde and Peiris (2008) L. Verde and H. V. Peiris, JCAP 0807, 009 (2008), eprint 0802.1219.
- Shafieloo et al. (2007) A. Shafieloo, T. Souradeep, P. Manimaran, P. K. Panigrahi, and R. Rangarajan, Phys. Rev. D75, 123502 (2007), eprint astro-ph/0611352.
- Shafieloo and Souradeep (2008) A. Shafieloo and T. Souradeep, Phys. Rev. D78, 023511 (2008), eprint 0709.1944.
- Matsumiya et al. (2002) M. Matsumiya, M. Sasaki, and J. Yokoyama, Phys. Rev. D65, 083007 (2002), eprint astro-ph/0111549.
- Kogo et al. (2004a) N. Kogo, M. Matsumiya, M. Sasaki, and J. Yokoyama, Astrophys. J. 607, 32 (2004a), eprint astro-ph/0309662.
- Kogo et al. (2004b) N. Kogo, M. Sasaki, and J. Yokoyama, Phys. Rev. D70, 103001 (2004b), eprint astro-ph/0409052.
- Nagata and Yokoyama (2008) R. Nagata and J. Yokoyama, Phys. Rev. D78, 123002 (2008), eprint 0809.4537.
- Nagata and Yokoyama (2009) R. Nagata and J. Yokoyama, Phys. Rev. D79, 043010 (2009), eprint 0812.4585.
- Hannestad (2001) S. Hannestad, Phys. Rev. D63, 043009 (2001), eprint astro-ph/0009296.
- Wang and Mathews (2002) Y. Wang and G. Mathews, Astrophys. J. 573, 1 (2002), eprint astro-ph/0011351.
- Mukherjee and Wang (2003a) P. Mukherjee and Y. Wang, Astrophys. J. 593, 38 (2003a), eprint astro-ph/0301058.
- Mukherjee and Wang (2003b) P. Mukherjee and Y. Wang, Astrophys. J. 599, 1 (2003b), eprint astro-ph/0303211.
- Hannestad (2004) S. Hannestad, JCAP 0404, 002 (2004), eprint astro-ph/0311491.
- Tocchini-Valentini et al. (2005) D. Tocchini-Valentini, M. Douspis, and J. Silk, Mon. Not. Roy. Astron. Soc. 359, 31 (2005), eprint astro-ph/0402583.
- Leach (2006) S. M. Leach, Mon. Not. Roy. Astron. Soc. 372, 646 (2006), eprint astro-ph/0506390.
- Paykari and Jaffe (2009) P. Paykari and A. H. Jaffe (2009), eprint astro-ph/0902.4399.
- Hu and Okamoto (2004) W. Hu and T. Okamoto, Phys. Rev. D69, 043004 (2004), eprint astro-ph/0308049.
- Reichardt et al. (2008) C. L. Reichardt et al. (2008), eprint astro-ph/0801.1491.
- Pryke et al. (2008) C. Pryke et al. (QUaD) (2008), eprint astro-ph/0805.1944.
- Jones et al. (2006) W. C. Jones et al., Astrophys. J. 647, 823 (2006), eprint astro-ph/0507494.
- Pearson et al. (2003) T. J. Pearson et al., Astrophys. J. 591, 556 (2003), eprint astro-ph/0205388.
- Kovac et al. (2002) J. Kovac et al., Nature 420, 772 (2002), eprint astro-ph/0209478.
- Montroy et al. (2006) T. E. Montroy et al., Astrophys. J. 647, 813 (2006), eprint astro-ph/0507514.
- Bischoff et al. (2008) . C. Bischoff et al. (CAPMAP) (2008), eprint 0802.0888.
- Nicholson (2009) G. Nicholson (2009), eprint in preparation.
- H. (1972) R. W. H., J. Opt. Soc. Am. 62, 55 (1972).
- Lucy (1974) L. B. Lucy, Astron. J. 79, 745 (1974).
- Jorissen et al. (2001) A. Jorissen, M. Mayor, and S. Udry (2001), eprint astro-ph/0105301.
- Surpi and Blandford (2003) G. Surpi and R. D. Blandford, Astrophys. J. 584, 100 (2003), eprint astro-ph/0111160.
- Helder and Vink (2008) E. A. Helder and J. Vink (2008), eprint 0806.3748.
- Shepp and Vardi (1982) L. A. Shepp and Y. Vardi, Medical Imaging, IEEE Transactions on 1, 113 (1982), ISSN 0278-0062.
- Lewis et al. (2000) A. Lewis, A. Challinor, and A. Lasenby, Astrophys. J. 538, 473 (2000), eprint astro-ph/9911177.
- Martin and Ringeval (2004) J. Martin and C. Ringeval, Phys. Rev. D69, 083515 (2004), eprint astro-ph/0310382.
- Martin and Ringeval (2005) J. Martin and C. Ringeval, JCAP 0501, 007 (2005), eprint hep-ph/0405249.
- Hamann et al. (2007) J. Hamann, L. Covi, A. Melchiorri, and A. Slosar, Phys. Rev. D76, 023503 (2007), eprint astro-ph/0701380.
- Liu et al. (2009) J. Liu, H. Li, J. Xia, and X. Zhang (2009), eprint 0901.2033.