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

Jun-Hwan Choi

Jun-Hwan Choi 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-phastro-ph.HEastro-ph.IM

Catalog footprint

What is connected

22works
5topics
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

22 published item(s)

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
preprint2016arXiv

Cosmic Dawn (CoDa): the First Radiation-Hydrodynamics Simulation of Reionization and Galaxy Formation in the Local Universe

Cosmic reionization by starlight from early galaxies affected their evolution, thereby impacting reionization, itself. Star formation suppression, for example, may explain the observed underabundance of Local Group dwarfs relative to N-body predictions for Cold Dark Matter. Reionization modelling requires simulating volumes large enough [~(100Mpc)^3] to sample reionization "patchiness", while resolving millions of galaxy sources above ~10^8 Msun , combining gravitational and gas dynamics with radiative transfer. Modelling the Local Group requires initial cosmological density fluctuations pre-selected to form the well-known structures of the local universe today. Cosmic Dawn ("CoDa") is the first such fully-coupled, radiation-hydrodynamics simulation of reionization of the local universe. Our new hybrid CPU-GPU code, RAMSES-CUDATON, performs hundreds of radiative transfer and ionization rate-solver timesteps on the GPUs for each hydro-gravity timestep on the CPUs. CoDa simulated (91Mpc)^3 with 4096^3 particles and cells, to redshift 4.23, on ORNL supercomputer Titan, utilizing 8192 cores and 8192 GPUs. Global reionization ended slightly later than observed. However, a simple temporal rescaling which brings the evolution of ionized fraction into agreement with observations also reconciles ionizing flux density, cosmic star formation history, CMB electron scattering optical depth and galaxy UV luminosity function with their observed values. Photoionization heating suppressed the star formation of haloes below ~2 x 10^9 Msun , decreasing the abun- dance of faint galaxies around MAB_1600 = [-10,-12]. For most of reionization, star formation was dominated by haloes between 10^10 - 10^11 Msun , so low-mass halo suppression was not reflected by a distinct feature in the global star formation history. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Baryon Cycle at High Redshifts: Effects of Galactic Winds on Galaxy Evolution in Overdense and Average Regions

We employ high-resolution cosmological zoom-in simulations focusing on a high-sigma peak and an average cosmological field at $z\sim 6-12$, in order to investigate the influence of environment and baryonic feedback on galaxy evolution in the reionization epoch. Strong feedback, e.g., galactic winds, caused by elevated star formation rates (SFRs) is expected to play an important role in this evolution. We compare different outflow prescriptions: (i) constant wind velocity (CW), (ii) variable wind scaling with galaxy properties (VW), and (iii) no outflows (NW). The overdensity leads to accelerated evolution of dark matter and baryonic structures, absent in the "normal" region, and to shallow galaxy stellar mass functions at the low-mass end. Although CW shows little dependence on both environments, the more physically motivated VW model does exhibit this effect. In addition, VW can reproduce the observed specific SFR (sSFR) and the sSFR-stellar mass relation, which CW and NW fail to satisfy simultaneously. Winds also differ substantially in affecting the state of the intergalactic medium (IGM). The difference lies in volume-filling factor of hot, high-metallicity gas which is near unity for CW, while it remains confined in massive filaments for VW, and locked up in galaxies for NW. Such gas is nearly absent in the normal region. Although all wind models suffer from deficiencies, the VW model seems to be promising in correlating the outflow properties to those of host galaxies. Further constraints on the state of the IGM at high-$z$ are needed to separate different wind models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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

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
preprint2014arXiv

The Gentle Growth of Galaxies at High Redshifts in Overdense Environments

We have explored prevailing modes of galaxy growth for redshifts z ~ 6-14, comparing substantially overdense and normal regions of the universe, using high-resolution zoom-in cosmological simulations. Such rare overdense regions have been projected to host high-z quasars. We demonstrate that galaxies in such environments grow predominantly by a smooth accretion from cosmological filaments which dominates the mass input from major, intermediate and minor mergers. We find that by z ~6, the accumulated galaxy mass fraction from mergers falls short by a factor of 10 of the cumulative accretion mass for galaxies in the overdense regions, and by a factor of 5 in the normal environments. Moreover, the rate of the stellar mass input from mergers also lies below that of an in-situ star formation (SF) rate. The fraction of stellar masses in galaxies contributed by mergers in overdense regions is ~12%, and ~33% in the normal regions, at these redshifts. Our median SF rates for ~few X 10^9 Mo galaxies agrees well with the recently estimated rates for z ~ 7 galaxies from Spitzer's SURF-UP survey. Finally, we find that the main difference between the normal and overdense regions lies in the amplified growth of massive galaxies in massive dark matter halos. This leads to the formation of >= 10^{10} Mo galaxies due to the ~100-fold increase in mass during the above time period. Such galaxies are basically absent in the normal regions at these redshifts.

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

Supermassive Black Hole Formation at High Redshifts via Direct Collapse: Physical Processes in the Early Stage

We use numerical simulations to explore whether direct collapse can lead to the formation of SMBH seeds at high-z. We follow the evolution of gas within DM halos of 2 x 10^8 Mo and 1 kpc. We adopt cosmological density profiles and j-distributions. Our goal is to understand how the collapsing flow overcomes the centrifugal barrier and whether it is subject to fragmentation. We find that the collapse leads either to a central runaway or to off-center fragmentation. A disk-like configuration is formed inside the centrifugal barrier. For more cuspy DM distribution, the gas collapses more and experiences a bar-like perturbation and a central runaway. We have followed this inflow down to ~10^{-4} pc. The flow remains isothermal and the specific angular momentum is efficiently transferred by gravitational torques in a cascade of nested bars. This cascade supports a self-similar, disk-like collapse. In the collapsing phase, virial supersonic turbulence develops and fragmentation is damped. For larger initial DM cores the timescales become longer. In models with more organized initial rotation, a torus forms and appears to be supported by turbulent motions. The evolution depends on the competition between two timescales, corresponding to the onset of the central runaway and off-center fragmentation. For less organized rotation, the torus is greatly weakened, the central accretion timescale is shortened, and off-center fragmentation is suppressed --- triggering the central runaway even in previously `stable' models. The resulting SMBH masses lie in the range 2 x 10^4 Mo - 2 x 10^6 Mo, much higher than for Population III remnants. We argue that the above upper limit appears to be more realistic mass. Corollaries of this model include a possible correlation between SMBH and DM halo masses, and similarity between the SMBH and halo mass functions, at time of formation.

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
preprint2011arXiv

Galaxy Formation in Heavily Overdense Regions at z~10: the Prevalence of Disks in Massive Halos

Using a high-resolution cosmological numerical simulation, we have analyzed the evolution of galaxies at z~10 in a highly overdense region of the universe. These objects could represent the high redshift galaxies recently observed by the Hubble's WFC3, and be as well possible precursors of QSOs at z~6-7. To overcome the sampling and resolution problems in cosmological simulations, we have used the Constrained Realizations method. Our main result for z~10 shows the region of 3.5h^{-1}Mpc radius in comoving coordinates completely dominated by disk galaxies in the total mass range of >=10^9h^{-1}Mo. We have verified that the gaseous and stellar disks we identify are robust morphological features, capable of surviving the ongoing merger process at these redshifts. Below this mass range, we find a sharp decline in the disk fraction to negligible numbers. At this redshift, the disks appear to be gas-rich and the dark matter halos baryon-rich, by a factor of ~2-3 above the average fraction of baryons in the universe. The prevalence of disk galaxies in the high density peaks during the epoch of reionization is contrary to the morphology-density trend observed at low redshifts.

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

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

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

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

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

The dynamics of satellite disruption in cold dark matter haloes

We investigate the physical mechanisms of tidal heating and satellite disruption in cold dark matter host haloes using N-body simulations based on cosmological initial conditions. We show the importance of resonant shocks and resonant torques with the host halo to satellite heating. A resonant shock (torque) couples the radial (tangential) motion of a satellite in its orbit to its phase space. For a satellite on a circular orbit, an ILR-like resonance dominates the heating and this heating results in continuous satellite mass loss. We estimate the requirements for simulations to achieve these dynamics using perturbation theory. Both resonant shocks and resonant torques affect satellites on eccentric orbits. We demonstrate that satellite mass loss is an outside-in process in energy space; a satellite's stars and gas are thus protected by their own halo against tidal stripping. We simulate the evolution of a halo similar to the Large Magellanic Cloud (LMC) in our Galactic dark matter halo and conclude that the LMC stars have not yet been stripped. Finally, we present a simple algorithm for estimating the evolution of satellite mass that includes both shock heating and resonant torques.

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