Researcher profile

Taysun Kimm

Taysun Kimm contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

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

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

14 published item(s)

preprint2022arXiv

Impact of Radiation Feedback on the Formation of Globular Cluster Candidates during Cloud-Cloud Collisions

To understand the impact of radiation feedback during the formation of a globular cluster (GC), we simulate a head-on collision of two turbulent giant molecular clouds (GMCs). A series of idealized radiation-hydrodynamic simulations is performed, with and without stellar radiation or Type II supernovae. We find that a gravitationally bound, compact star cluster of mass $M_{\rm GC} \sim 10^5\,M_\odot$ forms within $\approx 3\,{\rm Myr}$ when two GMCs with mass $M_{\rm GMC}=3.6\times10^5\,M_\odot$ collide. The GC candidate does not form during a single collapsing event but emerges due to the mergers of local dense gas clumps and gas accretion. The momentum transfer due to the absorption of the ionizing radiation is the dominant feedback process that suppresses the gas collapse and photoionization becomes efficient once a sufficient number of stars form. The cluster mass is larger by a factor of $\sim 2$ when the radiation feedback is neglected, and the difference is slightly more pronounced (16\%) when extreme $\rm Lyα$ feedback is considered in the fiducial run. In the simulations with radiation feedback, supernovae explode after the star-forming clouds are dispersed, and their metal ejecta are not instantaneously recycled to form stars.

preprint2022arXiv

LyC escape from SPHINX galaxies in the Epoch of Reionization

We measure escape fractions, $f_{\rm esc}$, of ionizing radiation from galaxies in the SPHINX suite of cosmological radiation-hydrodynamical simulations of reionization, resolving halos with $M_{\rm vir} \gtrapprox 7.5 \times 10^7 \ M_{\odot}$ with a minimum cell width of $\approx 10$ pc. Our new and largest $20$ co-moving Mpc wide volume contains tens of thousands of star-forming galaxies with halo masses up to a few times $10^{11} \ M_{\odot}$. The simulated galaxies agree well with observational constraints of the UV luminosity function in the Epoch of Reionization. The escape fraction fluctuates strongly in individual galaxies over timescales of a few Myrs, due to its regulation by supernova and radiation feedback, and at any given time a tiny fraction of star-forming galaxies emits a large fraction of the ionizing radiation escaping into the inter-galactic medium. Statistically, $f_{\rm esc}$ peaks in intermediate-mass, intermediate-brightness, and low-metallicity galaxies ($M_{*} \approx 10^7 \ M_{\odot}$, $M_{1500} \approx -17$, $Z\lesssim 5 \times 10^{-3} \ Z_{\odot}$), dropping strongly for lower and higher masses, brighter and dimmer galaxies, and more metal-rich galaxies. The escape fraction correlates positively with both the short-term and long-term specific star formation rate. According to SPHINX, galaxies too dim to be yet observed, with $M_{1500} \gtrapprox -17$, provide about $55$ percent of the photons contributing to reionization. The global averaged $f_{\rm esc}$ naturally decreases with decreasing redshift, as predicted by UV background models and low-redshift observations. This evolution is driven by decreasing specific star formation rates over cosmic time.

preprint2022arXiv

MgII in the JWST Era: a Probe of Lyman Continuum Escape?

Limited constraints on the evolution of the Lyman Continuum (LyC) escape fraction represent one of the primary uncertainties in the theoretical determination of the reionization history. Due to the intervening intergalactic medium (IGM), the possibility of observing LyC photons directly in the epoch of reionization is highly unlikely. For this reason, multiple indirect probes of LyC escape have been identified, some of which are used to identify low-redshift LyC leakers (e.g. O32), while others are primarily useful at $z>6$ (e.g. [OIII]/[CII] far infrared emission). The flux ratio of the resonant MgII doublet emission at 2796$\dot{\rm A}$ and 2803$\dot{\rm A}$ as well as the MgII optical depth have recently been proposed as ideal diagnostics of LyC leakage that can be employed at $z>6$ with JWST. Using state-of-the-art cosmological radiation hydrodynamics simulations post-processed with CLOUDY and resonant-line radiative transfer, we test whether MgII is indeed a useful probe of LyC leakage. Our simulations indicate that the majority of bright, star-forming galaxies with high LyC escape fractions are expected to be MgII emitters rather than absorbers at $z=6$. However, we find that the MgII doublet flux ratio is a more sensitive indicator of dust rather than neutral hydrogen, limiting its use as a LyC leakage indicator to only galaxies in the optically thin regime. Given its resonant nature, we show that MgII will be an exciting probe of the complex kinematics in high-redshift galaxies in upcoming JWST observations.

preprint2022arXiv

Predicting LyC emission of galaxies using their physical and Ly$α$ emission properties

The primary difficulty in understanding the sources and processes that powered cosmic reionization is that it is not possible to directly probe the ionizing Lyman Continuum (LyC) radiation at that epoch as those photons have been absorbed by the intervening neutral hydrogen in the IGM on their way to us. It is therefore imperative to build a model to accurately predict LyC emission using other properties of galaxies in the reionization era. In recent years, studies have shown that the LyC emission from galaxies may be correlated to their Lya emission. Here, we study this correlation by analyzing thousands of galaxies at high-z in the SPHINX cosmological simulation. We post-process these galaxies with the Lya radiative transfer code RASCAS and analyze the Lya - LyC connection. We find that the Lya and LyC luminosities are strongly correlated with each other, although with dispersion. There is a positive correlation between Lya and LyC escape fractions in the brightest Lya emitters (>$10^{41}$ erg/s), similar to the recent observational studies. However, when we also include fainter Lya emitters (LAEs), the correlation disappears, which suggests that the observed relationship may be driven by selection effects. We also find that bright LAEs are dominant contributors to reionization ($> 10^{40}$ erg/s galaxies contribute $> 90\%$ of LyC emission). Finally, we build predictive models using multivariate linear regression where we use the physical and the Lya properties of simulated galaxies to predict their intrinsic and escaping LyC luminosities with a high degree of accuracy. We find that the most important galaxy properties to predict the escaping LyC luminosity of a galaxy are its escaping Lya luminosity, gas mass, gas metallicity, and SFR. These models can be very useful to predict LyC emissions from galaxies and can help us identify the sources of reionization.

preprint2022arXiv

Simulating Jellyfish Galaxies: A Case Study for a Gas-Rich Dwarf Galaxy

We investigate the formation of jellyfish galaxies using radiation-hydrodynamic simulations of gas-rich dwarf galaxies with a multi-phase interstellar medium (ISM). We find that the ram-pressure-stripped (RPS) ISM is the dominant source of molecular clumps in the near wake within 10 kpc from the galactic plane, while in-situ formation is the major channel for dense gas in the distant tail of the gas-rich galaxy. Only 20% of the molecular clumps in the near wake originate from the intracluster medium (ICM); however, the fraction reaches 50% in the clumps located at $80\,{\rm kpc}$ from the galactic center since the cooling time of the RPS gas tends to be short due to the ISM--ICM mixing ($\lesssim$ 10 Myr). The tail region exhibits a star formation rate of $0.001-0.01\,{\rm M_{\odot}\,yr^{-1}}$, and most of the tail stars are born in the stripped wake within 10 kpc from the galactic plane. These stars induce bright H$α$ blobs in the tail, while H$α$ tails fainter than $6\times10^{38}\,{\rm erg\,s^{-1}\,kpc^{-2}}$ are mostly formed via collisional radiation and heating due to mixing. We also find that the stripped tails have intermediate X-ray to H$α$ surface brightness ratios (1.5$\lesssim F_{\rm X}/F_{\rm Hα}\lesssim$20), compared to the ISM ($\lesssim$1.5) or pure ICM ($\gg$20). Our results suggest that jellyfish features emerge when the ISM from gas-rich galaxies is stripped by strong ram pressure, mixes with the ICM, and enhances the cooling in the tail.

preprint2022arXiv

The Nature of High ${\rm [OIII]}_{\rm 88μm}$/${\rm [CII]}_{\rm 158μm}$ Galaxies in the Epoch of Reionization: Low Carbon Abundance and a Top-Heavy IMF?

ALMA observations of $z>6$ galaxies have revealed abnormally high [OIII]$_{\rm 88μm}$/[CII]$_{\rm 158μm}$ ratios and [CII]$_{\rm 158μm}$ deficits compared to local galaxies. Numerous solutions have been proposed including differences in C and O abundance ratios, observational bias, and differences in ISM properties. In order to elucidate the underlying physics that drives this high-redshift phenomenon, we employ SPHINX$^{20}$, a state-of-the-art, cosmological radiation-hydrodynamics simulation, that resolves detailed ISM properties of thousands of galaxies in the epoch of reionization. We find that the observed $z>6$ [OIII]-SFR and [CII]-SFR relations can only be reproduced when the C/O abundance ratio is $\sim8\times$ lower than Solar and the total metal production is $\sim4\times$ higher than that of a Kroupa IMF. This implies that high-redshift galaxies are potentially primarily enriched by low-metallicity core-collapse supernovae with a more top-heavy IMF. As AGB stars and type-Ia supernova begin to contribute to the galaxy metallicity, both the [CII]-SFR and [CII] luminosity functions are predicted to converge to observed values at $z\sim4.5$. While we demonstrate that ionisation parameter, LyC escape fraction, ISM gas density, and CMB attenuation all drive galaxies towards higher [OIII]/[CII], observed values at $z>6$ can only be reproduced with substantially lower C/O abundances compared to Solar. The combination of [CII] and [OIII] can be used to predict the values of ionisation parameter, ISM gas density, and LyC escape fraction and we provide estimates of these quantities for nine observed $z>6$ galaxies. Finally, we demonstrate that [OI]$_{\rm 63μm}$ can be used as a replacement for [CII] in high-redshift galaxies where [CII] is unobserved and argue that more observation time should be used to target [OI] at $z>6$. (Abridged)

preprint2022arXiv

The Obelisk simulation: galaxies contribute more than AGN to HI reionization of protoclusters

We present the Obelisk project, a cosmological radiation-hydrodynamics simulation following the assembly and reionization of a protocluster progenitor during the first two billions of years from the big bang, down to $z = 3.5$. The simulation resolves haloes down to the atomic cooling limit, and tracks the contribution of different sources of ionization: stars, active galactic nuclei, and collisions. The Obelisk project is designed specifically to study the coevolution of high redshift galaxies and quasars in an environment favouring black hole growth. In this paper, we establish the relative contribution of these two sources of radiation to reionization and their respective role in establishing and maintaining the high redshift ionizing background. Our volume is typical of an overdense region of the Universe and displays star formation rate and black hole accretion rate densities similar to high redshift protoclusters. We find that hydrogen reionization happens inside-out and is completed by $z \sim 6$ in our overdensity, and is predominantly driven by galaxies, while accreting black holes only play a role at $z \sim 4$.

preprint2021arXiv

A systematic study of the escape of LyC and Ly$α$ photons from star-forming, magnetized turbulent clouds

Understanding the escape of Lyman continuum (LyC) and Lyman $α$ (Ly$α$) photons from giant molecular clouds (GMCs) is crucial if we are to study the reionization of the Universe and to interpret spectra of observed galaxies at high redshift. To this end, we perform high-resolution, radiation-magneto-hydrodynamic simulations of GMCs with self-consistent star formation and stellar feedback. We find that a significant fraction (15-70%) of ionizing radiation escapes from the simulated GMCs with different masses ($10^5$ and $10^6\,M_\odot$), as the clouds are dispersed within about $2$-$5\,{\rm Myr}$ from the onset of star formation. The fraction of LyC photons leaked is larger when the GMCs are less massive, metal-poor, less turbulent, and less dense. The most efficient leakage of LyC radiation occurs when the total star formation efficiency of a GMC is about 20%. The escape of Ly$α$ shows a trend similar to that of LyC photons, except that the fraction of Ly$α$ photons escaping from the GMCs is larger ($f_{\rm esc}^{\rm Lyα}\approx f_{900}^{0.27}$) and that a GMC with strong turbulence shows larger $f_{\rm esc}^{\rm Lyα}$. The simulated GMCs show a characteristic velocity separation of $Δv\approx 120 \,{\rm km\,s^{-1}}$ in the time-averaged emergent Ly$α$ spectra, suggesting that Ly$α$ could be useful to infer the kinematics of the interstellar and circumgalactic medium. We show that Ly$α$ luminosities are a useful indicator of the LyC escape, provided the number of LyC photons can be deduced through stellar population modeling. Finally, we find that the correlations between the escape fractions of Ly$α$, ultraviolet photons at 1500A, and the Balmer $α$ line are weak.

preprint2021arXiv

Predicting the Observability of Population III Stars with ELT-HARMONI via the Helium $1640{\rmÅ}$ emission line

Population III (Pop. III) stars, as of yet, have not been detected, however as we move into the era of extremely large telescopes this is likely to change. One likely tracer for Pop. III stars is the HeII$\lambda1640$ emission line, which will be detectable by the HARMONI spectrograph on the European Extremely Large Telescope (ELT) over a broad range of redshifts ($2\leq z\leq14$). By post-processing galaxies from the cosmological, AMR-hydrodynamical simulation NewHorizon with theoretical spectral energy distributions (SED) for Pop. III stars and radiative transfer (i.e. the Yggdrasil Models and CLOUDY look-up tables respectively) we are able to compute the flux of HeII$\lambda1640$ for individual galaxies. From mock 10 hour observations of these galaxies we show that HARMONI will be able to detect Pop. III stars in galaxies up to $z\sim10$ provided Pop. III stars have a top heavy Initial Mass Function (IMF). Furthermore, we find that should Pop. III stars instead have an IMF similar to those of the Pop. I stars, the HeII$\lambda1640$ line would only be observable for galaxies with Pop. III stellar masses in excess of $10^{7}{\,M}_\odot\,{\rm yr}^{-1}$, average stellar age $<1{\rm Myr}$ at $z=4$. Finally, we are able to determine the minimal intrinsic flux required for HARMONI to detect Pop. III stars in a galaxy up to $z=10$.

preprint2020arXiv

How to Quench a Dwarf Galaxy: The Impact of Inhomogeneous Reionization on Dwarf Galaxies and Cosmic Filaments

We use the SPHINX suite of high-resolution cosmological radiation hydrodynamics simulations to study how spatially and temporally inhomogeneous reionization impacts the baryonic content of dwarf galaxies and cosmic filaments. The SPHINX simulations simultaneously model an inhomogeneous reionization, follow the escape of ionising radiation from thousands of galaxies, and resolve haloes well below the atomic cooling threshold. This makes them an ideal tool for examining how reionization impacts star formation and the gas content of dwarf galaxies. We compare simulations with and without stellar radiation to isolate the effects of radiation feedback from that of supernova, cosmic expansion, and numerical resolution. We find that the gas content of cosmic filaments can be reduced by more than 80% following reionization. The gas inflow rates into haloes with $M_{vir}<10^8M_{\odot}$ are strongly affected and are reduced by more than an order of magnitude compared to the simulation without reionization. A significant increase in gas outflow rates is found for halo masses $M_{vir}<7\times10^7M_{\odot}$. Our simulations show that inflow suppression, rather than photoevaporation, is the dominant mechanism by which the baryonic content of high-redshift dwarf galaxies is regulated. At fixed redshift and halo mass, there is a large scatter in the halo baryon fractions that is entirely dictated by the timing of reionization in the local region surrounding a halo which can change by $Δz>3$ at fixed mass. Finally, although the gas content of high-redshift dwarf galaxies is significantly impacted by reionization, we find that most haloes with $M_{vir}<10^8M_{\odot}$ can remain self-shielded and form stars long after reionization, until their local gas reservoir is depleted, suggesting that local group dwarf galaxies do not necessarily exhibit star formation histories that peak prior to $z=6$...

preprint2020arXiv

New Methods for Identifying Lyman Continuum Leakers and Reionization-Epoch Analogues

Identifying low-redshift galaxies that emit Lyman Continuum radiation (LyC leakers) is one of the primary, indirect methods of studying galaxy formation in the epoch of reionization. However, not only has it proved challenging to identify such systems, it also remains uncertain whether the low-redshift LyC leakers are truly &#34;analogues&#34; of the sources that reionized the Universe. Here, we use high-resolution cosmological radiation hydrodynamics simulations to examine whether simulated galaxies in the epoch of reionization share similar emission line properties to observed LyC leakers at $z\sim3$ and $z\sim0$. We find that the simulated galaxies with high LyC escape fractions ($f_{\rm esc}$) often exhibit high O32 and populate the same regions of the R23-O32 plane as $z\sim3$ LyC leakers. However, we show that viewing angle, metallicity, and ionisation parameter can all impact where a galaxy resides on the O32-$f_{\rm esc}$ plane. Based on emission line diagnostics and how they correlate with $f_{\rm esc}$, lower-metallicity LyC leakers at $z\sim3$ appear to be good analogues of reionization-era galaxies. In contrast, local [SII]-deficient galaxies do not overlap with the simulated high-redshift LyC leakers on the SII-BPT diagram; however, this diagnostic may still be useful for identifying leakers. We use our simulated galaxies to develop multiple new diagnostics to identify LyC leakers using IR and nebular emission lines. We show that our model using only [CII]$_{\rm 158μm}$ and [OIII]$_{\rm 88μm}$ can identify potential leakers from non-leakers from the local Dwarf Galaxy Survey. Finally, we apply this diagnostic to known high-redshift galaxies and find that MACS1149_JD1 at $z=9.1$ is the most likely galaxy to be actively contributing to the reionization of the Universe.

preprint2020arXiv

RASCAS: RAdiation SCattering in Astrophysical Simulations

Resonant lines are powerful probes of the interstellar and circumgalactic medium of galaxies. Their transfer in gas being a complex process, the interpretation of their observational signatures, either in absorption or in emission, is often not straightforward. Numerical radiative transfer simulations are needed to accurately describe the travel of resonant line photons in real and in frequency space, and to produce realistic mock observations. This paper introduces RASCAS, a new public 3D radiative transfer code developed to perform the propagation of any resonant line in numerical simulations of astrophysical objects. RASCAS was designed to be easily customisable and to process simulations of arbitrarily large sizes on large supercomputers. RASCAS performs radiative transfer on an adaptive mesh with an octree structure using the Monte Carlo technique. RASCAS features full MPI parallelisation, domain decomposition, adaptive load-balancing, and a standard peeling algorithm to construct mock observations. The radiative transport of resonant line photons through different mixes of species (e.g. \ion{H}{i}, \ion{Si}{ii}, \ion{Mg}{ii}, \ion{Fe}{ii}), including their interaction with dust, is implemented in a modular fashion to allow new transitions to be easily added to the code. RASCAS is very accurate and efficient. It shows perfect scaling up to a minimum of a thousand cores. It has been fully tested against radiative transfer problems with analytic solutions and against various test cases proposed in the literature. Although it was designed to describe accurately the many scatterings of line photons, RASCAS may also be used to propagate photons at any wavelength (e.g. stellar continuum or fluorescent lines), or to cast millions of rays to integrate the optical depths of ionising photons, making it highly versatile.

preprint2020arXiv

Reconstruction of Radio Relics and X-ray Tails in an Off-axis Cluster Merger: Hydrodynamical Simulations of A115

Although a merging galaxy cluster is a useful laboratory to study many important astrophysical processes and fundamental physics, only limited interpretations are possible without careful analysis of the merger scenario. However, the study is demanding because a thorough comparison of multi-wavelength observations with detailed numerical simulations is required. In this paper, we present such a study for the off-axis binary merger A115. The system possesses a number of remarkable observational features, but no convincing merger scenario, explaining the shape and location of the radio relic in harmony with the orientation of the cometary X-ray tails, has been presented. Our hydrodynamical simulation, with adaptive mesh refinement, suggests that the cometary X-ray tail of A115 might be a slingshot tail and can arise $\sim0.3$ Gyrs after the impact and before the two subclusters reach their apocenters. This scenario can predict the location and orientation of the giant radio relic, which is parallel to the northern X-ray tail. In addition, our study indicates that diffusive shock acceleration alone cannot generate the observed radio power unless aided by stronger magnetic fields and/or more significant presence of fossil electrons.

preprint2011arXiv

The angular momentum of baryons and dark matter halos revisited

Recent theoretical studies have shown that galaxies at high redshift are fed by cold, dense gas filaments, suggesting angular momentum transport by gas differs from that by dark matter. Revisiting this issue using high-resolution cosmological hydrodynamics simulations with adaptive mesh refinement, we find that at the time of accretion, gas and dark matter do carry a similar amount of specific angular momentum, but that it is systematically higher than that of the dark matter halo as a whole. At high redshift, freshly accreted gas rapidly streams into the central region of the halo, directly depositing this large amount of angular momentum within a sphere of radius r=0.1rvir. In contrast, dark matter particles pass through the central region unscathed, and a fraction of them ends up populating the outer regions of the halo (r/rvir>0.1), redistributing angular momentum in the process. As a result, large-scale motions of the cosmic web have to be considered as the origin of gas angular momentum rather than its virialised dark matter halo host. This generic result holds for halos of all masses at all redshifts, as radiative cooling ensures that a significant fraction of baryons remain trapped at the centre of the halos. Despite this injection of angular momentum enriched gas, we predict an amount for stellar discs which is in fair agreement with observations at z=0. This arises because the total specific angular momentum of the baryons remains close to that of dark matter halos. We propose a new scenario where gas efficiently carries the angular momentum generated by large-scale structure motions deep inside dark matter halos, redistributing it only in the vicinity of the disc.