BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Harley Katz

Harley Katz contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.GAastro-ph.COastro-ph.SR

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
3topics
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

12 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

RAMSES-RTZ: Non-Equilibrium Metal Chemistry and Cooling Coupled to On-The-Fly Radiation Hydrodynamics

Emission and absorption lines from elements heavier than helium (metals) represent one of our strongest probes of galaxy formation physics across nearly all redshifts accessible to observations. The vast majority of simulations that model these metal lines often assume either collisional or photoionisation equilibrium, or a combination of the two. For the few simulations that have relaxed these assumptions, a redshift-dependent meta-galactic UV background or fixed spectrum is often used in the non-equilibrium photoionisation calculation, which is unlikely to be accurate in the interstellar medium where the gas can self-shield as well as in the high-redshift circumgalactic medium where locally emitted radiation may dominate over the UV background. In this work, we relax this final assumption by coupling the ionisation states of individual metals to the radiation hydrodynamics solver present in RAMSES-RT. Our chemical network follows radiative recombination, dielectronic recombination, collisional ionisation, photoionisation, and charge transfer and we use the ionisation states to compute non-equilibrium optically-thin metal-line cooling. The fiducial model solves for the ionisation states of C, N, O, Mg, Si, S, Fe, and Ne in addition to H, He, and H$_2$, but can be easily extended for other ions. We provide interfaces to two different ODE solvers that are competitive in both speed and accuracy. The code has been benchmarked across a variety of gas conditions to reproduce results from CLOUDY when equilibrium is reached. We show an example isolated galaxy simulation with on-the-fly radiative transfer that demonstrates the utility of our code for translating between simulations and observations without the use of idealised photoionisation models.

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
preprint2022arXiv

Towards convergence of turbulent dynamo amplification in cosmological simulations of galaxies

Our understanding of the process through which magnetic fields reached their observed strengths in present-day galaxies remains incomplete. One of the advocated solutions is a turbulent dynamo mechanism that rapidly amplifies weak magnetic field seeds to the order of ${\sim}μ$G. However, simulating the turbulent dynamo is a very challenging computational task due to the demanding span of spatial scales and the complexity of the required numerical methods. In particular, turbulent velocity and magnetic fields are extremely sensitive to the spatial discretisation of simulated domains. To explore how refinement schemes affect galactic turbulence and amplification of magnetic fields in cosmological simulations, we compare two refinement strategies. A traditional quasi-Lagrangian adaptive mesh refinement approach focusing spatial resolution on dense regions, and a new refinement method that resolves the entire galaxy with a high resolution quasi-uniform grid. Our new refinement strategy yields much faster magnetic energy amplification than the quasi-Lagrangian method, which is also significantly greater than the adiabatic compressional estimate indicating that the extra amplification is produced through stretching of magnetic field lines. Furthermore, with our new refinement the magnetic energy growth factor scales with resolution following $\propto \Dres^{-1/2}$, in much better agreement with small-scale turbulent box simulations. Finally, we find evidence suggesting most magnetic amplification in our simulated galaxies occurs in the warm phase of their interstellar medium, which has a better developed turbulent field with our new refinement strategy.

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 &#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.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Black Hole Formation in the First Stellar Clusters

The early Universe was composed almost entirely of hydrogen and helium, with only trace amounts of heavy elements. It was only after the first generation of star formation that the Universe became sufficiently polluted to produce a second generation (Population II) of stars which are similar to those in our local Universe. Evidence of massive star cluster formation is nearly ubiquitous among the observed galaxy population and if this mode of star formation occurred at early enough epochs, the higher densities in the early Universe may have caused many of the stars in the cluster to strongly interact. In this scenario, it may be possible to form a very massive star by repeated stellar collisions that may directly collapse into a black hole and form a supermassive black hole seed. In this chapter, we will explore this scenario in detail to understand the dynamics which allow for this process to ensue and measure the probability for this type of seed to represent the supermassive black hole population observed at z > 6.

0 citations0 reviewsNo code linkedOpen access