A numerical code for solving various Lyman alpha (Lya) radiative transfer (RT) problems is presented. The code is suitable for an arbitrary, three-dimensional distribution of Lya emissivity, gas temperature, density, and velocity field. Capable of ha
ndling Lya RT in an adaptively refined grid-based structure, it enables detailed investigation of the effects of clumpiness of the interstellar (or intergalactic) medium. The code is tested against various geometrically and physically idealized configurations for which analytical solutions exist, and subsequently applied to three Lyman-break galaxies, extracted from high-resolution cosmological simulations at redshift z = 3.6. Proper treatment of the Lya scattering reveals a diversity of surface brightness (SB) and line profiles. Specifically, for a given galaxy the maximum observed SB can vary by an order of magnitude, and the total flux by a factor of 3 - 6, depending on the viewing angle. This may provide an explanation for differences in observed properties of high-redshift galaxies, and in particular a possible physical link between Lyman-break galaxies and regular Lya emitters.
We compute the escape of ionizing radiation from galaxies in the redshift interval z=4-10, i.e., during and after the epoch of reionization, using a high-resolution set of galaxies, formed in fully cosmological simulations. The simulations invoke ear
ly, energetic feedback, and the galaxies evolve into a realistic population at z=0. Our galaxies cover nearly four orders of magnitude in masses (10^{7.8}-10^{11.5}msun) and more than five orders in star formation rates (10^{-3.5}-10^{1.7}msunyr^{-1}), and we include an approximate treatment of dust absorption. We show that the source-averaged Lyman-limit escape fraction at z=10.4 is close to 80% declining monotonically with time as more massive objects build up at lower redshifts. Although the amount of dust absorption is uncertain to 1-1.5 dex, it is tightly correlated with metallicity; we find that dust is unlikely to significantly impact the observed UV output. These results support reionization by stellar radiation from low-luminosity dwarf galaxies and are also compatible with Lyman continuum observations and theoretical predictions at zsim3-4.
We present results from a numerical study of the multiphase interstellar medium in sub-Lyman-break galaxy protogalactic clumps. Such clumps are abundant at z=3 and are thought to be a major contributor to damped Ly-alpha absorption. We model the form
ation of winds from these clumps and show that during star formation episodes they feature outflows with neutral gas velocity widths up to several hundred km/s. Such outflows are consistent with the observed high-velocity dispersion in DLAs. In our models thermal energy feedback from winds and supernovae results in efficient outflows only when cold (~ 300 K), dense (> 100 msun/pc^3) clouds are resolved at grid resolution of 12 pc. At lower 24 pc resolution the first signs of the multiphase medium are spotted; however, at this low resolution thermal injection of feedback energy cannot yet create hot expanding bubbles around star-forming regions -- instead feedback tends to erase high-density peaks and suppress star formation. At 12 pc resolution feedback compresses cold clouds, often without disrupting the ongoing star formation; at the same time a larger fraction of feedback energy is channeled into low-density bubbles and winds. These winds often entrain compact neutral clumps which produce multi-component metal absorption lines.
The median observed velocity width v_90 of low-ionization species in damped Ly-alpha systems is close to 90 km/s, with approximately 10% of all systems showing v_90 > 210 km/s at z=3. We show that a relative shortage of such high-velocity neutral gas
absorbers in state-of-the-art galaxy formation models is a fundamental problem, present both in grid-based and particle-based numerical simulations. Using a series of numerical simulations of varying resolution and box size to cover a wide range of halo masses, we demonstrate that energy from gravitational infall alone is insufficient to produce the velocity dispersion observed in damped Ly-alpha systems, nor does this dispersion arise from an implementation of star formation and feedback in our highest resolution (~ 45 pc) models, if we do not put any galactic winds into our models by hand. We argue that these numerical experiments highlight the need to separate dynamics of different components of the multiphase interstellar medium at z=3.
