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.

We have developed a three-dimensional radiative transfer method designed specifically for use with parallel adaptive mesh refinement hydrodynamics codes. This new algorithm, which we call hybrid characteristics, introduces a novel form of ray tracing that can neither be classified as long, nor as short characteristics, but which applies the underlying principles, i.e. efficient execution through interpolation and parallelizability, of both. Primary applications of the hybrid characteristics method are radiation hydrodynamics problems that take into account the effects of photoionization and heating due to point sources of radiation. The method is implemented in the hydrodynamics package FLASH. The ionization, heating, and cooling processes are modelled using the DORIC ionization package. Upon comparison with the long characteristics method, we find that our method calculates the column density with a similarly high accuracy and produces sharp and well defined shadows. We show the quality of the new algorithm in an application to the photoevaporation of multiple over-dense clumps. We present several test problems demonstrating the feasibility of our method for performing high resolution three-dimensional radiation hydrodynamics calculations that span a large range of scales. Initial performance tests show that the ray tracing part of our method takes less time to execute than other parts of the calculation (e.g. hydrodynamics and adaptive mesh refinement), and that a high degree of efficiency is obtained in parallel execution. Although the hybrid characteristics method is developed for problems involving photoionization due to point sources, the algorithm can be easily adapted to the case of more general radiation fields.

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.

At the end of their lives low mass stars such as our Sun lose most of their mass. The resulting planetary nebulae show a wide variety of shapes, from spherical to highly bipolar. According to the generalized interacting stellar winds model, these shapes are due to an interaction between a very fast tenuous outflow, and a denser environment left over from an earlier slow phase of mass loss. Previous analytical and numerical work shows that this mechanism can explain cylindrically symmetric nebulae very well. However, many circumstellar nebulae have a multipolar or point-symmetric shape. With two-dimensional calculations, Icke showed that these seemingly enigmatic forms can be easily reproduced by a two-wind model in which the confining disk is warped, as is expected to occur in irradiated disks. Here, we present the extension to fully three-dimensional adaptive mesh refinement simulations of such an interaction.

Previous analytical and numerical work shows that the generalized interacting stellar winds model can explain the observed bipolar shapes of planetary nebulae very well. However, many circumstellar nebulae have a multipolar or point-symmetric shape. With two-dimensional calculations, Icke showed that these seemingly enigmatic forms can be easily reproduced by a two-wind model in which the confining disk is warped, as is expected to occur in irradiated disks. In this contribution we present the extension to fully three-dimensional adaptive mesh refinement simulations of such an interaction.

We have discovered a 7.06+-0.08 Hz quasi-periodic oscillation (QPO) in the X-ray flux of the low-luminosity low-mass X-ray binary (LMXB) and atoll source 4U 1820-30. This QPO was only observable at the highest observed mass accretion rate, when the source was in the uppermost part of the banana branch, at a 2-25 keV luminosity of 5.4x10^37 erg/s (for a distance of 6.4 kpc). The QPO had a FWHM of only 0.5+-0.2 Hz during small time intervals (32-s of data), and showed erratic shifts in the centroid frequency between 5.5 and 8 Hz. The rms amplitude over the energy range 2-60 keV was 5.6%+-0.2%. The amplitude increased with photon energy from 3.7%+-0.5% between 2.8 and 5.3 keV to 7.3%+-0.6% between 6.8 and 9.3 keV, above which it remained approximately constant at ~7%. The time lag of the QPO between 2.8-6.8 and 6.8-18.2 keV was consistent with being zero (-1.2+-3.4 ms). The properties of the QPO (i.e., its frequency and its presence only at the highest observed mass accretion rate) are similar to those of the 5-20 Hz QPO observed in the highest luminosity LMXBs (the Z sources) when they are accreting near the Eddington mass accretion limit. If this is indeed the same phenomenon, then models explaining the 5-20 Hz QPO in the Z sources, which require the near-Eddington accretion rates, will not hold. Assuming isotropic emission, the 2-25 keV luminosity of 4U 1820-30 at the time of the 7 Hz QPOs is at maximum only 40% (for a companion star with cosmic abundances), but most likely ~20% (for a helium companion star) of the Eddington accretion limit.