Measuring time delays from strongly lensed supernovae (SNe) is emerging as a novel and independent tool for estimating the Hubble constant $(H_0)$. This is very important given the recent discord in the value of $H_0$ from two methods that probe different distance ranges. The success of this technique will rely of our ability to discover strongly lensed SNe with measurable time delays. Here, we present the magnifications and the time delays for the multiply-imaged galaxies behind the Hubble Frontier Fields (HFF) galaxy clusters, by using recently published lensing models. Continuing on our previous work done for Abell 1689 (A1689) and Abell 370, we also show the prospects of observing strongly lensed SNe behind the HFF clusters with the upcoming James Webb Space Telescope (JWST). With four 1-hour visits in one year, the summed expectations of all six HFF clusters are $\sim0.5$ core-collapse (CC) SNe and $\sim0.06$ Type Ia SNe (SNe Ia) in F115W band, while with F150W the expectations are higher, $\sim0.9$ CC SNe and $\sim0.06$ SNe Ia. These estimates match those expected by only surveying A1689, proving that the performance of A1689 as gravitational telescope is superior. In the five HFF clusters presented here, we find that F150W will be able to detect SNe Ia (SNe IIP) exploding in 93 (80) pairs multiply-imaged galaxies with time delays of less than 5 years.

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.

The stellar disk size of a galaxy depends on the ratio of the disk stellar mass to the halo mass, $m_\star \equiv M_\star/M_{\rm dh}$, and the fraction of the dark halo angular momentum transferred to the stellar disk, $j_\star \equiv J_\star/J_{\rm dh}$. Since $m_{\star}$ and $j_{\star}$ are determined by many star-formation related processes, measuring $j_\star$ and $m_\star$ at various redshifts is essential to understand the formation history of disk galaxies. We use the 3D-HST GOODS-S, COSMOS, and AEGIS imaging data and photo-$z$ catalog to examine $j_\star$ and $m_\star$ for star-forming galaxies at $z \sim$ 2, 3, and 4, when disks are actively forming. We find that the $j_\star/m_\star$ ratio is $\simeq 0.77\pm 0.06$ for all three redshifts over the entire mass range examined, $8\times 10^{10} < M_{\rm dh}/h^{-1} M_\odot < 2\times 10^{12}$, with a possible ($<30\%$) decrease with mass. This high ratio is close to those of local disk galaxies, descendants of our galaxies in terms of $M_{\rm dh}$ growth, implying a nearly constant $j_\star/m_\star$ over past 12 Gyr. These results are remarkable because mechanisms controlling angular momentum transfer to disks such as inflows and feedbacks depend on both cosmic time and halo mass and indeed theoretical studies tend to predict $j_\star/m_\star$ changing with redshift and mass. It is found that recent theoretical galaxy formation simulations predict smaller $j_{\star}/m_{\star}$ than our values. We also find that a significant fraction of our galaxies appears to be unstable against bar formation.

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

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

POLARBEAR-2 is a next-generation receiver for precision measurements of the polarization of the cosmic microwave background (Cosmic Microwave Background (CMB)). Scheduled to deploy in early 2015, it will observe alongside the existing POLARBEAR-1 receiver, on a new telescope in the Simons Array on Cerro Toco in the Atacama desert of Chile. For increased sensitivity, it will feature a larger area focal plane, with a total of 7,588 polarization sensitive antenna-coupled Transition Edge Sensor (TES) bolometers, with a design sensitivity of 4.1 uKrt(s). The focal plane will be cooled to 250 milliKelvin, and the bolometers will be read-out with 40x frequency domain multiplexing, with 36 optical bolometers on a single SQUID amplifier, along with 2 dark bolometers and 2 calibration resistors. To increase the multiplexing factor from 8x for POLARBEAR-1 to 40x for POLARBEAR-2 requires additional bandwidth for SQUID readout and well-defined frequency channel spacing. Extending to these higher frequencies requires new components and design for the LC filters which define channel spacing. The LC filters are cold resonant circuits with an inductor and capacitor in series with each bolometer, and stray inductance in the wiring and equivalent series resistance from the capacitors can affect bolometer operation. We present results from characterizing these new readout components. Integration of the readout system is being done first on a small scale, to ensure that the readout system does not affect bolometer sensitivity or stability, and to validate the overall system before expansion into the full receiver. We present the status of readout integration, and the initial results and status of components for the full array.

We propose a high thermal conductivity infrared (IR) filter using alumina for use in millimeter wave detection systems. We constructed a prototype two-layer anti-reflection (AR) coated alumina filter with a diameter of 100mm and a thickness of 2mm, and characterized its thermal and optical properties. The transmittance of this filter at 95GHz and 150GHz is 97% and 95% while the estimated 3dB cutoff frequency is at 450GHz. The high thermal conductivity of alumina minimizes thermal gradients. We measure a differential temperature of only 0.21K between the center and the edge of the filter when it is mounted on a thermal anchor of 77K. We also constructed a thermal model based on the prototype filter and analyzed the scalability of the filter diameter. We conclude that temperature increase at the center of alumina IR filter is less than 6K even with a large diameter of 500mm, when the temperature at the edge of the filter is 50K. This is suitable for an application to a large-throughput next-generation cosmic microwave background (CMB) polarization experiment, such as POLARBEAR-2 (PB-2).

We examined the anisotropic point spread function (PSF) of Suprime-Cam data utilizing dense star field data. We decomposed the PSF ellipticities into three components, the optical aberration, atmospheric turbulence, and chip-misalignment in an empirical manner, and evaluated the amplitude of each component. We found that, for long-exposure data, the optical aberration has the largest contribution to the PSF ellipticities, which could be modeled well by a simple analytic function based on the lowest-order aberration theory. The statistical properties of PSF ellipticities resulting from the atmospheric turbulence were investigated by using the numerical simulations. The simulation results are in a reasonable agreement with the observed data. It is also found that the optical PSF can be well corrected by the standard correction method with a polynomial fitting function. However, for the atmospheric PSF, its correction is affected by the common limitation caused by sparse sampling of PSFs due to a limited number of stars. We also examined the effects of the residual PSF anisotropies on Suprime-Cam cosmic shear data. We found that the shape and amplitude of the B-mode shear variance are broadly consistent with those of the residual PSF ellipticities measured from the dense star field data. This indicates that most of the sources of residual systematic are understood, which is an important step for cosmic shear statistics to be a practical tool of the precision cosmology.

We present an improved prediction of the nonlinear perturbation theory (PT) via the Lagrangian picture, which was originally proposed by Matsubara (2008). Based on the relations between the power spectrum in standard PT and that in Lagrangian PT, we derive analytic expressions for the power spectrum in Lagrangian PT up to 2-loop order in both real and redshift spaces. Comparing the improved prediction of Lagrangian PT with $N$-body simulations in real space, we find that the 2-loop corrections can extend the valid range of wave numbers where we can predict the power spectrum within 1% accuracy by a factor of 1.0 ($z=0.5$), 1.3 (1), 1.6 (2) and 1.8 (3) vied with 1-loop Lagrangian PT results. On the other hand, in all redshift ranges, the higher-order corrections are shown to be less significant on the two-point correlation functions around the baryon acoustic peak, because the 1-loop Lagrangian PT is already accurate enough to explain the nonlinearity on those scales in $N$-body simulations.

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 investigate the distance-redshift relation in a realistic inhomogeneous universe where the mass distribution is described by the mass function of Sheth and Tormen. It is found that the derived distance deviates systematically from the standard distance up to 10% depending on the choice of the lowest halo mass in which baryonic matter condensed to form luminous object such as galaxies. Remarkably the derived distance is well approximated by the Dyer-Roeder distance if we choose the clumpiness parameter \alpha calculated by our model.We also discuss the effect of inhomogeneities in the determination of dark energy parameter in the supernovae observation, and find that this effect must be taken into account for the future high redshift supernovae observation.

We study the general theorem about gravitational lensing which states the relationship between the numbers of images with different parities. Our formulation allows an extension to the nontransparent and singular model.

Non-transparent models of multipole expansion model and two point-mass model are analyzed from the catastrophe theory. Singularity behaviours of $2^n$-pole moments are discussed. We apply these models to triple quasar PG1115+080 and compare with the typical transparent model, softened power law spheroids. Multipole expansion model gives the best fit among them.

Singularities of caustics appeared in gravitational lensing effect are discussed analytically. Multipole expansion model of lensing object is mainly studied since it is tractable and universal. Our analyses are confirmed by numerical calculations and applied to multiple quasar system of PG1115+080. Consistencies with elliptical lens models are also discussed.