We explore scenarios for reionizing the intergalactic medium with low galaxy ionizing photon escape fractions. We combine simulation-based halo-mass dependent escape fractions with an extrapolation of the observed galaxy rest-ultraviolet luminosity functions to solve for the reionization history from z=20 to z=4. We explore the posterior distributions for key unknown quantities, including the limiting halo mass for star-formation, the ionizing photon production efficiency, and a potential contribution from active galactic nuclei (AGN). We marginalize over the allowable parameter space using a Markov Chain Monte Carlo method, finding a solution which satisfies the most model-independent constraints on reionization. Our fiducial model can match observational constraints with an average escape fraction of <5% throughout the bulk of the epoch of reionization if: i) galaxies form stars down to the atomic cooling limit before reionization and a photosuppression mass of log(M_h/Msol)~9 during/after reionization (-13<M_UV,lim<-11); ii) galaxies become more efficient producers of ionizing photons at higher redshifts and fainter magnitudes, and iii) there is a significant, but sub-dominant, contribution by AGN at z < 7. In this model the faintest galaxies (M_UV>-15) dominate the ionizing emissivity, leading to an earlier start to reionization and a smoother evolution of the ionized volume filling fraction than models which assume a single escape fraction at all redshifts and luminosities. The ionizing emissivity from this model is consistent with observations at z=4-5 (and below, when extrapolated), in contrast to some models which assume a single escape fraction. Our predicted ionized volume filling fraction at z=7 of Q_HII=78% (+\- 8%) is in ~1-2 sigma tension with observations of Lya emitters at z~7 and the damping wing analyses of the two known z>7 quasars, which prefer Q_HII,z=7~40-50%.

The dark ages of the Universe end with the formation of the first generation of stars residing in primeval galaxies. These objects were the first to produce ultraviolet ionizing photons in a period when the cosmic gas changed from a neutral state to an ionized one, known as Epoch of Reionization (EoR). A pivotal aspect to comprehend the EoR is to probe the intertwined relationship between the fraction of ionizing photons capable to escape dark haloes, also known as the escape fraction ($f_{esc}$), and the physical properties of the galaxy. This work develops a sound statistical model suitable to account for such non-linear relationships and the non-Gaussian nature of $f_{esc}$. This model simultaneously estimates the probability that a given primordial galaxy starts the ionizing photon production and estimates the mean level of the $f_{esc}$ once it is triggered. The model was employed in the First Billion Years simulation suite, from which we show that the baryonic fraction and the rate of ionizing photons appear to have a larger impact on $f_{esc}$ than previously thought. A naive univariate analysis of the same problem would suggest smaller effects for these properties and a much larger impact for the specific star formation rate, which is lessened after accounting for other galaxy properties and non-linearities in the statistical model.

In this paper we introduce a new radiative transfer code SPRAI (Simplex Photon Radiation in the Arepo Implementation) based on the SimpleX radiation transfer method. This method, originally used only for post-processing, is now directly integrated into the Arepo code and takes advantage of its adaptive unstructured mesh. Radiated photons are transferred from the sources through the series of Voronoi gas cells within a specific solid angle. From the photon attenuation we derive corresponding photon fluxes and ionization rates and feed them to a primordial chemistry module. This gives us a self-consistent method for studying dynamical and chemical processes caused by ionizing sources in primordial gas. Since the computational cost of the SimpleX method does not scale directly with the number of sources, it is convenient for studying systems such as primordial star-forming halos that may form multiple ionizing sources.

The diffuse soft X-ray emissivity from galactic winds is computed during the Epoch of Reionization (EoR). We consider two analytic models, a pressure-driven wind and a superbubble model, and a 3D cosmological simulation including gas dynamics from the First Billion Years (FiBY) project. The analytic models are normalized to match the diffuse X-ray emissivity of star-forming galaxies in the nearby Universe. The cosmological simulation uses physically motivated star formation and wind prescriptions, and includes radiative transfer corrections. The models and the simulation all are found to produce sufficient heating of the Intergalactic Medium to be detectable by current and planned radio facilities through 21 cm measurements during the EoR. While the analytic models predict a 21 cm emission signal relative to the Cosmic Microwave Background sets in by $z_{\rm trans} \simeq 8 - 10$, the predicted signal in the FiBY simulation remains in absorption until reionization completes. The 21 cm absorption differential brightness temperature reaches a minimum of $\Delta T \simeq -130$ to $-200$ mK, depending on model. Allowing for additional heat from high mass X-ray binaries pushes the transition to emission to $z_{\rm trans} \simeq 10 - 12$, with shallower absorption signatures having a minimum of $\Delta T \simeq -110$ to $-140$ mK. The 21 cm signal may be a means of distinguishing between the wind models, with the superbubble model favouring earlier reheating. While an early transition to emission may indicate X-ray binaries dominate the reheating, a transition to emission as early as $z_{\rm trans} > 12$ would suggest the presence of additional heat sources.

We investigate the environment in which direct-collapse black holes may form by analysing a cosmological, hydrodynamical simulation that is part of the First Billion Years project. This simulation includes the most relevant physical processes leading to direct collapse of haloes, most importantly, molecular hydrogen depletion by dissociation of $H_2$ and $H^-$ from the evolving Lyman-Werner radiation field. We selected a sample of pristine atomic cooling haloes that have never formed stars in their past, have not been polluted with heavy elements and are cooling predominantly via atomic hydrogen lines. Amongst them we identified six haloes that could potentially harbour massive seed black holes formed via direct collapse (with masses in the range of $10^{4-6} M_{sun}$). These potential hosts of direct-collapse black holes form as satellites and are found within 15 physical kpc of proto-galaxies, with stellar masses in the range $10^{5-7} M_{sun}$ and maximal star formation rates of 0.1 Msun/yr over the past 5 Myr, and are exposed to the highest flux of Lyman-Werner radiation emitted from the neighbouring galaxies. It is the proximity to these proto-galaxies that differentiates these haloes from rest of the sample.

Using four different suites of cosmological simulations, we generate synthetic spectra for galaxies with different Lyman continuum escape fractions (fesc) at redshifts z=7-9, in the rest-frame wavelength range relevant for the James Webb Space Telescope (JWST) NIRSpec instrument. By investigating the effects of realistic star formation histories and metallicity distributions on the EW(Hb)-beta diagram (previously proposed as a tool for identifying galaxies with very high fesc), we find that neither of these effects are likely to jeopardize the identification of galaxies with extreme Lyman continuum leakage. Based on our models, we expect essentially all z=7-9 galaxies that exhibit rest-frame EW(Hb)< 30 {\AA} to have fesc>0.5. Incorrect assumptions concerning the ionizing fluxes of stellar populations or the dust properties of z>6 galaxies can in principle bias the selection, but substantial model deficiencies of this type will at the same time reveal themselves as an offset between the observed and simulated distribution of z>6 galaxies in the EW(Hb)-beta diagram. Such offsets would thereby allow JWST/NIRSpec measurements of these observables to serve as input for further model refinement.

Direct collapse black holes (DCBHs) formed from the collapse of atomically-cooled primordial gas in the early Universe are strong candidates for the seeds of supermassive BHs. DCBHs are thought to form in atomic cooling haloes in the presence of a strong molecule-dissociating, Lyman-Werner (LW) radiation field. Given that star forming galaxies are likely to be the source of the LW radiation in this scenario, ionizing radiation from these galaxies may accompany the LW radiation. We present cosmological simulations resolving the collapse of primordial gas into an atomic cooling halo, including the effects of both LW and ionizing radiation. We find that in cases where the gas is not self-shielded from the ionizing radiation, the collapse can be delayed by ~ 25 Myr. When the ionized gas does collapse, the free electrons that are present catalyze H2 formation. In turn, H2 cooling becomes efficient in the center of the halo, and DCBH formation is prevented. We emphasize, however, that in many cases the gas collapsing into atomic cooling haloes at high redshift is self-shielding to ionizing radiation. Therefore, it is only in a fraction of such haloes in which DCBH formation is prevented due to reionization.

We present results of full 3D hydrodynamical and radiative transfer simulations of the colliding stellar winds in the massive binary system Eta Carinae. We accomplish this by applying the SimpleX algorithm for 3D radiative transfer on an unstructured Voronoi-Delaunay grid to recent 3D smoothed particle hydrodynamics (SPH) simulations of the binary colliding winds. We use SimpleX to obtain detailed ionization fractions of hydrogen and helium, in 3D, at the resolution of the original SPH simulations. We investigate several computational domain sizes and Luminous Blue Variable primary star mass-loss rates. We furthermore present new methods of visualizing and interacting with output from complex 3D numerical simulations, including 3D interactive graphics and 3D printing. While we initially focus on Eta Car, the methods employed can be applied to numerous other colliding wind (WR 140, WR 137, WR 19) and dusty 'pinwheel' (WR 104, WR 98a) binary systems. Coupled with 3D hydrodynamical simulations, SimpleX simulations have the potential to help determine the regions where various observed time-variable emission and absorption lines form in these unique objects.

Proto-galaxies forming in low-mass dark matter haloes are thought to provide the majority of ionizing photons needed to reionize the Universe, due to their high escape fractions of ionizing photons. We study how the escape fraction in high-redshift galaxies relates to the physical properties of the halo in which the galaxies form, by computing escape fractions in more than 75000 haloes between redshifts 27 and 6 that were extracted from the First Billion Years project, high-resolution cosmological hydrodynamics simulations of galaxy formation. We find that the main constraint on the escape fraction is the gas column density in a radius of 10 pc around the stellar populations, causing a strong mass dependence of the escape fraction. The lower potential well in haloes with virial mass below 1e8 solar mass results in low column densities that can be penetrated by radiation from young stars (age < 5 Myr). In haloes with virial mass above 1e8 solar mass supernova feedback is important, but only 30% of the haloes in this mass range have an escape fraction higher than 1%. We find a large range of escape fractions in haloes with similar properties, caused by different distributions of the dense gas in the halo. This makes it very hard to predict the escape fraction on the basis of halo properties and results in a highly anisotropic escape fraction. The strong mass dependence, the large spread and the large anisotropy of the escape fraction may strongly affect the topology of reionization and is something current models of cosmic reionization should strive to take into account.

Spectral observations of the massive colliding wind binary Eta Carinae show phase-dependent variations, in intensity and velocity, of numerous helium emission and absorption lines throughout the entire 5.54-year orbit. Approaching periastron, the 3D structure of the wind-wind interaction region (WWIR) gets highly distorted due to the eccentric ($e \sim 0.9$) binary orbit. The secondary star ($\eta_{\mathrm{B}}$) at these phases is located deep within the primary's dense wind photosphere. The combination of these effects is thought to be the cause of the particularly interesting features observed in the helium lines at periastron. We perform 3D radiative transfer simulations of $\eta$ Car's interacting winds at periastron. Using the SimpleX radiative transfer algorithm, we post-process output from 3D smoothed particle hydrodynamic simulations of the inner 150 au of the $\eta$ Car system for two different primary star mass-loss rates ($\dot{M}_{\eta_{\mathrm{A}}}$). Using previous results from simulations at apastron as a guide for the initial conditions, we compute 3D helium ionization maps. We find that, for higher $\dot{M}_{\eta_{\mathrm{A}}}$, $\eta_{\mathrm{B}}$ He$^{0+}$-ionizing photons are not able to penetrate into the pre-shock primary wind. He$^{+}$ due to $\eta_{\mathrm{B}}$ is only present in a thin layer along the leading arm of the WWIR and in a small region close to the stars. Lowering $\dot{M}_{\eta_{\mathrm{A}}}$ allows $\eta_{\mathrm{B}}$'s ionizing photons to reach the expanding unshocked secondary wind on the apastron side of the system, and create a low fraction of He$^{+}$ in the pre-shock primary wind. With apastron on our side of the system, our results are qualitatively consistent with the observed variations in strength and radial velocity of $\eta$ Car's helium emission and absorption lines, which helps better constrain the regions where these lines arise.

We present the first 3D prints of output from a supercomputer simulation of a complex astrophysical system, the colliding stellar winds in the massive (>120 M_Sun), highly eccentric (e ~ 0.9) binary star system Eta Carinae. We demonstrate the methodology used to incorporate 3D interactive figures into a PDF journal publication and the benefits of using 3D visualization and 3D printing as tools to analyze data from multidimensional numerical simulations. Using a consumer-grade 3D printer (MakerBot Replicator 2X), we successfully printed 3D smoothed particle hydrodynamics (SPH) simulations of Eta Carinae's inner (r ~ 110 au) wind-wind collision interface at multiple orbital phases. The 3D prints and visualizations reveal important, previously unknown 'finger-like' structures at orbital phases shortly after periastron (phi ~ 1.045) that protrude radially outward from the spiral wind-wind collision region. We speculate that these fingers are related to instabilities (e.g. thin-shell, Rayleigh-Taylor) that arise at the interface between the radiatively-cooled layer of dense post-shock primary-star wind and the fast (3000 km/s), adiabatic post-shock companion-star wind. The success of our work and easy identification of previously unrecognized physical features highlight the important role 3D printing and interactive graphics can play in the visualization and understanding of complex 3D time-dependent numerical simulations of astrophysical phenomena.

The highly eccentric binary system Eta Carinae shows numerous time-variable emission and absorption features. These observational signatures are the result of interactions between the complex three-dimensional (3D) wind-wind collision regions and photoionization by the luminous stars. Specifically, helium presents several interesting spectral features that provide important clues on the geometry and physical properties of the system and the individual stars. We use the SimpleX algorithm to post-process 3D smoothed particle hydrodynamics simulation output of the interacting winds in Eta Car in order to obtain the fractions of ionized helium assuming three different primary star mass-loss rates. The resultant ionization maps constrain the regions where helium is singly- and doubly-ionized. We find that reducing the primary's mass-loss rate increases the volume of He+. Lowering the primary mass-loss rate produces large variations in the volume of He+ in the pre-shock primary wind on the periastron side of the system. Our results show that binary orientations in which apastron is on our side of the system are more consistent with available observations. We suggest that small variations in the primary's mass-loss rate might explain the observed increase in HeI absorption in recent decades, although numerous questions regarding this scenario remain open. We also propose that the absence of broad HeI lines in the spectra of Eta Car between its 1890's eruption and 1944 might be explained by the companion star's He0+ ionizing photons not being able to penetrate the wind-wind interaction region, due to a higher primary mass-loss rate at that time (by a factor >2, compared to the present value).

The contribution of stars in galaxies to cosmic reionisation depends on the star formation history in the Universe, the abundance of galaxies during reionisation, the escape fraction of ionising photons and the clumping factor of the inter-galactic medium (IGM). We compute the star formation rate and clumping factor during reionisation in a cosmological volume using a high-resolution hydrodynamical simulation. We post-process the output with detailed radiative transfer simulations to compute the escape fraction of ionising photons. Together, this gives us the opportunity to assess the contribution of galaxies to reionisation self-consistently. The strong mass and redshift dependence of the escape fraction indicates that reionisation occurred between z=15 and z=10 and was mainly driven by proto-galaxies forming in dark-matter haloes with masses between 1e7 and 1e8 solar mass. More massive galaxies that are rare at these redshifts and have significantly lower escape fractions contribute less photons to the reionisation process than the more-abundant low-mass galaxies. Star formation in the low-mass haloes is suppressed by radiative feedback from reionisation, therefore these proto-galaxies only contribute when the part of the Universe they live in is still neutral. After z~10, massive galaxies become more abundant and provide most of the ionising photons. In addition, we find that Population (Pop) III stars are too short-lived and not frequent enough to have a major contribution to reionisation. Although the stellar component of the proto-galaxies that produce the bulk of ionising photons during reionisation is too faint to be detected by the James Webb Space Telescope (JWST), these sources are brightest in the H-alpha and Ly-alpha recombination lines, which will likely be detected by JWST in deep surveys.

The UV escape fraction from high-redshift galaxies plays a key role in models of cosmic reionisation. Because it is currently not possible to deduce the escape fractions during the epoch of reionisation from observations, we have to rely on numerical simulations. Our aim is to better constrain the escape fraction from high-redshift dwarf galaxies, as these are the most likely sources responsible for reionising the Universe. We employ a N-body/SPH method that includes realistic prescriptions for the physical processes that are important for the evolution of dwarf galaxies. These models are post-processed with radiative transfer to determine the escape fraction of ionising radiation. We perform a parameter study to assess the influence of the spin parameter, gas fraction and formation redshift of the galaxy and study the importance of numerical parameters as resolution, source distribution and local gas clearing. We find that the UV escape fraction from high-redshift dwarf galaxies that have formed a rotationally supported disc lie between 1e-5 and 0.1. The mass and angular momentum of the galaxy are the most important parameters that determine the escape fraction. We compare our results to previous work and discuss the uncertainties of our models. The low escape fraction we find for high-redshift dwarf galaxies is balanced by their high stellar content, resulting in an efficiency parameter for stars that is only marginally lower than the values found by semi-analytic models of reionisation. We therefore conclude that dwarf galaxies play an important role in cosmic reionisation also after the initial starburst phase, when the gas has settled into a disc.

We present an improved version of the SimpleX method for radiative transfer on an unstructured Delaunay grid. The grid samples the medium through which photons are transported in an optimal way for fast radiative transfer calculations. In this paper we study the detailed working of SimpleX in test problems and show improvements over the previous version of the method. We have applied a direction conserving transport scheme that correctly transports photons in the optically thin regime, a regime where the original SimpleX algorithm lost its accuracy. In addition, a scheme of dynamic grid updates is employed to ensure correct photon transport when the optical depth changes during a simulation. For the application to large data sets, the method is parallellised for distributed memory machines using MPI. To test the new method, we have performed standard tests for cosmological radiative transfer. We show excellent correspondence with both the analytical solution (when available) and to the results of other codes compared to the former version of SimpleX, without sacrificing the benefits of the high computational speed of the method.

Context. Analytical and numerical analysis of the SimpleX radiative transfer algorithm, which features transport on a Delaunay triangulation. Aims. Verify whether the SimpleX radiative transfer algorithm conforms to mathematical expectations, to develop error analysis and present improvements upon earlier versions of the code. Methods. Voronoi-Delaunay tessellation, classical Markov theory. Results. Quantitative description of the error properties of the SimpleX method. Numerical validation of the method and verification of the analytical results. Improvements in accuracy and speed of the method. Conclusions. It is possible to transport particles such as photons in a physically correct manner with the SimpleX algorithm. This requires the use of weighting schemes or the modification of the point process underlying the transport grid. We have explored and applied several possibilities.