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

Taysun Kimm

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

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-ph.IMastro-ph.HEastro-ph.SR

Catalog footprint

What is connected

27works
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

27 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.

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

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

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

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

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

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

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

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

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

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

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Formation of globular clusters in atomic-cooling halos via rapid gas condensation and fragmentation during the epoch of reionization

We investigate the formation of metal-poor globular clusters (GCs) at the center of two dark matter halos with $M_{\rm halo} \sim4\times10^7 M_\odot$ at $z>10$ using cosmological radiation-hydrodynamics simulations. We find that very compact ($\lesssim$ 1 pc) and massive ($\sim6\times10^5 M_\odot$) clusters form rapidly when pristine gas collapses isothermally with the aid of efficient Ly$α$ emission during the transition from molecular-cooling halos to atomic-cooling halos. Because the local free-fall time of dense star-forming gas is very short ($\ll 1\,{\rm Myr}$), a large fraction of the collapsed gas is turned into stars before stellar feedback processes blow out the gas and shut down star formation. Although the early stage of star formation is limited to a small region of the central star-forming disk, we find that the disk quickly fragments due to metal enrichment from supernovae. Sub-clusters formed in the fragmented clouds eventually merge with the main cluster at the center. The simulated clusters closely resemble the local GCs in mass and size but show a metallicity spread that is much wider than found in the local GCs. We discuss a role of pre-enrichment by Pop III and II stars as a potential solution to the latter issue. Although not without shortcomings, it is encouraging that a naive blind (not tuned) cosmological simulation presents a possible channel for the formation of at least some massive GCs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantifying Distributions of Lyman Continuum Escape Fraction

Simulations have indicated that most of the escaped Lyman continuum photons escape through a minority of solid angles with near complete transparency, with the remaining majority of the solid angles largely opaque, resulting in a very broad and skewed probability distribution function (PDF) of the escape fraction when viewed at different angles. Thus, the escape fraction of Lyman continuum photons of a galaxy observed along a line of sight merely represents the properties of the interstellar medium along that line of sight, which may be an ill-representation of true escape fraction of the galaxy averaged over its full sky. Here we study how Lyman continuum photons escape from galaxies at $z=4-6$, utilizing high-resolution large-scale cosmological radiation-hydrodynamic simulations. We compute the PDF of the mean escape fraction ($\left<f_{\rm esc,1D}\right>$) averaged over mock observational samples, as a function of the sample size, compared to the true mean (had you an infinite sample size). We find that, when the sample size is small, the apparent mean skews to the low end. For example, for a true mean of 6.7%, an observational sample of (2,10,50) galaxies at $z=4$ would have have 2.5% probability of obtaining the sample mean lower than $\left<f_{\rm esc,1D}\right>=$(0.007%, 1.8%, 4.1%) and 2.5% probability of obtaining the sample mean being greater than (43%, 18%, 11%). Our simulations suggest that at least $\sim$ 100 galaxies should be stacked in order to constrain the true escape fraction within 20% uncertainty.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Towards simulating star formation in turbulent high-z galaxies with mechanical supernova feedback

Feedback from supernovae is essential to understanding the self-regulation of star formation in galaxies. However, the efficacy of the process in a cosmological context remains unclear due to excessive radiative losses during the shock propagation. To better understand the impact of SN explosions on the evolution of galaxies, we perform a suite of high-resolution (12 pc), zoom-in cosmological simulations of a Milky Way-like galaxy at z=3 with adaptive mesh refinement. We find that SN explosions can efficiently regulate star formation, leading to the stellar mass and metallicity consistent with the observed mass-metallicity relation and stellar mass-halo mass relation at z~3. This is achieved by making three important changes to the classical feedback scheme: i) the different phases of SN blast waves are modelled directly by injecting radial momentum expected at each stage, ii) the realistic time delay of SNe, commencing at as early as 3 Myr, is required to disperse very dense gas before a runaway collapse sets in at the galaxy centre via mergers of gas clumps, and iii) a non-uniform density distribution of the ISM is taken into account below the computational grid scale for the cell in which SN explodes. The last condition is motivated by the fact that our simulations still do not resolve the detailed structure of a turbulent ISM in which the fast outflows can propagate along low-density channels. The simulated galaxy with the SN feedback model shows strong outflows, which carry approximately ten times larger mass than star formation rate, as well as smoothly rising circular velocity. Other feedback models that do not meet the three conditions form too many stars, producing a peaked rotation curve. Our results suggest that understanding the structure of the turbulent ISM may be crucial to assess the role of SN and other feedback processes in galaxy formation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dancing in the dark: galactic properties trace spin swings along the cosmic web

A large-scale hydrodynamical cosmological simulation, Horizon-AGN, is used to investigate the alignment between the spin of galaxies and the cosmic filaments above redshift 1.2. The analysis of more than 150 000 galaxies per time step in the redshift range 1.2<z<1.8 with morphological diversity shows that the spin of low-mass blue galaxies is preferentially aligned with their neighbouring filaments, while high-mass red galaxies tend to have a perpendicular spin. The reorientation of the spin of massive galaxies is provided by galaxy mergers, which are significant in their mass build-up. We find that the stellar mass transition from alignment to misalignment happens around 3.10^10 M_sun. Galaxies form in the vorticity-rich neighbourhood of filaments, and migrate towards the nodes of the cosmic web as they convert their orbital angular momentum into spin. The signature of this process can be traced to the properties of galaxies, as measured relative to the cosmic web. We argue that a strong source of feedback such as active galactic nuclei is mandatory to quench in situ star formation in massive galaxies and promote various morphologies. It allows mergers to play their key role by reducing post-merger gas inflows and, therefore, keeping spins misaligned with cosmic filaments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Diverse Properties of Interstellar Medium Embedding Gamma-Ray Bursts at the Epoch of Reionization

Analysis is performed on ultra-high resolution large-scale cosmological radiation-hydrodynamic simulations to, for the first time, quantify the physical environment of long-duration gamma-ray bursts (GRBs) at the epoch of reionization. We find that, on parsec scales, 13% of GRBs remain in high density ($\ge 10^4$cm$^{-3}$) low-temperature star-forming regions, whereas 87% of GRBs occur in low-density ($\sim 10^{-2.5}$cm$^{-3}$) high temperature regions heated by supernovae. More importantly, the spectral properties of GRB afterglows, such as the neutral hydrogen column density, total hydrogen column density, dust column density, gas temperature and metallicity of intervening absorbers, vary strongly from sightline to sightline. Although our model explains extant limited observationally inferred values with respect to circumburst density, metallicity, column density and dust properties, a substantially larger sample of high-z GRB afterglows would be required to facilitate a statistically solid test of the model. Our findings indicate that any attempt to infer the physical properties (such as metallicity) of the interstellar medium of the host galaxy based on a very small number of (usually one) sightlines would be precarious. Utilizing high-z GRBs to probe interstellar medium and intergalactic medium should be undertaken properly taking into consideration the physical diversities of the interstellar medium.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Escape fraction of ionizing photons during reionization: effects due to supernova feedback and runaway OB stars

The fraction of hydrogen ionizing photons escaping from galaxies into the intergalactic medium is a critical ingredient in the theory of reionization. We use two zoomed-in, high-resolution (4 pc), cosmological radiation hydrodynamic simulations with adaptive mesh refinement to investigate the impact of two physical mechanisms (supernova feedback and runaway OB stars) on the escape fraction (f_esc) at the epoch of reionization (z>7). We implement a new, physically motivated supernova feedback model that can approximate the Sedov solutions at all (from the free expansion to snowplow) stages. We find that there is a significant time delay of about ten million years between the peak of star formation and that of escape fraction, due to the time required for the build-up and subsequent destruction of the star-forming cloud by supernova feedback. Consequently, the photon number-weighted mean escape fraction for dwarf galaxies in halos of mass 10^8-10^10.5 Msun is found to be <fesc>~11%, although instantaneous values of f_esc>20% are common when star formation is strongly modulated by the supernova explosions. We find that the inclusion of runaway OB stars increases the mean escape fraction by 22% to <fesc>~14%. As supernovae resulting from runaway OB stars tend to occur in less dense environments, the feedback effect is enhanced and star formation is further suppressed in halos with Mvir>10^9 Msun in the simulation with runaway OB stars compared with the model without them. While both our models produce enough ionizing photons to maintain a fully ionized universe at z>7 as observed, a still higher amount of ionizing photons at z>9 appears necessary to accommodate the high observed electron optical depth inferred from cosmic microwave background observations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Heavy dust obscuration of z=7 galaxies in a cosmological hydrodynamic simulation

Hubble Space Telescope observations with Wide Field Camera 3/IR reveal that galaxies at z~7 have very blue ultraviolet (UV) colors, consistent with these systems being dominated by young stellar populations with moderate or little attenuation by dust. We investigate UV and optical properties of the high-z galaxies in the standard cold dark matter model using a high-resolution adaptive mesh refinement cosmological hydrodynamic simulation. For this purpose, we perform panchromatic three-dimensional dust radiative transfer calculations on 198 galaxies of stellar mass 5x10^8-3x10^{10} Msun with three parameters, the dust-to-metal ratio, the extinction curve, and the fraction of directly escaped light from stars (\fesc). Our stellar mass function is found to be in broad agreement with Gonzalez et al., independent of these parameters. We find that our heavily dust-attenuated galaxies (A_V~1.8) can also reasonably match modest UV-optical colors, blue UV slopes, as well as UV luminosity functions, provided that a significant fraction (~10%) of light directly escapes from them. The observed UV slope and scatter are better explained with an Small Magellanic Cloud-type extinction curve, whereas Milky Way-type dust predicts too blue UV colors due to the 2175A bump. We expect that upcoming observations by ALMA will be able to test this heavily obscured model.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Infrared Properties of z=7 Galaxies from Cosmological Simulations

Three-dimensional panchromatic dust radiative transfer calculations are performed on a set of 198 galaxies of stellar masses in the range 5x10^8-3x10^10 Msun from a cosmological hydrodynamic simulation (resolved at 29pc/h) at z=7. In a companion paper (Kimm & Cen), the stellar mass and UV luminosity functions, and UV-optical and FUV-NUV colors are shown to be in good agreement with observations, if an SMC-type dust extinction curve is adopted. Here we make useful predictions, self-consistently, of the infrared properties of these z=7 simulated galaxies that can be confronted with upcoming ALMA data. Our findings are as follows. (1) The effective radius in the rest-frame MIPS 70 micron band is in the range of 80-400pc proper for z=7 galaxies with L_FIR=10^{11.3-12}Lsun. (2) The median of the peak wavelength of the far-infrared (FIR) spectral energy distribution is in the range of 45-60 micron, depending on the dust-to-metal ratio. (3) For star formation rate in the range 3-100 Msun/yr the median FIR to bolometric luminosity ratio is 60-90%. (4) The FIR luminosity function displays a power law in the high end with a slope of -3.1 +- 0.4, instead of the usual exponential decline.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Blowing cold flows away: the impact of early AGN activity on the formation of a brightest cluster galaxy progenitor

Supermassive black holes (BH) are powerful sources of energy that are already in place at very early epochs of the Universe (by z=6). Using hydrodynamical simulations of the formation of a massive M_vir=5 10^11 M_sun halo by z=6 (the most massive progenitor of a cluster of M_vir=2 10^15 M_sun at z=0), we evaluate the impact of Active Galactic Nuclei (AGN) on galaxy mass content, BH self-regulation, and gas distribution inside this massive halo. We find that SN feedback has a marginal influence on the stellar structure, and no influence on the mass distribution on large scales. In contrast, AGN feedback alone is able to significantly alter the stellar-bulge mass content by quenching star formation when the BH is self-regulating, and by depleting the cold gas reservoir in the centre of the galaxy. The growth of the BH proceeds first by a rapid Eddington-limited period fed by direct cold filamentary infall. When the energy delivered by the AGN is sufficiently large to unbind the cold gas of the bulge, the accretion of gas onto the BH is maintained both by smooth gas inflow and clump migration through the galactic disc triggered by merger-induced torques. The feedback from the AGN has also a severe consequence on the baryon mass content within the halo, producing large-scale hot superwinds, able to blow away some of the cold filamentary material from the centre and reduce the baryon fraction by more than 30 per cent within the halo's virial radius. Thus in the very young universe, AGN feedback is likely to be a key process, shaping the properties of the most massive galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constraining stellar assembly and AGN feedback at the peak epoch of star formation

We study stellar assembly and feedback from active galactic nuclei (AGN) around the epoch of peak star formation (1<z<2), by comparing hydrodynamic simulations to rest-frame UV-optical galaxy colours from the Wide Field Camera 3 (WFC3) Early-Release Science (ERS) Programme. Our Adaptive Mesh Refinement simulations include metal-dependent radiative cooling, star formation, kinetic outflows due to supernova explosions, and feedback from supermassive black holes. Our model assumes that when gas accretes onto black holes, a fraction of the energy is used to form either thermal winds or sub-relativistic momentum-imparting collimated jets, depending on the accretion rate. We find that the predicted rest-frame UV-optical colours of galaxies in the model that includes AGN feedback is in broad agreement with the observed colours of the WFC3 ERS sample at 1<z<2. The predicted number of massive galaxies also matches well with observations in this redshift range. However, the massive galaxies are predicted to show higher levels of residual star formation activity than the observational estimates, suggesting the need for further suppression of star formation without significantly altering the stellar mass function. We discuss possible improvements, involving faster stellar assembly through enhanced star formation during galaxy mergers while star formation at the peak epoch is still modulated by the AGN feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Feeding compact bulges and supermassive black holes with low angular-momentum cosmic gas at high redshift

We use cosmological hydrodynamical simulations to show that a significant fraction of the gas in high redshift rare massive halos falls nearly radially to their very centre on extremely short timescales. This process results in the formation of very compact bulges with specific angular momentum a factor 5-30$smaller than the average angular momentum of the baryons in the whole halo. Such low angular momentum originates both from segregation and effective cancellation when the gas flows to the centre of the halo along well defined cold filamentary streams. These filaments penetrate deep inside the halo and connect to the bulge from multiple rapidly changing directions. Structures falling in along the filaments (satellite galaxies) or formed by gravitational instabilities triggered by the inflow (star clusters) further reduce the angular momentum of the gas in the bulge. Finally, the fraction of gas radially falling to the centre appears to increase with the mass of the halo; we argue that this is most likely due to an enhanced cancellation of angular momentum in rarer halos which are fed by more isotropically distributed cold streams. Such an increasingly efficient funnelling of low-angular momentum gas to the centre of very massive halos at high redshift may account for the rapid pace at which the most massive supermassive black holes grow to reach observed masses around $10^9$M$_\odot$ at an epoch when the Universe is barely 1 Gyr old.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Are cold flows detectable with metal absorption lines?

[Abridged] Cold gas flowing within the "cosmic web" is believed to be an important source of fuel for star formation at high redshift. However, the presence of such filamentary gas has never been observationally confirmed. In this work, we investigate in detail whether such cold gas is detectable using low-ionisation metal absorption lines, such as CII \lambda1334 as this technique has a proven observational record for detecting gaseous structures. Using a large statistical sample of galaxies from the Mare Nostrum N-body+AMR cosmological simulation, we find that the typical covering fraction of the dense, cold gas in 10^12 Msun haloes at z~2.5 is lower than expected (~5%). In addition, the absorption signal by the interstellar medium of the galaxy itself turns out to be so deep and so broad in velocity space that it completely drowns that of the filamentary gas. A detectable signal might be obtained from a cold filament exactly aligned with the line of sight, but this configuration is so unlikely that it would require surveying an overwhelmingly large number of candidate galaxies to tease it out. Finally, the predicted metallicity of the cold gas in filaments is extremely low (\leq 0.001 Zsun). Should this result persist when higher resolution runs are performed, it would significantly increase the difficulty of detecting filamentary gas inflows using metal lines. However, even if we assume that filaments are enriched to Zsun, the absorption signal that we compute is still weak. We are therefore led to conclude that it is extremely difficult to observationally prove or disprove the presence of cold filaments as the favorite accretion mode of galaxies using low-ionisation metal absorption lines. The Ly-alpha emission route looks more promising but due to the resonant nature of the line, radiative transfer simulations are required to fully characterize the observed signal.

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Impact of Gas Stripping and Stellar Mass Loss on Satellite Galaxy Evolution

Current semi-analytic models (SAMs) of galaxy formation over-predict the fraction of passive small late-type satellite galaxies in dense environments by a factor of two to three. We hypothesize that this is due to inaccurate prescriptions on cold gas evolution. In the hope of solving this problem we apply detailed prescriptions on the evolution of diffuse hot gases in satellites and on stellar mass loss, both of which are critical to model cold gas evolution. We replace the conventional shock-heating motivated instant stripping with a realistic gradual prescription based on ram pressure and tidal stripping. We also carefully consider stellar mass loss in our model. When both mechanisms are included, the fraction of passive late types matches the data much more closely. The satellite over-quenching problem is still present in small galaxies in massive haloes, however. In terms of the detectable residual star formation rates, gradual diffuse gas stripping appears to be much more important than stellar mass loss in our model. The implications of these results and other possibilities, such as redshift-dependent merging geometry and tidal disruption, are also discussed.

0 citations0 reviewsNo code linkedOpen access