The Diffuse Nature of Stromgren Spheres

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.Comment: 4 pages, 2 figures, accepted for publication in A&A Letter

RADAMESH: Cosmological Radiative Transfer for Adaptive Mesh Refinement Simulations

We present a new three-dimensional radiative transfer (RT) code, RADAMESH, based on a ray-tracing, photon-conserving and adaptive (in space and time) scheme. RADAMESH uses a novel Monte Carlo approach to sample the radiation field within the computational domain on a "cell-by-cell" basis. Thanks to this algorithm, the computational efforts are now focused where actually needed, i.e. within the Ionization-fronts (I-fronts). This results in an increased accuracy level and, at the same time, a huge gain in computational speed with respect to a "classical" Monte Carlo RT, especially when combined with an Adaptive Mesh Refinement (AMR) scheme. Among several new features, RADAMESH is able to adaptively refine the computational mesh in correspondence of the I-fronts, allowing to fully resolve them within large, cosmological boxes. We follow the propagation of ionizing radiation from an arbitrary number of sources and from the recombination radiation produced by H and He. The chemical state of six species (HI, HII, HeI, HeII, HeIII, e) and gas temperatures are computed with a time-dependent, non-equilibrium chemistry solver. We present several validating tests of the code, including the standard tests from the RT Code Comparison Project and a new set of tests aimed at substantiating the new characteristics of RADAMESH. Using our AMR scheme, we show that properly resolving the I-front of a bright quasar during Reionization produces a large increase of the predicted gas temperature within the whole HII region. Also, we discuss how H and He recombination radiation is able to substantially change the ionization state of both species (for the classical Stroemgren sphere test) with respect to the widely used "on-the-spot" approximation.Comment: 19 pages, 24 figures; accepted for publication in MNRAS, version with high-resolution figures is avalaible at http://www.ast.cam.ac.uk/~cantal/Papers/CP10.pd

Diffuse continuum transfer in H II regions

We compare the accuracy of various methods for determining the transfer of the diffuse Lyman continuum in HII regions, by comparing them with a high-resolution discrete-ordinate integration. We use these results to suggest how, in multidimensional dynamical simulations, the diffuse field may be treated with acceptable accuracy without requiring detailed transport solutions. The angular distribution of the diffuse field derived from the numerical integration provides insight into the likely effects of the diffuse field for various material distributions.Comment: 12 pages, 7 figures, to be published in MNRA

Radiative transfer on hierarchial grids

We present new methods for radiative transfer on hierarchial grids. We develop a new method for calculating the scattered flux that employs the grid structure to speed up the computation. We describe a novel subiteration algorithm that can be used to accelerate calculations with strong dust temperature self-coupling. We compute two test models, a molecular cloud and a circumstellar disc, and compare the accuracy and speed of the new algorithms against existing methods. An adaptive model of the molecular cloud with less than 8 % of the cells in the uniform grid produced results in good agreement with the full resolution model. The relative RMS error of the surface brightness <4 % at all wavelengths, and in regions of high column density the relative RMS error was only 10^{-4}. Computation with the adaptive model was faster by a factor of ~5. The new method for calculating the scattered flux is faster by a factor of ~4 in large models with a deep hierarchy structure, when images of the scattered light are computed towards several observing directions. The efficiency of the subiteration algorithm is highly dependent on the details of the model. In the circumstellar disc test the speed-up was a factor of two, but much larger gains are possible. The algorithm is expected to be most beneficial in models where a large number of small, dense regions are embedded in an environment with a lower mean density.Comment: Accepted to A&A; 13 pages, 8 figures; (v2: minor typos corrected

LIME - a flexible, non-LTE line excitation and radiation transfer method for millimeter and far-infrared wavelengths

We present a new code for solving the molecular and atomic excitation and radiation transfer problem in a molecular gas and predicting emergent spectra. This code works in arbitrary three dimensional geometry using unstructured Delaunay latices for the transport of photons. Various physical models can be used as input, ranging from analytical descriptions over tabulated models to SPH simulations. To generate the Delaunay grid we sample the input model randomly, but weigh the sample probability with the molecular density and other parameters, and thereby we obtain an average grid point separation that scales with the local opacity. Our code does photon very efficiently so that the slow convergence of opaque models becomes traceable. When convergence between the level populations, the radiation field, and the point separation has been obtained, the grid is ray-traced to produced images that can readily be compared to observations. Because of the high dynamic range in scales that can be resolved using this type of grid, our code is particularly well suited for modeling of ALMA data. Our code can furthermore deal with overlapping lines of multiple molecular and atomic species.Comment: 13 pages, 12 figures, Accepted by A&A on 06/08/201

TRAPHIC - Radiative Transfer for Smoothed Particle Hydrodynamics Simulations

We present TRAPHIC, a novel radiative transfer scheme for Smoothed Particle Hydrodynamics (SPH) simulations. TRAPHIC is designed for use in simulations exhibiting a wide dynamic range in physical length scales and containing a large number of light sources. It is adaptive both in space and in angle and can be employed for application on distributed memory machines. The commonly encountered computationally expensive scaling with the number of light sources in the simulation is avoided by introducing a source merging procedure. The (time-dependent) radiative transfer equation is solved by tracing individual photon packets in an explicitly photon-conserving manner directly on the unstructured grid traced out by the set of SPH particles. To accomplish directed transport of radiation despite the irregular spatial distribution of the SPH particles, photons are guided inside cones. We present and test a parallel numerical implementation of TRAPHIC in the SPH code GADGET-2, specified for the transport of mono-chromatic hydrogen-ionizing radiation. The results of the tests are in excellent agreement with both analytic solutions and results obtained with other state-of-the-art radiative transfer codes.Comment: 31 pages, 20 figures. Accepted for publication in MNRAS. Revised version includes many clarifications and a new time-dependent radiative transfer calculation (fig. 19

The diffuse nature of Strömgren spheres

Wetensch. publicati