Artist: fast radiative transfer for large-scale simulations of the epoch of reionization

ABSTRACT We introduce the ‘Asymmetric Radiative Transfer In Shells Technique’ (artist), a new method for photon propagation on large scales that explicitly conserves photons, propagates photons at the speed of light, approximately accounts for photon directionality, and closely reproduces results of more detailed radiative transfer (RT) methods. Crucially, it is computationally fast enough to evolve the large cosmological volumes required to predict the 21cm power spectrum on scales that will be probed by future experiments targeting the epoch of reionization (EoR). Most seminumerical models aimed at predicting the EoR 21cm signal on these scales use an excursion set formalism (ESF) to model the gas ionization, which achieves computational viability by making a number of approximations. While artist is still roughly two orders of magnitude slower than ESF, it does allow to model the EoR without the need for such approximations. This is particularly important when considering a wide range of reionization scenarios for which artist would help limit the assumptions made. By implementing our RT method within the seminumerical code simfast21, we show that Artist predicts a significantly different evolution for the EoR ionization field compared to the code’s native ESF. In particular, artist predicts up to a factor of two difference in the power spectra, depending on the physical parameters assumed. Its application to large-scale EoR simulations will therefore allow more physically motivated constraints to be obtained for key EoR parameters. In particular, it will remove the need for the artificial rescaling of the escape fraction.

Cosmological filaments in the light of excursion set of saddle points

ABSTRACT The universe in large scales is structured as a network known as cosmic web. Filaments are one of the structural components of this web, which can be introduced as a novel probe to study the formation and evolution of structures and as a probe to study the cosmological models and to address the missing baryon problem. The aim of this work is to introduce an analytical framework to study the statistics of filaments such as number density of them and also to obtain the length-mass relation. For this objective, we model filaments as collapsed objects which have an extension in one direction, accordingly we use the ellipsoidal collapse to study the evolution of an over-dense region via gravitational instability. We find that the non-linear density of filaments in the epoch of formation is almost mass independent and is in order of ∼30. By introducing filament’s extended condition, we find a fitting function for length-mass relation. For the statistics of filaments, we propose a novel framework named excursion set of saddle points. In this approach, we count the saddle points of the density field Hessian matrix, and relate it to the count of filaments. In addition, we addressed the filament in filament problem with up-crossing approximation.

The Neutral Islands during the Late Epoch of Reionization

The large-scale structure of the ionization field during the epoch of reionization (EoR) can be modeled by the excursion set theory. While the growth of ionized regions during the early stage are described by the “bubble model”, the shrinking process of neutral regions after the percolation of the ionized region calls for an “island model”. An excursion set based analytical model and a semi-numerical code (islandFAST) have been developed. The ionizing background and the bubbles inside the islands are also included in the treatment. With two kinds of absorbers of ionizing photons, i.e. the large-scale under-dense neutral islands and the small-scale over-dense clumps, the ionizing background are self-consistently evolved in the model.

Excursion sets of infinitely divisible random fields with convolution equivalent Lévy measure

We consider a continuous, infinitely divisible random field in ℝd, d = 1, 2, 3, given as an integral of a kernel function with respect to a Lévy basis with convolution equivalent Lévy measure. For a large class of such random fields, we compute the asymptotic probability that the excursion set at level x contains some rotation of an object with fixed radius as x → ∞. Our main result is that the asymptotic probability is equivalent to the right tail of the underlying Lévy measure.

A precise mass function in the excursion set approach

In the present paper, using previous results from Del Popolo papers, we show how the mass function evolution can be obtained in the framework of a spherical collapse model, which has been modified to take account of dynamical friction, the cosmological constant, and angular momentum which proto-structures acquire through tidal interaction with neighbouring ones. We found an improved barrier which is in excellent agreement with simulations. The quoted barrier is used to calculated the mass function. In the case of the ΛCDM paradigm, our mass function is in good agreement (within some %) with the mass function of Klypin’s Bolshoi simulation for the virial mass range 5 × 10