Source author record

Naoki Yoshida

Naoki Yoshida appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA astro-ph.HE astro-ph.SR astro-ph gr-qc hep-ph astro-ph.IM astro-ph.EP physics.optics Cryptography and Security physics.comp-ph physics.data-an

Catalog footprint

What is connected

109works

13topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2024arXiv

A merger-driven scenario for clumpy galaxy formation in the epoch of reionization: Physical properties of clumps in the FirstLight simulation

Recent JWST observations with superb angular resolution have revealed the existence of clumpy galaxies at high redshift through the detection of rest-frame optical emission lines. We use the FirstLight simulation to study the properties of (sub-)galactic clumps that are bright in [OIII] 5007$\mathrm{\mathring{A}}$ line with flux greater than $\sim 10^{-18} \, {\rm erg\, s^{-1}\, cm^{-2}}$, to be detected by JWST. For 62 simulated galaxies that have stellar masses of $(0.5-6) \times 10^{10} \, M_\odot$ at $z=5$, we find clumps in 1828 snapshots in the redshift range $z = 9.5-5.5$. The clumps are identified by the surface density of star formation rate. About one-tenth of the snapshots show the existence of clumpy systems with two or more components. Most of the clumps are formed by mergers and can be characterized by their ages; central clumps dominated by stellar populations older than 50 Myr, and off-centered clumps dominated by younger stellar populations with specific star formation rates of $\sim 50 \, {\rm Gyr^{-1}}$. The latter type of young clumps is formed from gas debris in the tidal tails of major mergers with baryonic mass ratios of $1 \leq q < 4$. The merger-induced clumps are short-lived, and merge within a dynamical time of several tens million years. The number density of the clumpy systems is estimated to be $\sim 10^{-5}\, {\rm cMpc^{-3}}$, which is large enough to be detected in recent JWST surveys.

preprint2022arXiv

$\mathrm{H_2}$ cooling and gravitational collapse of supersonically induced gas objects

We study the formation and gravitational collapse of supersonically induced gas objects (SIGOs) in the early universe. We run cosmological hydrodynamics simulations of SIGOs, including relative streaming motions between baryons and dark matter. Our simulations also follow nonequilibrium chemistry and molecular hydrogen cooling in primordial gas clouds. A number of SIGOs are formed in the run with fast-streaming motions of 2 times the rms of the cosmological velocity fluctuations. We identify a particular gas cloud that condensates by H$_2$ cooling without being hosted by a dark matter halo. The SIGO remains outside the virial radius of its closest halo, and it becomes Jeans unstable when the central gas-particle density reaches $\sim 100~{\rm cm}^{-3}$ with a temperature of $\sim$ 200 K. The corresponding Jeans mass is $\sim 10^5 M_{\odot}$, and thus the formation of primordial stars or a star cluster is expected in the SIGO.

preprint2022arXiv

Deep-learning Real/Bogus classification for the Tomo-e Gozen transient survey

We present a deep neural network Real/Bogus classifier that improves classification performance in the Tomo-e Gozen transient survey by handling label errors in the training data. In the wide-field, high-frequency transient survey with Tomo-e Gozen, the performance of conventional convolutional neural network classifier is not sufficient as about $10^6$ bogus detections appear every night. In need of a better classifier, we have developed a new two-stage training method. In this training method, label errors in the training data are first detected by normal supervised learning classification, and then they are unlabeled and used for training of semi-supervised learning. For actual observed data, the classifier with this method achieves an area under the curve (AUC) of 0.9998 and a false positive rate (FPR) of 0.0002 at true positive rate (TPR) of 0.9. This training method saves relabeling effort by humans and works better on training data with a high fraction of label errors. By implementing the developed classifier in the Tomo-e Gozen pipeline, the number of transient candidates was reduced to $\sim$40 objects per night, which is $\sim$1/130 of the previous version, while maintaining the recovery rate of real transients. This enables more efficient selection of targets for follow-up observations.

preprint2022arXiv

Detecting preheating in proto-clusters with Lyman-$α$ Forest Tomography

Studies of low redshift galaxy clusters suggest the intra-cluster medium (ICM) has experienced non-gravitational heating during the formation phase of the clusters. Using simple phenomenological heating prescriptions, we simulate the effect of this preheating of the nascent ICM in galaxy proto-clusters and examine its effect on Lyman-$α$ (Ly$α$) forest tomographic maps. We analyse a series of cosmological zoom-in simulations of proto-clusters within the framework of the Ly$α$ transmission-dark matter (DM) density distribution. We find that the more energy is injected into the proto-ICM at $z$ = 3, the more the distribution at high DM density tilts towards higher Ly$α$ transmission. This effect has been confirmed in both low-resolution simulations adopting a preheating scheme based on entropy floors, as well as in higher-resolution simulations with another scheme based on energy floors. The evolution of the slope of this distribution is shown to vary with redshift. The methodology developed here can be applied to current and upcoming Ly$α$ forest tomographic survey data to help constrain feedback models in galaxy proto-clusters.

preprint2022arXiv

Possible Systematic Rotation in the Mature Stellar Population of a $z=9.1$ Galaxy

We present new observations with the Atacama Large Millimeter/submillimeter Array for a gravitationally-lensed galaxy at $z=9.1$, MACS1149-JD1. [O III] 88-$μ$m emission is detected at 10$σ$ with a spatial resolution of $\sim0.3$ kpc in the source plane, enabling the most distant morpho-kinematic study of a galaxy. The [O III] emission is distributed smoothly without any resolved clumps and shows a clear velocity gradient with $ΔV_{\rm obs}/2σ_{\rm tot}=0.84\pm0.23$, where $ΔV_{\rm obs}$ is the observed maximum velocity difference and $σ_{\rm tot}$ is the velocity dispersion measured in the spatially-integrated line profile, suggesting a rotating system. Assuming a geometrically thin self-gravitating rotation disk model, we obtain $V_{\rm rot}/σ_V=0.67_{-0.26}^{+0.73}$, where $V_{\rm rot}$ and $σ_V$ are the rotation velocity and velocity dispersion, respectively, still consistent with rotation. The resulting disk mass of $0.65_{-0.40}^{+1.37}\times10^{9}$ M$_\odot$ is consistent with being associated with the stellar mass identified with a 300 Myr-old stellar population independently indicated by a Balmer break in the spectral energy distribution. We conclude that the most of the dynamical mass is associated with the previously-identified mature stellar population that formed at $z\sim15$.

preprint2022arXiv

Public Release of A-SLOTH: Ancient Stars and Local Observables by Tracing Halos

The semi-analytical model A-SLOTH (Ancient Stars and Local Observables by Tracing Halos) is the first public code that connects the formation of the first stars and galaxies to observables. After several successful projects with this model, we publish the source code and describe the public version in this paper. The model is based on dark matter merger trees that can either be generated based on Extended Press-Schechter theory or that can be imported from dark matter simulations. On top of these merger trees, A-SLOTH applies analytical recipes for baryonic physics to model the formation of both metal-free and metal-poor stars and the transition between them with unprecedented precision and fidelity. A-SLOTH samples individual stars and includes radiative, chemical, and mechanical feedback. It is calibrated based on six observables, such as the optical depth to Thomson scattering, the stellar mass of the Milky Way and its satellite galaxies, the number of extremely-metal poor stars, and the cosmic star formation rate density at high redshift. A-SLOTH has versatile applications with moderate computational requirements. It can be used to constrain the properties of the first stars and high-z galaxies based on local observables, predicts properties of the oldest and most metal-poor stars in the Milky Way, can serve as a subgrid model for larger cosmological simulations, and predicts next-generation observables of the early Universe, such as supernova rates or gravitational wave events.

preprint2022arXiv

Stellar Wind Effect on the Atmospheric Escape of Hot Jupiters and their Ly-$α$ and H$α$ transits

Atmospheric escape of close-in exoplanets can be driven by high energy radiation from the host star. The planetary outflows interacting with the stellar wind may generate observable transit signals that depend on the strength of the stellar wind. We perform detailed radiation-hydrodynamics simulations of the atmospheric escape of hot Jupiters with including the wind from the host star in a self-consistent, dynamically coupled manner. We show that the planetary outflow is shaped by the balance between its thermal pressure and the ram pressure of the stellar wind. We use the simulation outputs to calculate the Lyman-$α$ and H$α$ transit signatures. Strong winds can confine the outflow and decrease the Lyman-$α$ transit depth. Contrastingly, the wind effect on H$α$ is weak because of the small contribution from the uppermost atmosphere of the planet. Observing both of the lines is important to understand the effect of the UV radiation and wind from the host. The atmospheric mass-loss rate is approximately independent of the strength of the wind. We also discuss the effect of the coronal mass ejections on the signatures. We argue that around M dwarfs the effect can be significant in every transit.

preprint2022arXiv

The three-year shear catalog of the Subaru Hyper Suprime-Cam SSP Survey

We present the galaxy shear catalog that will be used for the three-year cosmological weak gravitational lensing analyses using data from the Wide layer of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) Survey. The galaxy shapes are measured from the $i$-band imaging data acquired from 2014 to 2019 and calibrated with image simulations that resemble the observing conditions of the survey based on training galaxy images from the Hubble Space Telescope in the COSMOS region. The catalog covers an area of 433.48 deg$^2$ of the northern sky, split into six fields. The mean $i$-band seeing is 0.59 arcsec. With conservative galaxy selection criteria (e.g., $i$-band magnitude brighter than 24.5), the observed raw galaxy number density is 22.9 arcmin$^{-2}$, and the effective galaxy number density is 19.9 arcmin$^{-2}$. The calibration removes the galaxy property-dependent shear estimation bias to a level: $|δm|<9\times 10^{-3}$. The bias residual $δm$ shows no dependence on redshift in the range $0<z\leq 3$. We define the requirements for cosmological weak lensing science for this shear catalog, and quantify potential systematics in the catalog using a series of internal null tests for systematics related to point-spread function modelling and shear estimation. A variety of the null tests are statistically consistent with zero or within requirements, but (i) there is evidence for PSF model shape residual correlations; and (ii) star-galaxy shape correlations reveal additive systematics. Both effects become significant on $>1$ degree scales and will require mitigation during the inference of cosmological parameters using cosmic shear measurements.

preprint2022arXiv

Three-Dimensional Reconstruction of Weak Lensing Mass Maps with a Sparsity Prior. I. Cluster Detection

We propose a novel method to reconstruct high-resolution three-dimensional mass maps using data from photometric weak-lensing surveys. We apply an adaptive LASSO algorithm to perform a sparsity-based reconstruction on the assumption that the underlying cosmic density field is represented by a sum of Navarro-Frenk-White halos. We generate realistic mock galaxy shape catalogues by considering the shear distortions from isolated halos for the configurations matched to Subaru Hyper Suprime-Cam Survey with its photometric redshift estimates. We show that the adaptive method significantly reduces line-of-sight smearing that is caused by the correlation between the lensing kernels at different redshifts. Lensing clusters with lower mass limits of $10^{14.0} h^{-1}M_{\odot}$, $10^{14.7} h^{-1}M_{\odot}$, $10^{15.0} h^{-1}M_{\odot}$ can be detected with 1.5-$σ$ confidence at the low ($z<0.3$), median ($0.3\leq z< 0.6$) and high ($0.6\leq z< 0.85$) redshifts, respectively, with an average false detection rate of 0.022 deg$^{-2}$. The estimated redshifts of the detected clusters are systematically lower than the true values by $Δz \sim 0.03$ for halos at $z\leq 0.4$, but the relative redshift bias is below $0.5\%$ for clusters at $0.4<z\leq 0.85$. The standard deviation of the redshift estimation is $0.092$. Our method enables direct three-dimensional cluster detection with accurate redshift estimates.

preprint2021arXiv

Atmospheric Escape of Close-in Giants around Hot Stars: Far-Ultraviolet Radiation and Photoelectric Heating Effect

Atmospheric escape is an important process that controls the long-term evolution of close-in planets. We perform radiation hydrodynamics simulations of photo-evaporation of exoplanets' atmospheres to study the effect of photoelectric heating by far-ultraviolet (FUV) radiation. Specifically, we consider a close-in hot Jupiter around a hot A-star. Hot main-sequence stars emit not only extreme ultraviolet radiation but also FUV radiation, and thus can drive strong atmospheric escape by photoelectric heating. We show that the planetary atmosphere escapes at a rate as large as $\dot{M}\sim10^{14}\, \mathrm{g}~{\rm sec}^{-1}$ if the atmosphere contains a small amount of dust grains with the level of ten percent of the local interstellar medium. Close-in planets around hot stars can lose a significant fraction of the atmosphere during the long-term evolution. We quantify the amount of dust necessary for causing photoevaporation. The dust-to-gas mass ratio of $10^{-4}$ is sufficient to drive stronger atmospheric escape by FUV photoelectric heating than in the case with only extreme ultraviolet radiation. We also explore the metallicity dependence of the FUV-driven escape. The mass-loss rate increases with increasing the atmosphere's metallicity because of the enhanced photoelectric heating, but the stellar FUV flux decreases with increasing stellar metallicity. We derive an accurate estimate for the mass-loss rate as a function of FUV flux and metallicity, and of the planet's characteristics. The FUV driven atmospheric escape may be a key process to understand and explain the so-called sub-Jovian desert.

preprint2021arXiv

Capturing the inside-out quenching by black holes with far-infrared atomic line ratios

We propose to use relative strengths of far-infrared fine structure lines from galaxies to characterise early phases of the inside-out quenching by massive black holes (BHs). The BH feedback is thought to quench star formation by evacuating the ambient gas. In order to quantify the feedback effect on the gas density in the galactic centres, we utilise the outputs of IllustrisTNG and Illustris simulations, which implement different BH feedback models. We devise a physical model of H$_{\rm ~II}$ regions and compute the intensities of [O$_{\rm ~III}$] $52$ and $88~{\rm μm}$ lines. The line intensity ratio is sensitive to the local electron density, and thus can be used to measure the strength and physical extent of the BH quenching. If the BH feedback abruptly operates and expel the gas when it grows to a certain mass, as modelled in IllustrisTNG, the low-density gas yields relatively weak [O$_{\rm ~III}$] $52$ line with respect to $88~{\rm μm}$. In contrast, if the feedback strength and hence the local gas density are not strongly correlated with the BH mass, as in Illustris, the line ratio is not expected to vary significantly among galaxies with different evolutionary stages. We find these features are reproduced in the simulations. We also show that the line ratios are not sensitive to the aperture size for measurement, and thus observations do not need to resolve the galactic centres. We argue that the integrated line ratios can be used to capture the onset of the inside-out quenching by BHs.

preprint2021arXiv

Large-scale variation in reionization history caused by Baryon-dark matter streaming velocity

At cosmic recombination, there was supersonic relative motion between baryons and dark matter, which originated from the baryonic acoustic oscillations in the early universe. This motion has been considered to have a negligible impact on the late stage of cosmic reionization because the relative velocity quickly decreases. However, recent studies have suggested that the recombination in gas clouds smaller than the local Jeans mass ($\lesssim$ $10^8~M_\odot$) can affect the reionization history by boosting the number of ultraviolet photons required for ionizing the intergalactic medium. Motivated by this, we performed a series of radiation-hydrodynamic simulations to investigate whether the streaming motion can generate variation in the local reionization history by smoothing out clumpy small-scale structures and lowering the ionizing photon budget. We found that the streaming velocity can add a variation of $Δz_e$ $\sim$ $0.05$ $-$ $0.5$ in the end-of-reionization redshift, depending on the level of X-ray preheating and the time evolution of ionizing sources. The variation tends to be larger when the ionizing efficiency of galaxies decreases toward later times. Given the long spatial fluctuation scales of the streaming motion ($\gtrsim 100$ Mpc), it can help to explain the Ly$α$ opacity variation observed from quasars and leave large-scale imprints on the ionization field of the intergalactic medium during the reionization. The pre-reionization heating by X-ray sources is another critical factor that can suppress small-scale gas clumping and can diminish the variation in $z_e$ introduced by the streaming motion.

preprint2020arXiv

ALMA uncovers the [CII] emission and warm dust continuum in a z = 8.31 Lyman break galaxy

We report on the detection of the [CII] 157.7 $μ$m emission from the Lyman break galaxy (LBG) MACS0416_Y1 at z = 8.3113, by using the Atacama Large Millimeter/submillimeter Array (ALMA). The luminosity ratio of [OIII] 88 $μ$m (from previous campaigns) to [CII] is 9.31 $\pm$ 2.6, indicative of hard interstellar radiation fields and/or a low covering fraction of photo-dissociation regions. The emission of [CII] is cospatial to the 850 $μ$m dust emission (90 $μ$m rest-frame, from previous campaigns), however the peak [CII] emission does not agree with the peak [OIII] emission, suggesting that the lines originate from different conditions in the interstellar medium. We fail to detect continuum emission at 1.5 mm (160 $μ$m rest-frame) down to 18 $μ$Jy (3$σ$). This nondetection places a strong limit on the dust spectrum, considering the 137 $\pm$ 26 $μ$Jy continuum emission at 850 $μ$m. This suggests an unusually warm dust component (T $>$ 80 K, 90% confidence limit), and/or a steep dust-emissivity index ($β_{\rm dust}$ $>$ 2), compared to galaxy-wide dust emission found at lower redshifts (typically T $\sim$ 30 - 50 K, $β_{\rm dust}$ $\sim$ 1 - 2). If such temperatures are common, this would reduce the required dust mass and relax the dust production problem at the highest redshifts. We therefore warn against the use of only single-wavelength information to derive physical properties, recommend a more thorough examination of dust temperatures in the early Universe, and stress the need for instrumentation that probes the peak of warm dust in the Epoch of Reionization.

preprint2020arXiv

Cross Correlation of the Extragalactic Gamma-ray Background with Thermal Sunyaev-Zel'dovich Effect in the Cosmic Microwave Background

Cosmic rays in galaxy clusters are unique probes of energetic processes operating with large-scale structures in the Universe. Precise measurements of cosmic rays in galaxy clusters are essential for improving our understanding of non-thermal components in the intracluster-medium (ICM) as well as the accuracy of cluster mass estimates in cosmological analyses. In this paper, we perform a cross-correlation analysis with the extragalactic gamma-ray background and the thermal Sunyaev-Zel'dovich (tSZ) effect in the cosmic microwave background. The expected cross-correlation signal would contain rich information about the cosmic-ray-induced gamma-ray emission in the most massive galaxy clusters at $z\sim0.1-0.2$. We analyze the gamma-ray background map with 8 years of data taken by the Large Area Telescope onboard Fermi satellite and the publicly available tSZ map by Planck. We confirm that the measured cross-correlation is consistent with a null detection, and thus it enables us to put the tightest constraint on the acceleration efficiency of cosmic ray protons at shocks in and around galaxy clusters. We find the acceleration efficiency must be below 5\% with a $2σ$ confidence level when the hydrostatic mass bias of clusters is assumed to be 30\%, and our result is not significantly affected by the assumed value of the hydrostatic mass bias. Our constraint implies that the non-thermal cosmic-ray pressure in the ICM can introduce only a $\le 3\%$ level of the hydrostatic mass bias, highlighting that cosmic rays alone do not account for the mass bias inferred by the Planck analyses. Finally, we discuss future detectability prospects of cosmic-ray-induced gamma rays from the Perseus cluster for the Cherenkov Telescope Array.

preprint2020arXiv

Cross-correlation of the thermal Sunyaev-Zel'dovich effect and weak gravitational lensing: Planck and Subaru Hyper Suprime-Cam first-year data

Cross-correlation analysis of the thermal Sunyaev-Zel'dovich (tSZ) effect and weak gravitational lensing (WL) provides a powerful probe of cosmology and astrophysics of the intra-cluster medium. We present the measurement of the cross-correlation of tSZ and WL from Planck and Subaru Hyper-Suprime Cam. The combination enables us to study cluster astrophysics at high redshift. We use the tSZ-WL cross-correlation and the tSZ auto-power spectrum measurements to place a tight constraint on the hydrostatic mass bias, which is a measure of the degree of non-thermal pressure support in galaxy clusters. With the prior on cosmological parameters derived from the analysis of the cosmic microwave background anisotropies by Planck and taking into account foreground contributions both in the tSZ auto-power spectrum and the tSZ-WL cross-correlation, the hydrostatic mass bias is estimated to be $26.9^{+8.9}_{-4.4} \%$ ($68\%$ C.L.), which is consistent with recent measurements by mass calibration techniques.

preprint2020arXiv

Deep learning for intensity mapping observations: Component extraction

Line intensity mapping (LIM) is an emerging observational method to study the large-scale structure of the Universe and its evolution. LIM does not resolve individual sources but probes the fluctuations of integrated line emissions. A serious limitation with LIM is that contributions of different emission lines from sources at different redshifts are all confused at an observed wavelength. We propose a deep learning application to solve this problem. We use conditional generative adversarial networks to extract designated information from LIM. We consider a simple case with two populations of emission line galaxies; H$\rmα$ emitting galaxies at $z = 1.3$ are confused with [OIII] emitters at $z = 2.0$ in a single observed waveband at 1.5 $\rmμ$m. Our networks trained with 30,000 mock observation maps are able to extract the total intensity and the spatial distribution of H$\rmα$ emitting galaxies at $z = 1.3$. The intensity peaks are successfully located with 74% precision. The precision increases to 91% when we combine the results of 5 networks. The mean intensity and the power spectrum are reconstructed with an accuracy of $\sim$10%. The extracted galaxy distributions at a wider range of redshift can be used for studies on cosmology and on galaxy formation and evolution.

preprint2020arXiv

Deep Learning for Line Intensity Mapping Observations: Information Extraction from Noisy Maps

Line intensity mapping (LIM) is a promising observational method to probe large-scale fluctuations of line emission from distant galaxies. Data from wide-field LIM observations allow us to study the large-scale structure of the universe as well as galaxy populations and their evolution. A serious problem with LIM is contamination by foreground/background sources and various noise contributions. We develop conditional generative adversarial networks (cGANs) that extract designated signals and information from noisy maps. We train the cGANs using 30,000 mock observation maps with assuming a Gaussian noise matched to the expected noise level of NASA's SPHEREx mission. The trained cGANs successfully reconstruct Hα emission from galaxies at a target redshift from observed, noisy intensity maps. Intensity peaks with heights greater than 3.5 σ noise are located with 60 % precision. The one-point probability distribution and the power spectrum are accurately recovered even in the noise-dominated regime. However, the overall reconstruction performance depends on the pixel size and on the survey volume assumed for the training data. It is necessary to generate training mock data with a sufficiently large volume in order to reconstruct the intensity power spectrum at large angular scales. Our deep-learning approach can be readily applied to observational data with line confusion and with noise.

preprint2020arXiv

First Structure Formation under the Influence of Gas-Dark Matter Streaming Velocity and Density: Impact of the Baryons-trace-dark matter Approximation

The impact of the streaming between baryons and dark matter on the first structures has been actively explored by recent studies. We investigate how the key results are affected by two popular approximations. One is to implement the streaming by accounting for only the relative motion while assuming ``baryons trace dark matter" spatially at the initialization of simulation. This neglects the smoothing on the gas density taking place before the initialization. In our simulation initialized at $z_i=200$, it overestimates the gas density power spectrum by up to 40\% at $k\approx10^2~h~\mbox{Mpc}^{-1}$ at $z=20$. Halo mass ($M_h$) and baryonic fraction in halos ($f_{b,h}$) are also overestimated, but the relation between the two remains unchanged. The other approximation tested is to artificially amplify the density/velocity fluctuations in the cosmic mean density to simulate the first minihalos that form in overdense regions. This gives a head start to the halo growth while the subsequent growth is similar to that in the mean density. The growth in a true overdense region, on the other hand, is accelerated gradually in time. For example, raising $σ_8$ by 50\% effectively transforms $z\rightarrow\sqrt{1.5}z$ in the halo mass growth history while in 2-$σ$ overdensity, the growth is accelerated by a constant in redshift: $z\rightarrow{z+4.8}$. As a result, halos have grown more in the former than in the latter before $z\approx27$ and vice versa after. The $f_{b,h}$-$M_h$ relation is unchanged in those cases as well, suggesting that the Pop III formation rate for a given $M_h$ is insensitive to the tested approximations.

preprint2020arXiv

Photoevaporation of Minihalos during Cosmic Reionization: Primordial and Metal-Enriched Halos

The density distribution of the inter-galactic medium is an uncertain but highly important issue in the study of cosmic reionization. It is expected that there are abundant gas clouds hosted by low-mass "minihalos" in the early universe, which act as photon sinks until photoevaporated by the emerging ultra-violet background (UVB) radiation. We perform a suite of radiation hydrodynamics simulations to study the photoevaporation of minihalos. Our simulations follow hydrodynamics, non-equilibrium chemistry, and the associated cooling processes in a self-consistent manner. We conduct a parametric study by considering a wide range of gas metallicity ($0\,Z_\odot \leq Z \leq 10^{-3}\,Z_\odot$), halo mass ($10^3 M_\odot \leq M \leq 10^8 M_\odot$), UVB intensity ($0.01 \leq J_{21} \leq 1$), and turn-on redshift of ionizing sources ($10\leq z_{\rm IN} \leq 20$). We show that small halos are evaporated in a few tens million years, whereas larger mass halos survive for ten times longer. We show that the gas mass evolution of a minihalo can be characterized by a scaling parameter that is given by a combination of the halo mass, background radiation intensity, and redshift. Efficient radiative cooling in metal-enriched halos induces fast condensation of the gas to form a dense, self-shielded core. The cold, dense core can become gravitationally unstable in halos with high metallicities. Early metal enrichment may allow star formation in minihalos during cosmic reionization.

preprint2020arXiv

Photometric classification of HSC transients using machine learning

The advancement of technology has resulted in a rapid increase in supernova (SN) discoveries. The Subaru/Hyper Suprime-Cam (HSC) transient survey, conducted from fall 2016 through spring 2017, yielded 1824 SN candidates. This gave rise to the need for fast type classification for spectroscopic follow-up and prompted us to develop a machine learning algorithm using a deep neural network (DNN) with highway layers. This machine is trained by actual observed cadence and filter combinations such that we can directly input the observed data array into the machine without any interpretation. We tested our model with a dataset from the LSST classification challenge (Deep Drilling Field). Our classifier scores an area under the curve (AUC) of 0.996 for binary classification (SN Ia or non-SN Ia) and 95.3% accuracy for three-class classification (SN Ia, SN Ibc, or SN II). Application of our binary classification to HSC transient data yields an AUC score of 0.925. With two weeks of HSC data since the first detection, this classifier achieves 78.1% accuracy for binary classification, and the accuracy increases to 84.2% with the full dataset. This paper discusses the potential use of machine learning for SN type classification purposes.

preprint2020arXiv

R-process enrichment in ultrafaint dwarf galaxies

We study the enrichment and mixing of r-process elements in ultrafaint dwarf galaxies (UFDs). We assume that r-process elements are produced by neutron-star mergers (NSMs), and examine multiple models with different natal kick velocities and explosion energies. To this end, we perform cosmological simulations of galaxy formation to follow mixing of the dispersed r-process elements driven by star formation and the associated stellar feedback in progenitors of UFDs. We show that the observed europium abundance in Reticulum II is reproduced by our inner explosion model where a NSM is triggered at the centre of the galaxy, whereas the relatively low abundance in Tucana III is reproduced if a NSM occurs near the virial radius of the progenitor galaxy. The latter case is realised only if the neutron-star binary has a large natal kick velocity and travels over a long distance of a kilo-parsec before merger. In both the inner and outer explosion cases, it is necessary for the progenitor galaxy to sustain prolonged star formation over a few hundred million years after the NSM, so that the dispersed r-process elements are well mixed within the inter-stellar medium. Short-duration star formation results in inefficient mixing, and then a large variation is imprinted in the stellar europium abundances, which is inconsistent with the observations of Reticulum II and Tucana III.

preprint2020arXiv

Rapid Transients Originating from Thermonuclear Explosions in Helium White Dwarf Tidal Disruption Events

We study the emission properties of thermonuclear explosions in a helium white dwarf (WD) tidal disruption event (TDE). We consider a TDE where a 0.2 $M_{\odot}$ helium WD is disrupted by a $10^{2.5}\,M_{\odot}$ intermediate-mass black hole (IMBH). The helium WD is not only tidally disrupted but is also detonated by the tidal compression and by succeeding shocks. We focus on the emission powered by radioactive nuclei in the unbound TDE ejecta. We perform hydrodynamic simulations coupled with nuclear reactions, post-process detailed nucleosynthesis calculations, and then radiative transfer simulations. We thus derive multi-band light curves and spectra. The helium WD TDE shows rapid ($Δt_{1\mathrm{mag}}\simeq5\text{--}10$ days) and relatively faint ($L_{\mathrm{peak}}\simeq10^{42}\,\mathrm{erg}\,\mathrm{s}^{-1}$) light curves, because the ejecta mass and $^{56}$Ni mass are low ($0.12\,M_{\odot}$ and $0.03\,M_{\odot}$, respectively). The spectra show strong calcium and Fe-peak features and very weak silicon features, reflecting the peculiar elemental abundance. The key feature is the Doppler shift of the spectral lines up to $\simeq\pm12,000\,\mathrm{km}\,\mathrm{s}^{-1}$, depending on the viewing angle, due to the bulk motion of the ejecta. Our model matches well with two rapid and faint transients reported in Pursiainen et al. (2018). The particular model presented here does not match with observed SNe Iax, calcium-rich transients, or .Ia explosion candidates, either in the spectra or light curves. However, we expect a large variety of the observational signatures once a wide range of the WD/black hole masses and orbital parameters are considered. This study contributes to the search for WD TDEs with current and upcoming surveys, and to the identification of IMBHs as disrupters in the TDEs.

preprint2018arXiv

The onset of star formation 250 million years after the Big Bang

A fundamental quest of modern astronomy is to locate the earliest galaxies and study how they influenced the intergalactic medium a few hundred million years after the Big Bang. The abundance of star-forming galaxies is known to decline from redshifts of about 6 to 10, but a key question is the extent of star formation at even earlier times, corresponding to the period when the first galaxies might have emerged. Here we present spectroscopic observations of MACS1149-JD1, a gravitationally lensed galaxy observed when the Universe was less than four per cent of its present age. We detect an emission line of doubly ionized oxygen at a redshift of $9.1096\pm0.0006$, with an uncertainty of one standard deviation. This precisely determined redshift indicates that the red rest-frame optical colour arises from a dominant stellar component that formed about 250 million years after the Big Bang, corresponding to a redshift of about 15. Our results indicate the it may be possible to detect such early episodes of star formation in similar galaxies with future telescopes.

preprint2016arXiv

Cosmological Constraint on the Light Gravitino Mass from CMB Lensing and Cosmic Shear

Light gravitinos of mass $\lesssim \mathcal{O} (10)$ eV are of particular interest in cosmology, offering various baryogenesis scenarios without suffering from the cosmological gravitino problem. The gravitino may contribute considerably to the total matter content of the Universe and affect structure formation from early to present epochs. After the gravitinos decouple from other particles in the early Universe, they free-stream and consequently suppress density fluctuations of (sub-)galactic length scales. Observations of structure at the relevant length-scales can be used to infer or constrain the mass and the abundance of light gravitinos. We derive constraints on the light gravitino mass using the data of cosmic microwave background (CMB) lensing from Planck and of cosmic shear from the Canada France Hawaii Lensing Survey, combined with analyses of the primary CMB anisotropies and the signature of baryon acoustic oscillations in galaxy distributions. The obtained constraint on the gravitino mass is $m_{3/2} < 4.7$ eV (95% C.L.), which is substantially tighter than the previous constraint from clustering analysis of Ly-$α$ forests.

preprint2016arXiv

Cosmological constraints on dark matter annihilation and decay: Cross-correlation analysis of the extragalactic $γ$-ray background and cosmic shear

We derive constraints on dark matter (DM) annihilation cross section and decay lifetime from cross-correlation analyses of the data from Fermi-LAT and weak lensing surveys that cover a wide area of $\sim660$ squared degrees in total. We improve upon our previous analyses by using an updated extragalactic $γ$-ray background data reprocessed with the Fermi Pass 8 pipeline, and by using well-calibrated shape measurements of about twelve million galaxies in the Canada-France-Hawaii Lensing Survey (CFHTLenS) and Red-Cluster-Sequence Lensing Survey (RCSLenS). We generate a large set of full-sky mock catalogs from cosmological $N$-body simulations and use them to estimate statistical errors accurately. The measured cross correlation is consistent with null detection, which is then used to place strong cosmological constraints on annihilating and decaying DM. For leptophilic DM, the constraints are improved by a factor of $\sim100$ in the mass range of O(1) TeV when including contributions from secondary $γ$ rays due to the inverse-Compton upscattering of background photons. Annihilation cross-sections of $\langle σv \rangle \sim 10^{-23}\, {\rm cm}^3/{\rm s}$ are excluded for TeV-scale DM depending on channel. Lifetimes of $\sim 10^{25}$ sec are also excluded for the decaying TeV-scale DM. Finally, we apply this analysis to wino DM and exclude the wino mass around 200 GeV. These constraints will be further tightened, and all the interesting wino DM parameter region can be tested, by using data from future wide-field cosmology surveys.

preprint2016arXiv

Cosmological Simulations of Early Blackhole Formation: Halo Mergers, Tidal Disruption, and the Conditions for Direct Collapse

Gravitational collapse of a massive primordial gas cloud is thought to be a promising path for the formation of supermassive blackholes in the early universe. We study conditions for the so-called direct collapse (DC) blackhole formation in a fully cosmological context. We combine a semianalytic model of early galaxy formation with halo merger trees constructed from dark matter $N$-body simulations. We locate a total of 68 possible DC sites in a volume of $20\;h^{-1}\;\mathrm{Mpc}$ on a side. We then perform hydrodynamics simulations for 42 selected halos to study in detail the evolution of the massive clouds within them. We find only two successful cases where the gas clouds rapidly collapse to form stars. In the other cases, gravitational collapse is prevented by the tidal force exerted by a nearby massive halo, which otherwise should serve as a radiation source necessary for DC. Ram pressure stripping disturbs the cloud approaching the source. In many cases, a DC halo and its nearby light source halo merge before the onset of cloud collapse. Only when the DC halo is assembled through major mergers, the gas density increases rapidly to trigger gravitational instability. Based on our cosmological simulations, we conclude that the event rate of DC is an order of magnitude smaller than reported in previous studies, although the absolute rate is still poorly constrained. It is necessary to follow the dynamical evolution of a DC cloud and its nearby halo(s) in order to determine the critical radiation flux for DC.

preprint2016arXiv

Detection of an oxygen emission line from a high redshift galaxy in the reionization epoch

The physical properties and elemental abundances of the interstellar medium in galaxies during cosmic reionization are important for understanding the role of galaxies in this process. We report the Atacama Large Millimeter/submillimeter Array detection of an oxygen emission line at a wavelength of 88 micrometers from a galaxy at an epoch about 700 million years after the Big Bang. The oxygen abundance of this galaxy is estimated at about one-tenth that of the Sun. The non-detection of far-infrared continuum emission indicates a deficiency of interstellar dust in the galaxy. A carbon emission line at a wavelength of 158 micrometers is also not detected, implying an unusually small amount of neutral gas. These properties might allow ionizing photons to escape into the intergalactic medium.

preprint2016arXiv

Eigenanalysis of morphological diversity in silicon random nanostructures formed via resist collapse

This paper demonstrates eigenanalysis to quantitatively reveal the diversity and capacity of identities offered by the morphological diversity in silicon nanostructures formed via random collapse of resist. The analysis suggests that approximately 10^115 possible identities are provided per 0.18-um^2 area of nanostructures, indicating that nanoscale morphological signatures will be extremely useful for future information security applications where securing identities is critical. The eigenanalysis provides an intuitive physical picture and quantitative characterization of the diversity of structural fluctuations while unifying measurement stability concerns, which will be widely applicable to other materials, devices, and system architectures.

preprint2016arXiv

Formation of Massive Primordial Stars: Intermittent UV Feedback with Episodic Mass Accretion

We present coupled stellar evolution (SE) and 3D radiation-hydrodynamic (RHD) simulations of the evolution of primordial protostars, their immediate environment, and the dynamic accretion history under the influence of stellar ionizing and dissociating UV feedback. Our coupled SE-RHD calculations result in a wide diversity of final stellar masses covering 10 Msun $\lesssim M_* \lesssim$ 1000 Msun. The formation of very massive ($\gtrsim$ 250 Msun) stars is possible under weak UV feedback, whereas ordinary massive (a few x 10 Msun) stars form when UV feedback can efficiently halt the accretion. This may explain the peculiar abundance pattern of a Galactic metal-poor star recently reported by Aoki et al. (2014), possibly the observational signature of very massive precursor primordial stars. Weak UV feedback occurs in cases of variable accretion, in particular when repeated short accretion bursts temporarily exceed 0.01 Msun/yr, causing the protostar to inflate. In the bloated state, the protostar has low surface temperature and UV feedback is suppressed until the star eventually contracts, on a thermal adjustment timescale, to create an HII region. If the delay time between successive accretion bursts is sufficiently short, the protostar remains bloated for extended periods, initiating at most only short periods of UV feedback. Disk fragmentation does not necessarily reduce the final stellar mass. Quite the contrary, we find that disk fragmentation enhances episodic accretion as many fragments migrate inward and are accreted onto the star, thus allowing continued stellar mass growth under conditions of intermittent UV feedback. This trend becomes more prominent as we improve the resolution of our simulations. We argue that simulations with significantly higher resolution than reported previously are needed to derive accurate gas mass accretion rates onto primordial protostars.

preprint2016arXiv

Machine-learning Selection of Optical Transients in Subaru/Hyper Suprime-Cam Survey

We present an application of machine-learning (ML) techniques to source selection in the optical transient survey data with Hyper Suprime-Cam (HSC) on the Subaru telescope. Our goal is to select real transient events accurately and in a timely manner out of a large number of false candidates, obtained with the standard difference-imaging method. We have developed the transient selector which is based on majority voting of three ML machines of AUC Boosting, Random Forest, and Deep Neural Network. We applied it to our observing runs of Subaru-HSC in 2015 May and August, and proved it to be efficient in selecting optical transients. The false positive rate was 1.0% at the true positive rate of 90% in the magnitude range of 22.0--25.0 mag for the former data. For the latter run, we successfully detected and reported ten candidates of supernovae within the same day as the observation. From these runs, we learned the following lessons: (1) the training using artificial objects is effective in filtering out false candidates, especially for faint objects, and (2) combination of ML by majority voting is advantageous.

preprint2016arXiv

Matter Power Spectrum in Hidden Neutrino Interacting Dark Matter Models: A Closer Look at the Collision Term

Dark Matter (DM) models providing possible alternative solutions to the small- scale crisis of standard cosmology are nowadays of growing interest. We consider DM interacting with light hidden fermions via well motivated fundamental operators showing the resultant matter power spectrum is suppressed on subgalactic scales within a plausible parameter region. Our basic description of the evolution of cosmological perturbations relies on a fully consistent first principles derivation of a perturbed Fokker-Planck type equation, generalizing existing literature. The cosmological perturbation of the Fokker-Planck equation is presented for the first time in two different gauges, where the results transform into each other according to the rules of gauge transformation. Furthermore, our focus lies on a derivation of a broadly applicable and easily computable collision term showing important phenomenological differences to other existing approximations. As one of the main results and concerning the small-scale crisis, we show the equal importance of vector and scalar boson mediated interactions between DM and light fermions.

preprint2016arXiv

Nebular line emission from z > 7 galaxies in cosmological simulations: rest-frame UV to Optical lines

We have performed very large and high resolution cosmological hydrodynamic simulations in order to investigate detectability of nebular lines in the rest-frame ultraviolet (UV) to optical wavelength range from galaxies at $z > 7$. We find that the expected line fluxes are very well correlated with apparent UV magnitudes. The C IV $1549 \rm Å$ and the C III] $1909 \rm Å$ lines of galaxies brighter than $26~ \rm AB$ are detectable with current facilities such as the VLT/XShooter and the Keck/MOSFIRE. Metal lines such as C IV $1549 \rm Å$, C III] $1909 \rm Å$, [O II] $3727 \rm Å$ and [O III] $4959/5007 \rm Å$ can be good targets for the spectroscopic observation with Thirty Meter Telescope (TMT), European Extremely Large Telescope (E-ELT), Giant Magellan Telescope (GMT) and James Webb Space Telescope (JWST). H$α$ and H$β$ lines are also expected to be detectable with these telescopes. We also predict detectability of the nebular lines for $z > 10$ galaxies that will be found with JWST, Wide-Field Infrared Survey Telescope (WFIRST) and First Light And Reionization Explorer (FLARE) $(11 \leq z \leq 15$). We conclude that the C IV $1549 \rm Å$, C III] $1909 \rm Å$, [O III] $4959/5007 \rm Å$ and H$β$ lines even from $z \sim 15$ galaxies can be strong targets for TMT, ELT and JWST. We also find that the magnification by gravitational lensing is of great help to detect such high-$z$ galaxies. According to our model, the C III] $1909 \rm Å$ line in $z > 9$ galaxy candidates is detectable even using the current facilities.

preprint2016arXiv

Optical nano artifact metrics using silicon random nanostructures

Nano artifact metrics exploit unique physical attributes of nanostructured matter for authentication and clone resistance, which is vitally important in the age of Internet-of-Things where securing identities is critical. However, high-cost and huge experimental apparatuses, such as scanning electron microscopy, have been required in the former studies. Herein, we demonstrate an optical approach to characterise the nanoscale-precision signatures of silicon random structures towards realising low-cost and high-value information security technology. Unique and versatile silicon nanostructures are generated via resist collapse phenomena, which contains dimensions that are well below the diffraction limit of light. We exploit the nanoscale precision ability of confocal laser microscopy in the height dimension, and our experimental results demonstrate that the vertical precision of measurement is essential in satisfying the performances required for artifact metrics. Furthermore, by using state-of-the-art nanostructuring technology, we experimentally fabricate clones from the genuine devices. We demonstrate that the statistical properties of the genuine and clone devices are successfully exploited, showing that the liveness-detection-type approach, which is widely deployed in biometrics, is valid in artificially-constructed solid-state nanostructures. These findings pave the way for reasonable and yet sufficiently secure novel principles for information security based on silicon random nanostructures and optical technologies.

preprint2016arXiv

The Final Fates of Accreting Supermassive Stars

The formation of supermassive stars (SMSs) via rapid mass accretion and their direct collapse into black holes (BHs) is a promising pathway for sowing seeds of supermassive BHs in the early universe. We calculate the evolution of rapidly accreting SMSs by solving the stellar structure equations including nuclear burning as well as general relativistic (GR) effects up to the onset of the collapse. We find that such SMSs have less concentrated structure than fully-convective counterpart, which is often postulated for non-accreting ones. This effect stabilizes the stars against GR instability even above the classical upper mass limit $\gtrsim 10^5~M_\odot$ derived for the fully-convective stars. The accreting SMS begins to collapse at the higher mass with the higher accretion rate. The collapse occurs when the nuclear fuel is exhausted only for cases with $\dot M \lesssim 0.1~M_\odot~{\rm yr}^{-1}$. With $\dot{M} \simeq 0.3 - 1~M_\odot~{\rm yr}^{-1}$, the star becomes GR-unstable during the helium-burning stage at $M \simeq 2 - 3.5~\times 10^5~M_\odot$. In an extreme case with $10~M_\odot~{\rm yr}^{-1}$, the star does not collapse until the mass reaches $\simeq 8.0\times 10^5~M_\odot$, where it is still in the hydrogen-burning stage. We expect that BHs with roughly the same mass will be left behind after the collapse in all the cases.

preprint2016arXiv

The Hydrodynamic Feedback of Cosmic Reionization on Small-Scale Structures and Its Impact on Photon Consumption during the Epoch of Reionization

Density inhomogeneity in the intergalactic medium (IGM) can boost the recombination rate of ionized gas substantially, affecting the growth of HII regions during reionization. Previous attempts to quantify this effect typically failed to resolve down to the Jeans scale in the pre-ionization IGM, which is important in establishing this effect, along with the hydrodynamical back-reaction of reionization on it. Towards that end, we perform a set of fully-coupled, radiation-hydrodynamics simulations from cosmological initial conditions, extending the mass resolution of previous work to the scale of minihalos. Pre-reionization structure is evolved until a redshift $z_i$ at which the ionizing radiation from external sources arrives to sweep an R-type ionization front supersonically across the volume in a few Myr, until it is trapped on the surfaces of minihalos and converted to D-type, after which the minihalo gas is removed by photoevaporative winds. Small-scale density structures during this time lead to a high ($>$10) clumping factor for ionized gas, which hugely boosts the recombination rate until the structures are disrupted by the hydrodynamic feedback after $\sim 10-100~\rm{Myr}$. For incoming stellar radiation with intensity $J_{21}$ in a $200~h^{-1}~\rm{kpc}$ box with the mean density contrast $\barδ$, the number of extra recombinations per H atom, on top of what is expected from homogeneously distributed gas, is given by $0.32[J_{21}]^{0.12}[(1+z_i)/11]^{-1.7}[1+\barδ]^{2.5}$. In models in which most of the volume is ionized toward the end of reionization, this can add more than one recombination per H atom to the ionizing photon budget to achieve reionization.

preprint2016arXiv

The Subaru FMOS galaxy redshift survey (FastSound). IV. New constraint on gravity theory from redshift space distortions at $z\sim 1.4$

We measure the redshift-space correlation function from a spectroscopic sample of 2783 emission line galaxies from the FastSound survey. The survey, which uses the Subaru Telescope and covers the redshift ranges of $1.19<z<1.55$, is the first cosmological study at such high redshifts. We detect clear anisotropy due to redshift-space distortions (RSD) both in the correlation function as a function of separations parallel and perpendicular to the line of sight and its quadrupole moment. RSD has been extensively used to test general relativity on cosmological scales at $z<1$. Adopting a LCDM cosmology with the fixed expansion history and no velocity dispersion $σ_{\rm v}=0$, and using the RSD measurements on scales above 8Mpc/h, we obtain the first constraint on the growth rate at the redshift, $f(z)σ_8(z)=0.482\pm 0.116$ at $z\sim 1.4$ after marginalizing over the galaxy bias parameter $b(z)σ_8(z)$. This corresponds to $4.2σ$ detection of RSD. Our constraint is consistent with the prediction of general relativity $fσ_8\sim 0.392$ within the $1-σ$ confidence level. When we allow $σ_{\rm v}$ to vary and marginalize it over, the growth rate constraint becomes $fσ_8=0.494^{+0.126}_{-0.120}$. We also demonstrate that by combining with the low-z constraints on $fσ_8$, high-z galaxy surveys like the FastSound can be useful to distinguish modified gravity models without relying on CMB anisotropy experiments.

preprint2015arXiv

Cross-Correlation of the Extragalactic Gamma-ray Background with Luminous Red Galaxies

Measurements of the cross-correlation between the extragalactic gamma-ray background (EGB) and large-scale structure provide a novel probe of dark matter on extragalactic scales. We focus on luminous red galaxies (LRGs) as optimal targets to search for the signal of dark matter annihilation. We measure the cross-correlation function of the EGB taken from the Fermi Large Area Telescope with the LRGs from the Sloan Digital Sky Survey. Statistical errors are calculated using a large set of realistic mock LRG catalogs. The amplitude of the measured cross-correlation is consistent with null detection. Based on an accurate theoretical model of the distribution of dark matter associated with LRGs, we exclude dark matter annihilation cross-sections over $\langle σv\rangle =3\times10^{-25}-10^{-26}\, {\rm cm}^3 \,{\rm s}^{-1}$ for a 10 GeV dark matter. We further investigate systematic effects due to uncertainties in the Galactic gamma-ray foreground emission, which we find to be an order of magnitude smaller than the current statistical uncertainty. We also estimate the contamination from astrophysical sources in the LRGs by using known scaling relations between gamma-ray luminosity and star-formation rate, finding them to be negligibly small. Based on these results, we suggest that LRGs remain ideal targets for probing dark matter annihilation with future EGB measurement and galaxy surveys. Increasing the number of LRGs in upcoming galaxy surveys such as LSST would lead to big improvements of factors of several in sensitivity.

preprint2015arXiv

Early structure formation from primordial density fluctuations with a blue-tilted power spectrum

While observations of large-scale structure and the cosmic microwave background (CMB) provide strong constraints on the amplitude of the primordial power spectrum (PPS) on scales larger than 10~Mpc, the amplitude of the power spectrum on sub-galactic length scales is much more poorly constrained. We study early structure formation in a cosmological model with a blue-tilted PPS. We assume that the standard scale-invariant PPS is modified at small length scales as $P(k) \sim k^{m_{\rm s}}$ with $m_{\rm s} > 1$. We run a series of cosmological hydrodynamic simulations to examine the dependence of the formation epoch and the characteristic mass of primordial stars on the tilt of the PPS. In models with $m_{\rm s} > 1$, star-forming gas clouds are formed at $z > 100$, when formation of hydrogen molecules is inefficient because the intense CMB radiation destroys chemical intermediates. Without efficient coolant, the gas clouds gravitationally contract while keeping a high temperature. The protostars formed in such "hot" clouds grow very rapidly by accretion to become extremely massive stars that may leave massive black holes with a few hundred solar-masses at $z > 100$. The shape of the PPS critically affects the properties and the formation epoch of the first generation of stars. Future experiments of the CMB polarization and the spectrum distortion may provide important information on the nature of the first stars and their formation epoch, and hence on the shape of the small-scale power spectrum.

preprint2015arXiv

Impact of Baryonic Processes on Weak Lensing Cosmology: Power Spectrum, Non-Local Statistics, and Parameter Bias

We study the impact of baryonic physics on cosmological parameter estimation with weak lensing surveys. We run a set of cosmological hydrodynamics simulations with different galaxy formation models. We then perform ray-tracing simulations through the total matter density field to generate 100 independent convergence maps of 25 ${\rm deg}^2$ field-of-view, and use them to examine the ability of the following three lensing statistics as cosmological probes; power spectrum, peak counts, and Minkowski functionals. For the upcoming wide-field observations such as Subaru Hyper Suprime-Cam (HSC) survey with a sky coverage of 1400 ${\rm deg}^2$, these three statistics provide tight constraints on the matter density, density fluctuation amplitude, and dark energy equation of state, but parameter bias is induced by the baryonic processes such as gas cooling and stellar feedback. The bias can be significant when the statistical errors become small in future observations with a much larger survey area. We find the bias is induced in different directions in the parameter space depending on the statistics employed. While the two-point statistic, i.e. power spectrum, yields robust results against baryonic effects, the overall constraining power is weak compared with peak counts and Minkowski functionals. On the other hand, using one of peak counts or Minkowski functionals, or combined analysis with multiple statistics, results in biased parameter estimate. The bias can be as large as $1σ$ for HSC survey, and will be more significant for upcoming wider area surveys. We suggest to use an optimized combination so that the baryonic effects on parameter estimation are mitigated. Such `calibrated' combination can place stringent and robust constraints on cosmological parameters.

preprint2015arXiv

Particle splitting in smoothed particle hydrodynamics based on Voronoi diagram

We present a novel method for particle splitting in smoothed particle hydrodynamics simulations. Our method utilizes the Voronoi diagram for a given particle set to determine the position of fine daughter particles. We perform several test simulations to compare our method with a conventional splitting method in which the daughter particles are placed isotropically over the local smoothing length. We show that, with our method, the density deviation after splitting is reduced by a factor of about two compared with the conventional method. Splitting would smooth out the anisotropic density structure if the daughters are distributed isotropically, but our scheme allows the daughter particles to trace the original density distribution with length scales of the mean separation of their parent. We apply the particle splitting to simulations of the primordial gas cloud collapse. The thermal evolution is accurately followed to the hydrogen number density of 10^12 /cc. With the effective mass resolution of ~10^-4 Msun after the multi-step particle splitting, the protostellar disk structure is well resolved. We conclude that the method offers an efficient way to simulate the evolution of an interstellar gas and the formation of stars.

preprint2015arXiv

Primordial Star Formation under the Influence of Far Ultraviolet Radiation: 1540 Cosmological Halos and the Stellar Mass Distribution

We perform a large set of cosmological simulations of early structure formation and follow the formation and evolution of 1540 star-forming gas clouds to derive the mass distribution of primordial stars. The star formation in our cosmological simulations is characterized by two distinct populations, the so-called Population III.1 stars and primordial stars formed under the influence of far ultraviolet (FUV) radiation (Population III.2D stars). In this work, we determine the stellar masses by using the dependences on the physical properties of star-forming cloud and/or the external photodissociating intensity from nearby primordial stars, which are derived from the results of two-dimensional radiation hydrodynamic simulations of protostellar feedback. The characteristic mass of the Pop III stars is found to be a few hundred solar masses at z ~ 25, and it gradually shifts to lower masses with decreasing redshift. At high redshifts z > 20, about half of the star-forming gas clouds are exposed to intense FUV radiation and thus give birth to massive Pop III.2D stars. However, the local FUV radiation by nearby Pop III stars becomes weaker at lower redshifts, when typical Pop III stars have smaller masses and the mean physical separation between the stars becomes large owing to cosmic expansion. Therefore, at z < 20, a large fraction of the primordial gas clouds host Pop III.1 stars. At z =< 15, the Pop III.1 stars are formed in relatively cool gas clouds due to efficient radiative cooling by H_2 and HD molecules; such stars have masses of a few x 10 Msun. Since the stellar evolution and the final fate are determined by the stellar mass, Pop III stars formed at different epochs play different roles in the early universe.

preprint2015arXiv

Probing cosmology with weak lensing selected clusters I: Halo approach and all-sky simulations

We explore a variety of statistics of clusters selected with cosmic shear measurement by utilizing both analytic models and large numerical simulations. We first develop a halo model to predict the abundance and the clustering of weak lensing selected clusters. Observational effects such as galaxy shape noise are included in our model. We then generate realistic mock weak lensing catalogs to test the accuracy of our analytic model. To this end, we perform full-sky ray-tracing simulations that allow us to have multiple realizations of a large continuous area. We model the masked regions on the sky using the actual positions of bright stars, and generate 200 mock weak lensing catalogs with sky coverage of ~1000 squared degrees. We show that our theoretical model agrees well with the ensemble average of statistics and their covariances calculated directly from the mock catalogues. With a typical selection threshold, ignoring shape noise correction causes overestimation of the clustering of weak lensing selected clusters with a level of about 10%, and shape noise correction boosts the cluster abundance by a factor of a few. We calculate the cross-covariances using the halo model with accounting for the effective reduction of the survey area due to masks. The covariance of the cosmic shear auto power spectrum is affected by the mode-coupling effect that originates from sky masking. Our model and the results can be readily used for cosmological analysis with ongoing and future weak lensing surveys.

preprint2015arXiv

Probing cosmology with weak lensing selected clusters II: Dark energy and f(R) gravity models

Ongoing and future wide-field galaxy surveys can be used to locate a number of clusters of galaxies with cosmic shear measurement alone. We study constraints on cosmological models using statistics of weak lensing selected galaxy clusters. We extend our previous theoretical framework to model the statistical properties of clusters in variants of cosmological models as well as in the standard $Λ$CDM model. Weak lensing selection of clusters does not rely on the conventional assumption such as the relation between luminosity and mass and/or hydrostatic equilibrium, but a number of observational effects compromise robust identification. We use a large set of realistic mock weak-lensing catalogs as well as analytic models to perform a Fisher analysis and make forecast for constraining two competing cosmological models, $w$CDM model and $f(R)$ model proposed by Hu & Sawicki, with our lensing statistics. We show that weak lensing selected clusters are excellent probe of cosmology when combined with cosmic shear power spectrum even in presence of galaxy shape noise and masked regions. With the information of weak lensing selected clusters, the precision of cosmological parameter estimate can be improved by a factor of $\sim1.6$ and $\sim8$ for $w$CDM model and $f(R)$ model, respectively. Hyper Suprime-Cam survey with sky coverage of $1250$ squared degrees can constrain the equation of state of dark energy $w_{0}$ with a level of $Δw_0 \sim0.1$. It can also constrain the additional scalar degree of freedom in $f(R)$ model with a level of $|f_{R0}| \sim5\times10^{-6}$, when constraints from cosmic microwave background measurements are incorporated. Future weak lensing surveys with sky coverage of $20,000$ squared degrees will place tighter constraints on $w_{0}$ and $|f_{R0}|$ even without cosmic microwave background measurements.

preprint2015arXiv

Probing the Small Scale Matter Power Spectrum through Dark Matter Annihilation in the Early Universe

Recent observations of the cosmic microwave background (CMB) anisotropies and the distribution of galaxies, galaxy clusters, and the Lyman Alpha forest have constrained the shape of the power spectrum of matter fluctuations on large scales k < few h/Mpc. We explore a new technique to constrain the matter power spectrum on smaller scales, assuming the dark matter is a Weakly Interacting Massive Particle (WIMP) that annihilates at early epochs. Energy released by dark matter annihilation can modify the spectrum of CMB temperature fluctuations and thus CMB experiments such as Planck have been able to constrain the quantity f <sigma v> /m < 1/88 picobarn c / GeV, where f is the fraction of energy absorbed by gas, <sigma v> is the annihilation rate assumed constant, and m is the particle mass. We assume the standard scale-invariant primordial matter power spectrum of P_prim(k) ~ k^{n_s} at large scales k < k_p, while we adopt the modified power law of P_prim(k) ~ k_p^{n_s} (k/k_p)^{m_s} at small scales. We then aim at deriving constraints on m_s. For m_s > n_s, the excess small-scale power results in a much larger number of nonlinear small mass halos, particularly at high redshifts. Dark matter annihilation in these halos releases sufficient energy to partially ionize the gas, and consequently modify the spectrum of CMB fluctuations. We show that the recent Planck data can already be used to constrain the power spectrum on small scales. For a simple model with an NFW profile with halo concentration parameter c_200 = 5 and f <sigma v> / m = 1/100 picobarn c / GeV, we can limit the mass variance sigma_{max} < 100 at the 95% confidence level, corresponding to a power law index m_s < 1.43 (1.63) for k_p = 100 (1000) h/Mpc. Our results are also relevant to theories that feature a running spectral index.

preprint2015arXiv

Supermassive star formation via episodic accretion: protostellar disc instability and radiative feedback efficiency

The formation of SMSs is a potential pathway to seed SMBHs in the early universe. A critical issue for forming SMSs is stellar UV feedback, which may limit the stellar mass growth via accretion. In this paper we study the evolution of an accreting SMS and its UV emissivity under conditions of realistic variable accretion from a self-gravitating circumstellar disc. First we conduct a 2D hydrodynamical simulation to follow the long-term protostellar accretion until the stellar mass exceeds $10^4~M_\odot$. The disc fragments due to gravitational instability, creating a number of small clumps that rapidly migrate inward to fall onto the star. The resulting accretion history is thus highly time-dependent: short episodic accretion bursts are followed by longer, relative quiescent phases. We show that the circumstellar disc for the so-called direct collapse model is more unstable and generates greater variability over shorter timescales than normal Pop III cases. We conduct a post-process stellar evolution calculation using the obtained accretion history. Our results show that, regardless of the strong variability of the accretion rates, the stellar radius monotonically increases with almost constant effective temperature at $T_{\rm eff} \simeq 5000$ K as the stellar mass increases. The resulting UV feedback is too weak to hinder mass accretion due to the low flux of stellar UV photons, thus verifying our implicit assumption of no stellar feedback during the hydrodynamic simulations. The insensitivity of stellar evolution to variable accretion is attributed to the fact that typical timescales of variability, $\lesssim 10^3$ years, are too short to affect the stellar structure. We argue that this evolution will continue until the SMS eventually collapses to produce a massive black hole by the general relativistic instability after the stellar mass reaches $\gtrsim 10^5~M_\odot$.

preprint2014arXiv

Cross-Correlation of Cosmic Shear and Extragalactic Gamma-ray Background: Constraints on the Dark Matter Annihilation Cross-Section

We present the first measurement of the cross-correlation of weak gravitational lensing and the extragalactic gamma-ray background emission using data from the Canada-France-Hawaii Lensing Survey and the Fermi Large Area Telescope. The cross-correlation is a powerful probe of signatures of dark matter annihilation, because both cosmic shear and gamma-ray emission originate directly from the same DM distribution in the universe, and it can be used to derive constraints on dark matter annihilation cross-section. We show that the measured lensing-gamma correlation is consistent with a null signal. Comparing the result to theoretical predictions, we exclude dark matter annihilation cross sections of <sigma v> =10^{-24}-10^{-25} cm^3 s^-1 for a 100 GeV dark matter. If dark matter halos exist down to the mass scale of 10^-6 M_sun, we are able to place constraints on the thermal cross sections <sigma v> ~ 3 x 10^{-26} cm^3 s^-1 for a 10 GeV dark matter annihilation into tau^{+} tau^{-}. Future gravitational lensing surveys will increase sensitivity to probe annihilation cross sections of <sigma v> ~ 5 x 10^{-26} cm^3 s^-1 even for a 100 GeV dark matter. Detailed modeling of the contributions from astrophysical sources to the cross correlation signal could further improve the constraints by ~ 40-70 %.

preprint2014arXiv

Dust grain growth and the formation of the extremely primitive star SDSS J102915+172927

Dust grains in low-metallicity star-forming regions may be responsible for the formation of the first low-mass stars. The minimal conditions to activate dust-induced fragmentation require the gas to be pre-enriched above a critical dust-to-gas mass ratio Dcr=[2.6--6.3]x10^-9 with the spread reflecting the dependence on the grain properties. The recently discovered Galactic halo star SDSS J102915+172927 has a stellar mass of 0.8 Msun and a metallicity of Z=4.5x10^-5 Zsun and represents an optimal candidate for the dust-induced low-mass star formation. Indeed, for the two most plausible Population III supernova progenitors, with 20 Msun and 35 Msun, the critical dust-to-gas mass ratio can be overcome provided that at least 0.4 Msun of dust condenses in the ejecta, allowing for moderate destruction by the reverse shock. Here we show that even if dust formation in the first supernovae is less efficient or strong dust destruction does occur, grain growth during the collapse of the parent gas cloud is sufficiently rapid to activate dust cooling and likely fragmentation into low-mass and long-lived stars. Silicates and magnetite grains can experience significant grain growth in the density range 10^9 /cc < nH<10^12 /cc by accreting gas-phase species (SiO, SiO2, and Fe) until their gas-phase abundance drops to zero, reaching condensation efficiencies =1. The corresponding increase in the dust-to-gas mass ratio allows dust-induced cooling and fragmentation to be activated at 10^12 /cc < nH < 10^14 /cc, before the collapsing cloud becomes optically thick to continuum radiation. We show that for all the initial conditions that apply to the parent cloud of SDSS J102915+172927, dust-driven fragmentation is able to account for the formation of the star.

preprint2014arXiv

On the Systematic Errors of Cosmological-Scale Gravity Tests using Redshift Space Distortion: Non-linear Effects and the Halo Bias

Redshift space distortion (RSD) observed in galaxy redshift surveys is a powerful tool to test gravity theories on cosmological scales, but the systematic uncertainties must carefully be examined for future surveys with large statistics. Here we employ various analytic models of RSD and estimate the systematic errors on measurements of the structure growth-rate parameter, $fσ_8$, induced by non-linear effects and the halo bias with respect to the dark matter distribution, by using halo catalogues from 40 realisations of $3.4 \times 10^8$ comoving $h^{-3}$Mpc$^3$ cosmological N-body simulations. We consider hypothetical redshift surveys at redshifts z=0.5, 1.35 and 2, and different minimum halo mass thresholds in the range of $5.0 \times 10^{11}$ -- $2.0 \times 10^{13} h^{-1} M_\odot$. We find that the systematic error of $fσ_8$ is greatly reduced to ~5 per cent level, when a recently proposed analytical formula of RSD that takes into account the higher-order coupling between the density and velocity fields is adopted, with a scale-dependent parametric bias model. Dependence of the systematic error on the halo mass, the redshift, and the maximum wavenumber used in the analysis is discussed. We also find that the Wilson-Hilferty transformation is useful to improve the accuracy of likelihood analysis when only a small number of modes are available in power spectrum measurements.

preprint2014arXiv

One Hundred First Stars : Protostellar Evolution and the Final Masses

We perform a large set of radiation hydrodynamics simulations of primordial star formation in a fully cosmological context. Our statistical sample of 100 First Stars show that the first generation of stars have a wide mass distribution M_popIII = 10 ~ 1000 M_sun. We first run cosmological simulations to generate a set of primordial star-forming gas clouds. We then follow protostar formation in each gas cloud and the subsequent protostellar evolution until the gas mass accretion onto the protostar is halted by stellar radiative feedback. The accretion rates differ significantly among the primordial gas clouds which largely determine the final stellar masses. For low accretion rates the growth of a protostar is self-regulated by radiative feedback effects and the final mass is limited to several tens of solar masses. At high accretion rates the protostar's outer envelope continues to expand and the effective surface temperature remains low; such protostars do not exert strong radiative feedback and can grow in excess to one hundred solar masses. The obtained wide mass range suggests that the first stars play a variety of roles in the early universe, by triggering both core-collapse supernovae and pair-instability supernovae as well as by leaving stellar mass black holes. We find certain correlations between the final stellar mass and the physical properties of the star-forming cloud. These correlations can be used to estimate the mass of the first star from the properties of the parent cloud or of the host halo, without following the detailed protostellar evolution.

preprint2014arXiv

Physical Properties of UDF12 Galaxies in Cosmological simulations

We have performed a large cosmological hydrodynamics simulation tailored to the deep survey with the Hubble Space Telescope made in 2012, the so-called UDF12 campaign. After making a light-cone output, we have applied the same color selection criteria as the UDF12 campaign to select galaxies from our simulation, and then, have examined the physical properties of them as a proxy of the real observed UDF12 galaxies at $z > 7$. As a result, we find that the halo mass is almost linearly proportional to the observed ultraviolet (UV) luminosity ($4 \times 10^{11}~{\rm M_{\odot}}$ at $M_{\rm UV} = -21$). The dust attenuation and UV slope $β$ well correlates with the observed UV luminosity, which is consistent with observations quantitatively. The star formation rate (SFR) is also linearly proportional to the stellar mass and the specific SFR shows only a weak dependency on the mass. We also find an increasing star formation history with a time-scale of $\sim100$ Myr in the high-$z$ galaxies. An average metallicity weighted by the Lyman continuum luminosity reaches up to $>0.1$ Solar even at $z \sim 10$, suggesting a rapid metal enrichment. We also expect $\geq 0.1$ mJy at 350 GHz of the dust thermal emission from the galaxies with $H_{160} \leq 27$, which can be detectable with the Atacama Large Milimetre-submilimetre Array. The galaxies selected by the UDF12 survey contribute to only $52--12\%$ of the cosmic SFR density from $z \sim 7$ to $z \sim 10$, respectively. The James Webb Space Telescope will push the detection fraction up to $77--72\%$.

preprint2014arXiv

Reproducing cosmic evolution of galaxy population from $z = 4$ to $0$

We present cosmological hydrodynamic simulations performed to study evolution of galaxy population. The simulations follow timed release of mass, energy, and metals by stellar evolution and employ phenomenological treatments of supernova feedback, pre-supernova feedback modeled as feedback by radiation pressure from massive stars, and quenching of gas cooling in large halos. We construct the fiducial model so that it reproduces the observationally estimated galaxy stellar mass functions and the relationships between the galaxy stellar mass and the host halo mass from $z = 4$ to 0. We find that the fiducial model constructed this way naturally explains the cosmic star formation history, the galaxy downsizing, and the star formation rate and metallicity of the star-forming galaxies. The simulations without the quenching of the gas cooling in large halos overproduce massive galaxies at $z < 2$ and fail to reproduce galaxy downsizing. The simulations that do not employ the radiation pressure feedback from young stars predict too strong redshift evolution of the mass-metallicity relation. Furthermore, the slope of the relation becomes too steep at low redshift without the radiation pressure feedback. The metallicity dependence in the radiation pressure feedback is a key to explain the observed mass-metallicity relation. These facts indicate that these two processes in addition to supernova feedback are essential for galaxy evolution. Our simple phenomenological model is suitable to construct a mock galaxy sample to study physical properties of observed galaxy populations.

preprint2014arXiv

Statistical and Systematic Errors in Measurement of Weak-Lensing Minkowski Functionals: Application to Canada-France-Hawaii Lensing Survey

The measurement of cosmic shear using weak gravitational lensing is a challenging task that involves a number of complicated procedures. We study in detail the systematic errors in the measurement of weak lensing Minkowski Functionals (MFs). Specifically, we focus on systematics associated with galaxy shape measurements, photometric redshift errors, and shear calibration correction. We first generate mock weak-lensing catalogs that directly incorporate the actual observational characteristics of the Canada-France-Hawaii Lensing Survey (CFHTLenS). We then perform a Fisher analysis using the large set of mock catalogs for various cosmological models. We find that the statistical error associated with the observational effects degrades the cosmological parameter constraints by a factor of a few. The Subaru Hyper Suprime-Cam (HSC) survey with a sky coverage of ~1400 deg2 will constrain the dark energy equation of the state parameter with an error of Delta w_0 ~ 0.25 by the lensing MFs alone, but biases induced by the systematics can be comparable to the 1-sigma error. We conclude that the lensing MFs are powerful statistics beyond the two-point statistics, only if well-calibrated measurement of both the redshifts and the shapes of source galaxies is performed. Finally, we analyze the CFHTLenS data to explore the ability of the MFs to break degeneracies between a few cosmological parameters. Using a combined analysis of the MFs and the shear correlation function, we derive the matter density Omega_m0= 0.256+0.054-0.046.

preprint2014arXiv

Supernova dust formation and the grain growth in the early universe: The critical metallicity for low-mass star formation

We investigate the condition for the formation of low-mass second-generation stars in the early universe. It has been proposed that gas cooling by dust thermal emission can trigger fragmentation of a low-metallicity star-forming gas cloud. In order to determine the critical condition in which dust cooling induces the formation of low-mass stars, we follow the thermal evolution of a collapsing cloud by a one-zone semi-analytic collapse model. Earlier studies assume the dust amount in the local universe, where all refractory elements are depleted onto grains, and/or assume the constant dust amount during gas collapse. In this paper, we employ the models of dust formation and destruction in early supernovae to derive the realistic dust compositions and size distributions for multiple species as the initial conditions of our collapse calculations. We also follow accretion of heavy elements in the gas phase onto dust grains, i.e., grain growth, during gas contraction. We find that grain growth well alters the fragmentation property of the clouds, and that this still does not approach to the value in the local universe. The critical conditions can be written by the gas metallicity Zcr and the initial depletion efficiency fdep,0 of gas-phase metal onto grains, or dust-to-metal mass ratio, as (Zcr/10^{-5.5} Zsun) = (fdep,0/0.18)^{-0.44} with small scatters in the range of Zcr = [0.06--3.2]x10^{-5} Zsun. We also show that the initial dust composition and size distribution are important to determine Zcr.

preprint2014arXiv

The Dark Ages of the Universe and Hydrogen Reionization

One of the milestones in the cosmic history is the formation of the first luminous objects and Hydrogen reionization. The standard theory of cosmic structure formation predicts that the first generation of stars were born about a few hundred million years after the Big Bang. The dark Universe was then lit up once again, and eventually filled with ultraviolet photons emitted from stars, galaxies, and quasars. The exact epoch of the cosmic reionization and the details of the process, even the dominant sources, are not known except the fact that the universe was reionized early on. Signatures of reionization are expected to be imprinted in the cosmic microwave background radiation, especially in its large scale polarization. Future CMB experiments, together with other probes such as 21 cm surveys, will provide rich information on the process of reionization. We review recent studies on reionization. The implications from available observations in a wide range of wavelengths are discussed. Results from state-of-the-art computer simulations are presented. Finally, we discuss prospects for exploring the first few hundred million years of the cosmic history.

preprint2014arXiv

The origin of the most iron-poor star

We investigate the origin of carbon-enhanced metal-poor (CEMP) stars starting from the recently discovered $\rm [Fe/H]<-7.1$ star SMSS J031300 (Keller et al. 2014). We show that the elemental abundances observed on the surface of SMSS J031300 can be well fit by the yields of faint, metal free, supernovae. Using properly calibrated faint supernova explosion models, we study, for the first time, the formation of dust grains in such carbon-rich, iron-poor supernova ejecta. Calculations are performed assuming both unmixed and uniformly mixed ejecta and taking into account the partial destruction by the supernova reverse shock. We find that, due to the paucity of refractory elements beside carbon, amorphous carbon is the only grain species to form, with carbon condensation efficiencies that range between (0.15-0.84), resulting in dust yields in the range (0.025-2.25)M$_{\odot}$. We follow the collapse and fragmentation of a star forming cloud enriched by the products of these faint supernova explosions and we explore the role played by fine structure line cooling and dust cooling. We show that even if grain growth during the collapse has a minor effect of the dust-to-gas ratio, due to C depletion into CO molecules at an early stage of the collapse, the formation of CEMP low-mass stars, such as SMSS J031300, could be triggered by dust cooling and fragmentation. A comparison between model predictions and observations of a sample of C-normal and C-rich metal-poor stars supports the idea that a single common pathway may be responsible for the formation of the first low-mass stars.

preprint2013arXiv

ALMA will determine the spectroscopic redshift z>8 with FIR [O III] emission lines

We investigate the potential use of nebular emission lines in the rest-frame far-infrared (FIR) for determining spectroscopic redshift of z>8 galaxies with the Atacama Large Millimeter/sub-millimeter Array (ALMA). After making a line emissivity model as a function of metallicity, especially for the [O III] 88 micron line which is likely to be the strongest FIR line from H II regions, we predict the line fluxes from high-z galaxies based on a cosmological hydrodynamics simulation of galaxy formation. Since the metallicity of galaxies reaches at ~0.2 Zsun even at z>8 in our simulation, we expect the [O III] 88 micron line as strong as 1.3 mJy for 27 AB objects, which is detectable at a high significance by <1 hour integration with ALMA. Therefore, the [O III] 88 micron line would be the best tool to confirm the spectroscopic redshifts beyond z=8.

preprint2013arXiv

Constraints on Warm Dark Matter models from high-redshift long gamma-ray bursts

Structures in Warm Dark Matter (WDM) models are exponentially suppressed below a certain scale, characterized by the dark matter particle mass, $m_{\rm x}$. Since structures form hierarchically, the presence of collapsed objects at high-redshifts can set strong lower limits on $m_{\rm x}$. We place robust constraints on $m_{\rm x}$ using recent results from the {\it Swift} database of high-redshift gamma-ray bursts (GRBs). We parameterize the redshift evolution of the ratio between the cosmic GRB rate and star formation rate (SFR) as $\propto (1+z)^α$, thereby allowing astrophysical uncertainties to partially mimic the cosmological suppression of structures in WDM models. Using a maximum likelihood estimator on two different $z>4$ GRB subsamples (including two bursts at $z>8$), we constrain $m_{\rm x} \gtrsim 1.6$-1.8 keV at 95% C.L., when marginalized over a flat prior in $α$. We further estimate that 5 years of a SVOM-like mission would tighten these constraints to $m_{\rm x} \gtrsim 2.3 $ keV. Our results show that GRBs are a powerful probe of high-redshift structures, providing robust and competitive constraints on $m_{\rm x}$.

preprint2013arXiv

Detectability of High-Redshift Superluminous Supernovae with Upcoming Optical and Near-Infrared Surveys - II. Beyond z=6

Observational identification of the first stars is one of the great challenges in the modern astronomy. Although a single first star is too faint to be detected, supernova explosions of the first stars can be bright enough. An important question is whether such supernovae can be detected in the limited observational area with realistic observational resources. We perform detailed simulations to study the detectability of superluminous supernovae (SLSNe) at high redshifts, using the observationally-calibrated star formation rate density and supernova occurrence rate. We show that a 100 deg^2 survey with the limiting magnitude of 26 mag in near-infrared wavelengths will be able to discover about 10 SLSNe at z > 10. If the survey is extended to 200 deg^2 with 27 mag depth, about 10 SLSNe can be discovered at z > 15. We emphasize that the observations at > 3 um are crucial to detect and select SLSNe at z > 10. Our simulations are also applied to the planned survey with Euclid, WFIRST, and WISH. These surveys will be able to detect about 1000, 400, and 3000 SLSNe up to z ~ 5, 7, and 12, respectively. We conclude that detection of SLSNe at z > 10 is in fact achievable in the near future.

preprint2013arXiv

Effect of Mask Regions on Weak Lensing Statistics

Sky masking is unavoidable in wide-field weak lensing observations. We study how masks affect the measurement of statistics of matter distribution probed by weak gravitational lensing. We first use 1000 cosmological ray-tracing simulations to examine in detail the impact of masked regions on the weak lensing Minkowski Functionals (MFs). We consider actual sky masks used for a Subaru Suprime-Cam imaging survey. The masks increase the variance of the convergence field and the expected values of the MFs are biased. The bias then affects the non-Gaussian signals induced by the gravitational growth of structure. We then explore how masks affect cosmological parameter estimation. We calculate the cumulative signal-to-noise ratio S/N for masked maps to study the information content of lensing MFs. We show that the degradation of S/N for masked maps is mainly determined by the effective survey area. We also perform simple chi^2 analysis to show the impact of lensing MF bias due to masked regions. Finally, we compare ray-tracing simulations with data from a Subaru 2 deg^2 survey in order to address if the observed lensing MFs are consistent with those of the standard cosmology. The resulting chi^2/n_dof = 29.6/30 for three combined MFs, obtained with the mask effects taken into account, suggests that the observational data are indeed consistent with the standard LambdaCDM model. We conclude that the lensing MFs are powerful probe of cosmology only if mask effects are correctly taken into account.

preprint2013arXiv

Formation of Primordial Supermassive Stars by Rapid Mass Accretion

Supermassive stars (SMSs) forming via very rapid mass accretion (Mdot >~ 0.1 Msun/yr) could be precursors of supermassive black holes observed beyond redshift of about 6. Extending our previous work, we here study the evolution of primordial stars growing under such rapid mass accretion until the stellar mass reaches 10^{4 - 5} Msun. Our stellar evolution calculations show that a star becomes supermassive while passing through the "supergiant protostar" stage, whereby the star has a very bloated envelope and a contracting inner core. The stellar radius increases monotonically with the stellar mass, until =~ 100 AU for M_* >~ 10^4 Msun, after which the star begins to slowly contract. Because of the large radius the effective temperature is always less than 10^4 K during rapid accretion. The accreting material is thus almost completely transparent to the stellar radiation. Only for M_* >~ 10^5 Msun can stellar UV feedback operate and disturb the mass accretion flow. We also examine the pulsation stability of accreting SMSs, showing that the pulsation-driven mass loss does not prevent stellar mass growth. Observational signatures of bloated SMSs should be detectable with future observational facilities such as the James Webb Space Telescope. Our results predict that an inner core of the accreting SMS should suffer from the general relativistic instability soon after the stellar mass exceeds 10^5 Msun. An extremely massive black hole should form after the collapse of the inner core.

preprint2013arXiv

Growth of dust grains in a low-metallicity gas and its effect on the cloud fragmentation

In a low-metallicity gas, rapid cooling by dust thermal emission is considered to induce cloud fragmentation and play a vital role in the formation of low-mass stars (<~ 1 M_sun) in metal-poor environments. We investigate how the growth of dust grains through accretion of heavy elements in the gas phase onto grain surfaces alters the thermal evolution and fragmentation properties of a collapsing gas cloud. We calculate directly grain growth and dust emission cooling in a self-consistent manner. We show that MgSiO3 grains grow sufficiently at gas densities nH = 10^{10}, 10^{12}, and 10^{14} /cc for metallicities Z = 10^{-4}, 10^{-5}, and 10^{-6} Zsun, respectively, where the cooling of the collapsing gas cloud is enhanced. The condition for efficient dust cooling is insensitive to the initial condensation factor of pre-existing grains within the realistic range of 0.001--0.1, but sensitive to metallicity. The critical metallicity is Zcrit ~ 10^{-5.5} Zsun for the initial grain radius r_{MgSiO3,0} <~ 0.01 um and Zcrit ~ 10^{-4.5} Zsun for r_{MgSiO3,0} >~ 0.1 um. The formation of a recently discovered low-mass star with extremely low metallicity (<= 4.5x10^{-5} Zsun) could have been triggered by grain growth.

preprint2013arXiv

Modeling color-dependent galaxy clustering in cosmological simulations

We extend the subhalo abundance matching method to assign galaxy color to subhalos. We separate a luminosity-binned subhalo sample into two groups by a secondary subhalo property which is presumed to be correlated with galaxy color. The two subsamples then represent red and blue galaxy populations. We explore two models for the secondary property, namely subhalo assembly time and local dark matter density around each subhalo. The model predictions for the galaxy two-point correlation functions are compared with the recent results from the Sloan Digital Sky Survey. We show that the observed color dependence of galaxy clustering can be reproduced well by our method applied to cosmological N-body simulations without baryonic processes. We then compare the model predictions for the color-dependent galaxy-mass cross correlation functions with the results from gravitational lensing observations. The comparison allows us to distinguish the models, and also to discuss what subhalo property should be used to assign color to subhalos accurately. We show that the extended abundance matching method using the local dark matter density as a color proxy provides an accurate description of the galaxy populations in the local universe. We also study impacts of scatter in the local dark matter density-color relations. Introducing scatter improves agreements of our model predictions with the observed red and blue galaxy clustering and is needed to explain observed correlation functions in finer color bins. Finally, we study red galaxy fraction profiles in galaxy group- and cluster-sized halos and find the red fraction profiles have a relatively strong dependence on our model parameters. We argue that the red fraction profiles can be an important observational clue, in addition to galaxy clustering and galaxy-galaxy lensing, to explore the galaxy-(sub)halo connections.

preprint2013arXiv

On de-Sitter Geometry in Cosmic Void Statistics

Starting from the geometrical concept of a 4-dimensional de-Sitter configuration of spheres in Euclidean 3-space and modelling voids in the Universe as spheres, we show that a uniform distribution over this configuration space implies a power-law for the void number density which is consistent with results from the excursion set formalism and with data, for an intermediate range of void volumes. The scaling dimension of the large scale structure can be estimated as well. We also discuss the effect of restricting the survey geometry on the void statistics. This work is a new application of de-Sitter geometry to cosmology and also provides a new geometrical perspective on self-similarity in cosmology.

preprint2013arXiv

Structure of Dark Matter Halos in Warm Dark Matter models and in models with Long-Lived Charged Massive Particles

We study the formation of non-linear structures in Warm Dark Matter (WDM) models and in a Long-Lived Charged Massive Particle (CHAMP) model. CHAMPs with a decay lifetime of about 1 yr induce characteristic suppression in the matter power spectrum at subgalactic scales through acoustic oscillations in the thermal background. We explore structure formation in such a model. We also study three WDM models, where the dark matter particles are produced through the following mechanisms: i) WDM particles are produced in the thermal background and then kinematically decoupled; ii) WDM particles are fermions produced by the decay of thermal heavy bosons; and iii) WDM particles are produced by the decay of non-relativistic heavy particles. We show that the linear matter power spectra for the three models are all characterised by the comoving Jeans scale at the matter-radiation equality. Furthermore, we can also describe the linear matter power spectrum for the Long-Lived CHAMP model in terms of a suitably defined characteristic cut-off scale k_{Ch}, similarly to the WDM models. We perform large cosmological N-body simulations to study the non-linear growth of structures in these four models. We compare the halo mass functions, the subhalo mass functions, and the radial distributions of subhalos in simulated Milky Way-size halos. We study the models with k_{cut}=51, 410, 820 h/Mpc, and confirm that these statistics are indeed similar between the different WDM models and the Long-Lived CHAMP model. The result suggests that the cut-off scale k_{cut} not only characterises the linear power spectra but also can be used to predict the non-linear clustering properties. The radial distribution of subhalos in Milky Way-size halos is consistent with the observed distribution for k_{cut}~50-800 h/Mpc; such models resolve the so-called "missing satellite problem".

preprint2013arXiv

Weighing the Light Gravitino Mass with Weak Lensing Surveys

We explore the discovery potential of light gravitino mass m_{3/2} by combining future cosmology surveys and collider experiments. The former probe the imprint of light gravitinos in the cosmic matter density field, whereas the latter search signatures of a supersymmetry breaking mechanism. Free-streaming of light gravitinos suppresses the density fluctuations at galactic and sub-galactic length scales, where weak gravitational lensing can be used as a powerful probe. We perform numerical simulations of structure formation to quantify the effect. We then run realistic ray-tracing simulations of gravitational lensing to measure the cosmic shear in models with light gravitino. We forecast the possible reach of future wide-field surveys by Fisher analysis; the light gravitino mass can be determined with an accuracy of m_{3/2}=4\pm 1 eV by a combination of the Hyper Suprime Cam survey and cosmic microwave background anisotropy data obtained by Planck satellite. The corresponding accuracy to be obtained by the future Large Synoptic Survey Telescope is δm_{3/2}=0.6 eV. Data from experiments at Large Hadron Collider at 14 TeV will provide constraint at m_{3/2} \simeq 5 eV in the minimal framework of gauge-mediated supersymmetry breaking (GMSB) model. We conclude that a large class of the GMSB model can be tested by combining the cosmological observations and the collider experiments.

preprint2012arXiv

Detectability of High-Redshift Superluminous Supernovae with Upcoming Optical and Near-Infrared Surveys

Observations of high-redshift supernovae (SNe) open a novel opportunity to study the massive star population in the early Universe. We study the detectability of superluminous SNe with upcoming optical and near-infrared (NIR) surveys. Our calculations are based on the cosmic star formation history, the SN occurence rate, the characteristic colour and the light curve of the SNe that are all calibrated by available observations. We show that 15-150 SNe up to z ~ 4 will be discovered by the proposed Subaru/Hyper Suprime-Cam deep survey: 30 deg^2 survey with 24.5 AB mag depth in z-band for 3 months. With its ultra-deep layer (3.5 deg^2 with 25.6 AB mag depth in z-band for 4 months), the highest redshift can be extended to z ~ 5. We further explore the detectability by upcoming NIR survey utilizing future satellites such as Euclid, WFIRST, and WISH. The wide-field NIR surveys are very efficient to detect high-redshift SNe. With a hypothetical deep NIR survey for 100 deg^2 with 26 AB mag depth at 1-4 um, at least ~ 50 SNe will be discovered at z>3 in half a year. The number of the detected SNe can place a strong constraint on the stellar initial mass function or its slope especially at the high-mass end. Superluminous SNe at high redshifts can be distinguished from other types of SNe by the long time-scale of their light curves in the observer's frame, the optical colours redder than other core-collapse SNe and the NIR colours redder than any other types of SNe.

preprint2012arXiv

Direct Integration of the Collisionless Boltzmann Equation in Six-dimensional Phase Space: Self-gravitating Systems

We present a scheme for numerical simulations of collisionless self-gravitating systems which directly integrates the Vlasov--Poisson equations in six-dimensional phase space. By the results from a suite of large-scale numerical simulations, we demonstrate that the present scheme can simulate collisionless self-gravitating systems properly. The integration scheme is based on the positive flux conservation method recently developed in plasma physics. We test the accuracy of our code by performing several test calculations including the stability of King spheres, the gravitational instability and the Landau damping. We show that the mass and the energy are accurately conserved for all the test cases we study. The results are in good agreement with linear theory predictions and/or analytic solutions. The distribution function keeps the property of positivity and remains non-oscillatory. The largest simulations are run on 64^6 grids. The computation speed scales well with the number of processors, and thus our code performs efficiently on massively parallel supercomputers.

preprint2012arXiv

Distribution of Dust around Galaxies: An Analytic Model

We develop an analytic halo model for the distribution of dust around galaxies. The model results are compared with the observed surface dust density profile measured through reddening of background quasars in the Sloan Digital Sky Survey (SDSS) reported by Menard et al.(2010). We assume that the dust distribution around a galaxy is described by a simple power law, similarly to the mass distribution, but with a sharp cut-off at $αR_{\rm vir}$ where $R_{\rm vir}$ is the galaxy's virial radius and $α$ is a model parameter. Our model reproduces the observed dust distribution profile very well over a wide range of radial distance of $10 - 10^{4} h^{-1}$kpc. For the characteristic galaxy halo mass of $2\times 10^{12} h^{-1}M_{\odot}$ estimated for the SDSS galaxies, the best fit model is obtained if $α$ is greater than unity, which suggests that dust is distributed to over a few hundred kilo-parsecs from the galaxies. The observed large-scale dust distribution profile is reproduced if we assume the total amount of dust is equal to that estimated from the integrated stellar evolution over the cosmic time.

preprint2012arXiv

Formation and evolution of primordial protostellar systems

We investigate the formation of the first stars at the end of the cosmic dark ages with a suite of three-dimensional, moving mesh simulations that directly resolve the collapse of the gas beyond the formation of the first protostar at the centre of a dark matter minihalo. The simulations cover more than 25 orders of magnitude in density and have a maximum spatial resolution of 0.05 R_sun, which extends well below the radius of individual protostars and captures their interaction with the surrounding gas. In analogy to previous studies that employed sink particles, we find that the Keplerian disc around the primary protostar fragments into a number of secondary protostars, which is facilitated by H2 collisional dissociation cooling and collision-induced emission. The further evolution of the protostellar system is characterized by strong gravitational torques that transfer angular momentum between the secondary protostars formed in the disc and the surrounding gas. This leads to the migration of about half of the secondary protostars to the centre of the cloud in a free-fall time, where they merge with the primary protostar and enhance its growth to about five times the mass of the second most massive protostar. By the same token, a fraction of the protostars obtain angular momentum from other protostars via N-body interactions and migrate to higher orbits. On average, only every third protostar survives until the end of the simulation. However, the number of protostars present at any given time increases monotonically, suggesting that the system will continue to grow beyond the limited period of time simulated here.

preprint2012arXiv

Light Curve Modeling of Superluminous Supernova 2006gy: Collision between Supernova Ejecta and Dense Circumstellar Medium

We show model light curves of superluminous supernova 2006gy on the assumption that the supernova is powered by the collision of supernova ejecta and its dense circumstellar medium. The initial conditions are constructed based on the shock breakout condition, assuming that the circumstellar medium is dense enough to cause the shock breakout within it. We perform a set of numerical light curve calculations by using a one-dimensional multigroup radiation hydrodynamics code STELLA. We succeeded in reproducing the overall features of the early light curve of SN 2006gy with the circumstellar medium whose mass is about 15 Msun (the average mass-loss rate ~ 0.1 Msun/yr). Thus, the progenitor of SN 2006gy is likely a very massive star. The density profile of the circumstellar medium is not well constrained by the light curve modeling only, but our modeling disfavors the circumstellar medium formed by steady mass loss. The ejecta mass is estimated to be comparable to or less than 15 Msun and the explosion energy is expected to be more than 4e51 erg. No 56Ni is required to explain the early light curve. We find that the multidimensional effect, e.g., the Rayleigh-Taylor instability, which is expected to take place in the cool dense shell between the supernova ejecta and the dense circumstellar medium, is important in understanding supernovae powered by the shock interaction. We also show the evolution of the optical and near-infrared model light curves of high-redshift superluminous supernovae. They can be potentially used to identify SN 2006gy-like superluminous supernovae in the future optical and near-infrared transient surveys.

preprint2012arXiv

Low-mass star formation triggered by early supernova explosions

We study the formation of low-mass and extremely metal-poor stars in the early universe. Our study is motivated by the recent discovery of a low-mass (M < 0.8 Msun) and extremely metal-poor (Z <= 4.5 x 10^{-5} Zsun) star in the Galactic halo by Caffau et al. We propose a model that early supernova (SN) explosions trigger the formation of low-mass stars via shell fragmentation. We first perform one-dimensional hydrodynamic simulations of the evolution of an early SN remnant. We show that the shocked shell undergoes efficient radiative cooling and then becomes gravitationally unstable to fragment and collapse in about ten million years. We then follow the thermal evolution of the collapsing fragments using a one-zone code. Our one-zone calculation treats chemistry and radiative cooling self-consistently in low-metallicity gas. The collapsing gas cloud evolves roughly isothermally, until it cools rapidly by dust continuum emission at the density 10^{13}-10^{14} /cc. The cloud core then becomes thermally and gravitationally unstable and fragments. We argue that early SNe can trigger the formation of low-mass stars in the extremely metal-poor environment as Caffau et al. discovered recently.

preprint2012arXiv

Probing Primordial Non-Gaussianity with Weak Lensing Minkowski Functionals

We study the cosmological information contained in the Minkowski Functionals (MFs) of weak gravitational lensing convergence maps. We show that the MFs provide strong constraints on the local type primordial non-Gaussianity parameter f_NL. We run a set of cosmological N-body simulations and perform ray-tracing simulations of weak lensing, to generate 100 independent convergence maps of 25 deg^2 field-of-view for f_NL = -100, 0 and 100. We perform a Fisher analysis to study the degeneracy among other cosmological parameters such as the dark energy equation of state parameter w and the fluctuation amplitude sigma_8. We use fully nonlinear covariance matrices evaluated from 1000 ray-tracing simulations. For the upcoming wide-field observations such as Subaru Hyper Suprime-Cam survey with the proposed survey area of 1500 deg^2, the primordial non-Gaussianity can be constrained with a level of f_NL ~ 80 and w ~ 0.036 by weak lensing MFs. If simply scaled by the effective survey area, a 20000 deg^2 lensing survey using Large Synoptic Survey Telescope will give constraints of f_NL ~ 25 and w ~ 0.013. We show that these constraints can be further improved by a tomographic method using source galaxies in multiple redshift bins.

preprint2012arXiv

Protostellar Feedback and Final Mass of the Second-Generation Primordial Stars

The first stars in the universe ionized the ambient primordial gas through various feedback processes. "Second-generation" primordial stars potentially form from this disturbed gas after its recombination. In this Letter, we study the late formation stage of such second-generation stars, where a large amount of gas accretes onto the protostar and the final stellar mass is determined when the accretion terminates. We directly compute the complex interplay between the accretion flow and stellar ultraviolet (UV) radiation, performing radiation-hydrodynamic simulations coupled with stellar evolution calculations. Because of more efficient H2 and HD cooling in the pre-stellar stage, the accretion rates onto the star are ten times lower than in the case of the formation of the first stars. The lower accretion rates and envelope density result in the occurrence of an expanding bipolar HII region at a lower protostellar mass M_* \simeq 10Msun, which blows out the circumstellar material, thereby quenching the mass supply from the envelope to the accretion disk. At the same time the disk loses mass due to photoevaporation by the growing star. In our fiducial case the stellar UV feedback terminates mass accretion onto the star at M_* \simeq 17Msun. Although the derived masses of the second-generation primordial stars are systematically lower than those of the first generation, the difference is within a factor of only a few. Our results suggest a new scenario, whereby the majority of the primordial stars are born as massive stars with tens of solar masses, regardless of their generations.

preprint2012arXiv

Radiative cooling implementations in simulations of primordial star formation

We study the thermal evolution of primordial star-forming gas clouds using three-dimensional cosmological simulations. We critically examine how assumptions and approximations made in calculating radiative cooling rates affect the dynamics of the collapsing gas clouds. We consider two important molecular hydrogen cooling processes that operate in a dense primordial gas; H_2 line cooling and continuum cooling by H_2 collision-induced emission. To calculate the optically thick cooling rates, we follow the Sobolev method for the former, whereas we perform ray-tracing for the latter. We also run the same set of simulations using simplified fitting functions for the net cooling rates. We compare the simulation results in detail. We show that the time- and direction-dependence of hydrodynamic quantities such as gas temperature and local velocity gradients significantly affects the optically thick cooling rates. Gravitational collapse of the cloud core is accelerated when the cooling rates are calculated by using the fitting functions. The structure and evolution of the central pre-stellar disk are also affected. We conclude that physically motivated implementations of radiative transfer are necessary to follow accurately the thermal and chemical evolution of a primordial gas to high densities.

preprint2012arXiv

Simulations of Early Baryonic Structure Formation with Stream Velocity: I. Halo Abundance

It has been recently shown that the relative velocity between the dark matter and the baryons (vbc) at the time of recombination can affect the structure formation in the early universe (Tseliakhovich & Hirata 2010). We statistically quantify this effect using large cosmological simulations. We use three different high resolution sets of simulations (with separate transfer functions for baryons and dark matter) that vary in box size, particle number, and the value of the relative velocity between dark matter and baryons. We show that the total number density of halos is suppressed by ~ 20% at z = 25 for vbc = 1σ(vbc), where σ(vbc) is the variance of the relative velocity, while for vbc = 3.4σ(vbc) the relative suppression at the same redshift reaches 50%, remaining at or above the 30% level all the way to z = 11. We also find high abundance of "empty halos", i.e., halos that have gas fraction below half of the cosmic mean baryonic fraction fb. Specifically we find that for vbc = 1σ(vbc) all halos below 10^5M\odot are empty at z \geq 19. The high abundance of empty halos results in significant delay in the formation of gas rich mini-halos and the first galaxies.

preprint2012arXiv

Simulations of Early Baryonic Structure Formation with Stream Velocity: II. The Gas Fraction

Understanding the gas content of high redshift halos is crucial for studying the formation of the first generation of galaxies and reionization. Recently, Tseliakhovich & Hirata showed that the relative "stream" velocity between the dark matter and baryons at the time of recombination - formally a second order effect, but an unusually large one - can influence the later structure formation history of the Universe. We quantify the effect of the stream velocity on the so-called "characteristic mass" - the minimum mass of a dark matter halo capable of retaining most of its baryons throughout its formation epoch - using three different high-resolution sets of cosmological simulations (with separate transfer functions for baryons and dark matter) that vary in box size, particle number, and the value of the relative velocity between the dark matter and baryons. In order to understand this effect theoretically, we generalize the linear theory filtering mass to properly account for the difference between the dark matter and baryonic density fluctuation evolution induced by the stream velocity. We show that the new filtering mass provides an accurate estimate for the characteristic mass, while other theoretical ansatzes for the characteristic mass are substantially less precise.

preprint2012arXiv

Simulations of Wide-Field Weak Lensing Surveys I: Basic Statistics and Non-Gaussian Effects

We study the lensing convergence power spectrum and its covariance for a standard LCDM cosmology. We run 400 cosmological N-body simulations and use the outputs to perform a total of 1000 independent ray-tracing simulations. We compare the simulation results with analytic model predictions. The semi-analytic model based on Smith et al.(2003) fitting formula underestimates the convergence power by ~30% at arc-minute angular scales. For the convergence power spectrum covariance, the halo model reproduces the simulation results remarkably well over a wide range of angular scales and source redshifts. The dominant contribution at small angular scales comes from the sample variance due to the number fluctuations of halos in a finite survey volume. The signal-to-noise ratio for the convergence power spectrum is degraded by the non-Gaussian covariances by up to a factor 5 for a weak lensing survey to z_s ~1. The probability distribution of the convergence power spectrum estimators, among the realizations, is well approximated by a chi-square distribution with broadened variance given by the non-Gaussian covariance, but has a larger positive tail. The skewness and kurtosis have non-negligible values especially for a shallow survey. We argue that a prior knowledge on the full distribution may be needed to obtain an unbiased estimate on the ensemble averaged band power at each angular scale from a finite volume survey.

preprint2012arXiv

Sub-millimetre galaxies in cosmological hydrodynamic simulations: Source number counts and the spatial clustering

We use large cosmological Smoothed-Particle-Hydrodynamics simulations to study the formation and evolution of sub-millimetre galaxies (SMGs). In our previous work, we studied the statistical properties of ultra-violet selected star-forming galaxies at high redshifts. We populate the same cosmological simulations with SMGs by calculating the reprocess of stellar light by dust grains into far-infrared to millimetre wavebands in a self-consistent manner. We generate light-cone outputs to compare directly the statistical properties of the simulated SMGs with available observations. Our model reproduces the submm source number counts and the clustering amplitude. We show that bright SMGs with flux $S > 1$ mJy reside in halos with mass of $\sim 10^{13} M_{\odot}$ and have stellar masses greater than $10^{11}\sim \rm M_{\odot}$. The angular cross-correlation between the SMGs and Lyman-$α$ emitters is significantly weaker than that between the SMGs and Lyman-break galaxies. The cross-correlation is also weaker than the auto-correlation of the SMGs. The redshift distribution of the SMGs shows a broad peak at $z \sim 2$, where Bright SMGs contribute significantly to the global cosmic star formation rate density. Our model predicts that there are hundreds of SMGs with $S > 0.1$ mJy at $z > 5$ per 1 square degree field. Such SMGs can be detected by ALMA.

preprint2012arXiv

The nature of dark matter from the global high redshift HI 21 cm signal

We study the imprint of dark matter (DM) annihilation on the global 21 cm signal from the Dark Ages to Cosmic Reionization. Motivated by recent observations, we focus on three DM candidates: (i) a 10 GeV Bino-like neutralino (ii) a 200 GeV Wino and (iii) a 1 TeV heavier particle annihilating into leptons. For each DM candidate we assume two values for the thermally averaged annihilation cross section σv, the standard thermal value σv_th = 3 x 10^-26 cm^3 s^-1 and the maximum value allowed by WMAP7 data, σv_max. We include the enhancement of DM annihilations due to collapsed structures, detailed estimates of energy deposition into the intergalactic medium (IGM), as well realistic prescriptions for astrophysical sources of UV and X-ray radiation. In these models, the additional heat input from DM annihilation suppresses the mean 21cm brightness temperature offset by δTb ~ few-100 mK. In particular, the very deep δTb ~ -150 mK absorption feature at ~20 < z < 25 predicted by popular models of the first galaxies is considerably reduced or totally erased by some of the considered DM candidates. Such an enhancement in IGM heating could come from either DM annihilations or a stronger-than-expected astrophysical component (i.e. abundant early X-ray sources). However, we find that the two signatures are not degenerate, since the DM heating is dominated by halos several orders of magnitude smaller than those hosting galaxies, whose fractional abundance evolves more slowly resulting in a smaller gradient: d δTb/d ν< 4 mK/MHz in the range ν~60-80 MHz. The detection of such signals by future radio telescopes would be clear evidence of DM energy injection at high-redshifts.

preprint2012arXiv

WIMP DM and first stars: suppression of fragmentation in primordial star formation

We present the first 3D simulations to include the effects of dark matter annihilation feedback during the collapse of primordial mini-halos. We begin our simulations from cosmological initial conditions and account for dark matter annihilation in our treatment of the chemical and thermal evolution of the gas. The dark matter is modelled using an analytical density profile that responds to changes in the peak gas density. We find that the gas can collapse to high densities despite the additional energy input from the dark matter. No objects supported purely by dark matter annihilation heating are formed in our simulations. However, we find that the dark matter annihilation heating has a large effect on the evolution of the gas following the formation of the first protostar. Previous simulations without dark matter annihilation found that protostellar discs around Population III stars rapidly fragmented, forming multiple protostars that underwent mergers or ejections. When dark matter annihilation is included, however, these discs become stable to radii of 1000 AU or more. In the cases where fragmentation does occur, it is a wide binary that is formed.

preprint2011arXiv

Evolution of Primordial Stars Powered by Dark Matter Annihilation up to the Main-Sequence Stage

Primordial stars formed in the early universe are thought to be hosted by compact dark matter (DM) halos. If DM consists of Weakly Interacting Massive Particles (WIMPs), such stars may be powered by DM annihilation during the early phases of their evolutions. We study the pre-main sequence evolutions of the primordial star using a detailed stellar evolution code under the assumption that the annihilation of adiabatically contracted WIMPs DM within the star provides a sufficient energy to sustain the stellar equilibrium. We follow the evolution of accreting stars using several gas mass accretion rates derived from cosmological simulations. We show that the stellar mass becomes very large, up to 900 - 1000 M_sun when the star reaches the main-sequence phase for a reasonable set of model parameters such as DM particle mass and the annihilation cross section. During the dark star phase, the star expands over a thousand solar-radii, while the surface temperature remains below 10^4 K. The energy generated by nuclear reactions is not dominant during this phase. We also study models with different gas mass accretion rates and the DM particle masses. All our models for different DM particle masses pass the dark star phase. The final mass of the dark stars is essentially unchanged for DM mass of m_DM <= 10 GeV. Gravitational collapse of the massive dark stars will leave massive black holes with mass as large as 1000 M_sun in the early universe.

preprint2011arXiv

Lyman-alpha Emitters in Cosmological Simulations I: Lyman-alpha Escape Fraction and Statistical Properties at z=3.1

We use very large cosmological Smoothed-Particle-Hydrodynamics simulations to study the properties of high redshift Lyman-alpha emitters (LAEs). We identify star-forming galaxies at z=3.1 in a cosmological volume of 100 h^-1Mpc on a side. We develop a phenomenological model of absorption, scattering and escape of Lyman-alpha photons on the assumption that the clumpiness of the inter-stellar medium in a galaxy is correlated with the larger scale substructure richness. The radiative transfer effect proposed by Neufeld (1991) allows a large fraction of Lyman-alpha photons to escape from a clumpy galaxy even if it contains a substantial amount of dust. Our model reproduces, for the first time, all of the following observed properties of LAEs at z=3.1: the angular correlation function, ultra-violet and Lyman-alpha luminosity functions, and the equivalent width distribution. A simple model that takes only dust absorption into account fails in matching the observational data, suggesting that the kind of effect we consider is needed. Our model also predicts a bimodal age distribution for LAEs. There are old, massive and dusty LAEs, similar to recently found high redshift LAEs. The large LAEs have escape fractions of Lyman-alpha photons of f_esc ~ 0.05-0.1.

preprint2011arXiv

Matter Distribution around Galaxies

We explore the mass distribution of material associated with galaxies from the observation of gravitational weak lensing for the galaxy mass correlation function with the aid of $N$-body simulations of dark matter. The latter is employed to unfold various contributions that contribute to the integrated line of sight mass density. We conclude that galaxies have no definite edges of the matter distribution, extending to the middle to neighbouring galaxies with the density profile roughly $r^{-2.4}$ beyond the virial radius. The mass distributed beyond the virial radius (gravitationally bound radius) explains the gap seen in the mass density estimates, the global value $Ω_m\sim 0.27$ and typically $Ω_{\rm gal} \sim 0.15$ from the luminosity density multiplied by the mass to light ratio. We suggest to use a physical method of gravitational lensing to characterise galaxy samples rather than characterise them with photometric means.

preprint2011arXiv

Origin of multiple nuclei in ultraluminous infrared galaxies

Ultraluminous infrared galaxies (ULIRGs) with multiple ($\ge 3$) nuclei are frequently observed. It has been suggested that these nuclei are produced by multiple major mergers of galaxies. The expected rate of such mergers is, however, too low to reproduce the observed number of ULIRGs with multiple nuclei. We have performed high-resolution simulations of the merging of two gas-rich disk galaxies. We found that extremely massive and compact star clusters form from the strongly disturbed gas disks after the first or second encounter between the galaxies. The mass of such clusters reaches $\sim 10^8 M_{\odot}$, and their half-mass radii are $20-30 \rm{pc}$. Since these clusters consist of young stars, they appear to be several bright cores in the galactic central region ($\sim \rm{kpc}$). The peak luminosity of these clusters reaches $\sim 10%$ of the total luminosity of the merging galaxy. These massive and compact clusters are consistent with the characteristics of the observed multiple nuclei in ULIRGs. Multiple mergers are not necessary to explain multiple nuclei in ULIRGs.

preprint2011arXiv

Population III.1 and III.2 gamma-Ray Bursts: Constraints on the event rate for future radio and X-ray surveys

We calculate the theoretical event rate of gamma-ray bursts (GRBs) from the collapse of massive first-generation (Population III; Pop III) stars. The Pop III GRBs could be super-energetic with the isotropic energy up to $E_{\rm iso} \gtrsim 10^{55-57}$ ergs, providing a unique probe of the high-redshift Universe. We consider both the so-called Pop III.1 stars (primordial) and Pop III.2 stars (primordial but affected by radiation from other stars). We employ a semi-analytical approach that considers inhomogeneous hydrogen reionization and chemical evolution of the intergalactic medium. We show that Pop III.2 GRBs occur more than 100 times more frequently than Pop III.1 GRBs, and thus should be suitable targets for future GRB missions. Interestingly, our optimistic model predicts an event rate that is already constrained by the current radio transient searches. We expect $\sim10-10^4$ radio afterglows above $\sim$ 0.3 mJy on the sky with $\sim 1$ year variability and mostly without GRBs (orphans), which are detectable by ALMA, EVLA, LOFAR, and SKA, while we expect to observe maximum of $N < 20$ GRBs per year integrated over at $z>6$ for Pop III.2 and $N < 0.08$ per year integrated over at $z>10$ for Pop III.1 with EXIST, and $N < 0.2$ for Pop III.2 GRBs per year integrated over at $z > 6$ with \textit{Swift}.

preprint2011arXiv

Protostellar Feedback Halts the Growth of the First Stars in the Universe

The first stars fundamentally transformed the early universe by emitting the first light and by producing the first heavy elements. These effects were predetermined by the mass distribution of the first stars, which is thought to have been fixed by a complex interplay of gas accretion and protostellar radiation. We performed radiation-hydrodynamics simulations that followed the growth of a primordial protostar through to the early stages as a star with thermo-nuclear burning. The circumstellar accretion disk was evaporated by ultraviolet radiation from the star when its mass was 43 times that of the Sun. Such massive primordial stars, in contrast to the often postulated extremely massive stars, may help explain the fact that there are no signatures of the pair-instability supernovae in abundance patterns of metal-poor stars in our galaxy.

preprint2011arXiv

Shock-induced star cluster formation in colliding galaxies

We studied the formation process of star clusters using high-resolution N-body/smoothed particle hydrodynamcs simulations of colliding galaxies. The total number of particles is 1.2x10^8 for our high resolution run. The gravitational softening is 5 pc and we allow gas to cool down to \sim 10 K. During the first encounter of the collision, a giant filament consists of cold and dense gas found between the progenitors by shock compression. A vigorous starburst took place in the filament, resulting in the formation of star clusters. The mass of these star clusters ranges from 10^{5-8} Msun. These star clusters formed hierarchically: at first small star clusters formed, and then they merged via gravity, resulting in larger star clusters.

preprint2011arXiv

The First Galaxies

We review our current understanding of how the first galaxies formed at the end of the cosmic dark ages, a few 100 million years after the Big Bang. Modern large telescopes discovered galaxies at redshifts greater than seven, whereas theoretical studies have just reached the degree of sophistication necessary to make meaningful predictions. A crucial ingredient is the feedback exerted by the first generation of stars, through UV radiation, supernova blast waves, and chemical enrichment. The key goal is to derive the signature of the first galaxies to be observed with upcoming or planned next-generation facilities, such as the James Webb Space Telescope or Atacama Large Millimeter Array. From the observational side, ongoing deep-field searches for very high-redshift galaxies begin to provide us with empirical constraints on the nature of the first galaxies.

preprint2011arXiv

The nonlinear evolution of baryonic overdensities in the early universe: Initial conditions of numerical simulations

We run very large cosmological N-body hydrodynamical simulations in order to study statistically the baryon fractions in early dark matter halos. We critically examine how differences in the initial conditions affect the gas fraction in the redshift range z = 11--21. We test three different linear power spectra for the initial conditions: (1) A complete heating model, which is our fiducial model; this model follows the evolution of overdensities correctly, according to Naoz & Barkana (2005), in particular including the spatial variation of the speed of sound of the gas due to Compton heating from the CMB. (2) An equal-δ model, which assumes that the initial baryon fluctuations are equal to those of the dark matter, while conserving sigma8 of the total matter. (3) A mean cs model, which assumes a uniform speed of sound of the gas. The latter two models are often used in the literature. We calculate the baryon fractions for a large sample of halos in our simulations. Our fiducial model implies that before reionization and significant stellar heating took place, the minimum mass needed for a minihalo to keep most of its baryons throughout its formation was ~ 3 * 10^4 Msun. However, the alternative models yield a wrong (higher by about 50%) minimum mass, since the system retains a memory of the initial conditions. We also demonstrate this using the "filtering mass" from linear theory, which accurately describes the evolution of the baryon fraction throughout the simulated redshift range.

preprint2011arXiv

The pairwise velocity probability density function in models with local primordial non-Gaussianity

We study how primordial non-Gaussianity affects the pairwise velocity probability density function (PDF) using an analytical model and cosmological N-body simulations. We adopt the local type non-Gaussian models characterized by f_{nl}, and examine both the linear velocity difference PDF and the linear pairwise velocity PDF. We show explicitly how f_{nl} induces correlations between originally independent velocities along the parallel and the perpendicular to the line of separation directions. We compare the model results with measurements from N-body simulations of the non-Gaussian models. Linear theory fails to predict the PDF in the f_{nl} models. Therefore we develop an analytic model based on the Zeldovich approximation to describe the evolution of the velocity PDF. Our analytical model and simulation results show remarkably good agreement in both the parallel and the perpendicular directions for the PDF profiles, as well as the change in the PDF due to primordial non-Gaussianity. The agreement is particularly good for relatively small separations (< 10 h^{-1} Mpc). The inclusion of the evolution of the velocity PDF is important to obtain a good description on the signature of primordial non-Gaussianity in the PDF. Our model provides the foundation to constrain f_{nl} using the peculiar velocity in future surveys.

preprint2010arXiv

3D simulations of the accretion process in Kerr space-time with arbitrary value of the spin parameter

We present the results of three-dimensional general relativistic hydrodynamic simulations of adiabatic and spherically symmetric accretion in Kerr space-time. We consider compact objects with spin parameter $|a_*| \le 1$ (black holes) and with $|a_*| > 1$ (super-spinars). Our full three-dimensional simulations confirm the formation of equatorial outflows for high values of $|a_*|$, as found in our previous work in 2.5 dimensions. We show that the critical value of $|a_*|$ determining the onset of powerful outflows depends mainly on the radius of the compact object. The phenomenon of equatorial outflows can hardly occur around a black hole and may thus be used to test the bound $|a_*| \le 1$ for astrophysical black hole candidates.

preprint2010arXiv

Low-Metallicity Star Formation : Prestellar Collapse and Protostellar Accretion in the Spherical Symmetry

The collapse of dense cores with different metallicities is studied by hydrodynamical calculations coupled with detailed chemical and radiative processes. For this purpose, we construct a simple chemical network with non-equilibrium reactions among 15 chemical species, which reproduces the abundance of important molecular coolants by more detailed network very well. The evolution is followed until the formation of a hydrostatic protostar at the center. In a lower-metallicity gas cloud, the temperature during the collapse remains high owing to less efficient cooling. Using the temperature evolution at the center as a function the density, we discuss the possibility of fragmentation during the dust-cooling phase. The critical metallicity for the fragmentation is 10^{-5}Z_sun assuming moderate elongation of the cloud cores at the onset of this phase. From the density and velocity distributions at the time of protostar formation, we evaluate the mass accretion rate in the subsequent accretion phase. Using these accretion rates, we also calculate the evolution of the protostars under the assumption of stationary accretion flow. Finally, we discuss possible suppression of fragmentation by heating of the ambient gas by protostellar radiation, which is considered important in the contemporary star formation. We argue that it is negligible for <10^{-2}Zsun.

preprint2010arXiv

Non-Gaussian Error Contribution to Likelihood Analysis of the Matter Power Spectrum

We study the sample variance of the matter power spectrum for the standard Lambda Cold Dark Matter universe. We use a total of 5000 cosmological N-body simulations to study in detail the distribution of best-fit cosmological parameters and the baryon acoustic peak positions. The obtained distribution is compared with the results from the Fisher matrix analysis with and without including non-Gaussian errors. For the Fisher matrix analysis, we compute the derivatives of the matter power spectrum with respect to cosmological parameters using directly full nonlinear simulations. We show that the non-Gaussian errors increase the unmarginalized errors by up to a factor 5 for k_{max}=0.4h/Mpc if there is only one free parameter provided other parameters are well determined by external information. On the other hand, for multi-parameter fitting, the impact of the non-Gaussian errors is significantly mitigated due to severe parameter degeneracies in the power spectrum. The distribution of the acoustic peak positions is well described by a Gaussian distribution, with its width being consistent with the statistical interval predicted from the Fisher matrix. We also examine systematic bias in the best-fit parameter due to the non-Gaussian errors. The bias is found to be smaller than the 1 sigma statistical error for both the cosmological parameters and the acoustic scale positions.

preprint2010arXiv

Outflows from accreting super-spinars

In this paper we continue our study on the accretion process onto super-spinning Kerr objects with no event horizon (super-spinars). We discuss the counterpart of the Bondi accretion onto black holes. We first report the results of our numerical simulations. We found a quasi steady-state configuration for any choice of the parameters of our model. The most interesting feature is the presence of hot outflows. Unlike jets and outflows produced around black holes, which are thought to be powered by magnetic fields and emitted from the poles, here the outflows are produced by the repulsive gravitational force at a small distance from the super-spinar and are ejected around the equatorial plane. In some circumstances, the amount of matter in the outflow is considerable, which can indeed significantly reduce the gas mass accretion rate. Finally, we discuss a possible scenario of the accretion process in more realistic situations, which cannot be simulated by our code.

preprint2010arXiv

Shape and position of the shadow in the $δ= 2$ Tomimatsu-Sato space-time

Within 5-10 years, very long baseline interferometry facilities will be able to observe the "shadow" of super-massive black hole candidates. This will allow, for the first time, to test gravity in the strong field regime. In this paper, we study numerically the photon orbits in the $δ= 2$ Tomimatsu-Sato space-time. The $δ= 2$ Tomimatsu-Sato space-time is a stationary, axisymmetric, and asymptotically flat exact solution of the vacuum Einstein equations. We compare the associated shadow with the one of Kerr black holes. The shape of the shadow in the $δ= 2$ Tomimatsu-Sato space-time is oblate and the difference between the two axes can be as high as 6% when viewed on the equatorial plane. We argue that future space sub-mm interferometers (e.g. VSOP-3) may distinguish the two cases, and thus are able to test the Cosmic Censorship Conjecture.

preprint2010arXiv

Thick disk accretion in Kerr space-time with arbitrary spin parameters

In this paper we extend our previous works on spherically symmetric accretion onto black holes and super-spinars to the case in which the fluid has a finite angular momentum initially. We run 2.5D and 3D general relativistic hydrodynamic simulations of the accretion of a fat disk. We study how the accretion process changes by changing the values of the parameters of our model. We show that the value of the fluid angular momentum critically determines turn-on and off the production of powerful equatorial outflows around super-spinars. For corotating disks, equatorial outflows are efficiently generated, even for relatively low spin parameters or relatively large super-spinar radii. For counterrotating disks, equatorial outflows are instead significantly suppressed, and they are possible only in limited cases. We also study accretion around a tilted disk.

preprint2009arXiv

Dark Stars: Begynnelsen

The first phase of stellar evolution in the history of the universe may be Dark Stars, powered by dark matter heating rather than by fusion. Weakly interacting massive particles, which are their own antiparticles, can annihilate and provide an important heat source for the first stars in the universe. This and the following contribution present the story of Dark Stars. In this first part, we describe the conditions under which dark stars form in the early universe: 1) high dark matter densities, 2) the annihilation products get stuck inside the star, and 3) dark matter heating wins over all other cooling or heating mechanisms.

preprint2009arXiv

Dark Stars: Död och Återuppståndelse

The first phase of stellar evolution in the history of the universe may be Dark Stars, powered by dark matter heating rather than by fusion. Weakly interacting massive particles, which are their own antiparticles, can annihilate and provide an important heat source for the first stars in the universe. This and the previous contribution present the story of Dark Stars. In this second part, we describe the structure of Dark Stars and predict that they are very massive ($\sim 800 M_\odot$), cool (6000 K), bright ($\sim 10^6 L_\odot$), long-lived ($\sim 10^6$ years), and probable precursors to (otherwise unexplained) supermassive black holes. Later, once the initial dark matter fuel runs out and fusion sets in, dark matter annihilation can predominate again if the scattering cross section is strong enough, so that a Dark Star is born again.

preprint2009arXiv

Early Black Hole Formation by Accretion of Gas and Dark Matter

We propose a model in which intermediate-mass black holes (IMBHs) with mass of ~10000 Msun are formed in early dark matter halos. We carry out detailed stellar evolution calculations for accreting primordial stars including annihilation energy of dark matter particles. We follow the stellar core evolution consistently up to gravitational collapse. We show that very massive stars, as massive as 10000 Msun, can be formed in an early dark matter halo. Such stars are extremely bright with Log L/Lsun > 8.2. They gravitationally collapse to form IMBHs. These black holes could have seeded the formation of early super-massive blackholes.

preprint2009arXiv

The formation of the first stars and galaxies

Observations made using large ground-based and space-borne telescopes have probed cosmic history all the way from the present-day to a time when the Universe was less than a tenth of its present age. Earlier on lies the remaining frontier, where the first stars, galaxies, and massive black holes formed. They fundamentally transformed the early Universe by endowing it with the first sources of light and chemical elements beyond the primordial hydrogen and helium produced in the Big Bang. The interplay of theory and upcoming observations promises to answer the key open questions in this emerging field.

preprint2009arXiv

The onset of star formation in primordial haloes

Star formation remains an unsolved problem in astrophysics. Numerical studies of large-scale structure simulations cannot resolve the whole process and their approach usually assumes that only gas denser than a typical threshold can host and form stars. We investigate the onset of cosmological star formation and compare several very-high-resolution, three-dimensional, N-body/SPH simulations that include non-equilibrium, atomic and molecular chemistry, star formation prescriptions, and feedback effects. We study how primordial star formation depends on gas density thresholds, cosmological parameters and initial set-ups. For mean-density initial conditions, we find that standard low-density star-formation threshold (0.2 h^2/cm3) models predict the onset of star formation at z~25-31, depending on the adopted cosmology. In these models stars are formed regardless of the time between the moment when the threshold is reached and the effective runaway collapse. At high redshift, this time interval represents a significant fraction of the Hubble time and thus this assumption can induce large artificial off-sets to the onset of star formation. Choosing higher density thresholds (135 h^2/cm3) allows the entire cooling process to be followed, and the onset of star formation is then estimated to be at redshift z~12-16. When isolated, rare, high-density peaks are considered, the chemical evolution is much faster and the first star formation episodes occur at z > 40, almost regardless of the choice for the density threshold. These results could have implications for the formation redshift of the first cosmological objects, as inferred from direct numerical simulations of mean-density environments, and on the studies of the reionization history of the universe.

preprint2009arXiv

Toward First-Principle Simulations of Galaxy Formation: II. Shock-Induced Starburst at a Collision Interface During the First Encounter of Interacting Galaxies

We investigated the evolution of interacting disk galaxies using high-resolution $N$-body/SPH simulations, taking into account the multiphase nature of the interstellar medium (ISM). In our high-resolution simulations, a large-scale starburst occurred naturally at the collision interface between two gas disks at the first encounter, resulting in the formation of star clusters. This is consistent with observations of interacting galaxies. The probability distribution function (PDF) of gas density showed clear change during the galaxy-galaxy encounter. The compression of gas at the collision interface between the gas disks first appears as an excess at $n_{\rm H} \sim 10{\rm cm^{-3}}$ in the PDF, and then the excess moves to higher densities ($n_{\rm H} \gtrsim 100{\rm cm^{-3}}$) in a few times $10^7$ years where starburst takes place. After the starburst, the PDF goes back to the quasi-steady state. These results give a simple picture of starburst phenomena in galaxy-galaxy encounters.

preprint2008arXiv

Modeling Nonlinear Evolution of Baryon Acoustic Oscillations: Convergence Regime of N-body Simulations and Analytic Models

We use a series of cosmological N-body simulations and various analytic models to study the evolution of the matter power spectrum in real space in a ΛCold Dark Matter universe. We compare the results of N-body simulations against three analytical model predictions; standard perturbation theory, renormalized perturbation theory, and the closure approximation. We include the effects from finite simulation box size in the comparison. We determine the values of the maximum wavenumbers, k^{lim}_{1%} and k^{lim}_{3%}, below which the analytic models and the simulation results agree to within 1 and 3 percent, respectively. We then provide a simple empirical function which describes the convergence regime determined by comparison between our simulations and the analytical models. We find that if we use the Fourier modes within the convergence regime alone, the characteristic scale of baryon acoustic oscillations can be determined within 1% accuracy from future surveys with a volume of a few h^{-3}Gpc^3 at z\sim1 or z\sim3 in the absence of any systematic distortion of the power spectrum.

preprint2008arXiv

Toward first-principle simulations of galaxy formation: I. How should we choose star formation criteria in high-resolution simulations of disk galaxies?

We performed 3-dimensional N-body/SPH simulations to study how mass resolution and other model parameters such as the star formation efficiency parameter, C* and the threshold density, nth affect structures of the galactic gaseous/stellar disk in a static galactic potential. We employ 10^6 - 10^7 particles to resolve a cold and dense (T < 100 K & n_H > 100 cm^{-3}) phase. We found that structures of the ISM and the distribution of young stars are sensitive to the assumed nth. High-nth models with nth = 100 cm^{-3} yield clumpy multi-phase features in the ISM. Young stars are distributed in a thin disk of which half-mass scale height is 10 - 30 pc. In low-nth models with nth = 0.1 cm^{-3}, the stellar disk is found to be several times thicker, and the gas disk appears smoother than the high-nth models. A high-resolution simulation with high-nth is necessary to reproduce the complex structure of the gas disk. The global properties of the model galaxies in low-nth models, such as star formation histories, are similar to those in the high-nth models when we tune the value of C* so that they reproduce the observed relation between surface gas density and surface star formation rate density. We however emphasize that high-nth models automatically reproduce the relation, regardless of the values of C*. The ISM structure, phase distribution, and distributions of young star forming region are quite similar between two runs with values of C* which differ by a factor of 15. We also found that the timescale of the flow from n_H ~1 cm^{-3} to n_H > 100 cm^{-3} is about 5 times as long as the local dynamical time and is independent of the value of C*. The use of a high-nth criterion for star formation in high-resolution simulations makes numerical models fairy insensitive to the modelling of star formation. (Abridged)

preprint2006arXiv

Formation of Primordial Stars in a LCDM Universe

We study the formation of the first generation of stars in the standard cold dark matter model, using a very high-resolution hydordynamic simulations. Our simulation achieves a dynamic range of 10^{10} in length scale. With accurate treatment of atomic and molecular physics, it allows us to study the chemo-thermal evolution of primordial gas clouds to densities up to n = 10^{16}/cc without assuming any a priori equation of state; a six orders of magnitudes improvement over previous three-dimensional calculations. All the relevant atomic and molecular cooling and heating processes, including cooling by collision-induced continuum emission, are implemented. For calculating optically thick H2 cooling at high densities, we use the Sobolev method. To examine possible gas fragmentation owing to thermal instability, we compute explicitly the growth rate of isobaric perturbations. We show that the cloud core does not fragment in either the low-density or high-density regimes. We also show that the core remains stable against gravitational deformation and fragmentation. We obtain an accurate gas mass accretion rate within a 10 Msun innermost region around the protostar. The protostar is accreting the surrounding hot gas at a rate of 0.001-0.01 Msun/yr. From these findings we conclude that primordial stars formed in early minihalos are massive. We carry out proto-stellar evolution calculations using the obtained accretion rate. The resulting mass of the first star is M_ZAMS = 60-100 Msun, with the exact mass dependent on the actual accretion rate.

preprint2003arXiv

Early Structure Formation and Reionization in a Warm Dark Matter Cosmology

We study first structure formation in Lambda-dominated universes using large cosmological N-body/SPH simulations. We consider a standard LCDM model and a LWDM model in which the mass of the dark matter particles is taken to be m_X=10 keV. The linear power spectrum for the LWDM model has a characteristic cut-off at a wavenumber k=200 /Mpc, suppressing the formation of low mass (< 10^6 Msun) nonlinear objects early on. The absence of low mass halos in the WDM model makes the formation of primordial gas clouds with molecular hydrogen very inefficient at high redshifts. The first star-forming gas clouds form at z~21 in the WDM model, considerably later than in the CDM counterpart, and the abundance of these gas clouds differs by an order of magnitude between the two models. We carry out radiative transfer calculations by embedding massive Population III stars in the gas clouds. We show that the volume fraction of ionized gas rises up close to 100% by z=18 in the CDM case, whereas that of the WDM model remains extremely small at a level of a few percent. Thus the WDM model with m_X=10 keV is strongly inconsistent with the observed high optical depth by the WMAP satellite.

preprint2002arXiv

Improving the Accuracy of Cosmic Magnification Statistics

The systematic magnification of background sources by the weak gravitational-lensing effects of foreground matter, also called cosmic magnification, is becoming an efficient tool both for measuring cosmological parameters and for exploring the distribution of galaxies relative to the dark matter. We extend here the formalism of magnification statistics by estimating the contribution of second-order terms in the Taylor expansion of the magnification and show that the effect of these terms was previously underestimated. We test our analytical predictions against numerical simulations and demonstrate that including second-order terms allows the accuracy of magnification-related statistics to be substantially improved. We also show, however, that both numerical and analytical estimates can provide only lower bounds to real correlation functions, even in the weak lensing regime. We propose to use count-in-cells estimators rather than correlation functions for measuring cosmic magnification since they can more easily be related to correlations measured in numerical simulations.

preprint2001arXiv

Clustering of dark matter halos on the light-cone: scale-, time- and mass-dependence of the halo biasing in the Hubble volume simulations

We develop a phenomenological model to predict the clustering of dark matter halos on the light-cone by combining several existing theoretical models. Assuming that the velocity field of halos on large scales is approximated by linear theory, we propose an empirical prescription of a scale-, mass-, and time-dependence of halo biasing. We test our model against the Hubble Volume $N$-body simulation and examine its validity and limitations. We find a good agreement in two-point correlation functions of dark matter halos between the phenomenological model predictions and measurements from the simulation for $R>5h^{-1}$Mpc both in the real and redshift spaces. Although calibrated on the mass scale of groups and clusters and for redshifts up to $z\sim2$, the model is quite general and can be applied to a wider range of astrophysical objects, such as galaxies and quasars, if the relation between dark halos and visible objects is specified.

preprint2000arXiv

Weakly Self-Interacting Dark Matter and the Structure of Dark Halos

We study the formation of dark halos in a $Λ$CDM universe under the assumption that Cold Dark Matter particles have a finite cross-section for elastic collisions. We compare evolution when CDM mean free paths are comparable to halo sizes with the collisionless and fluid limits. We show that a few collisions per particle per Hubble time at halo centre can substantially affect the central density profile. Cross-sections an order of magnitude larger produce sufficient relaxation for rich clusters to develop core radii in the range 100-200 $h^{-1}$kpc. The structural evolution of halos is a competition between collisional relaxation caused by individual particle interactions and violent relaxation resulting from the infall and merging processes by which clusters grow. Although our simulations concentrate on systems of cluster size, we can scale our results to address the halo structure expected for dwarf galaxies. We find that collision cross-sections sufficiently large to significantly modify the cores of such galaxies produce cluster cores which are too large and/or too round to be consistent with observation. Thus the simplest model for self-interacting dark matter is unable to improve fits to published dwarf galaxy rotation curves without violating other observational constraints.

Naoki Yoshida

What is connected

Connect this record

See the researcher in context

Building this map preview

109 published item(s)

A merger-driven scenario for clumpy galaxy formation in the epoch of reionization: Physical properties of clumps in the FirstLight simulation

$\mathrm{H_2}$ cooling and gravitational collapse of supersonically induced gas objects

Deep-learning Real/Bogus classification for the Tomo-e Gozen transient survey

Detecting preheating in proto-clusters with Lyman-$α$ Forest Tomography

Possible Systematic Rotation in the Mature Stellar Population of a $z=9.1$ Galaxy

Public Release of A-SLOTH: Ancient Stars and Local Observables by Tracing Halos

Stellar Wind Effect on the Atmospheric Escape of Hot Jupiters and their Ly-$α$ and H$α$ transits

The three-year shear catalog of the Subaru Hyper Suprime-Cam SSP Survey

Three-Dimensional Reconstruction of Weak Lensing Mass Maps with a Sparsity Prior. I. Cluster Detection

Atmospheric Escape of Close-in Giants around Hot Stars: Far-Ultraviolet Radiation and Photoelectric Heating Effect

Capturing the inside-out quenching by black holes with far-infrared atomic line ratios

Large-scale variation in reionization history caused by Baryon-dark matter streaming velocity

ALMA uncovers the [CII] emission and warm dust continuum in a z = 8.31 Lyman break galaxy

Cross Correlation of the Extragalactic Gamma-ray Background with Thermal Sunyaev-Zel'dovich Effect in the Cosmic Microwave Background

Cross-correlation of the thermal Sunyaev-Zel'dovich effect and weak gravitational lensing: Planck and Subaru Hyper Suprime-Cam first-year data

Deep learning for intensity mapping observations: Component extraction

Deep Learning for Line Intensity Mapping Observations: Information Extraction from Noisy Maps

First Structure Formation under the Influence of Gas-Dark Matter Streaming Velocity and Density: Impact of the Baryons-trace-dark matter Approximation

Photoevaporation of Minihalos during Cosmic Reionization: Primordial and Metal-Enriched Halos

Photometric classification of HSC transients using machine learning

R-process enrichment in ultrafaint dwarf galaxies

Rapid Transients Originating from Thermonuclear Explosions in Helium White Dwarf Tidal Disruption Events

The onset of star formation 250 million years after the Big Bang

Cosmological Constraint on the Light Gravitino Mass from CMB Lensing and Cosmic Shear

Cosmological constraints on dark matter annihilation and decay: Cross-correlation analysis of the extragalactic $γ$-ray background and cosmic shear

Cosmological Simulations of Early Blackhole Formation: Halo Mergers, Tidal Disruption, and the Conditions for Direct Collapse

Detection of an oxygen emission line from a high redshift galaxy in the reionization epoch

Eigenanalysis of morphological diversity in silicon random nanostructures formed via resist collapse

Formation of Massive Primordial Stars: Intermittent UV Feedback with Episodic Mass Accretion

Machine-learning Selection of Optical Transients in Subaru/Hyper Suprime-Cam Survey

Matter Power Spectrum in Hidden Neutrino Interacting Dark Matter Models: A Closer Look at the Collision Term

Nebular line emission from z > 7 galaxies in cosmological simulations: rest-frame UV to Optical lines

Optical nano artifact metrics using silicon random nanostructures

The Final Fates of Accreting Supermassive Stars

The Hydrodynamic Feedback of Cosmic Reionization on Small-Scale Structures and Its Impact on Photon Consumption during the Epoch of Reionization

The Subaru FMOS galaxy redshift survey (FastSound). IV. New constraint on gravity theory from redshift space distortions at $z\sim 1.4$

Cross-Correlation of the Extragalactic Gamma-ray Background with Luminous Red Galaxies

Early structure formation from primordial density fluctuations with a blue-tilted power spectrum

Impact of Baryonic Processes on Weak Lensing Cosmology: Power Spectrum, Non-Local Statistics, and Parameter Bias

Particle splitting in smoothed particle hydrodynamics based on Voronoi diagram

Primordial Star Formation under the Influence of Far Ultraviolet Radiation: 1540 Cosmological Halos and the Stellar Mass Distribution

Probing cosmology with weak lensing selected clusters I: Halo approach and all-sky simulations

Probing cosmology with weak lensing selected clusters II: Dark energy and f(R) gravity models

Probing the Small Scale Matter Power Spectrum through Dark Matter Annihilation in the Early Universe

Supermassive star formation via episodic accretion: protostellar disc instability and radiative feedback efficiency

Cross-Correlation of Cosmic Shear and Extragalactic Gamma-ray Background: Constraints on the Dark Matter Annihilation Cross-Section

Dust grain growth and the formation of the extremely primitive star SDSS J102915+172927

On the Systematic Errors of Cosmological-Scale Gravity Tests using Redshift Space Distortion: Non-linear Effects and the Halo Bias

One Hundred First Stars : Protostellar Evolution and the Final Masses

Physical Properties of UDF12 Galaxies in Cosmological simulations

Reproducing cosmic evolution of galaxy population from $z = 4$ to $0$

Statistical and Systematic Errors in Measurement of Weak-Lensing Minkowski Functionals: Application to Canada-France-Hawaii Lensing Survey

Supernova dust formation and the grain growth in the early universe: The critical metallicity for low-mass star formation

The Dark Ages of the Universe and Hydrogen Reionization

The origin of the most iron-poor star

ALMA will determine the spectroscopic redshift z>8 with FIR [O III] emission lines

Constraints on Warm Dark Matter models from high-redshift long gamma-ray bursts

Detectability of High-Redshift Superluminous Supernovae with Upcoming Optical and Near-Infrared Surveys - II. Beyond z=6

Effect of Mask Regions on Weak Lensing Statistics

Formation of Primordial Supermassive Stars by Rapid Mass Accretion

Growth of dust grains in a low-metallicity gas and its effect on the cloud fragmentation

Modeling color-dependent galaxy clustering in cosmological simulations

On de-Sitter Geometry in Cosmic Void Statistics

Structure of Dark Matter Halos in Warm Dark Matter models and in models with Long-Lived Charged Massive Particles

Weighing the Light Gravitino Mass with Weak Lensing Surveys

Detectability of High-Redshift Superluminous Supernovae with Upcoming Optical and Near-Infrared Surveys

Direct Integration of the Collisionless Boltzmann Equation in Six-dimensional Phase Space: Self-gravitating Systems

Distribution of Dust around Galaxies: An Analytic Model

Formation and evolution of primordial protostellar systems

Light Curve Modeling of Superluminous Supernova 2006gy: Collision between Supernova Ejecta and Dense Circumstellar Medium

Low-mass star formation triggered by early supernova explosions

Probing Primordial Non-Gaussianity with Weak Lensing Minkowski Functionals

Protostellar Feedback and Final Mass of the Second-Generation Primordial Stars

Radiative cooling implementations in simulations of primordial star formation