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

Hyunbae Park

Hyunbae Park 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.IM

Catalog footprint

What is connected

14works
3topics
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

14 published item(s)

preprint2022arXiv

Scattering of Lyα Photons through the Reionizing Intergalactic Medium: I. Spectral Energy Distribution

During reionization, a fraction of galactic Ly$α$ emission is scattered in the intergalactic medium (IGM) and appears as a diffuse light extending megaparsecs from the source. We investigate how to probe the properties of the early galaxies and their surrounding IGM using this scattered light. We create a Monte Carlo algorithm to track individual photons and reproduce several test cases from previous literature. Then, we run our code on the simulated IGM of the CoDaII simulation. We find that the scattered light can leave an observable imprint on the emergent spectrum if collected over several square arcminutes. Scattering can redden the emission by increasing the path lengths of photons, but it can also make the photons bluer by upscattering them according to the peculiar motion of the scatterer. The photons emitted on the far blue side of the resonance appear more extended in both frequency and space compared to those emitted near the resonance. This provides a discriminating feature for the blueward emission, which cannot be constrained from the unscattered light coming directly from the source. The ionization state of the IGM also affects the scattered light spectrum. When the source is in a small HII region, the emission goes through more scatterings in the surrounding HI region regardless of the initial frequency and ends up more redshifted and spatially extended. This can result in a weakening of the scattered light toward high $z$ during reionization. Our results provide a framework for interpreting the scattered light to be measured by high-$z$ integral-field-unit surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Snowmass2021 Computational Frontier White Paper: Cosmological Simulations and Modeling

Powerful new observational facilities will come online over the next decade, enabling a number of discovery opportunities in the "Cosmic Frontier", which targets understanding of the physics of the early universe, dark matter and dark energy, and cosmological probes of fundamental physics, such as neutrino masses and modifications of Einstein gravity. Synergies between different experiments will be leveraged to present new classes of cosmic probes as well as to minimize systematic biases present in individual surveys. Success of this observational program requires actively pairing it with a well-matched state-of-the-art simulation and modeling effort. Next-generation cosmological modeling will increasingly focus on physically rich simulations able to model outputs of sky surveys spanning multiple wavebands. These simulations will have unprecedented resolution, volume coverage, and must deliver guaranteed high-fidelity results for individual surveys as well as for the cross-correlations across different surveys. The needed advances are as follows: (1) Development of scientifically rich and broadly-scoped simulations, which capture the relevant physics and correlations between probes (2) Accurate translation of simulation results into realistic image or spectral data to be directly compared with observations (3) Improved emulators and/or data-driven methods serving as surrogates for expensive simulations, constructed from a finite set of full-physics simulations (4) Detailed and transparent verification and validation programs for both simulations and analysis tools. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The short ionizing photon mean free path at z=6 in Cosmic Dawn III, a new fully-coupled radiation-hydrodynamical simulation of the Epoch of Reionization

Recent determinations of the mean free path of ionising photons (mfp) in the intergalactic medium (IGM) at $\rm z=6$ are lower than many theoretical predictions. To gain insight into this issue, we investigate the evolution of the mfp in our new massive fully coupled radiation hydrodynamics cosmological simulation of reionization: Cosmic Dawn III (CoDaIII). CoDaIII's scale ($\rm 94^3 \, cMpc^3$) and resolution ($\rm 8192^3$ grid) make it particularly suitable to study the evolution of the IGM during Reionization. The simulation was performed with RAMSES-CUDATON on Summit, and used 131072 processors coupled to 24576 GPUs, making it the largest Reionization simulation, and largest RAMSES simulation ever performed. A superior agreement with global constraints on Reionization is obtained in CoDaIII over CoDaII especially for the evolution of the neutral hydrogen fraction and the cosmic photo-ionization rate, thanks to an improved calibration, later end of reionization ($\rm z=5.6$), and higher spatial resolution. Analyzing the mfp, we find that CoDaIII reproduces the most recent observations very well, from $\rm z=6$ to $\rm z=4.6$. We show that the distribution of the mfp in CoDaIII is bimodal, with short (neutral) and long (ionized) mfp modes, respectively, due to the patchiness of reionization and the co-existence of neutral versus ionized regions during Reionization. The neutral mode peaks at sub-kpc to kpc scales of mfp, while the ionized mode peak evolves from $\rm 0.1 Mpc/h$ at $\rm z=7$ to $\sim 10$ Mpc/h at $\rm z=5.2$. Computing the mfp as the average of the ionized mode provides the best match to the recent observational determinations. The distribution reduces to a single neutral (ionized) mode at $\rm z>13$ ($\rm z<5$).

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Cosmological constraints from the redshift dependence of the Alcock-Paczynski effect: Possibility of estimateing the non-linear systematics using fast simulations

The tomographic AP method is so far the best method in separating the Alcock-Paczynski (AP) signal from the redshift space distortion (RSD) effects and deriving powerful constraints on cosmological parameters using the $\lesssim40h^{-1}\ \rm Mpc$ clustering region. To guarantee that the method can be easily applied to the future large scale structure (LSS) surveys, we study the possibility of estimating the systematics of the method using fast simulation method. The major contribution of the systematics comes from the non-zero redshift evolution of the RSD effects, which is quantified by $\hatξ_{Δs}(μ,z)$ in our analysis, and estimated using the BigMultidark exact N-body simulation and approximate COLA simulation samples. We find about 5\%/10\% evolution when comparing the $\hatξ_{Δs}(μ,z)$ measured as $z=0.5$/$z=1$ to the measurements at $z=0$. We checked the inaccuracy in the 2pCFs computed using COLA, and find it 5-10 times smaller than the intrinsic systematics of the tomographic AP method, indicating that using COLA to estimate the systematics is good enough. Finally, we test the effect of halo bias, and find $\lesssim$1.5\% change in $\hatξ_{Δs}$ when varying the halo mass within the range of $2\times 10^{12}$ to $10^{14}$ $M_{\odot}$. We will perform more studies to achieve an accurate and efficient estimation of the systematics in redshift range of $z=0-1.5$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Cosmological Information from the Small-scale Redshift Space Distortions

The redshift-space distortion (RSD) in the observed distribution of galaxies is known as a powerful probe of cosmology. Observations of large-scale RSD have given tight constraints on the linear growth rate of the large-scale structures in the universe. On the other hand, the small-scale RSD, caused by galaxy random motions inside clusters, has not been much used in cosmology, but also has cosmological information because universes with different cosmological parameters have different halo mass functions and virialized velocities. We focus on the projected correlation function $w(r_p)$ and the multipole moments $ξ_l$ on small scales ($1.4$ to $30\ h^{-1}\rm{Mpc}$). Using simulated galaxy samples generated from a physically motivated most bound particle (MBP)-galaxy correspondence scheme in the Multiverse Simulation, we examine the dependence of the small-scale RSD on the cosmological matter density parameter $Ω_m$, the satellite velocity bias with respect to MBPs, $b_v^s$, and the merger-time-scale parameter $α$. We find that $α=1.5$ gives an excellent fit to the $w(r_p)$ and $ξ_l$ measured from the SDSS-KIAS value added galaxy catalog. We also define the ``strength'' of Fingers-of-God as the ratio of the parallel and perpendicular size of the contour in the two-point correlation function set by a specific threshold value and show that the strength parameter helps constraining $(Ω_m, b_v^s, α)$ by breaking the degeneracy among them. The resulting parameter values from all measurements are $(Ω_m,b_v^s)=(0.272\pm0.013,0.982\pm0.040)$, indicating a slight reduction of satellite galaxy velocity relative to the MBP. However, considering that the average MBP speed inside haloes is $0.94$ times the dark matter velocity dispersion, the main drivers behind the galaxy velocity bias are gravitational interactions, rather than baryonic effects.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The impact of inhomogeneous subgrid clumping on cosmic reionization

Cosmic reionization was driven by the imbalance between early sources and sinks of ionizing radiation, both of which were dominated by small-scale structure and are thus usually treated in cosmological reionization simulations by subgrid modelling. The recombination rate of intergalactic hydrogen is customarily boosted by a subgrid clumping factor, ${\left<n^2\right>/\left<n\right>^2}$, which corrects for unresolved fluctuations in gas density ${n}$ on scales below the grid-spacing of coarse-grained simulations. We investigate in detail the impact of this inhomogeneous subgrid clumping on reionization and its observables, as follows: (1) Previous attempts generally underestimated the clumping factor because of insufficient mass resolution. We perform a high-resolution $N$-body simulation that resolves haloes down to the pre-reionization Jeans mass to derive the time-dependent, spatially-varying local clumping factor and a fitting formula for its correlation with local overdensity. (2) We then perform a large-scale $N$-body and radiative transfer simulation that accounts for this inhomogeneous subgrid clumping by applying this clumping factor-overdensity correlation. Boosting recombination significantly slows the expansion of ionized regions, which delays completion of reionization and suppresses 21 cm power spectra on large scales in the later stages of reionization. (3) We also consider a simplified prescription in which the globally-averaged, time-evolving clumping factor from the same high-resolution $N$-body simulation is applied uniformly to all cells in the reionization simulation, instead. Observables computed with this model agree fairly well with those from the inhomogeneous clumping model, e.g. predicting 21 cm power spectra to within 20% error, suggesting it may be a useful approximation.

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 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
preprint2015arXiv

The Impact of Nonlinear Structure Formation on the Power Spectrum of Transverse Momentum Fluctuations and the Kinetic Sunyaev-Zel'dovich Effect

Cosmological transverse momentum fields, whose directions are perpendicular to Fourier wave vectors, induce temperature anisotropies in the cosmic microwave background via the kinetic Sunyaev-Zeldovich (kSZ) effect. The transverse momentum power spectrum contains the four-point function of density and velocity fields, $\langleδδv v\rangle$. In the post-reionization epoch, nonlinear effects dominate in the power spectrum. We use perturbation theory and cosmological $N$-body simulations to calculate this nonlinearity. We derive the next-to-leading order expression for the power spectrum with a particular emphasis on the connected term that has been ignored in the literature. While the contribution from the connected term on small scales ($k>0.1\,h\,\rm{Mpc}^{-1}$) is subdominant relative to the unconnected term, we find that its contribution to the kSZ power spectrum at $\ell = 3000$ at $z<6$ can be as large as ten percent of the unconnected term, which would reduce the allowed contribution from the reionization epoch ($z>6$) by twenty percent. The power spectrum of transverse momentum on large scales is expected to scale as $k^2$ as a consequence of momentum conservation. We show that both the leading and the next-to-leading order terms satisfy this scaling. In particular, we find that both of the unconnected and connected terms are necessary to reproduce $k^2$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Kinetic Sunyaev-Zel'dovich effect as a probe of the physics of cosmic reionization: the effect of self-regulated reionization

We calculate the angular power spectrum of the Cosmic Microwave Background (CMB) temperature fluctuations induced by the kinetic Sunyaev-Zel'dovich (kSZ) effect from the epoch of reionization (EOR). We use detailed N-body+radiative transfer simulations to follow inhomogeneous reionization of the intergalactic medium (IGM). For the first time we take into account the "self-regulation" of reionization: star formation in low-mass dwarf galaxies (10^8 M_\sun \lesssim M \lesssim 10^9 M_\sun) or minihalos (10^5 M_\sun \lesssim M \lesssim 10^8 M_\sun) is suppressed if these halos form in the regions that were already ionized or Lyman-Werner dissociated. Some previous work suggested that the amplitude of the kSZ power spectrum from the EOR can be described by a two-parameter family: the epoch of half ionization and the duration of reionization. However, we argue that this picture applies only to simple forms of the reionization history which are roughly symmetric about the half-ionization epoch. In self-regulated reionization, the universe begins to be ionized early, maintains a low level of ionization for an extended period, and then finishes reionization as soon as high-mass atomically-cooling halos dominate. While inclusion of self-regulation affects the amplitude of the kSZ power spectrum only modestly (\sim 10 %), it can change the duration of reionization by a factor of more than two. We conclude that the simple two-parameter family does not capture the effect of a physical, yet complex, reionization history caused by self-regulation. When added to the post-reionization kSZ contribution, our prediction for the total kSZ power spectrum is below the current upper bound from the South Pole Telescope. Therefore, the current upper bound on the kSZ effect from the EOR is consistent with our understanding of the physics of reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simulating Cosmic Reionization and the Radiation Backgrounds from the Epoch of Reionization

Large-scale reionization simulations are described which combine the results of cosmological N-body simulations that model the evolving density and velocity fields and identify the galactic halo sources, with ray-tracing radiative transfer calculations which model the nonequilibrium ionization of the intergalactic medium. These simulations have been used to predict some of the signature effects of reionization on cosmic radiation backgrounds, including the CMB, near-IR, and redshifted 21cm backgrounds. We summarize some of our recent progress in this work, and address the question of whether observations of such signature effects can be used to distinguish the relative contributions of galaxies of different masses to reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gravitational Potential Environment of Galaxies: I. Simulation

We extend the concept of galaxy environment from the local galaxy number density to the gravitational potential and its functions like the shear tensor. For this purpose we examine whether or not one can make an accurate estimation of the gravitational potential from an observational sample which is finite in volume, biased due to galaxy biasing, and subject to redshift space distortion. Dark halos in a $Λ$CDM simulation are used in this test. We find that one needs to stay away from the sample boundaries by more than 30$h^{-1}$Mpc to reduce the error within 20% of the root mean square values of the potential or the shear tensor. The error due to the galaxy biasing can be significantly reduced by using the galaxy mass density field instead of the galaxy number density field. The error caused by the redshift space distortion can be effectively removed by correcting galaxy positions for the peculiar velocity effects. We inspect the dependence of dark matter halo properties on four environmental parameters; local density, gravitational potential, and the ellipticity and prolateness of the shear tensor. We find the local density has the strongest correlation with halo properties. This is evidence that the internal physical properties of dark halos are mainly controlled by small-scale physics. In high density regions dark halos are on average more massive and spherical, and have higher spin parameter and velocity dispersion. In high density regions dark halos are on average more massive and spherical, and have higher spin parameter and velocity dispersion. We also study the relation between the environmental parameters and the subtypes of dark halos. The spin parameter of satellite halos depends only weakly on the local density for all mass ranges studied while that of isloated or central halos depends more sensitively on the local density. (abridged)

0 citations0 reviewsNo code linkedOpen access