Search

The Ly$\alpha$ emission line is one of the most promising probes of cosmic reionisation but isolating the signature of a change in the ionisation state of the IGM is challenging because of intrinsic evolution and internal radiation transfer effects. We present the first study of the evolution of Ly$\alpha$ emitters (LAE) during the epoch of reionisation based on a full radiation-hydrodynamics cosmological simulation that is able to capture both the large-scale process of reionisation and the small-scale properties of galaxies. We predict the Ly$\alpha$ emission of galaxies in the $10^3$ cMpc$^3$ SPHINX simulation at $6\leq z\leq9$ by computing the full Ly$\alpha$ radiation transfer from ISM to IGM scales. SPHINX is able to reproduce many observational constraints such as the UV/Ly$\alpha$ luminosity functions and stellar mass functions at z $\geq$ 6 for the dynamical range probed by our simulation ($M_{\rm 1500}\gtrsim-18$, $L_{\rm Ly\alpha}\lesssim10^{42}$ erg/s, $M_{\star}\lesssim10^9$ M$_{\odot}$). As intrinsic Ly$\alpha$ emission and internal Ly$\alpha$ escape fractions barely evolve from $z=6$ to 9, the observed suppression of Ly$\alpha$ luminosities with increasing redshift is fully attributed to IGM absorption. For most observable galaxies ($M_{\rm 1500}\lesssim-16$), the Ly$\alpha$ line profiles are slightly shifted to the red due to internal radiative transfer effects which mitigates the effect of IGM absorption. Overall, the enhanced Ly$\alpha$ suppression during reionisation traces the IGM neutral fraction $x_{\rm HI}$ well but the predicted amplitude of this reduction is a strong function of the Ly$\alpha$ peak shift, which is set at ISM/CGM scales. We find that a large number of LAEs could be detectable in very deep surveys during reionisation when $x_{\rm HI}$ is still $\approx 50\%$., Comment: MNRAS accepted, 20 pages, 19 figures, plus appendices

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}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}, 21 pages, 11 figures, MNRAS Accepted

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 multi-frequency, on-the-fly radiation transfer and constrained transport ideal MHD in a cosmological context to simultaneously model the inhomogeneous process of reionization as well as the growth of PMFs. We run a series of $(5\,\text{cMpc})^3$ cosmological volumes, varying both the strength of the seed magnetic field ($B_0$) and its spectral index ($n_B$). We find that PMFs that have $n_B > -0.562\log_{10}\left(\frac{B_0}{1{\rm n}G}\right) - 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 Lya 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 sensitive indicator of the properties of the PMFs, even for those that satisfy our $\tau_e$ constraint. Furthermore, strong PMFs can marginally increase the escape of LyC photons by up to 25\% and shrink the effective radii of galaxies by $\sim44\%$ which could increase the completeness fraction of galaxy surveys. Finally, our simulations show that surveys with a magnitude limit of ${\rm M_{UV,1500}=-13}$ can probe the sources that provide the majority of photons for reionization out to $z=12$., Comment: 29 pages, 27 figures, accepted MNRAS

Infrared and nebular lines provide some of our best probes of the physics regulating the properties of the interstellar medium (ISM) at high-redshift. However, interpreting the physical conditions of high-redshift galaxies directly from emission lines remains complicated due to inhomogeneities in temperature, density, metallicity, ionisation parameter, and spectral hardness. We present a new suite of cosmological, radiation-hydrodynamics simulations, each centred on a massive Lyman-break galaxy that resolves such properties in an inhomogeneous ISM. Many of the simulated systems exhibit transient but well defined gaseous disks that appear as velocity gradients in [CII]~158.6$\mu$m emission. Spatial and spectral offsets between [CII]~158.6$\mu$m and [OIII]~88.33$\mu$m are common, but not ubiquitous, as each line probes a different phase of the ISM. These systems fall on the local [CII]-SFR relation, consistent with newer observations that question previously observed [CII]~158.6$\mu$m deficits. Our galaxies are consistent with the nebular line properties of observed $z\sim2-3$ galaxies and reproduce offsets on the BPT and mass-excitation diagrams compared to local galaxies due to higher star formation rate (SFR), excitation, and specific-SFR, as well as harder spectra from young, metal-poor binaries. We predict that local calibrations between H$\alpha$ and [OII]~3727$\AA$ luminosity and galaxy SFR apply up to $z>10$, as do the local relations between certain strong line diagnostics (R23 and [OIII]~5007$\AA$/H$\beta$) and galaxy metallicity. Our new simulations are well suited to interpret the observations of line emission from current (ALMA and HST) and upcoming facilities (JWST and ngVLA)., Comment: 21 pages, 16 figures, MNRAS accepted

Galaxy clusters constitute powerful cosmological probes thanks to comparisons between observed and simulated clusters. As such Virgo constitutes a formidable source for detailed observations facilitated by its proximity. However, the diversity of clusters complicates the comparisons on a one-to-one basis. Simulated clusters must be carefully selected, a daunting task since most properties are unknown. Alternatively, lookalikes produced in the proper large scale environment can be used. Additionally, their statistical study give access to the mean properties of the observed cluster including its most probable history as well as its deviation from an average cluster. This paper presents such a statistical study with 200 Virgo-like and 400+ cluster-size random dark matter halos. Only 18%(0.5%) of these random halos comply within 3(2)-sigma with the mean values (radius, velocity dispersion, number of substructures, spin, velocity, concentration, center of mass offset with respect to the spherical center) of Virgo halos at z=0 and abide by a similar merging history up to redshift 4. None are within 1-sigma because of environmentally induced properties (number of substructures and velocity). For further comparisons, random halos are selected to reproduce the mass distribution of the lookalikes to cancel mass bias effects. Redshift 1 appears then as a turning point: random to Virgo-like property ratios are alternatively smaller/larger than 1. This highlights the importance of studying clusters within their proper large scale environment: simulated galaxy population, grandly affected by the cluster history, can then be compared with the observed one in details. Direct lookalikes simplify grandly the challenge., Accepted for publication in MNRAS, 12 pages, 7 figures, 5 tables

At ~16-17Mpc from us, the Virgo cluster is a formidable source of information to study cluster formation and galaxy evolution in rich environments. Several observationally-driven formation scenarios arose within the past decade to explain the properties of galaxies that entered the cluster recently and the nature of the last significant merger that the cluster underwent. Confirming these scenarios requires extremely faithful numerical counterparts of the cluster. This paper presents the first CLONE, Constrained LOcal and Nesting Environment, simulation of the Virgo cluster within a ~15Mpc radius sphere. This cosmological hydrodynamical simulation, with feedback from supernovae and active galactic nuclei, with a ~3x10^7Msun dark matter particle mass and a minimum cell size of 350pc in the zoom region, reproduces Virgo within its large scale environment unlike a random cluster simulation. Overall the distribution of the simulated galaxy population matches the observed one including M87. The simulated cluster formation reveals exquisite agreements with observationally-driven scenarios: within the last Gigayear, about 300 small galaxies (M*>10^7Msun) entered the cluster, most of them within the last 500Myr. The last significant merger event occurred about 2 Gigayears ago: a group with a tenth of the mass of today's cluster entered from the far side as viewed from the Milky Way. This excellent numerical replica of Virgo will permit studying different galaxy type evolution (jellyfish, backsplash, etc.) as well as feedback phenomena in the cluster core via unbiased comparisons between simulated and observed galaxies and hot gas phase profiles to understand this great physics laboratory., Accepted for publication in MNRAS, 15 pages, 14 figures, 4 tables

HARMONI is the adaptive optics assisted, near-infrared and visible light integral field spectrograph for the Extremely Large Telescope (ELT). A first light instrument, it provides the work-horse spectroscopic capability for the ELT. As the project approaches its Final Design Review milestone, the design of the instrument is being finalized, and the plans for assembly, integration and testing are being detailed. We present an overview of the instrument's capabilities from a user perspective, provide a summary of the instrument's design, including plans for operations and calibrations, and provide a brief glimpse of the predicted performance for a specific observing scenario. The paper also provides some details of the consortium composition and its evolution since the project commenced in 2015.

The neutral intergalactic medium above redshift 6 is opaque to ionizing radiation, therefore one needs indirect measurements of the escape fraction of ionizing photons from galaxies of this epoch. Low-ionization state absorption lines are a common feature in the spectrum of galaxies, showing a diversity of strengths and shapes. Since these lines indicate the presence of neutral gas in front of the stars, they have been proposed to carry information on the escape of ionizing radiation from galaxies. We study which processes are responsible for the shape of the absorption lines, to better understand their origin. We then explore whether the absorption lines can be used to predict the escape fractions. Using a radiation-hydrodynamical zoom-in simulation and the radiative transfer code RASCAS, we generate mock CII 1334 and LyB lines of a virtual galaxy at redshift 3 as seen from many directions of observation. We also compute the escape fraction of ionizing photons in those directions and look for correlations between the lines and the escape fractions. We find that the resulting mock absorption lines are comparable to observations and that the lines and the escape fractions vary strongly depending on the direction of observation. Gas velocity and dust always affect the absorption profile significantly. We find no strong correlations between observable LyB or CII 1334 and the escape fraction. After correcting the continuum for attenuation by dust to recover the intrinsic continuum, the residual flux of CII 1334 correlates well with the escape fraction for directions with a dust corrected residual flux larger than 30%. For other directions, the relations have a strong dispersion, and the residual flux overestimates the escape fraction for most cases. Concerning LyB, the residual flux after dust correction does not correlate with the escape fraction but can be used as a lower limit. (abridged), 25 pages, 22 figures, accepted for publication in Astronomy and Astrophysics. New title, shorter abstract, corrected typos and url for Krome website

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., MNRAS accepted, 20 pages, 10 figures

The Lya 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 debated. Thanks to deep data obtained with MUSE, we can measure the evolution of X_LAE homogeneously over a wide redshift range of z~3-6 for UV-faint galaxies (down to M_1500~-17.75). This is significantly fainter than in former studies, and allows us to probe the bulk of the population of high-z star-forming galaxies. We construct a UV-complete photo-redshift sample following UV luminosity functions and measure the Lya emission with MUSE using the second data release from the MUSE HUDF Survey. We derive the redshift evolution of X_LAE for M_1500 in [-21.75;-17.75] for the first time with a equivalent width range EW(Lya)>=65 A and find low values of X_ LAE~50 A 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 Lya emission, hence X_LAE is biased towards higher values. 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 X_LAE much better than models where SF is a smooth function of time., Comment: 25 pages, 17 figures, 5 tables, accepted for publication in A&A

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., Comment: Accepted for publication in A&A. Code available at http://rascas.univ-lyon1.fr/

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\mu m}$ and [OIII]$_{\rm 88\mu 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., Comment: 18 pages, 11 figures, Accepted to MNRAS

We present velocities of galactic outflows in seven star-forming galaxies at $z=$5-6 with stellar masses of $M_* \sim10^{10.1}\ \rm{M_\odot}$. Although it is challenging to observationally determine the outflow velocities, we overcome this by using ALMA [CII]158 $\mu$m emission lines for systemic velocities and deep Keck spectra with metal absorption lines for velocity profiles available to date. We construct a composite Keck spectrum of the galaxies at $z=$5-6 with the [CII]-systemic velocities, and fit outflow-line profiles to the SiII1260, CII1335, and SiIV1394,1403 absorption lines in the composite spectrum. We measure the maximum (90%) and central outflow velocities to be $v_\rm{max}=700^{+180}_{-110}\ \rm{km\ s^{-1}}$ and $v_\rm{out}= 400^{+100}_{-150}\ \rm{km\ s^{-1}}$ on average, respectively, showing no significant differences between the outflow velocities derived with the low to high-ionization absorption lines. For $M_* \sim10^{10.1}\ \rm{M_\odot}$, we find that the $v_\rm{max}$ value of our $z=$5-6 galaxies is 3 times higher than those of $z\sim0$ galaxies and comparable to $z\sim2$ galaxies. Estimating the halo circular velocity $v_\rm{cir}$ from the stellar masses and the abundance matching results, we investigate a $v_\rm{max}$-$v_\rm{cir}$ relation. Interestingly, $v_\rm{max}$ for galaxies with $M_*=10^{10.0-10.8}\ \rm{M_\odot}$ shows a clear positive correlation with $v_\rm{cir}$ and/or the galaxy star formation rate over $z=$0-6 with a small scatter of $\simeq \pm 0.1$ dex, which is in good agreement with theoretical predictions. This positive correlation suggests that the outflow velocity is physically related to the halo circular velocity, and that the redshift evolution of $v_\rm{max}$ at fixed $M_*$ is explained by the increase in $v_\rm{cir}$ toward high redshift., Comment: 13 pages, 7 figures, 4 tables; Accepted for publication in ApJ

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., 9 pages, 4 figures, accepted for publication in MNRAS

We simulate the flux emitted from galaxy halos in order to quantify the brightness of the circumgalactic medium (CGM). We use dedicated zoom-in cosmological simulations with the hydrodynamical Adaptive Mesh Refinement code RAMSES, which are evolved down to z=0 and reach a maximum spatial resolution of 380 $h^{-1}$pc and a gas mass resolution up to 1.8$\times 10^{5} h^{-1} \rm{M}_{\odot}$ in the densest regions. We compute the expected emission from the gas in the CGM using CLOUDY emissivity models for different lines (e.g. Ly$\alpha$, CIV, OVI, CVI, OVIII) considering UV background fluorescence, gravitational cooling and continuum emission. In the case of Ly$\alpha$ we additionally consider the scattering of continuum photons. We compare our predictions to current observations and find them to be in good agreement at any redshift after adjusting the Ly$\alpha$ escape fraction. We combine our mock observations with instrument models for FIREBall-2 (UV balloon spectrograph) and HARMONI (visible and NIR IFU on the ELT) to predict CGM observations with either instrument and optimise target selections and observing strategies. Our results show that Ly$\alpha$ emission from the CGM at a redshift of 0.7 will be observable with FIREBall-2 for bright galaxies (NUV$\sim$18 mag), while metal lines like OVI and CIV will remain challenging to detect. HARMONI is found to be well suited to study the CGM at different redshifts with various tracers., Comment: accepted for publication in MNRAS

Understanding the escape of Lyman continuum (LyC) and Lyman alpha (Lya) photons from molecular clouds is one of the keys to constraining the reionization history of the Universe and the evolution of galaxies at high redshift. Using a set of radiation-hydrodynamic simulations with adaptive mesh refinement, we investigate how photons propagate and escape from turbulent clouds with different masses, star formation efficiencies (SFEs), and metallicities, as well as with different models of stellar spectra and supernova feedback. We find that the escape fractions in both LyC and Lya are generally increasing with time if the cloud is efficiently dispersed by radiation and supernova feedback. When the total SFE is low (1% of the cloud mass), 0.1-5% of LyC photons leave the metal-poor cloud, whereas the fractions increase to 20-70% in clouds with a 10% SFE. LyC photons escape more efficiently if gas metallicity is lower, if the upper mass limit in the stellar initial mass function is higher, if binary interactions are allowed in the evolution of stars, or if additional strong radiation pressure, such as Lya pressure, is present. As a result, the number of escaping LyC photons can easily vary by a factor of $\sim4$ on cloud scales. The escape fractions of Lya photons are systemically higher (60-80%) than those of LyC photons despite large optical depths at line centre ($\tau_0\sim10^6-10^9$). Scattering of Lya photons is already significant on cloud scales, leading to double-peaked profiles with peak separations of $v_{\rm sep}\sim400\,{\rm km\,s^{-1}}$ during the initial stage of the cloud evolution, while it becomes narrower than $v_{\rm sep} \le 150 \, {\rm km\,s^{-1}}$ in the LyC bright phase. Comparisons with observations of low-redshift galaxies suggest that Lya photons require further interactions with neutral hydrogen to reproduce their velocity offset for a given LyC escape fraction., Comment: 23 pages, 20 figures, accepted for publication in MNRAS

We use the photon flux and absorption tracer algorithm presented in Katz et al.~2018, to characterise the contribution of haloes of different mass and stars of different age and metallicity to the reionization of the Universe. We employ a suite of cosmological multifrequency radiation hydrodynamics AMR simulations that are carefully calibrated to reproduce a realistic reionization history and galaxy properties at $z \geq 6$. In our simulations, haloes with mass $10^9{\rm M_{\odot}}h^{-1}, Comment: 13 pages, 8 figures, Accepted to MNRAS

We investigate the ultraviolet (UV) spectral properties of faint Lyman-$\alpha$ emitters (LAEs) in the redshift range 2.93. We compute various averaged spectra of LAEs sub-sampled on the basis of their observational (e.g., Ly$\alpha$ strength, UV magnitude and spectral slope) and physical (e.g., stellar mass and star-formation rate) properties. We search for UV spectral features other than Ly$\alpha$, such as higher-ionization nebular emission lines and absorption features. We successfully observe the OIII]1666 and CIII]909 collisionally excited emission lines and the HeII1640 recombination feature, as well as the resonant CIV1550 doublet either in emission or P-Cygni. We compare the observed spectral properties of the different mean spectra and find the emission lines to vary with the observational and physical properties of the LAEs. In particular, the mean spectra of LAEs with larger Ly$\alpha$ equivalent widths, fainter UV magnitudes, bluer UV spectral slopes and lower stellar masses show the strongest nebular emission. The line ratios of these lines are similar to those measured in the spectra of local metal-poor galaxies, while their equivalent widths are weaker compared to the handful of extreme values detected in individual spectra of z>2 galaxies. This suggests that weak UV features are likely ubiquitous in high z, low-mass and faint LAEs. We publicly release the stacked spectra as they can serve as empirical templates for the design of future observations, such as those with the James Webb Space Telescope and the Extremely Large Telescope., Comment: 26 pages, 18 figures. Accepted for publication in A&A

We present a new, on-the-fly photon flux and absorption tracer algorithm designed to directly measure the contribution of different source populations to the metagalactic UV background and to the ionisation fraction of gas in the Universe. We use a suite of multifrequency radiation hydrodynamics simulations that are carefully calibrated to reproduce a realistic reionization history and galaxy properties at $z \ge 6$, to disentangle the contribution of photons emitted by different mass haloes and by stars with different metallicities and ages to the UV background during reionization. While at very early cosmic times low mass, metal poor haloes provide most of the LyC photons, their contribution decreases steadily with time. At $z = 6$ it is the photons emitted by massive systems (${\rm M_{halo}}/{\rm M_\odot} > 10^{10} \, {\rm h ^{-1}}$) and by the metal enriched stars ($10^{-3} < Z/Z_{\rm \odot} < 10^{-1.5}$) that provide the largest contribution to the ionising UV background. We demonstrate that there are large variations in the escape fraction depending on the source, with the escape fraction being highest ($\sim 45-60\%$) for photons emitted by the oldest stars that penetrate into the IGM via low opacity channels carved by the ionising photons and supernova from younger stars. Before HII regions begin to overlap, the photoionisation rate strongly fluctuates between different, isolated HII bubbles, depending on the embedded ionising source, which we suggest may result in spatial variations in the properties of dwarf galaxies., 19 pages, 14 figures, Accepted MNRAS

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 $10^{10}\,M_\odot$ halo, we show that the momentum transferred from resonantly scattered Lyman-$\alpha$ (Lya) photons is an important source of stellar feedback which can shape the evolution of galaxies. We find that Lya 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 Lya 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 $\sim 5$, compared to the model without the early feedback. More importantly, we find that galactic outflows become weaker in the presence of strong Lya 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 $\sim5-10$ near the mid plane, while it is reduced to $\sim1$ 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 Lya feedback., Comment: 19 pages, 16 figures, Accepted for publication in MNRAS

We present the SPHINX suite of cosmological adaptive mesh refinement simulations, the first radiation-hydrodynamical simulations to simultaneously capture large-scale reionization and the escape of ionizing radiation from thousands of resolved galaxies. Our $5$ and $10$ co-moving Mpc volumes resolve haloes down to the atomic cooling limit and model the inter-stellar medium with better than $\approx10$ pc resolution. The project has numerous goals in improving our understanding of reionization and making predictions for future observations. In this first paper we study how the inclusion of binary stars in computing stellar luminosities impacts reionization, compared to a model that includes only single stars. Owing to the suppression of galaxy growth via strong feedback, our galaxies are in good agreement with observational estimates of the galaxy luminosity function. We find that binaries have a significant impact on the timing of reionization: with binaries, our boxes are $99.9$ percent ionized by volume at $z\approx 7$, while without them our volumes fail to reionize by $z=6$. These results are robust to changes in volume size, resolution, and feedback efficiency. The escape of ionizing radiation from individual galaxies varies strongly and frequently. On average, binaries lead to escape fractions of $\approx 7-10$ percent, about $3.5$ times higher than with single stars only. The higher escape fraction is a result of a shallower decline in ionizing luminosity with age, and is the primary reason for earlier reionization, although the higher integrated luminosity with binaries also plays a sub-dominant role., Comment: 18 pages, 15 figures. Accepted for publication in MNRAS

We analyse the properties of circumgalactic gas around simulated galaxies in the redshift range z >= 3, utilising a new sample of cosmological zoom simulations. These simulations are intended to be representative of the observed samples of Lyman-alpha emitters recently obtained with the MUSE instrument (halo masses ~10^10-10^11 solar masses). We show that supernova feedback has a significant impact on both the inflowing and outflowing circumgalactic medium by driving outflows, reducing diffuse inflow rates, and by increasing the neutral fraction of inflowing gas. By temporally stacking simulation outputs we find that significant net mass exchange occurs between inflowing and outflowing phases: none of the phases are mass-conserving. In particular, we find that the mass in neutral outflowing hydrogen declines exponentially with radius as gas flows outwards from the halo centre. This is likely caused by a combination of both fountain-like cycling processes and gradual photo/collisional ionization of outflowing gas. Our simulations do not predict the presence of fast-moving neutral outflows in the CGM. Neutral outflows instead move with modest radial velocities (~ 50 kms^-1), and the majority of the kinetic energy is associated with tangential rather than radial motion., 24 pages, 17 figures, submitted to MNRAS after responding to a first referee report

The physical origin of the near-ultraviolet Mg II emission remains an underexplored domain, unlike more typical emission lines that are detected in the spectra of star-forming galaxies. We explore the nebular and physical properties of a sample of 381 galaxies between 0.70 < z < 2.34 drawn from the MUSE Hubble Ultra Deep Survey. The spectra of these galaxies show a wide variety of profiles of the Mg II λλ2796, 2803 resonant doublet, from absorption to emission. We present a study on the main drivers for the detection of Mg II emission in galaxy spectra. By exploiting photoionization models, we verified that the emission-line ratios observed in galaxies with Mg II in emission are consistent with nebular emission from HII regions. From a simultaneous analysis of MUSE spectra and ancillary Hubble Space Telescope information through spectral energy distribution fitting, we find that galaxies with Mg II in emission have lower stellar masses, smaller sizes, bluer spectral slopes, and lower optical depth than those with absorption. This leads us to suggest that Mg II emission is a potential tracer of physical conditions that are not merely related to those of the ionized gas. We show that these differences in Mg II emission and absorption can be explained in terms of a higher dust and neutral gas content in the interstellar medium (ISM) of galaxies showing Mg II in absorption, which confirms the extreme sensitivity of Mg II to the presence of the neutral ISM. We conclude with an analogy between the Mg II doublet and the Ly α line that lies in their resonant nature. Further investigations with current and future facilities, including the James Webb Space Telescope, are promising because the detection of Mg II emission and its potential connection with Lyα could provide new insights into the ISM content in the early Universe., Astronomy & Astrophysics, 617, ISSN:0004-6361, ISSN:1432-0746

We report the detection of extended Lyα haloes around 145 individual star-forming galaxies at redshifts 3 ≤ z ≤ 6 in the Hubble Ultra Deep Field observed with the Multi-Unit Spectroscopic Explorer (MUSE) at ESO-VLT. Our sample consists of continuum-faint (− 15 ≥ MUV ≥ −22) Lyα emitters (LAEs). Using a 2D, two-component (continuum-like and halo) decomposition of Lyα emission assuming circular exponential distributions, we measure scale lengths and luminosities of Lyα haloes. We find that 80% of our objects having reliable Lyα halo measurements show Lyα emission that is significantly more extended than the UV continuum detected by HST (by a factor ≈4 to >20). The median exponential scale length of the Lyα haloes in our sample is ≈4.5 kpc with a few haloes exceeding 10 kpc. By comparing the maximal detected extent of the Lyα emission with the predicted dark matter halo virial radii of simulated galaxies, we show that the detected Lyα emission of our selected sample of Lyα emitters probes a significant portion of the cold circum-galactic medium of these galaxies (>50% in average). This result therefore shows that there must be significant HI reservoirs in the circum-galactic medium and reinforces the idea that Lyα haloes are ubiquitous around high-redshift Lyα emitting galaxies. Our characterization of the Lyα haloes indicates that the majority of the Lyα flux comes from the halo (≈65%) and that their scale lengths seem to be linked to the UV properties of the galaxies (sizes and magnitudes). We do not observe a significant Lyα halo size evolution with redshift, although our sample for z> 5 is very small. We also explore the diversity of the Lyα line profiles in our sample and we find that the Lyα lines cover a large range of full width at half maximum (FWHM) from 118 to 512 km s-1. While the FWHM does not seem to be correlated to the Lyα scale length, most compact Lyα haloes and those that are not detected with high significance tend to have narrower Lyα profiles (-1). Finally, we investigate the origin of the extended Lyα emission but we conclude that our data do not allow us to disentangle the possible processes, i.e. scattering from star-forming regions, fluorescence, cooling radiation from cold gas accretion, and emission from satellite galaxies.

We present a clustering comparison of 12 galaxy formation models [including semi-analytic models (SAMs) and halo occupation distribution (HOD) models] all run on halo catalogues and merger trees extracted from a single Λ cold dark matter N-body simulation. We compare the results of the measurements of the mean halo occupation numbers, the radial distribution of galaxies in haloes and the two-point correlation functions (2PCF). We also study the implications of the different treatments of orphan (galaxies not assigned to any dark matter subhalo) and non-orphan galaxies in these measurements. Our main result is that the galaxy formation models generally agree in their clustering predictions but they disagree significantly between HOD and SAMs for the orphan satellites. Although there is a very good agreement between the models on the 2PCF of central galaxies, the scatter between the models when orphan satellites are included can be larger than a factor of 2 for scales smaller than 1 h(−1) Mpc. We also show that galaxy formation models that do not include orphan satellite galaxies have a significantly lower 2PCF on small scales, consistent with previous studies. Finally, we show that the 2PCF of orphan satellites is remarkably different between SAMs and HOD models. Orphan satellites in SAMs present a higher clustering than in HOD models because they tend to occupy more massive haloes. We conclude that orphan satellites have an important role on galaxy clustering and they are the main cause of the differences in the clustering between HOD models and SAMs., Monthly Notices of the Royal Astronomical Society, 469 (1), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

We present the deepest study to date of the Lya luminosity function (LF) in a blank field using blind integral field spectroscopy from MUSE. We constructed a sample of 604 Lya emitters (LAEs) across the redshift range 2.91 < z < 6.64 using automatic detection software in the Hubble Ultra Deep Field. We calculate accurate total Lya fluxes capturing low surface brightness extended Lya emission now known to be a generic property of high-redshift star-forming galaxies. We simulated realistic extended LAEs to characterise the selection function of our samples, and performed flux-recovery experiments to test and correct for bias in our determination of total Lya fluxes. We find an accurate completeness correction accounting for extended emission reveals a very steep faint-end slope of the LF, alpha, down to luminosities of log10 L erg s^-1< 41.5, applying both the 1/Vmax and maximum likelihood estimators. Splitting the sample into three broad redshift bins, we see the faint-end slope increasing from -2.03+1.42-inf at z ~ 3.44 to -2.86+0.76-inf at z ~ 5.48, however no strong evolution is seen between the 68% confidence regions in L*-alpha parameter space. Using the Lya line flux as a proxy for star formation activity, and integrating the observed LFs, we find that LAEs' contribution to the cosmic SFRD rises with redshift until it is comparable to that from continuum-selected samples by z ~ 6. This implies that LAEs may contribute more to the star-formation activity of the early Universe than previously thought - any additional interglactic medium correction would act to further boost the Lya luminosities. Finally, assuming fiducial values for the escape of Lya and LyC radiation, and the clumpiness of the IGM, we integrated the maximum likelihood LF at 5.00 < z < 6.64 and find we require only a small extrapolation beyond the data (< 1 dex in L) for LAEs alone to maintain an ionised IGM at z ~ 6., Comment: 16 pages, 10 figures. Accepted for publication in A&A special issue: "The MUSE Hubble Ultra Deep Field Survey" (MUSE UDF Series Paper VI)

Low mass galaxies are thought to provide the bulk of the ionizing radiation necessary to reionize the Universe. The amount of photons escaping the galaxies is poorly constrained theoretically, and difficult to measure observationally. Yet it is an essential parameter of reionization models. We study in detail how ionizing radiation can leak from high redshift galaxies. For this purpose, we use a series of high resolution radiation hydrodynamics simulations, zooming on three dwarf galaxies in a cosmological context. We find that the energy and momentum input from the supernova explosions has a pivotal role in regulating the escape fraction, by disrupting dense star forming clumps, and clearing sight lines in the halo. In the absence of supernovae, photons are absorbed very locally, within the birth clouds of massive stars. We follow the time evolution of the escape fraction, and find that it can vary by more than six orders of magnitude. This explains the large scatter in the value of the escape fraction found by previous studies. This fast variability also impacts the observability of the sources of reionization: a survey even as deep as $M_{\rm UV} = -14$ would miss about half of the underlying population of Lyman-continuum emitters., Comment: 18 pages, 17 figures. Accepted for publication in MNRAS

We present the first estimate of the Ly{\alpha} luminosity function using blind spectroscopy from the Multi Unit Spectroscopic Explorer, MUSE, in the Hubble Deep Field South. Using automatic source-detection software, we assemble a homogeneously-detected sample of 59 Ly{\alpha} emitters covering a flux range of -18.0 < log10 (F) < -16.3 (erg s^-1 cm^-2), corresponding to luminosities of 41.4 < log10 (L) < 42.8 (erg s^-1). As recent studies have shown, Ly{\alpha} fluxes can be underestimated by a factor of two or more via traditional methods, and so we undertake a careful assessment of each object's Ly{\alpha} flux using a curve-of-growth analysis to account for extended emission. We describe our self-consistent method for determining the completeness of the sample, and present an estimate of the global Ly{\alpha} luminosity function between redshifts 2.91 < z < 6.64 using the 1/Vmax estimator. We find the luminosity function is higher than many number densities reported in the literature by a factor of 2 - 3, although our result is consistent at the 1{\sigma} level with most of these studies. Our observed luminosity function is also in good agreement with predictions from semi-analytic models, and shows no evidence for strong evolution between the high- and low-redshift halves of the data. We demonstrate that one's approach to Ly{\alpha} flux estimation does alter the observed luminosity function, and caution that accurate flux assessments will be crucial in measurements of the faint end slope. This is a pilot study for the Ly{\alpha} luminosity function in the MUSE deep-fields, to be built on with data from the Hubble Ultra Deep Field which will increase the size of our sample by almost a factor of 10., Comment: 13 pages, 8 figures, submitted to MNRAS

We use the Ly-$\alpha$ 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-$\alpha$ forest, and compare to cross-correlations measured for quasars and damped Ly-$\alpha$ 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., Comment: Accepted by MNRAS, 23 pages, 20 figures. For a short video summarizing this paper, please visit: https://youtu.be/9ghGNtF16JA

We present the Mock Map Facility, a powerful tool to generate mock catalogues or images from semi-analytically post-processed snapshots of cosmological N-body simulations. The paper describes in detail an efficient technique to create such mocks from the GALICS semi-analytic model, providing the reader with an accurate quantification of the artifacts it introduces at every step. We show that replication effects introduce a negative bias on the clustering signal -- typically peaking at less than 10 percent around the correlation length. We also thoroughly discuss how the clustering signal is affected by finite volume effects, and show that it vanishes at scales larger than about a tenth of the simulation box size. For the purpose of analysing our method, we show that number counts and redshift distributions obtained with GALICS and MOMAF compare well to K-band observations and to the 2dFGRS. Given finite volume effects, we also show that the model can reproduce the APM angular correlation function. The MOMAF results discussed here are made publicly available to the astronomical community through a public database. Moreover, a user-friendly Web interface (http://galics.iap.fr) allows any user to recover her/his own favourite galaxy samples through simple SQL queries. The flexibility of this tool should permit a variety of uses ranging from extensive comparisons between real observations and those predicted by hierarchical models of galaxy formation, to the preparation of observing strategies for deep surveys and tests of data processing pipelines., 19 pages, 15 Figs, significantly modified version now accepted for publication in MNRAS. High-resolution version available at http://galics.cosmologie.fr/papers/momaf.ps.gz

Direct Lyα imaging of intergalactic gas at $z\sim 2$ has recently revealed giant cosmological structures around quasars, e.g., the Slug Nebula. Despite their high luminosity, the detection rate of such systems in narrow-band and spectroscopic surveys is less than 10%, possibly encoding crucial information on the distribution of gas around quasars and the quasar emission properties. In this study, we use the MUSE integral-field instrument to perform a blind survey for giant $\mathrm{Ly}\alpha $ nebulae around 17 bright radio-quiet quasars at $3\lt z\lt 4$ that does not suffer from most of the limitations of previous surveys. After data reduction and analysis performed with specifically developed tools, we found that each quasar is surrounded by giant $\mathrm{Ly}\alpha $ nebulae with projected sizes larger than 100 physical kiloparsecs and, in some cases, extending up to 320 kpc. The circularly averaged surface brightness profiles of the nebulae appear to be very similar to each other despite their different morphologies and are consistent with power laws with slopes $\approx -1.8$. The similarity between the properties of all these nebulae and the Slug Nebula suggests a similar origin for all systems and that a large fraction of gas around bright quasars could be in a relatively “cold” (T ∼ 10(4 )K) and dense phase. In addition, our results imply that such gas is ubiquitous within at least 50 kpc from bright quasars at $3\lt z\lt 4$ independently of the quasar emission opening angle, or extending up to 200 kpc for quasar isotropic emission.

Trabajo presentado en SPIE Astronomical Telescopes, celebrado en San Diego (California), del 26 de junio al 1 de julio de 2016, HARMONI is the E-ELT's first light visible and near-infrared integral field spectrograph. It will provide four different spatial scales, ranging from coarse spaxels of 60 × 30 mas best suited for seeing limited observations, to 4 mas spaxels that Nyquist sample the diffraction limited point spread function of the E-ELT at near-infrared wavelengths. Each spaxel scale may be combined with eleven spectral settings, that provide a range of spectral resolving powers (R 3500, 7500 and 20000) and instantaneous wavelength coverage spanning the 0.5 - 2.4 ¿m wavelength range of the instrument. In autumn 2015, the HARMONI project started the Preliminary Design Phase, following signature of the contract to design, build, test and commission the instrument, signed between the European Southern Observatory and the UK Science and Technology Facilities Council. Crucially, the contract also includes the preliminary design of the HARMONI Laser Tomographic Adaptive Optics system. The instrument's technical specifications were finalized in the period leading up to contract signature. In this paper, we report on the first activity carried out during preliminary design, defining the baseline architecture for the system, and the trade-off studies leading up to the choice of baseline.

Spatially resolved studies of high redshift galaxies, an essential insight into galaxy formation processes, have been mostly limited to stacking or unusually bright objects. We present here the study of a typical (L$^{*}$, M$_\star$ = 6 $\times 10^9$ $M_\odot$) young lensed galaxy at $z=3.5$, observed with MUSE, for which we obtain 2D resolved spatial information of Ly$\alpha$ and, for the first time, of CIII] emission. The exceptional signal-to-noise of the data reveals UV emission and absorption lines rarely seen at these redshifts, allowing us to derive important physical properties (T$_e\sim$15600 K, n$_e\sim$300 cm$^{-3}$, covering fraction f$_c\sim0.4$) using multiple diagnostics. Inferred stellar and gas-phase metallicities point towards a low metallicity object (Z$_{\mathrm{stellar}}$ = $\sim$ 0.07 Z$_\odot$ and Z$_{\mathrm{ISM}}$ $, Comment: 19 pages, 15 figures, accepted in MNRAS

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 f esc 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 new approach to constrain galaxy physical parameters from the combined interpretation of stellar and nebular emission in wide ranges of observations. This approach relies on the Bayesian analysis of any type of galaxy spectral energy distribution using a comprehensive library of synthetic spectra assembled using state-of-the-art models of star formation and chemical enrichment histories, stellar population synthesis, nebular emission and attenuation by dust. We focus on the constraints set by 5-band photometry and low- and medium-resolution spectroscopy at optical rest wavelengths on a set of physical parameters characterizing the stars and the interstellar medium. Since these parameters cannot be known a priori for any galaxy sample, we assess the accuracy to which they can be retrieved by simulating `pseudo-observations' using models with known parameters. Assuming that these models are good approximations of true galaxies, we find that the combined analysis of stellar and nebular emission in low-resolution galaxy spectra provides valuable constraints on all physical parameters. At higher resolution, the analysis of the combined stellar and nebular emission in 12,660 SDSS star-forming galaxies using our approach yields likelihood distributions of stellar mass, gas-phase oxygen abundance, optical depth of the dust and specific star formation rate similar to those obtained in previous separate analyses of the stellar and nebular emission at the original (twice higher) SDSS spectral resolution. We show that the constraints derived on galaxy physical parameters from these different types of observations depend sensitively on signal-to-noise ratio. Our approach can be extended to the analysis of any type of observation across the wavelength range covered by spectral evolution models. [abridged]

The VLT Multi Unit Spectroscopic Explorer (MUSE) integral-field spectrograph can detect Ly\alpha{} 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\alpha{} 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\alpha{} radiation transfer models in gas outflows. We expect $\approx$ 1500 bright LAEs ($F_{Ly\alpha}$ $\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\alpha}$ $\gtrsim$ $4 \times 10^{-19}$ erg s$^{-1}$ cm$^{-2}$) will be found. Our results suggest that faint Ly\alpha{} sources contribute significantly to the cosmic Ly\alpha{} 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., Comment: 18 pages, 13 figures, accepted for publication in MNRAS

We present numerical simulations of a 15 solar mass star in a suite of idealised environments in order to quantify the amount of energy transmitted to the interstellar medium (ISM). We include models of stellar winds, UV photoionisation and the subsequent supernova based on theoretical models and observations of stellar evolution. The system is simulated in 3D using RAMSES-RT, an Adaptive Mesh Refinement Radiation Hydrodynamics code. We find that stellar winds have a negligible impact on the system owing to their relatively low luminosity compared to the other processes. The main impact of photoionisation is to reduce the density of the medium into which the supernova explodes, reducing the rate of radiative cooling of the subsequent supernova. Finally, we present a grid of models quantifying the energy and momentum of the system that can be used to motivate simulations of feedback in the ISM unable to fully resolve the processes discussed in this work., 19 pages, 12 figures, accepted by MNRAS

We observed the Hubble Deep Field South with the new panoramic integral field spectrograph MUSE that we built and just commissioned at the VLT. The data cube resulting from 27 hours of integration covers one arcmin^2 field of view at an unprecedented depth with a 1 sigma emission line surface brightness limit of 1x$10^{-19}$ erg/s/cm$^2$/arcsec$^2$ and contains ~90,000 spectra. We present the combined and calibrated data cube, and we perform a first-pass analysis of the sources detected in the HDF-S imaging. We measured the redshifts of 189 sources up to a magnitude F814W = 29.5, increasing by more than an order of magnitude the number of known spectroscopic redshifts in this field. We also discovered 26 Lya emitting galaxies which are not detected in the HST WFPC2 deep broad band images. The intermediate spectral resolution of 2.3{\AA} allows us to separate resolved asymmetric Lya emitters, [O II] emitters, and C III] emitters and the large instantaneous wavelength range of 4500{\AA} helps to identify single emission lines. We also show how the three dimensional information of MUSE helps to resolve sources which are confused at ground-based image quality. Overall, secure identifications are provided for 83% of the 227 emission line sources detected in the MUSE data cube and for 32% of the 586 sources identified in the HST catalog of Casertano et al 2000. The overall redshift distribution is fairly flat to z=6.3, with a reduction between z=1.5 to 2.9, in the well-known redshift desert. The field of view of MUSE also allowed us to detect 17 groups within the field. We checked that the number counts of [O II] and Ly-a emitters are roughly consistent with predictions from the literature. Using two examples we demonstrate that MUSE is able to provide exquisite spatially resolved spectroscopic information on intermediate redshift galaxies present in the field., Comment: 21 pages, 20 figures. Note that the catalogues, data cube and associated spectra will be released upon acceptance of the paper

Number Counts of galaxies in two GALEX bands (FUV: 1530A and NUV: 2310A, both in AB magnitudes) are reported. They provide for the first time in the literature homogeneously calibrated number counts of UV galaxies covering continuously a very wide range of UV magnitude (14 -- 23.8). Both the FUV and NUV counts are inconsistent with a non-evolution model, while they are in good agreement with evolution models (essentially luminosity evolution) derived from the high-z UV luminosity functions of Arnouts et al. (2004). It is found that the contribution from galaxies detected by GALEX to the UV background is 0.68+-0.10 nW m-2 sr-1 at 1530A and 0.99+-0.15 nW m-2 sr-1 at 2310A. These are 66+-9% and 44+-6% of the total contributions of galaxies to the the UV background at 1530A, respectively, as estimated using the evolution models. ...

Non-resonant FeII* 2365, 2396, 2612, 2626 emission can potentially trace galactic winds in emission and provide useful constraints to wind models. From the 3'x3' mosaic of the Hubble Ultra Deep Field (UDF) obtained with the VLT/MUSE integral field spectrograph, we identify a statistical sample of 40 FeII* emitters and 50 MgII 2796, 2803 emitters from a sample of 271 [OII] 3726, 3729 emitters with reliable redshifts from z = 0.85 - 1.5 down to 2E-18 (3 sigma) ergs/s/cm^2 (for [OII]), covering the stellar mass range 10^8 - 10^11 Msun. The FeII* and MgII emitters follow the galaxy main sequence, but with a clear dichotomy. Galaxies with masses below 10^9 Msun and star formation rates (SFRs) of 10 Msun/year have FeII* emission without accompanying MgII emission. Between these two regimes, galaxies have both MgII and FeII* emission, typically with MgII P-Cygni profiles. Indeed, the MgII profile shows a progression along the main sequence from pure emission to P-Cygni profiles to strong absorption, due to resonant trapping. Combining the deep MUSE data with HST ancillary information, we find that galaxies with pure MgII emission profiles have lower star formation rate surface densities than those with either MgII P-Cygni profiles or FeII* emission. These spectral signatures produced through continuum scattering and fluorescence, MgII P-Cygni profiles and FeII* emission, are better candidates for tracing galactic outflows than pure MgII emission, which may originate from HII regions. We compare the absorption and emission rest-frame equivalent widths for pairs of FeII transitions to predictions from outflow models and find that the observations consistently have less total re-emission than absorption, suggesting either dust extinction or non-isotropic outflow geometries., 16 pages, 12 figures, 3 tables, accepted for publication in A&A (14 October 2017)

We present rest-frame Lyα equivalent widths (E W 0 ) of 417 Lyα emitters (LAEs) detected with Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) at 2.9 Ultra Deep Field. Based on the deep MUSE spectroscopy and ancillary Hubble Space Telescope (HST) photometry data, we carefully measured E W 0 values taking into account extended Lyα emission and UV continuum slopes (β ). Our LAEs reach unprecedented depths, both in Lyα luminosities and UV absolute magnitudes, from log (L Lyα /erg s-1 ) ~ 41.0 to 43.0 and from M UV ~ −16 to −21 (0.01−1.0 L * z =3 ). The E W 0 values span the range of ~ 5 to 240 A or larger, and their distribution can be well fitted by an exponential law N = N 0 exp(−E W 0 /w 0 ). Owing to the high dynamic range in M UV , we find that the scale factor, w 0 , depends on M UV in the sense that including fainter M UV objects increases w 0 , i.e., the Ando effect. The results indicate that selection functions affect the E W 0 scale factor. Taking these effects into account, we find that our w 0 values are consistent with those in the literature within 1σ uncertainties at 2.9 UV and L Lyα . Interestingly, we find 12 objects with E W 0 > 200 A above 1σ uncertainties. Two of these 12 LAEs show signatures of merger or AGN activity: the weak Civλ 1549 emission line. For the remaining 10 very large E W 0 LAEs, we find that the E W 0 values can be reproduced by young stellar ages ( ⊙ ). Otherwise, at least part of the Lyα emission in these LAEs needs to arise from anisotropic radiative transfer effects, fluorescence by hidden AGN or quasi-stellar object activity, or gravitational cooling.

International audience; Semi-analytical models (SAMs) are currently the best way to understand the formation of galaxies within the cosmic dark-matter structures. They are able to give a statistical view of the variety of the evolutionary histories of galaxies in terms of star formation and stellar mass assembly. While they reproduce the local stellar mass functions, correlation functions, and luminosity functions fairly well, they fail to match observations at high redshift (z \textgreater= 3) in most cases, particularly in the low-mass range. The inconsistency between models and CDM observations indicates that the history of gas accretion in galaxies, within their host dark-matter halo, and the transformation of gas into stars, are not followed well. We briefly present a new version of the GalICS semi-analytical model. With this new model, we explore the impact of classical mechanisms, such as supernova feedback or photoionization, on the evolution of the stellar mass assembly and the star formation rate. Even with strong efficiency, these two processes cannot explain the observed stellar mass function and star formation rate distribution or the stellar mass versus dark matter halo mass relation. We thus introduce an ad hoc modification of the standard paradigm, based on the presence of a no-star-forming gas component, and a concentration of the star-forming gas in galaxy discs. The main idea behind the existence of the no-star-forming gas reservoir is that only a fraction of the total gas mass in a galaxy is available to form stars. The reservoir generates a delay between the accretion of the gas and the star formation process. This new model is in much better agreement with the observations of the stellar mass function in the low-mass range than the previous models and agrees quite well with a large set of observations, including the redshift evolution of the specific star formation rate. However, it predicts a large amount of no-star-forming baryonic gas, potentially larger than observed, even if its nature has still to be examined in the context of the missing baryon problem. Outputs from all models are available at the CDS.

[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., Accepted for publication in ApJ (minor corrections from the previous version). Catalogs of mock BOSS spectra and relevant data can be found at: http://www2.iap.fr/users/peirani/lymas/lymas.htm

We present a new implementation of radiation hydrodynamics (RHD) in the adaptive mesh refinement (AMR) code RAMSES. The multi-group radiative transfer (RT) is performed on the AMR grid with a first-order Godunov method using the M1 closure for the Eddington tensor, and is coupled to the hydrodynamics via non-equilibrium thermochemistry of hydrogen and helium. This moment-based approach has the large advantage that the computational cost is independent of the number of radiative sources - it can even deal with continuous regions of emission such as bound-free emission from gas. As it is built directly into RAMSES, the RT takes natural advantage of the refinement and parallelization strategies already in place. Since we use an explicit advection solver for the radiative transport, the time step is restricted by the speed of light - a severe limitation that can be alleviated using the so--called "reduced speed of light" approximation. We propose a rigorous framework to assess the validity of this approximation in various conditions encountered in cosmology and galaxy formation. We finally perform with our newly developed code a complete suite of RHD tests, comparing our results to other RHD codes. The tests demonstrate that our code performs very well and is ideally suited for exploring the effect of radiation on current scenarios of structure and galaxy formation., Comment: 48 pages, 37 figures, accepted for publication in MNRAS

Simulations of galaxy evolution aim to capture our current understanding as well as to make predictions for testing by future experiments. Simulations and observations are often compared in an indirect fashion: physical quantities are estimated from the data and compared to models. However, many applications can benefit from a more direct approach, where the observing process is also simulated and the models are seen fully from the observer's perspective. To facilitate this, we have developed the Millennium Run Observatory (MRObs), a theoretical virtual observatory which uses virtual telescopes to `observe' semi-analytic galaxy formation models based on the suite of Millennium Run dark matter simulations. The MRObs produces data that can be processed and analyzed using the standard software packages developed for real observations. At present, we produce images in forty filters from the rest-frame UV to IR for two stellar population synthesis models, three different models of IGM absorption, and two cosmologies (WMAP1/7). Galaxy distributions for a large number of mock lightcones can be `observed' using models of major ground- and space-based telescopes. The data include lightcone catalogues linked to structural properties of galaxies, pre-observation model images, mock telescope images, and Source Extractor products that can all be traced back to the higher level dark matter, semi-analytic galaxy, and lightcone catalogues available in the Millennium database. Here, we describe our methods and announce a first public release of simulated surveys (e.g., SDSS, CFHT-LS, GOODS, GOODS/ERS, CANDELS, and HUDF). The MRObs browser, an online tool, further facilitates exploration of the simulated data. We demonstrate the benefits of a direct approach through a number of example applications (galaxy number counts in CANDELS, clusters, morphologies, and dropout selections)., MNRAS, in press. Millennium Run Observatory data products, online tools, and more available through http://galformod.mpa-garching.mpg.de/mrobs/

We develop a Principal Component Analysis aimed at classifying a sub-set of 27,350 spectra of galaxies in the range 0.4 < z < 1.0 collected by the VIMOS Public Extragalactic Redshift Survey (VIPERS). We apply an iterative algorithm to simultaneously repair parts of spectra affected by noise and/or sky residuals, and reconstruct gaps due to rest-frame transformation, and obtain a set of orthogonal spectral templates that span the diversity of galaxy types. By taking the three most significant components, we find that we can describe the whole sample without contamination from noise. We produce a catalogue of eigen-coefficients and template spectra that will be part of future VIPERS data releases. Our templates effectively condense the spectral information into two coefficients that can be related to the age and star formation rate of the galaxies. We examine the spectrophotometric types in this space and identify early, intermediate, late and starburst galaxies., 15 pages, 20 images, accepted for publication in MNRAS: MN-12-1739-MJ.R1

[Abridged] We perform on galaxy mock catalogues the same colour-density analysis made by Cucciati et al. (2006) on a 5 Mpc/h scale using the VVDS-Deep survey, and compare the results from mocks with observed data. We use mocks with the same flux limits (I=24) as the VVDS (CMOCKS), built using the semi- analytic model by De Lucia & Blaizot (2007) applied to the Millennium Simulation. From CMOCKS, we extracted samples of galaxies mimicking the VVDS observational strategy (OMOCKS). We computed the B-band Luminosity Function LF and the colour-density relation (CDR) in the mocks. We find that the LF in mocks roughly agrees with the observed LF, but at z, Comment: 15 pages, 12 figures, accepted for publication in A&A

Lyman-$\alpha$ nebulae are usually found in massive environments at high redshift ($z > 2$). The origin of their Lyman-$\alpha$ (Lya) emission remains debated. Recent polarimetric observations showed that at least some Lya sources are polarized. This is often interpreted as a proof that the photons are centrally produced, and opposed to the scenario in which the Lya emission is the cooling radiation emitted by gas heated during the accretion onto the halo. We suggest that this scenario is not incompatible with the polarimetric observations. In order to test this idea, we post-process a radiative hydrodynamics simulation of a blob with the MCLya Monte Carlo transfer code. We compute radial profiles for the surface brightness and the degree of polarization and compare them to existing observations. We find that both are consistent with a significant contribution of the extragalactic gas to the Lya emission. Most of the photons are centrally emitted and scattered inside the filament afterwards, producing the observed high level of polarization. We argue that the contribution of the extragalactic gas to the Lya emission does not prevent polarization to arise. On the contrary, we find that pure galactic emission causes the polarization profile to be too steep to be consistent with observations., Comment: 14 pages, 9 figures (+4 in appendices). Accepted for publication in Astronomy & Astrophysics

[abridged] Aims. The aim of this paper is to assess the impact of the interstellar medium (ISM) physics on Lyman-alpha (Lya) radiation transfer and to quantify how galaxy orientation with respect to the line of sight alters observational signatures. Methods. We compare the results of Lya radiation transfer calculations through the ISM of a couple of idealized galaxy simulations with different ISM models. Results. First, the small-scale structuration of the ISM plays a determinant role in shaping a galaxys Lya properties.The artificially warm, and hence smooth, ISM of G1 yields an escape fraction of 50 percent at the Lya line center, and produces symmetrical double-peak profiles. On the contrary, in G2, most young stars are embedded in thick star-forming clouds, and the result is a 10 times lower escape fraction. G2 also displays a stronger outflowing velocity field, which favors the escape of red-shifted photons, resulting in an asymmetric Lya line. Second, the Lya properties of G2 strongly depend on the inclination at which it is observed: From edge-on to face-on, the line goes from a double-peak profile with an equivalent width of -5 Angstrom to a 15 times more luminous red-shifted asymmetric line with EW 90 Angstrom. Conclusions. Lya radiation transfer calculations can only lead to realistic properties in simulations where galaxies are resolved into giant molecular clouds, putting these calculations out of reach of current large scale cosmological simulations. Finally, we find inclination effects to be much stronger for Lya photons than for continuum radiation. This could potentially introduce severe biases in the selection function of narrow-band Lya emitter surveys, which could indeed miss a significant fraction of the high-z galaxy population., Comment: 13 pages, 11 figures, accepted for publication in Astronomy&Astrophysics

[Abridged] We present a homogeneous and complete catalogue of optical groups identified in the purely flux limited (17.5, Submitted to A&A, revised version after referee comments, Table 5 added