Search

Galactic outflows are believed to play a critical role in the evolution of galaxies by regulating their mass build-up and star formation [1]. These outflows tend to form bipolar shapes above and below the galactic planes and extend well into the circumgalactic medium (CGM), up to tens of kpc perpendicular to the galaxy. They have been directly observed in the local Universe, e.g., around the Milky Way and M82 [2, 3]. At higher redshifts, cosmological simulations of galaxy formation predict an increase in the frequency and efficiency of galactic outflows due to the increasing star formation activity [4]. Outflows are responsible for removing potential fuel for star formation from the galaxy, while at the same enriching the CGM and the intergalactic medium. These feedback processes, although incorporated as key elements of cosmological simulations, are still poorly constrained on CGM scales. Here we present an ultra-deep MUSE image of the mean Mg II emission surrounding a sample of galaxies at z ≈ 1 that strongly suggests the presence of outflowing gas on physical scales of more than 10 kpc. We find a strong dependence of the detected signal on the inclination of the central galaxy, with edge-on galaxies clearly showing enhanced Mg II emission along the minor axis, while face-on galaxies display much weaker and more isotropic emission. We interpret these findings as supporting the idea that outflows typically have a bipolar cone geometry perpendicular to the galactic disk. We demonstrate that the signal is not dominated by a few outliers. After dividing the galaxy sample in subsamples by mass, the bipolar emission is only detected in galaxies with stellar mass M* >∼ 109.5 M⊙.

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

ALMA observations of z > 6 galaxies hav e rev ealed abnormally high [O III ] 88 μm /[C II ] 158 μm ratios and [C II ] 158 μm deficits compared to local galaxies. The origin of this behaviour is unknown. Numerous solutions have been proposed including differences in C and O abundance ratios, observational bias, and differences in ISM properties, including ionization parameter, gas density, or photodissociation region (PDR) covering fraction. In order to elucidate the underlying physics that drives this high- redshift phenomenon, we employ SPHINX 20 , a state-of-the-art, cosmological radiation-hydrodynamics simulation, that resolves detailed ISM properties of thousands of galaxies in the epoch of reionization which has been post-processed with CLOUDY to predict emission lines. We find that the observed z > 6 [O III ] 88 μm -SFR and [C II ] 158 μm -SFR relations can only be reproduced when the C/O abundance ratio is ∼8 ×lower than Solar and the total metal production is ∼4 ×higher than that of a Kroupa IMF. This implies that high-redshift galaxies are potentially primarily enriched by low-metallicity core-collapse supernovae with a more top-heavy IMF. As AGB stars and type-Ia supernova begin to contribute to the galaxy metallicity, both the [C II ] 158 μm -SFR and [C II ] 158 μm luminosity functions are predicted to converge to observed values at z ∼4.5. While we demonstrate that ionization parameter, LyC escape fraction, ISM gas density, and CMB attenuation all drive galaxies towards higher [O III ] 88 μm /[C II ] 158 μm , observed values at z > 6 can only be reproduced with substantially lower C/O abundances compared to Solar. The combination of [C II ] 158 μm and [O III ] 88 μm can be used to predict the values of ionization parameter, ISM gas density, and LyC escape fraction and we provide estimates of these quantities for nine observed z > 6 galaxies. Finally, we demonstrate that [O I ] 63 μm can be used as a replacement for [C II ] 158 μm in high-redshift galaxies where [C II ] 158 μm is unobserved and argue that more observation time should be used to target [O I ] 63 μm at z > 6. Future simulations will be needed to self-consistently address the numerous uncertainties surrounding a varying IMF at high redshift and the associated metal returns.

Limited constraints on the evolution of the Lyman Continuum (LyC) escape fraction represent one of the primary uncertainties in the theoretical determination of the reionization history. Due to the intervening intergalactic medium (IGM), the possibility of observing LyC photons directly in the epoch of reionization is highly unlikely. For this reason, multiple indirect probes of LyC escape have been identified, some of which are used to identify low-redshift LyC leakers (e.g. O32), while others are primarily useful at $z>6$ (e.g. [OIII]/[CII] far infrared emission). The flux ratio of the resonant MgII doublet emission at 2796$\dot{\rm A}$ and 2803$\dot{\rm A}$ as well as the MgII optical depth have recently been proposed as ideal diagnostics of LyC leakage that can be employed at $z>6$ with JWST. Using state-of-the-art cosmological radiation hydrodynamics simulations post-processed with CLOUDY and resonant-line radiative transfer, we test whether MgII is indeed a useful probe of LyC leakage. Our simulations indicate that the majority of bright, star-forming galaxies with high LyC escape fractions are expected to be MgII emitters rather than absorbers at $z=6$. However, we find that the MgII doublet flux ratio is a more sensitive indicator of dust rather than neutral hydrogen, limiting its use as a LyC leakage indicator to only galaxies in the optically thin regime. Given its resonant nature, we show that MgII will be an exciting probe of the complex kinematics in high-redshift galaxies in upcoming JWST observations., 25 Pages, 20 Figures, MNRAS Accepted

We measure escape fractions, $f_{\rm esc}$, of ionizing radiation from galaxies in the SPHINX suite of cosmological radiation-hydrodynamical simulations of reionization, resolving halos with $M_{\rm vir} \gtrapprox 7.5 \times 10^7 \ M_{\odot}$ with a minimum cell width of $\approx 10$ pc. Our new and largest $20$ co-moving Mpc wide volume contains tens of thousands of star-forming galaxies with halo masses up to a few times $10^{11} \ M_{\odot}$. The simulated galaxies agree well with observational constraints of the UV luminosity function in the Epoch of Reionization. The escape fraction fluctuates strongly in individual galaxies over timescales of a few Myrs, due to its regulation by supernova and radiation feedback, and at any given time a tiny fraction of star-forming galaxies emits a large fraction of the ionizing radiation escaping into the inter-galactic medium. Statistically, $f_{\rm esc}$ peaks in intermediate-mass, intermediate-brightness, and low-metallicity galaxies ($M_{*} \approx 10^7 \ M_{\odot}$, $M_{1500} \approx -17$, $Z\lesssim 5 \times 10^{-3} \ Z_{\odot}$), dropping strongly for lower and higher masses, brighter and dimmer galaxies, and more metal-rich galaxies. The escape fraction correlates positively with both the short-term and long-term specific star formation rate. According to SPHINX, galaxies too dim to be yet observed, with $M_{1500} \gtrapprox -17$, provide about $55$ percent of the photons contributing to reionization. The global averaged $f_{\rm esc}$ naturally decreases with decreasing redshift, as predicted by UV background models and low-redshift observations. This evolution is driven by decreasing specific star formation rates over cosmic time., Comment: 20 figures, 21 pages plus appendices. Accepted for publication in MNRAS

Understanding the escape of Lyman continuum (LyC) and Lyman $\alpha$ (Ly$\alpha$) photons from giant molecular clouds (GMCs) is crucial if we are to study the reionization of the Universe and to interpret spectra of observed galaxies at high redshift. To this end, we perform high-resolution, radiation-magneto-hydrodynamic simulations of GMCs with self-consistent star formation and stellar feedback. We find that a significant fraction (15-70%) of ionizing radiation escapes from the simulated GMCs with different masses ($10^5$ and $10^6\,M_\odot$), as the clouds are dispersed within about $2$-$5\,{\rm Myr}$ from the onset of star formation. The fraction of LyC photons leaked is larger when the GMCs are less massive, metal-poor, less turbulent, and less dense. The most efficient leakage of LyC radiation occurs when the total star formation efficiency of a GMC is about 20%. The escape of Ly$\alpha$ shows a trend similar to that of LyC photons, except that the fraction of Ly$\alpha$ photons escaping from the GMCs is larger ($f_{\rm esc}^{\rm Ly\alpha}\approx f_{900}^{0.27}$) and that a GMC with strong turbulence shows larger $f_{\rm esc}^{\rm Ly\alpha}$. The simulated GMCs show a characteristic velocity separation of $\Delta v\approx 120 \,{\rm km\,s^{-1}}$ in the time-averaged emergent Ly$\alpha$ spectra, suggesting that Ly$\alpha$ could be useful to infer the kinematics of the interstellar and circumgalactic medium. We show that Ly$\alpha$ luminosities are a useful indicator of the LyC escape, provided the number of LyC photons can be deduced through stellar population modeling. Finally, we find that the correlations between the escape fractions of Ly$\alpha$, ultraviolet photons at 1500A, and the Balmer $\alpha$ line are weak., Comment: 32 pages, 20 figures, Accepted for publication in the ApJS

Aims. The primary difficulty in understanding the sources and processes that powered cosmic reionization is that it is not possible to directly probe the ionizing Lyman-continuum (LyC) radiation at that epoch as those photons have been absorbed by the intervening neutral hydrogen. It is therefore imperative to build a model to accurately predict LyC emission using other properties of galaxies in the reionization era. Methods. In recent years, studies have shown that the LyC emission from galaxies may be correlated to their Lyman-alpha (Lyα) emission. In this paper we study this correlation by analyzing thousands of simulated galaxies at high redshift in the SPHINX cosmological simulation. We post-process these galaxies with the Lyα radiative transfer code RASCAS and analyze the Lyα – LyC connection. Results. We find that the Lyα and LyC luminosities are strongly correlated with each other, although with dispersion. There is a positive correlation between the escape fractions of Lyα and LyC radiations in the brightest Lyman-alpha emitters (LAEs; escaping Lyα luminosity LescLyα > 1041 erg s−1), similar to that reported by recent observational studies. However, when we also include fainter LAEs, the correlation disappears, which suggests that the observed relation may be driven by selection effects. We also find that the brighter LAEs are dominant contributors to reionization, with LescLyα > 1040 erg s−1 galaxies accounting for > 90% of the total amount of LyC radiation escaping into the intergalactic medium in the simulation. Finally, we build predictive models using multivariate linear regression, where we use the physical and Lyα properties of simulated reionization era galaxies to predict their LyC emission. We build a set of models using different sets of galaxy properties as input parameters and predict their intrinsic and escaping LyC luminosity with a high degree of accuracy (the adjusted R2 of these predictions in our fiducial model are 0.89 and 0.85, respectively, where R2 is a measure of how much of the response variance is explained by the model). We find that the most important galaxy properties for predicting the escaping LyC luminosity of a galaxy are its LescLyα, gas mass, gas metallicity, and star formation rate. Conclusions. These results and the predictive models can be useful for predicting the LyC emission from galaxies using their physical and Lyα properties and can thus help us identify the sources of reionization.

We investigate the formation of jellyfish galaxies using radiation-hydrodynamic simulations of gas-rich dwarf galaxies with a multi-phase interstellar medium (ISM). We find that the ram-pressure-stripped (RPS) ISM is the dominant source of molecular clumps in the near wake within 10 kpc from the galactic plane, while in-situ formation is the major channel for dense gas in the distant tail of the gas-rich galaxy. Only 20% of the molecular clumps in the near wake originate from the intracluster medium (ICM); however, the fraction reaches 50% in the clumps located at $80\,{\rm kpc}$ from the galactic center since the cooling time of the RPS gas tends to be short due to the ISM--ICM mixing ($\lesssim$ 10 Myr). The tail region exhibits a star formation rate of $0.001-0.01\,{\rm M_{\odot}\,yr^{-1}}$, and most of the tail stars are born in the stripped wake within 10 kpc from the galactic plane. These stars induce bright H$\alpha$ blobs in the tail, while H$\alpha$ tails fainter than $6\times10^{38}\,{\rm erg\,s^{-1}\,kpc^{-2}}$ are mostly formed via collisional radiation and heating due to mixing. We also find that the stripped tails have intermediate X-ray to H$\alpha$ surface brightness ratios (1.5$\lesssim F_{\rm X}/F_{\rm H\alpha}\lesssim$20), compared to the ISM ($\lesssim$1.5) or pure ICM ($\gg$20). Our results suggest that jellyfish features emerge when the ISM from gas-rich galaxies is stripped by strong ram pressure, mixes with the ICM, and enhances the cooling in the tail., Comment: 19 pages, 15 figures; Accepted for publication in ApJ

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

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

In the last decade we have seen an enormous increase in the size and quality of spectroscopic galaxy surveys, both at low and high redshift. New statistical techniques to analyse large portions of galaxy spectra are now finding favour over traditional index based methods. Here we will review a new robust and iterative Principal Component Analysis (PCA) algorithm, which solves several common issues with classic PCA. Application to the 4000 A break region of galaxies in the VIMOS VLT Deep Survey (VVDS) and Sloan Digital Sky Survey (SDSS) gives new high signal‐to‐noise ratio spectral indices easily interpretable in terms of recent star formation history. In particular, we identify a sample of post‐starburst galaxies at z∼0.7 and z∼0.07. We quantify for the first time the importance of post‐starburst galaxies, consistent with being descendents of gas‐rich major mergers, for building the red sequence. Finally, we present a comparison with new low and high redshift “mock spectroscopic surveys” derived from a Mi...

One open question in cosmology is the exact origin of cosmic Reionisation. Around redshift 6, or about one billion years after the Big Bang, the state of the hydrogen in the Intergalactic Medium (IGM) went from neutral to ionised, which constitutes the last phase transition of the Universe. The precise origin of this process is still debated, although the current consensus is that the ionising photons emitted by the young massive stars in early galaxies are responsible for Reionisation. However, most of the ionising photons, also called Lyman continuum (LyC), are absorbed by the Interstellar Medium (ISM) of the galaxies from which they are emitted. Therefore, to better understand which type of galaxy contributed the most to Reionisation, one must know how to measure the escape fraction of ionising photons (fesc) from high redshift galaxies. Since the ionising photons are absorbed by the IGM during Reionisation, they cannot be observed. Therefore, the measure of fesc has to be done with indirect tracers. The down-the-barrel absorption lines of low ionisation state (LIS) ions are promising tracers of fesc. They are absorption lines imprinted on the stellar continuum of galaxies by metallic ions such as C+ or Si+. Thanks to their wavelengths longer than the wavelength of the Lyman alpha transition of hydrogen, those absorption lines are not erased by the neutral hydrogen in the IGM. Furthermore, since the lines are from low ionisation elements, they trace the neutral medium of the ISM, which itself shapes fesc. To study how well we can infer fesc from absorption lines, I develop and exploit a unique post-processing pipeline to produce down-the-barrel LIS absorption line mock observations and to compute the escape fractions in our simulated galaxies. This allows me to compare the observables with fesc, and to investigate and quantify correlations and dispersion between them. After presenting the simulations that I am using in this work and the different steps of my post-processing pipeline, I begin by studying absorption lines themselves. I show that I reproduce realistic lines, which look like observed high redshift absorption lines. Then, I show that many complex physical processes affect the absorption profiles. The infilling effect, due to the scattering of line photons, changes the depth of the lines, and the large-scale velocity gradients of the gas in galaxies modify the lines significantly, due to the Doppler effect. Additionally, inhomogeneous dust screens also play an important role in the shaping of absorption lines. Then, I compare the properties of the lines with the escape fractions of ionising photons in selected directions of observations. Due to the complex processes affecting the lines, and to the imperfect co-location of LIS ions and neutral hydrogen, I find relations with a large scatter, and the absorption lines cannot in general trace the escape fraction of ionising photons. To investigate further the origin of the dispersion in the relations between escape fraction and absorption lines properties, I first analyse the complex but central role of dust in shaping the UV outputs of galaxies, and the impact of our dust modelling. By comparing observables such as the luminosity function, the beta slope and the infrared excess between simulations and observations, I assess how realistic our simulations are. I also test the observational methods used to compute the dust attenuation of high redshift galaxies. By using them on my mock spectra, and comparing with the exact dust attenuation known from the simulation, I deduce that the extinctions measured in observations are usually underestimated. Finally, I explore the link between absorption lines and gas kinematics. I compare the lines with mass outflow and inflow rates in the Circumgalactic Medium (CGM) of our simulations, to study how well they are correlated.; Une question ouverte en cosmologie est l’origine exacte de la Réionisation cosmique. Autour du redshift 6, soit environ un milliard d’années après le Big Bang, l’état de l’hydrogène dans le milieu intergalactique (IGM) est passé de neutre à ionisé, ce qui caractérise la dernière transition de phase de l’Univers. L’explication précise de ce processus est encore débattue, bien qu’il semble que les photons ionisants émis par les jeunes étoiles massives dans les premières galaxies soient responsables de la Réionisation. Cependant, la plupart des photons ionisants sont absorbés par le milieu interstellaire (ISM) des galaxies à partir desquelles ils sont émis. Par conséquent, pour mieux comprendre quel type de galaxie a le plus contribué à la Réionisation, il faut savoir mesurer la fraction d’échappement des photons ionisants (fesc) des galaxies à hauts redshifts. Les photons ionisants étant absorbés par l'IGM lors de la Réionisation, ils ne peuvent pas être observés depuis la Terre. Par conséquent, la mesure de fesc doit être faite avec des traceurs indirects. Les raies d’absorption des ions à faible ionisation sont des traceurs prometteurs de fesc. Ce sont des raies d’absorption créées par des ions métalliques tels que C+ ou Si+, apparaissant sur le continu stellaire des galaxies. Grâce à leur longueur d’onde plus longue que celle de la transition Lyman alpha de l’hydrogène, ces raies d’absorption ne sont pas effacées par l'IGM neutre de l’époque de la Réionisation. En outre, comme les raies viennent d'éléments à faible ionisation, elles tracent le milieu neutre de l’ISM, qui lui-même détermine fesc. Pour étudier dans quelle mesure nous pouvons déduire fesc des raies d’absorption, j’utilise mes codes de post-traitement uniques en leur genre pour produire des raies d’absorption et calculer fesc dans nos galaxies simulées. Cela me permet de comparer les deux quantités et de trouver des corrélations entre elles. Après avoir présenté les simulations que j’utilise dans ce travail et les étapes de mes codes de post-traitement, je commence par étudier les raies d’absorption elles-mêmes. Je montre que je reproduis des raies réalistes, ressemblant à des raies d’absorption de galaxies observées à redshift élevé. Ensuite, je montre que de nombreux processus physiques complexes affectent les profils d’absorption. L’effet de remplissage, dû à la diffusion des photons, modifie la profondeur des raies, et les gradients de vitesse du gaz dans les galaxies modifient considérablement les raies, en raison de l’effet Doppler. De plus, les écrans de poussière non homogènes jouent également un rôle important dans le façonnage des raies d’absorption. Ensuite, je compare les raies avec les fractions d’échappement des photons ionisants dans des directions d’observations choisies. En raison des processus complexes affectant les raies, et des différences entre la distribution des ions et celle de l’hydrogène neutre, je ne trouve que des relations dispersées, et les raies d’absorption ne peuvent en général pas prédire la fraction d’échappement des photons ionisants. Je continue l’analyse du continuum ultraviolet et des raies d’absorption pour étudier l'ISM et les flux de gaz dans les simulations. En comparant des quantités observables telles que la fonction de luminosité, la pente beta et l’excès infrarouge entre les simulations et les observations, j’évalue à quel point nos simulations sont réalistes. Je teste également les méthodes d’observation utilisées pour calculer l’atténuation par la poussière des galaxies à redshift élevé. En utilisant ces méthodes sur mes spectres simulés, et en comparant les résultats avec l’atténuation par la poussière exacte, connue dans la simulation, je déduis que les mesures sous-estiment la vraie atténuation. Enfin, je compare les raies d’absorption avec les débits de masse de gaz sortants et entrants dans le milieu circum-galactique de nos simulations, pour étudier dans quelle mesure ces quantités sont corrélées.

The Epoch of Reionisation, which spans during the first billion year of te Universe, corresponds to the period during which the first stars and galaxy form. In this context, the main topic of this thesis is to study the formation of those early structures and their radiative feedback to their environment. For this purpose, I use various numerical simulations tools designed to model the radiative transfer in galaxies in a cosmological framework.More specifically, I look at very small galaxies, which are believed to contribute the bulk of the photons required to reionise the Universe. I explore this idea using radiative hydrodynamics simulations performed with RAMSES-RT, focusing on three small galaxies with a very high spatial and temporal resolution. I first detail the mechanism that regulate te production and escape of ionising photons in galaxies, and I show that supernovae explosions are a crucial element for this regulation. I then started to investigate the observable properties of those galaxies.In a second part of my thesis, I developped an extension to the Monte Carlo radiative transfer code MCLya to take light polarisation into account and to model the diffuse emission. I applied this code to post-process the simulation of a Lyman-alpha blob (an extended Lyman-alpha source), and to study its polarisation properties with mock observations. Contrary to what was suggested before, I showed that polarisation is not a strong tracer of the origin of Lyman-alpha photons, L'époque de la réionisation, qui s'étend pendant le premier milliard d'années de l'Univers, correspond à la période où les premières étoiles et galaxies apparaissent. Dans ce contexte, l'enjeu majeur de cette thèse est d'étudier la formation de ces premières structures et leur rétroaction radiative sur leur environnement. Pour cela, j'utilise différentes méthodes de simulations numériques permettant de modéliser le transfert radiatif dans les galaxies, dans un cadre cosmologique.En particulier, les plus petites galaxies seraient celles qui apportent la majorité des photons nécessaires à ioniser l'Univers. J'explore cette hypothèse dans une première partie à l'aide de simulations radio-hydrodynamiques avec le code RAMSES-RT, ciblant trois petites galaxies avec une très haute résolution spatiale et temporelle. Je me suis d'abord intéressé aux mécanismes régulant la production et le transfert de photons ionisants dans les galaxies, et j'ai montré que les explosions de supernovae en sont un facteur essentiel. Ensuite, j'ai exploré les propriétés observables dérivées à partir de ces simulations.Dans un second temps, j'ai développé une extension au code Monte Carlo de transfert radiatif MCLya prenant en compte la polarisation du rayonnement et l'émission diffuse dans une simulation. J'ai utilisé ce code pour post-traiter une simulation d'un blob Lyman-alpha, une source étendue d'émission Lyman-alpha, et étudier ses propriétés de polarisation à l'aide de pseudo-observations. Contrairement à ce qui était proposé précédemment, j'ai pu montrer que la polarisation n'était pas un indicateur aussi utile qu'espéré pour tracer l'origine des photons Lyman-alpha

With the increasing sophistication and efficiency of cosmological hydrodynamics codes, ithas become viable to include ionizing radiative transfer (RT) in cosmological simulations,either in post-processing or in full-blown radiation-hydrodynamics (RHD) simulations. Inspite of the many hurdles involved, there has been much activity during the last decade or soon different strategies and implementations, because a number of interesting problems canbe addressed with RT and RHD, e.g. how and when the Universe became reionized, howradiation from stars and active galactic nuclei plays a part in regulating structure formationon small and large scales, and what predictions and interpretations we can make of observedphenomena such as the Lyman-alpha forest and diffuse sources of radiation.This coincides with the advent of the James Webb space telescope (JWST) and otherstate-of-the-art instruments which are about to give us an unprecedented glimpse into theend of the dark ages of the Universe, when the cosmos switched from a cold and neutralstate to a hot and ionized one, due to the turn-on of ionizing radiative sources.With a primary interest in the problem of radiative feedback in early structure formation,we have implemented an RHD version of the Ramses cosmological code we call RamsesRT,which is moment based and employs the local M1 Eddington tensor closure. This code allowsus to study the effects of ionizing radiation on-the-fly in cosmological RHD simulationsthat take full advantage of the adaptive mesh refinement and parallelization strategies ofRamses. For self-consistent RHD we have also implemented a non-equilibrium chemistry ofthe atomic hydrogen and helium species that interact with the transported radiation.I present in this thesis an extensive description of the RamsesRT implementation andnumerous tests to validate it.Thus far we have used the RHD implementation to study extended line emission fromaccretion streams, which are routinely predicted to exist at early redshift by cosmologicalsimulations but have never been unambiguously verified by observations, and to investigatewhether gravitational heating in those streams could be the dominant power source ofso-called Lyman-alpha blobs, an observed phenomenon which has been much studied anddebated during the last decade or two. Our conclusions from this investigation are thatLyman-alpha blobs can in principle be powered by gravitational heating, and furthermorethat accretion streams are on the verge of being directly detectable for the first time withupcoming instruments.My future intent is to use RamsesRT for high-resolution cosmological zoom simulations ofearly galaxy formation, up to the epoch of reionization, to study how radiative feedbackaffects the formation and evolution of those galaxies and to make observational predictionsthat can be tested with upcoming instruments such as the JWST., Avec l’essor actuel de la sophistication et de l’efficacité des codes de cosmologie hydrodynamique,il est devenu possible d’inclure le transfert radiatif (RT) des photons ionisants dansles simulations cosmologiques, soit en post-traitement, soit en simulations couplées rayonnement+hydrodynamique (RHD). Malgré de nombreux obstacles, il y a eu cette derniéredécennie beaucoup de recherches menées sur les différentes stratégies et implémentations,dû au fait qu’un nombre de problèmes intéressants peuvent être désormais abordés par laRT et RHD, par exemple comment et quand l’Univers s’est réionisé, comment l’émissionradiative des étoiles et des noyaux actifs de galaxies se comportent pour réguler la formationdes structures à des échelles petites et grandes, et quelles prédictions et interprétationsnous pouvons faire des phénomènes observés, tels que la forêt Lyman-alpha et des sourcesdiffuses de rayonnement.Cela coïncide avec l’avènement du télescope spatial James Webb (JWST) et d’autresinstruments de pointe qui sont sur le point de nous donner un aperçu sans précédent sur lafin des âges sombres de l’Univers, quand le cosmos est passé d’un état froid et neutre à unétat chaud et ionisé, à la suite de l’apparition des sources radiatives.Notre préoccupation principale étant les rétroactions radiatives des premieres structures,nous avons mis en place une version RHD du code cosmologique Ramses, que nous appelonsRamsesRT, basée sur la méthode des moments. Ce code nous permet d’étudier les effets durayonnement ionisant dans les simulations cosmologiques RHD qui tirent pleinement profitdes stratégies de raffinement adaptif de grille et de parallélisation de Ramses. Pour rendreauto-cohérent le RHD nous avons également mis en oeuvre une thermochimie hors-équilibreincluant des espèces de l’Hydrogène et de l’Hélium qui interagissent avec le rayonnementtransporté.Je présente dans cette thèse une description détaillée de RamsesRT et de nombreux testscontribuant à sa validation.Jusqu’à présent nous avons utilisé RamsesRT pour étudier l’émission Lyman-alpha decourants d’accrétion, qui sont prédits à grand redshift par les simulations cosmologiques,mais n’ont jamais été clairement identifiés par les observations. Nous avons également étudiéle chauffage gravitationnel dans ces courants pour déterminer si ce dernier pouvait être lasource motrice principale des Lyman-alpha blobs, un phénomène observé qui a été beaucoupétudié et débattu au cours de la dernière décennie. Cet étudie nous permet de conclure queles Lyman-alpha blobs peuvent, en principe, être alimentés par le chauffage gravitationnel,et que d’autre part, les courants d’accrétion sont sur le point d’être directement détectablesavec des instruments à venir.Mes intentions futures sont d’utiliser RamsesRT dans les simulations cosmologiques àhaute résolution, de la formation des premiéres galaxies jusqu’à l’époque de la réionisation,et ainsi étudier comment la rétroaction radiative affecte la formation et l’évolution de cesgalaxies et de faire des prévisions d’observation qui peuvent être testées avec des instrumentssophistiqués tels que le JWST.