A novel approach for accurate radiative transfer in cosmological hydrodynamic simulations

• Published in: Monthly notices of the Royal Astronomical Society Vol. 415; no. 4; pp. 3731 - 3749
• Main Authors: Petkova, Margarita, Springel, Volker
• Format: Journal Article
• Language: English
• Published: Oxford, UK Blackwell Publishing Ltd 01.08.2011; Wiley-Blackwell; Oxford University Press
• Subjects: Astronomy; Astrophysics; Cosmology; Earth, ocean, space; Exact sciences and technology; intergalactic medium; methods: numerical; Radiation; radiative transfer; Simulation; methods: numerical; intergalactic medium; radiative transfer; Intergalactic matter; Second order; Boltzmann equation; Light velocity; Radiation flux; Numerical method; Spheres; Numerical approximation; Discretization; Dynamics; Operator splitting; Radiative transfer; Diffusion; Intensity; Digital simulation; Shading; Transfer equation; Point sources; Advection; Solid angle; Local field; Hydrodynamic model; Field ionization
• DOI: 10.1111/j.1365-2966.2011.18986.x
• ISSN: 0035-8711

We present a numerical implementation of radiative transfer based on an explicitly photon-conserving advection scheme, where radiative fluxes over the cell interfaces of a structured or unstructured mesh are calculated with a second-order reconstruction of the intensity field. The approach employs a direct discretization of the radiative transfer equation in Boltzmann form with adjustable angular resolution that, in principle, works equally well in the optically-thin and optically-thick regimes. In our most general formulation of the scheme, the local radiation field is decomposed into a linear sum of directional bins of equal solid angle, tessellating the unit sphere. Each of these 'cone fields' is transported independently, with constant intensity as a function of the direction within the cone. Photons propagate at the speed of light (or optionally using a reduced speed of light approximation to allow larger time-steps), yielding a fully time-dependent solution of the radiative transfer equation that can naturally cope with an arbitrary number of sources, as well as with scattering. The method casts sharp shadows, subject to the limitations induced by the adopted angular resolution. If the number of point sources is small and scattering is unimportant, our implementation can alternatively treat each source exactly in angular space, producing shadows whose sharpness is only limited by the grid resolution. A third hybrid alternative is to treat only a small number of the locally most luminous point sources explicitly, with the rest of the radiation intensity followed in a radiative diffusion approximation. We have implemented the method in the moving-mesh code arepo, where it is coupled to the hydrodynamics in an operator-splitting approach that subcycles the radiative transfer alternately with the hydrodynamical evolution steps. We also discuss our treatment of basic photon sink processes relevant to cosmological reionization, with a chemical network that can accurately deal with non-equilibrium effects. We discuss several tests of the new method, including shadowing configurations in two and three dimensions, ionized sphere expansion in static and dynamic density fields and the ionization of a cosmological density field. The tests agree favourably with analytical expectations and results based on other numerical radiative transfer approximations.