Direct wavelet expansion of the primordial power spectrum

In order to constrain and possibly detect unusual physics during inflation, we allow the power spectrum of primordial matter density fluctuations, P-in(k), to be an arbitrary function in the estimation of cosmological parameters from data. The multiresolution and good localization properties of orthogonal wavelets make them suitable for detecting features in P-in(k). We expand P-in(k) directly in wavelet basis functions. The likelihood of the data is thus a function of the wavelet coefficients of P-in(k), as well as the Hubble constant H-0, baryon density Omega(b) h(2), cold dark matter density Omega(c) h(2), and the reionization optical depth tau(ri) in a flat LambdaCDM cosmology. We derive constraints on these parameters from cosmic microwave background anisotropy data (WMAP, CBI, and ACBAR) and large-scale structure data (2dFGRS and PSCZ) using the Markov chain Monte Carlo (MCMC) technique. The direct wavelet expansion method is different from and complementary to the wavelet band power method of Mukherjee & Wang, and results from the two methods are consistent. In addition, as we demonstrate, the direct wavelet expansion method has the advantage that once the wavelet coefficients have been constrained, the reconstruction of P-in(k) can be effectively denoised, i.e., P-in(k) can be reconstructed using only the coefficients that, say, deviate from zero at greater than 1 sigma. In doing so, we retain the essential properties of P-in(k). The reconstruction also suffers much less from the correlated errors of binning methods. The shape of the primordial power spectrum, as reconstructed in detail here, reveals an interesting new feature at 0.001 less than or similar to k/Mpc(-1) less than or similar to 0.005. It will be interesting to see whether this feature is confirmed by future data. The reconstructed and denoised P-in(k) is favored over the scale-invariant and power-law forms at greater than or similar to 1 sigma.