With a total integration time of 168 hours and a narrowband (NB) filter tuned to Lyman alpha at z = 8.8, the UltraVISTA survey has set out to find some of the most distant galaxies, on the verge of the Epoch of Reionization. Previous calculations of the expected number of detected Lya-emitting galaxies (LAEs) at this redshift did not explicitly take into account the radiative transfer (RT) of Lya. In this work we combine a theoretical model for the halo mass function with numerical results from high-res cosmological hydrosimulations with LyC+Lya RT, assessing the visibility of LAEs residing in these halos. Uncertainties such as cosmic variance and the anisotropic escape of Lya are taken into account, and it is predicted that once the survey has finished, the probabilities of detecting none, one, or more than one are ~90%, ~10%, and ~1%; a significantly smaller success rate compared to earlier predictions, due to the combined effect of a highly neutral IGM scattering Lya to such large distances from the galaxy that they fall outside the observational aperture, and to the actual depth of the survey being less than predicted. Because the IGM affects NB and broadband (BB) magnitudes differently, we argue for a relaxed color selection criterion of NB - BB ~ +0.85. But since the flux is continuum-dominated, even if a galaxy is detectable in the NB its probability of being selected as a NB excess object is <~35%. Various properties of galaxies at this redshift are predicted, e.g. UV and Lya LFs, M*-Mh relation, spectral shape, optimal aperture, and the anisotropic escape of Lya through both a dusty ISM and a partly neutral IGM. Finally, we describe and publish a fast numerical code for adding numbers with asymmetric uncertainties ("x_{-sigma_1}^{+sigma_2}") proving to be significantly better than the standard, but wrong, way of adding upper and lower uncertainties in quadrature separately.

The paper considers the evolution of the supernova envelopes produced by Population III stars with masses of $M_*\sim 25-200 M_\odot$ located in non-rotating protogalaxies with masses of $M\sim 10^7 M_\odot$ at redshifts $z=12$, with dark-matter density profiles in the form of modified isothermal spheres. The supernova explosion occurs in the ionization zone formed by a single parent star. The properties of the distribution of heavy elements (metals) produced by the parent star are investigated, as well as the efficiency with which they are mixed with the primordial gas in the supernova envelope. In supernovae with high energies ($E\simgt 5\times 10^{52}$ erg), an appreciable fraction of the gas can be ejected from the protogalaxy, but nearly all the heavy elements remain in the protogalaxy. In explosions with lower energies ($E\simlt 3\times 10^{52}$ erg), essentially no gas and heavy elements are lost from the protogalaxy: during the first one to threemillion years, the gas and heavy elements are actively carried from the central region of the protogalaxy ($r\sim 0.1 r_{vir}$, where $r_{vir}$ is the virial radius of the protogalaxy), but an appreciable fraction of the mass of metals subsequently returns when the hot cavity cools and the envelope collapses. Supernovae with high energies ($E\simgt 5\times 10^{52}$ erg) are characterized by a very low efficiency of mixing of metals; their heavy elements are located in the small volume occupied by the disrupted envelope (in a volume comparable with that of the entire envelope), with most of the metals remaining inside the hot, rarified cavity of the envelope. (abridged)

We consider the evolution of the ionization zone around Population III stars with $M_*\sim 25-200 M_\odot$ in protogalaxies with $M\sim 10^7 M_\odot$ at redshifts $z = 12$, assuming that the dark matter profile is a modified isothermal sphere. We study the conditions for the growth of instabilities in the ionization zones. The Rayleigh-Taylor and thermal instabilities develop efficiently in the ionization zones around 25-40 $M_\odot$ stars, while this efficiency is lower for stars with $\sim 120 M_\odot$. For more massive stars ($\sim 200 M_\odot$), the flux of ionizing photons is strong enough to considerably reduce the gas density in the ionization zone, and the typical lifetimes of stars ($\sim 2$ Myr) are insufficient for the growth of instabilities. The gas in a protogalaxy with $M\sim 10^7 M_\odot$ with a 200 $M_\odot$ central star is completely ionized by the end of the star's lifetime; in the case of a 120 $M_\odot$ central star, only one-third of the total mass of gas is ionized. Thus, ionizing photons from stars with $M_*\simlt 120 M_\odot$ cannot leave protogalaxies with $M\simgt 10^7 M_\odot$. If the masses of the central stars are 25 and 40 $M_\odot$, the gas in protogalaxies of this mass remains essentially neutral. We discuss the consequences of the evolution of the ionization zones for the propagation of the envelope after the supernova explosions of the stars and the efficiency of enrichment of the intergalactic medium in heavy elements.

We study the intergalactic transmission of radiation in the vicinity of the Ly{\alpha} wavelength. Simulating sightlines through the intergalactic medium (IGM) in detailed cosmological hydrosimulations, the impact of the IGM on the shape of the line profile from Ly{\alpha} emitting galaxies at redshifts 2.5 to 6.5 is investigated. In particular we show that taking into account the correlation of the density and velocity fields of the IGM with the galaxies, the blue part of the spectrum may be appreciably reduced, even at relatively low redshifts. This may in some cases provide an alternative to the often-invoked outflow scenario, although it is concluded that this model is still a plausible explanation of the many asymmetric Ly{\alpha} profiles observed. Applying the calculated wavelength dependent transmission to simulated spectra from Ly{\alpha} emitting galaxies, we derive the fraction of photons that are lost in the IGM, in addition to what is absorbed internally in the galaxies due to dust. Moreover, by comparing the calculated transmission of radiation blueward of the Ly{\alpha} line, the total optical depth to Thomson scattering of cosmic microwave background, with corresponding observations, we are able to constrain the epoch when the Universe was reionized to z <~ 8.5.

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 early, 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}\msun\yr^{-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 z\sim3-4.

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 handling 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 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 formation 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.

Based on cosmological simulations, we model Lyman continuum emission from a sample of 11 high-redshift star forming galaxies spanning a mass range of a factor 20. Each of the 11 galaxies has been simulated both with a Salpeter and a Kroupa initial mass function (IMF). We find that the Lyman continuum (LyC) luminosity of an average star forming galaxy in our sample declines from z=3.6 to 2.4 due to the steady gas infall and higher gas clumping at lower redshifts, increasingly hampering the escape of ionizing radiation. The galaxy-to-galaxy variation of apparent LyC emission at a fixed redshift is caused in approximately equal parts by the intrinsic variations in the LyC emission and by orientation effects. The combined scatter of an order of magnitude can explain the variance in the far-UV spectra of high-redshift galaxies detected by Shapley et al. (2006). Our results imply that the cosmic galactic ionizing UV luminosity would be monotonically decreasing from z=3.6 to 2.4, curiously anti-correlated with the star formation rate in the smaller galaxies, which on average rises during this redshift interval.

Using results from high-resolution galaxy formation simulations in a standard Lambda-CDM cosmology and a fully conservative multi-resolution radiative transfer code around point sources, we compute the energy-dependent escape fraction of ionizing photons from a large number of star forming regions in two galaxies at five different redshifts from z=3.8 to 2.39. All escape fractions show a monotonic decline with time, from (at the Lyman-limit) ~6-10% at z=3.6 to ~1-2% at z=2.39, due to higher gas clumping at lower redshifts. It appears that increased feedback can lead to higher f_esc at z>3.4 via evacuation of gas from the vicinity of star forming regions and to lower f_esc at z<2.39 through accumulation of swept-up shells in denser environments. Our results agree well with the observational findings of \citet{inoue..06} on redshift evolution of f_esc in the redshift interval z=2-3.6.

We use high resolution Eulerian hydrodynamics simulations to study kinematic properties of the low ionization species in damped Ly-alpha systems at redshift z=3. Our adaptive mesh refinement simulations include most key ingredients relevant for modeling neutral gas in high-column density absorbers: hydrodynamics, gravitational collapse, continuum radiative transfer and gas chemistry, but no star formation. We model high-resolution Keck spectra with unsaturated low ion transitions in two Si II lines (1526 and 1808 A), and compare simulated line profiles to the data from the SDSS DLA survey. We find that with increasing grid resolution the models show a trend in convergence towards the observed distribution of HI column densities. While in our highest resolution model we recover the cumulative number of DLAs per unit absorption distance, none of our models predicts DLA velocity widths as high as indicated by the data, suggesting that feedback from star formation might be important. At z=3 a non-negligible fraction of DLAs with column densities below 10^21 cm^-2 is caused by tidal tails due to galaxy-galaxy interactions in more massive halo environments. Lower column density absorbers with N_HI < 10^21.4 cm^-2 are sensitive to changes in the UV background resulting in a 10% reduction of the cumulative number of DLAs for twice the quasar background relative to the fiducial value. We find that the mass cut-off below which a large fraction of dwarf galaxies cannot retain gas after reionization is 7*10^7 msun, lower than the previous estimates. Finally, we show that models with self-shielding commonly used in the literature produce significantly lower DLA velocity widths than the full radiative transfer runs.

Radiative transfer simulations are now at the forefront of numerical astrophysics. They are becoming crucial for an increasing number of astrophysical and cosmological problems; at the same time their computational cost has come to the reach of currently available computational power. Further progress is retarded by the considerable number of different algorithms (including various flavours of ray-tracing and moment schemes) developed, which makes the selection of the most suitable technique for a given problem a non-trivial task. Assessing the validity ranges, accuracy and performances of these schemes is the main aim of this paper, for which we have compared 11 independent RT codes on 5 test problems: (0) basic physics, (1) isothermal H II region expansion and (2) H II region expansion with evolving temperature, (3) I-front trapping and shadowing by a dense clump, (4) multiple sources in a cosmological density field. The outputs of these tests have been compared and differences analyzed. The agreement between the various codes is satisfactory although not perfect. The main source of discrepancy appears to reside in the multi-frequency treatment approach, resulting in different thicknesses of the ionized-neutral transition regions and different temperature structure. The present results and tests represent the most complete benchmark available for the development of new codes and improvement of existing ones. To this aim all test inputs and outputs are made publicly available in digital form.

The computational difficulty of six-dimensional neutrino radiation hydrodynamics has spawned a variety of approximations, provoking a long history of uncertainty in the core-collapse supernova explosion mechanism. Under the auspices of the Terascale Supernova Initiative, we are honoring the physical complexity of supernovae by meeting the computational challenge head-on, undertaking the development of a new adaptive mesh refinement code for self-gravitating, six-dimensional neutrino radiation magnetohydrodynamics. This code--called {\em GenASiS}, for {\em Gen}eral {\em A}strophysical {\em Si}mulation {\em S}ystem--is designed for modularity and extensibility of the physics. Presently in use or under development are capabilities for Newtonian self-gravity, Newtonian and special relativistic magnetohydrodynamics (with `realistic' equation of state), and special relativistic energy- and angle-dependent neutrino transport--including full treatment of the energy and angle dependence of scattering and pair interactions.

The computational difficulty of six-dimensional neutrino radiation hydrodynamics has spawned a variety of approximations, provoking a long history of uncertainty in the core-collapse supernova explosion mechanism. Under the auspices of the Terascale Supernova Initiative, we are honoring the physical complexity of supernovae by meeting the computational challenge head-on, undertaking the development of a new adaptive mesh refinement code for self-gravitating, six-dimensional neutrino radiation magnetohydrodynamics. This code--called {\em GenASiS,} for {\em Gen}eral {\em A}strophysical {\em Si}mulation {\em S}ystem--is designed for modularity and extensibility of the physics. Presently in use or under development are capabilities for Newtonian self-gravity, Newtonian and special relativistic magnetohydrodynamics (with `realistic' equation of state), and special relativistic energy- and angle-dependent neutrino transport--including full treatment of the energy and angle dependence of scattering and pair interactions.

A new, very fast method for 3D radiative transfer on fully threaded grids with arbitrarily high angular resolution is presented. The method uses completely cell-based discretization, and is ideally suited for problems with diffuse background radiation, often encountered in cosmological and star formation models. We find that for accurate statistical study of intergalactic Ly-alpha absorption lines one needs of order of few hundred angular discretization elements even for models without radiative feedback from star forming galaxies.

We present new calculations of the inhomogeneous process of cosmological reionization by carefully following the radiative transfer in pre-computed hydrodynamical simulations of galaxy formation. These new computations represent an important step on the way towards fully self-consistent and adaptive calculations which will eventually cover the enormous range of scales from the sizes of individual mini-halos to the mean free path of ionizing photons in the post-overlap Universe. The goal of such simulations is to include enough realistic physics to accurately model the formation of early structures and the end of the `dark ages'. Our new calculations demonstrate that the process by which the ionized regions percolate the Universe is complex, and that the idea of voids being ionized before overdense regions is perhaps too simplistic. It seems that observational information pertaining to the reionization epoch may now be in our grasp, through the detection of Gunn-Peterson troughs at z~6. If so, then the comparison of information from many lines of sight with simulations such as ours may allow us to disentangle details of the ionization history and trace the early formation of structure.

A numerical scheme is proposed for the solution of the three-dimensional radiative transfer equation with variable optical depth. We show that time-dependent ray tracing is an attractive choice for simulations of astrophysical ionization fronts, particularly when one is interested in covering a wide range of optical depths within a 3D clumpy medium. Our approach combines the explicit advection of radiation variables with the implicit solution of local rate equations given the radiation field at each point. Our scheme is well suited to the solution of problems for which line transfer is not important, and could, in principle, be extended to those situations also. This scheme allows us to calculate the propagation of supersonic ionization fronts into an inhomogeneous medium. The approach can be easily implemented on a single workstation and also should be fully parallelizable.

Hydrodynamical calculations undertaken to simulate the collisional interaction between the Sgr dwarf galaxy and the Galactic outer HI disk are presented, constrained by recently derived orbital and mass parameters for this dwarf galaxy. It is found that a significant distortion to the structure of the Galactic HI disk will be induced by the collision if the mass of the dwarf exceeds ~10^9 M_Sun; this value is consistent with an estimate derived by requiring that the dwarf galaxy is sufficiently robust to survive tidal disruption until the present time. Though the precise details of the interaction are compromised in our simulations by the lack of a live Galactic halo, we find that for model masses > 5x10^9 M_Sun, prominent spiral arms and a substantial lopsidedness in the outer disk are produced. Furthermore, a noticeable warp-like structure is induced in the disk. Thus the Sgr dwarf may have significantly affected the star formation history and structure of the outer Galaxy. These simulations confirm the possibility of determining the current merging rate of low surface brightness, gas-poor dwarf galaxies of mass > 10^9 M_Sun onto giant spiral galaxies from careful analysis of observations of the structure of HI disks.