[Abridged] We use our most recent training set for the RICO code to estimate the impact of recombination uncertainties on the posterior probability distributions which will be obtained from future CMB experiments, and in particular the Planck satellite. Using a MCMC analysis to sample the posterior distribution of the cosmological parameters, we find that Planck will have biases of -0.7, -0.3 and -0.4 sigmas for n_S, Omega_b h2 and log(As), respectively, in the minimal 6-parameter LCDM model, if the description of the recombination history given by RICO is not used. The remaining parameters are not significantly affected. We also show, that the cosmology dependence of the corrections to the recombination history modeled with RICO has a negligible impact on the posterior distributions obtained for the case of the Planck satellite. In practice, this implies that the inclusion of additional corrections to existing recombination codes can be achieved using simple cosmology-independent `fudge functions'. Finally, we also investigated the impact of some recent improvements in the treatment of hydrogen recombination which are still not included in the current version of our training set for Rico, by assuming that the cosmology dependence of those corrections can be neglected. In summary, with our current understanding of the complete recombination process, the expected biases in the cosmological parameters inferred from Planck might be as large as -2.3, -1.7 and -1 sigmas for n_S, Omega_b h2 and log(As) respectively, if all those corrections are not taken into account. We note that although the list of physical processes that could be of importance for Planck seems to be nearly complete, still some effort has to be put in the validation of the results obtained by the different groups.

We present Rico, a code designed to compute the ionization fraction of the Universe during the epoch of hydrogen and helium recombination with an unprecedented combination of speed and accuracy. This is accomplished by training the machine learning code Pico on the calculations of a multi-level cosmological recombination code which self-consistently includes several physical processes that were neglected previously. After training, Rico is used to fit the free electron fraction as a function of the cosmological parameters. While, for example at low redshifts (z<~900), much of the net change in the ionization fraction can be captured by lowering the hydrogen fudge factor in Recfast by about 3%, Rico provides a means of effectively using the accurate ionization history of the full recombination code in the standard cosmological parameter estimation framework without the need to add new or refined fudge factors or functions to a simple recombination model. Within the new approach presented here it is easy to update Rico whenever a more accurate full recombination code becomes available. Once trained, Rico computes the cosmological ionization history with negligible fitting error in ~10 milliseconds, a speed-up of at least 10^6 over the full recombination code that was used here. Also Rico is able to reproduce the ionization history of the full code to a level well below 0.1%, thereby ensuring that the theoretical power spectra of CMB fluctuations can be computed to sufficient accuracy and speed for analysis from upcoming CMB experiments like Planck. Furthermore it will enable cross-checking different recombination codes across cosmological parameter space, a comparison that will be very important in order to assure the accurate interpretation of future cosmic microwave background data.

This paper presents the second release of Pico (Parameters for the Impatient COsmologist). Pico is a general purpose machine learning code which we have applied to computing the CMB power spectra and the WMAP likelihood. For this release, we have made improvements to the algorithm as well as the data sets used to train Pico, leading to a significant improvement in accuracy. For the 9 parameter nonflat case presented here Pico can on average compute the TT, TE and EE spectra to better than 1% of cosmic standard deviation for nearly all $\ell$ values over a large region of parameter space. Performing a cosmological parameter analysis of current CMB and large scale structure data, we show that these power spectra give very accurate 1 and 2 dimensional parameter posteriors. We have extended Pico to allow computation of the tensor power spectrum and the matter transfer function. Pico runs about 1500 times faster than CAMB at the default accuracy and about 250,000 times faster at high accuracy. Training Pico can be done using massively parallel computing resources, including distributed computing projects such as Cosmology@Home. On the homepage for Pico, located at http://cosmos.astro.uiuc.edu/pico, we provide new sets of regression coefficients and make the training code available for public use.

We present a fast, accurate, robust and flexible method of accelerating parameter estimation. This algorithm, called Pico, can compute the CMB power spectrum and matter transfer function as well as any computationally expensive likelihoods in a few milliseconds. By removing these bottlenecks from parameter estimation codes, Pico decreases their computational time by 1 or 2 orders of magnitude. Pico has several important properties. First, it is extremely fast and accurate over a large volume of parameter space. Furthermore, its accuracy can continue to be improved by using a larger training set. This method is generalizable to an arbitrary number of cosmological parameters and to any range of l-values in multipole space. Pico is approximately 3000 times faster than CAMB for flat models, and approximately 2000 times faster then the WMAP 3 year likelihood code. In this paper, we demonstrate that using Pico to compute power spectra and likelihoods produces parameter posteriors that are very similar to those using CAMB and the official WMAP3 code, but in only a fraction of the time. Pico and an interface to CosmoMC are made publicly available at http://www.astro.uiuc.edu/~bwandelt/pico/.