Model-independent constraints on dark energy density from flux-averaging analysis of type Ia supernova data

We reconstruct the dark energy density rho(X)(z) as a free function from current Type Ia supernova (SN Ia) data, together with the cosmic microwave background (CMB) shift parameter from CMB data from the Wilkinson Microwave Anisotropy Probe (WMAP), Cosmic Background Imager ( CBI), and Arcminute Cosmology Bolometer Array Receiver (ACBAR), and the large-scale structure (LSS) growth factor from the Two-Degree Field (2dF) galaxy survey data. We parameterize rho(X)(z) as a continuous function, given by interpolating its amplitudes at equally spaced z-values in the redshift range covered by SN Ia data, and a constant at larger z [where rho(X)(z) is only weakly constrained by CMB data]. We assume a flat universe and use the Markov Chain Monte Carlo (MCMC) technique in our analysis. We find that the dark energy density rho(X)(z) is constant for 0 less than or similar to z less than or similar to 0.5 and increases with redshift z for 0.5 less than or similar to z less than or similar to 1 at a 68.3% confidence level, but is consistent with a constant at a 95% confidence level. For comparison, we also give constraints on a constant equation of state for the dark energy. Flux averaging of SN Ia data is required to yield cosmological parameter constraints that are free of the bias induced by weak gravitational lensing. We set up a consistent framework for flux-averaging analysis of SN Ia data, based on the work of Wang. We find that flux averaging of SN Ia data leads to slightly lower Omega(m) and smaller time variation in rho(X)(z). This suggests that a significant increase in the number of SNe Ia from deep SN surveys on a dedicated telescope is needed to place a robust constraint on the time dependence of the dark energy density.