We report a measurement of the E-mode polarization power spectrum of the cosmic microwave background (CMB) using 150 GHz data taken from July 2014 to December 2016 with the POLARBEAR experiment. We reach an effective polarization map noise level of $32\,\mu\mathrm{K}$-$\mathrm{arcmin}$ across an observation area of 670 square degrees. We measure the EE power spectrum over the angular multipole range $500 \leq \ell <3000$, tracing the third to seventh acoustic peaks with high sensitivity. The statistical uncertainty on E-mode bandpowers is $\sim 2.3 \mu {\rm K}^2$ at $\ell \sim 1000$ with a systematic uncertainty of 0.5$\mu {\rm K}^2$. The data are consistent with the standard $\Lambda$CDM cosmological model with a probability-to-exceed of 0.38. We combine recent CMB E-mode measurements and make inferences about cosmological parameters in $\Lambda$CDM as well as in extensions to $\Lambda$CDM. Adding the ground-based CMB polarization measurements to the Planck dataset reduces the uncertainty on the Hubble constant by a factor of 1.2 to $H_0 = 67.20 \pm 0.57 {\rm km\,s^{-1} \,Mpc^{-1}}$. When allowing the number of relativistic species ($N_{eff}$) to vary, we find $N_{eff} = 2.94 \pm 0.16$, which is in good agreement with the standard value of 3.046. Instead allowing the primordial helium abundance ($Y_{He}$) to vary, the data favor $Y_{He} = 0.248 \pm 0.012$. This is very close to the expectation of 0.2467 from Big Bang Nucleosynthesis. When varying both $Y_{He}$ and $N_{eff}$, we find $N_{eff} = 2.70 \pm 0.26$ and $Y_{He} = 0.262 \pm 0.015$.

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.

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

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

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

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

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

In this paper, we show that single-field chaotic inflation on the brane with the potential $V=a \phi^p$ is compatible with the Swampland criteria. The spectral index and the running spectral index are within experimental bounds for $0<p \leq 2$. The tensor to scalar ratio is within observational bounds if $p \lesssim \mathcal{O}(1)$.

We examine the cosmology of warm dark matter (WDM), both stable and decaying, from the point of view of structure formation. We compare the matter power spectrum associated to WDM masses of 1.5 keV and 0.158 keV, with that expected for the stable cold dark matter $\Lambda$CDM$\equiv$SCDM paradigm, taken as our reference model. We scrutinize the effects associated to the warm nature of dark matter, as well as the fact that it decays. The decaying warm dark matter (DWDM) scenario is well-motivated, emerging in a broad class of particle physics theories where neutrino masses arise from the spontaneous breaking of a continuous global lepton number symmetry. The majoron arises as a Nambu-Goldstone boson, and picks up a mass from gravitational effects, that explicitly violate global symmetries. The majoron necessarily decays to neutrinos, with an amplitude proportional to their tiny mass, which typically gives it cosmologically long lifetimes. Using N-body simulations we show that our DWDM picture leads to a viable alternative to the $\Lambda$CDM scenario, with predictions that can differ substantially on small scales.

The recently announced results on the 21-cm absorption spectrum by the EDGES experiment can place very stringent limits on dark matter annihilation cross-sections. We properly take into account the heating energy released from dark matter annihilation from the radiation epoch to the 21-cm observation redshifts in the radiative transfer to compute the evolution of the gas temperature. Our results show that the global 21-cm absorption profile is a powerful cosmological probe of the dark matter interactions. For dark matter annihilating into electron-positron pairs, the EDGES results give a more stringent upper limit than the PLANCK result on the annihilation cross section at the lower dark matter mass region.

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

In this paper we study sneutrino dark matter in a recently proposed supersymmetric electroweak-scale inverse seesaw model, in which the majority of the sneutrino dark matter particle is a mixture of the right-handed sneutrino $\tilde{N}^c$ and the singlet field $\tilde{S}$. The scalar field $X$ responsible for the generation of neutrino masses can simultaneously play a crucial role for sneutrino annihilation in the early Universe via the pseudoscalar mediator $A_X$ into neutrinos. We focus here on the dominant annihilation channels and provide all the formulas together with analytic estimates in order to identify the relevant parameters. Furthermore, we show that the direct detection scattering cross section is many orders of magnitude below the current limits, and estimate the indirect detection annihilation rate, which is only a few orders of magnitude below the current limits.

With recent Lyman-alpha forest data from BOSS and XQ-100, some studies suggested that the lower mass limit on the fuzzy dark matter (FDM) particles is lifted up to $10^{-21}\,\mathrm{eV}$. However, such a limit was obtained by $\Lambda$CDM simulations with the FDM initial condition and the quantum pressure (QP) was not taken into account which could have generated non-trivial effects in large scales structures. We investigate the QP effects in cosmological simulations systematically, and find that the QP leads to further suppression of the matter power spectrum at small scales, as well as the halo mass function in the low mass end. We estimate the suppressing effect of QP in the 1D flux power spectrum of Lyman-alpha forest and compare it with data from BOSS and XQ-100. The rough uncertainties of thermal gas properties in the flux power spectrum model calculation were discussed. We conclude that more systematic studies, especially with QP taken into account, are necessary to constrain FDM particle mass using Lyman-alpha forest.

The Ultra-Light Axion (ULA) is a dark matter candidate with mass $\mathcal{O}(10^{-22})$eV and de Broglie wavelength of order kpc. Such an axion, also called the Fuzzy Dark Matter (FDM), thermalizes via the gravitational force and forms a Bose-Einstein condensate. Recent studies suggested that the quantum pressure from the FDM can significantly affect the structure formation in small scales, thus alleviating the so-called "small-scale crisis". In this paper, we develop a new technique to discretize the quantum pressure and illustrate the interactions among FDM particles in the $N$-body simulation, which accurately simulates the formation of the dark-matter halo and its inner structure in the region outside the softening length. In a self-gravitationally-bound virialized halo, we find a constant density, solitonic core, which is consistent with the theoretical prediction. The existence of the solitonic core reveals the non-linear effect of quantum pressure and impacts the structure formation in the FDM model.

We attempt to explain simultaneously the Galactic center gamma-ray excess and the 3.5 keV X-ray line from the Perseus cluster based on a class of non-abelian $SU(2)$ DM models, in which the dark matter and an excited state comprise a "dark" $SU(2)$ doublet. The non-abelian group kinetically mixes with the standard model gauge group via dimensions-5 operators. The dark matter particles annihilate into standard model fermions, followed by fragmentation and bremsstrahlung, and thus producing a continuous spectrum of gamma-rays. On the other hand, the dark matter particles can annihilate into a pair of excited states, each of which decays back into the dark matter particle and an X-ray photon, which has an energy equal to the mass difference between the dark matter and the excited state, which is set to be 3.5 keV. The large hierarchy between the required X-ray and $\gamma$-ray annihilation cross-sections can be achieved by a very small kinetic mixing between the SM and dark sector, which effectively suppresses the annihilation into the standard model fermions but not into the excited state.

In this paper, we consider large non-Gaussianity generated after D-term inflation in the case that the adiabatic curvature perturbation produced is dominated by a right-handed (RH) sneutrino curvaton. The cosmic string problem can also be evaded with around 10% contribution to the CMB power spectrum and the spectral index $n_s \simeq 1$. The non-Gaussianity produced can be as large as $10<f_{NL}<100$ for the right-handed sneutrino mass being around $m_\Phi=10 {GeV}$ with the Yukawa coupling $\lambda_\nu \sim 10^{-10}$. In this case, the lightest neutrino mass is $10^{-10}{eV}$ which favors RH sneutrino leptogenesis.

We present a global study of the simplest scalar phantom dark matter model. The best fit parameters of the model are determined by simultaneously imposing (i) relic density constraint from WMAP, (ii) 225 live days data from direct experiment XENON100, (iii) upper limit of gamma-ray flux from Fermi-LAT indirect detection based on dwarf spheroidal satellite galaxies, and (iv) the Higgs boson candidate with a mass about 125 GeV and its invisible branching ratio no larger than 40% if the decay of the Higgs boson into a pair of dark matter is kinematically allowed. The allowed parameter space is then used to predict annihilation cross sections for gamma-ray lines, event rates for three processes mono-b jet, single charged lepton and two charged leptons plus missing energies at the Large Hadron Collider, as well as to evaluate the muon anomalous magnetic dipole moment for the model.

An effective interaction approach is used to describe the interactions between the spin 0 or spin 1/2 dark matter particle and the degrees of freedom of the standard model. This approach is applicable to those models in which the dark matter particles do not experience the standard-model interactions, e.g., hidden-sector models. We explore the effects of these effective interaction operators on (i) dark matter relic density, (ii) spin-independent and spin-dependent dark matter-nucleon scattering cross sections, (iii) cosmic antiproton and gamma ray fluxes from the galactic halo due to dark matter annihilation, and (iv) monojet and monophoton production plus missing energy at the Tevatron and the Large Hadron Collider (LHC). We combine the experimental data of relic density from WMAP7, spin-independent cross section from XENON100, spin-dependent cross section from XENON10, ZEPLIN-III, and SIMPLE, cosmic antiproton flux from PAMELA, cosmic gamma-ray flux from ${\it Fermi}$-LAT, and the monojet and monophoton data from the Tevatron and the LHC, to put the most comprehensive limits on each effective operator.

Using an effective interaction approach to describe the interactions between the dark matter particle and the light degrees of freedom of the standard model, we calculate the gamma-ray flux due to the annihilation of the dark matter into quarks, followed by fragmentation into neutral pions which subsequently decay into photons. By comparison to the mid-latitude data released from the Fermi-LAT experiment, we obtain useful constraints on the size of the effective interactions and they are found to be comparable to those deduced from collider, gamma-ray line and anti-matter search experiments. However, the two operators induced by scalar and vector exchange among fermionic dark matter and light quarks that contribute to spin-independent cross sections are constrained more stringently by the recent XENON100 data.

In this paper we consider a class of inflation models on the brane where the dominant part of the inflaton scalar potential does not depend on the inflaton field value during inflation. In particular, we consider supernatural inflation, its hilltop version, A-term inflation, and supersymmetric (SUSY) D- and F-term hybrid inflation on the brane. We show that the parameter space can be broadened, the inflation scale generally can be lowered, and still possible to have the spectral index $n_s=0.96$.

Following \cite{Lin:2009yt}, we explore the parameter space of the case when the supersymmetry (SUSY) breaking scale is lower, for example, in gauge mediated SUSY breaking model. During inflation, the form of the potential is $V_0$ plus MSSM (MSSM stands for Minimal Supersymmetric Standard Model) (or A term (It is called A-term inflation when the inflaton field is any direction (gauge or singlet) that generates an A-term.)) inflation. We show that the model works for a wide range of the potential $V_0$ with the soft SUSY breaking mass $m\sim O(1)$ TeV. The implication to MSSM (or A-term) inflation is that the flat directions which is lifted by the non-renormalizable terms described by the superpotential $W=\lambda_p \phi^{p-1}/M^{p-3}_P$ with $p=4$ and $p=5$ are also suitable to be an inflaton field for $\lambda_p=O(1)$ provided there is an additional false vacuum term $V_0$ with appropriate magnitude. The flat directions corresponds to $p=6$ also works for $0 \lae V_0/M_P^4 \lae 10^{-40}$.

Motivated by a two-bump (or 1-peak plus 1-hump) structure in the ATIC data, we perform a statistical analysis fitting the PAMELA and ATIC data to a dark matter model, in which the dark matter particle can undergo both annihilation and decay. Using a chi-square analysis we show that both data can be simultaneously fitted better with such a double-action dark matter particle. We use an existing neutrino mass model in literature to illustrate the idea.

Supernatural inflation is an attractive model based just on a flat direction with soft SUSY breaking mass terms in the framework of supersymmetry. The beauty of the model is inferred from its name that the model needs no fine-tuning. However, the prediction of the spectral index is $n_s \gae 1$, in contrast to experimental data. In this paper, we show that the beauty of supernatural inflation with the spectral index reduced to $n_s=0.96$ without any fine-tuning, by considering the general feature that a flat direction is lifted by a non-renormalizable term with an A-term.

In this paper, we consider logarithmic radiative corrections and higher order terms to the supersymmetric hilltop F- and D-term hybrid inflation models. Conventional F- and D-term hybrid inflation only predicts $n_s \gae 0.98$. We show that via a positive quadratic and a negative quartic correction the spectral index can be reduced to $n_s=0.96$ suggested from latest WMAP result and also cosmic string problem appeared in SUSY hybrid inflation can be solved with mild tuning of the parameters if $\kappa \lae 0.01$ for F-term inflation and $g \lae 0.05$ for D-term inflation.

One of the most important and surprising discoveries in cosmology in recent years is the realization that our Universe is dominated by a mysterious dark energy, which leads to an accelerating expansion of space-time. A simple generalization of the standard Friedmann-Robertson-Walker equations based on General Relativity and the Cosmological Principles with the inclusion of a number of closed extra dimensions reproduces the currently observed data on dark energy, without the introduction of any cosmological constant or new particles. In particular, with a few extra dimensions, we obtain the redshift dependence of the deceleration parameter, the dark energy equation of state, as well as the age of the Universe in agreement with data, with essentially no free parameter. This model predicts that the extra dimensions have been compactified throughout the cosmic history, and it therefore suggests that signals from early Universe may give promising signatures of extra dimensions.

In a model of TeV right-handed (RH) neutrino by Krauss, Nasri, and Trodden, the sub-eV scale neutrino masses are generated via a 3-loop diagram with the vanishing see-saw mass forbidden by a discrete symmetry, and the TeV mass RH neutrino is simultaneously a novel candidate for the cold dark matter. However, we show that with a single RH neutrino it is not possible to generate two mass-square differences as required by the oscillation data. We extend the model by introducing one more TeV RH neutrino and show that it is possible to satisfy the oscillation pattern within the modified model. After studying in detail the constraints coming from the dark matter, lepton flavor violation and the muon anomalous magnetic moment, and the neutrinoless double beta decay, we explore the parameter space and derive predictions of the model. Finally, we study the production and decay signatures of the TeV RH neutrinos at TeV $e^+ e^-/\mu^+ \mu^-$ colliders.

We compare the tau neutrino flux arising from the galaxy and the earth atmosphere for $10^{3} \leq E/{GeV} \leq 10^{11}$. The intrinsic and oscillated tau neutrino fluxes from both sources are considered. We find that, for $E\geq 10^3$ GeV, the oscillated $\nu_{\tau}$ flux along the galactic plane dominates over the maximal intrinsic atmospheric $\nu_{\tau}$ flux, i.e., the flux along the horizontal direction. We also briefly comment on the prospects for observing these high-energy tau neutrinos.

We compare the tau neutrino flux arising from the galaxy and the earth atmosphere for 10^3 < E/GeV < 10^11. The intrinsic and oscillated tau neutrino fluxes from both sources are calculated. The intrinsic galactic \nu_{\tau} flux (E > 10^3 GeV) is calculated by considering the interactions of high-energy cosmic-rays with the matter present in our galaxy, whereas the oscillated galactic \nu_{\tau} flux is coming from the oscillation of the galactic \nu_{\mu} flux. For the intrinsic atmospheric \nu_{\tau} flux, we extend the validity of a previous calculation from E < 10^6 GeV up to E < 10^11 GeV. The oscillated atmospheric \nu_{\tau} flux is, on the other hand, rather suppressed. We find that, for 10^3 < E/GeV < 5\cdot 10^7, the oscillated \nu_{\tau} flux along the galactic plane dominates over the maximal intrinsic atmospheric \nu_{\tau} flux, i.e., the flux along the horizontal direction. We also briefly mention the presently envisaged prospects for observing these high-energy tau neutrinos.