Source author record

Hidenobu Yajima

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

29works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

29 published item(s)

preprint2022arXiv

Far and extreme UV radiation feedback in molecular clouds and its influence on the mass and size of star clusters

We study the formation of star clusters in molecular clouds by performing three-dimensional radiation hydrodynamics simulations with far ultraviolet (FUV; $6 ~{\rm eV} \leqq h ν\leqq 13.6 ~{\rm eV}$) and extreme ultraviolet (EUV; $hν\geqq 13.6~{\rm eV}$) radiative feedback. We find that the FUV feedback significantly suppresses the star formation in diffuse clouds with the initial surface densities of $Σ_{\rm cl} \lesssim \rm 50~M_{\odot} \; pc^{-2}$. In the cases of clouds with $Σ_{\rm cl} \sim \rm 100-200~M_{\odot} \; pc^{-2}$, the EUV feedback plays a main role and decrease the star formation efficiencies less than $0.3$. We show that thermal pressure from PDRs or H{\sc ii} regions disrupts the clouds and makes the size of the star clusters larger. Consequently, the clouds with the mass $M_{\rm cl} \lesssim 10^{5}~\rm M_{\odot}$ and the surface density $Σ_{\rm cl} \lesssim 200~\rm M_{\odot}\; pc^{-2}$ remain the star clusters with the stellar densities of $\sim 100~\rm M_{\odot}\; pc^{-3}$ that nicely match the observed open clusters in the Milky Way. If the molecular clouds are massive ($M_{\rm cl} \gtrsim 10^{5}~\rm M_{\odot}$) and compact ($Σ\gtrsim 400~\rm M_{\odot}\; pc^{-2}$), the radiative feedback is not effective and they form massive dense cluster with the stellar densities of $\sim 10^{4}~\rm M_{\odot}\; pc^{-3}$ like observed globular clusters or young massive star clusters. Thus, we suggest that the radiative feedback and the initial conditions of molecular clouds are key factors inducing the variety of the observed star clusters.

preprint2022arXiv

FOREVER22: galaxy formation in protocluster regions

We present results from a new cosmological hydrodynamics simulation campaign of protocluster (PC) regions, FOREVER22: FORmation and EVolution of galaxies in Extremely-overdense Regions motivated by SSA22. The simulations cover a wide range of cosmological scales using three different zoom set-ups in a parent volume of $(714.2~\rm cMpc)^{3}$: PCR (Proto-Cluster Region; $V= (28.6~{\rm cMpc})^{3} $, SPH particle mass, $m_{\rm{SPH}} = 4.1 \times 10^{6}~\rm M_{\odot}$ and final redshift, $z_{\rm end}=2.0$), BCG (Brightest proto-Cluster Galaxy; $V \sim (10~{\rm cMpc})^{3} $, $m_{\rm SPH} = 5.0\times10^{5}~\rm M_{\odot}$ and $z_{\rm end}=4.0$ ), and First ( $V \sim (3~{\rm cMpc})^{3} $, $m_{\rm SPH} = 7.9 \times 10^{3}~\rm M_{\odot}$ and $z_{\rm end}=9.5$) runs, that allow to focus on different aspects of galaxy formation. In the PCR runs, we follow 10 PCs, each harbouring 1 - 4 SMBHs with $M_{\rm BH} \ge 10^{9}~\rm M_{\odot}$. One of the PC cores shows a spatially close arrangement of seven starburst galaxies with ${\rm SFR} \gtrsim 100~\rm M_{\odot}~yr^{-1}$ each, that are dust-obscured and would appear as submillimeter galaxies with flux $\gtrsim 1~$ mJy at $1.1~ \rm mm$ in observations. The BCG runs show that the total SFRs of haloes hosting BCGs are affected by AGN feedback, but exceed $1000~\rm M_{\odot}~yr^{-1}$ at $z \lesssim 6$. The First runs resolve mini-haloes hosting population (Pop) III stars and we show that, in PC regions, the dominant stellar population changes from Pop III to Pop II at $z \gtrsim 20$, and the first galaxies with ${\rm SFR} \gtrsim 18~\rm M_{\odot}~yr^{-1}$ form at $z \sim 10$. These can be prime targets for future observations with the James Webb Space Telescope. Our simulations successfully reproduce the global star formation activities in observed PCs and suggest that PCs can kickstart cosmic reionization.

preprint2022arXiv

TRINITY: a three-dimensional time-dependent radiative transfer code for in-vivo near-infrared imaging

We develop a new three-dimensional time-dependent radiative transfer code, TRINITY (Time-dependent Radiative transfer In Near-Infrared TomographY), for in-vivo diffuse optical tomography (DOT). The simulation code is based on the design of long radiation rays connecting boundaries of a computational domain, which allows us to calculate light propagation with little numerical diffusion. We parallelize the code with Message Passing Interface (MPI) using the domain decomposition technique and confirm the high parallelization efficiency, so that simulations with a spatial resolution of $\sim 1$ millimeter can be performed in practical time. As a first application, we study the light propagation for a pulse collimated within $θ\sim 15^\circ$ in a phantom, which is a uniform medium made of polyurethane mimicking biological tissue. We show that the pulse spreads in all forward directions over $\sim$ a few millimeters due to the multiple scattering process of photons. Our simulations successfully reproduce the time-resolved signals measured with eight detectors for the phantom. We also introduce the effects of reflection and refraction at the boundary of medium with a different refractive index and demonstrate the faster propagation of photons in an air hole that is an analogue for the respiratory tract.

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.

preprint2020arXiv

ART$^2$: A 3D Parallel Multi-wavelength Radiative Transfer Code for Continuum and Atomic and Molecular Lines

ART$^2$ is a 3D multi-wavelength Monte Carlo radiative transfer (RT) code that couples continuum and emission lines to track the propagation of photons and their interactions with the interstellar medium (ISM). The previous version of ART$^2$, which included continuum and Ly$α$ line, has been extensively applied to hydrodynamics simulations to study multi-band properties of galaxies and ISM. Here, we describe new implementations of non-local thermodynamic equilibrium RT of molecular and atomic fine structure emission lines, and the parallelization of the code using a number of novel methods. The new ART$^2$ can efficiently and self-consistently produce a full spectrum that includes both continuum and lines such as [CII], [NII], [OIII], Ly$α$, and CO. These essential features, together with the multi-phase ISM model and the adaptive grid, make ART$^2$ a multi-purpose code to study multi-wavelength properties of a wide range of astrophysical systems from planetary disks to large-scale structures. To demonstrate the capability of the new ART$^2$, we applied it to two hydrodynamics simulations: the zoom-in Milky Way Simulation to obtain panchromatic properties of individual galaxies, and the large-scale IllustrisTNG100 Simulation to obtain global properties such as the line intensity mappings. These products are essential for a broad array of studies such as the correlations between physical and panchromatic properties and their evolution. By enabling direct comparison between numerical simulations and multi-band observations, ART$^2$ provides a crucial theoretical framework for the interpretations of existing observations, the plan for future surveys, and the synergy between multi-band galaxy surveys and line intensity mappings. Therefore, ART$^2$ is a powerful and versatile tool to bridge the gap between theories and observations of cosmic structures.

preprint2020arXiv

FIR-luminous [CII] emitters in the ALMA-SCUBA-2 COSMOS survey (AS2COSMOS): The nature of submillimeter galaxies in a 10 comoving Mpc-scale structure at z~4.6

We report the discovery of a 10 comoving Mpc-scale structure traced by massive submillimeter galaxies (SMGs) at z~4.6. These galaxies are selected from an emission line search of ALMA Band 7 observations targeting 184 luminous submillimeter sources ($S_{850μ{\rm m}}\geq$ 6.2 mJy) across 1.6 degrees$^2$ in the COSMOS field. We identify four [CII] emitting SMGs and two probable [CII] emitting SMG candidates at z=4.60-4.64 with velocity-integrated signal-to-noise ratio of SNR>8. Four of the six emitters are near-infrared blank SMGs. After excluding one SMG whose emission line is falling at the edge of the spectral window, all galaxies show clear velocity gradients along the major axes that are consistent with rotating gas disks. The estimated rotation velocities of the disks are 330-550 km s$^{-1}$ and the inferred host dark-matter halo masses are ~2-8 $\times$ 10$^{12}$M$_{\odot}$. From their estimated halo masses and [CII] luminosity function, we suggest that these galaxies have a high (50-100%) duty cycle and high (~0.1) baryon conversion efficiency (SFR relative to baryon accretion rate), and that they contribute $\simeq$2% to the total star-formation rate density at z=4.6. These SMGs are concentrated within just 0.3% of the full survey volume, suggesting they are strongly clustered. The extent of this structure and the individual halo masses suggest that these SMGs will likely evolve into members of a ~10$^{15}$M$_{\odot}$ cluster at z=0. This survey reveals synchronized dusty starburst in massive halos at z>4, which could be driven by mergers or fed by smooth gas accretion.

preprint2020arXiv

High circular polarization of near infrared light induced by micron-size dust grains

We explore the induction of circular polarization (CP) of near-infrared light in star-forming regions using three-dimensional radiative transfer calculations. The simulations trace the change of Stokes parameters at each scattering/absorption process in a dusty gas slab composed of aligned grains. We find that the CP degree enlarges significantly according as the size of dust grains increases and exceeds $\sim 20$ percent for micron-size grains. Therefore, if micron-size grains are dominant in a dusty gas slab, the high CP observed around luminous young stellar objects can be accounted for. The distributions of CP show the asymmetric quadrupole patters regardless of the grain sizes. Also, we find that the CP degree depends on the relative position of a dusty gas slab. If a dusty gas slab is located behind a star-forming region, the CP reaches $\sim 60$ percent in the case of $1.0~{\rm μm}$ dust grains. Hence, we suggest that the observed variety of CP maps can be explained by different size distributions of dust grains and the configuration of aligned grains.

preprint2020arXiv

Star cluster formation and cloud dispersal by radiative feedback: dependence on metallicity and compactness

We study star cluster formation in various environments with different metallicities and column densities by performing a suite of three-dimensional radiation hydrodynamics simulations. We find that the photoionization feedback from massive stars controls the star formation efficiency (SFE) in a star-forming cloud, and its impact sensitively depends on the gas metallicity $Z$ and initial cloud surface density $Σ$. At $Z=1~Z_{\odot}$, SFE increases as a power law from 0.03 at $Σ= 10~M_{\odot}{\rm pc^{-2}}$ to 0.3 at $Σ= 300~M_{\odot}{\rm pc^{-2}}$. In low-metallicity cases $10^{-2}- 10^{-1} Z_{\odot}$, star clusters form from atomic warm gases because the molecule formation time is not short enough with respect to the cooling or dynamical time. In addition, the whole cloud is disrupted more easily by expanding H{\sc ii} bubbles which have higher temperature owing to less efficient cooling. With smaller dust attenuation, the ionizing radiation feedback from nearby massive stars is stronger and terminate star formation in dense clumps. These effects result in inefficient star formation in low-metallicity environments: the SFE drops by a factor of $\sim 3$ at $Z=10^{-2}~Z_{\odot}$ compared to the results for $Z=1~Z_{\odot}$, regardless of $Σ$. Newborn star clusters are also gravitationally less bound. We further develop a new semi-analytical model that can reproduce the simulation results well, particularly the observed dependencies of the SFEs on the cloud surface densities and metallicities.

preprint2019arXiv

CHORUS. III. Photometric and Spectroscopic Properties of Ly$α$ Blobs at $z=4.9-7.0$

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

preprint2019arXiv

First Identification of 10-kpc Scale [CII] 158um Halos around Star-Forming Galaxies at z=5-7

We report the discovery of 10-kpc scale [CII] 158um halos surrounding star-forming galaxies in the early Universe. We choose deep ALMA data of 18 galaxies each with a star-formation rate of ~ 10-70 Msun with no signature of AGN whose [CII] lines are individually detected at z=5.153-7.142, and conduct stacking of the [CII] lines and dust-continuum in the uv-visibility plane. The radial profiles of the surface brightnesses show a 10-kpc scale [CII] halo at the 9.2sigma level significantly extended more than the HST stellar continuum data by a factor of ~5 on the exponential-profile basis, as well as the dust continuum. We also compare the radial profiles of [CII] and Lya halos universally found in star-forming galaxies at this epoch, and find that the scale lengths agree within the 1sigma level. While two independent hydrodynamical zoom-in simulations match the dust and stellar continuum properties, the simulations cannot reproduce the extended [CII] line emission. The existence of the extended [CII] halo is the evidence of outflow remnants in the early galaxies and suggest that the outflows may be dominated by cold-mode outflows expelling the neutral gas.

preprint2019arXiv

Galaxy evolution and radiative properties in the early Universe: multi-wavelength analysis in cosmological simulations

Recent observations have successfully detected UV or infrared flux from galaxies at the epoch of reionization. However, the origin of their radiative properties has not been fully understood yet. Combining cosmological hydrodynamic simulations and radiative transfer calculations, we present theoretical predictions of multi-wavelength radiative properties of the first galaxies at z=6-15. We find that most of the gas and dust are ejected from star-forming regions due to supernova (SN) feedback, which allows UV photons to escape. We show that the peak of SED rapidly shifts between UV and infrared wavelengths on a timescale of 100 Myr due to intermittent star formation and feedback. When dusty gas covers the star-forming regions, the galaxies become bright in the observed-frame sub-millimeter wavelengths. In addition, we find that the escape fraction of ionizing photons also changes between 1-40% at z>10. The mass fraction of HII region changes with the star formation history, resulting in the fluctuations of metal lines and Lyman-alpha line luminosities. In the starbursting phase of galaxies with the halo mass $\sim 10^{11}\,{\rm M_{\odot}}$ ($10^{12}\,{\rm M_{\odot}}$), the simulated galaxy has $L_{\rm [O\,III]} \sim 10^{42}\, (10^{43})\, {\rm erg\,s^{-1}}$, which is consistent with the observed star-forming galaxies at z>7. Our simulations suggest that deep [C II] observation with ALMA can trace the distribution of neutral gas extending over $\sim 20$ physical kpc. We also find that the luminosity ratio $L_{\rm [O\,III]}/L_{\rm [C\,II]}$ decreases with bolometric luminosity due to metal enrichment. Our simulations show that the combination of multi-wavelength observations by ALMA and JWST will be able to reveal the multi-phase ISM structure and the transition from starbursting to outflowing phases of high-z galaxies.

preprint2016arXiv

The role of stellar relaxation in the formation and evolution of the first massive black holes

We present calculations on the formation of massive black holes with 10^5 Msun at z > 6 that can be the seeds of supermassive black holes at z > 6. Under the assumption of compact star cluster formation in merging galaxies, star clusters in haloes of 10^8 ~ 10^9 Msun can undergo rapid core-collapse leading to the formation of very massive stars (VMSs) with ~1000 Msun which directly collapse into black holes with similar masses. Star clusters in halos of > 10^9 Msun experience type-II supernovae before the formation of VMSs due to long core-collapse time scales. We also model the subsequent growth of black holes via accretion of residual stars in clusters. 2-body relaxation efficiently re-fills the loss cones of stellar orbits at larger radii and resonant relaxation at small radii is the main driver for accretion of stars onto black holes. As a result, more than ninety percent of stars in the initial cluster are swallowed by the central black holes before z=6. Using dark matter merger trees we derive black hole mass functions at z=6-20. The mass function ranges from 10^3 to 10^5 Msun at z <~ 15. Major merging of galaxies of >~ 4*10^8 Msun at z ~ 20 successfully leads to the formation of >~ 10^5 Msun BHs by z >~ 10 which can be the potential seeds of supermassive black holes seen today.

preprint2015arXiv

Can the 21 cm signal probe Population III and II star formation?

Using varying models for the star formation rate (SFR) of Population (Pop) III and II stars at z>6 we derive the expected redshift history of the global 21 cm signal from the inter-galactic medium (IGM). To recover the observed Thomson scattering optical depth of the cosmic microwave background (CMB) requires SFRs at the level of ~ 10^-3 Msun/yr/Mpc^3 at z ~ 15 from Pop III stars, or ~ 10^-1Msun/yr/Mpc^3 at z ~ 7 from Pop II stars. In the case the SFR is dominated by Pop III stars, the IGM quickly heats above the CMB at z > 12 due to heating from supernovae. In addition, Lyman-alpha photons from haloes hosting Pop III stars couple the spin temperature to that of the gas, resulting in a deep absorption signal. If the SFR is dominated by Pop II stars, the IGM slowly heats and exceeds the CMB temperature at z ~ 10. However, the larger and varying fraction of Pop III stars is able to break this degeneracy. We find that the impact of the initial mass function (IMF) of Pop III stars on the 21 cm signal results in an earlier change to a positive signal if the IMF slope is ~ -1.2. Measuring the 21 cm signal at z > 10 with next generation radio telescopes such as the Square Kilometre Array will be able to investigate the contribution from Pop III and Pop II stars to the global star formation rate.

preprint2015arXiv

Cold accretion in early galaxy formation and its Lyman-alpha signatures

The Lyman-alpha (Lya) emission has played an important role in detecting high-redshift galaxies, including recently distant ones at redshift z > 7. It may also contain important information on the origin of these galaxies. Here, we investigate the formation of a typical L* galaxy and its observational signatures at the earliest stage, by combining a cosmological hydrodynamic simulation with three-dimensional radiative transfer calculations using the newly improved ART^2 code. Our cosmological simulation uses the Aquila initial condition which zooms in onto a Milky Way-like halo with high resolutions, and our radiative transfer couples multi-wavelength continuum, Lya line, and ionization of hydrogen. We find that the modeled galaxy starts to form at redshift z ~ 24 through efficient accretion of cold gas, which produces a strong Lya line with a luminosity of L(Lya) ~ 10^42 erg/s as early as z ~ 14. The Lya emission appears to trace the cold, dense gas. The lines exhibit asymmetric, single-peak profiles, and are shifted to the blue wing, a characteristic feature of gas inflow. Moreover, the contribution to the total Lya luminosity by excitation cooling increases with redshift, and it becomes dominant at z >~ 6. We predict that L* galaxies such as the modeled one may be detected at z <~ 8 by JWST and ALMA with a reasonable integration time. Beyond redshift 12, however, only Lya line may be observable by spectroscopic surveys. Our results suggest that Lya line is one of the most powerful tools to detect the first generation of galaxies, and to decipher their formation mechanism.

preprint2015arXiv

Observational Properties of Simulated Galaxies in Overdense and Average Regions at High Redshifts z= 6-12

We use high-resolution zoom-in cosmological simulations of galaxies of Romano-Diaz et al., post-processing them with a panchromatic three-dimensional radiation transfer code to obtain the galaxy UV luminosity function (LF) at z ~ 6-12. The galaxies are followed in a rare, heavily overdense region within a ~ 5-sigma density peak, which can host high-z quasars, and in an average density region, down to the stellar mass of M_star ~ 4* 10^7 Msun. We find that the overdense regions evolve at a substantially accelerated pace --- the most massive galaxy has grown to M_star ~ 8.4*10^10 Msun by z = 6.3, contains dust of M_dust~ 4.1*10^8 Msun, and is associated with a very high star formation rate, SFR ~ 745 Msun/yr.The attained SFR-M_star correlation results in the specific SFR slowly increasing with M_star. Most of the UV radiation in massive galaxies is absorbed by the dust, its escape fraction f_esc is low, increasing slowly with time. Galaxies in the average region have less dust, and agree with the observed UV LF. The LF of the overdense region is substantially higher, and contains much brighter galaxies. The massive galaxies are bright in the infrared (IR) due to the dust thermal emission, with L_IR~ 3.7*10^12 Lsun at z = 6.3, while L_IR < 10^11 Lsun for the low-mass galaxies. Therefore, ALMA can probe massive galaxies in the overdense region up to z ~ 10 with a reasonable integration time. The UV spectral properties of disky galaxies depend significantly upon the viewing angle.The stellar and dust masses of the most massive galaxy in the overdense region are comparable to those of the sub-millimetre galaxy (SMG) found by Riechers et al. at z = 6.3, while the modelled SFR and the sub-millimetre flux fall slightly below the observed one. Statistical significance of these similarities and differences will only become clear with the upcoming ALMA observations.

preprint2014arXiv

Angular momentum loss of primordial gas in Lyman-alpha radiation field

We present results on the radiation drag exerted by an isotropic and homogenous background of Lya photons on neutral gas clouds orbiting within HII regions around Population III stars of different masses. The Doppler shift causes a frequency difference between photons moving in the direction of the cloud and opposite to it resulting in a net momentum loss of the cloud in the direction of motion. We find that half of the angular momentum of gas with v_theta <~ 20 km/s near (r <~ 3 kpc) a Population III star of 120 Msun at z=20 is lost within ~ 10^6 yr. The radiation drag is a strong function of cloud velocity that peaks at v ~ 20 km/s reflecting the frequency dependence of the photon cross section. Clouds moving with velocities larger than ~ 100 km/s loose their angular momentum on time scales of ~ 10^8 yr. At lower redshifts radiation drag becomes inefficient as the Lya photon density in HII regions decreases by a factor (1+z)^3 and angular momentum is lost on time scales > 10^8 yr even for low velocity clouds. Our results suggest that a sweet spot exists for the loss of angular momentum by radiation drag for gas clouds at z > 10 and with v ~ 20 km/s. Comparison to dynamical friction forces acting on typical gas clouds suggest that radiation drag is the dominant effect impacting the orbit. We propose that this effect can suppress the formation of extended gas discs in the first galaxies and help gas accretion near galactic centres and central black holes.

preprint2014arXiv

Distinctive 21 cm structures of the first stars, galaxies, and quasars

Observations of the redshifted 21 cm line with upcoming radio telescopes promise to transform our understanding of the cosmic reionization. To unravel the underlying physical process, we investigate the 21 cm structures of three different ionizing sources, Population (Pop) III stars, the first galaxies, and the first quasars, by using radiative transfer simulations that include both ionization of neutral hydrogen and resonant scattering of Lya photons. We find that Pop III stars and quasars produce a smooth transition from an ionized and hot state to a neutral and cold one, owing to their hard spectral energy distribution with abundant ionizing photons, in contrast to the sharp transition in galaxies. Furthermore, Lya scattering plays a dominant role in producing the 21 cm signal as it determines the relation between hydrogen spin temperature and gas kinetic temperature. This effect, also called Wouthuysen-Field coupling, depends strongly on the ionizing source. It is the strongest around galaxies, where the spin temperature is highly coupled to that of the gas, resulting in extended absorption troughs in the 21 cm brightness temperature. On the other hand, in the case of Pop III stars, the 21 cm signal shows both emission and absorption regions around a small HII bubble. For quasars, a large emission region in the 21 cm signal is produced, and the absorption region decreases as the size of the HII bubble becomes large due to the limited traveling time of photons. We predict that future surveys from large radio arrays such as MWA, LOFAR and SKA may be able to detect the 21 cm signals of primordial galaxies and quasars, but not likely Pop III stars due to its small angular diameter.

preprint2014arXiv

Dust properties of Lyman break galaxies in cosmological simulations

Recent observations have indicated the existence of dust in high-redshift galaxies, however, the dust properties in them are still unknown. Here we present theoretical constraints on dust properties in Lyman break galaxies (LBGs) at z=3 by post-processing a cosmological smoothed particle hydrodynamics simulation with radiative transfer calculations. We calculate the dust extinction in 2800 dark matter halos using the metallicity information of individual gas particles in our simulation. We use only bright galaxies with rest-frame UV magnitude M_1700 < -20 mag, and study the dust size, dust-to-metal mass ratio, and dust composition. From the comparison of calculated color excess between B and V-band (i.e., E(B-V)) and the observations, we constrain the typical dust size, and show that the best-fitting dust grain size is ~ 0.05 micron, which is consistent with the results of theoretical dust models for Type-II supernova. Our simulation with the dust extinction effect can naturally reproduce the observed rest-frame UV luminosity function of LBGs at z=3 without assuming an ad hoc constant extinction value. In addition, in order to reproduce the observed mean E(B-V), we find that the dust-to-metal mass ratio needs to be similar to that of the local galaxies, and that the graphite dust is dominant or at least occupy half of dust mass.

preprint2014arXiv

Escape of Lyman-alpha and continuum photons from star-forming galaxies

A large number of high-redshift galaxies have been discovered via their narrow-band Lya line or broad-band continuum colors in recent years. The nature of the escaping process of photons from these early galaxies is crucial to understanding galaxy evolution and the cosmic reionization. Here, we investigate the escape of Lya, non-ionizing UV-continuum (l = 1300 - 1600 angstrom in rest frame), and ionizing photons (l < 912 angstrom) from galaxies by combining a cosmological hydrodynamic simulation with three-dimensional multi-wavelength radiative transfer calculations. The galaxies are simulated in a box of 5^3 h^-3 Mpc^3 with high resolutions using the Aquila initial condition which reproduces a Milky Way-like galaxy at redshift z=0. We find that the escape fraction (fesc) of these different photons shows a complex dependence on redshift and galaxy properties: fesc(Lya) and fesc(UV) appear to evolve with redshift, and they show similar, weak correlations with galaxy properties such as mass, star formation, metallicity, and dust content, while fesc(Ion) remains roughly constant at ~ 0.2 from z ~ 0 - 10, and it does not show clear dependence on galaxy properties. fesc(Lya) correlates more strongly with fesc(UV) than with fesc(Ion). In addition, we find a relation between the emergent Lya luminosity and the ionizing photon emissivity of Lyman Alpha Emitters (LAEs). By combining this relation with the observed luminosity functions of LAEs at different redshift, we estimate the contribution from LAEs to the reionization of intergalactic medium (IGM). Our result suggests that ionizing photons from LAEs alone are not sufficient to ionize IGM at z > 6, but they can maintain the ionization of IGM at z ~ 0 - 5.

preprint2013arXiv

Effects of photon trapping on the Lyman-alpha properties of star-forming galaxies

Recent observations show that a large number of Lyman-alpha emitters (LAEs) at high redshift z >= 3 have unusually high Lya equivalent widths (EW > 400 angstrom). However, the origin of these high EWs is an open question. Here, we investigate the impacts of photon trapping on the Lya EW and other properties by tracking the Lya radiative transfer in spherical galactic clouds. We find that the delayed escape of the Lya photons can change the Lya properties significantly. During the transition phase from optically thick to optically thin where the Lya photons can escape simultaneously, the EW can be boosted to ~ 1000 angstrom, the Lya luminosity can be increased by a factor of a few, and the line profile can be significantly broadened. The boost factor appears to depend on the galaxy properties such as mass and star formation rate and timescale, therefore future investigation combing 3D Lya RT calculations with cosmological simulations of galaxy formation and evolution is needed to fully understand the Lya properties of early star-forming galaxies.

preprint2013arXiv

Extended Lyman-alpha emission from interacting galaxies at high redshifts

Recent observations have discovered a population of extended Lya sources, dubbed Lyman-alpha blobs (LABs), at high redshift z ~ 3 - 6.6. These LABs typically have a luminosity of L ~ 10^42-10^44 erg/s, and a size of tens of kiloparsecs, with some giant ones reaching up to D ~ 100 kpc. However, the origin of these LABs is not well understood. In this paper, we investigate a merger model for the formation of LABs by studying Lya emission from interacting galaxies at high redshifts by means of a combination of hydrodynamics simulations with three dimensional radiative transfer calculations. Our galaxy simulations focus on a set of binary major mergers of galaxies with a mass range of 3-7 *10^12 Msun in the redshift range of z ~ 3 -7, and we use the newly improved ART^2 code to perform the radiative transfer calculations which couple multi-wavelength continuum, ionization of hydrogen, and Lya line emission. We find that intense star formation and enhanced cooling induced by gravitational interaction produce strong Lya emission from these merging galaxies. The Lya emission appears to be extended due to the extended distribution of sources and gas. During the close encounter of galaxy progenitors when the star formation rate peaks at ~ 10^3 Msun/yr, our model produces Lya blobs with luminosity of L ~ 10^42-10^44 erg/s, and size of D ~ 10-20 kpc at z>6 and D ~ 20-50 kpc at z ~ 3, in broad agreement with observations in the same redshift range. Our results suggest that merging galaxies may produce some typical LABs as observed, but the giant ones may be produced by mergers more massive than those in our model, or a combination of mergers and cold accretion from filaments on a large scale.

preprint2013arXiv

Luminous supernova-like UV/optical/infrared transients associated with ultra-long gamma-ray bursts from metal-poor blue supergiants

Metal-poor massive stars may typically end up their lives as blue supergiants (BSGs). Gamma-ray bursts (GRBs) from such progenitors could have ultra-long duration of relativistic jets. For example Population III (Pop III) GRBs at z ~ 10-20 might be observable as X-ray rich events with a typical duration of T_90 ~ 10^4(1+z) sec. Recent GRB111209A at z = 0.677 has an ultra long duration of T_90 ~ 2.5*10^4 sec so that it have been suggested that the progenitor might be a metal-poor BSGs in the local universe. Here, we suggest luminous UV/optical/infrared emissions associated with such a new class of GRB from metal poor BSGs. Before the jet head breaks out the progenitor envelope, the energy injected by the jet is stored in a hot-plasma cocoon, which finally emerges and expands as a baryon-loaded fireball. We show that the photospheric emissions from the cocoon fireball could be intrinsically very bright (L_peak ~ 10^(42-44) erg/sec) in UV/optical bands (E_peak ~ 10 eV) with a typical duration of ~ 100 days in the rest frame. Such cocoon emissions from Pop III GRB might be detectable in infrared bands at ~ years after Pop III GRBs at up to z ~ 15 by up-coming facilities like JWST. We also suggest that GRB111209A might have been rebrightening in UV/optical bands up to an AB magnitude of < 26. The cocoon emissions from local metal-poor BSGs might have been already observed as luminous supernovae without GRB since they can be seen from the off-axis direction of the jet.

preprint2012arXiv

ART^2 : Coupling Lyman-alpha Line and Multi-wavelength Continuum Radiative Transfer

Narrow-band Lya line and broad-band continuum have played important roles in the discovery of high-redshift galaxies in recent years. Hence, it is crucial to study the radiative transfer of both Lya and continuum photons in the context of galaxy formation and evolution in order to understand the nature of distant galaxies. Here, we present a three-dimensional Monte Carlo radiative transfer code, All-wavelength Radiative Transfer with Adaptive Refinement Tree (ART^2), which couples Lya line and multi-wavelength continuum, for the study of panchromatic properties of galaxies and interstellar medium. This code is based on the original version of Li et al., and features three essential modules: continuum emission from X-ray to radio, Lya emission from both recombination and collisional excitation, and ionization of neutral hydrogen. The coupling of these three modules, together with an adaptive refinement grid, enables a self-consistent and accurate calculation of the Lya properties. As an example, we apply ART^2 to a cosmological simulation that includes both star formation and black hole growth, and study in detail a sample of massive galaxies at redshifts z=3.1 - 10.2. We find that these galaxies are Lya emitters (LAEs), whose Lya emission traces the dense gas region, and that their Lya lines show a shape characteristic of gas inflow. Furthermore, the Lya properties, including photon escape fraction, emergent luminosity, and equivalent width, change with time and environment. Our results suggest that LAEs evolve with redshift, and that early LAEs such as the most distant one detected at z ~ 8.6 may be dwarf galaxies with a high star formation rate fueled by infall of cold gas, and a low Lya escape fraction.

preprint2012arXiv

Effect of radiative transfer on damped Lyman-alpha and Lyman limit systems in cosmological SPH simulations

We study the effect of local stellar radiation and UVB on the physical properties of DLAs and LLSs at z=3 using cosmological SPH simulations. We post-process our simulations with the ART code for radiative transfer of local stellar radiation and UVB. We find that the DLA and LLS cross sections are significantly reduced by the UVB, whereas the local stellar radiation does not affect them very much except in the low-mass halos. This is because clumpy high-density clouds near young star clusters effectively absorb most of the ionizing photons from young stars. We also find that the UVB model with a simple density threshold for self-shielding effect can reproduce the observed column density distribution function of DLAs and LLSs very well, and we validate this model by direct radiative transfer calculations of stellar radiation and UVB with high angular resolution. We show that, with a self-shielding treatment, the DLAs have an extended distribution around star-forming regions typically on ~ 10-30 kpc scales, and LLSs are surrounding DLAs on ~ 30-60 kpc scales. Our simulations suggest that the median properties of DLA host haloes are: Mh = 2.4*10^10 Msun, SFR = 0.3 Msun/yr, M* = 2.4*10^8 Msun, and Z/Zsun = 0.1. About 30 per cent of DLAs are hosted by haloes having SFR = 1 - 20 Msun/yr, which is the typical SFR range for LBGs. More than half of DLAs are hosted by the LBGs that are fainter than the current observational limit. Our results suggest that fractional contribution to LLSs from lower mass haloes is greater than for DLAs. Therefore the median values of LLS host haloes are somewhat lower with Mh = 9.6*10^9 Msun, SFR = 0.06 Msun/yr, M* = 6.5*10^7 Msun and Z/Zsun = 0.08. About 80 per cent of total LLS cross section are hosted by haloes with SFR < 1 Msun/yr, hence most LLSs are associated with low-mass halos with faint LBGs below the current detection limit.

preprint2012arXiv

Were progenitors of local L* galaxies Lyman-alpha emitters at high redshift?

The Lya emission has been observed from galaxies over a redshift span z ~ 0 - 8.6. However, the evolution of high-redshift Lya emitters (LAEs), and the link between these populations and local galaxies, remain poorly understood. Here, we investigate the Lya properties of progenitors of a local L* galaxy by combining cosmological hydrodynamic simulations with three-dimensional radiative transfer calculations using the new ART^2 code. We find that the main progenitor (the most massive one) of a Milky Way-like galaxy has a number of Lya properties close to those of observed LAEs at z ~ 2 - 6, but most of the fainter ones appear to fall below the detection limits of current surveys. The Lya photon escape fraction depends sensitively on a number of physical properties of the galaxy, such as mass, star formation rate, and metallicity, as well as galaxy morphology and orientation. Moreover, we find that high-redshift LAEs show blue-shifted Lya line profiles characteristic of gas inflow, and that the Lya emission by excitation cooling increases with redshift, and becomes dominant at z > 6. Our results suggest that some observed LAEs at z ~ 2-6 with luminosity of L_Lya ~ 10^{42-43} ergs/s may be similar to the main progenitor of the Milky Way at high redshift, and that they may evolve into present-day L* galaxies.

preprint2011arXiv

Sub-millimeter brightness of early star-forming galaxies

Based on a three-dimensional model of an early star-forming galaxy, we explore the evolution of the sub-millimeter brightness. The model galaxy is employed from an ultra-high-resolution chemodynamic simulation of a primordial galaxy by Mori & Umemura, where the SFR is ~ 10 Msun/yr at t<0.3 Gyr and several Msun/yr at t>0.3 Gyr. The former phase well reproduces the observed properties of LAEs and the latter does LBGs. We solve the three-dimensional radiative transfer in the clumpy interstellar media in this model galaxy, taking the size distributions of dust grains into account, and calculate the dust temperature as a function of galactic evolutionary time. We find that the clumpiness of interstellar media plays an important role for the sub-millimeter brightness. In the LAE phase, dust grains are concentrated on clumpy star-forming regions that are distributed all over the galaxy, and the grains can effectively absorb UV radiation from stars. As a result, the dust is heated up to T>35 K. In the LBG phase, the continuous supernovae drive dust grains far away from star-forming regions. Then, the grains cannot absorb much radiation from stars, and becomes into a cold state close to the CMB temperature. Consequently, the dust temperature decreases with the evolutionary time, where the mass-weighted mean temperature is T=26 K at t=0.1 Gyr and T=21 K at t=1.0 Gyr. By this analysis, it turns out that the sub-millimeter brightness is higher in the LAE phase than that in the LBG phase, although the dust-to-gas ratio increases monotonically as a function of time. We derive the spectral energy distributions by placing the model galaxy at a given redshift. The peak flux at 850 micron is found to be S_850 ~ 0.2 - 0.9 mJy if the model galaxy is placed at 6>z>2. This means that ALMA can detect an early star-forming galaxy with SFR of ~ 10 Msun/yr by less than one hour integration with 16 antennas.

preprint2010arXiv

Effects of Ultraviolet Background and Local Stellar Radiation on the H_I Column Density Distribution

We study the impact of ultraviolet background (UVB) radiation field and the local stellar radiation on the H_I column density distribution f(N_HI) of damped Ly-alpha systems (DLAs) and sub-DLAs at z=3 using cosmological smoothed particle hydrodynamics simulations. We find that, in the previous simulations with an optically thin approximation, the UVB was sinking into the H_I cloud too deeply, and therefore we underestimated the f(N_HI) at 19 < log(N_HI) < 21.2 compared to the observations. If the UVB is shut off in the high-density regions with n_gas > 6 x 10^{-3} cm^{-3}, then we reproduce the observed f(N_HI) at z=3 very well. We also investigate the effect of local stellar radiation by post-processing our simulation with a radiative transfer code, and find that the local stellar radiation does not change the f(N_HI) very much. Our results show that the shape of f(N_HI) is determined primarily by the UVB with a much weaker effect by the local stellar radiation and that the optically thin approximation often used in cosmological simulation is inadequate to properly treat the ionization structure of neutral gas in and out of DLAs. Our result also indicates that the DLA gas is closely related to the transition region from optically-thick neutral gas to optically-thin ionized gas within dark matter halos.

preprint2010arXiv

Escape fraction of ionizing photons from high-redshift galaxies in cosmological SPH simulations

Combing the three-dimensional radiative transfer (RT) calculation and cosmological SPH simulations, we study the escape fraction of ionizing photons (f_esc) of high-redshift galaxies at z=3-6. Our simulations cover the halo mass range of M_h = 10^9 - 10^12 M_sun. We postprocess several hundred simulated galaxies with the Authentic Radiative Transfer (ART) code to study the halo mass dependence of f_esc. In this paper, we restrict ourselves to the transfer of stellar radiation from local stellar population in each dark matter halo. We find that the average f_esc steeply decreases as the halo mass increases, with a large scatter for the lower mass haloes. The low mass haloes with M_h ~ 10^9 M_sun have large values of f_esc (with an average of ~ 0.4), whereas the massive haloes with M_h ~ 10^11 M_sun show small values of f_esc (with an average of ~ 0.07). This is because in our simulations, the massive haloes show more clumpy structure in gas distribution, and star-forming regions are embedded inside these clumps, making it more difficult for the ionizing photons to escape. On the other hand, in low mass haloes, there are often conical regions of highly ionized gas due to the shifted location of young star clusters from the center of dark matter halo, which allows the ionizing photons to escape more easily than in the high-mass haloes. By counting the number of escaped ionizing photons, we show that the star-forming galaxies can ionize the intergalactic medium at z=3-6. The main contributor to the ionizing photons is the haloes with M_h < 10^10 M_sun owing to their high f_esc. The large dispersion in f_esc suggests that there may be various sizes of H{\sc ii} bubbles around the haloes even with the same mass in the early stages of reionization. We also examine the effect of UV background radiation field on f_esc using simple, four different treatment of UV background.

preprint2009arXiv

The escape of ionizing photons from supernova-dominated primordial galaxies

In order to assess the contribution of Lyman break galaxies (LBGs) and Lyman alpha emitters (LAEs) at redshifts 3<z<7 to the ionization of intergalactic medium (IGM), we investigate the escape fractions of ionizing photons from supernova-dominated primordial galaxies by solving the three-dimensional radiative transfer. The model galaxy is employed from an ultra-high-resolution chemodynamic simulation of a primordial galaxy by Mori & Umemura (2006),which well reproduces the observed properties of LAEs and LBGs. The total mass of model galaxy is 10^{11}M_sun. We solve not only photo-ionization but also collisional ionization by shocks. In addition, according to the chemical enrichment, we incorporate the effect of dust extinction, taking the size distributions of dust into account. As a result, we find that dust extinction reduces the escape fractions by a factor 1.5-8.5 in the LAE phase and by a factor 2.5-11 in the LBG phase, while the collisional ionization by shocks increases the escape fractions by a factor ~ 2. The resultant escape fractions are 0.07-0.47 in the LAE phase and 0.06-0.17 in the LBG phase. These results are well concordant with the recent estimations derived from the flux ratio at 1500 angstrom to 900 angstrom of LAEs and LBGs. Combining the resultant escape fractions with the luminosity functions of LAEs and LBGs, we find that high-z LAEs and LBGs can ionize the IGM at z=3-5. However, ionizing radiation from LAEs as well as LBGs falls short to ionize the IGM at z>6. That implies that additional ionization sources may required at z>6.