A consensus view on the formation of planetesimals is now exposed to a threat, since recent numerical studies on the mechanical properties of dust aggregates tend to dispute the conceptual picture that submicrometer-sized grains conglomerate into planetesimals in protoplanetary disks. With the advent of precise laboratory experiments and extensive computer simulations on the interaction between elastic spheres comprising dust aggregates, we revisit a model for the tensile strength of dust aggregates consisting of small elastic grains. In the framework of contact mechanics and fracture mechanics, we examine outcomes of computer simulations and laboratory experiments on the tensile strength of dust aggregates. We provide a novel analytical formula that explicitly incorporates the volume effect on the tensile strength, namely, the dependence of tensile strength on the volume of dust aggregates. We find that our model for the tensile strength of dust aggregates well reproduces results of computer simulations and laboratory experiments, if appropriate values are adopted for the elastic parameters used in the model. Moreover, the model with dust aggregates of submicrometer-sized grains is in good harmony with the tensile strength of cometary dust and meteoroids derived from astronomical observations. Therefore, we reaffirm the commonly believed idea that the formation of planetesimals begins with conglomeration of submicrometer-sized grains condensed in protoplanetary disks.

We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

We report results from searches for new physics with low-energy electronic recoil data recorded with the XENON1T detector, with an exposure of 0.65 t-y and an unprecedentedly low background rate of $76 \pm 2_{stat}$ events/(t y keV) between 1-30 keV. An excess over known backgrounds is observed below 7 keV, rising towards lower energies and prominent between 2-3 keV. The solar axion model has a 3.5$\sigma$ significance, and a three-dimensional 90% confidence surface is reported for axion couplings to electrons, photons, and nucleons. This surface is inscribed in the cuboid defined by $g_{ae} < 3.7 \times 10^{-12}$, $g_{ae}g_{an}^{eff} < 4.6 \times 10^{-18}$, and $g_{ae}g_{a\gamma} < 7.6\times10^{-22}~{GeV}^{-1}$, and excludes either $g_{ae}=0$ or $g_{ae}g_{a\gamma}=g_{ae}g_{an}^{eff}=0$. The neutrino magnetic moment signal is similarly favored over background at 3.2$\sigma$ and a confidence interval of $\mu_{\nu} \in (1.4,~2.9)\times10^{-11}\mu_B$ (90% C.L.) is reported. Both results are in strong tension with stellar constraints. The excess can also be explained by $\beta$ decays of tritium at 3.2$\sigma$ significance with a corresponding tritium concentration in xenon of $(6.2 \pm 2.0) \times 10^{-25}$ mol/mol. Such a trace amount can be neither confirmed nor excluded with current knowledge of production and reduction mechanisms. The significances of the solar axion and neutrino magnetic moment hypotheses are decreased to 2.1$\sigma$ and 0.9$\sigma$, respectively, if an unconstrained tritium component is included in the fitting. With respect to bosonic dark matter, the excess favors a monoenergetic peak at ($2.3 \pm 0.2$) keV (68% C.L.) with a 3.0$\sigma$ global (4.0$\sigma$ local) significance over background. This analysis sets the most restrictive direct constraints to date on pseudoscalar and vector bosonic dark matter for most masses between 1 and 210 keV/c$^2$.

We have searched for exotic neutrino-electron interactions that could be produced by a neutrino millicharge, by a neutrino magnetic moment, or by dark photons using solar neutrinos in the XMASS-I liquid xenon detector. We observed no significant signals in 711 days of data. We obtain an upper limit for neutrino millicharge of 5.4$\times$10$^{-12} e$ at 90\% confidence level assuming all three species of neutrino have common millicharge. We also set flavor dependent limits assuming the respective neutrino flavor is the only one carrying a millicharge, $7.3 \times 10^{-12} e$ for $\nu_e$, $1.1 \times 10^{-11} e$ for $\nu_{\mu}$, and $1.1 \times 10^{-11} e$ for $\nu_{\tau}$. These limits are the most stringent yet obtained from direct measurements. We also obtain an upper limit for the neutrino magnetic moment of 1.8$\times$10$^{-10}$ Bohr magnetons. In addition, we obtain upper limits for the coupling constant of dark photons in the $U(1)_{B-L}$ model of 1.3$\times$10$^{-6}$ if the dark photon mass is 1$\times 10^{-3}$ MeV$/c^{2}$, and 8.8$\times$10$^{-5}$ if it is 10 MeV$/c^{2}$.

We construct a model of H$\alpha$ emitters (HAEs) based on a semi-analytic galaxy formation model, the New Numerical Galaxy Catalog ($\nu^2$GC). In this paper, we report our estimate for the field variance of the HAE distribution. By calculating the H$\alpha$ luminosity from the star-formation rate of galaxies, our model well reproduces the observed H$\alpha$ luminosity function (LF) at $z=0.4$. The large volume of the $\nu^2$GC makes it possible to examine the spatial distribution of HAEs over a region of (411.8 Mpc)$^3$ in the comoving scale. The surface number density of $z=0.4$ HAEs with $L_{\rm H\alpha} \geq 10^{40}$ erg s$^{-1}$ is 308.9 deg$^{-2}$. We have confirmed that the HAE is a useful tracer for the large-scale structure of the Universe because of their significant overdensity ($>$ 5$\sigma$) at clusters and the filamentary structures. The H$\alpha$ LFs within a survey area of $\sim$2 deg$^2$ (typical for previous observational studies) show a significant field variance up to $\sim$1 dex. Based on our model, one can estimate the variance on the H$\alpha$ LFs within given survey areas.

Xenon dual-phase time projection chambers designed to search for Weakly Interacting Massive Particles have so far shown a relative energy resolution which degrades with energy above $\sim$200 keV. This has limited their sensitivity in the search for rare events like the neutrinoless double-beta decay of $^{136}$Xe at its $Q$-value, $Q_{\beta\beta}\simeq$ 2.46 MeV. For the XENON1T dual-phase time projection chamber, we demonstrate that the relative energy resolution at 1 $\sigma/\mu$ is as low as (0.79$\pm$0.02) % in its one-ton fiducial mass, and for single-site interactions at $Q_{\beta\beta}$. We achieve this by using a signal correction method to rectify the saturation effects of the signal readout system. The very good energy resolution from keV to MeV energies demonstrated in XENON1T opens up new windows for the xenon dual-phase dark matter detectors to simultaneously search for other rare events.

We present spatial variations of the CO J=2-1/1-0 line ratio in M83 using Total Power array data from ALMA. While the intensities of these two lines correlate tightly, the ratio varies over the disk, with a disk average ratio of 0.69, and shows the galactic center and a two-arm spiral pattern. It is high (>0.7) in regions of high molecular gas surface density, but ranges from low to high ratios in regions of low surface density. The ratio correlates well with the spatial distributions and intensities of FUV and IR emissions, with FUV being the best correlated. It also correlates better with the ratio of IR intensities (70/350mic), a proxy for dust temperature, than with the IR intensities. Taken together, these results suggest either a direct or indirect link between the dust heating by the interstellar radiation field and the condition of GMCs, even though no efficient mechanism is known for a thermal coupling of dust and bulk gas in GMCs.

We present our ALMA Band 8 observations of a damped Ly${\alpha}$ absorption (DLA) system at $z$=3.150 observed in the spectrum of the quasar Q2233+131 at $z$=3.295. The optical counterpart of this DLA has been identified and it shows a double-peaked Ly${\alpha}$ emission line. Since one possible origin of DLAs at high redshift is an outflowing gas from star-forming galaxies, DLA2233+131 provides a good laboratory to investigate the nature of high-$z$ DLAs. Motivated by this, we have carried out ALMA band 8 observations to study the [C II] line in this system. However, we do not detect any significant emission line in the observed pass bands. Instead, we have serendipitously found three submm continuum sources in the observed sky area. One appears to be the quasar Q2233+131 itself while the other two sources are newly identified submm galaxies (SMGs), called SMG1 and SMG2 in this paper. They are located at a separation of 400.7 and 800.1 from Q2233+131, respectively. Their 646 ${\mu}$m fluxes are 6.35 mJy and 6.43 mJy, respectively, being higher than that of Q2233+131, 3.62 mJy. Since these two SMGs are not detected in the optical images obtained with the Hubble Space Telescope and the Subaru Telescope, they have a very red spectral energy distribution. It is, therefore,suggested that they are high-redshift galaxies or very dusty galaxies at intermediate redshift although we cannot rule out the possibility that they are optically very faint SMG analogs at low redshift. Follow up observations will be necessary to explore the nature of this interesting region.

Direct dark matter detection experiments based on a liquid xenon target are leading the search for dark matter particles with masses above $\sim$ 5 GeV/c$^2$, but have limited sensitivity to lighter masses because of the small momentum transfer in dark matter-nucleus elastic scattering. However, there is an irreducible contribution from inelastic processes accompanying the elastic scattering, which leads to the excitation and ionization of the recoiling atom (the Migdal effect) or the emission of a Bremsstrahlung photon. In this letter, we report on a probe of low-mass dark matter with masses down to about 85 MeV/c$^2$ by looking for electronic recoils induced by the Migdal effect and Bremsstrahlung, using data from the XENON1T experiment. Besides the approach of detecting both scintillation and ionization signals, we exploit an approach that uses ionization signals only, which allows for a lower detection threshold. This analysis significantly enhances the sensitivity of XENON1T to light dark matter previously beyond its reach.

We report constraints on light dark matter (DM) models using ionization signals in the XENON1T experiment. We mitigate backgrounds with strong event selections, rather than requiring a scintillation signal, leaving an effective exposure of $(22 \pm 3)$ tonne-days. Above $\sim\!0.4\,\mathrm{keV}_\mathrm{ee}$, we observe $<1 \, \text{event}/(\text{tonne} \times \text{day} \times \text{keV}_\text{ee})$, which is more than one thousand times lower than in similar searches with other detectors. Despite observing a higher rate at lower energies, no DM or CEvNS detection may be claimed because we cannot model all of our backgrounds. We thus exclude new regions in the parameter spaces for DM-nucleus scattering for DM masses $m_\chi$ within $3-6\,\mathrm{GeV}/\mathrm{c}^2$, DM-electron scattering for $m_\chi > 30\,\mathrm{MeV}/\mathrm{c}^2$, and absorption of dark photons and axion-like particles for $m_\chi$ within $0.186 - 1 \, \mathrm{keV}/\mathrm{c}^2$.

The Ginzburg-Landau (GL) effective theory is a useful tool to study a superconductivity or superfluidity near the critical temperature, and usually the expansion up to the 4th order in terms of order parameters is sufficient for the description of the second-order phase transition. In this paper, we discuss the GL equation for the neutron $^{3}P_{2}$ superfluidity relevant for interior of neutron stars. We derive the GL expansion up to the 8th order in the condensates and find that this order is necessary for the system to have the unique ground state, unlike the ordinary cases. Starting from the $LS$ potential, which provides the dominant attraction between two neutrons at the high density, we derive the GL equation in the path-integral formalism, where the auxiliary field method and the Nambu-Gor'kov representation are used. We present the detailed description for the trace calculation necessary in the derivation of the GL equation. As numerical results, we show the phase diagram of the neutron $^{3}P_{2}$ superfluidity on the plane spanned by the temperature and magnetic field, and find that the 8th order terms lead to a first-order phase transition, whose existence was predicted in the Bogoliubov-de Gennes equation but has not been found thus far within the framework of the GL expansion up to the 6th order.The first-order phase transition will affect the interior structures inside the neutron stars.

We report the Subaru Hyper Suprime-Cam (HSC) discovery of two Ly$\alpha$ blobs (LABs), dubbed z70-1 and z49-1 at $z=6.965$ and $z=4.888$ respectively, that are Ly$\alpha$ emitters with a bright ($\log L_{\rm Ly\alpha}/{\rm [erg\ s^{-1}]}>43.4$) and spatially-extended Ly$\alpha$ emission, and present the photometric and spectroscopic properties of a total of seven LABs; the two new LABs and five previously-known LABs at $z=5.7-6.6$. The z70-1 LAB shows the extended Ly$\alpha$ emission with a scale length of $1.4\pm 0.2$ kpc, about three times larger than the UV continuum emission, making z70-1 the most distant LAB identified to date. All of the 7 LABs, except z49-1, exhibit no AGN signatures such as X-ray emission, {\sc Nv}$\lambda$1240 emission, or Ly$\alpha$ line broadening, while z49-1 has a strong {\sc Civ}$\lambda$1548 emission line indicating an AGN on the basis of the UV-line ratio diagnostics. We carefully model the point-spread functions of the HSC images, and conduct two-component exponential profile fitting to the extended Ly$\alpha$ emission of the LABs. The Ly$\alpha$ scale lengths of the core (star-forming region) and the halo components are $r_{\rm c}=0.6-1.2$ kpc and $r_{\rm h}=2.0-13.8$ kpc, respectively. The average $r_{\rm h}$ of the LABs falls on the extrapolation of the $r_{\rm h}$-Ly$\alpha$ luminosity relation of the Ly$\alpha$ halos around VLT/MUSE star-forming galaxies at the similar redshifts, suggesting that typical LABs at $z\gtrsim5$ are not special objects, but star-forming galaxies at the bright end.

We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

The JAXA/ISAS spacecraft DESTINY$^+$ will be launched to the active asteroid (3200) Phaethon in 2022. Among the proposed core payload is the DESTINY+ Dust Analyzer (DDA) which is an upgrade of the Cosmic Dust Analyzer flown on the Cassini spacecraft to Saturn (Srama et al. 2011). We use two up-to-date computer models, the ESA Interplanetary Meteoroid Engineering Model (IMEM, Dikarev et al. 2005), and the interstellar dust module of the Interplanetary Meteoroid environment for EXploration model (IMEX; Sterken2013 et al., Strub et al. 2019) to study the detection conditions and fluences of interplanetary and interstellar dust with DDA. Our results show that a statistically significant number of interplanetary and interstellar dust particles will be detectable with DDA during the 4-years interplanetary cruise of DESTINY+. The particle impact direction and speed can be used to descriminate between interstellar and interplanetary particles and likely also to distinguish between cometary and asteroidal particles.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

A search for Weakly Interacting Massive Particles (WIMPs) was conducted with the single-phase liquid-xenon detector XMASS through inelastic scattering in which $^{129}$Xe nuclei were excited, using an exposure ($\rm 327\; kg \times 800.0 \; days$) 48 times larger than that of our previous study. The inelastic excitation sensitivity was improved by detailed evaluation of background, event classification based on scintillation timing that distinguished $\gamma$-rays and $\beta$-rays, and simultaneous fitting of the energy spectra of $\gamma$-like and $\beta$-like samples. No evidence of a WIMP signal was found. Thus, we set the upper limits of the inelastic channel cross section at 90\% confidence level, for example, $4.1\times 10^{-39} \;{\rm cm^2}$ for a $200\; {\rm GeV}/c^2$ WIMP. This result provides the most stringent limits on the SD WIMP-neutron interaction and is better by a factor of 7.7 at $200\;{\rm GeV}/c^2$ than the existing experimental limit.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

We examined radial variations in molecular-gas based star formation efficiency (SFE), which is defined as star formation rate per unit molecular gas mass, for 80 galaxies selected from the CO Multi-line Imaging of Nearby Galaxies project (Sorai et al. 2019). The radial variations in SFE for individual galaxies are typically a factor of 2 -- 3, which suggests that SFE is nearly constant along galactocentric radius. We found the averaged SFE in 80 galaxies of $(1.69 \pm 1.1) \times 10^{-9}$ yr$^{-1}$, which is consistent with Leroy et al. 2008 if we consider the contribution of helium to the molecular gas mass evaluation and the difference in the assumed initial mass function between two studies. We compared SFE among different morphological (i.e., SA, SAB, and SB) types, and found that SFE within the inner radii ($r/r_{25} < 0.3$, where $r_{25}$ is $B$-band isophotal radius at 25 mag arcsec$^{-2}$) of SB galaxies is slightly higher than that of SA and SAB galaxies. This trend can be partly explained by the dependence of SFE on global stellar mass, which probably relates to the CO-to-H$_2$ conversion factor through the metallicity. For two representative SB galaxies in our sample, NGC 3367 and NGC 7479, the ellipse of $r/r_{25}$ = 0.3 seems to cover not only the central region but also the inner part of the disk, mainly the bar. These two galaxies show higher SFE in the bar than in spiral arms. However, we found an opposite trend in NGC 4303; SFE is lower in the bar than in spiral arms, which is consistent with earlier studies (e.g., Momose et al. 2010). These results suggest diversity of star formation activities in the bar.

XMASS-I is a single-phase liquid xenon detector whose purpose is direct detection of dark matter. To achieve the low background requirements necessary in the detector, a new model of photomultiplier tubes (PMTs), R10789, with a hexagonal window was developed based on the R8778 PMT used in the XMASS prototype detector. We screened the numerous component materials for their radioactivity. During development, the largest contributions to the reduction of radioactivity came from the stem and the dynode support. The glass stem was exchanged to the Kovar alloy one and the ceramic support were changed to the quartz one. R10789 is the first model of Hamamatsu Photonics K. K. that adopted these materials for low background purposes and provided a groundbreaking step for further reductions of radioactivity in PMTs. Measurements with germanium detectors showed 1.2$\pm$0.3 mBq/PMT of $^{226}$Ra, less than 0.78 mBq/PMT of $^{228}$Ra, 9.1$\pm$2.2 mBq/PMT of $^{40}$K, and 2.8$\pm$0.2 mBq/PMT of $^{60}$Co. In this paper, the radioactive details of the developed R10789 are described together with our screening methods and the components of the PMT.

A search for dark matter using an underground single-phase liquid xenon detector was conducted at the Kamioka Observatory in Japan, particularly for Weakly Interacting Massive Particles (WIMPs). We have used 705.9 live days of data in a fiducial volume containing 97 kg of liquid xenon at the center of the detector. The event rate in the fiducial volume after the data reduction was ${\rm (4.2 \pm 0.2) \times 10^{-3} \, day^{-1}kg^{-1} keV_{ee}^{-1}}$ at ${\rm 5 \, keV_{ee}}$, with a signal efficiency of ${\rm 20\%}$. All the remaining events are consistent with our background evaluation, mostly of the "mis-reconstructed events" originated from $^{210}$Pb in the copper plates lining the detector's inner surface. The obtained upper limit on a spin-independent WIMP-nucleon cross section was ${\rm 2.2 \times 10^{-44} \, cm^{2}}$ for a WIMP mass of ${\rm 60 \, GeV/c^{2}}$ at the $90\%$ confidence level, which was the most stringent limit among results from single-phase liquid xenon detectors.

A search for dark matter (DM) with mass in the sub-GeV region (0.32-1 GeV) was conducted by looking for an annual modulation signal in XMASS, a single-phase liquid xenon detector. Inelastic nuclear scattering accompanied by bremsstrahlung emission was used to search down to an electron equivalent energy of 1 keV. The data used had a live time of 2.8 years (3.5 years in calendar time), resulting in a total exposure of 2.38 ton-years. No significant modulation signal was observed and 90% confidence level upper limits of $1.6 \times 10^{-33}$ cm$^2$ at 0.5 GeV was set for the DM-nucleon cross section. This is the first experimental result of a search for DM mediated by the bremsstrahlung effect. In addition, a search for DM with mass in the multi-GeV region (4-20 GeV) was conducted with a lower energy threshold than previous analysis of XMASS. Elastic nuclear scattering was used to search down to a nuclear recoil equivalent energy of 2.3 keV, and upper limits of 2.9 $\times$10$^{-42}$ cm$^2$ at 8 GeV was obtained.

Hidden photons and axion-like particles are candidates for cold dark matter if they were produced non-thermally in the early universe. We conducted a search for both of these bosons using 800 live-days of data from the XMASS detector with 327 kg of liquid xenon in the fiducial volume. No significant signal was observed, and thus we set constraints on the $\alpha' / \alpha$ parameter related to kinetic mixing of hidden photons and the coupling constant $g_{Ae}$ of axion-like particles in the mass range from 40 to 120 keV/$c^2$, resulting in $\alpha' / \alpha < 6 \times 10^{-26}$ and $g_{Ae} < 4 \times 10^{-13}$. These limits are the most stringent derived from both direct and indirect searches to date.

We present the latest results of a semi-analytic model of galaxy formation, "New Numerical Galaxy Catalogue", which is combined with large cosmological N-body simulations. This model can reproduce statistical properties of galaxies at z < 6.0. We focus on the properties of active galactic nuclei (AGNs) and supermassive black holes, especially on the accretion timescale onto black holes. We find that the number density of AGNs at z < 1.5 and at hard X-ray luminosity 10^{ 44 }< erg/s is underestimated compared with recent observational estimates when we assume the exponentially decreasing accretion rate and the accretion timescale which is proportional to the dynamical time of the host halo or the bulge, as is often assumed in semi-analytic models. We show that to solve this discrepancy, the accretion timescale of such less luminous AGNs instead should be a function of the black hole mass and the accreted gas mass. This timescale can be obtained from a phenomenological modelling of the gas angular momentum loss in the circumnuclear torus and/or the accretion disc. Such models predict a longer accretion timescale for less luminous AGNs at z < 1.0 than bright QSOs whose accretion timescale would be 10^{ 7-8 } yr. With this newly introduced accretion timescale, our model can explain the observed luminosity functions of AGNs at z < 6.0.

The observed value $\Lambda_{\rm obs}$ of the cosmological constant $\Lambda$ is extremely smaller than theoretical expectations, and the anthropic argument has been proposed as a solution to this problem because galaxies do not form when $\Lambda \gg \Lambda_{\rm obs}$. However, the contemporary galaxy formation theory predicts that stars form even with a high value of $\Lambda / \Lambda_{\rm obs} \sim$ 50, which makes the anthropic argument less persuasive. Here we calculate the probability distribution of $\Lambda$ using a model of cosmological galaxy formation, considering extinction of observers caused by radiation from nearby supernovae. The life survival probability decreases in a large $\Lambda$ universe because of higher stellar density. Using a reasonable rate of lethal supernovae, we find that the mean expectation value of $\Lambda$ can be close to $\Lambda_{\rm obs}$, and hence this effect may be essential to understand the small but nonzero value of $\Lambda$. It is predicted that we are located on the edge of habitable regions about stellar density in the Galaxy, which may be tested by future exoplanet studies.

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

An annual modulation signal due to the Earth orbiting around the Sun would be one of the strongest indications of the direct detection of dark matter. In 2016, we reported a search for dark matter by looking for this annual modulation with our single-phase liquid xenon XMASS-I detector. That analysis resulted in a slightly negative modulation amplitude at low energy. In this work, we included more than one year of additional data, which more than doubles the exposure to 800 live days with the same 832 kg target mass. When we assume weakly interacting massive particle (WIMP) dark matter elastically scattering on the xenon target, the exclusion upper limit for the WIMP-nucleon cross section was improved by a factor of 2 to 1.9$\times$10$^{-41}$cm$^2$ at 8 GeV/c$^2$ at 90\% confidence level with our newly implemented data selection through a likelihood method. For the model-independent case, without assuming any specific dark matter model, we obtained more consistency with the null hypothesis than before with a $p$-value of 0.11 in the 1$-$20 keV energy region. This search probed this region with an exposure that was larger than that of DAMA/LIBRA. We also did not find any significant amplitude in the data for periodicity with periods between 50 and 600 days in the energy region between 1 to 6 keV.

Recent radio observations towards nearby galaxies started to map the whole disk and to identify giant molecular clouds (GMCs) even in the regions between galactic spiral structures. Observed variations of GMC mass functions in different galactic environment indicates that massive GMCs preferentially reside along galactic spiral structures whereas inter-arm regions have many small GMCs. Based on the phase transition dynamics from magnetized warm neutral medium to molecular clouds, Kobayashi et al. 2017 proposes a semi-analytical evolutionary description for GMC mass functions including cloud-cloud collision (CCC) process. Their results show that CCC is less dominant in shaping the mass function of GMCs compared with the accretion of dense HI gas driven by the propagation of supersonic shock waves. However, their formulation does not take into account the possible enhancement of star formation by CCC. Radio observations within the Milky Way indicate the importance of CCC for the formation of star clusters and massive stars. In this article, we reformulate the time evolution equation largely modified from Kobayashi et al. 2017 so that we additionally compute star formation subsequently taking place in CCC clouds. Our results suggest that, although CCC events between smaller clouds outnumber the ones between massive GMCs, CCC-driven star formation is mostly driven by massive GMCs > 10^5.5 Msun (where Msun is the solar mass). The resultant cumulative CCC-driven star formation may amount to a few 10 percent of the total star formation in the Milky Way and nearby galaxies.

The 21-cm signal in the vicinity of the first stars is expected to reflect properties of the first stars. In this study we pay special attention to tracing the time evolution of the ionizing photons' escape fraction, which affects the distribution of neutral hydrogen, by performing radiation hydrodynamics (RHD) simulations resolving dense gas in a halo. We find that the radial profile of 21-cm differential brightness temperature is quite sensitive to the stellar and halo masses, which reflects the time evolution of the escape fraction. In the case of a less massive star, ionizing photons hardly escape from its host halo due to the absorption by dense halo gas, thus an deep 21-cm absorption feature at just outside the halo lasts a long time. Whereas photons from a massive star well working to heat the ambient intergalactic medium turn out to cause a spatially extended 21-cm emission signature. Although individual signals are found to be undetectable with the Square Kilometre Array, our analysis using the results from the RHD simulations indicates that the properties of the first stars are imprinted on the 21-cm global signal: its amplitude depends not only on the cosmic star formation rate density, but also on the typical mass of the first stars due to the stellar-mass-dependent heating rate. Thus, we suggest that the initial mass function of the first stars is an essential factor in understanding the global signal.

We present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

In theories with the large extra dimensions beyond the standard 4-dimensional spacetime, axions could propagate in such extra dimensions, and acquire Kaluza-Klein (KK) excitations. These KK axions are produced in the Sun and could solve unexplained heating of the solar corona. While most of the solar KK axions escape from the solar system, a small fraction is gravitationally trapped in orbits around the Sun. They would decay into two photons inside a terrestrial detector. The event rate is expected to modulate annually depending on the distance from the Sun. We have searched for the annual modulation signature using $832\times 359$ kg$\cdot$days of XMASS-I data. No significant event rate modulation is found, and hence we set the first experimental constraint on the KK axion-photon coupling of $4.8 \times 10^{-12}\, \mathrm{GeV}^{-1}$ at 90% confidence level for a KK axion number density of $\bar{n}_\mathrm{a} = 4.07 \times 10^{13}\, \mathrm{m}^{-3}$, the total number of extra dimensions $n = 2$, and the number of extra dimensions $\delta = 2$ that axions can propagate in.

We present the Ly$\alpha$ luminosity functions (LFs) at $z=$5.7 and 6.6 derived from a new large sample of 1,266 Ly$\alpha$ emitters (LAEs) identified in total areas of 14 and 21 deg$^2$, respectively, based on the early narrowband data of the Subaru/Hyper Suprime-Cam (HSC) survey. Together with careful Monte-Carlo simulations that account for the incompleteness of the LAE selection and the flux estimate systematics in the narrowband imaging, we have determined the Ly$\alpha$ LFs with the unprecedentedly small statistical and systematic uncertainties in a wide Ly$\alpha$ luminosity range of $10^{42.8-43.8}$ erg s$^{-1}$. We obtain the best-fit Schechter parameters of $L^{*}_{{\rm Lya}}=1.6^{+2.2}_{-0.6} (1.7^{+0.3}_{-0.7}) \times10^{43}$ erg s$^{-1}$, $\phi^{*}_{{\rm Lya}}=0.85^{+1.87}_{-0.77}\ (0.47^{+1.44}_{-0.44})\times10^{-4}$ Mpc$^{-3}$, and $\alpha=-2.6^{+0.6}_{-0.4}\ (-2.5^{+0.5}_{-0.5})$ at $z=5.7$ ($6.6$). We confirm that our best-estimate Ly$\alpha$ LFs are consistent with the majority of the previous studies, but find that our Ly$\alpha$ LFs do not agree with the high number densities of LAEs recently claimed by Matthee/Santos et al.'s studies that may overcorrect the incompleteness and the flux systematics. Our Ly$\alpha$ LFs at $z=5.7$ and $6.6$ show an indication that the faint-end slope is very steep ($\alpha \simeq -2.5$), although it is also possible that the bright-end LF results are enhanced by systematic effects such as the contribution from AGNs, blended merging galaxies, and/or large ionized bubbles around bright LAEs. Comparing our Ly$\alpha$ LF measurements with four independent reionization models, we estimate the neutral hydrogen fraction of the IGM to be $x_{\rm HI}=0.3\pm0.2$ at $z=$6.6 that is consistent with the small Thomson scattering optical depth obtained by Planck 2016.

We present an unprecedentedly large catalog consisting of 2,230 > L* Lya emitters (LAEs) at z=5.7 and 6.6 on the 13.8 and 21.2 deg2 sky, respectively, that are identified by the SILVERRUSH program with the first narrowband imaging data of the Hyper Suprime-Cam (HSC) survey. We confirm that the LAE catalog is reliable on the basis of 96 LAEs whose spectroscopic redshifts are already determined by this program and the previous studies. This catalogue is also available on-line. Based on this catalogue, we derive the rest-frame Lya equivalent-width distributions of LAEs at z~5.7-6.6 that are reasonably explained by the exponential profiles with the scale lengths of ~120-170A, showing no significant evolution from z~5.7 to z~6.6. We find that 275 LAEs with a large equivalent width (LEW) of >240A are candidates of young-metal poor galaxies and AGNs. We also find that the fraction of LEW LAEs to all ones is 4% and 21% at z~5.7 and z~6.6, respectively. Our LAE catalog includes 11 Lya blobs (LABs) that are LAEs with spatially extended Lya emission whose profile is clearly distinguished from those of stellar objects at the >~ 3sigma level. The number density of the LABs at z=6-7 is ~ 10^-7-10^-6 Mpc^-3, being ~ 10-100 times lower than those claimed for LABs at z~ 2-3, suggestive of disappearing LABs at z>~6, albeit with the different selection methods and criteria for the low and high-z LABs.

We present Lya and UV-nebular emission line properties of bright Lya emitters (LAEs) at z=6-7 with a luminosity of log L_Lya/[erg s-1] = 43-44 identified in the 21-deg2 area of the SILVERRUSH early sample developed with the Subaru Hyper Suprime-Cam (HSC) survey data. Our optical spectroscopy newly confirm 21 bright LAEs with clear Lya emission, and contribute to make a spectroscopic sample of 96 LAEs at z=6-7 in SILVERRUSH. From the spectroscopic sample, we select 7 remarkable LAEs as bright as Himiko and CR7 objects, and perform deep Keck/MOSFIRE and Subaru/nuMOIRCS near-infrared spectroscopy reaching the 3sigma-flux limit of ~ 2x10^{-18} erg s-1 for the UV-nebular emission lines of He II1640, C IV1548,1550, and O III]1661,1666. Except for one tentative detection of C IV, we find no strong UV-nebular lines down to the flux limit, placing the upper limits of the rest-frame equivalent widths (EW_0) of ~2-4 A for He II, C IV, and O III] lines. Here we also investigate the VLT/X-SHOOTER spectrum of CR7 whose 6 sigma detection of He II is claimed by Sobral et al. Although two individuals and the ESO-archive service carefully re-analyze the X-SHOOTER data that are used in the study of Sobral et al., no He II signal of CR7 is detected, supportive of weak UV-nebular lines of the bright LAEs even for CR7. Spectral properties of these bright LAEs are thus clearly different from those of faint dropouts at z~7 that have strong UV-nebular lines shown in the various studies. Comparing these bright LAEs and the faint dropouts, we find anti-correlations between the UV-nebular line EW_0 and UV-continuum luminosity, which are similar to those found at z~2-3.

We present the SILVERRUSH program strategy and clustering properties investigated with $\sim 2,000$ Ly$\alpha$ emitters at $z=5.7$ and $6.6$ found in the early data of the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey exploiting the carefully designed narrowband filters. We derive angular correlation functions with the unprecedentedly large samples of LAEs at $z=6-7$ over the large total area of $14-21$ deg$^2$ corresponding to $0.3-0.5$ comoving Gpc$^2$. We obtain the average large-scale bias values of $b_{\rm avg}=4.1\pm 0.2$ ($4.5\pm 0.6$) at $z=5.7$ ($z=6.6$) for $\gtrsim L^*$ LAEs, indicating the weak evolution of LAE clustering from $z=5.7$ to $6.6$. We compare the LAE clustering results with two independent theoretical models that suggest an increase of an LAE clustering signal by the patchy ionized bubbles at the epoch of reionization (EoR), and estimate the neutral hydrogen fraction to be $x_{\rm HI}=0.15^{+0.15}_{-0.15}$ at $z=6.6$. Based on the halo occupation distribution models, we find that the $\gtrsim L^*$ LAEs are hosted by the dark-matter halos with the average mass of $\log (\left < M_{\rm h} \right >/M_\odot) =11.1^{+0.2}_{-0.4}$ ($10.8^{+0.3}_{-0.5}$) at $z=5.7$ ($6.6$) with a Ly$\alpha$ duty cycle of 1 % or less, where the results of $z=6.6$ LAEs may be slightly biased, due to the increase of the clustering signal at the EoR. Our clustering analysis reveals the low-mass nature of $\gtrsim L^*$ LAEs at $z=6-7$, and that these LAEs probably evolve into massive super-$L^*$ galaxies in the present-day universe.

Super-Eddington mass accretion has been suggested as an efficient mechanism to grow supermassive black holes (SMBHs). We investigate the imprint left by the radiative efficiency of the super-Eddington accretion process on the clustering of quasars using a new semi-analytic model of galaxy and quasar formation based on large-volume cosmological $N$-body simulations. Our model includes a simple model for the radiative efficiency of a quasar, which imitates the effect of photon trapping for a high mass accretion rate. We find that the model of radiative efficiency affects the relation between the quasar luminosity and the quasar host halo mass. The quasar host halo mass has only weak dependence on quasar luminosity when there is no upper limit for quasar luminosity. On the other hand, it has significant dependence on quasar luminosity when the quasar luminosity is limited by its Eddington luminosity. In the latter case, the quasar bias also depends on the quasar luminosity, and the quasar bias of bright quasars is in agreement with observations. Our results suggest that the quasar clustering studies can provide a constraint on the accretion disc model.

We detect 20 $z=7.0$ Ly$\alpha$ emitter (LAE) candidates to $L({\rm Ly}\alpha) \geq 2 \times 10^{42}$ erg s$^{-1}$ or $0.3$ $L^*_{z=7}$ and in $6.1\times 10^5$ Mpc$^3$ volume in the Subaru Deep Field and the Subaru/XMM-Newton Deep Survey field by 82 and 37 hours of Subaru Suprime-Cam narrowband NB973 and reddest optical $y$-band imaging. We compare their Ly$\alpha$ and UV luminosity functions (LFs) and densities and Ly$\alpha$ equivalent widths (EWs) to those of $z=5.7$, 6.6 and 7.3 LAEs from previous Suprime-Cam surveys. The Ly$\alpha$ LF (density) rapidly declines by a factor of $\times$1.5 (1.9) in $L({\rm Ly}\alpha)$ at $z=5.7-6.6$ (160 Myr), $\times$1.5 (1.6) at $z=6.6-7.0$ (60 Myr) at the faint end and $\times$2.0 (3.8) at $z=7.0-7.3$ (40 Myr). Also, in addition to the systematic decrease in EW at $z=5.7-6.6$ previously found, 2/3 of the $z=7.0$ LAEs detected in the UV continuum exhibit lower EWs than the $z=6.6$ ones. Moreover, while the UV LF and density do not evolve at $z=5.7-6.6$, they modestly decline at $z=6.6-7.0$, implying galaxy evolution contributing to the decline of the Ly$\alpha$ LF. Comparison of the $z=7.0$ Ly$\alpha$ LF to the one predicted by an LAE evolution model further reveals that galaxy evolution alone cannot explain all the decline of Ly$\alpha$ LF. If we attribute the discrepancy to Ly$\alpha$ attenuation by neutral hydrogen, the intergalactic medium transmission of Ly$\alpha$ photons at $z=7.0$ would be $T_{{\rm Ly}\alpha}^{\rm IGM} \leq 0.6-0.7$. It is lower (higher) than $T_{{\rm Ly}\alpha}^{\rm IGM}$ at $z=6.6$ (7.3) derived by previous studies, suggesting rapid increase in neutral fraction at $z > 6$.

In order to investigate the physical relationship between strong Ly-alpha absorbers (logN_HI>20.0 cm^-2) such as damped Ly-alpha absorption systems (DLAs) and young star-forming galaxies at high redshift, we have conducted narrow-band observations of Ly-alpha emitters (LAEs) in a concentrated region of strong Ly-alpha absorbers at z= 2.3, the J1230+34 field. Using a catalog of Ly-alpha absorbers with logN_HI>20.0 cm^-2 based on the baryon oscillation spectroscopic survey (BOSS), we found 6 fields where 3 or more absorbers are concentrated within a (50 Mpc)^3 cubic box in the comoving scale. Among them, we focus on the J1230+34 field, where 2 DLAs and 2 sub-DLAs present. Our narrow-band imaging observations with Subaru/Suprime-Cam using a custom-made filter, NB400 (lambda_c=4003 A and FWHM=92 A) yield a sample of 149 LAEs in this field. In the large scale (~50 Mpc), we have found no differences between the obtained Ly-alpha luminosity function and those in the blank fields at similar redshifts. We also compare the frequency distribution of the Ly-alpha rest-frame equivalent width (EW_0) in the target field and other fields including both overdensity region and blank field, but find no differences. On the other hand, in the small scale (~10 Mpc), we have found a possible overdensity of LAEs around a DLA with the highest HI column density (N_HI=21.08 cm^-2) in the target field while there are no density excess around the other absorbers with a lower N_HI.

Recent radio observations show that the giant molecular cloud (GMC) mass functions noticeably vary across galactic disks. High-resolution magnetohydrodynamics simulations show that multiple episodes of compression are required for creating a molecular cloud in the magnetized interstellar medium. In this article, we formulate the evolution equation for the GMC mass function to reproduce the observed profiles, for which multiple compression are driven by the network of expanding shells due to HII regions and supernova remnants. We introduce the cloud-cloud collision (CCC) terms in the evolution equation in contrast to the previous work (Inutsuka et al. 2015). The computed time evolution suggests that the GMC mass function slope is governed by the ratio of GMC formation timescale to its dispersal timescale, and that the CCC effect is limited only in the massive-end of the mass function. In addition, we identify a gas resurrection channel that allows the gas dispersed by massive stars to regenerate GMC populations or to accrete onto the pre-existing GMCs. Our results show that almost all of the dispersed gas contribute to the mass growth of pre-existing GMCs in arm regions whereas less than 60 per cent in inter-arm regions. Our results also predict that GMC mass functions have a single power-law exponent in the mass range < 10^5.5 Msun (where Msun represents the solar mass), which is well characterized by GMC self-growth and dispersal timescales. Measurement of the GMC mass function slope provides a powerful method to constrain those GMC timescales and the gas resurrecting factor in various environment across galactic disks.

We present a new cosmological galaxy formation model, $\nu^2$GC, as an updated version of our previous model $\nu$GC. We adopt the so-called "semi-analytic" approach, in which the formation history of dark matter halos is computed by ${\it N}$-body simulations, while the baryon physics such as gas cooling, star formation and supernova feedback are simply modeled by phenomenological equations. Major updates of the model are as follows: (1) the merger trees of dark matter halos are constructed in state-of-the-art ${\it N}$-body simulations, (2) we introduce the formation and evolution process of supermassive black holes and the suppression of gas cooling due to active galactic nucleus (AGN) activity, (3) we include heating of the intergalactic gas by the cosmic UV background, and (4) we tune some free parameters related to the astrophysical processes using a Markov chain Monte Carlo method. Our ${\it N}$-body simulations of dark matter halos have unprecedented box size and mass resolution (the largest simulation contains 550 billion particles in a 1.12 Gpc/h box), enabling the study of much smaller and rarer objects. The model was tuned to fit the luminosity functions of local galaxies and mass function of neutral hydrogen. Local observations, such as the Tully-Fisher relation, size-magnitude relation of spiral galaxies and scaling relation between the bulge mass and black hole mass were well reproduced by the model. Moreover, the model also well reproduced the cosmic star formation history and the redshift evolution of rest-frame ${\it K}$-band luminosity functions. The numerical catalog of the simulated galaxies and AGNs is publicly available on the web.

Recent studies have derived quasar luminosity functions (QLFs) at various redshifts. However, the faint side of the QLF at high redshifts is still too uncertain. An accurate estimate of the survey completeness is essential to derive an accurate QLF for use in studying the luminosity-dependent density evolution of the quasar population. Here we investigate how the luminosity dependence of quasar spectra (the Baldwin effect) and the attenuation model for the inter-galactic medium (IGM) affect the completeness estimates. For this purpose, we revisit the completeness of quasar surveys specifically at $z\sim4-5$, using the COSMOS images observed with Subaru/Suprime-Cam. As the result, we find that the completeness estimates are sensitive to the luminosity dependence of the quasar spectrum and difference in the IGM attenuation models. At $z\sim4$, the number density of quasars when we adopt the latest IGM model and take the luminosity dependence of spectra into account are $(3.49\pm1.62)\times10^{-7}$ Mpc$^{-3}$ mag$^{-1}$ for $-24.09<M_{1450}<-23.09$ and $(5.24\pm2.13)\times10^{-7}$ Mpc$^{-3}$ mag$^{-1}$ for $-23.09<M_{1450}<-22.09$, respectively, which are $\sim$$24$$\%$ lower than that estimated by the conventional method. On the other hand, at $z\sim5$, a $1\sigma$ confidence upper limit of the number density at $-24.5<M_{1450}<-22.5$ in our new estimates is $\sim$$43$$\%$ higher than that estimated previously. The results suggest that the luminosity dependence of the quasar spectrum and the IGM model are important for deriving accurate number density of high-$z$ quasars. Even taking these effects into account, the inferred luminosity-dependent density evolution of quasars is consistent with the AGN down-sizing evolutionary picture.

We report the measurement of the emission time profile of scintillation from gamma-ray induced events in the XMASS-I 832 kg liquid xenon scintillation detector. Decay time constant was derived from a comparison of scintillation photon timing distributions between the observed data and simulated samples in order to take into account optical processes such as absorption and scattering in liquid xenon. Calibration data of radioactive sources, $^{55}$Fe, $^{241}$Am, and $^{57}$Co were used to obtain the decay time constant. Assuming two decay components, $\tau_1$ and $\tau_2$, the decay time constant $\tau_2$ increased from 27.9 ns to 37.0 ns as the gamma-ray energy increased from 5.9 keV to 122 keV. The accuracy of the measurement was better than 1.5 ns at all energy levels. A fast decay component with $\tau_1 \sim 2$ ns was necessary to reproduce data. Energy dependencies of $\tau_2$ and the fraction of the fast decay component were studied as a function of the kinetic energy of electrons induced by gamma-rays. The obtained data almost reproduced previously reported results and extended them to the lower energy region relevant to direct dark matter searches.

We explore the effect of varying the mass of the seed black hole on the resulting black hole mass - bulge mass relation at z ~ 0, using a semi-analytic model of galaxy formation combined with large cosmological N-body simulations. We constrain our model by requiring the observed properties of galaxies at z ~ 0 are reproduced. In keeping with previous semi-analytic models, we place a seed black hole immediately after a galaxy forms. When the mass of the seed is set at 10^5 M_sun, we find that the model results become inconsistent with recent observational results of the black hole mass - bulge mass relation for dwarf galaxies. In particular, the model predicts that bulges with ~ 10^9 M_sun harbour larger black holes than observed. On the other hand, when we employ seed black holes with 10^3 M_sun, or randomly select their mass within a 10^(3-5) M_sun range, the resulting relation is consistent with observation estimates, including the observed dispersion. We find that to obtain stronger constraints on the mass of seed black holes, observations of less massive bulges at z ~ 0 are a more powerful comparison than the relations at higher redshifts.

A search for dark matter was conducted by looking for an annual modulation signal due to the Earth's rotation around the Sun using XMASS, a single phase liquid xenon detector. The data used for this analysis was 359.2 live days times 832 kg of exposure accumulated between November 2013 and March 2015. When we assume Weakly Interacting Massive Particle (WIMP) dark matter elastically scattering on the target nuclei, the exclusion upper limit of the WIMP-nucleon cross section 4.3$\times$10$^{-41}$cm$^2$ at 8 GeV/c$^2$ was obtained and we exclude almost all the DAMA/LIBRA allowed region in the 6 to 16 GeV/c$^2$ range at $\sim$10$^{-40}$cm$^2$. The result of a simple modulation analysis, without assuming any specific dark matter model but including electron/$\gamma$ events, showed a slight negative amplitude. The $p$-values obtained with two independent analyses are 0.014 and 0.068 for null hypothesis, respectively. we obtained 90\% C.L. upper bounds that can be used to test various models. This is the first extensive annual modulation search probing this region with an exposure comparable to DAMA/LIBRA.

We study the extent to which very bright (-23.0 < MUV < -21.75) Lyman-break selected galaxies at redshifts z~7 display detectable Lya emission. To explore this issue, we have obtained follow-up optical spectroscopy of 9 z~7 galaxies from a parent sample of 24 z~7 galaxy candidates selected from the 1.65 sq.deg COSMOS-UltraVISTA and SXDS-UDS survey fields using the latest near-infrared public survey data, and new ultra-deep Subaru z'-band imaging (which we also present and describe in this paper). Our spectroscopy has yielded only one possible detection of Lya at z=7.168 with a rest-frame equivalent width EW_0 = 3.7 (+1.7/-1.1) Angstrom. The relative weakness of this line, combined with our failure to detect Lya emission from the other spectroscopic targets allows us to place a new upper limit on the prevalence of strong Lya emission at these redshifts. For conservative calculation and to facilitate comparison with previous studies at lower redshifts, we derive a 1-sigma upper limit on the fraction of UV bright galaxies at z~7 that display EW_0 > 50 Angstrom, which we estimate to be < 0.23. This result may indicate a weak trend where the fraction of strong Lya emitters ceases to rise, and possibly falls between z~6 and z~7. Our results also leave open the possibility that strong Lya may still be more prevalent in the brightest galaxies in the reionization era than their fainter counterparts. A larger spectroscopic sample of galaxies is required to derive a more reliable constraint on the neutral hydrogen fraction at z~7 based on the Lya fraction in the bright galaxies.

We present the evolution of dark matter halos in six large cosmological N-body simulations, called the $\nu^2$GC (New Numerical Galaxy Catalog) simulations on the basis of the LCDM cosmology consistent with observational results obtained by the Planck satellite. The largest simulation consists of $8192^3$ (550 billion) dark matter particles in a box of $1.12 \, h^{-1} \rm Gpc$ (a mass resolution of $2.20 \times 10^{8} \, h^{-1} M_{\odot}$). Among simulations utilizing boxes larger than $1 \, h^{-1} \rm Gpc$, our simulation yields the highest resolution simulation that has ever been achieved. A $\nu^2$GC simulation with the smallest box consists of eight billions particles in a box of $70 \, h^{-1} \rm Mpc$ (a mass resolution of $3.44 \times 10^{6} \, h^{-1} M_{\odot}$). These simulations can follow the evolution of halos over masses of eight orders of magnitude, from small dwarf galaxies to massive clusters. Using the unprecedentedly high resolution and powerful statistics of the $\nu^2$GC simulations, we provide statistical results of the halo mass function, mass accretion rate, formation redshift, and merger statistics, and present accurate fitting functions for the Planck cosmology. By combining the $\nu^2$GC simulations with our new semi-analytic galaxy formation model, we are able to prepare mock catalogs of galaxies and active galactic nuclei, which will be made publicly available in the near future.

In this paper, we describe 220Rn calibration source that was developed for liquid noble gas detectors. The key advantage of this source is that it can provide 212Bi-212Po consecutive events, which enables us to evaluate the vertex resolution of a detector at low energy by comparing low-energy events of 212Bi and corresponding higher-energy alpha-rays from 212Po. Since 220Rn is a noble gas, a hot metal getter can be used when introduced using xenon as the carrier gas. In addition, no long-life radioactive isotopes are left behind in the detector after the calibration is complete; this has clear advantage over the use of 222Rn which leaves long- life radioactivity, i.e., 210Pb. Using a small liquid xenon test chamber, we developed a system to introduce 220Rn via the xenon carrier gas; we demonstrated the successful introduction of 6 times 10^2 220Rn atoms in our test environment.

We investigate morphological properties of 61 Lyman-alpha emitters (LAEs) at z = 4.86 identified in the COSMOS field, based on Hubble Space Telescope Advanced Camera for Surveys (ACS) imaging data in the F814W-band. Out of the 61 LAEs, we find the ACS counterparts for the 54 LAEs. Eight LAEs show double-component structures with a mean projected separation of 0."63 (~ 4.0 kpc at z = 4.86). Considering the faintness of these ACS sources, we carefully evaluate their morphological properties, that is, size and ellipticity. While some of them are compact and indistinguishable from the PSF half-light radius of 0."07 (~ 0.45 kpc), the others are clearly larger than the PSF size and spatially extended up to 0."3 (~ 1.9 kpc). We find that the ACS sources show a positive correlation between ellipticity and size and that the ACS sources with large size and round shape are absent. Our Monte Carlo simulation suggests that the correlation can be explained by (1) the deformation effects via PSF broadening and shot noise or (2) the source blending in which two or more sources with small separation are blended in our ACS image and detected as a single elongated source. Therefore, the 46 single-component LAEs could contain the sources which consist of double (or multiple) components with small spatial separation (i.e., < 0."3 or 1.9 kpc). Further observation with high angular resolution at longer wavelengths (e.g., rest-frame wavelengths of > 4000 A) is inevitable to decipher which interpretation is adequate for our LAE sample.

We investigate clustering properties of quasars using a new version of our semi-analytic model of galaxy and quasar formation with state-of-the-art cosmological N-body simulations. In this study, we assume that a major merger of galaxies triggers cold gas accretion on to a supermassive black hole and quasar activity. Our model can reproduce the downsizing trend of the evolution of quasars. We find that the median mass of quasar host dark matter haloes increases with cosmic time by an order of magnitude from z=4 (a few 1e+11 Msun) to z=1 (a few 1e+12 Msun), and depends only weakly on the quasar luminosity. Deriving the quasar bias through the quasar--galaxy cross-correlation function in the model, we find that the quasar bias does not depend on the quasar luminosity, similar to observed trends. This result reflects the fact that quasars with a fixed luminosity have various Eddington ratios and thus have various host halo masses that primarily determine the quasar bias. We also show that the quasar bias increases with redshift, which is in qualitative agreement with observations. Our bias value is lower than the observed values at high redshifts, implying that we need some mechanisms that make quasars inactive in low-mass haloes and/or that make them more active in high-mass haloes.

The Rosetta lander Philae successfully landed on the nucleus of comet 67P/Churyumov-Gerasimenko on 12 November 2014. Philae carries the Dust Impact Monitor (DIM) on board, which is part of the Surface Electric Sounding and Acoustic Monitoring Experiment (SESAME). DIM employs piezoelectric PZT sensors to detect impacts by sub-millimeter and millimeter-sized ice and dust particles that are emitted from the nucleus and transported into the cometary coma. The DIM sensor measures dynamical data like flux and the directionality of the impacting particles. Mass and speed of the particles can be constrained assuming density and elastic particle properties. DIM was operated during three mission phases of Philae at the comet: (1) Before Philae's separation from Rosetta at distances of about 9.6 km, 11.8 km, and 25.3 km from the nucleus barycenter. In this mission phase particles released from the nucleus on radial trajectories remained undetectable because of significant obscuration by the structures of Rosetta, and no dust particles were indeed detected. (2) During Philae's descent to its nominal landing site Agilkia, DIM detected one approximately millimeter-sized particle at a distance of 5.0 km from the nucleus' barycenter, corresponding to an altitude of 2.4 km from the surface. This is the closest ever dust detection at a cometary nucleus by a dedicated in-situ dust detector. (3) At Philae's final landing site, Abydos, DIM detected no dust impact which may be due to low cometary activity in the vicinity of Philae, or due to shading by obstacles close to Philae, or both. Laboratory calibration experiments showed that the material properties of the detected particle are compatible with a porous particle having a bulk density of approximately $250\, \mathrm{kg\,m^{-3}}$. The particle could have been lifted off from the comet's surface by sublimating water ice.

We present both the observations and the data reduction procedures of the Subaru COSMOS 20 project that is an optical imaging survey of the HST COSMOS field, carried out by using Suprime-Cam on the Subaru Telescope with the following 20 optical filters: 6 broad-band (B, g', V, r', i', and z'), 2 narrow-band (NB711 and NB816), and 12 intermediate-band filters (IA427, IA464, IA484, IA505, IA527, IA574, IA624, IA679, IA709, IA738, IA767, and IA827). A part of this project is described in Taniguchi et al. (2007) and Capak et al. (2007) for the six broad-band and one narrow-band (NB816) filter data. In this paper, we present details of the observations and data reduction for remaining 13 filters (the 12 IA filters and NB711). In particular, we describe the accuracy of both photometry and astrometry in all the filter bands. We also present optical properties of the Suprime-Cam IA filter system in Appendix.

We report a discovery of 6 massive galaxies with both extremely large Lya equivalent width and evolved stellar population at z ~ 3. These MAssive Extremely STrong Lya emitting Objects (MAESTLOs) have been discovered in our large-volume systematic survey for strong Lya emitters (LAEs) with twelve optical intermediate-band data taken with Subaru/Suprime-Cam in the COSMOS field. Based on the SED fitting analysis for these LAEs, it is found that these MAESTLOs have (1) large rest-frame equivalent width of EW_0(Lya) ~ 100--300 A, (2) M_star ~ 10^10.5--10^11.1 M_sun, and (3) relatively low specific star formation rates of SFR/M_star ~ 0.03--1 Gyr^-1. Three of the 6 MAESTLOs have extended Ly$\alpha$ emission with a radius of several kpc although they show very compact morphology in the HST/ACS images, which correspond to the rest-frame UV continuum. Since the MAESTLOs do not show any evidence for AGNs, the observed extended Lya emission is likely to be caused by star formation process including the superwind activity. We suggest that this new class of LAEs, MAESTLOs, provides a missing link from star-forming to passively evolving galaxies at the peak era of the cosmic star-formation history.

We have investigated effects of dust attenuation on quasar luminosity functions using a semi-analytic galaxy formation model combined with a large cosmological N-body simulation. We estimate the dust attenuation of quasars self-consistently with that of galaxies by considering the dust in their host bulges. We find that the luminosity of the bright quasars is strongly dimmed by the dust attenuation, about 2 mag in the B-band. Assuming the empirical bolometric corrections for active galactic nuclei (AGNs) by Marconi et al., we find that this dust attenuation is too strong to explain the B-band and X-ray quasar luminosity functions simultaneously. We consider two possible mechanisms that weaken the dust attenuation. As such a mechanism, we introduce a time delay for AGN activity, that is, gas fueling to a central black hole starts some time after the beginning of the starburst induced by a major merger. The other is the anisotropy in the dust distribution. We find that in order to make the dust attenuation of the quasars negligible, either the gas accretion into the black holes has to be delayed at least three times the dynamical timescale of their host bulges or the dust covering factor is as small as 0.1.

Current constraints on dark matter density profiles from weak lensing are typically limited to radial scales greater than 50-100 kpc. In this paper, we explore the possibility of probing the very inner regions of galaxy/halo density profiles by measuring stacked weak lensing on scales of only a few tens of kpc. Our forecasts focus on scales smaller than the equality radius (Req) where the stellar component and the dark matter component contribute equally to the lensing signal. We compute the evolution of Req as a function of lens stellar mass and redshift and show that Req=7-34 kpc for galaxies with the stellar mass of 10^{9.5}-10^{11.5} solar masses. Unbiased shear measurements will be challenging on these scales. We introduce a simple metric to quantify how many source galaxies overlap with their neighbours and for which shear measurements will be challenging. Rejecting source galaxies with close-by companions results in about a 20 per cent decrease in the overall source density. Despite this decrease, we show that Euclid and WFIRST will be able to constrain galaxy/halo density profiles at Req with signal-to-noise ratio >20 for the stellar mass of >10^{10} solar masses. Weak lensing measurements at Req, in combination with stellar kinematics on smaller scales, will be a powerful means by which to constrain both the inner slope of the dark matter density profile as well as the mass and redshift dependence of the stellar initial mass function.

Using the Hubble Space Telescope/Wide Field Camera 3 imaging data and multi-wavelength photometric catalog, we investigated the dust temperature of passively evolving and star-forming galaxies at 0.2<z<1.0 in the CANDELS fields. We estimated the stellar radiation field by low-mass stars from the stellar mass and surface brightness profile of these galaxies and then calculated their steady-state dust temperature. At first, we tested our method using nearby early-type galaxies with the deep FIR data by the Herschel Virgo cluster survey and confirmed that the estimated dust temperatures are consistent with the observed temperatures within the uncertainty. We then applied the method to galaxies at 0.2<z<1.0, and found that most of passively evolving galaxies with Mstar > 10^{10} Msun have a relatively high dust temperature of Tdust > 20 K, for which the formation efficiency of molecular hydrogen on the surface of dust grains in the diffuse ISM is expected to be very low from the laboratory experiments. The fraction of passively evolving galaxies strongly depends on the expected dust temperature at all redshifts and increases rapidly with increasing the temperature around Tdust ~ 20 K. These results suggest that the dust heating by low-mass stars in massive galaxies plays an important role for the continuation of their passive evolution, because the lack of the shielding effect of the molecular hydrogen on the UV radiation can prevent the gas cooling and formation of new stars.

We present a new chemical evolution model for dwarf spheroidal galaxies (dSphs) in the Local Universe. Our main aim is to explain both their observed star formation histories and metallicity distribution functions simultaneously. Applying our new model for the four local dSphs, that is, Fornax, Sculptor, Leo II, and Sextans, we find that our new model reproduces the observed chemical properties of the dSphs consistently. Our results show that the dSphs have evolved with both a low star formation efficiency and a large gas outflow efficiency compared with the Milky Way, as suggested by previous works. Comparing the observed [alpha/Fe]--[Fe/H] relation of the dSphs with the model predictions, we find that our model favors a longer onset time of Type Ia supernovae (i.e., 0.5 Gyr) than that suggested in previous studies (i.e., 0.1 Gyr). We discuss the origin of this discrepancy in detail.

Recent observations show that the space density of luminous active galactic nuclei (AGNs) peaks at higher redshifts than that of faint AGNs. This downsizing trend in the AGN evolution seems to be contradictory to the hierarchical structure formation scenario. In this study, we present the AGN space density evolution predicted by a semi-analytic model of galaxy and AGN formation based on the hierarchical structure formation scenario. We demonstrate that our model can reproduce the downsizing trend of the AGN space density evolution. The reason for the downsizing trend in our model is a combination of the cold gas depletion as a consequence of star formation, the gas cooling suppression in massive halos and the AGN lifetime scaling with the dynamical timescale. We assume that a major merger of galaxies causes a starburst, spheroid formation, and cold gas accretion onto a supermassive black hole (SMBH). We also assume that this cold gas accretion triggers AGN activity. Since the cold gas is mainly depleted by star formation and gas cooling is suppressed in massive dark halos, the amount of cold gas accreted onto SMBHs decreases with cosmic time. Moreover, AGN lifetime increases with cosmic time. Thus, at low redshifts, major mergers do not always lead to luminous AGNs. Because the luminosity of AGNs is correlated with the mass of accreted gas onto SMBHs, the space density of luminous AGNs decreases more quickly than that of faint AGNs. We conclude that the anti-hierarchical evolution of the AGN space density is not contradictory to the hierarchical structure formation scenario.

We present a new version of a semi-analytic model of cosmological galaxy formation, incorporating a star formation law with a feedback depending on the galaxy-scale mean dust opacity and metallicity, motivated by recent observations of star formation in nearby galaxies and theoretical considerations. This new model is used to investigate the effect of such a feedback on shaping the galaxy luminosity function and its evolution. Star formation activity is significantly suppressed in dwarf galaxies by the new feedback effect, and the faint-end slope of local luminosity functions can be reproduced with a reasonable strength of supernova feedback, which is in contrast to the previous models that require a rather extreme strength of supernova feedback. Our model can also reproduce the early appearance of massive galaxies manifested in the bright-end of high redshift K-band luminosity functions. Though some of the previous models also succeeded in reproducing this, they assumed a star formation law depending on the galaxy-scale dynamical time, which is not supported by observations. We argue that the feedback depending on dust opacity (or metal column density) is essential, rather than that simply depending on gas column density, to get these results.

Using the Hubble Space Telescope/Advanced Camera for Surveys data in the COSMOS field, we systematically searched clumpy galaxies at 0.2<z<1.0 and investigated the fraction of clumpy galaxies and its evolution as a function of stellar mass, star formation rate (SFR), and specific SFR (SSFR). The fraction of clumpy galaxies in star-forming galaxies with Mstar > 10^9.5 Msun decreases with time from ~0.35 at 0.8<z<1.0 to ~0.05 at 0.2<z<0.4 irrespective of the stellar mass, although the fraction tends to be slightly lower for massive galaxies with Mstar > 10^10.5 Msun at each redshift. On the other hand, the fraction of clumpy galaxies increases with increasing both SFR and SSFR in all the redshift ranges we investigated. In particular, we found that the SSFR dependences of the fractions are similar among galaxies with different stellar masses, and the fraction at a given SSFR does not depend on the stellar mass in each redshift bin. The evolution of the fraction of clumpy galaxies from z~0.9 to z~0.3 seems to be explained by such SSFR dependence of the fraction and the evolution of SSFRs of star-forming galaxies. The fraction at a given SSFR also appears to decrease with time, but this can be due to the effect of the morphological K-correction. We suggest that these results are understood by the gravitational fragmentation model for the formation of giant clumps in disk galaxies, where the gas mass fraction is a crucial parameter.

In order to understand nature of building blocks of galaxies in the early universe, we investigate "genuine irregular galaxies (GIGs)" in the nearby universe. Here, GIGs are defined as isolated galaxies without regular structures (spheroid, bulge, disk, bar, spiral arm, and nucleus). Using the results of two excellent studies on galaxy morphology based on the Sloan Digital Sky Survey (SDSS), we obtain a sample of 66 irregular galaxies. We carry out new classification of them into GIGs and non-GIGs which have regular structure or show evidence for galaxy interaction, by using the SDSS Data Release 10 images. We then find that a half of these irregular galaxies (33/66) are GIGs and obtain an unambiguous sample of 33 GIGs for the first time. We discuss their observational properties by comparing them with those of elliptical, S0, spiral galaxies, and irregular galaxies without the GIGs. We find that our GIGs have smaller sizes, lower optical luminosities, bluer rest-frame optical colors, lower surface stellar mass densities, and lower gas metallicity than normal galaxies. All these properties suggest that they are in chemically and dynamically younger phases even in the nearby universe.

We present a new model of the extragalactic background light (EBL) and corresponding gamma-gamma opacity for intergalactic gamma-ray absorption from z = 0 up to z = 10, based on a semi-analytical model of hierarchical galaxy formation that reproduces key observed properties of galaxies at various redshifts. Including the potential contribution from Population III stars and following the cosmic reionization history in a simplified way, the model is also broadly consistent with available data concerning reionization, particularly the Thomson scattering optical depth constraints from WMAP. In comparison with previous EBL studies up to z ~ 3-5, our predicted gamma-gamma opacity is in general agreement for observed gamma-ray energy below 400/(1 + z) GeV, whereas it is a factor of ~ 2 lower above this energy because of a correspondingly lower cosmic star formation rate, even though the observed UV luminosity is well reproduced by virtue of our improved treatment of dust obscuration and direct estimation of star formation rate. The horizon energy at which the gamma-ray opacity is unity does not evolve strongly beyond z ~ 4 and approaches ~ 20 GeV. The contribution of Population III stars is a minor fraction of the EBL at z = 0, and is also difficult to distinguish through gamma-ray absorption in high-z objects, even at the highest levels allowed by the WMAP constraints. Nevertheless, the attenuation due to Population II stars should be observable in high-z gamma-ray sources by telescopes such as Fermi or CTA and provide a valuable probe of the evolving EBL in the rest-frame UV. The detailed results of our model are publicly available in numerical form at the URL http://www.slac.stanford.edu/~yinoue/Download.html.

During the course of our Lyman continuum imaging survey, we found that the spectroscopically confirmed Lya emitter LAE 221724+001716 at z = 3.10 in the SSA 22 field shows strong Lyman continuum emission (lambda_rest ~ 900 A) that escapes from this galaxy. However, another recent spectroscopic survey revealed that the supposed Lyman continuum emission could arise from a foreground galaxy at z = 1.76 if the emission line newly detected from the galaxy at lambda_obs ~ 3360 A is Lya. If this is the case, as the angular separation between these two galaxies is very small (~ 0.6"), LAE 221724+001716 at z = 3.10 could be amplified by the gravitational lensing caused by this intervening galaxy. Here we present a possible gravitational lens model for the system of LAE 221724+001716. First, we estimate the stellar mass of the intervening galaxy as Mstar ~ 3.5x10^9 Msun from its UV luminosity and ~ 3.0x10^7--2.4x10^9 Msun through the SED fitting. Then, we find that the gravitational magnification factor ranges from 1.01 to 1.16 using the so-called singular isothermal sphere model for strong lensing. While LAE 221724+001716 is the first system of an LAE-LAE lensing reported so far, the estimated magnification factor is not so significant because the stellar mass of the intervening galaxy is small.

The cosmic star formation rate density (CSFRD) has been observationally investigated out to redshift z~10. However, most of theoretical models for galaxy formation underpredict the CSFRD at z>1. Since the theoretical models reproduce the observed luminosity functions (LFs), luminosity densities (LDs), and stellar mass density at each redshift, this inconsistency does not simply imply that theoretical models should incorporate some missing unknown physical processes in galaxy formation. Here, we examine the cause of this inconsistency in UV wavelengths by using a mock catalog of galaxies generated by a semi-analytic model of galaxy formation. We find that this inconsistency is due to two observational uncertainties: dust obscuration correction and conversion from UV luminosity to star formation rate (SFR). The methods for correction of obscuration and SFR conversion used in observational studies result in the overestimation of CSFRD by ~ 0.1-0.3 dex and ~ 0.1-0.2 dex, respectively, compared to the results obtained directly from our mock catalog. We present new empirical calibrations for dust attenuation and conversion from observed UV LFs and LDs into CSFRD.

We have studied the properties of more than 1600 low-redshift galaxies by utilizing high-quality infrared flux measurements of the AKARI All-Sky Survey and physical quantities based on optical and 21-cm observations. Our goal is to understand the physics determining the infrared spectral energy distribution (SED). The ratio of the total infrared luminosity L_TIR, to the star-formation rate (SFR) is tightly correlated by a power-law to specific SFR (SSFR), and L_TIR is a good SFR indicator only for galaxies with the largest SSFR. We discovered a tight linear correlation for normal galaxies between the radiation field strength of dust heating, estimated by infrared SED fits (U_h), and that of galactic-scale infrared emission (U_TIR ~ L_TIR/R^2), where R is the optical size of a galaxy. The dispersion of U_h along this relation is 0.3 dex, corresponding to 13% dispersion in the dust temperature. This scaling and the U_h/U_TIR ratio can be explained physically by a thin layer of heating sources embedded in a thicker, optically-thick dust screen. The data also indicate that the heated fraction of the total dust mass is anti-correlated to the dust column density, supporting this interpretation. In the large U_TIR limit, the data of circumnuclear starbursts indicate the existence of an upper limit on U_h, corresponding to the maximum SFR per gas mass of ~ 10 Gyr^{-1}. We find that the number of galaxies sharply drops when they become optically thin against dust-heating radiation, suggesting that a feedback process to galaxy formation (likely by the photoelectric heating) is working when dust-heating radiation is not self-shielded on a galactic scale. Implications are discussed for the M_HI-size relation, the Kennicutt-Schmidt relation, and galaxy formation in the cosmological context.

We present a new theoretical calculation of the contribution to the extragalactic gamma-ray background radiation (EGRB) from star-forming galaxies, based on a state-of-the-art model of hierarchical galaxy formation that is in quantitative agreement with a variety of observations of local and high-redshift galaxies. Gamma-ray luminosity ($L_\gamma$) and spectrum of galaxies are related to star formation rate ($\psi$), gas mass ($M_gas$), and star formation mode (quiescent or starburst) of model galaxies using latest observed data of nearby galaxies. We try the two limiting cases about gamma-ray production: the escape limit ($L_\gamma \propto \psi M_gas$) and the calorimetric limit ($L_\gamma \propto \psi$), and our standard model predicts 7 and 4% contribution from star-forming galaxies to the total EGRB flux (including bright resolved sources) recently reported by the Fermi Gamma-Ray Space Telescope. Systematic uncertainties do not allow us to determine the EGRB flux better than by a factor of ~2. The predicted number of nearby galaxies detectable by Fermi is consistent with the observation. Intergalactic absorption by pair-production attenuates the EGRB flux only by a modest factor of ~1.3 at the highest Fermi energy band, and the reprocessed cascade emission does not significantly alter EGRB at lower photon energies. The sum of the known contributions from AGNs and star-forming galaxies can explain a large part of EGRB, with a remarkable agreement between the predicted model spectrum and observation.

We study the relationship between the metallicity of gamma-ray burst (GRB) progenitors and the probability distribution function (PDF) of GRB host galaxies as a function of luminosity using cosmological hydrodynamic simulations of galaxy formation. We impose a maximum limit to the gas metallicity in which GRBs can occur, and examine how the predicted luminosity PDF of GRB host galaxies changes in the simulation. We perform the Kolmogorov-Smirnov test, and show that the result from our simulation agrees with the observed luminosity PDF of core-collapse supernovae (SNe) host galaxies when we assume that the core-collapse SNe trace star formation. When we assume that GRBs occur only in a low-metallicity environment with $Z\lesssim 0.1 \Zsun$, GRBs occur in lower luminosity galaxies, and the simulated luminosity PDF becomes quantitatively consistent with the observed luminosity PDF. The observational bias against the host galaxies of optically dark GRBs owing to dust extinction may be another reason for the lower luminosities of GRB host galaxies, but the observed luminosity PDF of GRB host galaxies cannot be reproduced solely by the dust bias in our simulation.

We investigate the prospects for discovering blazars at very high-redshifts (z>3-6) with the Fermi Gamma-Ray Space Telescope (Fermi), employing a model for the evolving gamma-ray luminosity function (GLF) of the blazar population. Our previous GLF model is used as a basis, which features luminosity-dependent density evolution implied from X-ray data on active galactic nuclei, as well as the blazar sequence paradigm for their spectral energy distribution, and which is consistent with EGRET and current Fermi observations of blazars.Here we augment the high-redshift evolution of this model by utilizing the luminosity function of quasars from the Sloan Digital Sky Survey (SDSS), which is well-constrained up to z~5. We find that Fermi may discover a few blazars up to z~6 in the entire sky during its 5-year survey. We further discuss how such high-redshift blazar candidates may be efficiently selected in future Fermi data.

We conducted a deep narrowband NB973 (FWHM = 200 A centered at 9755 A) survey of z=7 Lyman alpha emitters (LAEs) in the Subaru/XMM-Newton Deep Survey Field, using the fully depleted CCDs newly installed on the Subaru Telescope Suprime-Cam, which is twice more sensitive to z=7 Lyman alpha at ~ 1 micron than the previous CCDs. Reaching the depth 0.5 magnitude deeper than our previous survey in the Subaru Deep Field that led to the discovery of a z=6.96 LAE, we detected three probable z=7 LAE candidates. Even if all the candidates are real, the Lyman alpha luminosity function (LF) at z=7 shows a significant deficit from the LF at z=5.7 determined by previous surveys. The LAE number and Lyman alpha luminosity densities at z=7 is ~ 7.7-54% and ~5.5-39% of those at z=5.7 to the Lyman alpha line luminosity limit of L(Ly-alpha) >~ 9.2 x 10^{42} erg s^{-1}. This could be due to evolution of the LAE population at these epochs as a recent galaxy evolution model predicts that the LAE modestly evolves from z=5.7 to 7. However, even after correcting for this effect of galaxy evolution on the decrease in LAE number density, the z=7 Lyman alpha LF still shows a deficit from z=5.7 LF. This might reflect the attenuation of Lyman alpha emission by neutral hydrogen remaining at the epoch of reionization and suggests that reionization of the universe might not be complete yet at z=7. If we attribute the density deficit to reionization, the intergalactic medium (IGM) transmission for Lyman alpha photons at z=7 would be 0.4 <= T_{Ly-alpha}^{IGM} <= 1, supporting the possible higher neutral fraction at the earlier epochs at z > 6 suggested by the previous surveys of z=5.7-7 LAEs, z ~ 6 quasars and z > 6 gamma-ray bursts.

A rich variety of order parameter manifolds of multicomponent Bose-Einstein condensates (BECs) admit various kinds of topological excitations, such as fractional vortices, monopoles, skyrmions, and knots. In this paper, we discuss two topological excitations in spinor BECs: non-Abelian vortices and knots. Unlike conventional vortices, non-Abelian vortices neither reconnect themselves nor pass through each other, but create a rung between them in a topologically stable manner. We discuss the collision dynamics of non-Abelian vortices in the cyclic phase of a spin-2 BEC. In the latter part, we show that a knot, which is a unique topological object characterized by a linking number or a Hopf invariant [$\pi_3 (S^2)=Z$], can be created using a conventional quadrupole magnetic field in a cold atomic system.

The latest determination of the extragalactic gamma-ray background (EGRB) radiation by Fermi is compared with the theoretical prediction of the blazar component by Inoue & Totani (2009; hereafter IT09). The Fermi EGRB spectrum is in excellent agreement with IT09, indicating that blazars are the dominant component of the EGRB, and contributions from any other sources (e.g., dark matter annihilations) are minor. It also indicates that the blazar SED (spectral energy distribution) sequence taken into account in IT09 is a valid description of mean blazar SEDs. The possible contribution of MeV blazars to the EGRB in the MeV band is also discussed. In five total years of observations, we predict that Fermi will detect ~1200 blazars all sky down to the corresponding sensitivity limit. We also address the detectability of the highest-redshift blazars. Updating our model with regard to high-redshift evolution based on SDSS quasar data, we show that Fermi may find some blazars up to z~6 during the five-year survey. Such blazars could provide a new probe of early star and galaxy formation through GeV spectral attenuation signatures induced by high-redshift UV background radiation.

We present theoretical predictions of UV continuum luminosity function (UV LF) and Lya equivalent width (EW) distribution of Lyman alpha emitters (LAEs) in the framework of the hierarchical clustering model of galaxy formation. The model parameters about LAEs were determined by fitting to the observed Lya LF at z=5.7 in our previous study, and the fit indicates that extinction of Lya photons by dust is significantly less effective than that of UV continuum photons, implying clumpy dust distribution in interstellar medium. We then compare the predictions about UV LFs and EW distributions with a variety of observations at z ~ 3-6, allowing no more free parameters and paying careful attention to the selection conditions of LAEs in each survey. We find that the predicted UV LFs and EW distributions are in nice agreement with observed data, and especially, our model naturally reproduces the existence of large EW LAEs (> 240 A) without introducing Pop III stars or top-heavy initial mass function. We show that both the stellar population (young age and low metallicity) and extinction by clumpy dust are the keys to reproduce large EW LAEs. The evidence of EW enhancement by clumpy dust is further strengthened by the quantitative agreement between our model and recent observations about a positive correlation between EW and extinction. The observed trend that brighter LAEs in UV continuum tend to have smaller mean EW is also reproduced, and the clumpy dust is playing an important role again for this trend. We suggested in our previous study that the transmission of intergalactic medium for Lya emission rapidly decreases from z ~ 6 to 7 by the fitting to Lya LFs, and this evidence is quantitatively strengthened by the comparison with the UV LF and EW distribution at z ~ 6.6.

Long gamma-ray bursts (GRBs) have been suggested to occur preferentially in low metallicity environment. We discuss the possibility and theoretical aspects of using Lyman alpha emission properties of long GRB host galaxies as a metallicity indicator of high redshift GRB environments, where direct metallicity measurements are not easy. We propose to use the fraction of Ly-alpha emitters (LAEs) in long GRB host galaxies as a function of UV luminosity, which can be compared with star-formation-rate weighted LAE fraction of Lyman-break galaxies as the standard in the case of no metallicity dependence. There are two important effects of metallicity dependence of long GRB rate to change the LAE fraction of host galaxies. One is the enhancement of intrinsic Ly-alpha equivalent width (EW) by stronger ionizing UV luminosity of low metallicity stellar population, and the other is extinction by interstellar dust to change the observable EW. Based on a latest theoretical model of LAEs that reproduce observations, we argue that the latter is likely to work in the opposite direction to the former, i.e., to decrease LAE fraction if GRBs preferentially occur in low-metallicity environments, because of the clumpy inter-stellar medium effect. The high LAE fraction of GRB host galaxies indicated by observations is quantitatively be explained by the LAE model if GRBs occur when Z < 0.1Z_sun, although this result is still indicative because of the limited statistics and theoretical uncertainties. This result demonstrates that the LAE statistics of GRB hosts may give us useful information in the future.

We investigate the collision dynamics of two non-Abelian vortices and find that, unlike Abelian vortices, they neither reconnect themselves nor pass through each other, but create a rung between them in a topologically stable manner. Our predictions are verified using the model of the cyclic phase of a spin-2 spinor Bose-Einstein condensate.

We made a narrowband NB973 (bandwidth of 200A at 9755A) imaging of the Subaru Deep Field (SDF) and found two z=7 Lyman alpha emitter (LAE) candidates down to NB973=24.9. Carrying out deep follow-up spectroscopy, we identified one of them as a real z=6.96 LAE. This has shown that galaxy formation was in progress just 750 Myr after the Big Bang. Meanwhile, the Lyman alpha line luminosity function of LAE is known to decline from z=5.7 to 6.6 in the SDF. L* at z=6.6 is 40-60% of that at z=5.7. We also confirm that the number density of z=7 LAE is only 17% of the density at z=6.6 comparing the latest SDF LAE samples. This series of significant decreases in LAE density with increasing redshift can be the result of galaxy evolution during these epochs. However, using the UV continuum luminosity functions of LAEs, those of Lyman break galaxies, and a LAE evolution model based on the hierarchical clustering, we find that galaxy evolution alone cannot explain all the decrease in density. This extra density deficit can be interpreted as the attenuation of the Lyman alpha photons from LAEs due to a rapid evolution of neutral hydrogen fraction during the ongoing cosmic reionization at z~6.6-7.

We present a new theoretical model for the luminosity functions (LFs) of Lyman alpha (Lya) emitting galaxies in the framework of hierarchical galaxy formation. We extend a semi-analytic model of galaxy formation that reproduces a number of observations for local and high-z galaxies, without changing the original model parameters but introducing a physically-motivated modelling to describe the escape fraction of Lya photons from host galaxies (f_esc). Though a previous study using a hierarchical clustering model simply assumed a constant and universal value of f_esc, we incorporate two new effects on f_esc: extinction by interstellar dust and galaxy-scale outflow induced as a star formation feedback. It is found that the new model nicely reproduces all the observed Lya LFs of the Lya emitters (LAEs) at different redshifts in z ~ 3-6. Especially, the rather surprisingly small evolution of the observed LAE Lya LFs compared with the dark halo mass function is naturally reproduced. Our model predicts that galaxies with strong outflows and f_esc ~ 1 are dominant in the observed LFs. This is also consistent with available observations, while the simple universal f_esc model requires f_esc << 1 not to overproduce the brightest LAEs. On the other hand, we found that our model significantly overpredicts LAEs at z > 6, and absorption of Lya photons by neutral hydrogen in intergalactic medium (IGM) is a reasonable interpretation for the discrepancy. This indicates that the IGM neutral fraction x_HI rapidly evolves from x_HI << 1 at z < 6 to a value of order unity at z ~ 6-7, which is broadly consistent with other observational constraints on the reionization history.

LXeGRIT - Liquid Xenon Gamma-Ray Imaging Telescope - is the first prototype of a Compton telescope for \MeV \g-ray astrophysics based on a LXe time projection chamber. One of the most relevant figures of merit for a Compton telescope is the detection efficiency for \g-rays, which depends on diverse contributions such as detector geometry and passive materials, trigger efficiency, dead time, etc. A detailed study of the efficiency of the LXeGRIT instrument, based both on laboratory measurements and Monte Carlo simulations, is presented in this paper.

To explore the possibility that an inflationary universe can be created out of a stable particle in the laboratory, we consider the classical and quantum dynamics of a magnetic monopole in the thin-shell approximation. Classically there are three types of solutions: stable, collapsing and inflating monopoles. We argue that the transition from a stable monopole to an inflating one could occur either by collision with a domain wall or by quantum tunneling.

The first galaxies in the Universe are built up where cold dark matter (CDM) forms large scale filamentary structure. Although the galaxies are expected to emit numerous Lya photons, they are surrounded by plentiful neutral hydrogen with a typical optical depth for Lya of ~10^5 (HI halos) before the era of cosmological reionization. The HI halo almost follows the cosmological Hubble expansion with some anisotropic corrections around the galaxy because of the gravitational attraction by the underlying CDM filament. In this paper, we investigate the detectability of the Lya emissions from the first galaxies, examining their dependence on viewing angles. Solving the Lya line transfer problem in an anisotropically expanding HI halo, we show that the escape probability from the HI halo is the largest in direction along the filament axis. If the Lya source is observed with a narrow-band filter, the difference of apparent Lya line luminosities among viewing angles can be a factor of > 40 at an extreme case. Furthermore, we evaluate the predicted physical features of the Lya sources and flux magnification by gravitational lensing effect due to clusters of galaxies along the filament. We conclude that, by using the next generation space telescopes like the JWST, the Lya emissions from the first galaxies whose CDM filament axes almost face to us can be detected with the S/N of > 10.

The expected episodic or intermittent star formation histories (SFHs) of gas rich dwarf irregular galaxies (dIrrs) are the longstanding puzzles to understand their whole evolutional history. Solving this puzzle, we should grasp what physical mechanism causes the quiescent phase of star formation under the very gas rich condition after the first starburst phase. We consider that this quiescent phase is kept by lack of H2, which can be important coolant to generate the next generation of stars in the low-metal environment like dIrrs. Furthermore, in dIrrs, H2 formation through gas-phase reactions may dominate the one on dust-grain surfaces because their interstellar medium (ISM) are very plentiful and the typical dust-to-gas ratio of dIrrs (D_dIrrs = 1.31 x 10^-2 D_MW, where D_MW is its value for the local ISM) is on the same order with a critical value D_cr ~ 10^-2 D_MW. We show that the lack of H2 is mainly led by H- destruction when gas-phase H2 formation dominates since H- is important intermediary of gas-phase H2 formation. H- is destroyed by the radiation from all stars born in the previous starburst phase because H- destroying infrared photon can penetrate the whole ISM of dIrrs. Considering the physical process which timescale is the longest as main process in regulating global star formation, we can show this lack of H2 leads the quiescent phase of star formation. Hence, we can say that the stellar radiation which destroys H- and leads low H2 abundance should be properly treated in studying SFHs of dIrrs.

Ly$\alpha$ line emission from the Loeb-Rybicki (LR) halo, which is the expanding HI IGM (intergalactic medium) around the first star clusters and the ionized interstellar medium, is investigated by solving a radiative transfer problem. While the initial scattering optical depth is $\sim 10^5$ for the Ly$\alpha$ photons, most of the Ly$\alpha$ photons can escape when the cumulative frequency-shift due to the expansion of the HI IGM becomes significantly large. The current paper improves upon previous treatments of the scattering processes and the opacity for the Ly$\alpha$ transfer. Confirming the previous results of the LR halo, we investigate the effect of the aspherical expansion of the IGM. The asphericity is hypothesized to follow the initial stage of the gravitational deceleration to form the large scale filamentary structure of the Universe. According to our results, the effect of the asphericity lets the peak wavelength of the line profile shift to longer wavelengths and the FWHM of the profile become wider than those of the spherically expanding model. To detect these features is meaningful if we are interested in the initial evolution of the large scale structure, since they reflect the dynamical properties of the IGM at that time. Furthermore, given the recent discovery of the high redshift cosmological reionization, we briefly comment on the effects of the redshift and the cosmological parameters on the line profile.

LXeGRIT is the first prototype of a novel Compton telescope for MeV gamma-ray astrophysics based on a Liquid Xenon Time Projection Chamber (LXeTPC), sensitive in the energy band of 0.15-10 MeV. In this homogeneous, 3D position sensitive detector, gamma rays with at least two interactions in the sensitive volume of 2800 cm$^{3}$, are imaged as in a standard Compton telescope. Gamma-rays with a single interaction cannot be imaged and constitute a background which can be easily identified and rejected. Charged particles and localized beta-particles background is also easily suppressed based on the TPC localization capability with millimeter resolution. A measurement of the total gamma-ray background rate in near space conditions and the background rejection power of the LXeTPC was a primary goal of the LXeGRIT balloon flight program. We present here a preliminary analysis addressing this question, based on balloon flight data acquired during the Oct 4-5, 2000 LXeGRIT balloon flight from Ft. Sumner, NM. In this long duration (27 hr) balloon experiment, the LXeGRIT TPC was not surrounded by any gamma-ray or charged particle shield. Single site events and charged particles were mostly rejected on-line at the first and second trigger level. The remaining count rate of single-site \g-ray events, at an average atmospheric depth of 3.2 g cm$^{-2}$, is consistent with that expected from atmospheric and diffuse gamma-ray background, taking into account the instrument mass model and response.

LXeGRIT is the first prototype of a novel concept of Compton telescope, based on the complete 3-dimensional reconstruction of the sequence of interactions of individual gamma rays in one position sensitive detector. This balloon-borne telescope consists of an unshielded time projection chamber with an active volume of 400 cm$^2 \times 7$ cm filled with high purity liquid xenon. Four VUV PMTs detect the fast xenon scintillation light signal, providing the event trigger. 124 wires and 4 anodes detect the ionization signals, providing the event spatial coordinates and total energy. In the period 1999 -- 2001, LXeGRIT has been extensively tested both in the laboratory and at balloon altitude, and its response in the MeV region has been thoroughly characterized. Here we summarize some of the results on pre-flight calibration, event reconstruction techniques, and performance during a 27 hour balloon flight on October 4 -- 5. We further present briefly the on-going efforts directed to improve the performance of this prototype towards the requirements for a base module of a next-generation Compton telescope.

XENON is a novel liquid xenon experiment concept for a sensitive dark matter search using a 1-tonne active target, distributed in an array of ten independent time projection chambers. The design relies on the simultaneous detection of ionization and scintillation signals in liquid xenon, with the goal of extracting as much information as possible on an event-by-event basis, while maintaining most of the target active. XENON is expected to have effective and redundant background identification and discrimination power, higher than 99.5%, and to achieve a very low threshold, on the order of 4 keV visible recoil energy. Based on this expectation and the 1-tonne mass of active xenon, we project a sensitivity of 0.0001 events/kg/day, after 3 yr operation in an appropriate underground location. The XENON experiment has been recently proposed to the National Science Foundation (NSF) for an initial development phase leading to the development of the 100 kg unit module.