Researcher profile

Naoki Yoshida

Naoki Yoshida contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
37works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

37 published item(s)

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&#39; 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&#39;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&#39;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&#39;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&#39;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&#39;dovich effect and weak gravitational lensing: Planck and Subaru Hyper Suprime-Cam first-year data

Cross-correlation analysis of the thermal Sunyaev-Zel&#39;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&#39;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&#34; 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 &#34;minihalos&#34; 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.

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.

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

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.

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

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 &#34;shadow&#34; 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.

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.

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.