We develop a framework for joint constraints on the CO luminosity function based on power spectra (PS) and voxel intensity distributions (VID), and apply this to simulations of COMAP, a CO intensity mapping experiment. This Bayesian framework is based on a Markov chain Monte Carlo (MCMC) sampler coupled to a Gaussian likelihood with a joint PS + VID covariance matrix computed from a large number of fiducial simulations, and re-calibrated with a small number of simulations per MCMC step. The simulations are based on dark matter halos from fast peak patch simulations combined with the $L_\text{CO}(M_\text{halo})$ model of Li et al. (2016). We find that the relative power to constrain the CO luminosity function depends on the luminosity range of interest. In particular, the VID is more sensitive at both small and large luminosities, while the PS is more sensitive at intermediate luminosities. The joint analysis is superior to using either observable separately. When averaging over CO luminosities ranging between $L_\text{CO} = 10^4-10^7L_\odot$, and over 10 cosmological realizations of COMAP Phase 2, the uncertainties (in dex) are larger by 58 % and 30 % for the PS and VID, respectively, when compared to the joint analysis (PS + VID). This method is generally applicable to any other random field, with a complicated likelihood, as long a fast simulation procedure is available.

Line-intensity mapping (LIM or IM) is an emerging field of observational work, with strong potential to fit into a larger effort to probe large-scale structure and small-scale astrophysical phenomena using multiple complementary tracers. Taking full advantage of such complementarity means, in part, undertaking line-intensity surveys with galaxy surveys in mind. We consider the potential for detection of a cross-correlation signal between COMAP and blind surveys based on photometric redshifts (as in COSMOS) or based on spectroscopic data (as with the HETDEX survey of Lyman-$\alpha$ emitters). We find that obtaining $\sigma_z/(1+z)\lesssim0.003$ accuracy in redshifts and $\gtrsim10^{-4}$ sources per Mpc$^3$ with spectroscopic redshift determination should enable a CO-galaxy cross spectrum detection significance at least twice that of the CO auto spectrum. Either a future targeted spectroscopic survey or a blind survey like HETDEX may be able to meet both of these requirements.

Warm dark matter cosmologies have been widely studied as an alternative to the cold dark matter paradigm, the characteristic feature being a suppression of structure formation on small cosmological scales. A very similar situation occurs if standard cold dark matter particles are kept in local thermal equilibrium with a, possibly dark, relativistic species until the universe has cooled down to keV temperatures. We perform a systematic phenomenological study of this possibility, and classify all minimal models containing dark matter and an arbitrary radiation component that allow such a late kinetic decoupling. We recover explicit cases recently discussed in the literature and identify new classes of examples that are very interesting from a model-building point of view. In some of these models dark matter is inevitably self-interacting, which is remarkable in view of recent observational support for this possibility. Hence, dark matter models featuring late kinetic decoupling have the potential not only to alleviate the missing satellites problem but also to address other problems of the cosmological concordance model on small scales, in particular the cusp-core and too-big-too-fail problems, in some cases without invoking any additional input.

We test for foreground residuals in the foreground cleaned Planck Cosmic Microwave Background (CMB) maps outside and inside U73 mask commonly used for cosmological analysis. The aim of this paper is to introduce a new method to validate masks by looking at the differences in cleaned maps obtained by different component separation methods. By analyzing the power spectrum as well as the mean, variance and skewness of needlet coefficients on bands outside and inside the U73 mask we first confirm that the pixels already masked by U73 are highly contaminated and cannot be used for cosmological analysis. We further find that the U73 mask needs extension in order to reduce large scale foreground residuals to a level of less than $20\%$ of the standard deviation of CMB fluctuations within the bands closest to the galactic equator. We also find 276 point sources in the cleaned foreground maps which are currently not masked by the U73 mask. Our final publicly available extended mask leaves $65.9\%$ of the sky for cosmological analysis. Note that this extended mask may be important for analyses on local sky patches; in full sky analyses the additional residuals near the galactic equator may average out.