We present a new model for high redshift Lyman α emitters (LAEs) in the cosmological context which takes into account the resonant scattering of Lyα 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, in post-processing, a gas outflow model for each galaxy based on simple scaling arguments. The coupling with a library of numerical experiments of Lyα transfer through expanding (or static) dusty shells of gas allows us to derive the Lyα escape fraction and profile of each galaxy. Results obtained with this new approach are compared with simpler models often used in the literature. The predicted distribution of Lyα photons escape fraction shows that galaxies with a low star formation rate (SFR) have a fesc of the order of unity, suggesting that, for those objects, Lyα may be used to trace the SFR 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 ultraviolet (UV) and Lyα luminosity function data at 3 < z < 5. We find that we are in good agreement with both the bright Lyα data and the faint LAE population observed by Rauch et al. at z= 3 whereas a simpler constant Lyαescape fraction model fails to do so. Most of the Lyα profiles of our LAEs are redshifted by the diffusion in the expanding gas which suppresses intergalactic medium absorption and scattering. The bulk of the observed Lyα equivalent width (EW) distribution is recovered by our model, but we fail to obtain the very large values sometimes detected. Our predictions for stellar masses and UV luminosity functions of LAEs show a satisfactory agreement with observational estimates. The UV-brightest galaxies are found to show only low Lyα 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 H i column densities that we predict for these objects, which quench preferentially resonant Lyα photons via dust extinction.
We present a clustering analysis of a sample of 238 Ly alpha emitters at redshift 3 less than or similar to z less than or similar to 6 from the MUSE-Wide survey. This survey mosaics extragalactic ...
We present a study of the galaxy Lyman-alpha luminosity function (LF) using a sample of 17 lensing clusters observed by the MUSE/VLT. Magnification from strong gravitational lensing by clusters of galaxies and MUSE apabilities allow us to blindly detect LAEs without any photometric pre-selection, reaching the faint luminosity regime. 600 lensed LAEs were selected behind these clusters in the redshift range 2.9<$z$< 6.7, covering four orders of magnitude in magnification-corrected Lyman-alpha luminosity (39.042 are consistent with those obtained from blank field observations. In the faint luminosity regime, the density of sources is well described by a steep slope, $α\sim-2$ for the global redshift range. Up to log(L)$\sim$41, the steepening of the faint end slope with redshift, suggested by the earlier work of DLV19 is observed, but the uncertainties remain large. A significant flattening is observed towards the faintest end, for the highest redshift bins (log$L$<41). Using face values, the steep slope at the faint-end causes the SFRD to dramatically increase with redshift, implying that LAEs could play a major role in the process of cosmic reionization. The flattening observed towards the faint end for the highest redshift bins still needs further investigation. This turnover is similar to the one observed for the UV LF at $z\geq6$ in lensing clusters, with the same conclusions regarding the reliability of current results (e.g.arXiv:1803.09747(N); arXiv:2205.11526(N)).
Context. The Ly α emitter (LAE) fraction, X LAE , is a potentially powerful probe of the evolution of the intergalactic neutral hydrogen gas fraction. However, uncertainties in the measurement of X LAE are still under debate. Aims. Thanks to deep data obtained with the integral field spectrograph Multi Unit Spectroscopic Explorer (MUSE), we can measure the evolution of the LAE fraction homogeneously over a wide redshift range of z ≈ 3–6 for UV-faint galaxies (down to UV magnitudes of M 1500 ≈ −17.75). This is a significantly fainter range than in former studies ( M 1500 ≤ −18.75) and it allows us to probe the bulk of the population of high-redshift star-forming galaxies. Methods. We constructed a UV-complete photometric-redshift sample following UV luminosity functions and measured the Ly α emission with MUSE using the latest (second) data release from the MUSE Hubble Ultra Deep Field Survey. Results. We derived the redshift evolution of X LAE for M 1500 ∈ [ − 21.75; −17.75] for the first time with a equivalent width range EW (Ly α ) ≥ 65 Å and found low values of X LAE ≲ 30% at z ≲ 6. The best-fit linear relation is X LAE = 0.07 +0.06 −0.03 z − 0.22 +0.12 −0.24 . For M 1500 ∈ [ − 20.25; −18.75] and EW (Ly α ) ≥ 25 Å, our X LAE values are consistent with those in the literature within 1 σ at z ≲ 5, but our median values are systematically lower than reported values over the whole redshift range. In addition, we do not find a significant dependence of X LAE on M 1500 for EW (Ly α ) ≥ 50 Å at z ≈ 3–4, in contrast with previous work. The differences in X LAE mainly arise from selection biases for Lyman Break Galaxies (LBGs) in the literature: UV-faint LBGs are more easily selected if they have strong Ly α emission, hence X LAE is biased towards higher values when those samples are used. Conclusions. Our results suggest either a lower increase of X LAE towards z ≈ 6 than previously suggested, or even a turnover of X LAE at z ≈ 5.5, which may be the signature of a late or patchy reionization process. We compared our results with predictions from a cosmological galaxy evolution model. We find that a model with a bursty star formation (SF) can reproduce our observed LAE fractions much better than models where SF is a smooth function of time.
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.
The Lyman continuum (LyC) cannot be observed at the epoch of reionization (z {\gtrsim} 6) due to intergalactic H I absorption. To identify Lyman continuum emitters (LCEs) and infer the fraction of escaping LyC, astronomers have developed various indirect diagnostics of LyC escape. Using measurements of the LyC from the Low-redshift Lyman Continuum Survey (LzLCS), we present the first statistical test of these diagnostics. While optical depth indicators based on Lyα, such as peak velocity separation and equivalent width, perform well, we also find that other diagnostics, such as the [O III]/[O II] flux ratio and star formation rate surface density, predict whether a galaxy is a LCE. The relationship between these galaxy properties and the fraction of escaping LyC flux suggests that LyC escape depends strongly on H I column density, ionization parameter, and stellar feedback. We find LCEs occupy a range of stellar masses, metallicities, star formation histories, and ionization parameters, which may indicate episodic and/or different physical causes of LyC escape.
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.
We present an analysis of HI Lyman-alpha emission in deep VLT/MUSE observations of two highly magnified and extended galaxies at z=3.5 and 4.03, including a newly discovered, almost complete Einstein ring. While these Lyman-alpha haloes are intrinsically similar to the ones typically seen in other MUSE deep fields, the benefits of gravitational lensing allows us to construct exceptionally detailed maps of Lyman-alpha line properties at sub-kpc scales. By combining all multiple images, we are able to observe complex structures in the Lyman-alpha emission and uncover small (~ 120 km/s in Lyman-alpha peak shift), but significant at > 4 sigma, systematic variations in the shape of the Lyman-alpha line profile within each halo. Indeed, we observe a global trend for the line peak shift to become redder at large radii, together with a strong correlation between the peak wavelength and line width. This systematic intrahalo variation is markedly similar to the object-to-object variations obtained from the integrated properties of recent large samples. Regions of high surface brightness correspond to relatively small line shifts, which could indicate that Lyman-alpha emission escapes preferentially from regions where the line profile has been less severely affected by scattering of Lyman-alpha photons.
Abstract Understanding the origin of strong galactic outflows and the suppression of star formation in dwarf galaxies is a key problem in galaxy formation. Using a set of radiation-hydrodynamic simulations of an isolated dwarf galaxy embedded in a 1010 M⊙ halo, we show that the momentum transferred from resonantly scattered Lyman-α (Lyα) photons is an important source of stellar feedback which can shape the evolution of galaxies. We find that Lyα feedback suppresses star formation by a factor of two in metal-poor galaxies by regulating the dynamics of star-forming clouds before the onset of supernova explosions (SNe). This is possible because each Lyα photon resonantly scatters and imparts ∼10–300 times greater momentum than in the single scattering limit. Consequently, the number of star clusters predicted in the simulations is reduced by a factor of ∼5, compared to the model without the early feedback. More importantly, we find that galactic outflows become weaker in the presence of strong Lyα radiation feedback, as star formation and associated SNe become less bursty. We also examine a model in which radiation field is arbitrarily enhanced by a factor of up to 10, and reach the same conclusion. The typical mass-loading factors in our metal-poor dwarf system are estimated to be ∼5–10 near the mid-plane, while it is reduced to ∼1 at larger radii. Finally, we find that the escape of ionizing radiation and hence the reionization history of the Universe is unlikely to be strongly affected by Lyα feedback.
'%3e%3cpath fill='%23012169' d='M0 0v30h60V0z'/%3e%3cpath stroke='%23fff' stroke-width='6' d='m0 0 60 30m0-30L0 30'/%3e%3cpath stroke='%23C8102E' stroke-width='4' d='m0 0 60 30m0-30L0 30' clip-path='url(%23b)'/%3e%3cpath stroke='%23fff' stroke-width='10' d='M30 0v30M0 15h60'/%3e%3cpath stroke='%23C8102E' stroke-width='6' d='M30 0v30M0 15h60'/%3e%3c/g%3e%3c/svg%3e)
