Using only cosmic microwave background polarization data from the POLARBEAR experiment, we measure $B$-mode polarization delensing on subdegree scales at more than $5\sigma$ significance. We achieve a 14% $B$-mode power variance reduction, the highest to date for internal delensing, and improve this result to 2% by applying for the first time an iterative maximum a posteriori delensing method. Our analysis demonstrates the capability of internal delensing as a means of improving constraints on inflationary models, paving the way for the optimal analysis of next-generation primordial $B$-mode experiments.

Deep convolutional neural networks have been a popular tool for image generation and restoration. The performance of these networks is related to the capability of learning realistic features from a large dataset. In this work, we applied the problem of inpainting non-Gaussian signal, in the context of Galactic diffuse emissions at the millimetric and sub-millimetric regimes, specifically Synchrotron and Thermal Dust emission. Both of them are affected by contamination at small angular scales due to extra-galactic radio sources (the former) and to dusty star-forming galaxies (the latter). We consider the performances of a nearest-neighbors inpainting technique and compare it with two novels methodologies relying on generative Neural Networks. We show that the generative network is able to reproduce the statistical properties of the ground truth signal more consistently with high confidence level. The Python Inpainter for Cosmological and AStrophysical SOurces (PICASSO) is a package encoding a suite of inpainting methods described in this work and has been made publicly available.

We present a measurement of the gravitational lensing deflection power spectrum reconstructed with two seasons cosmic microwave background polarization data from the POLARBEAR experiment. Observations were taken at 150 GHz from 2012 to 2014 which survey three patches of sky totaling 30 square degrees. We test the consistency of the lensing spectrum with a Cold Dark Matter (CDM) cosmology and reject the no-lensing hypothesis at a confidence of 10.9 sigma including statistical and systematic uncertainties. We observe a value of A_L = 1.33 +/- 0.32 (statistical) +/- 0.02 (systematic) +/- 0.07 (foreground) using all polarization lensing estimators, which corresponds to a 24% accurate measurement of the lensing amplitude. Compared to the analysis of the first year data, we have improved the breadth of both the suite of null tests and the error terms included in the estimation of systematic contamination.

We describe the Fourier Transform Spectrometer (FTS) used for in-field testing of the POLARBEAR receiver, an experiment located in the Atacama Desert of Chile which measures the cosmic microwave background (CMB) polarization. The POLARBEAR-FTS (PB-FTS) is a Martin-Puplett interferometer designed to couple to the Huan Tran Telescope (HTT) on which the POLARBEAR receiver is installed. The PB-FTS measured the spectral response of the POLARBEAR receiver with signal-to-noise ratio (SNR) $>20$ for $\sim$69% of the focal plane detectors due to three features: a high throughput of 15.1 steradian cm$^{2}$, optimized optical coupling to the POLARBEAR optics using a custom designed output parabolic mirror, and a continuously modulated output polarizer. The PB-FTS parabolic mirror is designed to mimic the shape of the 2.5 m-diameter HTT primary reflector which allows for optimum optical coupling to the POLARBEAR receiver, reducing aberrations and systematics. One polarizing grid is placed at the output of the PB-FTS, and modulated via continuous rotation. This modulation allows for decomposition of the signal into different harmonics that can be used to probe potentially pernicious sources of systematic error in a polarization-sensitive instrument. The high throughput and continuous output polarizer modulation features are unique compared to other FTS calibrators used in the CMB field. In-field characterization of the POLARBEAR receiver was accomplished using the PB-FTS in April 2014. We discuss the design, construction, and operation of the PB-FTS and present the spectral characterization of the POLARBEAR receiver. We introduce future applications for the PB-FTS in the next-generation CMB experiment, the Simons Array.

Recent developments of transition-edge sensors (TESs), based on extensive experience in ground-based experiments, have been making the sensor techniques mature enough for their application on future satellite CMB polarization experiments. LiteBIRD is in the most advanced phase among such future satellites, targeting its launch in Japanese Fiscal Year 2027 (2027FY) with JAXA's H3 rocket. It will accommodate more than 4000 TESs in focal planes of reflective low-frequency and refractive medium-and-high-frequency telescopes in order to detect a signature imprinted on the cosmic microwave background (CMB) by the primordial gravitational waves predicted in cosmic inflation. The total wide frequency coverage between 34GHz and 448GHz enables us to extract such weak spiral polarization patterns through the precise subtraction of our Galaxy's foreground emission by using spectral differences among CMB and foreground signals. Telescopes are cooled down to 5Kelvin for suppressing thermal noise and contain polarization modulators with transmissive half-wave plates at individual apertures for separating sky polarization signals from artificial polarization and for mitigating from instrumental 1/f noise. Passive cooling by using V-grooves supports active cooling with mechanical coolers as well as adiabatic demagnetization refrigerators. Sky observations from the second Sun-Earth Lagrangian point, L2, are planned for three years. An international collaboration between Japan, USA, Canada, and Europe is sharing various roles. In May 2019, the Institute of Space and Astronautical Science (ISAS), JAXA selected LiteBIRD as the strategic large mission No. 2.

We report a 4.8$\sigma$ measurement of the cross-correlation signal between the cosmic microwave background (CMB) lensing convergence reconstructed from measurements of the CMB polarization made by the POLARBEAR experiment and the infrared-selected galaxies of the Herschel-ATLAS survey. This is the first measurement of its kind. We infer a best-fit galaxy bias of $b = 5.76 \pm 1.25$, corresponding to a host halo mass of $\log_{10}(M_h/M_\odot) =13.5^{+0.2}_{-0.3}$ at an effective redshift of $z \sim 2$ from the cross-correlation power spectrum. Residual uncertainties in the redshift distribution of the sub-mm galaxies are subdominant with respect to the statistical precision. We perform a suite of systematic tests, finding that instrumental and astrophysical contaminations are small compared to the statistical error. This cross-correlation measurement only relies on CMB polarization information that, differently from CMB temperature maps, is less contaminated by galactic and extra-galactic foregrounds, providing a clearer view of the projected matter distribution. This result demonstrates the feasibility and robustness of this approach for future high-sensitivity CMB polarization experiments.

We present the first measurement of cross-correlation between the lensing potential, reconstructed from cosmic microwave background (CMB) {\it polarization} data, and the cosmic shear field from galaxy shapes. This measurement is made using data from the POLARBEAR CMB experiment and the Subaru Hyper Suprime-Cam (HSC) survey. By analyzing an 11~deg$^2$ overlapping region, we reject the null hypothesis at 3.5$\sigma$\ and constrain the amplitude of the {\bf cross-spectrum} to $\widehat{A}_{\rm lens}=1.70\pm 0.48$, where $\widehat{A}_{\rm lens}$ is the amplitude normalized with respect to the Planck~2018{} prediction, based on the flat $\Lambda$ cold dark matter cosmology. The first measurement of this {\bf cross-spectrum} without relying on CMB temperature measurements is possible due to the deep POLARBEAR map with a noise level of ${\sim}$6\,$\mu$K-arcmin, as well as the deep HSC data with a high galaxy number density of $n_g=23\,{\rm arcmin^{-2}}$. We present a detailed study of the systematics budget to show that residual systematics in our results are negligibly small, which demonstrates the future potential of this cross-correlation technique.

We present a measurement of the $B$-mode polarization power spectrum of the cosmic microwave background (CMB) using taken from July 2014 to December 2016 with the POLARBEAR experiment. The CMB power spectra are measured using observations at 150 GHz with an instantaneous array sensitivity of $\mathrm{NET}_\mathrm{array}=23\, \mu \mathrm{K} \sqrt{\mathrm{s}}$ on a 670 square degree patch of sky centered at (RA, Dec)=($+0^\mathrm{h}12^\mathrm{m}0^\mathrm{s},-59^\circ18^\prime$). A continuously rotating half-wave plate is used to modulate polarization and to suppress low-frequency noise. We achieve $32\,\mu\mathrm{K}$-$\mathrm{arcmin}$ effective polarization map noise with a knee in sensitivity of $\ell = 90$, where the inflationary gravitational wave signal is expected to peak. The measured $B$-mode power spectrum is consistent with a $\Lambda$CDM lensing and single dust component foreground model over a range of multipoles $50 \leq \ell \leq 600$. The data disfavor zero $C_\ell^{BB}$ at $2.2\sigma$ using this $\ell$ range of POLARBEAR data alone. We cross-correlate our data with Planck high frequency maps and find the low-$\ell$ $B$-mode power in the combined dataset to be consistent with thermal dust emission. We place an upper limit on the tensor-to-scalar ratio $r < 0.90$ at 95% confidence level after marginalizing over foregrounds.

The Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

We present Atacama Large Millimeter/submillimiter Array (ALMA) high sensitivity ($\sigma_P \simeq 0.4\,$mJy) polarimetric observations at $97.5\,$GHz (Band 3) of a complete sample of $32$ extragalactic radio sources drawn from the faint Planck-ATCA Co-eval Observations (PACO) sample ($b<-75^\circ$, compact sources brighter than $200\,$mJy at $20\,$GHz). We achieved a detection rate of $~97\%$ at $3\,\sigma$ (only $1$ non-detection). We complement these observations with new Australia Telescope Compact Array (ATCA) data between $2.1$ and $35\,$GHz obtained within a few months and with data published in earlier papers from our collaboration. Adding the co-eval GaLactic and Extragalactic All-sky Murchison widefield array (GLEAM) survey detections between $70\,$ and $230\,$MHz for our sources, we present spectra over more than $3$ decades in frequency in total intensity and over about $1.7$ decades in polarization. The spectra of our sources are smooth over the whole frequency range, with no sign of dust emission from the host galaxy at mm wavelengths nor of a sharp high frequency decline due, for example, to electron ageing. We do however find indications of multiple emitting components and present a classification based on the number of detected components. We analyze the polarization fraction behaviour and distributions up to $97\,$GHz for different source classes. Source counts in polarization are presented at $95\,$GHz.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

While Type Ia Supernovae (SNe Ia) are one of the most mature cosmological probes, the next era promises to be extremely exciting in the number of different ways SNe Ia are used to measure various cosmological parameters. Here we review the experiments in the 2020s that will yield orders of magnitudes more SNe Ia, and the new understandings and capabilities to constrain systematic uncertainties at a level to match these statistics. We then discuss five different cosmological probes with SNe Ia: the conventional Hubble diagram for measuring dark energy properties, the distance ladder for measuring the Hubble constant, peculiar velocities and weak lensing for measuring sigma8 and strong-lens measurements of H0 and other cosmological parameters. For each of these probes, we discuss the experiments that will provide the best measurements and also the SN Ia-related systematics that affect each one.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

The Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $\mu$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $\mu$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5\sigma)$, or will rule out with more than $5\sigma$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4\sigma$, constrain the effective number of light particle species with $\Delta N_{\rm eff}<0.06~(2\sigma)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $\sigma_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science.

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.

The Simons Observatory (SO) is a set of cosmic microwave background instruments that will be deployed in the Atacama Desert in Chile. The key science goals include setting new constraints on cosmic inflation, measuring large scale structure with gravitational lensing, and constraining neutrino masses. Meeting these science goals with SO requires high sensitivity and improved calibration techniques. In this paper, we highlight a few of the most important instrument calibrations, including spectral response, gain stability, and polarization angle calibrations. We present their requirements for SO and experimental techniques that can be employed to reach those requirements.

The polarization of the atmosphere has been a long-standing concern for ground-based experiments targeting cosmic microwave background (CMB) polarization. Ice crystals in upper tropospheric clouds scatter thermal radiation from the ground and produce a horizontally-polarized signal. We report the detailed analysis of the cloud signal using a ground-based CMB experiment, POLARBEAR, located at the Atacama desert in Chile and observing at 150 GHz. We observe horizontally-polarized temporal increases of low-frequency fluctuations ("polarized bursts," hereafter) of $\lesssim$0.1 K when clouds appear in a webcam monitoring the telescope and the sky. The hypothesis of no correlation between polarized bursts and clouds is rejected with $>$24$\sigma$ statistical significance using three years of data. We consider many other possibilities including instrumental and environmental effects, and find no other reasons other than clouds that can explain the data better. We also discuss the impact of the cloud polarization on future ground-based CMB polarization experiments.

In this proceeding, we present studies of instrumental systematic effects for the Simons Obsevatory (SO) that are associated with the detector system and its interaction with the full SO experimental systems. SO will measure the Cosmic Microwave Background (CMB) temperature and polarization anisotropies over a wide range of angular scales in six bands with bandcenters spanning from 27 GHz to 270 GHz. We explore effects including intensity-to-polarization leakage due to coupling optics, bolometer nonlinearity, uncalibrated gain variations of bolometers, and readout crosstalk. We model the level of signal contamination, discuss proposed mitigation schemes, and present instrument requirements to inform the design of SO and future CMB projects.

The Simons Observatory (SO) will make precision temperature and polarization measurements of the cosmic microwave background (CMB) using a series of telescopes which will cover angular scales between one arcminute and tens of degrees and sample frequencies between 27 and 270 GHz. Here we present the current design of the large aperture telescope receiver (LATR), a 2.4 m diameter cryostat that will be mounted on the SO 6 m telescope and will be the largest CMB receiver to date. The cryostat size was chosen to take advantage of the large focal plane area having high Strehl ratios, which is inherent to the Cross-Dragone telescope design. The LATR will be able to accommodate thirteen optics tubes, each having a 36 cm diameter aperture and illuminating several thousand transition-edge sensor (TES) bolometers. This set of equipment will provide an opportunity to make measurements with unparalleled sensitivity. However, the size and complexity of the LATR also pose numerous technical challenges. In the following paper, we present the design of the LATR and include how we address these challenges. The solutions we develop in the process of designing the LATR will be informative for the general CMB community, and for future CMB experiments like CMB-S4.

The Simons Observatory (SO) will provide precision polarimetry of the cosmic microwave background (CMB) using a series of telescopes which will cover angular scales from arc-minutes to tens of degrees, contain over 60,000 detectors, and observe in frequency bands between 27 GHz and 270 GHz. SO will consist of a six-meter-aperture telescope initially coupled to ~35,000 detectors along with an array of 0.5m aperture refractive cameras, coupled to an additional 30,000+ detectors. The large aperture telescope receiver (LATR) is coupled to a six-meter crossed Dragone telescope and will be 2.4m in diameter, weigh over 3 tons, and have five cryogenic stages (80 K, 40 K, 4 K, 1 K and 100 mK). The LATR is coupled to the telescope via 13 independent optics tubes containing cryogenic optical elements and detectors. The cryostat will be cooled by by two Cryomech PT90 (80 K) and three Cryomech PT420 (40 K and 4 K) pulse tube cryocoolers, with cooling of the 1 K and 100 mK stages by a commercial dilution refrigerator system. The second component, the small aperture telescope (SAT), is a single optics tube refractive cameras of 42cm diameter. Cooling of the SAT stages will be provided by two Cryomech PT420, one of which is dedicated to the dilution refrigeration system which will cool the focal plane to 100 mK. SO will deploy a total of three SATs. In order to estimate the cool down time of the camera systems given their size and complexity, a finite difference code based on an implicit solver has been written to simulate the transient thermal behavior of both cryostats. The result from the simulations presented here predict a 35 day cool down for the LATR. The simulations suggest additional heat switches between stages would be effective in distribution cool down power and reducing the time it takes for the LATR to cool. The SAT is predicted to cool down in one week, which meets the SO design goals.

The Simons Observatory (SO) will observe the temperature and polarization anisotropies of the cosmic microwave background (CMB) over a wide range of frequencies (27 to 270 GHz) and angular scales by using both small (0.5 m) and large (6 m) aperture telescopes. The SO small aperture telescopes will target degree angular scales where the primordial B-mode polarization signal is expected to peak. The incoming polarization signal of the small aperture telescopes will be modulated by a cryogenic, continuously-rotating half-wave plate (CRHWP) to mitigate systematic effects arising from slowly varying noise and detector pair-differencing. In this paper, we present an assessment of some systematic effects arising from using a CRHWP in the SO small aperture systems. We focus on systematic effects associated with structural properties of the HWP and effects arising when operating a HWP, including the amplitude of the HWP synchronous signal (HWPSS), and I -> P (intensity to polarization) leakage that arises from detector non-linearity in the presence of a large HWPSS. We demonstrate our ability to simulate the impact of the aforementioned systematic effects in the time domain. This important step will inform mitigation strategies and design decisions to ensure that SO will meet its science goals.

The Simons Observatory (SO) will make precision temperature and polarization measurements of the cosmic microwave background (CMB) using a series of telescopes which will cover angular scales between one arcminute and tens of degrees, contain over 60,000 detectors, and sample frequencies between 27 and 270 GHz. SO will consist of a six-meter-aperture telescope coupled to over 30,000 detectors along with an array of half-meter aperture refractive cameras, which together couple to an additional 30,000+ detectors. SO will measure fundamental cosmological parameters of our universe, find high redshift clusters via the Sunyaev-Zeldovich effect, constrain properties of neutrinos, and seek signatures of dark matter through gravitational lensing. In this paper we will present results of the simulations of the SO large aperture telescope receiver (LATR). We will show details of simulations performed to ensure the structural integrity and thermal performance of our receiver, as well as will present the results of finite element analyses (FEA) of designs for the structural support system. Additionally, a full thermal model for the LATR will be described. The model will be used to ensure we meet our design requirements. Finally, we will present the results of FEA used to identify the primary vibrational modes, and planned methods for suppressing these modes. Design solutions to each of these problems that have been informed by simulation will be presented.

The Simons Observatory (SO) is a new experiment that aims to measure the cosmic microwave background (CMB) in temperature and polarization. SO will measure the polarized sky over a large range of microwave frequencies and angular scales using a combination of small ($\sim0.5 \, \rm m$) and large ($\sim 6\, \rm m $) aperture telescopes and will be located in the Atacama Desert in Chile. This work is part of a series of papers studying calibration, sensitivity, and systematic errors for SO. In this paper, we discuss current efforts to model optical systematic effects, how these have been used to guide the design of the SO instrument, and how these studies can be used to inform instrument design of future experiments like CMB-S4. While optical systematics studies are underway for both the small aperture and large aperture telescopes, we limit the focus of this paper to the more mature large aperture telescope design for which our studies include: pointing errors, optical distortions, beam ellipticity, cross-polar response, instrumental polarization rotation and various forms of sidelobe pickup.

New telescopes are being built to measure the Cosmic Microwave Background (CMB) with unprecedented sensitivity, including Simons Observatory (SO), CCAT-prime, the BICEP Array, SPT-3G, and CMB Stage-4. We present observing strategies for telescopes located in Chile that are informed by the tools used to develop recent Atacama Cosmology Telescope (ACT) and Polarbear surveys. As with ACT and Polarbear, these strategies are composed of scans that sweep in azimuth at constant elevation. We explore observing strategies for both small (0.42 m) aperture telescopes (SAT) and a large (6 m) aperture telescope (LAT). We study strategies focused on small sky areas to search for inflationary gravitational waves as well as strategies spanning roughly half the low-foreground sky to constrain the effective number of relativistic species and measure the sum of neutrino masses via the gravitational lensing signal due to large scale structure. We present these strategies specifically considering the telescope hardware and science goals of the SO, located at 23 degrees South latitude, 67.8 degrees West longitude. Observations close to the Sun and the Moon can introduce additional systematics by applying additional power to the instrument through telescope sidelobes. Significant side lobe contamination in the data can occur even at tens of degrees or more from bright sources. Therefore, we present several strategies that implement Sun and Moon avoidance constraints into the telescope scheduling. Strategies for resolving conflicts between simultaneously visible fields are discussed. We focus on maximizing telescope time spent on science observations. It will also be necessary to schedule calibration measurements, however that is beyond the scope of this work. The outputs of this study are algorithms that can generate specific schedule commands for the Simons Observatory instruments.

The Simons Observatory (SO) is an upcoming experiment that will study temperature and polarization fluctuations in the cosmic microwave background (CMB) from the Atacama Desert in Chile. SO will field both a large aperture telescope (LAT) and an array of small aperture telescopes (SATs) that will observe in six bands with center frequencies spanning from 27 to 270~GHz. Key considerations during the SO design phase are vast, including the number of cameras per telescope, focal plane magnification and pixel density, in-band optical power and camera throughput, detector parameter tolerances, and scan strategy optimization. To inform the SO design in a rapid, organized, and traceable manner, we have created a Python-based sensitivity calculator with several state-of-the-art features, including detector-to-detector optical white-noise correlations, a handling of simulated and measured bandpasses, and propagation of low-level parameter uncertainties to uncertainty in on-sky noise performance. We discuss the mathematics of the sensitivity calculation, the calculator's object-oriented structure and key features, how it has informed the design of SO, and how it can enhance instrument design in the broader CMB community, particularly for CMB-S4.

An estimation of the sky signal from streams of Time Ordered Data (TOD) acquired by Cosmic Microwave Background (\cmb) experiments is one of the most important steps in the context of \cmb data analysis referred to as the map-making problem. The continuously growing \cmb data sets render the \cmb map-making problem more challenging in terms of computational cost and memory in particular in the context of ground based experiments. In this context, we study a novel class of the Preconditioned Conjugate Gradient (PCG) solvers which invoke two-level preconditioners. We compare them against PCG solvers commonly used in the map-making context considering their precision and time-to-solution. We compare these new methods on realistic, simulated data sets reflecting the characteristics of current and forthcoming \cmb ground-based experiment. We develop an embarrassingly parallel implementation of the approach where each processor performs a sequential map-making for a subset of the TOD. We find that considering the map level residuals the new class of solvers permits achieving tolerance of up to 3 orders of magnitude better than the standard approach, where the residual level often saturates before convergence is reached. This corresponds to an important improvement in the precision of recovered power spectra in particular on the largest angular scales. The new method also typically requires fewer iterations to reach a required precision and thus shorter runtimes for a single map-making solution. However, the construction of an appropriate two-level preconditioner can be as costly as a single standard map-making run. Nevertheless, if the same problem needs to be solved multiple times, e.g., as in Monte Carlo simulations, this cost has to be incurred only once, and the method should be competitive not only as far as its precision but also its performance is concerned.

We combine the latest datasets obtained with different surveys to study the frequency dependence of polarized emission coming from Extragalactic Radio Sources (ERS). We consider data over a very wide frequency range starting from $1.4$ GHz up to $217$ GHz. This range is particularly interesting since it overlaps the frequencies of the current and forthcoming Cosmic Microwave Background (\cmb) experiments. Current data suggest that at high radio frequencies, ($ \nu \geq 20$ GHz) the fractional polarization of ERS does not depend on the total flux density. Conversely, recent datasets indicate a moderate increase of polarization fraction as a function of frequency, physically motivated by the fact that Faraday depolarization is expected to be less relevant at high radio-frequencies. We compute ERS number counts using updated models based on recent data, and we forecast the contribution of unresolved ERS in CMB polarization spectra. Given the expected sensitivities and the observational patch sizes of forthcoming \cmb\ experiments about $\sim 200 $ ( up to $\sim 2000 $ ) polarized ERS are expected to be detected. Finally, we assess that polarized ERS can contaminate the cosmological B-mode polarization if the tensor-to-scalar ratio is $r< 0.05$ and they have to be robustly controlled to de-lens \cmb\ B-modes at the arcminute angular scales.

We report an improved measurement of the cosmic microwave background (CMB) $B$-mode polarization power spectrum with the POLARBEAR experiment at 150 GHz. By adding new data collected during the second season of observations (2013-2014) to re-analyzed data from the first season (2012-2013), we have reduced twofold the band-power uncertainties. The band powers are reported over angular multipoles $500 \leq \ell \leq 2100$, where the dominant $B$-mode signal is expected to be due to the gravitational lensing of $E$-modes. We reject the null hypothesis of no $B$-mode polarization at a confidence of 3.1$\sigma$ including both statistical and systematic uncertainties. We test the consistency of the measured $B$-modes with the $\Lambda$ Cold Dark Matter ($\Lambda$CDM) framework by fitting for a single lensing amplitude parameter $A_L$ relative to the Planck best-fit model prediction. We obtain $A_L = 0.60 ^{+0.26} _{-0.24} ({\rm stat}) ^{+0.00} _{-0.04}({\rm inst}) \pm 0.14 ({\rm foreground}) \pm 0.04 ({\rm multi})$, where $A_{L}=1$ is the fiducial $\Lambda$CDM value, and the details of the reported uncertainties are explained later in the manuscript.

A continuously rotating half-wave plate (CRHWP) is a promising tool to improve the sensitivity to large angular scales in cosmic microwave background (CMB) polarization measurements. With a CRHWP, single detectors can measure three of the Stokes parameters, $I$, $Q$ and $U$, thereby avoiding the set of systematic errors that can be introduced by mismatches in the properties of orthogonal detector pairs. We focus on the implementation of CRHWPs in large aperture telescopes (i.e. the primary mirror is larger than the current maximum half-wave plate diameter of $\sim$0.5 m), where the CRHWP can be placed between the primary mirror and focal plane. In this configuration, one needs to address the intensity to polarization ($I{\rightarrow}P$) leakage of the optics, which becomes a source of 1/f noise and also causes differential gain systematics that arise from CMB temperature fluctuations. In this paper, we present the performance of a CRHWP installed in the POLARBEAR experiment, which employs a Gregorian telescope with a 2.5 m primary illumination pattern. The CRHWP is placed near the prime focus between the primary and secondary mirrors. We find that the $I{\rightarrow}P$ leakage is larger than the expectation from the physical properties of our primary mirror, resulting in a 1/f knee of 100 mHz. The excess leakage could be due to imperfections in the detector system, i.e. detector non-linearity in the responsivity and time-constant. We demonstrate, however, that by subtracting the leakage correlated with the intensity signal, the 1/f noise knee frequency is reduced to 32 mHz ($\ell \sim$39 for our scan strategy), which is very promising to probe the primordial B-mode signal. We also discuss methods for further noise subtraction in future projects where the precise temperature control of instrumental components and the leakage reduction will play a key role.

We present a model for simulating Carbon Monoxide (CO) rotational line emission in molecular clouds, taking account of their 3D spatial distribution in galaxies with different geometrical properties. The model implemented is based on recent results in the literature and has been designed for performing Monte-Carlo simulations of this emission. We compare the simulations produced with this model and calibrate them, both on the map level and on the power spectrum level, using the second release of data from the Planck satellite for the Galactic plane, where the signal-to-noise ratio is highest. We use the calibrated model to extrapolate the CO power spectrum at low Galactic latitudes where no high sensitivity observations are available yet. We then forecast the level of unresolved polarized emission from CO molecular clouds which could contaminate the power spectrum of Cosmic Microwave Background (CMB) polarization B-modes away from the Galactic plane. Assuming realistic levels of the polarization fraction, we show that the level of contamination is equivalent to a cosmological signal with $r \lesssim 0.02$. The Monte-Carlo MOlecular Line Emission (MCMole3D) Python package, which implements this model, is being made publicly available.

Analysis of cosmic microwave background (CMB) datasets typically requires some filtering of the raw time-ordered data. Filtering is frequently used to minimize the impact of low frequency noise, atmospheric contributions and/or scan synchronous signals on the resulting maps. In this work we explicitly construct a general filtering operator, which can unambiguously remove any set of unwanted modes in the data, and then amend the map-making procedure in order to incorporate and correct for it. We show that such an approach is mathematically equivalent to the solution of a problem in which the sky signal and unwanted modes are estimated simultaneously and the latter are marginalized over. We investigate the conditions under which this amended map-making procedure can render an unbiased estimate of the sky signal in realistic circumstances. We then study the effects of time-domain filtering on the noise correlation structure in the map domain, as well as impact it may have on the performance of the popular pseudo-spectrum estimators. We conclude that although maps produced by the proposed estimators arguably provide the most faithful representation of the sky possible given the data, they may not straightforwardly lead to the best constraints on the power spectra of the underlying sky signal and special care may need to be taken to ensure this is the case. By contrast, simplified map-makers which do not explicitly correct for time-domain filtering, but leave it to subsequent steps in the data analysis, may perform equally well and be easier and faster to implement. We focus on polarization-sensitive measurements targeting the B-mode component of the CMB signal and apply the proposed methods to realistic simulations based on characteristics of an actual CMB polarization experiment, POLARBEAR.

POLARBEAR-2 (PB-2) is a cosmic microwave background (CMB) polarization experiment that will be located in the Atacama highland in Chile at an altitude of 5200 m. Its science goals are to measure the CMB polarization signals originating from both primordial gravitational waves and weak lensing. PB-2 is designed to measure the tensor to scalar ratio, r, with precision {\sigma}(r) < 0.01, and the sum of neutrino masses, {\Sigma}m{\nu}, with {\sigma}({\Sigma}m{\nu}) < 90 meV. To achieve these goals, PB-2 will employ 7588 transition-edge sensor bolometers at 95 GHz and 150 GHz, which will be operated at the base temperature of 250 mK. Science observations will begin in 2017.

The E and B Experiment (EBEX) is a balloon-borne polarimeter designed to measure the polarization of the cosmic microwave background radiation and to characterize the polarization of galactic dust. EBEX was launched December 29, 2012 and circumnavigated Antarctica observing $\sim$6,000 square degrees of sky during 11 days at three frequency bands centered around 150, 250 and 410 GHz. EBEX was the first experiment to operate a kilo-pixel array of transition-edge sensor bolometers and a continuously rotating achromatic half-wave plate aboard a balloon platform. It also pioneered the use of detector readout based on digital frequency domain multiplexing. We describe the temperature calibration of the experiment. The gain response of the experiment is calibrated using a two-step iterative process. We use signals measured on passes across the Galactic plane to convert from readout-system counts to power. The effective smoothing scale of the EBEX optics and the star camera-to-detector offset angles are determined through \c{hi}2 minimization using the compact HII region RCW 38. This two-step process is initially performed with parameters measured before the EBEX 2013 flight and then repeated until the calibration factor and parameters converge.

We constrain anisotropic cosmic birefringence using four-point correlations of even-parity $E$-mode and odd-parity $B$-mode polarization in the cosmic microwave background measurements made by the POLARization of the Background Radiation (POLARBEAR) experiment in its first season of observations. We find that the anisotropic cosmic birefringence signal from any parity-violating processes is consistent with zero. The Faraday rotation from anisotropic cosmic birefringence can be compared with the equivalent quantity generated by primordial magnetic fields if they existed. The POLARBEAR nondetection translates into a 95% confidence level (C.L.) upper limit of 93 nanogauss (nG) on the amplitude of an equivalent primordial magnetic field inclusive of systematic uncertainties. This four-point correlation constraint on Faraday rotation is about 15 times tighter than the upper limit of 1380 nG inferred from constraining the contribution of Faraday rotation to two-point correlations of $B$-modes measured by Planck in 2015. Metric perturbations sourced by primordial magnetic fields would also contribute to the $B$-mode power spectrum. Using the POLARBEAR measurements of the $B$-mode power spectrum (two-point correlation), we set a 95% C.L. upper limit of 3.9 nG on primordial magnetic fields assuming a flat prior on the field amplitude. This limit is comparable to what was found in the Planck 2015 two-point correlation analysis with both temperature and polarization. We perform a set of systematic error tests and find no evidence for contamination. This work marks the first time that anisotropic cosmic birefringence or primordial magnetic fields have been constrained from the ground at subdegree scales.

Atmosphere is one of the most important noise sources for ground-based cosmic microwave background (CMB) experiments. By increasing optical loading on the detectors, it amplifies their effective noise, while its fluctuations introduce spatial and temporal correlations between detected signals. We present a physically motivated 3d-model of the atmosphere total intensity emission in the millimeter and sub-millimeter wavelengths. We derive a new analytical estimate for the correlation between detectors time-ordered data as a function of the instrument and survey design, as well as several atmospheric parameters such as wind, relative humidity, temperature and turbulence characteristics. Using an original numerical computation, we examine the effect of each physical parameter on the correlations in the time series of a given experiment. We then use a parametric-likelihood approach to validate the modeling and estimate atmosphere parameters from the POLARBEAR-I project first season data set. We derive a new 1.0% upper limit on the linear polarization fraction of atmospheric emission. We also compare our results to previous studies and weather station measurements. The proposed model can be used for realistic simulations of future ground-based CMB observations.