Search

Understanding the survival of gas within subhalos under various astrophysical processes is crucial for elucidating cosmic structure formation and evolution. We study the resilience of gas in subhalos, focusing on the impact of tidal and ram pressure stripping through hydrodynamic simulations. Our results uncover significant gas stripping primarily driven by ram pressure effects, which also profoundly influence the gas distribution within these subhalos. Notably, despite their vulnerability to ram pressure effects, the low-mass subhalos can play a pivotal role in influencing the observable characteristics of cosmic structures due to their large abundance. Specifically, we explore the application of our findings to the 21 cm forest, showing how the survival dynamics of gas in subhalos can modulate the 21 cm optical depth, a key probe for detecting minihalos in the pre-reionization era. (abridged) In this work, we further investigate the contribution of subhalos to the 21 cm optical depth with hydrodynamics simulations, particularly highlighting the trajectories and fates of subhalos within mass ranges of \(10^{4-6} M_{\odot}h^{-1}\) in a host halo of \(10^7 M_{\odot}h^{-1}\). Despite their susceptibility to ram pressure stripping, the contribution of abundant low-mass subhalos to the 21-cm optical depth is more significant than that of their massive counterparts primarily due to their greater abundance. We find that the 21-cm optical depth can be increased by a factor of approximately two due to the abundant low-mass subhalos. (abridged) Our work provides critical insights into the gas dynamics within subhalos in the early Universe, highlighting their resilience against environmental stripping effects, and their impact on observable 21-cm signals.

On the basis of standard wavelength conversion by the use of angular phase matching of nonlinear optical crystals arranged in free space, applicable to a wide range of wavelengths and laser intensities, we both numerically and experimentally present an optimal design for achieving near-full energy conversion while maintaining good single-mode properties of fundamental laser radiation in the nanosecond regime., Comment: 9 pages, including 4 figures

The dark age of the universe, when no luminous object had existed, ended with the birth of the first stars, galaxies, and blackholes. This epoch is called cosmic dawn. Cosmic reionization is the major transition of the intergalactic medium (IGM) in the universe driven by ionizing photons emitted from luminous objects. Although the epoch through the dark age to reionization is a milestone in the universe, our knowledge of this epoch has not been sufficient yet. Cosmic 21cm signal, which is emitted from neutral hydrogen, is expected to open a new window for this epoch. In this review paper, we first introduce the basic physics of the 21cm line and how first stars impact on the 21cm line signal. Next, we briefly summarize how we extract astrophysical information from the 21cm line signal by means of statistical and machine learning approaches. We also discuss the synergy between the 21cm line signal and other emission lines. Finally, we summarize the current status of 21cm experiments., Comment: 30 pages, 19 figures. Published in PASJ as Special Issue "Metre and Centimetre Radio Astronomy in the Next Decade"

The tomography of the polarized Sunyaev-Zeldvich effect due to free electrons of galaxy clusters can be used to constrain the nature of dark energy because CMB quadrupoles at different redshifts as the polarization source are sensitive to the integrated Sachs-Wolfe effect. Here we show that the low multipoles of the temperature and E-mode polarization anisotropies from the all-sky CMB can improve the constraint further through the correlation between them and the CMB quadrupoles viewed from the galaxy clusters. Using a Monte-Carlo simulation, we find that low multipoles of the temperature and E-mode polarization anisotropies potentially improve the constraint on the equation of state of dark energy parameter by $\sim 17$ percent.

We study the contribution of subhalos to the 21 cm forest signal. The halos can host the substructures and including the effects of those small scale clumps can potentially boost the 21 cm optical depth in favor of detecting the 21 cm forest signals. We estimate the boost factor representing the ratio of the optical depth due to the subhalo contribution and that due to the host halo alone (without subhalos). Even though the optical depth boost factor is negligible for a small host halo with the mass of order $10^5 M_{\odot}$, the subhalo contribution can enhance the optical depth by an order of magnitude for a host halo of order $10^7 M_{\odot}$. The resultant 21 cm absorption line abundance which is obtained by integrating over the halo mass range relevant for the 21 cm forest signal can be enhanced by up to of order $10\%$ due to the substructures. The larger boost factor for a larger host halo would be of particular interest for the 21 cm forest detection because the the contribution of the larger host halos to the 21 cm forest signals is smaller due to their higher temperature and less abundance than the smaller host halos. The subhalos hence can well help the larger host halos more important for the signal estimation which, without considering the subhalos, may not give appreciable contribution to 21 cm forest signals., Comment: 14 pages, 7 figures

Cosmic Reionization imprints its signature on the temperature and polarization anisotropies of the cosmic microwave background (CMB). Advances in CMB telescopes have already placed a significant constraint on the history of reionization. As near-future CMB telescopes target the maximum sensitivity, or observations limited only by the cosmic variance (CV), we hereby forecast the potential of future CMB observations in constraining the history of reionization. In this study, we perform Markov Chain Monte Carlo analysis for CV-limited E-mode polarization observations such as the LiteBIRD (Light satellite for the studies of B-mode polarization and Inflation from cosmic background Radiation Detection), based on a few different methods that vary in the way of sampling reionization histories. We focus especially on estimating the very early history of reionization that occurs at redshifts $z>15$, which is quantified by the partial CMB optical depth due to free electrons at $z>15$, $\tau_{z>15}$. We find that reionization with $\tau_{z>15} \sim 0.008$, which are well below the current upper limit $\tau_{z>15} \sim 0.02$, are achievable by reionization models with minihalo domination in the early phase and can be distinguished from those with $\tau_{z>15} \lesssim 5\times 10^{-4}$ through CV-limited CMB polarization observations. An accurate estimation of $\tau_{z>15}$, however, remains somewhat elusive. We investigate whether resampling the E-mode polarization data with limited spherical-harmonic modes may resolve this shortcoming., Comment: 13 pages, 18figures

Theoretically modelling the 21-cm signals caused by Population III stars (Pop III stars) is the key to extracting fruitful information on Pop III stars from current and forthcoming 21-cm observations. In this work we develop a new module of Pop III stars in which the escape fractions of ionizing photons and Lyman-Werner (LW) photons, photo-heating by UV radiation, and LW feedback are consistently incorporated. By implementing the module into a public 21-cm semi-numerical simulation code, 21CMFAST, we demonstrate 21-cm signal calculations and investigate the importance of Pop III star modelling. What we find is that the contribution from Pop III stars to cosmic reionization significantly depends on the treatment of the escape fraction. With our escape fraction model, Pop III stars hardly contribute to reionization because less massive halos, whose escape fraction are high, cannot host Pop III stars due to LW feedback. On the other hand, Pop III stars well contribute to reionization with the conventional constant escape fraction. We also find that UV photo-heating has non-negligible impact on the 21-cm global signal and the 21-cm power spectrum if the ionization fraction of the Universe is higher than roughly 1 percent. In this case, the strength of the 21-cm global signal depends on the photo-heating efficiency and thus on the Pop III star mass. We conclude that detailed modelling of Pop III stars is imperative to predict 21-cm observables accurately for future observations., Comment: 9 pages, 8 figures, submitted to MNRAS

Precise instrumental calibration is of crucial importance to 21-cm cosmology experiments. The Murchison Widefield Array's (MWA) Phase II compact configuration offers us opportunities for both redundant calibration and sky-based calibration algorithms; using the two in tandem is a potential approach to mitigate calibration errors caused by inaccurate sky models. The MWA Epoch of Reionization (EoR) experiment targets three patches of the sky (dubbed EoR0, EoR1, and EoR2) with deep observations. Previous work in \cite{Li_2018} and \cite{Wenyang_2019} studied the effect of tandem calibration on the EoR0 field and found that it yielded no significant improvement in the power spectrum over sky-based calibration alone. In this work, we apply similar techniques to the EoR1 field and find a distinct result: the improvements in the power spectrum from tandem calibration are significant. To understand this result, we analyze both the calibration solutions themselves and the effects on the power spectrum over three nights of EoR1 observations. We conclude that the presence of the bright radio galaxy Fornax A in EoR1 degrades the performance of sky-based calibration, which in turn enables redundant calibration to have a larger impact. These results suggest that redundant calibration can indeed mitigate some level of model-incompleteness error., Comment: 20 pages, 11 figures, accepted for publication in PASA

The radio observation of 21\,cm-line signal from the Epoch of Reionization (EoR) enables us to explore the evolution of galaxies and intergalactic medium in the early universe. However, the detection and imaging of the 21\,cm-line signal are tough due to the foreground and instrumental systematics. In order to overcome these obstacles, as a new approach, we propose to take a cross correlation between observed 21\,cm-line data and 21\,cm-line images generated from the distribution of the Lyman-$\alpha$ emitters (LAEs) through machine learning. In order to create 21\,cm-line maps from LAE distribution, we apply conditional Generative Adversarial Network (cGAN) trained with the results of our numerical simulations. We find that the 21\,cm-line brightness temperature maps and the neutral fraction maps can be reproduced with correlation function of 0.5 at large scales $k<0.1~{\rm Mpc}^{-1}$. Furthermore, we study the detectability of the the cross correlation assuming the the LAE deep survey of the Subaru Hyper Suprime Cam, the 21\,cm observation of the MWA Phase II and the presence of the foreground residuals. We show that the signal is detectable at $k < 0.1~{\rm Mpc}^{-1}$ with 1000 hours of MWA observation even if the foreground residuals are 5 times larger than the 21\,cm-line power spectrum. Our new approach of cross correlation with image construction using the cGAN can not only boost the detectability of EoR 21\,cm-line signal but also allow us to estimate the 21\,cm-line auto-power spectrum., Comment: 16 pages, 12 figures, 4 tables, revised version

Ultracompact minihalos~(UCMHs) can form after the epoch of matter-radiation equality, if the density fluctuations of dark matter have significantly large amplitude on small scales. The constraint on the UCMH abundance allows us to access such small-scale fluctuations. In this paper, we present that, through the measurement of 21-cm fluctuations before the epoch of reionization~ we can obtain a constraint on the UCMH abundance. We calculate the 21-cm signal from UCMHs and show that UCMHs provide the enhancement of the 21-cm fluctuations. We also investigate the constraint on the UCMH abundance and small-scale curvature perturbations. Our results indicate that the upcoming 21-cm observation, the Square Kilometre Array (SKA), provides the constraint on amplitude of primordial curvature power spectrum, ${\cal A}_{\zeta} \lesssim 10^{-6}$ on $100~{\rm Mpc}^{-1} \lesssim k \lesssim 1000~{\rm Mpc}^{-1}$. Although it is not stronger than the one from the non-detection of gamma rays induced by dark matter annihilation in UCMHs, the constraint by the SKA will be important because this constraint is independent of the dark matter particle model., Comment: 8pages, 6 figures, published in MNRAS

Detecting $\rm H_I$ 21cm line in the intergalactic medium (IGM) during the Epoch of Reionization (EoR) suffers from foreground contamination such as Galactic synchrotron and extragalactic radio sources. Cross-correlation between the 21cm line and Lyman-$\alpha$ emitter (LAE) galaxies is a powerful tool to identify the 21cm signal since the 21cm line emission has correlation with LAEs while the LAEs are statistically independent of the foregrounds. So far, the detectability of 21cm-LAE cross-power spectrum has been investigated with simple LAE models where the observed Ly$\alpha$ luminosity is proportional to the dark matter halo mass. However, the previous models were inconsistent with the latest observational data of LAEs obtained with Subaru/Hyper Suprime-Cam (HSC). Here, we revisit the detectability of 21cm-LAE cross-power spectrum adopting a state-of-the-art LAE model consistent with all Subaru/HSC observations such as the Ly$\alpha$ luminosity function, LAE angular auto-correlation, and the LAE fractions in the continuum selected galaxies. We find that resultant cross-power spectrum with the updated LAE model is reduced at small scales ($k\sim 1\ \rm Mpc^{-1}$) compared to the simple models, while the amplitudes at large scales ($k \lesssim 0.2 \ \rm Mpc^{-1}$) are not affected so much. We conclude that the large-scale signal would be detectable with Square Kilometre Array (SKA) and HSC LAE cross-correlation but detecting the small scale signal would require an extended HSC LAE survey with an area of $\sim 75\ \rm deg^2$ or 3000 hrs observation time of 21cm line with SKA., Comment: 11 pages, 10 figures, accepted to Monthly Notices of the Royal Astronomical Society

We investigate the distribution of metals in the cosmological volume at $z\sim3$, in particular, provided by massive population III (Pop III) stars using a cosmological $N$-body simulation in which a model of Pop III star formation is implemented. Owing to the simulation, we can choose minihaloes where Pop III star formation occurs at $z>10$ and obtain the spatial distribution of the metals at lower-redshifts. To evaluate the amount of heavy elements provided by Pop III stars, we consider metal yield of pair-instability or core-collapse supernovae (SNe) explosions of massive stars. By comparing our results to the Illustris-1 simulation, we find that heavy elements provided by Pop III stars often dominate those from galaxies in low density regions. The median value of the volume averaged metallicity is $Z\sim 10^{-4.5 - -2} Z_{\odot}$ at the regions. Spectroscopic observations with the next generation telescopes are expected to detect the metals imprinted on quasar spectra., Comment: 10 pages, 8 figures, accepted for publication in MNRAS

In the present universe, magnetic fields exist with various strengths and on various scales. One possible origin of these cosmic magnetic fields is the primordial magnetic fields (PMFs) generated in the early universe. PMFs are considered to contribute to matter density evolution via Lorentz force and the thermal history of intergalactic medium (IGM) gas due to ambipolar diffusion. Therefore, information about PMFs should be included in the temperature anisotropy of the Cosmic Microwave Background through the thermal Sunyaev-Zel'dovich (tSZ) effect in IGM. In this article, given an initial power spectrum of PMFs, we show the spatial fluctuation of mass density and temperature of the IGM and tSZ angular power spectrum created by the PMFs. Finally, we find that the tSZ angular power spectrum induced by PMFs becomes significant on small scales, even with PMFs below the observational upper limit. Therefore, we conclude that the measurement of tSZ anisotropy on small scales will provide the most stringent constraint on PMFs., Comment: 8 pages, 2 figures, Conference report of "The Power of Faraday Tomography---Towards 3D Mapping of Cosmic Magnetic Fields---", published in Galaxies, special issue "The Power of Faraday Tomography"

We measure the scale dependence and redshift dependence of 21 cm line emitted from the neutral hydrogen gas at redshift 13 on scales of k>1h/Mpc, but it is roughly constant at lower redshift z<3. The redshift evolution of HI bias is relatively slow compared to that of QSOs at similar redshift range. We also measure a redshift space distortion (RSD) of HI gas to explore the properties of HI clustering. Fitting to a widely applied theoretical prediction, we find that the constant bias is consistent with that measured directly from the real-space power spectra, and the velocity dispersion is marginally consistent with the linear perturbation prediction. Finally we compare the results obtained from our simulation and the Illustris simulation, and conclude that the detailed astrophysical effects do not affect the scale dependence of HI bias very much, which implies that the cosmological analysis using 21 cm line of HI will be robust against the uncertainties arising from small-scale astrophysical processes such as star formation and supernova feedback., Comment: 10 pages, 6 figures, submitted to MNRAS

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., Comment: 13 pages, 11 figures, submitted to MNRAS

We present a novel Lyman alpha (Ly$\alpha$) radiative transfer code, SEURAT, where line scatterings are solved adaptively with the resolution of the smoothed particle hydrodynamics (SPH). The radiative transfer method implemented in SEURAT is based on a Monte Carlo algorithm in which the scattering and absorption by dust are also incorporated. We perform standard test calculations to verify the validity of the code; (i) emergent spectra from a static uniform sphere, (ii) emergent spectra from an expanding uniform sphere, and (iii) escape fraction from a dusty slab. Thereby we demonstrate that our code solves the Ly$\alpha$ radiative transfer with sufficient accuracy. We emphasise that SEURAT can treat the transfer of Ly$\alpha$ photons even in highly complex systems that have significantly inhomogeneous density fields. The high adaptivity of SEURAT is desirable to solve the propagation of Ly$\alpha$ photons in the interstellar medium of young star-forming galaxies like Ly$\alpha$ emitters (LAEs). Thus, SEURAT provides a powerful tool to model the emergent spectra of Ly$\alpha$ emission, which can be compared to the observations of LAEs., Comment: 11 pages, 10 figures, accepted for publication in MNRAS

We report fourteen and twenty-eight protocluster candidates at z=5.7 and 6.6 over 14 and 19 deg^2 areas, respectively, selected from 2,230 (259) Lya emitters (LAEs) photometrically (spectroscopically) identified with Subaru/Hyper Suprime-Cam (HSC) deep images (Keck, Subaru, and Magellan spectra and the literature data). Six out of the 42 protocluster candidates include 1-12 spectroscopically confirmed LAEs at redshifts up to z=6.574. By the comparisons with the cosmological Lya radiative transfer (RT) model reproducing LAEs with the reionization effects, we find that more than a half of these protocluster candidates are progenitors of the present-day clusters with a mass of > 10^14 M_sun. We then investigate the correlation between LAE overdensity delta and Lya rest-frame equivalent width EW_Lya^rest, because the cosmological Lya RT model suggests that a slope of EW_Lya^rest-delta relation is steepened towards the epoch of cosmic reionization (EoR), due to the existence of the ionized bubbles around galaxy overdensities easing the escape of Lya emission from the partly neutral intergalactic medium (IGM). The available HSC data suggest that the slope of the EW_Lya^rest-delta correlation does not evolve from the post-reionization epoch z=5.7 to the EoR z=6.6 beyond the moderately large statistical errors. There is a possibility that we would detect the evolution of the EW_Lya^rest - delta relation from z=5.7 to 7.3 by the upcoming HSC observations providing large samples of LAEs at z=6.6-7.3.

The survey of Lyman $\alpha$ emitters (LAEs) with Subaru Hyper Suprime-Cam, called SILVERRUSH (Ouchi et al.), is producing massive data of LAEs at $z\gtrsim6$. Here we present LAE simulations to compare the SILVERRUSH data. In 162$^3$ comoving Mpc$^3$ boxes, where numerical radiative transfer calculations of reionization were performed, LAEs have been modeled with physically motivated analytic recipes as a function of halo mass. We have examined $2^3$ models depending on the presence or absence of dispersion of halo Ly$\alpha$ emissivity, dispersion of the halo Ly$\alpha$ optical depth, $\tau_\alpha$, and halo mass dependence of $\tau_\alpha$. The unique free parameter in our model, a pivot value of $\tau_\alpha$, is calibrated so as to reproduce the $z=5.7$ Ly$\alpha$ luminosity function (LF). We compare our model predictions with Ly$\alpha$ LFs at $z=6.6$ and $7.3$, LAE angular auto-correlation functions (ACFs) at $z=5.7$ and $6.6$, and LAE fractions in Lyman break galaxies at $5

The Lyman continuum (LyC) escape fraction is a key quantity to determine the contribution of galaxies to cosmic reionization. It has been known that the escape fractions estimated by observations and numerical simulations show a large diversity. However, the origins of the diversity are still uncertain. In this work, to understand what quantities of galaxies are responsible for controlling the escape fraction, we numerically evaluate the escape fraction by performing ray-tracing calculation with simplified disc galaxy models. With a smooth disc model, we explore the dependence of the escape fraction on the disposition of ionizing sources, and find that the escape fraction varies up to $\sim 3$ orders of magnitude. It is also found that the halo mass dependence of disc scale height determines whether the escape fraction increases or decreases with halo mass. With a clumpy disc model, it turns out that the escape fraction increases as the clump mass fraction increases because the density in the inter-clump region decreases. In addition, we find that clumpiness regulates the escape fraction via two ways when the total clump mass dominates the total gas mass; the escape fraction is controlled by the covering factor of clumps if the clumps are dense sufficient to block LyC photons, otherwise the clumpiness works to reduce the escape fraction by increasing the total number of recombination events in a galaxy., Comment: 14 pages, 14 figures, 1 table

Detection of the redshifted 21cm-line signal from neutral hydrogen in the intergalactic medium (IGM) during the Epoch of Reionization (EoR) is complicated by intense foregrounds such as galactic synchrotron and extragalactic radio galaxies. The 21cm-Lyman-$\alpha$ emitter(LAE) cross-correlation is one of the tools available to reduce the foreground effects because the foreground emission from such radio sources is statistically independent of LAE distribution. LAE surveys during the EoR at redshifts $z=6.6$ and $7.3$ are ongoing by the Subaru Hyper Suprime-Cam (HSC). Additionally, Prime Focus Spectrograph (PFS) will provide precise redshift information of the LAEs discovered by the HSC survey. In this paper, we investigate the detectability of the 21cm signal with the 21cm-LAE cross-correlation by using our improved reionization simulations. We also focus on the error budget and evaluate it quantitatively in order to consider a strategy to improve the signal-to-noise ratio. In addition, we explore an expansion of the LAE survey to suggest optimal survey parameters and show a potential to measure a characteristic size of ionized bubbles via the turnover scale of the cross-power spectrum. As a result, we find that the Murchison Widefield Array (MWA) has ability to detect the cross-power spectrum signal on large scales by combining LAE Deep field survey of HSC. We also show that the sensitivity is improved dramatically at small scales by adding redshift information from the PFS measurements. The Square Kilometre Array (SKA) has a potential to measure the turnover scale with an accuracy of $6\times10^{-3}~{\rm Mpc^{-1}}$., Comment: 13pages, 13figures, 3tables, accepted to MNRAS

The presence of ubiquitous magnetic fields in the universe is suggested from observations of radiation and cosmic ray from galaxies or the intergalactic medium (IGM). One possible origin of cosmic magnetic fields is the magnetogenesis in the primordial universe. Such magnetic fields are called primordial magnetic fields (PMFs), and are considered to affect the evolution of matter density fluctuations and the thermal history of the IGM gas. Hence the information of PMFs is expected to be imprinted on the anisotropies of the cosmic microwave background (CMB) through the thermal Sunyaev-Zel'dovich (tSZ) effect in the IGM. In this study, given an initial power spectrum of PMFs as $P(k)\propto B_{\rm 1Mpc}^2 k^{n_{B}}$, we calculate dynamical and thermal evolutions of the IGM under the influence of PMFs, and compute the resultant angular power spectrum of the Compton $y$-parameter on the sky. As a result, we find that two physical processes driven by PMFs dominantly determine the power spectrum of the Compton $y$-parameter; (i) the heating due to the ambipolar diffusion effectively works to increase the temperature and the ionization fraction, and (ii) the Lorentz force drastically enhances the density contrast just after the recombination epoch. These facts result in making the tSZ angular power spectrum induced by the PMFs more remarkable at $\ell >10^4$ than that by galaxy clusters even with $B_{\rm 1Mpc}=0.1$ nG and $n_{B}=-1.0$ because the contribution from galaxy clusters decreases with increasing $\ell$. The measurement of the tSZ angular power spectrum on high $\ell$ modes can provide the stringent constraint on PMFs., Comment: 10 pages, 3 figures, 1 table, published in PRD

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., Comment: 16 pages, 10 figures. Accepted for publication in PASJ special issue

A large fraction of baryons predicted from the standard cosmology has been missing observationally. Although previous numerical simulations have indicated that most of the missing baryons reside in large-scale filaments in the form of Warm Hot Intergalactic Medium (WHIM), it is generally very difficult to detect signatures from such a diffuse gas. In this work, we focus on the hyperfine transition of neutral hydrogen (HI) called 21-cm line as a tool to trace the WHIM. For the purpose, we first construct the map of the 21-cm signals by using the data provided by the state-of-the-art cosmological hydrodynamics simulation project, Illustris, in which detailed processes affecting the dynamical and thermal evolution of the WHIM are implemented. From the comparison with the constructed 21-cm signal map with the expected noise level of the Square Kilometre Array phase 1 mid-frequency instrument (SKA1-mid), we find that the 21-cm signals from the filamentary structures at redshifts z=0.5-3 are detectable with the SKA1-mid if we assume the angular resolution of \Delta\theta > 10 arcmin and the observing time of t_obs > 100 hours. However, it also turns out that the signals mainly come from galaxies residing in the filamentary structures and the contribution from the WHIM is too small to detect with the SKA1-mid. Our results suggest that about 10 times higher sensitivity than the SKA1-mid is possibly enough to detect the WHIM at z=0.5-3., Comment: 9 pages, 7 figures, 1 table, submitted to MNRAS

Understanding Ly$\alpha$ emitting galaxies (LAEs) can be a key to reveal cosmic reionization and galaxy formation in the early Universe. Based on halo merger trees and Ly$\alpha$ radiation transfer calculations, we model redshift evolution of LAEs and their observational properties at $z \ge 6$. We consider ionized bubbles associated with individual LAEs and IGM transmission of Ly$\alpha$ photons. We find that Ly$\alpha$ luminosity tightly correlates with halo mass and stellar mass, while the relation with star formation rate has a large dispersion. Comparing our models with the observed luminosity function by Konno et al., we suggest that LAEs at $z \sim 7$ have galactic wind of $V_{\rm out} \gtrsim 100~\rm km\, s^{-1}$ and HI column density of $N_{\rm HI} \gtrsim 10^{20}~\rm cm^{-2}$. Number density of bright LAEs rapidly decreases as redshift increases, due to both lower star formation rate and smaller HII bubbles. Our model predicts future wide deep surveys with next generation telescopes, such as JWST, E-ELT and TMT, can detect LAEs at $z \sim 10$ with a number density of $n_{\rm LAE} \sim {\rm a~few~} \times10^{-6} ~\rm Mpc^{-3}$ for the flux sensitivity of $10^{-18} ~\rm erg\, cm^{-2}\, s^{-1}$. By combining these surveys with future 21-cm observations, it could be possible to detect both LAEs with $L_{Ly\alpha} \gtrsim 10^{42}~\rm erg~s^{-1}$ and their associated giant HII bubbles with the size $\gtrsim 250 ~\rm kpc$ at $z \sim 10$., Comment: 17 pages, 13 figures, submitted to MNRAS

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., Comment: 18 pages, 15 figures, 3 tables, accepted to ApJ; A typo in the title of the version 1 is corrected and replaced as the version 2

We explore the possibility of the formation of globular clusters under ultraviolet (UV) background radiation. One-dimensional spherical symmetric radiation hydrodynamics (RHD) simulations by Hasegawa et al. have demonstrated that the collapse of low-mass (10^6-10^7 solar masses) gas clouds exposed to intense UV radiation can lead to the formation of compact star clusters like globular clusters (GCs) if gas clouds contract with supersonic infall velocities. However, three-dimensional effects, such as the anisotropy of background radiation and the inhomogeneity in gas clouds, have not been studied so far. In this paper, we perform three-dimensional RHD simulations in a semi-cosmological context, and reconsider the formation of compact star clusters in strong UV radiation fields. As a result, we find that although anisotropic radiation fields bring an elongated shadow of neutral gas, almost spherical compact star clusters can be procreated from a "supersonic infall" cloud, since photo-dissociating radiation suppresses the formation of hydrogen molecules in the shadowed regions and the regions are compressed by UV heated ambient gas. The properties of resultant star clusters match those of GCs. On the other hand, in weak UV radiation fields, dark matter-dominated star clusters with low stellar density form due to the self-shielding effect as well as the positive feedback by ionizing photons. Thus, we conclude that the "supersonic infall" under a strong UV background is a potential mechanism to form GCs., Comment: 16 pages, 12 figures, 3 tables, accepted for publication in MNRAS

Cosmic reionization is known to be a major phase transition of the gas in the Universe. Since astronomical objects formed in the early Universe, such as the first stars, galaxies and black holes, are expected to have caused cosmic reionization, the formation history and properties of such objects are closely related to the reionization process. In spite of the importance of exploring reionization, our understandings regarding reionization is not sufficient yet. Square Kilometre Array (SKA) is a next-generation large telescope that will be operated in the next decade. Although several programs of next-generation telescopes are currently scheduled, the SKA will be the unique telescope with a potential to directly observe neutral hydrogen up to z~30, and provide us with valuable information on the Cosmic Dawn (CD) and the Epoch of Reionization (EoR). The early science with the SKA will start in a few years; it is thus the time for us to elaborate a strategy for CD/EoR Science with the SKA. The purpose of this document is to introduce Japanese scientific interests in the SKA project and to report results of our investigation., Comment: 14 pages, 12 figures, Japan SKA consortium (SKA-JP) Cosmic Dawn/Epoch of Reionization Science book

Understanding the detailed process of cosmic reionization is one of the remaining problems in astrophysics and cosmology. Here we construct a model of cosmic reionization that includes contribution from high-$z$ galaxies and active galactic nuclei (AGNs), and calculate the reionization and thermal histories with the model. To keep the model to be general and realistic, we vary the escape fraction of ionizing photons, $f_{\rm esc}$, and the faint-end slope of the AGN luminosity function at high redshifts, $\alpha_{\rm hz}$, within the constraints from the observed cosmic star formation history and the observed bright-end UV luminosity functions at $z \leq 6$. Additionally, we model the spectral energy distribution (SED) of AGNs which depends on the Eddington ratio and the black-hole mass. By comparing the computed reionization histories with the observed \ion{H}{i} fractions and the optical depth for Thomson scattering from Planck, we find that $\alpha_{\rm hz} > -1.5$ and $f_{\rm esc} < 0.15$ are favored when we employ the bright-end luminosity function obtained by Giallongo et al. Our result suggests that an AGN-dominated model with the abundance of faint-AGNs being as large as the estimate by Giallongo et al. is allowed only if the contribution from high-z galaxies is almost negligible, while a galaxy-dominant model is also allowed. We also find that the shape of SED has a significant impact on the thermal history. Therefore it is expected that measurements of the thermal state of the IGM provide useful information on properties of ionizing sources., Comment: 14 pages, 14 figures, accepted to MNRAS

We study the number and the distribution of low mass Pop III stars in the Milky Way. In our numerical model, hierarchical formation of dark matter minihalos and Milky Way sized halos are followed by a high resolution cosmological simulation. We model the Pop III formation in H2 cooling minihalos without metal under UV radiation of the Lyman-Werner bands. Assuming a Kroupa IMF from 0.15 to 1.0 Msun for low mass Pop III stars, as a working hypothesis, we try to constrain the theoretical models in reverse by current and future observations. We find that the survivors tend to concentrate on the center of halo and subhalos. We also evaluate the observability of Pop III survivors in the Milky Way and dwarf galaxies, and constraints on the number of Pop III survivors per minihalo. The higher latitude fields require lower sample sizes because of the high number density of stars in the galactic disk, the required sample sizes are comparable in the high and middle latitude fields by photometrically selecting low metallicity stars with optimized narrow band filters, and the required number of dwarf galaxies to find one Pop III survivor is less than ten at <100 kpc for the tip of redgiant stars. Provided that available observations have not detected any survivors, the formation models of low mass Pop III stars with more than ten stars per minihalo are already excluded. Furthermore, we discuss the way to constrain the IMF of Pop III star at a high mass range of > 10 Msun., Comment: accepted by ApJ, 11 pages, 8 figures

We present a new numerical scheme to solve the transfer of diffuse radiation on three-dimensional mesh grids which is efficient on processors with highly parallel architecture such as recently popular GPUs and CPUs with multi- and many-core architectures. The scheme is based on the ray-tracing method and the computational cost is proportional to $N_{\rm m}^{5/3}$ where $N_{\rm m}$ is the number of mesh grids, and is devised to compute the radiation transfer along each light-ray completely in parallel with appropriate grouping of the light-rays. We find that the performance of our scheme scales well with the number of adopted CPU cores and GPUs, and also that our scheme is nicely parallelized on a multi-node system by adopting the multiple wave front scheme, and the performance scales well with the amount of the computational resources. As numerical tests to validate our scheme and to give a physical criterion for the angular resolution of our ray-tracing scheme, we perform several numerical simulations of the photo-ionization of neutral hydrogen gas by ionizing radiation sources without the "on-the-spot" approximation, in which the transfer of diffuse radiation by radiative recombination is incorporated in a self-consistent manner., Comment: 26 pages, 15 figures, accepted for publication in PASJ

We perform cosmological hydrodynamics simulations with non-equilibrium primordial chemistry to obtain 59 minihalos that host first stars. The obtained minihalos are used as initial conditions of local three dimensional radiation hydrodynamics simulations to investigate the formation of the first stars. We find two-thirds of the minihalos host multiple stars, while the rest of them have single stars. The mass of the stars found in our simulations are in the range of 1 Msun \la M \la 300 Msun, peaking at several x 10 Msun. Most of the very massive stars of \ga 140 Msun are born as single stars, although not all of the single stars are very massive. We also find a few stars of \la 1 Msun that are kicked by the gravitational three body interactions to the position distant from the center of mass. The frequency that a star forming minihalo contains a binary system is \sim 50%. We also investigate the abundance pattern of the stellar remnants by summing up the contributions from the first stars in the simulations. Consequently, the pattern is compatible with that of the low metallicity Damped Lyman-alpha systems or the Extremely Metal Poor (EMP) stars, if the mass spectrum obtained in our experiment is shifted to the low mass side by 0.2 dex. If we consider the case that an EMP star is born in the remnant of the individual minihalo without mixing with others, the chemical signature of the pair instability supernova is more prominent, because most of them are born as single stars., Comment: 18 pages ApJ accepted

We perform three-dimensional radiation hydrodynamic simulations of uniform dusty gas clouds irradiated by an active galactic nucleus (AGN) to investigate the dependence of evolution of clouds on the ionization parameter $\mathcal{U}$ and the Str{\"o}mgren number $\mathcal{N}_{S}$. We find that the evolution can be classified into two cases depending on $\mathcal{U}$. In low $\mathcal{U}$ cases ($\mathcal{U}\approx 10^{-2}$), the evolution is mainly driven by photo-evaporation. A approximately spherically-symmetric evaporation flow with velocity of $100\operatorname{-}150\;\mathrm{km\;s^{-1}}$ is launched from the irradiated face. The cloud is compressed by a D-type shock with losing its mass due to photo-evaporation and is finally turned into a dense filament by $t\lesssim 1.5t_{\mathrm{sc}}$. In high $\mathcal{U}$ cases ($\mathcal{U}\approx 5\times 10^{-2}$), radiation pressure suppresses photo-evaporation from the central part of the irradiated face, reducing photo-evaporation rate. A evaporation flow from the outskirts of the irradiated face is turned into a high velocity ($\lesssim 500\;\mathrm{km\;s^{-1}}$) gas wind because of radiation pressure on dust. The cloud is swept by a radiation pressure-driven shock and becomes a dense gas disk by $t\approx t_{\mathrm{sweep}}$. Star formation is expected in these dense regions for both cases of $\mathcal{U}$. We discuss the influences of the AGN radiation on the clumpy torus. A simple estimate suggests that the clumps are destroyed in timescales shorter than their orbital periods. For the clumpy structure to be maintained over long period, the incident radiation field needs to be sufficiently weaken for most of the clumps, or, some mechanism that creates the clumps continuously is needed., Comment: 33 pages, 36 figures, 6 tables; accepted for publication in Monthly Notices of the Royal Astronomical Society

First, the formation of first objects driven by dark matter is revisited by high-resolution hydrodynamic simulations. It is revealed that dark matter haloes of ~10^4M_sun can produce first luminous objects with the aid of dark matter cusps. Therefore, the mass of first objects is smaller by roughly two orders of magnitude than in the previous prediction. This implies that the number of Pop III stars formed in the early universe could be significantly larger than hitherto thought. Secondly, the feedback by photo-ionization and photo-dissociation photons in the first objects is explored with radiation hydrodynamic simulations, and it is demonstrated that multiple stars can form in a 10^5M_sun halo. Thirdly, the fragmentation of an accretion disk around a primordial protostar is explored with photo-dissociation feedback. As a result, it is found that the photo-dissociation can reduce the mass accretion rate onto protostars. Also, protostars as small as 0.8M_sun may be ejected and evolve with keeping their mass, which might be detected as "real first stars" in the Galactic halo. Finally, state-of-the-art radiation hydrodynamic simulations are performed to investigate the internal ionization of first galaxies and the escape of ionizing photons. We find that UV feedback by forming massive stars enhances the escape fraction even in a halo as massive as > 6* 10^9M_sun, while it reduces the star formation rate significantly. This may have a momentous impact on the cosmic reionization., Comment: 26 pages, 14 figures

We explore the impacts of ultraviolet (UV) radiation feedback on galaxies during the epoch of reionisation by cosmological simulations in which hydrodynamics and the transfer of the H and He ionising photons are consistently coupled. Moreover we take into account H_2 non-equilibrium chemistry, including photo-dissociation. The most striking feature of the simulations is a high spatial resolution for the radiative transfer (RT) calculation which enables us to start considering not only external UV feedback processes but also internal UV feedback processes in each galaxy. We find that the star formation is significantly suppressed due to the internal UV and supernova (SN) feedback. In low mass galaxies with M<10^9Msun, a large amount of gas is evacuated by photo-evaporation as previous studies have shown, which results in the suppression of star formation. Surprisingly, star formation in massive halos is also strongly suppressed despite the fact that these halos hardly lose any gas by photo-evaporation. The suppression of star formation in massive halos is mainly caused by following two factors; (i) small scale clumpy structures in the galaxies are smoothened by the internal feedback, (ii) although the dense gas in the galaxies is mostly neutral, the H_2 formation and cooling processes are disturbed by mild photo-heating. Photo-dissociating radiations actually suppress star formation, but the magnitude of the effect is not so large in massive galaxies. Even though our simulation volume is too small to be a representative patch of the Universe during reionisation, we find that our simulated star formation rate densities and HI fractions at z=6-7 are consistent with those found in observations., Comment: 14 pages, 13 figures, 1 table, accepted for publication in MNRAS

We present a novel radiation hydrodynamics code, START, which is a smoothed particle hydrodynamics (SPH) scheme coupled with accelerated radiative transfer. The basic idea for the acceleration of radiative transfer is parallel to the tree algorithm that is hitherto used to speed up the gravitational force calculation in an N-body system. It is demonstrated that the radiative transfer calculations can be dramatically accelerated, where the computational time is scaled as Np log Ns for Np SPH particles and Ns radiation sources. Such acceleration allows us to readily include not only numerous sources but also scattering photons, even if the total number of radiation sources is comparable to that of SPH particles. Here, a test simulation is presented for a multiple source problem, where the results with START are compared to those with a radiation SPH code without tree-based acceleration. We find that the results agree well with each other if we set the tolerance parameter as < 1.0, and then it demonstrates that START can solve radiative transfer faster without reducing the accuracy. One of important applications with START is to solve the transfer of diffuse ionizing photons, where each SPH particle is regarded as an emitter. To illustrate the competence of START, we simulate the shadowing effect by dense clumps around an ionizing source. As a result, it is found that the erosion of shadows by diffuse recombination photons can be solved. Such an effect is of great significance to reveal the cosmic reionization process., Comment: 14 pages, 23 figures, accepted for publication in MNRAS

We explore the formation of secondary Population III (Pop III) stars under radiation hydrodynamic (RHD) feedback by a preformed massive star. To properly treat RHD feedback, we perform three-dimensional RHD simulations incorporating the radiative transfer of ionizing photons as well as H_2 dissociating photons from a preformed star. A collapsing gas cloud is settled at a given distance from a 120Msun Pop III star, and the evolution of the cloud is pursued including RHD feedback. We derive the threshold density depending on the distance, above which the cloud can keep collapsing owing to the shielding of H_2 dissociating radiation. We find that an H_2 shell formed ahead of an ionizing front works effectively to shield the H_2 dissociating radiation, leading to the positive feedback for the secondary Pop III star formation. Also, near the threshold density, the envelope of gas cloud is stripped significantly by a shock associated with an ionizing front. By comparing the mass accretion timescale with the Kelvin-Helmholtz timescale, we estimate the mass of secondary Pop III stars. It turns out that the stripping by a shock can reduce the mass of secondary Pop III stars down to \approx 20Msum., Comment: ApJ accepted, 11 pages, 11 figures

We present a novel scenario for globular cluster (GC) formation, where the ultraviolet (UV) background radiation effectively works so as to produce compact star clusters. Here, we explore the formation of GCs in UV radiation fields. For this purpose, we calculate baryon and dark matter (DM) dynamics in spherical symmetry, incorporating the self-shielding effects by solving the radiative transfer of UV radiation. In addition, we prescribe the star formation in cooled gas components and pursue the dynamics of formed stars. As a result, we find that the evolution of subgalactic objects in UV background radiation are separated into three types, that is, (1) prompt star formation, where less massive clouds ~10^{5-8} M_sun are promptly self-shielded and undergo star formation, (2) delayed star formation, where photoionized massive clouds >10^8 M_sun collapse despite high thermal pressure and are eventually self-shielded to form stars in a delayed fashion, and (3) supersonic infall, where photoionized less massive clouds ~10^{5-8} M_sun contract with supersonic infall velocity and are self-shielded when a compact core forms. In particular, the type (3) is a novel type found in the present simulations, and eventually produces a very compact star cluster. The resultant mass-to-light ratios, half-mass radii, and velocity dispersions for the three types are compared to the observations of GCs, dwarf spheroidals (dSphs), and ultra-compact dwarfs (UCDs). It turns out that the properties of star clusters resulting from supersonic infall match well with those of observed GCs, whereas the other two types are distinct from GCs. Hence, we conclude that supersonic infall in a UV background is a promising mechanism to form GCs., Comment: 11pages, 9figures, accepted for publication in MNRAS, a few typos are fixed

We explore the impact of ultraviolet (UV) radiation from massive Population III (Pop III) stars of 25, 40, 80, and 120 M_sun on the subsequent Pop III star formation. In this paper, particular attention is paid to the dependence of radiative feedback on the mass of source Pop III star. UV radiation from the source star can work to impede the secondary star formation through the photoheating and photodissociation processes. Recently, Susa & Umemura (2006) have shown that the ionizing radiation alleviates the negative effect by H_2-dissociating radiation from 120$M_sun PopIII star, since an H_2 shell formed ahead of an ionizing front can effectively shield H_2-dissociating radiation. On the other hand, it is expected that the negative feedback by H_2-dissociating radiation can be predominant if a source star is less massive, since a ratio of the H_2-dissociating photon number to the ionizing photon number becomes higher. In order to investigate the radiative feedback effects from such less massive stars, we perform three-dimensional radiation hydrodynamic simulations, incorporating the radiative transfer effect of ionizing and H_2-dissociating radiation. As a result, we find that if a source star is less massive than ~25M_sun, the ionizing radiation cannot suppress the negative feedback of H_2-dissociating radiation. Therefore, the fate of the neighboring clouds around such less massive stars is determined solely by the flux of H_2-dissociating radiation from source stars. With making analytic estimates of H_2 shell formation and its shielding effect, we derive the criteria for radiation hydrodynamic feedback depending on the source star mass., Comment: 8 pages, 9 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

We study (not necessarily connected) Z-graded A-infinity-algebras and their A-infinity-modules. Using the cobar and the bar construction and Quillen's homotopical algebra, we describe the localisation of the category of A-infinity-algebras with respect to A-infintity-quasi-isomorphisms. We then adapt these methods to describe the derived category of an augmented A-infinity-algebra A. The case where A is not endowed with an augmentation is treated differently. Nevertheless, when A is strictly unital, its derived category can be described in the same way as in the augmented case. Next, we compare two different notions of A-infinity-unitarity : strict unitarity and homological unitarity. We show that, up to homotopy, there is no difference between these two notions. We then establish a formalism which allows us to view A-infini-categories as A-infinity-algebras in suitable monoidal categories. We generalize the fundamental constructions of category theory to this setting : Yoneda embeddings, categories of functors, equivalences of categories... We show that any algebraic triangulated category T which admits a set of generators is A-infinity-pretriangulated, that is to say, T is equivalent to $H^0 tw A$, where $tw A$ is the A-infinity-category of twisted objets of a certain A-infinity-category A., Comment: 230 pages, 11 figures, ph. D. thesis in french