Source author record

Jérémy Blaizot

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.CO astro-ph.IM Distributed, Parallel, and Cluster Computing

Catalog footprint

Research overview

13works

4topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

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

Radiation-MagnetoHydrodynamics simulations of cosmic ray feedback in disc galaxies

Cosmic rays (CRs) are thought to play an important role in galaxy evolution. We study their effect when coupled to other important sources of feedback, namely supernovae and stellar radiation, by including CR anisotropic diffusion and radiative losses but neglecting CR streaming. Using the RAMSES-RT code, we perform the first radiation-magnetohydrodynamics simulations of isolated disc galaxies with and without CRs. We study galaxies embedded in dark matter haloes of $10^{10}$ , $10^{11}$ and $10^{12}\, \rm M_{\odot}$ with a maximum resolution of $9 \,\rm pc$ . We find that CRs reduce star formation rate in our two dwarf galaxies by a factor 2, with decreasing efficiency with increasing galaxy mass. They increase significantly the outflow mass loading factor in all our galaxies and make the outflows colder. We study the impact of the CR diffusion coefficient, exploring values from $κ= 10^{27}$ to $\rm 3\times 10^{29}\, cm^2\, s^{-1}$ . With lower $κ$ , CRs remain confined for longer on small scales and are consequently efficient in suppressing star formation, whereas a higher diffusion coefficient reduces the effect on star formation and increases the generation of cold outflows. Finally, we compare CR feedback to a calibrated 'strong' supernova feedback model known to sufficiently regulate star formation in high-redshift cosmological simulations. We find that CR feedback is not sufficiently strong to replace this strong supernova feedback. As they tend to smooth out the ISM and fill it with denser gas, CRs also lower the escape fraction of Lyman continuum photons from galaxies.

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)

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

Tracing the simulated high-redshift circum-galactic medium with Lyman alpha emission

With the Multi Unit Spectroscopic Explorer (MUSE), it is now possible to detect spatially extended Lyman alpha emission from individual faint (M_UV ~ -18) galaxies at redshifts, 3 < z < 6, tracing gas out to circum-galactic scales comparable to the dark matter halo virial radius. To explore the implications of such observations, we present a cosmological radiation hydrodynamics simulation of a single galaxy, chosen to be typical of the Lyman alpha-emitting galaxies detected by MUSE in deep fields. We use this simulation to study the origin and dynamics of the high-redshift circum-galactic medium (CGM). We find that the majority of the mass in the diffuse CGM is comprised of material infalling for the first time towards the halo center, but with the inner CGM also containing a comparable amount of mass that has moved past first-pericentric passage, and is in the process of settling into a rotationally supported configuration. Making the connection to Lyman alpha emission, we find that the observed extended surface brightness profile is due to a combination of three components: scattering of galactic Lyman alpha emission in the CGM, in-situ emission of CGM gas (mostly infalling), and Lyman alpha emission from small satellite galaxies. The weight of these contributions vary with distance from the galaxy such that (1) scattering dominates the inner regions (r < 7 kpc), at surface brightness larger than a few 10^-19 cgs, (2) all components contribute equally around r ~ 10 kpc (or SB ~10^-19), and (3) the contribution of small satellite galaxies takes over at large distances (or SB ~10^-20). Our simulation fails to reproduce the characteristic observed Lyman alpha spectral morphology that is red-shifted with respect to the systemic velocity, with the implication that the simulation is missing an important component of neutral outflowing gas.

preprint2020arXiv

The MUSE Hubble Ultra Deep Field Survey XIII. Spatially resolved spectral properties of Lyman alpha haloes around star-forming galaxies at z > 3

We present spatially resolved maps of six individually-detected Lyman alpha haloes (LAHs) as well as a first statistical analysis of the Lyman alpha (Lya) spectral signature in the circum-galactic medium of high-redshift star-forming galaxies using MUSE. Our resolved spectroscopic analysis of the LAHs reveals significant intrahalo variations of the Lya line profile. Using a three-dimensional two-component model for the Lya emission, we measure the full width at half maximum (FWHM), the peak velocity shift and the asymmetry of the Lya line in the core and in the halo of 19 galaxies. We find that the Lya line shape is statistically different in the halo compared to the core for ~40% of our galaxies. Similarly to object-by-object based studies and a recent resolved study using lensing, we find a correlation between the peak velocity shift and the width of the Lya line both at the interstellar and circum-galactic scales. While there is a lack of correlation between the spectral properties and the spatial scale lengths of our LAHs, we find a correlation between the width of the line in the LAH and the halo flux fraction. Interestingly, UV bright galaxies show broader, more redshifted and less asymmetric Lya lines in their haloes. The most significant correlation found is for the FWHM of the line and the UV continuum slope of the galaxy, suggesting that the redder galaxies have broader Lya lines. The generally broad and red line shapes found in the halo component suggests that the Lya haloes are powered either by scattering processes through an outflowing medium, fluorescent emission from outflowing cold clumps of gas, or a mix of both. Considering the large diversity of the Lya line profiles observed in our sample and the lack of strong correlation, the interpretation of our results is still broadly open and underlines the need for realistic spatially resolved models of the LAHs.

preprint2016arXiv

Modeling Lyman-α Forest Cross-Correlations with LyMAS

We use the Ly- $α$ Mass Association Scheme (LyMAS; Peirani et al. 2014) to predict cross-correlations at $z=2.5$ between dark matter halos and transmitted flux in the Ly- $α$ forest, and compare to cross-correlations measured for quasars and damped Ly- $α$ systems (DLAs) from the Baryon Oscillation Spectroscopic Survey (BOSS) by Font-Ribera et al. (2012, 2013). We calibrate LyMAS using Horizon-AGN hydrodynamical cosmological simulations of a $(100\ h^{-1}\ \mathrm{Mpc})^3$ comoving volume. We apply this calibration to a $(1\ h^{-1}\ \mathrm{Gpc})^3$ simulation realized with $2048^3$ dark matter particles. In the 100 $h^{-1}$ Mpc box, LyMAS reproduces the halo-flux correlations computed from the full hydrodynamic gas distribution very well. In the 1 $h^{-1}$ Gpc box, the amplitude of the large scale cross-correlation tracks the halo bias $b_h$ as expected. We provide empirical fitting functions that describe our numerical results. In the transverse separation bins used for the BOSS analyses, LyMAS cross-correlation predictions follow linear theory accurately down to small scales. Fitting the BOSS measurements requires inclusion of random velocity errors; we find best-fit RMS velocity errors of 399 km s $^{-1}$ and 252 km s $^{-1}$ for quasars and DLAs, respectively. We infer bias-weighted mean halo masses of $M_h/10^{12}\ h^{-1}M_\odot=2.19^{+0.16}_{-0.15}$ and $0.69^{+0.16}_{-0.14}$ for the host halos of quasars and DLAs, with $\sim 0.2$ dex systematic uncertainty associated with redshift evolution, IGM parameters, and selection of data fitting range.

preprint2015arXiv

Lyman-α Emitters in the context of hierarchical galaxy formation: predictions for VLT/MUSE surveys

The VLT Multi Unit Spectroscopic Explorer (MUSE) integral-field spectrograph can detect Lyα emitters (LAE) in the redshift range $2.8 \lesssim z \lesssim 6.7$ in a homogeneous way. Ongoing MUSE surveys will notably probe faint Lyα sources that are usually missed by current narrow-band surveys. We provide quantitative predictions for a typical wedding-cake observing strategy with MUSE based on mock catalogs generated with a semi-analytic model of galaxy formation coupled to numerical Lyα radiation transfer models in gas outflows. We expect $\approx$ 1500 bright LAEs ( $F_{Lyα}$ $\gtrsim$ $10^{-17}$ erg s $^{-1}$ cm $^{-2}$ ) in a typical Shallow Field (SF) survey carried over $\approx$ 100 arcmin $^2$ , and $\approx$ 2,000 sources as faint as $10^{-18}$ erg s $^{-1}$ cm $^{-2}$ in a Medium-Deep Field (MDF) survey over 10 arcmin $^2$ . In a typical Deep Field (DF) survey of 1 arcmin $^2$ , we predict that $\approx$ 500 extremely faint LAEs ( $F_{Lyα}$ $\gtrsim$ $4 \times 10^{-19}$ erg s $^{-1}$ cm $^{-2}$ ) will be found. Our results suggest that faint Lyα sources contribute significantly to the cosmic Lyα luminosity and SFR budget. While the host halos of bright LAEs at z $\approx$ 3 and 6 have descendants with median masses of $2 \times 10^{12}$ and $5 \times 10^{13}$ $M_{\odot}$ respectively, the faintest sources detectable by MUSE at these redshifts are predicted to reside in halos which evolve into typical sub- $L^{*}$ and $L^{*}$ galaxy halos at z = 0. We expect typical DF and MDF surveys to uncover the building blocks of Milky Way-like objects, even probing the bulk of the stellar mass content of LAEs located in their progenitor halos at z $\approx$ 3.

preprint2014arXiv

LyMAS: Predicting Large-Scale Lyman-alpha Forest Statistics from the Dark Matter Density Field

[abridged] We describe LyMAS (Ly-alpha Mass Association Scheme), a method of predicting clustering statistics in the Ly-alpha forest on large scales from moderate resolution simulations of the dark matter distribution, with calibration from high-resolution hydrodynamic simulations of smaller volumes. We use the "Horizon MareNostrum" simulation, a 50 Mpc/h comoving volume evolved with the adaptive mesh hydrodynamic code RAMSES, to compute the conditional probability distribution P(F_s|delta_s) of the transmitted flux F_s, smoothed (1-dimensionally) over the spectral resolution scale, on the dark matter density contrast delta_s, smoothed (3-dimensionally) over a similar scale. In this study we adopt the spectral resolution of the SDSS-III BOSS at z=2.5, and we find optimal results for a dark matter smoothing length sigma=0.3 Mpc/h (comoving). In extended form, LyMAS exactly reproduces both the 1-dimensional power spectrum and 1-point flux distribution of the hydro simulation spectra. Applied to the MareNostrum dark matter field, LyMAS accurately predicts the 2-point conditional flux distribution and flux correlation function of the full hydro simulation for transverse sightline separations as small as 1 Mpc/h, including redshift-space distortion effects. It is substantially more accurate than a deterministic density-flux mapping ("Fluctuating Gunn-Peterson Approximation"), often used for large volume simulations of the forest. With the MareNostrum calibration, we apply LyMAS to 1024^3 N-body simulations of a 300 Mpc/h and 1.0 Gpc/h cube to produce large, publicly available catalogs of mock BOSS spectra that probe a large comoving volume. LyMAS will be a powerful tool for interpreting 3-d Ly-alpha forest data, thereby transforming measurements from BOSS and other massive quasar absorption surveys into constraints on dark energy, dark matter, space geometry, and IGM physics.

preprint2012arXiv

Modelling high redshift Lyman-alpha Emitters

We present a new model for high redshift Lyman-Alpha Emitters (LAEs) in the cosmological context which takes into account the resonant scattering of Ly-a photons through expanding gas. The GALICS semi-analytic model provides us with the physical properties of a large sample of high redshift galaxies. We implement a gas outflow model for each galaxy based on simple scaling arguments. The coupling with a library of numerical experiments of Ly-a transfer through expanding or static dusty shells of gas allows us to derive the Ly-a escape fractions and profiles. The predicted distribution of Ly-a photons escape fraction shows that galaxies with a low star formation rate have a f_esc of the order of unity, suggesting that, for those objects, Ly-a may be used to trace the star formation rate assuming a given conversion law. In galaxies forming stars intensely, the escape fraction spans the whole range from 0 to 1. The model is able to get a good match to the UV and Ly-a luminosity function (LF) data at 3 < z < 5. We find that we are in good agreement with both the bright Ly-a data and the faint population observed by Rauch et al. (2008) at z=3. Most of the Ly-a profiles of our LAEs are redshifted by the diffusion in the outflow which suppresses IGM absorption. The bulk of the observed Ly-a equivalent width (EW) distribution is recovered by our model, but we fail to obtain the very large values sometimes detected. Predictions for stellar masses and UV LFs of LAEs show a satisfactory agreement with observational estimates. The UV-brightest galaxies are found to show only low Ly-a EWs in our model, as it is reported by many observations of high redshift LAEs. We interpret this effect as the joint consequence of old stellar populations hosted by UV-bright galaxies, and high HI column densities that we predict for these objects, which quench preferentially resonant Ly-a photons via dust extinction.

preprint2010arXiv

Cosmological Simulations on a Grid of Computers

The work presented in this paper aims at restricting the input parameter values of the semi-analytical model used in GALICS and MOMAF, so as to derive which parameters influence the most the results, e.g., star formation, feedback and halo recycling efficiencies, etc. Our approach is to proceed empirically: we run lots of simulations and derive the correct ranges of values. The computation time needed is so large, that we need to run on a grid of computers. Hence, we model GALICS and MOMAF execution time and output files size, and run the simulation using a grid middleware: DIET. All the complexity of accessing resources, scheduling simulations and managing data is harnessed by DIET and hidden behind a web portal accessible to the users.

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.

Jérémy Blaizot

Research overview

Is this your work?

13 publications

LyC escape from SPHINX galaxies in the Epoch of Reionization

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

Radiation-MagnetoHydrodynamics simulations of cosmic ray feedback in disc galaxies

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

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?

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

Tracing the simulated high-redshift circum-galactic medium with Lyman alpha emission

The MUSE Hubble Ultra Deep Field Survey XIII. Spatially resolved spectral properties of Lyman alpha haloes around star-forming galaxies at z > 3

Modeling Lyman-α Forest Cross-Correlations with LyMAS

Lyman-α Emitters in the context of hierarchical galaxy formation: predictions for VLT/MUSE surveys

LyMAS: Predicting Large-Scale Lyman-alpha Forest Statistics from the Dark Matter Density Field

Modelling high redshift Lyman-alpha Emitters

Cosmological Simulations on a Grid of Computers

See the researcher in context

Building this map preview

Jérémy Blaizot

Research overview

Is this your work?

13 publications

LyC escape from SPHINX galaxies in the Epoch of Reionization

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

Radiation-MagnetoHydrodynamics simulations of cosmic ray feedback in disc galaxies

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

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

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

Tracing the simulated high-redshift circum-galactic medium with Lyman alpha emission

The MUSE Hubble Ultra Deep Field Survey XIII. Spatially resolved spectral properties of Lyman alpha haloes around star-forming galaxies at z > 3

Modeling Lyman-α Forest Cross-Correlations with LyMAS

Lyman-α Emitters in the context of hierarchical galaxy formation: predictions for VLT/MUSE surveys

LyMAS: Predicting Large-Scale Lyman-alpha Forest Statistics from the Dark Matter Density Field

Modelling high redshift Lyman-alpha Emitters

Cosmological Simulations on a Grid of Computers

See the researcher in context

Building this map preview

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?

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