Foreground Separation and Constraints on Primordial Gravitational Waves with the PICO Space Mission by Ragnhild Aurlien et al. on Monday 28 NovemberPICO is a concept for a NASA probe-scale mission aiming to detect orconstrain the tensor to scalar ratio $r$, a parameter that quantifies theamplitude of inflationary gravity waves. We carry out map-based componentseparation on simulations with five foreground models and input $r$ values$r_{in}=0$ and $r_{in} = 0.003$. We forecast $r$ determinations using aGaussian likelihood assuming either no delensing or a residual lensing factor$A_{\rm lens}$ = 27%. By implementing the first full-sky, postcomponent-separation, map-domain delensing, we show that PICO should be able toachieve $A_{\rm lens}$ = 22% - 24%. For four of the five foreground models wefind that PICO would be able to set the constraints $r < 1.3 \times 10^{-4}\,\, \mbox{to} \,\, r <2.7 \times 10^{-4}\, (95\%)$ if $r_{in}=0$, thestrongest constraints of any foreseeable instrument. For these models,$r=0.003$ is recovered with confidence levels between $18\sigma$ and$27\sigma$. We find weaker, and in some cases significantly biased, upperlimits when removing few low or high frequency bands. The fifth model gives a$3\sigma$ detection when $r_{in}=0$ and a $3\sigma$ bias with $r_{in} = 0.003$.However, by correlating $r$ determinations from many small 2.5% sky areas withthe mission's 555 GHz data we identify and mitigate the bias. This analysisunderscores the importance of large sky coverage. We show that when only lowmultipoles $\ell \leq 12$ are used, the non-Gaussian shape of the truelikelihood gives uncertainties that are on average 30% larger than a Gaussianapproximation.arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.14342v1