Maximum likelihood search results without 85GHz and 145GHz

  B. Racine

In this posting we run the maximum likelihood search on the DC4, removing the 85GHz and 145GHz channels.
The goal is to assess some benefits of the split-band strategy. Since we did not create a new dataset redistribution the effort in joint bands, we will have a decrease in sensitivity. This posting is mostly focusing on the possible bias increase when dropping the split band.


This posting is a modification of the analysis of CMB-S4 Data Challenge 04 using a BICEP/Keck-style parametric multicomponent analysis, reported in this posting, where we also described the parameterization, studied the dependence on priors, presented the simulations etc.

In the DC4 challenge, we used 9 bands (20, 30, 40, 85, 95, 145, 155, 220 and 270 GHz), specified here, and plotted here. Splitting the "90" and "150" bands should help to check the possible foreground residuals, but it's importance will strongly depend on the foreground model complexity.
The band splitting was studied with the performance-based Fisher forecasting code in this posting. Switching between the "optimal" and "force-split" path, you can see that it prefers split band when we reach the sensitivity \(\sigma(r) = 5\times 10^{-4}\).
In the current posting, we study the effect of the band splitting on the bias on the recovered r for 10 different sky models.

1: Summary Plot

In Figure 1, we summarize the r results, for the full frequency coverage and for the case where we drop the 85 and 145GHz channel. We show it for different values of \(A_L\), which is a crude approximation to delensing performances.
We also show the L=-2 log(Likelihood)/dof, where the Likelihood is the one returned by the minuit minimizer and dof is the number of degrees of freedom, i.e. \(N_{\rm bandpowers} * N_{\rm channels} * (N_{\rm channels}+1)/2 - N_{\rm param}\). Here \(N_{\rm bandpowers}=9\), \(N_{\rm params}=10\) and \(N_{\rm channels}\) is either 9 or 7. It is not exactly a reduced \(\chi^2\), since we are not using a Gaussian likelihood (we use Hammimeche-Lewis), but it is a good approximation.
For a description of the sky models used, see this posting.
For a table with the numbers plotted here, see this table.

Comments on the figure:

Figure 1: Summary of the results for different sky models (from 04.00 to 04.09, as described here).
In red, the analysis of the simulations with r=0.003, in green, the r=0 case. On the left, the case without decorrelation in the parametrization, on the right, with linear-\(\ell\) decorrelation. The outer error bars show the standard deviation \(\sigma\) of the \(N_{sims}\) simulations' ML results ( \(N_{sims}\simeq 500\) for models 00 to 06, 150 for 7, 8 and 9), and the inner error bars show \(\sigma/\sqrt{N_{sims}}\). We report both the r values, as well as L=-2log(Likelihood)/dof.


For the 10 foreground models considered in our data challenge, we don't see significant bias increase when removing the 85 and 145GHz channels. Note that except for model 5, most models don't have a strong decorrelation, as can be seen here. It is hard to assess how helpful the band splitting would be in the case of strong decorrelation, since we are not comparing cases with a similar effort.