Dr Harley Katz

Abstract:

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.Comment: Accepted to Astronomy and Astrophysics (A&A) Journal. 27 pages, 21 Figure

Abstract:

We investigate the formation of jellyfish galaxies using radiation-hydrodynamic simulations of gas-rich dwarf galaxies with a multiphase 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 kpc from the galactic center since the cooling time of the RPS gas tends to be short owing to the ISM–ICM mixing (≲10 Myr). The tail region exhibits a star formation rate of 0.001–0.01 M⊙ 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 × 1038 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 ≲ FX/FHα ≲ 20), compared to the ISM (≲1.5) or pure ICM (≫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.

Abstract:

ALMA observations of z > 6 galaxies hav e rev ealed abnormally high [O III ] 88 μm /[C II ] 158 μm ratios and [C II ] 158 μm deficits compared to local galaxies. The origin of this behaviour is unknown. Numerous solutions have been proposed including differences in C and O abundance ratios, observational bias, and differences in ISM properties, including ionization parameter, gas density, or photodissociation region (PDR) covering fraction. 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 which has been post-processed with CLOUDY to predict emission lines. We find that the observed z > 6 [O III ] 88 μm -SFR and [C II ] 158 μm -SFR relations can only be reproduced when the C/O abundance ratio is ∼8 ×lower than Solar and the total metal production is ∼4 ×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 [C II ] 158 μm -SFR and [C II ] 158 μm luminosity functions are predicted to converge to observed values at z ∼4.5. While we demonstrate that ionization parameter, LyC escape fraction, ISM gas density, and CMB attenuation all drive galaxies towards higher [O III ] 88 μm /[C II ] 158 μm , observed values at z > 6 can only be reproduced with substantially lower C/O abundances compared to Solar. The combination of [C II ] 158 μm and [O III ] 88 μm can be used to predict the values of ionization 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 [O I ] 63 μm can be used as a replacement for [C II ] 158 μm in high-redshift galaxies where [C II ] 158 μm is unobserved and argue that more observation time should be used to target [O I ] 63 μm at z > 6. Future simulations will be needed to self-consistently address the numerous uncertainties surrounding a varying IMF at high redshift and the associated metal returns.

Abstract:

We present the first results from SPHINX-MHD, a suite of cosmological radiation-magnetohydrodynamics simulations designed to study the impact of primordial magnetic fields (PMFs) on galaxy formation and the evolution of the intergalactic medium during the epoch of reionization. The simulations are among the first to employ on-the-fly radiation transfer and constrained transport ideal MHD in a cosmological context to simultaneously model the inhomogeneous process of reionization and the growth of PMFs. We run a series of $(5{\rm Mpc})^3$ cosmological volumes, varying both the strength of the seed magnetic field and its spectral index. We find that PMFs with a spectral index ($n_B$) and a comoving amplitude ($B_0$) that have $n_B>-0.562\log_{10}(B_0/1{\rm n}G) - 3.35$ produce electron optical depths ($\tau_e$) that are inconsistent with CMB constraints due to the unrealistically early collapse of low-mass dwarf galaxies. For $n_B\geq-2.9$, our constraints are considerably tighter than the $\sim{\rm n}G$ constraints from Planck. PMFs that do not satisfy our constraints have little impact on the reionization history or the shape of the UV luminosity function. Likewise, detecting changes in the Ly$\alpha$ forest due to PMFs will be challenging because photoionisation and photoheating efficiently smooth the density field. However, we find that the first absorption feature in the global 21cm signal is a particularly sensitive indicator of the properties of the PMFs, even for those that satisfy our $\tau_e$ constraint. Furthermore, strong PMFs can increase the escape of LyC photons by up to 25% and shrink the effective radii of galaxies by 44% which could increase the completeness fraction of galaxy surveys. Finally, our simulations show that surveys with a magnitude limit of ${\rm M_{UV,1500{\rm A}}=-13}$ can probe the sources that provide the 50% of photons for reionization out to $z=12$.