Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA radiative transfer scheme for cosmological reionization based on a local Eddington tensor
101
0
0.0
(
0
)
Authors
Dominique Aubert
Ask ChatGPT about the research
No Arabic abstract
A radiative transfer scheme is presented, based on a moment description of the equation of radiative transfer and the so-called ``M1 closure model for the Eddington tensor. This model features a strictly hyperbolic transport step for radiation: it has been implemented using standard Godunov--like techniques in a new code called ATON. Coupled to simple models of ionization chemistry and photo-heating, ATON is able to reproduce the results of other schemes on a various set of standard tests such as the expansion of a HII region, the shielding of the radiation by dense clumps and cosmological ionization by multiple sources. Being simple yet robust, such a scheme is intended to be naturally and easily included in grid--based cosmological fluid solvers.
rate research
Read More
We describe Rabacus, a Python package for calculating the transfer of hydrogen ionizing radiation in simplified geometries relevant to astronomy and cosmology. We present example solutions for three specific cases: 1) a semi-infinite slab gas distribution in a homogeneous isotropic background, 2) a spherically symmetric gas distribution with a point source at the center, and 3) a spherically symmetric gas distribution in a homogeneous isotropic background. All problems can accommodate arbitrary spectra and density profiles as input. The solutions include a treatment of both hydrogen and helium, a self-consistent calculation of equilibrium temperatures, and the transfer of recombination radiation. The core routines are written in Fortran 90 and then wrapped in Python leading to execution speeds thousands of times faster than equivalent routines written in pure Python. In addition, all variables have associated units for ease of analysis. The software is part of the Python Package Index and the source code is available on Bitbucket at https://bitbucket.org/galtay/rabacus . In addition, installation instructions and a detailed users guide are available at http://pythonhosted.org//rabacus .
We show that the algorithm based on the weighted essentially nonoscillatory (WENO) scheme with anti-diffusive flux corrections can be used as a solver of the radiative transfer equations. This algorithm is highly stable and robust for solving problems with both discontinuities and smooth solution structures. We test this code with the ionized sphere around point sources. It shows that the WENO scheme can reveal the discontinuity of the radiative or ionizing fronts as well as the evolution of photon frequency spectrum with high accuracy on coarse meshes and for a very wide parameter space. This method would be useful to study the details of the ionized patch given by individual source in the epoch of reionization. We demonstrate this method by calculating the evolution of the ionized sphere around point sources in physical and frequency spaces. It shows that the profile of the fraction of neutral hydrogen and the ionized radius are sensitively dependent on the intensity of the source.
A three dimensional parallel Monte Carlo (MC) dust radiative transfer code is presented. To overcome the huge computing time requirements of MC treatments, the computational power of vectorized hardware is used, utilizing either multi-core computer power or graphics processing units. The approach is a self-consistent way to solve the radiative transfer equation in arbitrary dust configurations. The code calculates the equilibrium temperatures of two populations of large grains and stochastic heated polycyclic aromatic hydrocarbons (PAH). Anisotropic scattering is treated applying the Heney-Greenstein phase function. The spectral energy distribution (SED) of the object is derived at low spatial resolution by a photon counting procedure and at high spatial resolution by a vectorized ray-tracer. The latter allows computation of high signal-to-noise images of the objects at any frequencies and arbitrary viewing angles. We test the robustness of our approach against other radiative transfer codes. The SED and dust temperatures of one and two dimensional benchmarks are reproduced at high precision. We utilize the Lucy-algorithm for the optical thin case where the Poisson noise is high, the iteration free Bjorkman & Wood method to reduce the calculation time, and the Fleck & Canfield diffusion approximation for extreme optical thick cells. The code is applied to model the appearance of active galactic nuclei (AGN) at optical and infrared wavelengths. The AGN torus is clumpy and includes fluffy composite grains of various sizes made-up of silicates and carbon. The dependence of the SED on the number of clumps in the torus and the viewing angle is studied. The appearance of the 10 micron silicate features in absorption or emission is discussed. The SED of the radio loud quasar 3C 249.1 is fit by the AGN model and a cirrus component to account for the far infrared emission.
The linear perturbation theory of inhomogeneous reionization (LPTR) has been developed as an analytical tool for predicting the global ionized fraction and large-scale power spectrum of ionized density fluctuations during reionization. In the original formulation of the LPTR, the ionization balance and radiative transfer equations are linearized and solved in Fourier space. However, the LPTRs approximation to the full solution of the radiative transfer equation is not straightforward to interpret, since the latter is most intuitively conceptualized in position space. To bridge the gap between the LPTR and the language of numerical radiative transfer, we present a new, equivalent, position-space formulation of the LPTR that clarifies the approximations it makes and facilitates its interpretation. We offer a comparison between the LPTR and the excursion-set model of reionization (ESMR), and demonstrate the built-in capability of the LPTR to explore a wide range of reionization scenarios, and to go beyond the ESMR in exploring scenarios involving X-rays.
One of the methodologies that carry out the division of the electrical grid into zones is based on the aggregation of nodes characterized by similar Power Transfer Distribution Factors (PTDFs). Here, we point out that satisfactory clustering algorithm should take into account two aspects. First, nodes of similar impact on cross-border lines should be grouped together. Second, cross-border power flows should be relatively insensitive to differences between real and assumed Generation Shift Key matrices. We introduce a theoretical basis of a novel clustering algorithm (BubbleClust) that fulfills these requirements and we perform a case study to illustrate social welfare consequences of the division.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
