We present spectroscopic confirmation of three new two-image gravitationally lensed quasars, compiled from existing strong lens and X-ray catalogs. Images of HSC J091843.27-022007.5 show a red galaxy with two blue point sources at either side, separated by 2.26 arcsec. This system has a source and a lens redshifts $z_s = 0.804$ and $z_l = 0.459$, respectively, as obtained by our follow-up spectroscopic data. CXCO J100201.50+020330.0 shows two point sources separated by 0.85 arcsec on either side of an early-type galaxy. The follow-up spectroscopic data confirm the fainter quasar has the same redshift with the brighter quasar from the SDSS fiber spectrum at $z_s = 2.016$. The deflecting foreground galaxy is a typical early-type galaxy at a redshift of $z_l = 0.439$. SDSS J135944.21+012809.8 has two point sources with quasar spectra at the same redshift $z_s = 1.096$, separated by 1.05 arcsec, and fits to the HSC images confirm the presence of a galaxy between these. These discoveries demonstrate the power of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP)'s deep imaging and wide sky coverage. Combined with existing X-ray source catalogues and follow-up spectroscopy, the HSC-SSP provides us unique opportunities to find multiple-image quasars lensed by a foreground galaxy.

We report the largest sample of candidate strong gravitational lenses belonging to the Survey of Gravitationally-lensed Objects in HSC Imaging for group-to-cluster scale (SuGOHI-c) systems. These candidates are compiled from the S18A data release of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) Survey. We visually inspect $\sim39,500$ galaxy clusters, selected from several catalogs, overlapping with the Wide, Deep, and UltraDeep fields, spanning the cluster redshift range $0.05<z_{cl}<1.38$. We discover 641 candidate lens systems, of which 536 are new. From the full sample, 47 are almost certainly bonafide lenses, 181 of them are highly probable lenses and 413 are possible lens systems. Additionally, we present 131 lens candidates at galaxy-scale serendipitously discovered during the inspection. We obtained spectroscopic follow-up of 10 candidates using the X-shooter. With this follow-up, we confirm 8 systems as strong gravitational lenses. Of the remaining two, one of the sources is too faint to detect any emission, and the other has a tentative redshift close to the lens redshift, but additional arcs in this system are yet to be observed spectroscopically. Since the HSC-SSP is an ongoing survey, we expect to find $\sim600$ definite or probable lenses using this procedure and even more if combined with other lens finding methods.

We report our high-resolution (0.03 arcsec to 0.07 arcsec) ALMA imaging of the quadruply lensed radio-loud quasar MG\,J0414+0534 at redshift $z=2.639$ in the continuum and the broad CO(11-10) line at $\sim 340\,$GHz. With the help of strong lensing magnification and ALMA's high-resolution, we succeeded in resolving the jet/dust and CO gas in the quasar host galaxy both extending up to $\sim 1\,$kpc, with a resolution of $\sim 50\,$pc for the first time. Both the continuum emission and the CO(11-10) line have a similar bimodal structure aligned with the quasar jets ($\sim 200\,$pc) observed by Very Long Baseline Interferometry (VLBI) at $5\,$GHz and $8.4\,$GHz. The CO gas in the vicinity of both the eastern and western jet components at the location of $\sim 80\,$pc from the quasar core are moving at high velocities up to $\pm 600\,$$\textrm{km}\,\textrm{s}^{-1}$ relative to the core. The observed features show clear evidence of strong interaction between the jets and interstellar medium (ISM). High temperature and high density environments in the ISM of the quasar host galaxy, as suggested from CO line ratios, also support this result. The small scale of the jets, the jet-ISM interaction, and the continuum spectral energy distribution of this source indicate that we are watching the infancy stage of quasar radio activity.

We report the serendipitous discovery of HSC J0904$-$0102, a quadruply-lensed Lyman break galaxy (LBG) in the Survey of Gravitationally-lensed Objects in Hyper Suprime-Cam Imaging (SuGOHI). Owing to its point-like appearance, the source was thought to be a lensed active galactic nucleus. We obtained follow-up spectroscopic data with the Gemini Multi-Object Spectrographs on the Gemini South Telescope, which confirmed this to be a lens system. The deflecting foreground galaxy is a typical early-type galaxy at a high redshift of $z_{\ell} = 0.957$ with stellar velocity dispersion $\sigma_v=259\pm56$ km~s$^{-1}$. The lensed source is identified as an LBG at $z_{\rm s} = 3.403$, based on the sharp drop bluewards of Ly$\alpha$ and other absorption features. A simple lens mass model for the system, assuming a singular isothermal ellipsoid, yields an Einstein radius of $\theta_{\rm Ein} = 1. 23^{\prime\prime}$ and a total mass within the Einstein radius of $M_{\rm Ein} = (5.55\pm 0.24) \times 10^{11}M_{\odot}$ corresponding to a velocity dispersion of $\sigma_{\rm SIE}= 283\pm 3$ km~s$^{-1}$, which is in good agreement with the value derived spectroscopically. The most isolated lensed LBG image has a magnification of $\sim 6.5$. In comparison with other lensed LBGs and typical $z\sim4$ LBG populations, HSC J0904$-$0102 is unusually compact, an outlier at $>2\sigma$ confidence. Together with a previously discovered SuGOHI lens, HSC J1152$+$0047, that is similarly compact, we believe that the HSC Survey is extending LBG studies down to smaller galaxy sizes.

The cosmic microwave background (CMB) contains an anomalous cold spot with a surrounding hot ring, known as the Cold Spot. Inoue & Silk (2006) proposed that this feature could be explained by postulating a supervoid: if such a large underdense region exists, then the growth of matter perturbing around the spot might differ from the average value in the Universe and the differences might affect weak lensing analysis of peak statistics. To investigate environmental effects on halo number count and peak statistics, we used a publicly available ray-tracing simulation for a box size of 2250$h^{-1}$Mpc on a side (Takahashi et al. 2017). We found that the number counts for massive haloes toward the largest underdense region in the simulation decreases and the corresponding significance of the difference, based on a cosmic average, is $\geq3\sigma$. On the basis of the results of peak statistics analysis, the number of high peaks decreases with the decrement of massive haloes, but the number of low peaks increases with the lack of matter in the line of sight. The highest significance of the decrement in peak counts in large underdense regions is $5\sigma$ in the total signal-to-noise ratio. Our result implies that environmental effects on halo abundance and weak lensing peak statistic can be used to probe the presence and properties of supervoids.

On the strength of a double Nobel prize winning experiment (Super)Kamiokande and an extremely successful long baseline neutrino programme, the third generation Water Cherenkov detector, Hyper-Kamiokande, is being developed by an international collaboration as a leading worldwide experiment based in Japan. The Hyper-Kamiokande detector will be hosted in the Tochibora mine, about 295 km away from the J-PARC proton accelerator research complex in Tokai, Japan. The currently existing accelerator will be steadily upgraded to reach a MW beam by the start of the experiment. A suite of near detectors will be vital to constrain the beam for neutrino oscillation measurements. A new cavern will be excavated at the Tochibora mine to host the detector. The experiment will be the largest underground water Cherenkov detector in the world and will be instrumented with new technology photosensors, faster and with higher quantum efficiency than the ones in Super-Kamiokande. The science that will be developed will be able to shape the future theoretical framework and generations of experiments. Hyper-Kamiokande will be able to measure with the highest precision the leptonic CP violation that could explain the baryon asymmetry in the Universe. The experiment also has a demonstrated excellent capability to search for proton decay, providing a significant improvement in discovery sensitivity over current searches for the proton lifetime. The atmospheric neutrinos will allow to determine the neutrino mass ordering and, together with the beam, able to precisely test the three-flavour neutrino oscillation paradigm and search for new phenomena. A strong astrophysical programme will be carried out at the experiment that will detect supernova neutrinos and will measure precisely solar neutrino oscillation.

The cosmic microwave background (CMB) has non-Gaussian features in the temperature fluctuations. An anomalous cold spot surrounded with a hot ring, called the Cold Spot is one of such features. If a large underdence region (supervoid) resides towards the Cold Spot, we would be able to detect a systematic shape distortion in the images of background source galaxies via weak lensing effect. In order to estimate the detectability of such signals, we used the data of $N$-body simulations to simulate full-sky ray-tracing of source galaxies. We searched for a most prominent underdense region using the simulated convergence maps smoothed at a scale of 20 degree and obtained tangential shears around it. The lensing signal expected in a concordant $\Lambda$CDM model can be detected at a signal-to-noise ratio $S/N\sim3$. If a supervoid with a radius of $\sim 200\,h^{-1}\,\textrm{Mpc}$ and a density contrast $\delta_0 \sim -0.3$ at the centre resides at a redshift $z\sim 0.2$, on-going and near-future weak gravitational lensing surveys would detect a lensing signal with $S/N\gtrsim4$ without resorting to stacking. From the tangential shear profile, we can obtain a constraint on the projected mass distribution of the supervoid.

By monitoring a large number of stars in the Local Group galaxies, we can detect nanolensing events by sub-lunar mass compact objects (SULCOs) such as primordial black holes (PBHs) and rogue (free-floating) dwarf planets in the Milky Way halo. In contarst to microlensing by stellar-mass objects, the finite-source size effect becomes important and the lensing time duration becomes shorter ($\sim 10^{1-4}\,\textrm{s}$). Using stars with $V<26$ in M33 as sources, for one-night observation, we would be able to detect $10^{3-4}$ nanolensing events caused by SULCOs in the Milky Way halo with a mass of $10^{-9}M_{\odot}$ to $10^{-7}M_{\odot}$ for sources with S/N$>5$ if SULCOs constitute all the dark matter components. Moreover, we expect $10^{1-2}$ events in which bright blue stars with S/N$>100$ are weakly amplified due to lensing by SULCOs with a mass range of $10^{-11}M_{\odot}$ to $10^{-9}M_{\odot}$. Thus the method would open a new window on SULCOs in the Milky Way halo that would otherwise not be observable.

We present monitoring campaign observations at optical and near-infrared (NIR) wavelengths for a radio-loud active galactic nucleus (AGN) at z=0.840, SDSS~J110006.07+442144.3 (hereafter, J1100+4421), which was identified during a flare phase in late February, 2014. The campaigns consist of three intensive observing runs from the discovery to March, 2015, mostly within the scheme of the OISTER collaboration. Optical-NIR light curves and simultaneous spectral energy distributions (SEDs) are obtained. Our measurements show the strongest brightening in March, 2015. We found that the optical-NIR SEDs of J1100+4421 show an almost steady shape despite the large and rapid intranight variability. This constant SED shape is confirmed to extend to $\sim5~\mu$m in the observed frame using the archival WISE data. Given the lack of absorption lines and the steep power-law spectrum of $\alpha_{\nu}\sim-1.4$, where $f_{\nu}\propto\nu^{\alpha_{\nu}}$, synchrotron radiation by a relativistic jet with no or small contributions from the host galaxy and the accretion disk seems most plausible as an optical-NIR emission mechanism. The steep optical-NIR spectral shape and the large amplitude of variability are consistent with this object being a low $\nu_{\rm{peak}}$ jet-dominated AGN. In addition, sub-arcsec resolution optical imaging data taken with Subaru Hyper Suprime-Cam does not show a clear extended component and the spatial scales are significantly smaller than the large extensions detected at radio wavelengths. The optical spectrum of a possible faint companion galaxy does not show any emission lines at the same redshift and hence a merging hypothesis for this AGN-related activity is not supported by our observations.

The dark matter search project by means of ultra high purity NaI(Tl) scintillator is now underdevelopment. An array of large volume NaI(Tl) detectors whose volume is 12.7 cm$\phi\times$12.7 cm is applied to search for dark matter signal. To remove radioactive impurities in NaI(Tl) crystal is one of the most important task to find small number of dark matter signals. We have developed high purity NaI(Tl) crystal which contains small amounts of radioactive impurities, $<4$ ppb of $^{nat}$K, 0.3 ppt of Th chain, 58 $\mu$Bq/kg of $^{226}$Ra and 30 $\mu$Bq/kg of $^{210}$Pb. Future prospects to search for dark matter by means of a large volume and high purity NaI(Tl) scintillator is discussed.

We investigate anomalous strong lens systems, particularly the effects of weak lensing by structures in the line of sight, in models with long-lived electrically charged massive particles (CHAMPs). In such models, matter density perturbations are suppressed through the acoustic damping and the flux ratio of lens systems are impacted, from which we can constrain the nature of CHAMPs. For this purpose, first we perform $N$-body simulations and develop a fitting formula to obtain non-linear matter power spectra in models where cold neutral dark matter and CHAMPs coexist in the early Universe. By using the observed anomalous quadruple lens samples, we obtained the constraints on the lifetime ($\tau_{\rm Ch}$) and the mass density fraction ($r_{\rm Ch}$) of CHAMPs. We show that, for $r_{\rm Ch}=1$, the lifetime is bounded as $\tau_{\rm Ch} < 0.96\,$yr (95% confidence level), while a longer lifetime $\tau_{\rm Ch} = 10\,$yr is allowed when $r_{\rm Ch} < 0.5$ at the 95% confidence level. Implications of our result for particle physics models are also discussed.

We report the $4 \, \sigma$ detection of a faint object with a flux of ~ 0.3 mJy, in the vicinity of the quadruply lensed QSO MG0414+0534 using the Atacama Large Millimeter/submillimeter array (ALMA) Band 7. The object is most probably a dusty dark dwarf galaxy, which has not been detected in either the optical, near-infrared (NIR) or radio (cm) bands. An anomaly in the flux ratio of the lensed images observed in Band 7 and the mid-infrared (MIR) band and the reddening of the QSO light color can be simultaneously explained if we consider the object as a lensing substructure with an ellipticity ~ 0.7 at a redshift of $0.5 \lesssim z \lesssim 1$. Using the best-fit lens models with three lenses, we find that the dark matter plus baryon mass associated with the object is $\sim 10^9\, M_{\odot}$, the dust mass is $\sim 10^7\,M_{\odot}$ and the linear size is $\gtrsim 5\,$kpc. Thus our findings suggest that the object is a dusty dark dwarf galaxy. A substantial portion of faint submillimeter galaxies (SMGs) in the universe may be attributed to such dark objects.

We present a search for low energy antineutrino events coincident with the gravitational wave events GW150914 and GW151226, and the candidate event LVT151012 using KamLAND, a kiloton-scale antineutrino detector. We find no inverse beta-decay neutrino events within $\pm 500$ seconds of either gravitational wave signal. This non-detection is used to constrain the electron antineutrino fluence and the luminosity of the astrophysical sources.

In the cosmic microwave background or galaxy density maps, missing fluctuations in masked regions can be reconstructed from fluctuations in the surrounding unmasked regions if the original fluctuations are sufficiently smooth. One reconstruction method involves applying a harmonic expansion iteratively to fluctuations in the unmasked region. In this paper, we discuss how well this reconstruction method can recover the original fluctuations depending on the prior of fluctuations and property of the masked region. The reconstruction method is formulated with an asymptotic expansion in terms of the size of mask for a fixed iteration number. The reconstruction accuracy depends on the mask size, the spectrum of the underlying density fluctuations, the scales of the fluctuations to be reconstructed and the number of iterations. For Gaussian fluctuations with the Harrison--Zel'dovich spectrum, the reconstruction method provides more accurate restoration than naive methods based on brute--forth matrix inversion or the singular value decomposition. We also demonstrate that an isotropic non-Gaussian prior does not change the results but an anisotropic non-Gaussian prior can yield a higher reconstruction accuracy compared to the Gaussian prior case.

We explore the origin of flux ratio anomaly in quadruple lens systems. Using a semi-analytic method based on $N$-body simulations, we estimate the effect of possible magnification perturbation caused by subhaloes with a mass scale of <$\sim 10^9\,h^{-1} \textrm{M}_\odot$ in lensing galaxy haloes. Taking into account astrometric shifts, assuming that the primary lens is described by a singular isothermal ellipsoid, the expected change to the flux ratios per a multiply lensed image is just a few percent and the mean of the expected convergence perturbation at the effective Einstein radius of the lensing galaxy halo is $ 0.003$, corresponding to the mean of the ratio of a projected dark matter mass fraction in subhaloes at the effective Einstein radius $0.006$. In contrast, the expected change to the flux ratio caused by line-of-sight structures is typically $\sim 10$ percent and the mean of the convergence perturbation is $0.008$, corresponding to $0.017$. The contribution of magnification perturbation caused by subhaloes is $\sim 40$ percent of the total at a source redshift $z_S= 0.7$ and decreases monotonically in $z_S$ to $\sim 20$ percent at $z_S= 3.6$. Assuming statistical isotropy, the convergence perturbation estimated from observed 11 quadruple lens systems has a positive correlation with the source redshift $z_S$, which is much stronger than that with the lens redshift $z_L$. This feature also supports an idea that the flux ratio anomaly is caused mainly by line-of-sight structures rather than subhaloes. We also discuss about a possible imprint of line-of-sight structures in demagnification of minimum images due to locally underdense structures in the line of sight.

A high purity and large volume NaI(Tl) scintillator was developed to search for cosmic dark matter. The required densities of radioactive impurities (RIs) such as U-chain, Th-chain are less than a few ppt to establish high sensitivity to dark matter. The impurity of RIs were effectively reduced by selecting raw materials of crucible and by performing chemical reduction of lead ion in NaI raw powder. The impurity of $^{226}$Ra was reduced less than 100 $\mu$Bq/kg in NaI(Tl) crystal. It should be remarked that the impurity of $^{210}$Pb, which is difficult to reduce, is effectively reduced by chemical processing of NaI raw powder down to less than 30 $\mu$Bq/kg. The expected sensitivity to cosmic dark matter by using 250 kg of the high purity and large volume NaI(Tl) scintillator (PICO-LON; Pure Inorganic Crystal Observatory for LOw-background Neutr(al)ino) is 7$\times$10$^{-45}$ cm$^{2}$ for 50 GeV$/c^{2}$ WIMPs.

Recently it has been claimed that the warm dark matter (WDM) model cannot at the same time reproduce the observed Lyman-{\alpha} forests in distant quasar spectra and solve the small-scale issues in the cold dark matter (CDM) model. As an alternative candidate, it was shown that the mixed dark matter (MDM) model that consists of WDM and CDM can satisfy the constraint from Lyman-{\alpha} forests and account for the "missing satellite problem" as well as the reported 3.5 keV anomalous X-ray line. We investigate observational constraints on the MDM model using strong gravitational lenses. We first develop a fitting formula for the nonlinear power spectra in the MDM model by performing N-body simulations and estimate the expected perturbations caused by line-of-sight structures in four quadruply lensed quasars that show anomaly in the flux ratios. Our analysis indicates that the MDM model is compatible with the observed anomaly if the mass fraction of the warm component is smaller than 0.47 at the 95% confidence level. The MDM explanation to the anomalous X-ray line and the small-scale issues is still viable even after this constraint is taken into account.

We present an analysis of the ALMA long baseline science verification data of the gravitational lens system SDP.81. We fit the positions of the brightest clumps at redshift z=3.042 and a possible AGN component of the lensing galaxy at redshift z=0.2999 in the band 7 continuum image using a canonical lens model, a singular isothermal ellipsoid plus an external shear. Then, we measure the ratio of fluxes in some apertures at the source plane where the lensed images are inversely mapped. We find that the aperture flux ratios of band 7 continuum image are perturbed by 10-20 percent with a significance at 2 ~ 3 sigma level. Moreover, we measure the astrometric shifts of multiply lensed images near the caustic using the CO(8-7) line. Using a lens model best-fitted to the band 7 continuum image, we reconstruct the source image of the CO(8-7) line by taking linear combination of inverted quadruply lensed images. At the 50th channel (rest-frame velocity 28.6 km/s) of the CO(8-7) line, we find an imprint of astrometric shifts of the order of 0.01 arcsec in the source image. Based on a semi-analytic calculation, we find that the observed anomalous flux ratios and the astrometric shifts can be explained by intergalactic dark structures in the line of sight. A compensated homogeneous spherical clump with a mean surface mass density of the order of 10^8 solar mass h^-1 arcsec^-2 can explain the observed anomaly and astrometric shifts simultaneously.

In the late stages of nuclear burning for massive stars ($M>8~M_{\sun}$), the production of neutrino-antineutrino pairs through various processes becomes the dominant stellar cooling mechanism. As the star evolves, the energy of these neutrinos increases and in the days preceding the supernova a significant fraction of emitted electron anti-neutrinos exceeds the energy threshold for inverse beta decay on free hydrogen. This is the golden channel for liquid scintillator detectors because the coincidence signature allows for significant reductions in background signals. We find that the kiloton-scale liquid scintillator detector KamLAND can detect these pre-supernova neutrinos from a star with a mass of $25~M_{\sun}$ at a distance less than 690~pc with 3$\sigma$ significance before the supernova. This limit is dependent on the neutrino mass ordering and background levels. KamLAND takes data continuously and can provide a supernova alert to the community.

The PICO-LON project aims at search for cold dark matter by means of highly radio-pure and large volume NaI(Tl) scintillator. The NaI powder was purifed by chemical processing to remove lead isotopes and selecting a high purity graphite crucible. The concentrations of radioactive impurities of $^{226}$Ra and $^{228}$Th were effectively reduced to 58$\pm$4 $\mu$Bq/kg and 1.5$\pm$1.9 $\mu$Bq/kg, respectively. It should be remarked that the concentration of $^{210}$Pb, which is crucial for the sensitivity to dark matter, was reduced to 24$\pm$2 $\mu$Bq/kg. The total background rate at 10 keVee was as low as 8 keV$^{-1}$kg$^{-1}$day$^{-1}$, which was sufficiently low to search for dark matter. Further purification of NaI(Tl) ingot and future prospect of PICO-LON project is discussed.

We report a measurement of the neutrino-electron elastic scattering rate of 862 keV 7Be solar neutrinos based on a 165.4 kton-day exposure of KamLAND. The observed rate is 582 +/- 90 (kton-day)^-1, which corresponds to a 862 keV 7Be solar neutrino flux of (3.26 +/- 0.50) x 10^9 cm^-2s^-1, assuming a pure electron flavor flux. Comparing this flux with the standard solar model prediction and further assuming three flavor mixing, a nu_e survival probability of 0.66 +/- 0.14 is determined from the KamLAND data. Utilizing a global three flavor oscillation analysis, we obtain a total 7Be solar neutrino flux of (5.82 +/- 0.98) x 10^9 cm^-2s^-1, which is consistent with the standard solar model predictions.

We search for electron anti-neutrinos ($\overline{\nu}_e$) from long and short-duration gamma-ray bursts~(GRBs) using data taken by the KamLAND detector from August 2002 to June 2013. No statistically significant excess over the background level is found. We place the tightest upper limits on $\overline{\nu}_e$ fluence from GRBs below 7 MeV and place first constraints on the relation between $\overline{\nu}_e$ luminosity and effective temperature.

We investigate the weak lensing effect by line-of-sight structures with a surface mass density of <~10^8 solar mass/arcsec^2 in QSO-galaxy quadruple lens systems. Using high-resolution N-body simulations in warm dark matter (WDM) models and observed four quadruple lenses that show anomalies in the flux ratios, we obtain constraints on the mass of thermal WDM, m_WDM>= 1.3keV(95%CL) assuming that the density of the primary lens is described by a singular isothermal ellipsoid (SIE). The obtained constraint is consistent with those from Lyman-$\alpha$ forests and the number counts of high-redshift galaxies at z>4. Our results show that WDM with a free-streaming comoving wavenumber k_{fs} <= 27 h/Mpc is disfavored as the major component of cosmological density at redshifts 0.5 <~ z <~ 4 provided that the SIE models describe the gravitational potentials of the primary lenses correctly.

We explore the weak lensing effects by ministructures in the line-of-sight in a quadruply lensed quasar MG0414+0534 that shows an anomaly in the flux-ratios. We find that the observed flux-ratio anomaly can be explained by a presence of either a minifilament or a minivoid in the line-of-sight with a surface mass density of the order of 10^(8-9) h^(-1) solar mass /arcsec^2 without taking into account any subhalos in the lensing galaxy. The astrometric perturbation by a possible minifilament/minivoid is <~ 0.001 arcsec and the amplitudes of convergence perturbations due to these perturbers are ~ 0.004-0.008 at the place of an image that shows anomaly. In order to discriminate models with the line-of-sight ministructures from those with a subhalo(s) in the lensing galaxy, we need to precisely measure the projected convergence and shear around the lensing galaxy. The differential magnification effect could break the model degeneracy if the source size is > ~100 pc. Observation at the submillimeter band using interferometers will enable us to determine the origin of anomalies in the flux ratios.

KamLAND-PICO project aims to search for WIMPs dark matter by means of NaI(Tl) scintillator. To investigate the WIMPs candidate whose cross section is as small as $10^{-9}$ pb, a pure NaI(Tl) crystal was developed by chemical processing and taking care of surroundings. The concentration of U and Th chain was reduced to $5.4\pm0.9$ ppt and $3.3\pm2.2$ ppt, respectively. It should be remarked that the concentration of $^{210}$Pb which was difficult to reduce reached to the high purity as $58\pm26$ $\mu$Bq/kg.

We investigate the weak lensing effects of line-of-sight structures on quadruple images in quasar-galaxy strong lens systems based on N-body and ray-tracing simulations that can resolve halos with a mass of 10^5 solar mass. The intervening halos and voids disturb the magnification ratios of lensed images as well as their relative positions due to lensing. The magnification ratios typically change by O(10%) when the shifts of relative angular positions of lensed images are constrained to <0.004 arcsec. The constrained amplitudes of projected density perturbations due to line-of-sight structures are O(10^8) solar mass per arcsec^2. These results are consistent with our new analytical estimate based on the two-point correlation of density fluctuations. The observed mid-infrared (MIR) flux ratios for 6 quasar-galaxy lens systems with quadruple images agree well with the numerically estimated values without taking into account of subhalos residing in the lensing galaxies. We find that the constrained mean amplitudes of projected density perturbations in the line-of-sight are negative, which suggests that the fluxes of lensed images are perturbed mainly by minivoids and minihalos in underdense regions. We derive a new fitting formula for estimating the probability distribution function of magnification perturbation. We also find that the mean amplitude of magnification perturbation roughly equals the standard deviation regardless of the model parameters.

We study the effect of global topology of the spatial geometry on the cosmic microwave background (CMB) for closed flat and closed hyperbolic models in which the spatial hypersurface is multiply connected. If the CMB temperature fluctuations were entirely produced at the last scattering, then the large-angle fluctuations would be much suppressed in comparison with the simply connected counterparts which is at variance with the observational data. However, as we shall show in this thesis, for low matter density models the observational constraints are less stringent since a large amount of large-angle fluctuations could be produced at late times. On the other hand, a slight suppression in large-angle temperature correlations in such models explains rather naturally the observed anomalously low quadrupole which is incompatible with the prediction of the "standard" Friedmann-Robertson-Walker-Lemaitre models. Interestingly, moreover, the development in the astronomical observation technology has made it possible to directly explore the imprint of the non-trivial topology by looking for identical objects so called "ghosts" in wide separated directions. For the CMB temperature fluctuations identical patterns would appear on a pair of circles in the sky. Another interesting feature is the non-Gaussianity in the temperature fluctuations. Inhomogeneous and anisotropic Gaussian fluctuations for a particular choice of position and orientation are regarded as non-Gaussian fluctuations for a homogeneous and isotropic ensemble.

A galaxy halo may contain a large number of intermediate mass black holes (IMBHs) with masses in the range of 10^{2-6} solar mass. We propose to directly detect these IMBHs by observing multiply imaged QSO-galaxy or galaxy-galaxy strong lens systems in the submillimeter bands with high angular resolution. The silhouette of an IMBH in the lensing galaxy halo would appear as either a monopole-like or a dipole-like variation at the scale of the Einstein radius against the Einstein ring of the dust-emitting region surrounding the QSO. We use a particle tagging technique to dynamically populate a Milky Way-sized dark matter halo with black holes, and show that the surface mass density and number density of IMBHs have power-law dependences on the distance from the center of the host halo if smoothed on a scale of ~ 1 kpc. Most of the black holes orbiting close to the center are freely roaming as they have lost their dark matter hosts during infall due to tidal stripping. Next generation submillimeter telescopes with high angular resolution (< 0.3 mas) will be capable of directly mapping such off-nuclear freely roaming IMBHs with a mass of ~ 10^6 solar mass in a lensing galaxy that harbours a O(10^9) solar mass supermassive black hole in its nucleus.

We explore the weak lensing effect by line-of-sight halos and sub-halos with a mass of M < 10^7 solar mass in QSO-galaxy strong lens systems with quadruple images in a concordant LCDM universe. Using a polynomially fitted non-linear power spectrum P(k) obtained from N-body simulations that can resolve halos with a mass of M ~ 10^5 solar mass, or structures with a comoving wavenumber of k ~ 3*10^2 h/Mpc, we find that the ratio of magnification perturbation due to intervening halos to that of a primary lens is typically ~10 per cent and the predicted values agree well with the estimated values for 6 observed QSO-galaxy lens systems with quadruple images in the mid-infrared band without considering the effects of substructures inside a primary lens. We also find that the estimated amplitudes of convergence perturbation for the 6 lenses increase with the source redshift as predicted by theoretical models. Using an extrapolated matter power spectrum, we demonstrate that small halos or sub-halos in the line-of-sight with a mass of M=10^3-10^7 solar mass, or structures with a comoving wavenumber of k=3*10^2-10^4 h/Mpc can significantly affect the magnification ratios of the lensed images. Flux ratio anomalies in QSO-galaxy strong lens systems offer us a unique probe into clustering property of mini-halos with a mass of M < 10^6 solar mass.

In a concordant $\Lambda$ Cold Dark Matter ($\Lambda$CDM) model, large-angle Cosmic Microwave Background (CMB) temperature anisotropy due to linear perturbations in the local universe is not negligible. We explore a possible role of an underdense region (void) that may cause an anomalous Cold Spot (CS) in the CMB sky. Although the observed anomalous cold region with a surrounding hot ring can be produced by an underdense region surrounded by a massive wall, a decrement in the CMB temperature in the line-of-sight is suppressed because of blueshift of CMB photons that pass the wall. Therefore, undercompensated models give better agreement with the observed data in comparison with overcompensated or compensated models. We find that it is likely that $\sim$90 per cent of the CMB fluctuation is generated due to an overdense region surrounded by an underdense region at the last scattering surface, and the remaining $\sim 10$ per cent is produced due to a single spherical underdense region with a radius $r\sim 6\times 10^2 h^{-1}$Mpc and a density contrast $\delta_m\sim -0.009$ ($2 \sigma$) at redshift $z\sim 1$ in the line-of-sight to the CS. The probability of chance alignment of such two structures is $\sim 0.7$ per cent if the perturbation with an underdense region at $z\sim 1$ is moderately undercompensated.

We propose the Hyper-Kamiokande (Hyper-K) detector as a next generation underground water Cherenkov detector. It will serve as a far detector of a long baseline neutrino oscillation experiment envisioned for the upgraded J-PARC, and as a detector capable of observing -- far beyond the sensitivity of the Super-Kamiokande (Super-K) detector -- proton decays, atmospheric neutrinos, and neutrinos from astronomical origins. The baseline design of Hyper-K is based on the highly successful Super-K, taking full advantage of a well-proven technology. (to be continued)

We report a measurement of the neutrino-electron elastic scattering rate from 8B solar neutrinos based on a 123 kton-day exposure of KamLAND. The background-subtracted electron recoil rate, above a 5.5 MeV analysis threshold is 1.49+/-0.14(stat)+/-0.17(syst) events per kton-day. Interpreted as due to a pure electron flavor flux with a 8B neutrino spectrum, this corresponds to a spectrum integrated flux of 2.77+/-0.26(stat)+/-0.32(syst) x 10^6 cm^-2s^-1. The analysis threshold is driven by 208Tl present in the liquid scintillator, and the main source of systematic uncertainty is due to background from cosmogenic 11Be. The measured rate is consistent with existing measurements and with Standard Solar Model predictions which include matter enhanced neutrino oscillation.

Recent measurements of hot and cold spots on the cosmic microwave background (CMB) sky suggest a presence of super-structures on (>100 h^{-1}Mpc) scales. We develop a new formalism to estimate the expected amplitude of temperature fluctuations due to the integrated Sachs-Wolfe (ISW) effect from prominent quasi-linear structures. Applying the developed tools to the observed ISW signals from voids and clusters in catalogs of galaxies at redshifts z<1, we find that they indeed imply a presence of quasi-linear super-structures with a comoving radius 100~300 h^{-1}Mpc and a density contrast ~O(0.1). We find that the observed ISW signals are at odd with the concordant \Lambda cold dark matter (CDM) model that predicts Gaussian primordial perturbations at equal to or larger than 3 sigma level. We also confirm that the mean temperature around the CMB cold spot in the southern Galactic hemisphere filtered by a compensating top-hat filter deviates from a mean value at ~3 sigma level, implying that a quasi-linear supervoid or an underdensity region surrounded by a massive wall may reside at low redshifts z<0.3 and the actual angular size (16^\circ-17^\circ) may be larger than the apparent size (4^\circ-10^\circ) discussed in literature. Possible solutions are briefly discussed.

Because accretion and merger shocks in clusters of galaxies may accelerate particles to high energies, clusters are candidate sites for the origin of ultra-high-energy (UHE) cosmic-rays. A prediction was presented for gamma-ray emission from a cluster of galaxies at a detectable level with the current generation of imaging atmospheric Cherenkov telescopes. The gamma-ray emission was produced via inverse Compton upscattering of cosmic microwave background (CMB) photons by electron-positron pairs generated by collisions of UHE cosmic rays in the cluster. We observed two clusters of galaxies, Abell 3667 and Abell 4038, searching for very-high-energy gamma-ray emission with the CANGAROO-III atmospheric Cherenkov telescope system in 2006. The analysis showed no significant excess around these clusters, yielding upper limits on the gamma-ray emission. From a comparison of the upper limit for the north-west radio relic region of Abell 3667 with a model prediction, we derive a lower limit for the magnetic field of the region of ~0.1 micro G. This shows the potential of gamma-ray observations in studies of the cluster environment. We also discuss the flux upper limit from cluster center regions using a model of gamma-ray emission from neutral pions produced in hadronic collisions of cosmic-ray protons with the intra-cluster medium (ICM). The derived upper limit of the cosmic-ray energy density within this framework is an order of magnitude higher than that of our Galaxy.

Observation by the CANGAROO-III stereoscopic system of the Imaging Cherenkov Telescope has detected extended emission of TeV gamma rays in the vicinity of the pulsar PSR B1706$-$44. The strength of the signal observed as gamma-ray-like events varies when we apply different ways of emulating background events. The reason for such uncertainties is argued in relevance to gamma-rays embedded in the "off-source data", that is, unknown sources and diffuse emission in the Galactic plane, namely, the existence of a complex structure of TeV gamma-ray emission around PSR B1706$-$44.

We present mid-infrared imaging at 11.7 $\mu$m for the quadruple lens systems, MG0414+0534 and Q2237+030, using the cooled mid-infrared camera and spectrometer (COMICS) attached on the Subaru telescope. MG0414+0534 is characterized by a bright pair of lensed images (A1, A2) and their optical flux ratio A2/A1 deviates significantly from the prediction of a smooth lens model. Q2237+030 is comprised of four lensed images, which are significantly affected by microlensing in a foreground lensing galaxy. Our mid-infrared observations of these lensed images have revealed that the mid-infrared flux ratio for A2/A1 of MG0414+0534 is nearly unity (0.90 +- 0.04). We find that this flux ratio is systematically small, at 4 to 5 sigma level, compared with the prediction of a best smooth lens model (1.09) represented by a singular isothermal ellipsoid and external shear, thereby suggesting the presence of substructures to explain our observational result. In contrast, for Q2237+030, our high signal-to-noise observation indicates that the mid-infrared flux ratios between all the four images of Q2237+030 are virtually consistent with the prediction of a smooth lens model. Based on the size estimate of the dust torus surrounding the nuclei of these QSOs, we set limits on the mass of a substructure in these lens systems, which can cause anomalies in the flux ratios. For MG0414+0534, since the required mass of a substructure inside its Einstein radius is > 360 M_sun, millilensing by a CDM substructure is most likely. If it is modeled as a singular isothermal sphere, the mass inside radius of 100 pc is given as > 1.0x10^5 M_sun. For Q2237+030, there is no significant evidence of millilensing, so the reported anomalous flux ratios in shorter wavelengths are entirely caused by microlensing by stars (abridged).

Recent observational data of supernovae indicate that we may live in an underdense region, which challenges the Copernican principle. We show that the integrated Sachs-Wolfe (ISW) effect is an excellent discriminator between anti-Copernican inhomogeneous models and the standard Copernican models. As a reference model, we consider an anti-Copernican inhomogeneous model that consists of two inner negatively curved underdense regions and an outer flat Einstein-de Sitter region. We assume that these regions are connected by two thin-walls at redshifts z = 0.067 and z=0.45. In the inner two regions, the first-order ISW effect is dominant and comparable to that in the concordant flat-Lambda models. In the outer Einstein-de Sitter region, the first-order ISW effect vanishes but the second-order ISW effect plays a dominant role, while the first-order ISW effect is dominant in the flat-Lambda models at moderate redshifts. This difference can discrimate the anti-Copernican models from the concordant flat-Lambda model. At high redshits, the second-order ISW effect is dominant both in our inhomogeneous model and the concordant model. In the outer region, moreover, the ISW effect due to large-scale density perturbations with a present matter density contrast much less than 0.37 is negligible, while the effect due to small-scale density perturbations (such as clusters of galaxies, superclusters and voids) with matter density contrast much larger than 0.37 would generate anisotropies which are larger than those generated by the ISW effect in the concordant model.

We study the cosmic microwave background (CMB) anisotropy due to spherically symmetric nonlinear structures in flat universes with dust and a cosmological constant. By modeling a time-evolving spherical compensated void/lump by Lemaitre-Tolman-Bondi spacetimes, we numerically solve the null geodesic equations with the Einstein equations. We find that a nonlinear void redshifts the CMB photons that pass through it regardless of the distance to it. In contrast, a nonlinear lump blueshifts (or redshifts) the CMB photons if it is located near (or sufficiently far from) us. The present analysis comprehensively covers previous works based on a thin-shell approximation and a linear/second order perturbation method and the effects of shell thickness and full nonlinearity. Our results indicate that, if quasi-linear and large ($>100$Mpc) voids/lumps would exist, they could be observed as cold or hot spots with temperature variance $>10^{-5}$K in the CMB sky.

We develop a new interpolation scheme, based on harmonic inpainting, for reconstructing the cosmic microwave background temperature data within the Galaxy mask from the data outside the mask. We find that, for scale-invariant isotropic random Gaussian fluctuations, the developed algorithm reduces the errors in the reconstructed map for the odd-parity modes significantly for azimuthally symmetric masks with constant galactic latitudes. For a more realistic Galaxy mask, we find a modest improvement in the even-parity modes as well.

We study second order gravitational effects of local inhomogeneities on the cosmic microwave background radiation in flat universes with matter and a cosmological constant $\Lambda$. We find that the general relativistic correction to the Newtonian approximation is negligible at second order provided that the size of the inhomogeneous region is sufficiently smaller than the horizon scale. For a spherically symmetric top-hat type quasi-linear perturbation, the first order temperature fluctuation corresponding to the linear integrated Sachs-Wolfe (ISW) effect is enhanced(suppressed) by the second order one for a compensated void(lump). As a function of redshift of the local inhomogeneity, the second order temperature fluctuations due to evolution of the gravitational potential have a peak before the matter-$\Lambda$ equality epoch for a fixed comoving size and a density contrast. The second order gravitational effects from local quasi-linear inhomogeneities at a redshift z~1 may significantly affect the cosmic microwave background.

To set useful limits on the abundance of small-scale dark matter halos (subhalos) in a galaxy scale, we have carried out mid-infrared imaging and integral-field spectroscopy for a sample of quadruple lens systems showing anomalous flux ratios. These observations using Subaru have been successful for distinguishing millilensing by subhalos from microlensing by stars. Current status for our lensing analysis of dark matter substructure is reported.

Hierarchical clustering models of cold dark matter (CDM) predict that about 5% - 10% of a galaxy-sized halo with mass ~ 10^12 solar masses (M_sun) resides in substructures (CDM subhalos) with masses <= 10^8 M_sun. To directly identify such substructures, we propose to observe radio continuum emission from multiply imaged QSOs using VSOP-2 with a high angular resolution.

It has been argued that the large-angle cosmic microwave background anisotropy has anomalies at 3-sigma level. We review various proposed ideas to explain the origin of the anomalies and discuss how we can constrain the proposed models using future observational data.

We explore the large angular scale temperature anisotropies in the cosmic microwave background (CMB) due to homogeneous local dust-filled voids in a flat Friedmann-Robertson-Walker universe with a cosmological constant. In comparison with the equivalent dust-filled void model in the Einstein-de Sitter background, we find that the anisotropy for compensated asymptotically expanding local voids can be larger because second-order effects enhance the linear integrated Sachs-Wolfe (ISW) effect. However, for local voids that expand sufficiently faster than the asymptotic velocity of the wall, the second-order effect can suppress the fluctuation due to the linear ISW effect. A pair of quasi-linear compensated asymptotic local voids with radius (2-3)*10^2 ~h^{-1} Mpc and a matter density contrast ~-0.3 can be observed as cold spots with a temperature anisotropy Delta T/T~O(10^{-5}) that might help explain the observed large-angle CMB anomalies. We predict that the associated anisotropy in the local Hubble constant in the direction of the voids could be as large as a few percent.

We have observed the quadruply lensed quasar 1RXS J1131-1231 with the integral field spectrograph mode of the Kyoto Tridimensional Spectrograph II mounted on the Subaru telescope. Its field of view has covered simultaneously the three brighter lensed images A, B, and C, which are known to exhibit anomalous flux ratios in their continuum emission. We have found that the [OIII] line flux ratios among these lensed images are consistent with those predicted by smooth-lens models. The absence of both microlensing and millilensing effects on this [OIII] narrow line region sets important limits on the mass of any substructures along the line of sight, which is expressed as M_E < 10^5 M_solar for the mass inside an Einstein radius. In contrast, the H_beta line emission, which originates from the broad line region, shows an anomaly in the flux ratio between images B and C, i.e., a factor two smaller C/B ratio than predicted by smooth-lens models. The ratio of A/B in the H_beta line is well reproduced. We show that the anomalous C/B ratio for the H_beta line is caused most likely by micro/milli-lensing of image C. This is because other effects, such as the differential dust extinction and/or arrival time difference between images B and C, or the simultaneous lensing of another pair of images A and B, are all unlikely. In addition, we have found that the broad H_beta line of image A shows a slight asymmetry in its profile compared with those in the other images, which suggests the presence of a small microlensing effect on this line emitting region of image A.

We explore the large angular scale temperature anisotropies in the cosmic microwave background due to expanding homogeneous local voids at redshift z~1. A compensated spherically symmetric homogeneous dust-filled void with radius \~3*10^2 h^{-1}Mpc, and density contrast ~-0.3 can be observed as a cold spot with a temperature anisotropy -1*10^{-5} surrounded by a slightly hotter ring. We find that a pair of these circular cold spots separated by ~50 degree can account both for the planarity of the octopole and for the alignment between the quadrupole and the octopole in the cosmic microwave background (CMB) anisotropy. The cold spot in the Galactic southern hemisphere which is anomalous at the ~3sigma level can be explained by such a large void at z~1. The observed north-south asymmetry in the large-angle CMB power can be attributed to the asymmetric distribution of these local voids between the two hemispheres. The statistical significance of the low quadrupole is further reduced in this interpretation of the large angular scale CMB anomalies.

A galaxy-sized halo may contain a large number of intermediate mass (10^{2-4} solar mass) compact objects (IMCOs), which can be intermediate mass black holes (IMBHs) or the CDM subhalos. We propose to directly detect the IMBHs by observing multiply imaged QSO-galaxy lens systems with a high angular resolution (~ 0.03mas), which would be achieved by the next-VLBI space missions. The silhouette of the IMBHs would appear as an either monopole-like or dipole-like variation at the scale of the Einstein radius against the QSO jets. As a byproduct, we can also directly detect the 10^{4-5} solar mass CDM subhalos. From a measurement of the local distortion in the surface brightness of the QSO jet, we can make a distinction between a point mass (corresponding to an IMBH) and an extended structure (corresponding to a CDM subhalo). It would be a unique probe of the IMCOs whose nature has been under the veil of mistery.

We investigate the extended source size effects on gravitational lensing in which a lens consists of a smooth potential and small mass clumps (``substructure lensing''). We first consider a lens model that consists of a clump modeled as a singular isothermal sphere (SIS) and a primary lens modeled as an external background shear and convergence. For this simple model, we derive analytic formulae for (de)magnification of circularly symmetric top-hat sources with three types of parity for their lensed images, namely, positive, negative, and doubly negative parities. Provided that the source size is sufficiently larger than the Einstein radius of the SIS, we find that in the positive (doubly negative) parity case, an extended source is always magnified (demagnified) in comparison with the unperturbed macrolens system, whereas in the negative parity case, the (de)magnification effect, which depends on the sign of convergence minus unity is weaker than those in other parities. It is shown that a measurement of the distortion pattern in a multiply lensed image enables us to break the degeneracy between the lensing effects of clump mass and those of clump distance if lensing parameters of the relevant macrolens model are determined from the position and flux of multiple images. We also show that an actual density profile of a clump can be directly measured by analyzing the ``fine structure'' in a multiply lensed image within the Einstein radius of the clump.

The cold dark matter (CDM) structure formation model predicts that about 5-10 percent of a typical galactic halo of mass $\sim 10^{12} \ms$ is in substructures with masses $\lesssim 10^8 \ms$. To directly detect such substructures, we propose to observe dust continuum emission from a strongly lensed QSO-host galaxy using a large submillimeter interferometer array with a high angular resolution of $\sim 0.01$arcsec such as the planned Atacama Large Submillimeter Array (ALMA). To assess their observational feasibility, we numerically simulate millilensing of an extended circular source by a CDM substructure modeled as a tidally truncated singular isothermal sphere (SIS) embedded in a typical QSO-galaxy lens system, B1422+231, modeled as a singular isothermal ellipsoid (SIE) with an external constant shear and a constant convergence. Assuming an angular resolution of 0.01arcsec, we find that the angular positions of $\sim 10^8 \ms$ substructures at several kpc from the center of the macrolens halo can be directly measured if the size of the dust continuum emission region and the gradient of the surface brightness at the position of the perturber are sufficiently large. From the astrometric shift on a scale of a few times $10~$mas of an image perturbed by a subhalo with respect to an unperturbed macrolensed image, we can break the degeneracy between subhalo mass and distance provided that macrolensing parameters are determined from positions and fluxes of multiple images.

We present mid-infrared imaging at 11.7 mu m for the quadruple lens systems, PG1115+080 and B1422+231, using the cooled mid-infrared camera and spectrometer (COMICS) attached on the Subaru telescope. These lensed QSOs are characterized by their anomalous optical and radio flux ratios as obtained for (A1, A2) images of PG1115+080 and (A, B, C) images of B1422+231, respectively, i.e., such flux ratios are hardly reproduced by lens models with smooth mass distribution. Our mid-infrared observations for these images have revealed that the mid-infrared flux ratio A2/A1 of PG1115+080 is virtually consistent with smooth lens models (but inconsistent with the optical flux ratio), whereas for B1422+231, the mid-infrared flux ratios among (A, B, C) are in good agreement with the radio flux ratios. We also identify a clear infrared bump in the spectral energy distributions of these QSOs, thereby indicating that the observed mid-infrared fluxes originate from a hot dust torus around a QSO nucleus. Based on the size estimate of the dust torus, we place limits on the mass of a substructure in these lens systems, causing the anomalous optical or radio flux ratios. For PG1115+080, the mass of a substructure inside an Einstein radius, M_E, is < 16 Msun, corresponding to either a star or a low-mass CDM subhalo having the mass of M_{100}^{SIS} < 2.2 * 10^4 Msun inside radius of 100 pc if modeled as a singular isothermal sphere (SIS). For B1422+231, we obtain M_E > 209 Msun, indicating that a CDM subhalo is more likely, having the mass of M_{100}^{SIS} > 7.4 * 10^4 Msun

We investigate the expected gravitational wave emission from coalescing supermassive black hole (SMBH) binaries resulting from mergers of their host galaxies. We employ a semi-analytic model of galaxy and quasar formation based on the hierarchical clustering scenario to estimate the amplitude of the expected stochastic gravitational wave background owing to inspiraling SMBH binaries and bursts rates owing to the SMBH binary coalescence events. We find that the characteristic strain amplitude of the background radiation is $h_c(f) \sim 10^{-16} (f/1 \mu {\rm Hz})^{-2/3}$ for $f \lesssim 1 \mu {\rm Hz}$. The main contribution to the total strain amplitude of the background radiation comes from SMBH coalescence events at $0<z<1$. We also find that a future space-based gravitational wave interferometer such as the planned \textit{Laser Interferometer Space Antenna} ({\sl LISA}) might detect intense gravitational wave bursts associated with coalescence of SMBH binaries with total mass $M_{\rm tot} < 10^7 M_{\odot}$ at $z \gtrsim 2$ at a rate $ \sim 1.0 {\rm yr}^{-1}$. Our model predicts that burst signals with a larger amplitude $h_{\rm burst} \sim 10^{-15}$ correspond to coalescence events of massive SMBH binary with total mass $M_{\rm tot} \sim 10^8 M_{\odot}$ at low redshift $ z \lesssim 1$ at a rate $ \sim 0.1 {\rm yr}^{-1}$ whereas those with a smaller amplitude $h_{\rm burst} \sim 10^{-17}$ correspond to coalescence events of less massive SMBH binary with total mass $M_{\rm tot} \sim 10^6 M_{\odot}$ at high redshift $ z \gtrsim 3$.

We present the results of indirect searches for Weakly Interacting Massive Particles (WIMPs) with 1679.6 live days of data from the Super-Kamiokande detector using neutrino-induced upward through-going muons. The search is performed by looking for an excess of high energy muon neutrinos from WIMP annihilations in the Sun, the core of the Earth, and the Galactic Center, as compared to the number expected from the atmospheric neutrino background. No statistically significant excess was seen. We calculate flux limits in various angular cones around each of the above celestial objects. We obtain conservative model-independent upper limits on WIMP-nucleon cross-section as a function of WIMP mass and compare these results with the corresponding results from direct dark matter detection experiments.

We investigate the expected gravitational wave emission from coalescing supermassive black hole (SMBH) binaries resulting from mergers of their host galaxies. When galaxies merge, the SMBHs in the host galaxies sink to the center of the new merged galaxy and form a binary system. We employ a semi-analytic model of galaxy and quasar formation based on the hierarchical clustering scenario to estimate the amplitude of the expected stochastic gravitational wave background owing to inspiraling SMBH binaries and bursts owing to the SMBH binary coalescence events. We find that the characteristic strain amplitude of the background radiation is $h_c(f) \sim 10^{-16} (f/1 \mu {\rm Hz})^{-2/3}$ for $f \lesssim 1 \mu {\rm Hz}$ just below the detection limit from measurements of the pulsar timing provided that SMBHs coalesce simultaneously when host galaxies merge. The main contribution to the total strain amplitude of the background radiation comes from SMBH coalescence events at $0<z<1$. We also find that a future space-based gravitational wave interferometer such as the planned \textit{Laser Interferometer Space Antenna} ({\sl LISA}) might detect intense gravitational wave bursts associated with coalescence of SMBH binaries with total mass $M_{\rm tot} < 10^7 M_{\odot}$ at $z \gtrsim 2$ at a rate $ \sim 1.0 {\rm yr}^{-1}$. Our model predicts that burst signals with a larger amplitude $h_{\rm burst} \sim 10^{-15}$ correspond to coalescence events of massive SMBH binary with total mass $M_{\rm tot} \sim 10^8 M_{\odot}$ at low redshift $ z \lesssim 1$ at a rate $ \sim 0.1 {\rm yr}^{-1}$ whereas those with a smaller amplitude $h_{\rm burst} \sim 10^{-17}$ correspond to coalescence events of less massive SMBH binary with total mass $M_{\rm tot} \sim 10^6 M_{\odot}$ at high redshift $ z \gtrsim 3$.

In many braneworld models, gravity is largely modified at the electro-weak scale ~ 1TeV. In such models, primordial black holes (PBHs) with lunar mass M ~ 10^{-7}M_sun might have been produced when the temperature of the universe was at ~ 1TeV. If a significant fraction of the dark halo of our galaxy consists of these lunar mass PBHs, a huge number of BH binaries will exist in our neighborhood. Third generation detectors such as EURO can detect gravitational waves from these binaries, and can also determine their chirp mass. With a new detector designed to be sensitive at high frequency bands greater than 1 kHz, the existence of extradimensions could be confirmed.

A significant fraction of non-baryonic or baryonic dark matter in galactic halos may consist of MASsive Compact Objects (MASCOs) with mass M=10^{1-4}M_{sun}. Possible candidates for such compact objects include primordial black holes or remnants of primordial (Population III) stars. We propose a method for directly detecting MASCOs in extragalactic halos, using the VLBI techniques with extremely high resolution. If a galactic halo comprising a large number of MASCOs produces multiple images of a background radio-loud QSO by gravitational lensing, then a high-resolution radio map of each macro-lensed image should reveal microlensing effects by MASCOs. To assess their observational feasibility, we simulate microlensing of the radio-loud, four-image lensed QSO, B1422+231, assuming angular resolution of ~0.01 mas. MASCOs are represented by point masses. For comparison, we also simulate microlensing of B1422+231 by singular isothermal spheres. We find that the surface brightness of the macro-lensed images shows distinct spatial patterns on the scale of the Einstein radius of the perturbers. In the case of point-mass perturbers, many tiny dark spots also appear in the macro-lensed images associated with a decrease in the surface brightness toward the fringe of the original QSO image, whereas no such spots are available in the SIS models. Based on the size, position and magnified or demagnified patterns of images, we shall be able to determine the mass and density profile of a MASCO as well as its spatial distribution and abundance in a galactic halo.

We reexamine constraints on the spatial size of closed toroidal models with cold dark matter and the cosmological constant from cosmic microwave background. We carry out Bayesian analyses using the Cosmic Background Explorer (COBE) data properly taking into account the statistically anisotropic correlation, i.e., off-diagonal elements in the covariance. We find that the COBE constraint becomes more stringent in comparison with that using only the angular power spectrum, if the likelihood is marginalized over the orientation of the observer. For some limited choices of orientations, the fit to the COBE data is considerably better than that of the infinite counterpart. The best-fit matter normalization is increased because of large-angle suppression in the power and the global anisotropy of the temperature fluctuations. We also study several deformed closed toroidal models in which the fundamental cell is described by a rectangular box. In contrast to the cubic models, the large-angle power can be enhanced in comparison with the infinite counterparts if the cell is sufficiently squashed in a certain direction. It turns out that constraints on some slightly deformed models are less stringent. We comment on how these results affect our understanding of the global topology of our universe.

We revisit the cloud-in-cloud problem for non-Gaussian density fluctuations. We show that the extended Press-Schechter (EPS) formalism for non-Gaussian fluctuations has a flaw in describing mass functions regardless of type of filtering. As an example, we consider non-Gaussian models in which density fluctuations at a point obeys a \chi^2 distribution with \nu degrees of freedom. We find that mass functions predicted by using an integral formula proposed by Jedamzik, and Yano, Nagashima and Gouda, properly taking into account correlation between objects at different scales, deviate from those predicted by using the EPS formalism, especially for strongly non-Gaussian fluctuations. Our results for the mass function at large mass scales are consistent with those by Avelino and Viana obtained from numerical simulations.

We revisit the observational constraints on compact(closed) hyperbolic(CH) models from cosmic microwave background(CMB). We carry out Bayesian analyses for CH models with volume comparable to the cube of the present curvature radius using the COBE-DMR data and show that a slight suppression in the large-angle temperature correlations owing to the non-trivial topology explains rather naturally the observed anomalously low quadrupole which is incompatible with the prediction of the standard infinite Friedmann-Robertson-Walker models. While most of positions and orientations are ruled out, the likelihoods of CH models are found to be much better than those of infinite counterparts for some specific positions and orientations of the observer, leading to less stringent constraints on the volume of the manifolds. Even if the spatial geometry is nearly flat as $\Omega_{tot}=0.9-0.95$, suppression of the angular power on large angular scales is still prominent for CH models with volume much less than the cube of the present curvature radius if the cosmological constant is dominant at present.

In this paper, the cosmic microwave background (CMB) anisotropy in a multiply-connected compact flat 3-torus model with the cosmological constant is investigated. Using the COBE-DMR 4-year data, a full Bayesian analysis revealed that the constraint on the topology of the flat 3-torus model with low-matter-density is less stringent. As in compact hyperbolic models, the large-angle temperature fluctuations can be produced as the gravitational potential decays at the $\Lambda$-dominant epoch well after the last scattering. The maximum allowed number $N$ of images of the cell (fundamental domain) within the observable region at present is approximately 49 for $\Omega_m=0.1$ and $\Omega_\Lambda=0.9$ whereas $N\sim8$ for $\Omega_m=1.0$ and $\Omega_\Lambda=0$.

We have investigated the cosmic microwave background (CMB) anisotropy in closed multiply connected universes (flat and hyperbolic) with low matter density. We show that the COBE constraints on these low matter density models with non-trivial topology are less stringent since a large amount of CMB anisotropy on large angular scales can be produced due to the decay of the gravitational potential at late time.

The effect of a non-trivial topology on the temperature correlations on the cosmic microwave background (CMB) in a small compact hyperbolic universe with volume comparable to the cube of the curvature radius is investigated. Because the bulk of large-angular CMB fluctuations is produced at the late epoch in low $\Omega_0$ models, the effect of a long wavelength cut-off due to the periodic structure does not lead to the significant suppression of large-angular power as in compact flat models. The angular power spectra are consistent with the COBE data for $\Omega_0 \ge 0.1$.

In this paper, Gaussianity of eigenmodes and non-Gaussianity in the Cosmic Microwave Background (CMB) temperature fluctuations in two smallest compact hyperbolic (CH) models are investigated. First, it is numerically found that the expansion coefficients of low-lying eigenmodes on the two CH manifolds behave as if they are Gaussian random numbers at almost all the places. Next, non-Gaussianity of the temperature fluctuations in the (l,m) space in these models is studied. Assuming that the initial fluctuations are Gaussian, the real expansion coefficients b_{l m} of the temperature fluctuations in the sky are found to be distinctively non-Gaussian. In particular, the cosmic variances are found to be much larger than that for Gaussian models. On the other hand, the anisotropic structure is vastly erased if one averages the fluctuations at a number of different observing points because of the Gaussian pseudo-randomness of the eigenmodes. Thus the dominant contribution to the two-point correlation functions comes from the isotropic terms described by the angular power spectra C_l. Finally, topological quantities: the total length and the genus of isotemperature contours are investigated. The variances of total length and genus at high and low threshold levels are found to be considerably larger than that of Gaussian models while the means almost agree with them.

Measurements of cosmic microwave background (CMB) anisotropy are ideal experiments for discovering the non-trivial global topology of the universe. To evaluate the CMB anisotropy in multiply-connected compact cosmological models, one needs to compute the eigenmodes of the Laplace-Beltrami operator. Using the direct boundary element method, we numerically obtain the low-lying eigenmodes on a compact hyperbolic 3-space called the Thurston manifold which is the second smallest in the known compact hyperbolic 3-manifolds. The computed eigenmodes are expanded in terms of eigenmodes on the unit three-dimensional pseudosphere. We numerically find that the expansion coefficients behave as Gaussian pseudo-random numbers for low-lying eigenmodes. The observed gaussianity in the CMB fluctuations can partially be attributed to the Gaussian pseudo-randomness of the expansion coefficients assuming that the Gaussian pseudo-randomness is the universal property of the compact hyperbolic spaces.

Measurements of CMB anisotropy are ideal experiments for discovering the non-trivial global topology of the universe. To evaluate the CMB anisotropy in multiply-connected compact cosmological models, one needs to compute eigenmodes of the Laplace-Beltrami operator. We numerically obtain the eigenmodes on a compact 3-hyperbolic space cataloged as $m003(-2,3)$ in SnapPea (computer program by Jeff Weeks) using the direct boundary element method, which enables one to simulate the CMB in multiply-connected compact models with high precision. The angular power spectra $C_{l}$'s ($2 \leq l\leq18$) are calculated using computed eigenmodes for $5.4\leqk<10$ and Gaussian random approximation for the expansion coefficients for $10\leqk<20$. Assuming that the initial power spectrum is the Harrison-Zeldovich spectrum, the computed $C_l$'s are consistent with the COBE data for $0.1 \leq\Omega_o \leq 0.6$. In low $\Omega_o$ models, the large-angular fluctuations can be produced at periods after the last scattering as the curvature perturbations decay in the curvature dominant era.