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.

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 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.

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.

LiteBIRD is a next-generation satellite mission to measure the polarization of the cosmic microwave background (CMB) radiation. On large angular scales the B-mode polarization of the CMB carries the imprint of primordial gravitational waves, and its precise measurement would provide a powerful probe of the epoch of inflation. The goal of LiteBIRD is to achieve a measurement of the characterizing tensor to scalar ratio $r$ to an uncertainty of $\delta r=0.001$. In order to achieve this goal we will employ a kilo-pixel superconducting detector array on a cryogenically cooled sub-Kelvin focal plane with an optical system at a temperature of 4~K. We are currently considering two detector array options; transition edge sensor (TES) bolometers and microwave kinetic inductance detectors (MKID). In this paper we give an overview of LiteBIRD and describe a TES-based polarimeter designed to achieve the target sensitivity of 2~$\mu$K$\cdot$arcmin over the frequency range 50 to 320~GHz.

We investigate constraints on the time variation of the fine structure constant between the recombination epoch and the present epoch, \Delta\alpha/\alpha \equiv (\alpha_{rec} - \alpha_{now})/\alpha_{now}, from cosmic microwave background (CMB) taking into account simultaneous variation of other physical constants, namely the electron mass m_{e} and the proton mass m_{p}. In other words, we consider the variation of Yukawa coupling and the QCD scale \Lambda_{QCD} in addition to the electromagnetic coupling. We clarify which parameters can be determined from CMB temperature anisotropy in terms of singular value decomposition. Assuming a relation among variations of coupling constants governed by a single scalar field (the dilaton), the 95% confidence level (C.L.) constraint on \Delta\alpha/\alpha is found to be -8.28 \times 10^{-3} < \Delta\alpha/\alpha < 1.81 \times 10^{-3}, which is tighter than the one obtained by considering only the change of \alpha and m_{e}. We also obtain the constraint on the time variation of the proton-to-electron mass ratio \mu \equiv m_{p}/m_{e} to be -0.52 < \Delta\mu/\mu < 0.17 (95% C.L.) under the same assumption. Finally, we also implement a forecast for constraints from the PLANCK survey.

As a test of the standard inflationary cosmology, which generically predicts nearly scale-invariant spectrum of primordial curvature fluctuations, we perform Markov-Chain Monte-Carlo analysis to search for possible modulations in the power spectrum and determine its shape together with the cosmological parameters using cosmic microwave background radiation data. By incorporating various three-parameter features on the simple power-law spectrum, we find an oscillatory modulation localized around the comoving wavenumber $k\simeq 0.009\mathrm{Mpc}^{-1}$ at 99.995% confidence level which improves the log-likelihood as much as $-\Delta 2\ln {\cal L}\equiv \Delta \chi^2_{\rm eff}=-22$. This feature can be detected even if we use only the cross correlation between the temperature and the E-mode polarization anisotropies.

The primordial curvature fluctuation spectrum is reconstructed by the maximum likelihood reconstruction method using the five-year Wilkinson Microwave Anisotropy Probe data of the cosmic microwave background temperature anisotropy. We apply the covariance matrix analysis and decompose the reconstructed spectrum into statistically independent band-powers. The prominent peak off a simple power-law spectrum found in our previous analysis turn out to be a $3.3\sigma$ deviation. From the statistics of primordial spectra reconstructed from mock observations, the probability that a primordial spectrum including such excess is realized in a power-law model is estimated to be about 2%.

The primordial curvature fluctuation spectrum is reconstructed by the cosmic inversion method using the five-year WMAP data of the cosmic microwave background temperature anisotropy. We apply the covariance matrix analysis and decompose the reconstructed spectrum into statistically independent band-powers. The statistically significant deviation from a simple power-law spectrum suggested by the analysis of the first-year data is not found in the five-year data except possibly at one point near the border of the wavenumber domain where accurate reconstruction is possible.

The constraints on the time variation of the fine structure constant at recombination epoch relative to its present value, $\Delta\alpha/\alpha \equiv (\alpha_{\mathrm{rec}} - \alpha_{\mathrm{now}})/\alpha_{\mathrm{now}}$, are obtained from the analysis of the 5-year WMAP cosmic microwave background data. As a result of Markov-Chain Monte-Carlo analysis, it is found that, contrary to the analysis based on the previous WMAP data, the mean value of $\Delta\alpha/\alpha=-0.0009$ does not change significantly whether we use the Hubble Space Telescope (HST) measurement of the Hubble parameter as a prior or not. The resultant 95% confidence ranges of $\Delta\alpha/\alpha$ are $-0.028 < \Delta\alpha/\alpha < 0.026$ with HST prior and $-0.050 < \Delta\alpha/\alpha < 0.042$ without HST prior.

We show a text-book potential for single-field inflation, namely, the Coleman-Weinberg model can induce double inflation and formation of primordial black holes (PBHs), because fluctuations that leave the horizon near the end of first inflation are anomalously enhanced at the onset of second inflation when the time-dependent mode turns to a growing mode rather than a decaying mode. The mass of PBHs produced in this mechanism lies in several discrete ranges depending on the model parameters. We also calculate the effects of non-Gaussian statistics due to higher-order interactions on the abundance of PBHs, which turns out to be small.

We present observational constraints on a scalar-tensor gravity theory by $\chi^2$ test for CMB anisotropy spectrum. We compare the WMAP temperature power spectrum with the harmonic attractor model, in which the scalar field has its harmonic effective potential with curvature $\beta$ in the Einstein conformal frame and the theory relaxes toward Einstein gravity with time. We found that the present value of the scalar coupling, i.e. the present level of deviation from Einstein gravity $(\alpha_0^2)$, is bounded to be smaller than $5\times 10^{-4-7\beta}$ ($2\sigma$), and $10^{-2-7\beta}$ ($4\sigma$) for $0< \beta<0.45$. This constraint is much stronger than the bound from the solar system experiments for large $\beta$ models, i.e., $\beta> 0.2$ and 0.3 in $2\sigma$ and $4\sigma$ limits, respectively. Furthermore, within the framework of this model, the variation of the gravitational constant at the recombination epoch is constrained as $|G(z=z_{rec})-G_0|/G_0 < 0.05(2\sigma)$, and $0.23(4\sigma)$.

Evolution of primordial fluctuations in a Brans-Dicke type scalar-tensor gravity theory is comprehensively investigated. The harmonic attractor model, in which the scalar field has its harmonic effective potential in the Einstein conformal frame and the theory relaxes toward Einstein gravity with time, is considered. The evolution of adiabatic initial perturbations in flat SCDM models is examined from the radiation-dominated epoch up to the present. We discuss how the scalar-tensor gravity affects the evolution of metric and matter perturbations, mainly focusing on the observational consequences, i.e., the matter power spectrum and the power spectrum of cosmic microwave background temperature. We find that the early time deviation is characterized only by the large static gravitational constant while the late time behavior is qualitatively different from that in Einstein gravity because the time variation of the gravitational constant and its fluctuation have non-negligible effects. The attracting scalar-tensor gravity affects only small scale modes due to its attracting nature, the degree of which is far beyond the post-Newtonian deviation at the present epoch.