We present a formulation of the observed number density fluctuations of gravitational-wave (GW) sources in a three dimensional space. In GW observations, redshift identification for each GW source is a challenging issue, in particular, for high redshift sources. The use of observed luminosity distance as a distance indicator will be a simple yet optimal way for measuring the clustering signal. We derive the density fluctuations of GW sources estimated from observed luminosity distance and sky position of each source. The density fluctuations are distorted as similar to the so-called redshift space distortions in galaxy surveys but with several differences. We then show the two-point correlation function and multipole power spectrum in the presence of the distortion effect. We find that the line-of-sight derivative of the lensing convergence, which does not appear in the redshift-space distortions, leads to a significant impact on the observed correlation function and higher-order multipole power spectrum.

We introduce the skew-spectrum statistic for weak lensing convergence $\kappa$ maps and test it against state-of-the-art high-resolution all-sky numerical simulations. We perform the analysis as a function of source redshift and smoothing angular scale for individual tomographic bins. We also analyse the cross-correlation between different tomographic bins. We compare the numerical results to fitting-functions used to model the bispectrum of the underlying density field as a function of redshift and scale. We derive a closed form expression for the skew-spectrum for gravity-induced secondary non-Gaussianity. We also compute the skew-spectrum for the projected $\kappa$ inferred from Cosmic Microwave Background (CMB) studies. As opposed to the low redshift case we find the post-Born corrections to be important in the modelling of the skew-spectrum for such studies. We show how the presence of a mask and noise can be incorporated in the estimation of a skew-spectrum.

We provide a new fitting formula of the matter bispectrum in the nonlinear regime calibrated by high-resolution cosmological $N$-body simulations of $41$ cold dark matter ($w$CDM, $w=$ constant) models around the Planck 2015 best-fit parameters. As the parameterization in our fitting function is similar to that in Halofit, our fitting is named BiHalofit. The simulation volume is sufficiently large ($> 10 \, {\rm Gpc}^3$) to cover almost all measurable triangle bispectrum configurations in the universe. The function is also calibrated using one-loop perturbation theory at large scales ($k<0.3 \, h \, {\rm Mpc}^{-1}$). Our formula reproduced the matter bispectrum to within $10 \, (15) \, \%$ accuracy in the Planck 2015 model at wavenumber $k< 3 \, (10) \, h \, {\rm Mpc}^{-1}$ and redshifts $z=0-3$. The other $40$ $w$CDM models obtained poorer fits, with accuracy approximating $20 \, \%$ at $k<3 \, h \, {\rm Mpc}^{-1}$ and $z=0-1.5$ (the deviation includes the $10 \, \%$-level sample variance of the simulations). We also provide a fitting formula that corrects the baryonic effects such as radiative cooling and active galactic nucleus feedback, using the latest hydrodynamical simulation IllustrisTNG. We demonstrate that our new formula more accurately predicts the weak-lensing bispectrum than the existing fitting formulas. This formula will assist current and future weak-lensing surveys and cosmic microwave background lensing experiments. Numerical codes of the formula are available, written in Python, C and Fortran.

We present new constraints on anisotropic birefringence of the cosmic microwave background polarization using two seasons of data from the Atacama Cosmology Telescope covering $456$ square degrees of sky. The birefringence power spectrum, measured using a curved-sky quadratic estimator, is consistent with zero. Our results provide the tightest current constraint on birefringence over a range of angular scales between $5$ arcminutes and $9$ degrees. We improve previous upper limits on the amplitude of a scale-invariant birefringence power spectrum by a factor of between $2$ and $3$. Assuming a nearly-massless axion field during inflation, our result is equivalent to a $2\,\sigma$ upper limit on the Chern-Simons coupling constant between axions and photons of $g_{\alpha\gamma}<4.0\times 10^{-2}/H_I$ where $H_I$ is the inflationary Hubble scale.

We construct cosmic microwave background lensing mass maps using data from the 2014 and 2015 seasons of observations with the Atacama Cosmology Telescope (ACT). These maps cover 2100 square degrees of sky and overlap with a wide variety of optical surveys. The maps are signal dominated on large scales and have fidelity such that their correlation with the cosmic infrared background is clearly visible by eye. We also create lensing maps with thermal Sunyaev-Zel'dovich contamination removed using a novel cleaning procedure that only slightly degrades the lensing signal-to-noise ratio. The cross-spectrum between the cleaned lensing map and the BOSS CMASS galaxy sample is detected at $10$-$\sigma$ significance, with an amplitude of $A=1.02 \pm 0.10$ relative to the Planck best-fit LCDM cosmological model with fiducial linear galaxy bias. Our measurement lays the foundation for lensing cross-correlation science with current ACT data and beyond.

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.

CMB-HD is a proposed ultra-deep (0.5 uk-arcmin), high-resolution (15 arcseconds) millimeter-wave survey over half the sky that would answer many outstanding questions in both fundamental physics of the Universe and astrophysics. This survey would be delivered in 7.5 years of observing 20,000 square degrees, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.

Families of cosmic inflation models predict a primordial gravitational-wave background that imprints B-mode polarization pattern in the Cosmic Microwave Background (CMB). High sensitivity instruments with wide frequency coverage and well-controlled systematic errors are needed to constrain the faint B-mode amplitude. We have developed antenna-coupled Transition Edge Sensor (TES) arrays for high-sensitivity polarized CMB observations over a wide range of millimeter-wave bands. BICEP Array, the latest phase of the BICEP/Keck experiment series, is a multi-receiver experiment designed to search for inflationary B-mode polarization to a precision $\sigma$(r) between 0.002 and 0.004 after 3 full years of observations, depending on foreground complexity and the degree of lensing removal. We describe the electromagnetic design and measured performance of BICEP Array low-frequency 40-GHz detector, their packaging in focal plane modules, and optical characterization including efficiency and beam matching between polarization pairs. We summarize the design and simulated optical performance, including an approach to improve the optical efficiency due to mismatch losses. We report the measured beam maps for a new broad-band corrugation design to minimize beam differential ellipticity between polarization pairs caused by interactions with the module housing frame, which helps minimize polarized beam mismatch that converts CMB temperature to polarization ($T \rightarrow P$) anisotropy in CMB maps.

Branches of cosmic inflationary models, such as slow-roll inflation, predict a background of primordial gravitational waves that imprints a unique odd-parity B-mode pattern in the Cosmic Microwave Background (CMB) at amplitudes that are within experimental reach. The BICEP/Keck (BK) experiment targets this primordial signature, the amplitude of which is parameterized by the tensor-to-scalar ratio r, by observing the polarized microwave sky through the exceptionally clean and stable atmosphere at the South Pole. B-mode measurements require an instrument with exquisite sensitivity, tight control of systematics, and wide frequency coverage to disentangle the primordial signal from the Galactic foregrounds. BICEP Array represents the most recent stage of the BK program, and comprises four BICEP3-class receivers observing at 30/40, 95, 150 and 220/270 GHz. The 30/40 GHz receiver will be deployed at the South Pole during the 2019/2020 austral summer. After 3 full years of observations with 30,000+ detectors, BICEP Array will measure primordial gravitational waves to a precision $\sigma(r)$ between 0.002 and 0.004, depending on foreground complexity and the degree of lensing removal. In this paper we give an overview of the instrument, highlighting the design features in terms of cryogenics, magnetic shielding, detectors and readout architecture as well as reporting on the integration and tests that are ongoing with the first receiver at 30/40 GHz.

The BICEP/Keck (BK) experiment aims to detect the imprint of primordial gravitational waves in the Cosmic Microwave Background polarization, which would be direct evidence of the inflation theory. While the tensor-to-scalar ratio has been constrained to be r_0.05 < 0.06 at 95% c.l., further improvements on this upper limit are hindered by polarized Galactic foreground emissions and removal of gravitational lensing polarization. The 30/40 GHz receiver of the BICEP Array (BA) will deploy at the end of 2019 and will constrain the synchrotron foreground with unprecedented accuracy within the BK sky patch. We will show the design of the 30/40 GHz detectors and test results summarizing its performance. The low optical and atmospheric loading at these frequencies requires our TES detectors to have low saturation power in order to be photon-noise dominated. To realize the low thermal conductivity required from a 250 mK base temperature, we developed new bolometer leg designs. We will present the relevant measured detector parameters: G, Tc, Rn, Psat , and spectral bands, and noise spectra. We achieved a per bolometer NEP including all noise components of 2.07E-17 W/sqrt(Hz), including an anticipated photon noise level 1.54E-17 W/sqrt(Hz).

The BICEP/Keck experiment (BK) is a series of small-aperture refracting telescopes observing degree-scale Cosmic Microwave Background (CMB) polarization from the South Pole in search of a primordial $B$-mode signature. This $B$-mode signal arises from primordial gravitational waves interacting with the CMB, and has amplitude parametrized by the tensor-to-scalar ratio $r$. Since 2016, BICEP3 and the Keck Array have been observing with 4800 total antenna-coupled transition-edge sensor detectors, with frequency bands spanning 95, 150, 220, and 270 GHz. Here we present the optical performance of these receivers from 2016 to 2019, including far-field beams measured in situ with an improved chopped thermal source and instrument spectral response measured with a field-deployable Fourier Transform Spectrometer. As a pair differencing experiment, an important systematic that must be controlled is the differential beam response between the co-located, orthogonally polarized detectors. We generate per-detector far-field beam maps and the corresponding differential beam mismatch that is used to estimate the temperature-to-polarization leakage in our CMB maps and to give feedback on detector and optics fabrication. The differential beam parameters presented here were estimated using improved low-level beam map analysis techniques, including efficient removal of non-Gaussian noise as well as improved spatial masking. These techniques help minimize systematic uncertainty in the beam analysis, with the goal of constraining the bias on $r$ induced by temperature-to-polarization leakage to be subdominant to the statistical uncertainty. This is essential as we progress to higher detector counts in the next generation of CMB experiments.

We describe an on-sky demonstration of a microwave-multiplexing readout system in one of the receivers of the Keck Array, a polarimetry experiment observing the cosmic microwave background at the South Pole. During the austral summer of 2018-2019, we replaced the time-division multiplexing readout system with microwave-multiplexing components including superconducting microwave resonators coupled to radio-frequency superconducting quantum interference devices at the sub-Kelvin focal plane, coaxial-cable plumbing and amplification between room temperature and the cold stages, and a SLAC Microresonator Radio Frequency system for the warm electronics. In the range 5-6 GHz, a single coaxial cable reads out 528 channels. The readout system is coupled to transition-edge sensors, which are in turn coupled to 150-GHz slot-dipole phased-array antennas. Observations began in April 2019, and we report here on an initial characterization of the system performance.

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 study the reconstruction of the cosmic rotation power spectrum produced by parity-violating physics, with an eye to ongoing and near future cosmic microwave background (CMB) experiments such as BICEP Array, CMBS4, LiteBIRD and Simons Observatory. In addition to the inflationary gravitational waves and gravitational lensing, measurements of other various effects on CMB polarization open new window into the early Universe. One of these is anisotropies of the cosmic polarization rotation which probes the Chern-Simons term generally predicted by string theory. The anisotropies of the cosmic rotation are also generated by the primordial magnetism and in the Standard Model extention framework. The cosmic rotation anisotropies can be reconstructed as quadratic in CMB anisotropies. However, the power of the reconstructed cosmic rotation is a CMB four-point correlation and is not directly related to the cosmic-rotation power spectrum. Understanding all contributions in the four-point correlation is required to extract the cosmic rotation signal. Assuming a scale-invariant rotation spectrum motivated by the inflationary cosmic-rotation models, we employ simulation to quantify each contribution to the four-point correlation and find that (1) a secondary contraction of the trispectrum increases the total signal-to-noise, (2) a bias from the lensing-induced trispectrum is significant compared to the statistical errors in, e.g., LiteBIRD and CMBS4-like experiments, (3) the use of a realization-dependent estimator decreases the statistical errors by 10%-20%, depending on experimental specifications, and (4) other higher-order contributions are negligible at least for near future 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.

Recent studies have demonstrated that {\em secondary} non-Gaussianity induced by gravity will be detected with a high signal-to-noise (S/N) by future and even by on-going weak lensing surveys. One way to characterise such non-Gaussianity is through the detection of a non-zero three-point correlation function of the lensing convergence field, or of its harmonic transform, the bispectrum. A recent study analysed the properties of the squeezed configuration of the bispectrum, when two wavenumbers are much larger than the third one. We extend this work by estimating the amplitude of the (reduced) bispectrum in four generic configurations, i.e., {\em squeezed, equilateral, isosceles} and {\em folded}, and for four different source redshifts $z_s=0.5,1.0,1.5,2.0$, by using an ensemble of all-sky high-resolution simulations. We compare these results against theoretical predictions. We find that, while the theoretical expectations based on widely used fitting functions can predict the general trends of the reduced bispectra, a more accurate theoretical modelling will be required to analyse the next generation of all-sky weak lensing surveys. The disagreement is particularly pronounced in the squeezed limit.

Precision measurements of cosmic microwave background (CMB) polarization require extreme control of instrumental systematics. In a companion paper we have presented cosmological constraints from observations with the BICEP2 and Keck Array experiments up to and including the 2015 observing season (BK15), resulting in the deepest CMB polarization maps to date and a statistical sensitivity to the tensor-to-scalar ratio of $\sigma(r) = 0.020$. In this work we characterize the beams and constrain potential systematic contamination from main beam shape mismatch at the three BK15 frequencies (95, 150, and 220 GHz). Far-field maps of 7,360 distinct beam patterns taken from 2010-2015 are used to measure differential beam parameters and predict the contribution of temperature-to-polarization leakage to the BK15 B-mode maps. In the multifrequency, multicomponent likelihood analysis that uses BK15, Planck, and WMAP maps to separate sky components, we find that adding this predicted leakage to simulations induces a bias of $\Delta r = 0.0027 \pm 0.0019$. Future results using higher-quality beam maps and improved techniques to detect such leakage in CMB data will substantially reduce this uncertainty, enabling the levels of systematics control needed for BICEP Array and other experiments that plan to definitively probe large-field inflation.

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 robust constraints on the stochastic gravitational waves (GWs) at Mpc scales from the cosmic microwave background (CMB) data. CMB constraints on GWs are usually characterized as the tensor-to-scalar ratio, assuming specifically a power-law form of the primordial spectrum, and are obtained from the angular spectra of CMB. Here, we relax the assumption of the power-law form, and consider to what extent one can constrain a monochromatic GW at shorter wavelengths. Previously, such a constraint has been derived at the wavelengths larger than the resolution scale of the CMB measurements, typically above $100$Mpc (below $10^{-16}$Hz in frequency). However, GWs whose wavelength is much shorter than $100$Mpc can imprint a small but non-negligible signal on CMB anisotropies at observed angular scales, $\ell<1000$. Here, using the CMB temperature, polarization, and lensing data set, we obtain the best constraints to date at $10^{-16}-10^{-14}$Hz of the GWs produced before the time of decoupling, which are tighter than those derived from the astrometric measurements and upper bounds on extra radiations. In the future, the constraints on GWs at Mpc scales will be further improved by several orders of magnitude with the precision $B$-mode measurement on large scales, $\ell<100$.

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

A millimeter-wave survey over half the sky, that spans frequencies in the range of 30 to 350 GHz, and that is both an order of magnitude deeper and of higher-resolution than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. By providing such a deep, high-resolution millimeter-wave survey (about 0.5 uK-arcmin noise and 15 arcsecond resolution at 150 GHz), CMB-HD will enable major advances. It will allow 1) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10/hMpc), which probes dark matter particle properties. It will also allow 2) measurements of the thermal and kinetic Sunyaev-Zel'dovich effects on small scales to map the gas density and gas pressure profiles of halos over a wide field, which probes galaxy evolution and cluster astrophysics. In addition, CMB-HD would allow us to cross critical thresholds in fundamental physics: 3) ruling out or detecting any new, light (< 0.1eV), thermal particles, which could potentially be the dark matter, and 4) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe. Such a survey would also 5) monitor the transient sky by mapping the full observing region every few days, which opens a new window on gamma-ray bursts, novae, fast radio bursts, and variable active galactic nuclei. Moreover, CMB-HD would 6) provide a census of planets, dwarf planets, and asteroids in the outer Solar System, and 7) enable the detection of exo-Oort clouds around other solar systems, shedding light on planet formation. CMB-HD will deliver this survey in 5 years of observing half the sky, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. The telescopes will field about 2.4 million detectors (600,000 pixels) in total.

We present the strongest constraints to date on anisotropies of cosmic microwave background (CMB) polarization rotation derived from 150 GHz data taken by the BICEP2/Keck Array CMB experiments up to and including the 2014 observing season (BK14). The definition of the polarization angle in BK14 maps has gone through self-calibration in which the overall angle is adjusted to minimize the observed TB and EB power spectra. After this procedure, the QU maps lose sensitivity to a uniform polarization rotation but are still sensitive to anisotropies of polarization rotation. This analysis places constraints on the anisotropies of polarization rotation, which could be generated by CMB photons interacting with axionlike pseudoscalar fields or Faraday rotation induced by primordial magnetic fields. The sensitivity of BK14 maps ($\sim 3\mu$K-arcmin) makes it possible to reconstruct anisotropies of the polarization rotation angle and measure their angular power spectrum much more precisely than previous attempts. Our data are found to be consistent with no polarization rotation anisotropies, improving the upper bound on the amplitude of the rotation angle spectrum by roughly an order of magnitude compared to the previous best constraints. Our results lead to an order of magnitude better constraint on the coupling constant of the Chern-Simons electromagnetic term $g_{a\gamma}\leq 7.2\times 10^{-2}/H_I$ (95% confidence) than the constraint derived from the B-mode spectrum, where $H_I$ is the inflationary Hubble scale. This constraint leads to a limit on the decay constant of $10^{-6}\lesssim f_a/M_{\rm pl}$ at mass range of $10^{-33}< m_a< 10^{-28}$ eV for $r=0.01$, assuming $g_{a\gamma}\sim\alpha/(2\pi f_a)$ with $\alpha$ denoting the fine structure constant. The upper bound on the amplitude of the primordial magnetic fields is 30nG (95% confidence) from the polarization rotation anisotropies.

Polarization data will soon provide the best avenue for measurements of the CMB lensing potential, although it is potentially sensitive to several instrumental effects including beam asymmetry, polarization angle uncertainties, sky coverage, as well as analysis choices such as masking. We derive "bias-hardened" lensing estimators to mitigate these effects, at the expense of somewhat larger reconstruction noise, and test them numerically on simulated data. We find that the mean-field bias from masking is significant for the EE quadratic lensing estimator, however the bias-hardened estimator combined with filtering techniques can mitigate the mean field. On the other hand, the EB estimator does not significantly suffer from the mean-field from the point source masking and survey window function. The contamination from beam asymmetry and polarization angle uncertainties, however, can generate mean-field biases for the EB estimator. These can also be mitigated using bias-hardened estimators, with at most a factor of ~ 3 degradation of noise level compared to the conventional approach.

We study the impact of the nonlinear growth of the large-scale structure (LSS) on the removal of the gravitational lensing effect (delensing) in cosmic microwave background (CMB) $B$ modes. The importance of the nonlinear growth of the LSS in the gravitational lensing analysis of CMB has been recently recognized by several works, while its impact on delensing is not yet explored. The delensing using mass-tracers such as galaxies and cosmic infrared background (CIB) could be also affected by the nonlinear growth. We find that the nonlinear growth of the LSS leads to $\sim 0.3\%$ corrections to $B$-mode spectrum after delensing with a high-$z$ mass tracer ($z_m\sim 2$) at $\ell=1000$-$2000$. The off-diagonal correlation coefficients of the lensing $B$-mode template spectrum become significant for delensing with low-$z$ tracers ($z_m\lesssim 0.5$), but are negligible with high-$z$ tracers (such as CIB). On the other hand, the power spectrum covariance of the delensed $B$ mode is not significantly affected by the nonlinear growth of the LSS, and the delensing efficiency is not significantly changed even if we use low-$z$ tracers. The CMB $B$-mode internal delensing is also not significantly affected by the nonlinear growth.

Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey would allow for major advances in measuring the distribution of dark matter and gas on small-scales, and yield needed insight on 1.) dark matter particle properties, 2.) the evolution of gas and galaxies, 3.) new light particle species, 4.) the epoch of inflation, and 5.) the census of bodies orbiting in the outer Solar System.

Cosmic microwave background (CMB) lensing is an integrated effect whose kernel is greater than half the peak value in the range $1<z<5$. Measuring this effect offers a powerful tool to probe the large-scale structure of the Universe at high redshifts. With the increasing precision of ongoing CMB surveys, other statistics than the lensing power spectrum, in particular the lensing bi-spectrum, will be measured at high statistical significance. This will provide ways to improve the constraints on cosmological models and lift degeneracies. Following on an earlier paper, we test analytical predictions of the CMB lensing bi-spectrum against full-sky lensing simulations, and discuss their validity and limitation in detail. The tree-level prediction of perturbation theory agrees with the simulation only up to $\ell\sim 200$, but the one-loop order allows capturing the simulation results up to $\ell\sim 600$. We also show that analytical predictions based on fitting formulas for the matter bi-spectrum agree reasonably well with simulation results, although the precision of the agreement depends on the configurations and scales considered. For instance, the agreement is at the $10\%$-level for the equilateral configuration at multipoles up to $\ell\sim2000$, but the difference in the squeezed limit raises to more than a factor of two at $\ell\sim2000$. This discrepancy appears to come from limitations in the fitting formula of the matter bi-spectrum. We also find that the analytical prediction for the post-Born correction to the bi-spectrum is in good agreement with the simulation. We conclude by discussing the bi-spectrum prediction in some theories of modified gravity.

The BICEP and Keck Array experiments are a suite of small-aperture refracting telescopes observing the microwave sky from the South Pole. They target the degree-scale B-mode polarization signal imprinted in the Cosmic Microwave Background (CMB) by primordial gravitational waves. Such a measurement would shed light on the physics of the very early universe. While BICEP2 observed for the first time a B-mode signal at 150 GHz, higher frequencies from the Planck satellite showed that it could be entirely due to the polarized emission from Galactic dust, though uncertainty remained high. Keck Array has been observing the same region of the sky for several years, with an increased detector count, producing the deepest polarized CMB maps to date. New detectors at 95 GHz were installed in 2014, and at 220 GHz in 2015. These observations enable a better constraint of galactic foreground emissions, as presented here. In 2015, BICEP2 was replaced by BICEP3, a 10 times higher throughput telescope observing at 95 GHz, while Keck Array is now focusing on higher frequencies. In the near future, BICEP Array will replace Keck Array, and will allow unprecedented sensitivity to the gravitational wave signal. High resolution observations from the South Pole Telescope (SPT) will also be used to remove the lensing contribution to B-modes.

We present results from an analysis of all data taken by the BICEP2/Keck CMB polarization experiments up to and including the 2015 observing season. This includes the first Keck Array observations at 220 GHz and additional observations at 95 & 150 GHz. The $Q/U$ maps reach depths of 5.2, 2.9 and 26 $\mu$K$_{cmb}$ arcmin at 95, 150 and 220 GHz respectively over an effective area of $\approx 400$ square degrees. The 220 GHz maps achieve a signal-to-noise on polarized dust emission approximately equal to that of Planck at 353 GHz. We take auto- and cross-spectra between these maps and publicly available WMAP and Planck maps at frequencies from 23 to 353 GHz. We evaluate the joint likelihood of the spectra versus a multicomponent model of lensed-$\Lambda$CDM+$r$+dust+synchrotron+noise. The foreground model has seven parameters, and we impose priors on some of these using external information from Planck and WMAP derived from larger regions of sky. The model is shown to be an adequate description of the data at the current noise levels. The likelihood analysis yields the constraint $r_{0.05}<0.07$ at 95% confidence, which tightens to $r_{0.05}<0.06$ in conjunction with Planck temperature measurements and other data. The lensing signal is detected at $8.8 \sigma$ significance. Running maximum likelihood search on simulations we obtain unbiased results and find that $\sigma(r)=0.020$. These are the strongest constraints to date on primordial gravitational waves.

BICEP Array is the newest multi-frequency instrument in the BICEP/Keck Array program. It is comprised of four 550 mm aperture refractive telescopes observing the polarization of the cosmic microwave background (CMB) at 30/40, 95, 150 and 220/270 GHz with over 30,000 detectors. We present an overview of the receiver, detailing the optics, thermal, mechanical, and magnetic shielding design. BICEP Array follows BICEP3's modular focal plane concept, and upgrades to 6" wafer to reduce fabrication with higher detector count per module. The first receiver at 30/40 GHz is expected to start observing at the South Pole during the 2019-20 season. By the end of the planned BICEP Array program, we project $\sigma(r) \sim 0.003$, assuming current modeling of polarized Galactic foreground and depending on the level of delensing that can be achieved with higher resolution maps from the South Pole Telescope.

BICEP3 is a 520mm aperture on-axis refracting telescope observing the polarization of the cosmic microwave background (CMB) at 95GHz in search of the B-mode signal originating from inflationary gravitational waves. BICEP3's focal plane is populated with modularized tiles of antenna-coupled transition edge sensor (TES) bolometers. BICEP3 was deployed to the South Pole during 2014-15 austral summer and has been operational since. During the 2016-17 austral summer, we implemented changes to optical elements that lead to better noise performance. We discuss this upgrade and show the performance of BICEP3 at its full mapping speed from the 2017 and 2018 observing seasons. BICEP3 achieves an order-of-magnitude improvement in mapping speed compared to a Keck 95GHz receiver. We demonstrate $6.6\mu K\sqrt{s}$ noise performance of the BICEP3 receiver.

Bicep Array is a cosmic microwave background (CMB) polarization experiment that will begin observing at the South Pole in early 2019. This experiment replaces the five Bicep2 style receivers that compose the Keck Array with four larger Bicep3 style receivers observing at six frequencies from 30 to 270GHz. The 95GHz and 150GHz receivers will continue to push the already deep Bicep/Keck CMB maps while the 30/40GHz and 220/270GHz receivers will constrain the synchrotron and galactic dust foregrounds respectively. Here we report on the design and performance of the Bicep Array instruments focusing on the mount and cryostat systems.

Targeting faint polarization patterns arising from Primordial Gravitational Waves in the Cosmic Microwave Background requires excellent observational sensitivity. Optical elements in small aperture experiments such as Bicep3 and Keck Array are designed to optimize throughput and minimize losses from transmission, reflection and scattering at millimeter wavelengths. As aperture size increases, cryostat vacuum windows must withstand larger forces from atmospheric pressure and the solution has often led to a thicker window at the expense of larger transmission loss. We have identified a new candidate material for the fabrication of vacuum windows: with a tensile strength two orders of magnitude larger than previously used materials, woven high-modulus polyethylene could allow for dramatically thinner windows, and therefore significantly reduced losses and higher sensitivity. In these proceedings we investigate the suitability of high-modulus polyethylene windows for ground-based CMB experiments, such as current and future receivers in the Bicep/Keck Array program. This includes characterizing their optical transmission as well as their mechanical behavior under atmospheric pressure. We find that such ultra-thin materials are promising candidates to improve the performance of large-aperture instruments at millimeter wavelengths, and outline a plan for further tests ahead of a possible upcoming field deployment of such a science-grade window.

BICEP Array is a degree-scale Cosmic Microwave Background (CMB) experiment that will search for primordial B-mode polarization while constraining Galactic foregrounds. BICEP Array will be comprised of four receivers to cover a broad frequency range with channels at 30/40, 95, 150 and 220/270 GHz. The first low-frequency receiver will map synchrotron emission at 30 and 40 GHz and will deploy to the South Pole at the end of 2019. In this paper, we give an overview of the BICEP Array science and instrument, with a focus on the detector module. We designed corrugations in the metal frame of the module to suppress unwanted interactions with the antenna-coupled detectors that would otherwise deform the beams of edge pixels. This design reduces the residual beam systematics and temperature-to-polarization leakage due to beam steering and shape mismatch between polarized beam pairs. We report on the simulated performance of single- and wide-band corrugations designed to minimize these effects. Our optimized design alleviates beam differential ellipticity caused by the metal frame to about 7% over 57% bandwidth (25 to 45 GHz), which is close to the level due the bare antenna itself without a metal frame. Initial laboratory measurements are also presented.

We investigate the impact of non-Gaussian lensing deflections on measurements of the CMB lensing power spectrum. We find that the false assumption of their Gaussianity significantly biases these measurements in current and future experiments at the percent level. The bias is detected by comparing CMB lensing reconstructions from simulated CMB data lensed with Gaussian deflection fields to reconstructions from simulations lensed with fully non-Gaussian deflection fields. The non-Gaussian deflections are produced by ray-tracing through snapshots of an N-body simulation and capture both the non-Gaussianity induced by non-linear structure formation and by multiple correlated deflections. We find that the amplitude of the measured bias is in agreement with analytical predictions by B\"ohm et al. 2016. The bias is largest in temperature-based measurements and we do not find evidence for it in measurements from a combination of polarization fields ($EB,EB$). We argue that the non-Gaussian bias should be even more important for measurements of cross-correlations of CMB lensing with low-redshift tracers of large-scale structure.

Cosmological structures grow differently in theories of gravity which are modified as compared to Einstein's General relativity (GR). Cosmic microwave background (CMB) fluctuation patterns at the last scattering surface are lensed by these structures along the photon path to the observer. The observed CMB pattern therefore keeps trace of the growth history of structures. We show that observations of the CMB lensing bi-spectrum offer an interesting way to constrain deviations from GR in a broad class of scalar-tensor theories of gravity called "beyond Horndeski". We quantify how the constraints on generic parameters describing the deviations from GR depend on the effective multipole range of the analysis. Our results further indicate that an accurate nonlinear correction of the matter bi-spectrum in the modified gravity considered is necessary when the bi-spectrum is used to probe scales beyond a multipole $\ell_{\rm max} \gtrsim 1500$. We also found that the results are insensitive to details of the implementation of the screening mechanism, at very small scales. We finally demonstrate the potential of the lensing bi-spectrum to provide a blind reconstruction of the redshift evolution of our modified gravity parameters by combining the analysis of CMB and low-z source lensing data.

Inflation is the leading theory of the first instant of the universe. Inflation, which postulates that the universe underwent a period of rapid expansion an instant after its birth, provides convincing explanation for cosmological observations. Recent advancements in detector technology have opened opportunities to explore primordial gravitational waves generated by the inflation through B-mode (divergent-free) polarization pattern embedded in the Cosmic Microwave Background anisotropies. If detected, these signals would provide strong evidence for inflation, point to the correct model for inflation, and open a window to physics at ultra-high energies. LiteBIRD is a satellite mission with a goal of detecting degree-and-larger-angular-scale B-mode polarization. LiteBIRD will observe at the second Lagrange point with a 400 mm diameter telescope and 2,622 detectors. It will survey the entire sky with 15 frequency bands from 40 to 400 GHz to measure and subtract foregrounds. The U.S. LiteBIRD team is proposing to deliver sub-Kelvin instruments that include detectors and readout electronics. A lenslet-coupled sinuous antenna array will cover low-frequency bands (40 GHz to 235 GHz) with four frequency arrangements of trichroic pixels. An orthomode-transducer-coupled corrugated horn array will cover high-frequency bands (280 GHz to 402 GHz) with three types of single frequency detectors. The detectors will be made with Transition Edge Sensor (TES) bolometers cooled to a 100 milli-Kelvin base temperature by an adiabatic demagnetization refrigerator.The TES bolometers will be read out using digital frequency multiplexing with Superconducting QUantum Interference Device (SQUID) amplifiers. Up to 78 bolometers will be multiplexed with a single SQUID amplidier. We report on the sub-Kelvin instrument design and ongoing developments for the LiteBIRD mission.

We present constraints on the patchy reionization by measuring the trispectrum of the Planck 2015 cosmic microwave background (CMB) temperature anisotropies. The patchy reionization leads to anisotropies in the CMB optical depth, and the statistics of the observed CMB anisotropies is altered. We estimate the trispectrum of the CMB temperature anisotropies to constrain spatial variation of the optical depth. We show that the measured trispectrum is consistent with that from the standard lensed CMB simulation at $2\sigma$. While no evidence of the patchy reionization is found in the Planck 2015 temperature trispectrum, the CMB constraint on the patchy reionization is significantly improved from previous works. Assuming the analytic bubble-halo model of Wang and Hu (2006), the constraint obtained in this work rules out the typical bubble size at the ionization fraction of $\sim0.5$ as $R\gtrsim 10$ Mpc. Further, our constraint implies that large-scale $B$-modes from the patchy reionization are not a significant contamination in detecting the primordial gravitational waves of $r\gtrsim0.001$ if the $B$ mode induced by the patchy reionization is described by Dvorkin et al. (2009). The CMB trispectrum data starts to provide meaningful constraints on the patchy reionization.

We present 108 full-sky gravitational lensing simulation data sets generated by performing multiple-lens plane ray-tracing through high-resolution cosmological $N$-body simulations. The data sets include full-sky convergence and shear maps from redshifts $z=0.05$ to $5.3$ at intervals of $150 \, h^{-1}{\rm Mpc}$ comoving radial distance (corresponding to a redshift interval of $\Delta z \simeq 0.05$ at the nearby universe), enabling the construction of a mock shear catalog for an arbitrary source distribution up to $z=5.3$. The dark matter halos are identified from the same $N$-body simulations with enough mass resolution to resolve the host halos of the Sloan Digital Sky Survey (SDSS) CMASS and Luminous Red Galaxies (LRGs). Angular positions and redshifts of the halos are provided by a ray-tracing calculation, enabling the creation of a mock halo catalog to be used for galaxy-galaxy and cluster-galaxy lensing. The simulation also yields maps of gravitational lensing deflections for a source redshift at the last scattering surface, and we provide 108 realizations of lensed cosmic microwave background (CMB) maps in which the post-Born corrections caused by multiple light scattering are included. We present basic statistics of the simulation data, including the angular power spectra of cosmic shear, CMB temperature and polarization anisotropies, galaxy-galaxy lensing signals for halos, and their covariances. The angular power spectra of the cosmic shear and CMB anisotropies agree with theoretical predictions within $5\%$ up to $\ell = 3000$ (or at an angular scale $\theta > 0.5$ arcsin). The simulation data sets are generated primarily for the ongoing Subaru Hyper Suprime-Cam survey but are freely available for download at http://cosmo.phys.hirosaki-u.ac.jp/takahasi/allsky_raytracing.

We explore the sensitivity of weak lensing observables to the expansion history of the universe and to the growth of cosmic structures, as well as the relative contribution of both effects to constraining cosmological parameters. We utilize ray-tracing dark-matter-only N-body simulations and validate our technique by comparing our results for the convergence power spectrum with analytic results from past studies. We then extend our analysis to non-Gaussian observables which cannot be easily treated analytically. We study the convergence (equilateral) bispectrum and two topological observables, lensing peaks and Minkowski functionals, focusing on their sensitivity to the matter density $\Omega_m$ and the dark energy equation of state $w$. We find that a cancelation between the geometry and growth effects is a common feature for all observables, and exists at the map level. It weakens the overall sensitivity by up to a factor of 3 and 1.5 for $w$ and $\Omega_m$, respectively, with the bispectrum worst affected. However, combining geometry and growth information alleviates the degeneracy between $\Omega_m$ and $w$ from either effect alone. As a result, the magnitude of marginalized errors remain similar to those obtained from growth-only effects, but with the correlation between the two parameters switching sign. These results shed light on the origin of cosmology-sensitivity of non-Gaussian statistics, and should be useful in optimizing combinations of observables.

We present a new method for delensing $B$ modes of the cosmic microwave background (CMB) using a lensing potential reconstructed from the same realization of the CMB polarization (CMB internal delensing). The $B$-mode delensing is required to improve sensitivity to primary $B$ modes generated by, e.g., the inflationary gravitational waves, axion-like particles, modified gravity, primordial magnetic fields, and topological defects such as cosmic strings. However, the CMB internal delensing suffers from substantial biases due to correlations between observed CMB maps to be delensed and that used for reconstructing a lensing potential. We construct a realization-dependent (RD) estimator for correcting these biases by deriving a general optimal estimator for higher-order correlations. The RD estimator is less sensitive to simulation uncertainties. Compared to the previous methods, we find that the RD estimator corrects the biases without substantial degradation of the delensing efficiency.

In the past several decades, the standard cosmological model has been established and its parameters have been measured to a high precision, while there are still many of the fundamental questions in cosmology; such as the physics in the very early Universe, the origin of the cosmic acceleration and the nature of the dark matter. The future world's largest radio telescope, Square Kilometre Array (SKA), will be able to open the new frontier of cosmology and will be one of the most powerful tools for cosmology in the next decade. The cosmological surveys conducted by the SKA would have the potential not only to answer these fundamental questions but also deliver the precision cosmology. In this article we briefly review the role of the SKA from the view point of the modern cosmology. The cosmology science led by the SKA-Japan Consortium (SKA-JP) Cosmology Science Working Group is also discussed.

BICEP3, the latest telescope in the BICEP/Keck program, started science observations in March 2016. It is a 550mm aperture refractive telescope observing the polarization of the cosmic microwave background at 95 GHz. We show the focal plane design and detector performance, including spectral response, optical efficiency and preliminary sensitivity of the upgraded BICEP3. We demonstrate 9.72$\mu$K$\sqrt{\textrm{s}}$ noise performance of the BICEP3 receiver.

BICEP3 is a small-aperture refracting cosmic microwave background (CMB) telescope designed to make sensitive polarization maps in pursuit of a potential B-mode signal from inflationary gravitational waves. It is the latest in the BICEP/Keck Array series of CMB experiments at the South Pole, which has provided the most stringent constraints on inflation to date. For the 2016 observing season, BICEP3 was outfitted with a full suite of 2400 optically coupled detectors operating at 95 GHz. In these proceedings we report on the far field beam performance using calibration data taken during the 2015-2016 summer deployment season in situ with a thermal chopped source. We generate high-fidelity per-detector beam maps, show the array-averaged beam profile, and characterize the differential beam response between co-located, orthogonally polarized detectors which contributes to the leading instrumental systematic in pair differencing experiments. We find that the levels of differential pointing, beamwidth, and ellipticity are similar to or lower than those measured for BICEP2 and Keck Array. The magnitude and distribution of BICEP3's differential beam mismatch - and the level to which temperature-to-polarization leakage may be marginalized over or subtracted in analysis - will inform the design of next-generation CMB experiments with many thousands of detectors.

BICEP3 is a 520 mm aperture, compact two-lens refractor designed to observe the polarization of the cosmic microwave background (CMB) at 95 GHz. Its focal plane consists of modularized tiles of antenna-coupled transition edge sensors (TESs), similar to those used in BICEP2 and the Keck Array. The increased per-receiver optical throughput compared to BICEP2/Keck Array, due to both its faster f/1.7 optics and the larger aperture, more than doubles the combined mapping speed of the BICEP/Keck program. The BICEP3 receiver was recently upgraded to a full complement of 20 tiles of detectors (2560 TESs) and is now beginning its second year of observation (and first science season) at the South Pole. We report on its current performance and observing plans. Given its high per-receiver throughput while maintaining the advantages of a compact design, BICEP3-class receivers are ideally suited as building blocks for a 3rd-generation CMB experiment, consisting of multiple receivers spanning 35 GHz to 270 GHz with total detector count in the tens of thousands. We present plans for such an array, the new "BICEP Array" that will replace the Keck Array at the South Pole, including design optimization, frequency coverage, and deployment/observing strategies.

The first discovery of the gravitational wave (GW) event, GW150914, suggests a higher merger rate of black-hole (BH) binaries. If this is true, a number of BH binaries will be observed via the second-generation GW detectors, and the statistical properties of the observed BH binaries can be scrutinized. A naive but important question to ask is whether the spatial distribution of BH binaries faithfully traces the matter inhomogeneities in the Universe or not. Although the BH binaries are thought to be formed inside the galaxies in most of the scenarios, there is no observational evidence to confirm such a hypothesis. Here, we estimate how well the second-generation GW detectors can statistically confirm the BH binaries to be a tracer of the large-scale structure by looking at the auto- and cross-correlation of BH binaries with photometric galaxies and weak lensing measurements, finding that, with a 3 year observation, the $>3\sigma$ detection of non-zero signal is possible if the BH merger rate today is $\dot{n}_0\gtrsim100$ Gpc$^{-3}$yr$^{-1}$ and the clustering bias of BH binaries is $b_{\rm BH,0}\gtrsim1.5$.

A linear polarization field on the sphere can be uniquely decomposed into an E-mode and a B-mode component. These two components are analytically defined in terms of spin-2 spherical harmonics. Maps that contain filtered modes on a partial sky can also be decomposed into E-mode and B-mode components. However, the lack of full sky information prevents orthogonally separating these components using spherical harmonics. In this paper, we present a technique for decomposing an incomplete map into E and B-mode components using E and B eigenmodes of the pixel covariance in the observed map. This method is found to orthogonally define E and B in the presence of both partial sky coverage and spatial filtering. This method has been applied to the BICEP2 and the Keck Array maps and results in reducing E to B leakage from LCDM E-modes to a level corresponding to a tensor-to-scalar ratio of $r<1\times10^{-4}$.

The B-mode polarization in the cosmic microwave background (CMB) anisotropies at large angular scales provides a smoking-gun evidence for the primordial gravitational waves (GWs). It is often stated that a discovery of the GWs establishes the quantum fluctuation of vacuum during the cosmic inflation. Since the GWs could also be generated by source fields, however, we need to check if a sizable signal exists due to such source fields before reaching a firm conclusion when the B-mode is discovered. Source fields of particular types can generate non-Gaussianity (NG) in the GWs. Testing statistics of the B-mode is a powerful way of detecting such NG. As a concrete example, we show a model in which a gauge field sources chiral GWs via a pseudoscalar coupling, and forecast the detection significance at the future CMB satellite LiteBIRD. Effects of residual foregrounds and lensing B-mode are both taken into account. We find the B-mode bispectrum "BBB" is in particular sensitive to the source-field NG, which is detectable at LiteBIRD with a $> 3 \sigma$ significance. Therefore the search for the "BBB" will be indispensable toward unambiguously establishing quantum fluctuation of vacuum when the B-mode is discovered. We also introduced the Minkowski functional to detect the NGs. While we find that the Minkowski functional is less efficient than the harmonic-space bispectrum estimator, it still serves as a useful cross check. Finally, we also discuss the possibility of extracting clean information on parity violation of GWs, and new types of parity-violating observables induced by lensing.

We present measurements of polarization lensing using the 150 GHz maps which include all data taken by the BICEP2 & Keck Array CMB polarization experiments up to and including the 2014 observing season (BK14). Despite their modest angular resolution ($\sim 0.5^\circ$), the excellent sensitivity ($\sim 3\mu$K-arcmin) of these maps makes it possible to directly reconstruct the lensing potential using only information at larger angular scales ($\ell\leq 700$). From the auto-spectrum of the reconstructed potential we measure an amplitude of the spectrum to be $A^{\phi\phi}_{\rm L}=1.15\pm 0.36$ (Planck $\Lambda$CDM prediction corresponds to $A^{\phi\phi}_{\rm L}=1$), and reject the no-lensing hypothesis at 5.8$\sigma$, which is the highest significance achieved to date using an EB lensing estimator. Taking the cross-spectrum of the reconstructed potential with the Planck 2015 lensing map yields $A^{\phi\phi}_{\rm L}=1.13\pm 0.20$. These direct measurements of $A^{\phi\phi}_{\rm L}$ are consistent with the $\Lambda$CDM cosmology, and with that derived from the previously reported BK14 B-mode auto-spectrum ($A^{\rm BB}_{\rm L}=1.20\pm 0.17$). We perform a series of null tests and consistency checks to show that these results are robust against systematics and are insensitive to analysis choices. These results unambiguously demonstrate that the B-modes previously reported by BICEP / Keck at intermediate angular scales ($150\lesssim\ell\lesssim 350$) are dominated by gravitational lensing. The good agreement between the lensing amplitudes obtained from the lensing reconstruction and B-mode spectrum starts to place constraints on any alternative cosmological sources of B-modes at these angular scales.

We discuss detectability of the nonlinear growth of the large-scale structure in the cosmic microwave background (CMB) lensing. Lensing signals involved in CMB anisotropies have been measured from multiple CMB experiments, such as Atacama Cosmology Telescope (ACT), Planck, POLARBEAR, and South Pole Telescope (SPT). Reconstructed lensing signals are useful to constrain cosmology via their angular power spectrum, while detectability and cosmological application of their bispectrum induced by the nonlinear evolution are not well studied. Extending the analytic estimate of the galaxy lensing bispectrum presented by Takada and Jain (2004) to the CMB case, we show that even near term CMB experiments such as Advanced ACT, Simons Array and SPT3G could detect the CMB lensing bispectrum induced by the nonlinear growth of the large-scale structure. In the case of the CMB Stage-IV, we find that the lensing bispectrum is detectable at $\gtrsim 50\sigma$ statistical significance. This precisely measured lensing bispectrum has rich cosmological information, and could be used to constrain cosmology, e.g., the sum of the neutrino masses and the dark-energy properties.

We investigate the statistical power of higher-order statistics and cross-correlation statistics to constrain the primordial non-Gaussianity from the imaging surveys. In particular, we consider the local-type primordial non- Gaussianity and discuss how well one can tightly constrain the higher-order non-Gaussian parameters ($g_{\rm NL}$ and $\tau_{\rm NL}$) as well as the leading order parameter $f_{\rm NL}$ from the halo/galaxy clustering and weak gravitational lensing measurements. Making use of a strong scale-dependent behavior in the galaxy/halo clustering, Fisher matrix analysis reveals that the bispectra can break the degeneracy between non-Gaussian parameters ($f_{\rm NL}$, $g_{\rm NL}$ and $\tau_{\rm NL}$) and this will give simultaneous constraints on those three parameters. The combination of cross-correlation statistics further improves the constraints by factor of 2. As a result, upcoming imaging surveys like the Large Synoptic Survey Telescope have the potential to improve the constraints on the primordial non-Gaussianity much tighter than those obtained from the CMB measurement by Planck, giving us an opportunity to test the single-sourced consistency relation, $\tau_{\rm NL} \ge (36/25) f_{\rm NL}^2$.

Gravitational waves (GWs) from compact binary stars at cosmological distances are promising and powerful cosmological probes, referred to as the GW standard sirens. With future GW detectors, we will be able to precisely measure source luminosity distances out to a redshift $z\sim5$. To extract cosmological information, previously proposed cosmological studies using the GW standard sirens rely on source redshift information obtained through an extensive electromagnetic follow-up campaign. However, the redshift identification is typically time consuming and rather challenging. Here, we propose a novel method for cosmology with the GW standard sirens free from the redshift measurements. Utilizing the anisotropies of the number density and luminosity distances of compact binaries originated from the large-scale structure, we show that, once GW observations will be well established in the future, (i) these anisotropies can be measured even at very high redshifts ($z\geq 2$), where the identification of the electromagnetic counterpart is difficult, (ii) the expected constraints on the primordial non-Gaussianity with the Einstein Telescope would be comparable to or even better than the other large-scale structure probes at the same epoch, and (iii) the cross-correlation with other cosmological observations is found to have high-statistical significance, providing additional cosmological information at very high redshifts.

We explore the potential use of the Radio Continuum (RC) survey conducted by the Square Kilometre Array (SKA) to remove (delens) the lensing-induced B-mode polarization and thus enhance future cosmic microwave background (CMB) searches for inflationary gravitational waves. Measurements of large-scale B-modes of the CMB are considered to be the best method for probing gravitational waves from the cosmic inflation. Future CMB experiments will, however, suffer from contamination by non-primordial B-modes, one source of which is the lensing B-modes. Delensing, therefore, will be required for further improvement of the detection sensitivity for gravitational waves. Analyzing the use of the two-dimensional map of galaxy distribution provided by the SKA RC survey as a lensing mass tracer, we find that joint delensing using near future CMB experiments and the SKA phase 1 will improve the constraints on the tensor-to-scalar ratio by more than a factor of $\sim 2$ compared to those without the delensing analysis. Compared to the use of CMB data alone, the inclusion of the SKA phase 1 data will increase the significance of the constraints on the tensor-to-scalar ratio by a factor $1.2$-$1.6$. For LiteBIRD combined with a ground-based experiment such as Simons Array and Advanced ACT, the constraint on the tensor-to-scalar ratio when adding SKA phase 2 data is improved by a factor of $2.3$-$2.7$, whereas delensing with CMB data alone improves the constraints by only a factor $1.3$-$1.7$. We conclude that the use of SKA data is a promising method for delensing upcoming CMB experiments such as LiteBIRD.

The B-mode polarization of the cosmic microwave background on large scales has been considered as a probe of gravitational waves from the cosmic inflation. Ongoing and future experiments will, however, suffer from contamination due to the B-modes of non-primordial origins, one of which is the lensing induced B-mode polarization. Subtraction of the lensing B-modes, usually referred to as delensing, will be required for further improvement of detection sensitivity of the gravitational waves. In such experiments, knowledge of statistical properties of the B-modes after delensing is indispensable to likelihood analysis particularly because the lensing B-modes are known to be non-Gaussian. In this paper, we study non-Gaussian structure of the delensed B-modes on large scales, comparing them with those of the lensing B-modes. In particular, we investigate the power spectrum correlation matrix and the probability distribution function (PDF) of the power spectrum amplitude. Assuming an experiment in which the quadratic delensing is an almost optimal method, we find that delensing reduces correlations of the lensing B-mode power spectra between different multipoles, and that the PDF of the power spectrum amplitude is well described as a normal distribution function with a variance larger than that in the case of a Gaussian field. These features are well captured by an analytic model based on the 4th order Edgeworth expansion. As a consequence of the non-Gaussianity, the constraint on the tensor-to-scalar ratio after delensing is degraded within approximately a few percent, which depends on the multipole range included in the analysis.

We present a method of lensing reconstruction on the full sky, by extending the optimal quadratic estimator proposed by Okamoto & Hu (2003) to the case including the curl mode of deflection angle. The curl mode is induced by the vector and tensor metric perturbations, and the reconstruction of the curl mode would be a powerful tool to not only check systematics in the estimated gradient mode but also probe any vector and tensor sources. We find that the gradient and curl modes can be reconstructed separately, thanks to the distinctive feature in the parity symmetry between the gradient and curl modes. We compare our estimator with the flat-sky estimator proposed by Cooray et al (2005). Based on the new formalism, the expected signal-to-noise ratio of the curl mode produced by the primordial gravitational-waves and a specific model of cosmic strings are estimated, and prospects for future observations are discussed.

We discuss the future detectability of gravitational-wave induced lensing from high-sensitivity cosmic microwave background (CMB) experiments. Gravitational waves can induce a rotational component of the weak-lensing deflection angle, usually referred to as the curl mode, which would be imprinted on the CMB maps. Using the technique of reconstructing lensing signals involved in CMB maps, this curl mode can be measured in an unbiased manner, offering an independent confirmation of the gravitational waves complementary to the B-mode polarization experiments. Based on the Fisher matrix analysis, we first show that with the noise levels necessary to confirm the consistency relation for the primordial gravitational waves, the future CMB experiments will be able to detect the gravitational-wave induced lensing signals. For a tensor-to-scalar ratio of $r < 0.1$, even if the consistency relation is difficult to confirm with a high significance, the gravitational-wave induced lensing would be detected at more than $3\,\sigma$ significance level. Further, we point out that high-sensitivity experiments will be also powerful to constrain the gravitational waves generated after the recombination epoch. Compared to the B-mode polarization, the curl mode is particularly sensitive to gravitational waves generated at low redshifts ($z < 10$) with a low frequency ($k < 10^{-3}$ Mpc$^{-1}$), and it could give a much tighter constraint on their energy density $\Omega_{\rm GW}$ by more than three orders of magnitude.

The lensing signals involved in CMB polarization maps have already been measured with ground-based experiments such as SPTpol and POLARBEAR, and would become important as a probe of cosmological and astrophysical issues in the near future. Sizes of polarization maps from ground-based experiments are, however, limited by contamination of long wavelength modes of observational noise. To further extract the lensing signals, we explore feasibility of measuring lensing signals from a collection of small sky maps each of which is observed separately by a ground-based large telescope, i.e., lensing reconstruction from a patchwork map of large sky coverage organized from small sky patches. We show that, although the B-mode power spectrum obtained from the patchwork map is biased due to baseline uncertainty, bias on the lensing potential would be negligible if the B-mode on scales larger than the blowup scale of $1/f$ noise is removed in the lensing reconstruction. As examples of cosmological applications, we also show 1) the cross-correlations between the reconstructed lensing potential and full-sky temperature/polarization maps from satellite missions such as PLANCK and LiteBIRD, and 2) the use of the reconstructed potential for delensing B-mode polarization of LiteBIRD observation.

The weak lensing effect on the cosmic microwave background (CMB) induces distortions in spatial pattern of CMB anisotropies, and statistical properties of CMB anisotropies become a weakly non-Gaussian field. We first summarize the weak lensing effect on the CMB (CMB lensing) in the presence of scalar, vector and tensor perturbations. Then we focus on the lensing effect on CMB statistics and methods to estimate deflection angles and their power spectrum. We end by summarizing recent observational progress and future prospect.

We present constraints on a cosmic string network with a measurement of weak gravitational lensing from CMB temperature map. The cosmic string network between observer and last scattering surface of CMB photons generates vector and/or tensor metric perturbations, and the deflection of CMB photons by these gravitational fields has curl mode which is not produced by the scalar metric perturbations. In this paper, we use the power spectrum of curl mode obtained from Planck to constrain the string tension, G\mu, and the reconnection probability, P. In demonstrating the parameter constraints with Planck curl mode, we also measure the lensing power spectrum from the Atacama Cosmology Telescope (ACT) 2008 season data, which have better angular resolution with lower instrumental noise on a much smaller chunk of the sky. Assuming P=1, the upper bound on tension is G\mu = 6.6\times 10^{-5} with 2\sigma (95% C.L.), using curl mode from Planck, which is weaker than that from the small-scale temperature power spectrum. For small values of P, however, the constraint from curl mode becomes tighter compared to that from temperature power spectrum. For P\lsim 10^{-2}, we obtain the constraint on the combination of the string parameters as G\mu P^{-1}\leq 3.4\times 10^{-5} at more than 2\sigma (95% C.L.).

We systematically derive full-sky formulae for the weak lensing power spectra generated by scalar, vector and tensor perturbations from the total angular momentum (TAM) method. Based on both the geodesic and geodesic deviation equations, we first give the gauge-invariant expressions for the deflection angle and Jacobi map as observables of the CMB lensing and cosmic shear experiments. We then apply the TAM method, originally developed in the theoretical studies of CMB, to a systematic derivation of the angular power spectra. The TAM representation, which characterizes the total angular dependence of the spatial modes projected along a line-of-sight, can carry all the information of the lensing modes generated by scalar, vector, and tensor metric perturbations. This greatly simplifies the calculation, and we present a complete set of the full-sky formulae for angular power spectra in both the E-/B-mode cosmic shear and gradient-/curl-mode lensing potential of deflection angle. Based on the formulae, we give illustrative examples of non-vanishing B-mode cosmic shear and curl-mode of deflection angle in the presence of the vector and tensor perturbations, and explicitly compute the power spectra.

We present new methods for lensing reconstruction from CMB temperature fluctuations which have smaller mean-field and reconstruction noise bias corrections than current lensing estimators, with minimal loss of signal-to-noise. These biases are usually corrected using Monte Carlo simulations, and to the extent that these simulations do not perfectly mimic the underlying sky there are uncertainties in the bias corrections. The bias-hardened estimators which we present can have reduced sensitivity to such uncertainties, and provide a desirable cross-check on standard results. To test our approach, we also show the results of lensing reconstruction from simulated temperature maps given on 100 deg^2, and confirm that our approach works well to reduce biases for a typical masked map in which 70 square masks each having 10 arcminute on a side exist, covering 2% of the simulated map, which is similar to the masks used in the current SPT lensing analysis.

We study the observational signature of vector metric perturbations through the effect of weak gravitational lensing. In the presence of vector perturbations, the non-vanishing signals for B-mode cosmic shear and curl-mode deflection angle, which have never appeared in the case of scalar metric perturbations, naturally arise. Solving the geodesic and geodesic deviation equations, we drive the full-sky formulas for angular power spectra of weak lensing signals, and give the explicit expressions for E-/B-mode cosmic shear and gradient-/curl-mode deflection angle. As a possible source for seeding vector perturbations, we then consider a cosmic string network, and discuss its detectability from upcoming weak lensing and CMB measurements. Based on the formulas and a simple model for cosmic string network, we calculate the angular power spectra and expected signal-to-noise ratios for the B-mode cosmic shear and curl-mode deflection angle. We find that the weak lensing signals are enhanced for a smaller intercommuting probability of the string network, $P$, and they are potentially detectable from the upcoming cosmic shear and CMB lensing observations. For $P\sim 10^{-1}$, the minimum detectable tension of the cosmic string will be down to $G\mu\sim 5\times 10^{-8}$. With a theoretically inferred smallest value $P\sim 10^{-3}$, we could even detect the string with $G\mu\sim 5\times 10^{-10}$.

We present forecast results for constraining the primordial non-Gaussianity from photometric surveys through a large-scale enhancement of the galaxy clustering amplitude. In photometric surveys, the distribution of observed galaxies at high redshifts suffers from the gravitational-lensing magnification, which systematically alters the number density for magnitude-limited galaxy samples. We estimate size of the systematic bias in the best-fit cosmological parameters caused by the magnification effect, particularly focusing on the primordial non-Gaussianity. For upcoming deep and/or wide photometric surveys like HSC, DES and LSST, the best-fit value of the non-Gaussian parameter, fNL, obtained from the galaxy count data is highly biased, and the true values of fNL would typically go outside the 3-sigma error of the biased confidence region, if we ignore the magnification effect in the theoretical template of angular power spectrum. The additional information from cosmic shear data helps not only to improve the constraint, but also to reduce the systematic bias. As a result, the size of systematic bias on fNL would become small enough compared to the expected 1-sigma error for HSC and DES, but it would be still serious for deep surveys with z_m > 1.5, like LSST. Tomographic technique improves the constraint on fNL by a factor of 2-3 compared to the one without tomography, but the systematic bias would increase.

We discuss the synergy of the cosmic shear and CMB lensing experiments to simultaneously constrain the neutrino mass and dark energy properties. Taking fully account of the CMB lensing, cosmic shear, CMB anisotropies, and their cross correlation signals, we clarify a role of each signal, and investigate the extent to which the upcoming observations by a high-angular resolution experiment of CMB and deep galaxy imaging survey can tightly constrain the neutrino mass and dark energy equation-of-state parameters. Including the primary CMB information as a prior cosmological information, the Fisher analysis reveals that the time varying equation-of-state parameters, given by w(a)=w_0+w_a(1-a), can be tightly constrained with the accuracies of 5% for w_0 and 15% for w_a, which are comparable to or even better than those of the stage-III type surveys neglecting the effect of massive neutrinos. In other words, including the neutrino mass in the parameter estimation would not drastically alter the Figure-of-Merit estimates of dark energy parameters from the weak lensing measurements. For the neutrino mass, a clear signal for total neutrino mass with 0.1 eV can be detected with 2-sigma significance. The robustness and sensitivity of these results are checked in detail by allowing the setup of cosmic shear experiment to vary as a function of observation time or exposure time, showing that the improvement of the constraints very weakly depends on the survey parameters, and the results mentioned above are nearly optimal for the dark energy parameters and the neutrino mass.