In this paper, we present a new method for the solution of those linear transport processes that may be described by a Master Equation, such as electron, neutron and photon transport, and more exotic variants thereof. We base our algorithm on a Markov process on a Voronoi-Delaunay grid, a nonperiodic lattice which is derived from a random point process that is chosen to optimally represent certain properties of the medium through which the transport occurs. Our grid is locally translation and rotation invariant in the mean. We illustrate our approach by means of a particular example, in which the expectation value of the length of a grid line corresponds to the local mean free path. In this example, the lattice is a direct representation of the `free path space' of the medium. Subsequently, transport is defined as simply moving particles from one node to the next, interactions taking place at each point. We derive the statistical properties of such lattices, describe the limiting behavior, and show how interactions are incorporated as global coefficients. Two elementary linear transport problems are discussed: that of free ballistic transport, and the transport of particles through a scattering medium. We also mention a combination of these two. We discuss the efficiency of our method, showing that it is much faster than most other methods, because the operation count does not scale with the number of sources. (ABRIDGED)

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.

In this Letter, we argue that the standard analytical derivations of properties of HII regions, such as the speed, shape and asymptotic position of ionisation fronts require a more precise treatment. These derivations use the on the spot approximation, which in effect ignores the diffuse component of the radiation field. We show that, in fact, HII regions are diffusion dominated. This has as a result that the morphology of inhomogeneous HII regions will be drastically different, because shadowing effects have a less profound impact on the apparent shape. Moreover, it will have influence on the propagation speed of ionisation fronts. We quantify our claims by analytically deriving the internal radiation structure of HII regions, taking diffusion fully into account for several different cosmologically relevant density distributions.

We present a new numerical approach that is able to solve the multi-dimensional radiative transfer equations in all opacity regimes on a Lagrangian, unstructured network of characteristics based on a stochastic point process. Our method reverses the limiting procedure used to derive the transfer equations, by going back to the original Markov process. Thus, we reduce this highly complex system of coupled differential equations to a simple one-dimensional random walk on a graph, which is shown to be computationally very efficient. Specifically, we use a Delaunay graph, which makes it possible to combine our scheme with a new smoothed particle hydrodynamics (SPH) variant proposed by Pelupessy et al.(2003). We show that the results of applying a two-dimensional implementation of our method with various suitable test cases agree with the analytical results, and we point out the advantages of using our method with inhomogeneous point distributions, showing examples in the progress. Hereafter, we present a supplement to our method, which can be useful in cases where the medium is optically very thin, and we conclude by stating some anticipated properties of this method in three dimensions, and announce future extensions.