BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Kentaro Nagamine

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

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-ph.HEastro-phastro-ph.SRastro-ph.IM

Catalog footprint

What is connected

44works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

44 published item(s)

preprint2022arXiv

Dust diffusion in SPH simulations of an isolated galaxy

We compute the evolution of the grain size distribution (GSD) in a suite of numerical simulations of an isolated Milky-Way-like galaxy using the $N$-body/smoothed-particle-hydrodynamics code {\sc Gadget-4}. The full GSD is sampled on a logarithmically spaced grid with 30 bins, and its evolution is calculated self-consistently with the hydrodynamical and chemical evolution of the galaxy using a state-of-the-art star formation and feedback model. In previous versions of this model, the GSD tended to be slightly biased towards larger grains and the extinction curve had a tendency to be flatter than the observations. This work addresses these issues by considering the diffusion of dust and metals through turbulence on subgrid scales and introducing a multi-phase subgrid model that enables a smoother transition from diffuse to dense gas. We show that diffusion can significantly enhance the production of small grains and improve the agreement with the observed dust extinction curve in the Milky Way.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Osaka Feedback Model II: Modeling Supernova Feedback Based on High-Resolution Simulations

Feedback from supernovae (SNe) is an essential mechanism that self-regulates the growth of galaxies, and a better model of SN feedback is still needed in galaxy formation simulations. In the first part of this paper, using an Eulerian hydrodynamic code Athena++, we find universal scaling relations for the time evolution of momentum and radius for a superbubble, when the momentum and time are scaled by those at the shell-formation time. In the second part of this paper, we develop an SN feedback model based on the Athena++ simulation results utilizing Voronoi tessellation around each star particle, and implement it into the GADGET3-Osaka smoothed particle hydrodynamic code. Our feedback model was demonstrated to be isotropic and conservative in terms of energy and momentum. We examined the mass/energy/metal loading factors and find that our stochastic thermal feedback model produced galactic outflow that carries metals high above the galactic plane but with weak suppression of star formation. Additional mechanical feedback further suppressed star formation and brought the simulation results in better agreement with the observations of the Kennicutt--Schmidt relation, with all the results being within the uncertainties of observed data. We argue that both thermal and mechanical feedback are necessary for the SN feedback model of galaxy evolution when an individual SN bubble is unresolved.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The co-evolution of molecular hydrogen and the grain size distribution in an isolated galaxy

Understanding the evolution of dust and molecular hydrogen (H$_2$) is a critical aspect of galaxy evolution, as they affect star formation and the spectral energy distribution of galaxies. We use the $N$-body/smoothed-particle-hydrodynamics code {\sc Gadget-4} to compute the evolution of dust and H$_2$ in a suite of numerical simulations of an isolated Milky-Way-like galaxy. The evolution of the full grain size distribution (GSD) is solved by sampling the grain size on a logarithmically spaced grid with 30 bins. The evolution of a primordial chemistry network with twelve species is solved consistently with the hydrodynamic evolution of the system, including star formation, metal and energy ejections from stars into the interstellar medium through supernova feedback and stellar winds. The formation model for H$_2$ considers the GSD and photo-dissociation through the UV radiation of young stars. We identify the processes needed for producing a sizeable amount of H$_2$, verify that the resulting star formation law in the later stages of galaxy evolution is consistent with observations of local spirals, and show that our model manages to produce a galactic molecular gas fraction in line with observations of Milky-Way-like galaxies. We stress the importance of the co-evolution of the GSD and H$_2$, as models assuming a fixed MRN shape for the GSD overestimate the production of H$_2$ in regimes where the dust abundance is dominated by large grains and underestimate it in the regime where the dust is dominated by small grains, both of which are realized in simulations of dust evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

3D Distribution Map of HI Gas and Galaxies Around an Enormous Ly$α$ Nebula and Three QSOs at $z=2.3$ Revealed by the HI Tomographic Mapping Technique

We present an IGM HI tomography map in a survey volume of $16 \times 19 \times 131 \ h^{-3} {\rm comoving \ Mpc}^{3}$ (cMpc$^3$) centered at MAMMOTH-1 nebula and three neighbouring quasars at $z=2.3$. MAMMOTH-1 nebula is an enormous Ly$α$ nebula (ELAN), hosted by a type-II quasar dubbed MAMMOTH1-QSO, that extends over $1\ h^{-1}$ cMpc with not fully clear physical origin. Here we investigate the HI-gas distribution around MAMMOTH1-QSO with the ELAN and three neighbouring type-I quasars, making the IGM HI tomography map with a spatial resolution of $2.6\ h^{-1}$ cMpc. Our HI tomography map is reconstructed with HI Ly$α$ forest absorption of bright background objects at $z=2.4-2.9$: one eBOSS quasar and 16 Keck/LRIS galaxy spectra. We estimate the radial profile of HI flux overdensity for MAMMOTH1-QSO, and find that MAMMOTH1-QSO resides in a volume with significantly weak HI absorption. This suggests that MAMMOTH1-QSO has a proximity zone where quasar illuminates and photo-ionizes the surrounding HI gas and suppresses HI absorption, and that the ELAN is probably a photo-ionized cloud embedded in the cosmic web. The HI radial profile of MAMMOTH1-QSO is very similar to those of three neighbouring type-I quasars at $z=2.3$, which is compatible with the AGN unification model. We compare the distributions of the HI absorption and star-forming galaxies in our survey volume, and identify a spatial offset between density peaks of star-forming galaxies and HI gas. This segregation may suggest anisotropic UV background radiation created by star-forming galaxy density fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Direct collapse to supermassive black hole seeds: the critical conditions for suppression of $\rm H_2$ cooling

Observations of high-redshift quasars imply the presence of supermassive black holes already at z~ 7.5. An appealing and promising pathway to their formation is the direct collapse scenario of a primordial gas in atomic-cooling haloes at z ~ 10 - 20, when the $\rm H_2$ formation is inhibited by a strong background radiation field, whose intensity exceeds a critical value, $J_{\rm crit}$. To estimate $J_{\rm crit}$, typically, studies have assumed idealized spectra, with a fixed ratio of $\rm H_{2}$ photo-dissociation rate $k_{\rm H_2}$ to the $\rm H^-$ photo-detachment rate $k_{\rm H^-}$. This assumption, however, could be too narrow in scope as the nature of the background radiation field is not known precisely. In this work we argue that the critical condition for suppressing the $\rm H_2$ cooling in the collapsing gas could be described in a more general way by a combination of $k_{\rm H_2}$ and $k_{\rm H^-}$ parameters. By performing a series of cosmological zoom-in simulations with an encompassing set of $k_{\rm H_2}$ and $k_{\rm H^-}$, we examine the gas flow by following evolution of basic parameters of the accretion flow. We test under what conditions the gas evolution is dominated by $\rm H_{2}$ and/or atomic cooling. We confirm the existence of a critical curve in the $k_{\rm H_2}-k_{\rm H^-}$ plane, and provide an analytical fit to it. This curve depends on the conditions in the direct collapse, and reveals domains where the atomic cooling dominates over the molecular cooling. Furthermore, we have considered the effect of $\rm H_{2}$ self-shielding on the critical curve, by adopting three methods for the effective column density approximation in $\rm H_{2}$. We find that the estimate of the characteristic length-scale for shielding can be improved by using $λ_{\rm Jeans25}$, which is 0.25 times that of the local Jeans length.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Statistical correlation between the distribution of Ly$α$ emitters and IGM HI at $z\sim2.2$ mapped by Subaru/Hyper Suprime-Cam

The correlation between neutral Hydrogen (HI) in the intergalactic medium (IGM) and galaxies now attracts great interests. We select four fields which include several coherently strong Ly$α$ absorption systems at $z\sim2.2$ detected by using background quasars from the whole SDSS/(e)BOSS database. Deep narrow-band and $g$-band imaging are performed using the Hyper Suprime-Cam on the Subaru Telescope. We select out 2,642 Ly$α$ emitter (LAE) candidates at $z=2.177\pm0.023$ down to the Ly$α$ luminosity of $L_{\text{Ly}α}\approx 2 \times 10^{42} {\rm erg~s}^{-1}$ to construct the galaxy overdensity maps, covering an effective area of 5.39 deg$^2$. Combining the sample with the Ly$α$ absorption estimated from 64 (e)BOSS quasar spectra, we find a moderate to strong correlation between the LAE overdensity $δ_{\rm LAE}$ and the effective optical depth $τ_{\rm LoS}$ in line-of-sights, with $P$-value$=0.09\%$ ($<0.01\%$) when the field that contains a significant quasar overdensity is in(ex)cluded. The cross-correlation analysis also clearly suggests that up to $4\pm1$ pMpc, LAEs tend to cluster in the regions rich in HI gas, indicated by the high $τ_{\rm LoS}$, and avoid the low $τ_{\rm LoS}$ region where the HI gas is deficient. By averaging the $τ_{\rm LoS}$ as a function of the projected distance ($d$) to LAEs, we find a $30\%$ excess signal at $2σ$ level at $d<200$ pkpc, indicating the dense HI in circumgalactic medium, and a tentative excess at $400<d<600$ pkpc in IGM regime, corroborating the cross-correlation signal detected at about $0.5$ pMpc. These statistical analyses indicate that galaxy$-$IGM HI correlations exist on scales ranging from several hundred pkpc to several pMpc at $z\sim2.2$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The AGORA high-resolution galaxy simulations comparison project: Public data release

As part of the AGORA High-resolution Galaxy Simulations Comparison Project (Kim et al. 2014, 2016) we have generated a suite of isolated Milky Way-mass galaxy simulations using 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical galaxy formation community. In these simulations we adopted identical galactic disk initial conditions, and common physics models (e.g., radiative cooling and ultraviolet background by a standardized package). Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production were carefully constrained. Here we release the simulation data to be freely used by the community. In this release we include the disk snapshots at 0 and 500Myr of evolution per each code as used in Kim et al. (2016), from simulations with and without star formation and feedback. We encourage any member of the numerical galaxy formation community to make use of these resources for their research - for example, compare their own simulations with the AGORA galaxies, with the common analysis yt scripts used to obtain the plots shown in our papers, also available in this release.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Galaxy simulation with the evolution of grain size distribution

We compute the evolution of interstellar dust in a hydrodynamic simulation of an isolated disc galaxy. We newly implement the evolution of full grain size distribution by sampling 32 grid points on the axis of the grain radius. We solve it consistently with the chemical enrichment and hydrodynamic evolution of the galaxy. This enables us to theoretically investigate spatially resolved evolution of grain size distribution in a galaxy. The grain size distribution evolves from a large-grain-dominated ($\gtrsim 0.1~μ$m) phase to a small-grain production phase, eventually converging to a power-law-like grain size distribution similar to the so-called MRN distribution. We find that the small-grain abundance is higher in the dense ISM in the early epoch ($t\lesssim 1$ Gyr) because of efficient dust growth by accretion, while coagulation makes the small-grain abundance less enhanced in the dense ISM later. This leads to steeper extinction curves in the dense ISM than in the diffuse ISM in the early phase, while they show the opposite trend later. The radial trend is also described by faster evolution in the inner part. We also confirm that the simulation reproduces the observed trend in the relation between dust-to-gas ratio and metallicity, and in the radial gradients of dust-to-gas ratio and dust-to-metal ratio. Since the above change in the grain size distribution occurs at $t\sim 1$ Gyr, the age and density dependence of grain size distribution has a significant impact on the extinction curves even at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Connecting the Dots: Tracking Galaxy Evolution Using Constant Cumulative Number Density at 3<z<7

Using the cosmological smoothed particle hydrodynamical code GADGET-3 we make a realistic assessment of the technique of using constant cumulative number density as a tracer of galaxy evolution at high redshift. We find that over a redshift range of $3\leq z \leq7$ one can on average track the growth of the stellar mass of a population of galaxies selected from the same cumulative number density bin to within $\sim 0.20$ dex. Over the stellar mass range we probe ($10^{10.39}\leq M_s/M_\odot \leq 10^{10.75}$ at $z =$ 3 and $10^{8.48}\leq M_s/M_\odot \leq 10^{9.55}$ at $z =$ 7) one can reduce this bias by selecting galaxies based on an evolving cumulative number density. We find the cumulative number density evolution exhibits a trend towards higher values which can be quantified by simple linear formulations going as $-0.10Δz$ for descendants and $0.12Δz$ for progenitors. Utilizing such an evolving cumulative number density increases the accuracy of descendant/progenitor tracking by a factor of $\sim2$. This result is in excellent agreement, within $0.10$ dex, with abundance matching results over the same redshift range. However, we find that our more realistic cosmological hydrodynamic simulations produce a much larger scatter in descendant/progenitor stellar masses than previous studies, particularly when tracking progenitors. This large scatter makes the application of either the constant cumulative number density or evolving cumulative number density technique limited to average stellar masses of populations only, as the diverse mass assembly histories caused by stochastic physical processes such as gas accretion, mergers, and star formation of individual galaxies will lead to a larger scatter in other physical properties such as metallicity and star-formation rate.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Direct Collapse to Supermassive Black Hole Seeds: Comparing the AMR and SPH Approaches

We provide detailed comparison between the AMR code Enzo-2.4 and the SPH/N- body code GADGET-3 in the context of isolated or cosmological direct baryonic collapse within dark matter (DM) halos to form supermassive black holes. Gas flow is examined by following evolution of basic parameters of accretion flows. Both codes show an overall agreement in the general features of the collapse, however, many subtle differences exist. For isolated models, the codes increase their spatial and mass resolutions at different pace, which leads to substantially earlier collapse in SPH than in AMR cases due to higher gravitational resolution in GADGET-3. In cosmological runs, the AMR develops a slightly higher baryonic resolution than SPH during halo growth via cold accretion permeated by mergers. Still, both codes agree in the buildup of DM and baryonic structures. However, with the onset of collapse, this difference in mass and spatial resolution is amplified, so evolution of SPH models begins to lag behind. Such a delay can have effect on formation/destruction rate of H2 due to UV background, and on basic properties of host halos. Finally, isolated non-cosmological models in spinning halos, with spin parameter λ ~ 0.01 - 0.07, show delayed collapse for greater λ, but pace of this increase is faster for AMR. Within our simulation setup, GADGET-3 requires significantly larger computational resources than Enzo- 2.4 during collapse, and needs similar resources, during the pre-collapse, cosmological structure formation phase. Yet it benefits from substantially higher gravitational force and hydrodynamic resolutions, except at the end of collapse.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Evolution of binary seeds in collapsing protostellar gas clouds

We perform three dimensional smoothed particle hydrodynamics (SPH) simulations of gas accretion onto the seeds of binary stars to investigate their short-term evolution. Our simulation setup is more realistic compared to the previous works by taking into account of dynamically evolving envelope with non-uniform distribution of gas density and angular momentum of accreting flow. Our initial condition includes a seed binary and a surrounding gas envelope, modelling the phase of core collapse of gas cloud when the fragmentation has already occurred. We assume that the seed binary has no eccentricity and no growth by gas accretion. The envelope is assumed to be an isothermal gas with no self-gravity. We run multiple simulations with different values of initial mass ratio $q_0$ (the ratio of secondary over primary mass) and gas temperature, and find a critical value of $q_{\rm c} = 0.25$ which distinguishes the later evolution of mass ratio $q$ as a function of time. If $q_0$ >~ $q_{\rm c}$, the secondary seed grows faster and $q$ increases monotonically towards unity. If $q_0$ <~ $q_{\rm c}$, on the other hand, the primary seed grows faster and $q$ is lower than $q_0$ at the end of the simulation. Based on our numerical results, we analytically calculate the long-term evolution of the seed binary including the growth of binary by gas accretion. We find that the seed binary with $q_0$ >~ $q_{\rm c}$ evolves towards an equal-mass binary star, and that with $q_0$ <~ $q_{\rm c}$ evolves to a binary with an extreme value of $q$. Binary separation is a monotonically increasing function of time for any $q_0$, suggesting that the binary growth by accretion does not lead to the formation of close binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Metallicity Measurements of Gamma-Ray Burst and Supernova Explosion Sites: Lessons from HII regions in M31

We examine how the small-scale ($<$\,kpc) variation of metallicity within a galaxy, which is found in nearby galaxies, affect the observational estimates of metallicity in the explosion sites of transient events such as core-collapse supernovae (CC SNe) and gamma-ray bursts (GRBs). Assuming the same luminosity, metallicity, and spatial distributions of \ion{H}{ii}\ regions (hereafter HIIRs) as observed in M31, we compute the apparent metallicities that we would obtain when the spectrum of a target region is blended with those of surrounding HIIRs within the length scale of typical spatial resolution. When the spatial resolution of spectroscopy is $\lesssim$ 0.5\,kpc, which is typical for the existing studies of CC SN sites, we find that the apparent metallicities reflect the metallicities of target regions, but with significant systematic uncertainties in some cases. When the spatial resolution is $\gtrsim$ 1.0 kpc, regardless of the target regions (which has a wide range of metallicity that spans $\sim$ 0.6 dex for the M31 HIIRs), we always obtain the apparent metallicities similar to the average metallicity of the M31 HIIRs. Given that the apparent metallicities measured with $\gtrsim$\,kpc scale resolution do not necessarily reflect the immediate environment of the stellar explosions, the current observational estimates of high metallicities for some of the long GRB host galaxies do {\em not} rule out the hypothesis that the long GRBs are exclusively born in a low-metallicity environment.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The rise and fall of a challenger: the Bullet Cluster in $Λ$ Cold Dark Matter simulations

The Bullet Cluster has provided some of the best evidence for the $Λ$ cold dark matter ($Λ\rm{CDM}$) model via direct empirical proof of the existence of collisionless dark matter, while posing a serious challenge owing to the unusually high inferred pairwise velocities of its progenitor clusters. Here we investigate the probability of finding such a high-velocity pair in large-volume N-body simulations, particularly focusing on differences between halo finding algorithms. We find that algorithms that do not account for the kinematics of infalling groups yield vastly different statistics and probabilities. When employing the ROCKSTAR halo finder that considers particle velocities, we find numerous Bullet-like pair candidates that closely match not only the high pairwise velocity, but also the mass, mass ratio, separation distance, and collision angle of the initial conditions that have been shown to produce the Bullet Cluster in non-cosmological hydrodynamic simulations. The probability of finding a high pairwise velocity pair among haloes with $M_{\rm halo}\geq10^{14} M_{\odot}$ is $4.6\times 10^{-4}$ using ROCKSTAR, while it is $\approx 34\times$ lower using a friends-of-friends (FOF) based approach as in previous studies. This is because the typical spatial extent of Bullet progenitors is such that FOF tends to group them into a single halo despite clearly distinct kinematics. Further requiring an appropriately high average mass among the two progenitors, we find the comoving number density of potential Bullet-like candidates to be on the order of $\approx10^{-10} \,{\rm Mpc}^{-3}$. Our findings suggest that $Λ\rm{CDM}$ straightforwardly produces massive, high relative velocity halo pairs analogous to Bullet Cluster progenitors, and hence the Bullet Cluster does not present a challenge to the $Λ\rm{CDM}$ model.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Impact of H_2-based star formation model on the z>=6 luminosity function and the ionizing photon budget for reionization

We present the results of a numerical study examining the effect of H_2-based star formation (SF) model on the rest-frame UV luminosity function and star formation rate function (SFRF) of z>=6 galaxies, and the implications for reionization. Using cosmological hydrodynamical simulations outfitted with an H_2-SF model, we find good agreement with our previous results (non-H_2 SF model) and observations at Muv<=-18. However at Muv>-18, we find that the LF deviates from both our previous work and current observational extrapolations, producing significantly fewer low-luminosity galaxies and exhibiting additional turnover at the faint end. We constrain the redshift evolution of this turnover point using a modified Schechter function that includes additional terms to quantify the turnover magnitude (Muv^t) and subsequent slope (Beta). We find that Muv^t evolves from Muv^t=-17.33 (at z=8) to -15.38 (z=6), while Beta becomes shallower by ΔBeta=0.22 during the same epoch. This occurs in an Muv range which will be observable by JWST. By integrating the SFRF, we determine that even though H_2-SF model significantly reduces the number density of low luminosity galaxies at Muv>-18, it does not suppress the total SFR density enough to affect the capability of SF to maintain reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Molecular hydrogen regulated star formation in cosmological SPH simulations

It has been shown observationally that star formation (SF) correlates tightly with the presence of molecular hydrogen (H$_2$). Therefore it would be important to investigate its implication on galaxy formation in a cosmological context. In the present work, we track the H$_2$ mass fraction within our cosmological smoothed particle hydrodynamics (SPH) code GADGET-3 using an equilibrium analytic model by Krumholz et al. This model allows us to regulate the star formation in our simulation by the local abundance of H$_2$ rather than the total cold gas density, and naturally introduce the dependence of star formation on metallicity. We investigate implications of the equilibrium H$_2$-based SF model on galaxy population properties, such as the stellar-to-halo mass ratio (SHMR), baryon fraction, cosmic star formation rate density (SFRD), galaxy specific SFR, galaxy stellar mass functions (GSMF), and Kennicutt-Schmidt (KS) relationship. The advantage of our work over the previous ones is having a large sample of simulated galaxies in a cosmological volume from high-redshift to $z=0$. We find that low-mass halos with M$_{DM}<10^{10.5}$ M$_\odot$ are less efficient in producing stars in the H$_2$-based SF model at $z\geq6$, which brings the simulations to a better agreement with observational estimates of SHMR and GSMF at the low-mass end. This is particularly evident by a reduction in the number of low-mass galaxies at M$_\star\leq10^{8}$ M$_\odot$ in the GSMF. The overall SFRD is also reduced at high-$z$ in the H$_2$ run, which results in slightly higher SFRD at low-redshift due to more abundant gas available for star formation at later times. This new H$_2$ model is able to reproduce the empirical KS relationship at $z=0$ naturally without the need for setting its normalization by hand, and overall it seems to have more advantages than the previous pressure-based SF model.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the Virialization of Disk Winds: Implications for the Black Hole Mass Estimates in AGN

Estimating the mass of a supermassive black hole (SMBH) in an active galactic nucleus (AGN) usually relies on the assumption that the broad line region (BLR) is virialized. However, this assumption seems invalid in BLR models that consists of an accretion disk and its wind. The disk is likely Keplerian and therefore virialized. However, the wind material must, beyond a certain point, be dominated by an outward force that is stronger than gravity. Here, we analyze hydrodynamic simulations of four different disk winds: an isothermal wind, a thermal wind from an X-ray heated disk, and two line-driven winds, one with and the other without X-ray heating and cooling. For each model, we check whether gravity governs the flow properties, by computing and analyzing the volume-integrated quantities that appear in the virial theorem: internal, kinetic, and gravitational energies, We find that in the first two models, the winds are non-virialized whereas the two line-driven disk winds are virialized up to a relatively large distance. The line-driven winds are virialized because they accelerate slowly so that the rotational velocity is dominant and the wind base is very dense. For the two virialized winds, the so-called projected virial factor scales with inclination angle as $1/ \sin^2{i}$. Finally, we demonstrate that an outflow from a Keplerian disk becomes unvirialized more slowly when it conserves the gas specific angular momentum -- as in the models considered here, than when it conserves the angular velocity -- as in the so-called magneto-centrifugal winds.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Duty Cycle and the Increasing Star Formation History of z>=6 Galaxies

We examine the duty cycle and the history of star formation (SFH) for high-redshift galaxies at z>=6 using cosmological hydrodynamic simulations. We find that, even though individual galaxies have bursty SFH, the averaged SFH between z~15 to z=6 can be characterized well by either an exponentially increasing functional form with characteristic time-scales of 70 Myr to 200 Myr for galaxies with stellar masses Ms~10^6 Msun to >10^10 Msun respectively, or by a simple power-law form which exhibits a similar mass dependent time-scales. Using the SFH of individual galaxies, we measure the duty cycle of star formation (DC_SFH); i.e., the fraction of time a galaxy of a particular mass spends above a star formation rate (SFR) threshold which would make it observable to the Hubble Space Telescope (HST) during a given epoch. We also examine the fraction of galaxies at a given redshift that are brighter than a rest-frame UV magnitude (Muv ~ -18), which is sufficient enough to make them observable (DC_Muv). We find that both DC_SFH and DC_Muv make a sharp transition from zero (for galaxies with Ms <= 10^7 Msun) to unity (for Ms > 10^9 Msun). The measured duty cycle is also manifested in the intrinsic scatter in the Ms-SFR relationship (~ 1 dex) and Ms-Muv relationship (ΔMuv ~ +-1 mag). We provide analytic fits to the DC as a function of Ms using a sigmoid function, which can be used to correct for catalogue incompleteness. We consider the effects of duty cycle to the observational estimate of galaxy stellar mass functions (GSMF) and the star formation rate density (SFRD), and find that it results in a much shallower low-mass end slopes of the GSMF and a reduction of >~ 70% of our intrinsic SFRD, making our simulation results more compatible with observational estimates.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gas Accretion onto a Supermassive Black Hole: a step to model AGN feedback

We study the gas accretion onto a supermassive black hole (SMBH) using the 3D SPH code GADGET-3 on scales of 0.1-200 pc. First we test our code with spherically symmetric, adiabatic Bondi accretion problem. We find that our simulation can reproduce the expected Bondi accretion flow very well for a limited amount of time until the effect of outer boundary starts to be visible. We also find artificial heating of gas near the inner accretion boundary due to the artificial viscosity of SPH. Second, we implement radiative cooling and heating due to X-rays, and examine the impact of thermal feedback by the central X-ray source. The accretion flow roughly follows the Bondi solution for low central X-ray luminosities, however, the flow starts to exhibit non-spherical fragmentation due to thermal instability for a certain range of central L_X, and a strong overall outflow develops for greater L_X. The cold gas develops filamentary structures that fall into the central SMBH, whereas the hot gas tries to escape through the channels in-between the cold filaments. Such fragmentation of accreting gas can assist in the formation of clouds around AGN, induce star-formation, and contribute to the observed variability of narrow-line regions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Multiphase, non-spherical gas accretion onto a black hole

(Abridged) We investigate non-spherical behavior of gas accreting onto a central supermassive black hole performing simulations using the SPH code GADGET-3 including radiative cooling and heating by the central X-ray source. As found in earlier 1D studies, our 3D simulations show that the accretion mode depends on the X-ray luminosity (L_X) for a fixed density at infinity and accretion efficiency. In the low L_X limit, gas accretes in a stable, spherically symmetric fashion. In the high L_X limit, the inner gas is significantly heated up and expands, reducing the central mass inflow rate. The expanding gas can turn into a strong enough outflow capable of expelling most of the gas at larger radii. For some intermediate L_X, the accretion flow becomes unstable developing prominent non-spherical features, the key reason for which is thermal instability (TI) as shown by our analyses. Small perturbations of the initially spherically symmetric accretion flow that is heated by the intermediate L_X quickly grow to form cold and dense clumps surrounded by overheated low density regions. The cold clumps continue their inward motion forming filamentary structures; while the hot infalling gas slows down because of buoyancy and can even start outflowing through the channels in between the filaments. We found that the ratio between the mass inflow rates of the cold and hot gas is a dynamical quantity depending on several factors: time, spatial location, and L_X; and ranges between 0 and 4. We briefly discuss astrophysical implications of such TI-driven fragmentation of accreting gas on the formation of clouds in narrow and broad line regions of AGN, the formation of stars, and the observed variability of the AGN luminiosity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gamma-ray burst rate: high-redshift excess and its possible origins

Prompted by various analyses of long (Type II) GRB rates and their relationship to the cosmic star-formation history, metallicity and luminosity function evolution, we systematically analyze these effects with a Monte Carlo code. We test various cosmic star-formation history models including analytical and empirical models as well as those derived from cosmological simulations. We also explore expressions for metallicity enhancement of the GRB rate with redshift, as presented in the literature, and discuss improvements to these analytic expressions from the point of view of galactic evolution. These are also compared to cosmological simulations on metal enrichment. Additionally we explore possible evolutionary effects of the GRB rate and luminosity function with redshift. The simulated results are tested with the observed Swift sample including the L, z, and peak flux (log N-log P) distributions. The observational data imply that an increase in the GRB rate is necessary to account for the observations at high redshift, although the form of this enhancement is unclear. A rate increase due to lower metallicity at higher redshift may not be the singular cause and is subject to a variety of uncertainties. Alternatively, evolution of the GRB luminosity function break with redshift shows promise as a possible alternative.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the inconsistency between the estimates of cosmic star formation rate and stellar mass density of high redshift galaxies

There are mainly two different approaches to measure the cosmic star formation history: direct star formation rate density (SFRD) and stellar mass density rhostar as functions of redshift. Compilations of current observations seem to show a disparity in the two quantities, in the sense that the integral of SFRD is higher than the observed rhostar (after considering gas recycling). Using cosmological smoothed particle hydrodynamics simulations based on the concordance Lambda cold dark matter model, we show that the two quantities become more consistent with each other when we consider the observed galaxy mass limit. The comparison between simulations and (dust corrected) observed cosmic SFRD shows a good agreement, while the observed rhostar is significantly lower than the simulation results. This can be reconciled if the current high-$z$ galaxy surveys are missing faint low-mass galaxies due to their flux limit. Our simulated GSMFs have steep low-mass end slopes of alpha < -2 at z>3, and when these numerous low-mass galaxies are included, the total rhostar matches with the integral of SFRD.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Pairwise Velocities of Dark Matter Halos: a Test for the Lambda Cold Dark Matter Model using the Bullet Cluster

The existence of 1E0657-56 poses a challenge to the concordance Lambda cold dark matter model. Here we investigate the velocity distribution of dark matter halo pairs in large N-body simulations with differing box sizes (250Mpc/h-2 Gpc/h) and resolutions. We examine statistics such as the halo masses, pairwise halo velocities (v_{12}), and pair separation distances. We then compare our results to the initial conditions (ICs) required to reproduce the observational properties of 1E0657-56 in non-cosmological hydrodynamical simulations. We find that the high velocity tail of the v_{12} distribution extends to greater velocities as we increase the box size. We also find that the number of high-v_{12} pairs increases as we increase the particle count and resolution with a fixed box size, however, this increase is mostly due to lower mass halos which do not match the observed masses of 1E0657-56. We find that the redshift evolution is not very strong for the v_{12} distribution function between z=0.0 and z~0.5. We identify some pairs whose v_{12} resemble the required ICs, however, even the best candidates have either wrong halo mass ratios, or too large separations. Our simulations suggest that it is very difficult to produce such ICs at z=0.0, 0.296,& 0.489 in comoving volumes as large as (2Gpc/h)^3. Based on the extrapolation of our cumulative v_{12} function, we find that one needs a simulation with a comoving box size of (4.48Gpc/h)^3 and 2240^3 DM particles in order to produce at least one pair of halos that resembles the required v_{12} and observed masses of 1E0657-56. We find that the probability of finding a halo pair with v_{12}>=3000km/s and masses >=10^{14}Msun to be 2.76x10^{-8} at z=0.489. We conclude that either 1E0657-56 is incompatible with the concordance LCDM universe, or the ICs suggested by the non-cosmological simulations must be revised to give a lower value of v_{12}.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

SPH Simulations of Black Hole Accretion: A Step to Model Black Hole Feedback in Galaxies

(Abridged) We test how accurately the smoothed particle hydrodynamics (SPH) numerical technique can follow spherically-symmetric Bondi accretion. Using the 3D SPH code GADGET-3, we perform simulations of gas accretion onto a central supermassive black hole (SMBH) of mass 10^8 M_sun within the radial range of 0.1 - 200 pc. We carry out simulations without and with radiative heating by a central X-ray corona and radiative cooling. For an adiabatic case, the radial profiles of hydrodynamical properties match the Bondi solution, except near the inner and outer radius of the computational domain. We find that adiabatic Bondi accretion can be reproduced for durations of a few dynamical times at the Bondi radius, and for longer times if the outer radius is increased. With radiative heating and cooling included, the spherically accreting gas takes a longer time to reach a steady-state than the adiabatic Bondi accretion runs, and in some cases does not reach a steady-state even within several hundred dynamical times. We find that artificial viscosity in the GADGET code causes excessive heating near the inner radius, making the thermal properties of the gas inconsistent with a physical solution. This overheating occurs typically only in the supersonic part of the flow, so that it does not affect the mass accretion rate. We see that increasing the X-ray luminosity produces a lower central mass inflow rate, implying that feedback due to radiative heating is operational in our simulations. With a sufficiently high X-ray luminosity, the inflowing gas is radiatively heated up, and an outflow develops. We conclude that the SPH simulations can capture the gas dynamics needed to study radiative feedback provided artificial viscosity alters only highly supersonic part of the inflow.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Steep Faint-end Slopes of Galaxy Mass and Luminosity Functions at z>=6 and the Implications for Reionisation

We present the results of a numerical study comparing photometric and physical properties of simulated z=6-9 galaxies to the observations taken by the WFC3 instrument aboard the Hubble Space Telescope. Using cosmological hydrodynamical simulations we find good agreement with observations in color-color space at all studied redshifts. We also find good agreement between observations and our Schechter luminosity function fit in the observable range, Muv<= -18, provided that a moderate dust extinction effect exists for massive galaxies. However beyond what currently can be observed, simulations predict a very large number of low-mass galaxies and evolving steep faint-end slopes from alpha_L = -2.15 at z=6 to alpha_L = -2.64 at z=9, with a dependence of |alpha_L| \propto (1+z)^0.59. During the same epoch, the normalization phi* increases and the characteristic magnitude Muv* becomes moderately brighter with decreasing redshift. We find similar trends for galaxy stellar mass function with evolving low-mass end slope from alpha_M = - 2.26 at z=6 to alpha_M = -2.87 at z=9, with a dependence of |alpha_M| \propto (1+z)^0.65. Together with our recent result on the high escape fraction of ionizing photons for low-mass galaxies, our results suggest that the low-mass galaxies are important contributor of ionizing photons for the reionisation of the Universe at z>=6.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Cross-Correlation between Damped Ly-alpha Systems and Ly-break Galaxies in Cosmological SPH Simulations

We calculate the cross-correlation function (CCF) between damped Ly-a systems (DLAs) and Lyman break galaxies (LBGs) using cosmological hydrodynamic simulations at z=3. We compute the CCF with two different methods. First, we assume that there is one DLA in each dark matter halo if its DLA cross section is non-zero. In our second approach we weight the pair-count by the DLA cross section of each halo, yielding a cross-section-weighted CCF. We also compute the angular CCF for direct comparison with observations. Finally, we calculate the auto-correlation functions of LBGs and DLAs, and their bias against the dark matter distribution. For these different approaches, we consistently find that there is good agreement between our simulations and observational measurements by Cooke et al. and Adelberger et al. Our results thus confirm that the spatial distribution of LBGs and DLAs can be well described within the framework of the concordance Lambda CDM model. We find that the correlation strengths of LBGs and DLAs are consistent with the actual observations, and in the case of LBGs it is higher than would be predicted by low-mass galaxy merger models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Effects of cosmological parameters and star formation models on the cosmic star formation history in LambdaCDM cosmological simulations

We investigate the effects of the change of cosmological parameters and star formation (SF) models on the cosmic SF history using cosmological smoothed particle hydrodynamics (SPH) simulations based on the cold dark matter (CDM) model. We vary the cosmological parameters within 1-sigma error from the WMAP best-fit parameters, and find that such changes in cosmological parameters mostly affect the amplitude of the cosmic SF history. At high redshift (hereafter high-z), the star formation rate (SFR) is sensitive to the cosmological parameters that control the small-scale power of the primordial power spectrum, while the cosmic matter content becomes important at lower redshifts. We also test two new SF models: 1) the `Pressure' model based on the work by Schaye & Dalla Vecchia (2008), and 2) the `Blitz' model that takes the effect of molecular hydrogen formation into account, based on the work by Blitz & Rosolowsky (2006). Compared to the previous conventional SF model, the Pressure model reduces the SFR in low-density regions and shows better agreement with the observations of the Kennicutt-Schmidt law. This model also suppresses the early star formation and shifts the peak of the cosmic SF history toward lower redshift, more consistently with the recent observational estimates of cosmic SFR density. The simulations with the new SF model also predict lower global stellar mass densities at high-z, larger populations of low-mass galaxies and a higher gas fraction in high-z galaxies. Our results suggest that there is room left in the model uncertainties to reconcile the discrepancy that was found between the theory and observations of cosmic SF history and stellar mass density. Nevertheless, our simulations still predict higher stellar mass densities than most of the observational estimates.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Infrared Emission from High-Redshift Galaxies in Cosmological SPH Simulations

We compute the infrared (IR) emission from high-redshift galaxies in cosmological smoothed particle hydrodynamics simulations by coupling the output of the simulation with the population synthesis code `GRASIL' by Silva et al. Based on the stellar mass, metallicity and formation time of each star particle, we estimate the full spectral energy distribution of each star particle from ultraviolet to IR, and compute the luminosity function of simulated galaxies in the Spitzer broadband filters for direct comparison with the available Spitzer observations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Large-Scale Outflows from AGN: a link between central black holes and galaxies

We summarize the results from numerical simulations of mass outflows from AGN. We focus on simulations of outflows driven by radiation from large-scale inflows. We discuss the properties of these outflows in the context of the so-called AGN feedback problem. Our main conclusion is that this type of outflows are efficient in removing matter but inefficient in removing energy.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Luminosity Distribution of Gamma-Ray Burst Host Galaxies at redshift z=1 in Cosmological Smoothed Particle Hydrodynamic Simulations: Implications for the Metallicity Dependence of GRBs

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Lyman-alpha Emitters and Lyman-break Galaxies at z=3-6 in Cosmological SPH Simulations

We study the properties of Lyman-alpha emitters (LAEs) and Lyman-break galaxies (LBGs) at z=3-6 using cosmological SPH simulations. We investigate two simple scenarios for explaining the observed Ly-a and rest-frame UV luminosity functions (LFs) of LAEs: (i) the "escape fraction" scenario, in which the "effective" escape fraction (including the IGM attenuation) of Ly-a photons is f_Lya ~0.1 (0.15) at z=3 (6), and (ii) the "stochastic" scenario, in which the fraction of LAEs that are turned on at z=3 (6) is \Cstoc ~0.07 (0.2) after correcting for the IGM attenuation. Our comparisons with a number of different observations suggest that the stochastic scenario is preferred over the escape fraction scenario. We find that the mean values of stellar mass, metallicity and black hole mass hosted by LAEs are all smaller in the stochastic scenario than in the escape fraction scenario. In our simulations, the galaxy stellar mass function evolves rapidly, as expected in hierarchical structure formation. However, its evolution is largely compensated by a beginning decline in the specific star formation rate, resulting in little evolution of the rest-frame UV LF from z=6 to 3. The rest-frame UV LF of both LAEs and LBGs at z=3 & 6 can be described well by the stochastic scenario provided the extinction is moderate, E(B-V) ~0.15, for both populations, although our simulation might be overpredicting the number of bright LBGs at z=6. We also discuss the correlation function and bias of LAEs. The Ly-a LFs at z=6 in a field-of-view of 0.2 deg^2 show a significantly larger scatter owing to cosmic variance relative to that in a 1 deg^2 field, and the scatter seen in the current observational estimates of the Ly-a LF can be accounted for by cosmic variance.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multicomponent and Variable Velocity Galactic Outflow in Cosmological Hydrodynamic Simulations

We develop a new ``Multicomponent and Variable Velocity'' (MVV) galactic outflow model for cosmological smoothed particle hydrodynamic (SPH) simulations. The MVV wind model reflects the fact that the wind material can arise from different phases in the interstellar medium (ISM), and the mass-loading factor in the MVV model is a function of galaxy stellar mass. We find that the simulation with the MVV outflow has the following characteristics: (i) the intergalactic medium (IGM) is hardly heated up, and the mean IGM temperature is almost the same as in the no-wind run; (ii) it has lower cosmic star formation rates (SFRs) compared to the no-wind run, but higher SFRs than the constant velocity wind run; (iii) it roughly agrees with the observed IGM metallicity, and roughly follows the observed evolution of Omega(Civ); (iv) the lower mass galaxies have larger mass-loading factors, and the low-mass end of galaxy stellar mass function is flatter than in the previous simulations. Therefore, the MVV outflow model mildly alleviates the problem of too steep galaxy stellar mass function seen in the previous SPH simulations. In summary, the new MVV outflow model shows reasonable agreement with observations, and gives better results than the constant velocity wind model.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On the Feedback Efficiency of Active Galactic Nuclei

(Abridged) We measure and analyze the energy, momentum, and mass feedback efficiencies due to radiation from AGN in relatively large scale outflows. Our measurements are based on the axisymmetric and time-dependent radiation-hydrodynamical simulations recently presented in Kurosawa & Proga. These simulations follow dynamics of gas under the influences of the AGN radiation and the gravity of the central 10^8 Msun black hole on scales from ~0.01 to ~10pc. We compare our Mdot_a-rho_o relation with that predicted by the Bondi accretion model. For high luminosities comparable to the Eddington limit, the power-law fit (Mdot_a ~ rho_o^q) to our models yields q=~0.5 instead of q=1.0 which is predicted by the Bondi model. This difference is caused by the outflows which are important for the overall mass budget at high luminosities. The maximum momentum and mass feedback efficiencies found in our models are ~10^-2 and ~10^-1, respectively. However, the outflows are much less important energetically: their thermal and kinetic powers in units of the radiative luminosity are ~10^-5 and ~10^-4, respectively. The efficiencies do not increase monotonically with the accretion luminosity but rather peak around the Eddington limit beyond which a steady state disk-wind-like solution exists. Our energy feedback efficiencies are significantly lower than 0.05, which is required in some cosmological and galaxy merger simulations. The low feedback efficiencies found here could have significant implications on the mass growth of super massive black holes in the early universe. We stress however that we have not considered the innermost parts of the accretion and outflow where radiation and matter interact most strongly. The feedback from this region could have efficiencies significantly above the low values found here.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Star Formation History of Dwarf Galaxies in Cosmological Hydrodynamic Simulations

We examine the past and current work on the star formation (SF) histories of dwarf galaxies in cosmological hydrodynamic simulations. The results obtained from different numerical methods are still somewhat mixed, but the differences are understandable if we consider the numerical and resolution effects. It remains a challenge to simulate the episodic nature of SF history in dwarf galaxies at late times within the cosmological context of a cold dark matter model. More work is needed to solve the mysteries of SF history of dwarf galaxies employing large-scale hydrodynamic simulations on the next generation of supercomputers.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Effects of metal enrichment and metal cooling in galaxy growth and cosmic star formation history

We present the results of a numerical study on the effects of metal enrichment and metal cooling on galaxy formation and cosmic star formation (SF) history using cosmological hydrodynamic simulations. We find following differences in the simulation with metal cooling when compared to the run without it: (1) the cosmic star formation rate (SFR) is enhanced by about 50 & 20% at z=1 & 3, respectively; (2) the gas mass fraction in galaxies is lower; (3) the total baryonic mass function (gas + star) at z=3 does not differ significantly, but shows an increase in the number of relatively massive galaxies at z=1; (4) the baryonic mass fraction of intergalactic medium (IGM) is reduced at z<3 due to more efficient cooling and gas accretion onto galaxies. Our results suggest that the metal cooling enhances the galaxy growth by two different mechanisms: (1) increase of SF efficiency in the local interstellar medium (ISM), and (2) increase of IGM accretion onto galaxies. The former process is effective throughout most of the cosmic history, while the latter is effective only at z<3 when the IGM is sufficiently enriched by metals owing to feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Incidence Rate of GRB-host-DLAs at High Redshift

We study the incidence rate of damped Ly-a systems associated with the host galaxies of gamma-ray bursts (GRB-host-DLAs) as functions of neutral hydrogen column density (N_HI) and projected star formation rate (SFR) using cosmological SPH simulations. Assuming that the occurrence of GRBs is correlated with the local SFR, we find that the median N_HI of GRB-host-DLAs progressively shifts to lower N_HI values with increasing redshift, and the incidence rate of GRB-host-DLAs with log N_HI > 21.0 decreases rapidly at z>=6. Our results suggest that the likelihood of observing the signature of IGM attenuation in GRB afterglows increases towards higher redshift, because it will not be blocked by the red damping wing of DLAs in the GRB host galaxies. This enhances the prospects of using high-redshift GRBs to probe the reionization history of the Universe. The overall incidence rate of GRB-host-DLAs decreases monotonically with increasing redshift, whereas that of QSO-DLAs increases up to z=6. A measurement of the difference between the two incidence rates would enable an estimation of the value of η_grb, which is the mass fraction of stars that become GRBs for a given amount of star formation. Our predictions can be tested by upcoming high-z GRB missions, including JANUS (Joint Astrophysics Nascent Universe Scout) and SVOM (Space multi-band Variable Object Monitor).

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Massive galaxies at redshift 2 in cosmological hydrodynamic simulations

We study the properties of galaxies at z=2 in a Lambda cold dark matter universe, using two different types of hydrodynamic simulation methods -- Eulerian TVD and smoothed particle hydrodynamics (SPH) -- and a spectrophotometric analysis in the U_n, G, R filter set. The simulated galaxies at z=2 satisfy the color-selection criteria proposed by Adelberger et al. (2004) and Steidel et al. (2004) when we assume Calzetti extinction with E(B-V)=0.15. We find that the number density of simulated galaxies brighter than R<25.5 at z=2 is about 1e-2 h^3 Mpc^-3 for E(B-V)=0.15, which is roughly twice that of the number density found by Erb et al. (2004) for the UV bright sample. This suggests that roughly half of the massive galaxies with M*>10^{10} Msun/h at z=2 are UV bright population, and the other half is bright in the infra-red wavelengths. The most massive galaxies at z=2 have stellar masses >= 10^{11-12} Msun. They typically have been continuously forming stars with a rate exceeding 30 Msun/yr over a few Gyrs from z=10 to z=2, together with significant contribution by starbursts reaching up to 1000 Msun/yr which lie on top of the continuous component. TVD simulations indicate a more sporadic star formation history than the SPH simulations. Our results do not imply that hierarchical galaxy formation fails to account for the observed massive galaxies at z>=1. The global star formation rate density in our simulations peaks at z>=5, a much higher redshift than predicted by the semianalytic models. This star formation history suggests early build-up of the stellar mass density, and predicts that 70 (50, 30)% of the total stellar mass at z=0 had already been formed by z=1 (2, 3). Upcoming observations by Spitzer and Swift might help to better constrain the star formation history at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Massive galaxies in cosmological simulations: UV-selected sample at redshift z=2

We study the properties of galaxies at z=2 in a Lambda CDM universe, using two different types of hydrodynamic simulation methods (Eulerian TVD and SPH) and a spectrophotometric analysis in the Un, G, R filter set. The simulated galaxies at z=2 satisfy the color-selection criteria proposed by Adelberger et al. (2004) when we assume Calzetti extinction with E(B-V)=0.15. We find that the number density of simulated galaxies brighter than R<25.5 at z=2 is about 2e-2 h^3/Mpc^3, roughly one order of magnitude larger than that of Lyman break galaxies at z=3. The most massive galaxies at z=2 have stellar masses >~1e11 Msun, and their observed-frame G-R colors lie in the range 0.0<G-R<1.0. They typically have been continuously forming stars with a rate exceeding 30 Msun/yr over a few Gyrs from z=10 to z=2, although the TVD simulation indicates a more sporadic star formation history than the SPH simulations. Of order half of their stellar mass was already assembled by z~4. The reddest massive galaxies at z=2 with G-R >= 1.0 and Mstar>1e10 Msun/h finished the build-up of their stellar mass by z~3. Interestingly, our study suggests that the majority of the most massive galaxies at z=2 should be detectable at rest-frame UV wavelengths, contrary to some recent claims made on the basis of near-IR studies of galaxies at the same epoch, provided the median extinction is less than E(B-V)<0.3. However, our results also suggest that the fraction of stellar mass contained in galaxies that pass the color-selection criteria could be as low as 50% of the total stellar mass in the Universe at z=2. Our simulations suggest that the missing stellar mass is contained in fainter (R>25.5) and intrinsically redder galaxies. Our results do not suggest that hierarchical galaxy formation fails to account for the massive galaxies at z>=1. (abridged)

0 citations0 reviewsNo code linkedOpen access