# Dark Energy and the CMB

\[

## Abstract

We find that current Cosmic Microwave Background (CMB) anisotropy data strongly constrain the mean spatial curvature of the Universe to be near zero, or, equivalently, the total energy density to be near critical—as predicted by inflation. This result is robust to editing of data sets, and variation of other cosmological parameters (totaling seven, including a cosmological constant). Other lines of argument indicate that the energy density of non-relativistic matter is much less than critical. Together, these results are evidence, independent of supernovae data, for dark energy in the Universe.

\] Introduction. Cosmologists have long realized that there is more to the Universe than meets the eye. A wide variety of evidence points to the existence of dark matter in the Universe, matter which cannot be seen, but which can be indirectly detected by its contribution to the gravitational field. As observations have improved, the phenomenology of the “dark sector” has become richer. While dark matter was originally posited to explain what would otherwise be excessively attractive gravity, dark energy explains the accelerating expansion—an apparently repulsive gravitational effect. The most well-known argument for this additional dark component is based on inferences of the luminosity distances to high–$`z`$ supernovae. The anomalously large distances indicate that the Universe was expanding more slowly in the past than it is now; i.e., the expansion rate is accelerating. Acceleration only occurs if the bulk pressure is negative, and this could only be due to a previously undetected component.

Here we argue for dark energy based on another gravitational effect: its influence on the mean spatial curvature. This argument  does not rely on the supernovae observations and therefore avoids the systematic uncertainties in the inferred luminosity distances. It is based on a lower limit to the total density, and a smaller upper limit on the density of non-relativistic matter. The lower limit comes from measurements of the anisotropy of the cosmic microwave background (CMB) whose statistical properties depend on the mean spatial curvature, which in turn depends on the mean total density. We find that the CMB strongly indicates that $`\mathrm{\Omega }>0.4`$, where $`\mathrm{\Omega }`$ is the ratio of the total mean density to the critical density (that for which the mean curvature would be zero). Upper limits to the density of non-relativistic matter come from a variety of sources which quite firmly indicate $`\mathrm{\Omega }_m<0.4`$.

The CMB sensitivity to curvature is due to the dependence on curvature of the angular extent of objects of known size, at known redshifts. CMB photons that are penetrating our galaxy today, were emitted from a thin shell at a redshift of $`z1100`$ (called the “last-scattering surface”) during the transition from an ionized plasma to a neutral medium. The “object” of known size at known redshift is the sound-horizon of the plasma at the epoch of last-scattering. Its observational signature is the location of a series of peaks in the angular power spectrum of the CMB.

One must be careful about using current CMB data to determine $`\mathrm{\Omega }`$ or any other cosmological parameters for several reasons. First, these are very difficult experiments, and the data sets they produce have low signal-to-noise ratios and limited frequency ranges, complicating the detection of systematic errors. Use of different calibration standards further increases the risk of underestimated systematic error. To counter these problems, we examine the robustness of our results to editing of data sets, and check that the distribution of model residuals is consistent with the stated measurement uncertainties.

Second, the CMB angular power spectra depend on a number of parameters other than the curvature. To some degree, a change in curvature can be mimicked by changes in other parameters. We therefore vary six parameters besides the curvature, placing mild prior constraints on some of these so as not to explore unrealistic regions of the parameter space.

Finally, existing data are insufficient to firmly establish the paradigm for structure formation which we have assumed: structure grew via gravitational instability from primordial adiabatic perturbations. Our conclusions depend on this assumption. At present, this counts as a possible source of systematic error. Fortunately, future CMB data will verify (or refute) the paradigm and will also allow for the determination of $`\mathrm{\Omega }`$ with greatly reduced model dependence.

The data. Present data are already so abundant that it must be compressed before it can serve as the basis for a multi-dimensional parameter search. Fortunately, all data sets have been compressed to constraints on the angular power spectrum, $`C_l2\pi C(\theta )P_l(\mathrm{cos}\theta )d(\mathrm{cos}\theta )`$ where $`C(\theta )`$ is the correlation function. Because of the tremendous reduction in the size of the data sets, this data compression is called “radical compression” .

Here we use the radically compressed data from http://www.cita.utoronto.ca/k̃nox/radical.html. In this compilation the non-Gaussianity of the power spectrum uncertainties has been characterized for a number of experiments with a lot of the weight (including all those plotted with large symbols in Fig. 1); assuming Gaussianity leads to biases.

The Search Method. We search over a seven-dimensional parameter space specified by $`\mathrm{\Omega }`$, $`\mathrm{\Omega }_bh^2`$, $`\mathrm{\Omega }_{\mathrm{cdm}}h^2`$, $`\mathrm{\Omega }_\mathrm{\Lambda }h^2`$, $`\tau `$, $`n_s`$ and $`C_{10}`$, where $`\mathrm{\Omega }_i=\rho _i/\rho _c`$ and $`i=b,\mathrm{cdm},\mathrm{\Lambda }`$ is for baryons, cold dark matter and a cosmological constant respectively, $`\rho _c3H_0^2/(8\pi G)`$ is the critical density, $`\tau `$ is the optical depth to Thomson scattering, $`n_s`$ is the power-law index of the primordial matter power spectrum, and $`C_{10}`$ serves as the normalization parameter. The Hubble constant, $`H_0100h\mathrm{km}\mathrm{sec}^1\mathrm{Mpc}^1`$, is a dependent variable in this space, due to the sum rule: $`\mathrm{\Omega }_\mathrm{\Lambda }+\mathrm{\Omega }_b+\mathrm{\Omega }_{\mathrm{cdm}}=\mathrm{\Omega }`$. Note that, for specificity and simplicity, we have chosen the dark energy to be a cosmological constant; other choices (e.g., qunitessence) would not significantly affect our curvature constraints.

For each value of $`\mathrm{\Omega }`$ we vary the 23 other parameters (six cosmological and 17 calibration—one for each experiment) to find the minimum value of $`\chi ^2=\chi _d^2+\chi _p^2`$. Here $`\chi _d^2`$ is the offset log-normal form explicitly given in Eq. 39-43 of , which was shown to be a good approximation to the log of the likelihood function. Information from non-CMB observations is included as a prior contribution, $`\chi _p^2`$. Unless otherwise stated, we assume that $`h=0.65\pm 0.1`$ (a reasonable interpretation of several measurements) and $`\mathrm{\Omega }_bh^2=0.019\pm 0.003`$ (from  but with a 40% increase in their uncertainty). We use the Levenberg-Marquardt method to find the minimum value of $`\chi ^2`$ for each value of $`\mathrm{\Omega }`$. We stop the hunt when the new $`\chi ^2`$ is within $`0.1`$ of the old value. We tested this method on simulated data and recovered the correct results.

The likelihood of the best-fit model, $`(\mathrm{\Omega })`$ is proportional to $`\mathrm{exp}(\chi ^2/2)`$. Ideally we would marginalize over the non-$`\mathrm{\Omega }`$ parameters rather than maximizing over them. However, we note that in the limit that the likelihood is Gaussian, these two procedures are equivalent. More generally, in order for marginalization to give qualitatively different answers there would have to be, with decreasing $`\mathrm{\Omega }`$, a very rapid increase in the volume of parameter space in the non-$`\mathrm{\Omega }`$ direction with $`\chi ^2`$’s comparable to the minimum $`\chi ^2`$. Inspection of the Fisher matrix leads us to believe this is not the case.

The Results. Our main results are shown in Fig. 2: the relative likelihood ($`\mathrm{exp}(\chi ^2/2)`$) of the different values of $`\mathrm{\Omega }`$. Including all the data, the best-fit (minimum $`\chi ^2`$) $`\mathrm{\Omega }=1`$ model is $`2\times 10^7`$ times more probable than the best-fit $`\mathrm{\Omega }=0.4`$ model. $`\mathrm{\Omega }<0.7`$ is ruled out at the $`95\%`$ confidence level.

To test the robustness of this result, we edited out single data sets suspected of providing the most weight. Most of these editings produced little change. Only the omission of TOCO changes things substantially, and even then, the best-fit $`\mathrm{\Omega }=1`$ model is $`150`$ times more probable than the best-fit $`\mathrm{\Omega }=0.4`$ model. We also edited pairs of data sets: for no CAT and TOCO, no MSAM and CAT, and no MSAM and TOCO, we find $`\mathrm{\Omega }=1`$ to be $`120`$, $`2.5\times 10^6`$ and $`8`$ times more likely than $`\mathrm{\Omega }=0.4`$. Also shown, as measures of goodness-of-fit, are $`\chi ^2`$ and the degrees of freedom. The $`\chi ^2`$ value for the “All” case is a bit high, but one expects even higher ones over $`8\%`$ of the time, so there is no strong evidence for inconsistencies in the data. As further indication of the robustness of the result, one can see from the “TOCO” panel of Fig. 2 that it persists even when all but a single data set is removed.

For the “All” case, the best-fit $`\mathrm{\Omega }=1`$ model has $`\mathrm{\Omega }_bh^2=0.019`$, $`h=0.65`$, $`\mathrm{\Omega }_\mathrm{\Lambda }=0.69`$, $`\tau =0.17`$ and $`n=1.12`$ and is plotted in Fig. 1. There are degeneracies among these parameters though and none of them is strongly constrained on its own. For example, an equivalently good fit (to just the CMB data) is given by the following model with no tilt or reionization: $`\mathrm{\Omega }=1`$, $`\mathrm{\Omega }_bh^2=0.021`$, $`h=0.65`$, $`\mathrm{\Omega }_\mathrm{\Lambda }=0.65`$, $`\tau =0`$ and $`n=1`$.

We also covered the $`\mathrm{\Omega }_m`$, $`\mathrm{\Omega }_\mathrm{\Lambda }`$ plane, at each point finding the minimum $`\chi ^2`$ possible with variation of the remaining 5 parameters. Figure 3 shows the $`\mathrm{\Delta }\chi ^2=`$ 1, 4 and 9 contours in this plane (which, for a Gaussian, correspond to the 40%, 87% and 99% confidence regions, respectively).

Discussion. Figure 3 also shows constraints on $`\mathrm{\Omega }_m`$ from clusters. Although constraints on $`\mathrm{\Omega }_m`$ arise from a variety of techniques (for reviews see ) perhaps the most reliable are those based on the determination of the ratio of baryonic matter to dark matter in clusters of galaxies . With the assumption that the cluster ratio is the mean ratio (reasonable due to the large size of the clusters), and the baryonic mean density from nucleosynthesis, one can constrain the range of allowable values of $`\mathrm{\Omega }_m`$. Since only the baryonic intracluster gas is detected, the upper limits on $`\mathrm{\Omega }_m`$ from this method are better understood than the lower limits. Mohr et al.  find, from a sample of 27 X-ray clusters, that (including corrections for clumping and depletion of the gas) $`\mathrm{\Omega }_m<(0.32\pm 0.03)/\sqrt{h/0.65}`$. Including the Hubble constant uncertainty ($`h=0.65\pm 0.1`$) this becomes $`\mathrm{\Omega }_m<0.32\pm 0.05`$. Assuming 10% of the baryons to be in galaxies as opposed to the gas, as estimated by , we find $`\mathrm{\Omega }_m=0.29\pm 0.05`$. Results from observations of the Sunyaev-Zeldovich effect in clusters are consistent, though less restrictive: $`\mathrm{\Omega }_m=0.31\pm 0.1`$ . Most other methods (those that do not rely on the cluster baryon fraction) generally result in formally stronger upper limits to $`\mathrm{\Omega }_m`$. This increases our confidence in the Mohr et al. $`\mathrm{\Omega }_m`$ upper limit, but we do not quote these stronger constraints due to our concerns that they are affected by systematic uncertainties that are more difficult to quantify than those in the baryon fraction method.

There have been a number of other analyses of CMB anisotropy data which generally obtained weaker constraints on $`\mathrm{\Omega }`$. There are technical differences between our work and previous work: we account for the non-Gaussianity of the likelihood function, allow for calibration uncertainties, place “sanity” priors on the Hubble constant and the baryon density, and vary six parameters in addition to the curvature. Also, much of the strength of our argument comes from data reported within the last year.

The verdict from the CMB is now in. It does not depend on any one, or even any two, experiments. It clearly points towards a flat Universe and, together with cluster data, strongly indicates the existence of dark energy. These conclusions are consistent with, and independent of, the supernovae results. The completely different sets of systematic uncertainties in the two arguments further strengthen the case. Other constraints in the $`\mathrm{\Omega }_m,\mathrm{\Omega }_\mathrm{\Lambda }`$ plane were recently obtained  by combining cosmic flow data with supernovae observations.

We have neglected several data sets, all of which, if included, would only strengthen our conclusions. Two of these are PythonV and Viper. PythonV and Viper together trace out a peak with centroid near $`l=200`$, and a significant drop in power by $`l=400`$. They have not been included because of the strong correlations in the existing reductions of the data; a new reduction with all correlations specified will soon be available for PythonV.

Any model without a drop in power from $`l=200`$ to $`l=400`$ has difficulties agreeing with all the data. Models fitting this description include the adiabatic models considered here with $`\mathrm{\Omega }<0.4`$ and also topological defect models, whose breadth is a consequence of the loss of the coherent peak structure .

We have been concentrating on implications of the peak location, but the height is also of interest. With fixed $`h`$, it is additional evidence for low $`\mathrm{\Omega }_m`$. The lower $`\mathrm{\Omega }_mh^2`$, the later the transition from a radiation-dominated Universe to a matter-dominated Universe and the larger the early ISW effect, which contributes in the region of the first peak . For flat models, the best fit is at $`\mathrm{\Omega }_m=0.4`$ with $`\mathrm{\Omega }_m=1`$ four times less likely.

Conclusions. We have shown that $`\mathrm{\Omega }=1`$ is strongly favored over $`\mathrm{\Omega }=0.4`$. This result is interesting for two reasons. First, $`\mathrm{\Omega }=1`$ is a prediction of the simplest models of inflation. Second, together with the constraint $`\mathrm{\Omega }_m<0.4`$, it is evidence for dark energy.

The CMB can say little about the nature of the dark energy. A cosmological constant fits the current data, but then so would many of the other forms of dark energy proposed over the past few years. Generation and exploration of new theoretical ideas as to the nature of this dark energy is clearly warranted.

Measurements of CMB anisotropy have already delivered on their promise to provide new clues towards an improved understanding of cosmological structure formation and fundamental physics. We look forward to greater clarification of the dark energy problem, as well as possibly new surprises, from improved CMB anisotropy measurements in the near future.

###### Acknowledgements.

We are grateful to J. Mohr for useful conversations and A. Jaffe for supplyings us with the SN data. We used CMBFAST many, many times. SD is supported by the DOE and by NASA Grant NAG 5-7092. LK is supported by the DOE, NASA grant NAG5-7986 and NSF grant OPP-8920223.