اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدCRASH_alpha: coupling continuum and line radiative transfer
31
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
In this paper we present CRASH_alpha, the first radiative transfer code for cosmological application that follows the parallel propagation of Ly_alpha and ionizing photons. CRASH_alpha is a version of the continuum radiative transfer code CRASH with a new algorithm to follow the propagation of Ly_alpha photons through a gas configuration whose ionization structure is evolving. The implementation introduces the time evolution for Ly_alpha photons (a feature commonly neglected in line radiative transfer codes) and, to reduce the computational time needed to follow each scattering, adopts a statistical approach to the Ly_alpha treatment by making extensive use of pre-compiled tables. With this statistical approach we experience a drastic increase of the computational speed and, at the same time, an excellent agreement with the full Ly_alpha radiative transfer computations of the code MCLy_alpha. We find that the emerging spectra keep memory of the ionization history which generates a given ionization configuration of the gas and, to properly account for this effect, a self-consistent joint evolution of line and ionizing continuum radiation as implemented in CRASH_alpha is necessary. A comparison between the results from our code and from Ly_alpha scattering alone on a fixed HI density field shows that the extent of the difference between the emerging spectra depends on the particular configuration considered, but it can be substantial and can thus affect the physical interpretation of the problem at hand. These differences should furthermore be taken into account when computing the impact of the Ly_alpha radiation on e.g. the observability of the 21 cm line from neutral hydrogen at epochs preceeding complete reionization.
قيم البحث
اقرأ أيضاً
Thermal dust emission carries information on physical conditions and dust properties in many astronomical sources. Because observations represent a sum of emission along the line of sight, their interpretation often requires radiative transfer modell
ing. We describe a new radiative transfer program SOC for computations of dust emission and examine its performance in simulations of interstellar clouds with external and internal heating. SOC implements the Monte Carlo radiative transfer method as a parallel program for shared memory computers. It can be used to study dust extinction, scattering, and emission. We tested SOC with realistic cloud models and examined the convergence and noise of the dust temperature estimates and of the resulting surface brightness maps. SOC has been demonstrated to produces accurate estimates for dust scattering and for thermal dust emission. It performs well with both with CPUs and with GPUs, the latter providing up to an order of magnitude speed-up. In the test cases, ALI improved the convergence rates but also was sensitive to Monte Carlo noise. Run-time refinement of the hierarchical-grid models did not help in reducing the run times required for a given accuracy of solution. The use of a reference field, without ALI, works more robustly. It also allows the run time to be optimised if the number of photon packages is increased only as the iterations progress. The use of GPUs in radiative transfer computations should be investigated further.
Radiative transfer modelling is part of many astrophysical simulations and is used to make synthetic observations and to assist analysis of observations. We concentrate on the modelling of the radio lines emitted by the interstellar medium. In connec
tion with high-resolution models, this can be significant computationally challenge. Our goal is a line radiative transfer (RT) program that makes good use of multi-core CPUs and GPUs. Parallelisation is essential to speed up computations and to enable the tackling of large modelling tasks with personal computers. The program LOC is based on ray-tracing and uses standard accelerated lambda iteration (ALI) methods for faster convergence. The program works on 1D and 3D grids. The 1D version makes use of symmetries to speed up the RT calculations. The 3D version works with octree grids and, to enable calculations with large models, is optimised for low memory usage. Tests show that LOC gives results that are in agreement with other RT codes to within ~2%. This is typical of code-to-code differences, which often are related to different interpretations of the model set-up. LOC run times compare favourably with those of Monte Carlo codes. In 1D tests, LOC runs were by up to a factor ~20 faster on a GPU than on a single CPU core. In spite of the complex path calculations, up to ~10 speed-up was observed also for 3D models using octree discretisation. GPUs enable calculations of models with hundreds of millions of cells, as encountered in the context of large-scale simulations of interstellar clouds. LOC shows good performance and accuracy and and is able to handle many RT modelling tasks on personal computers. Being written in Python, with the computing-intensive parts implemented as compiled OpenCL kernels, it can also a serve as a platform for further experimentation with alternative RT implementations.
Context. Magnetic fields are important to the dynamics of many astrophysical processes and can typically be studied through polarization observations. Polarimetric interferometry capabilities of modern (sub)millimeter telescope facilities have made i
t possible to obtain detailed velocity resolved maps of molecular line polarization. To properly analyze these for the information they carry regarding the magnetic field, the development of adaptive three-dimensional polarized line radiative transfer models is necessary. Aims. We aim to develop an easy-to-use program to simulate the polarization maps of molecular and atomic (sub)millimeter lines in magnetized astrophysical regions, such as protostellar disks, circumstellar envelopes, or molecular clouds. Methods. By considering the local anisotropy of the radiation field as the only alignment mechanism, we can model the alignment of molecular or atomic species inside a regular line radiative transfer simulation by only making use of the converged output of this simulation. Calculations of the aligned molecular or atomic states can subsequently be used to ray trace the polarized maps of the three-dimensional simulation. Results. We present a three-dimensional radiative transfer code, POlarized Radiative Transfer Adapted to Lines (PORTAL), that can simulate the emergence of polarization in line emission through a magnetic field of arbitrary morphology. Our model can be used in stand-alone mode, assuming LTE excitation, but it is best used when processing the output of regular three-dimensional (nonpolarized) line radiative transfer modeling codes. We present the spectral polarization map of test cases of a collapsing sphere and protoplanetary disk for multiple three-dimensional magnetic field morphologies.
We present a novel Lyman alpha (Ly$alpha$) radiative transfer code, SEURAT, where line scatterings are solved adaptively with the resolution of the smoothed particle hydrodynamics (SPH). The radiative transfer method implemented in SEURAT is based on
a Monte Carlo algorithm in which the scattering and absorption by dust are also incorporated. We perform standard test calculations to verify the validity of the code; (i) emergent spectra from a static uniform sphere, (ii) emergent spectra from an expanding uniform sphere, and (iii) escape fraction from a dusty slab. Thereby we demonstrate that our code solves the Ly$alpha$ radiative transfer with sufficient accuracy. We emphasise that SEURAT can treat the transfer of Ly$alpha$ photons even in highly complex systems that have significantly inhomogeneous density fields. The high adaptivity of SEURAT is desirable to solve the propagation of Ly$alpha$ photons in the interstellar medium of young star-forming galaxies like Ly$alpha$ emitters (LAEs). Thus, SEURAT provides a powerful tool to model the emergent spectra of Ly$alpha$ emission, which can be compared to the observations of LAEs.
Aims. We present MCFOST-art, a new non-local thermodynamic equilibrium radiative transfer solver for multilevel atomic systems. The code is embedded in the 3D radiative transfer code MCFOST and is compatible with most of the MCFOST modules. The code
is versatile and designed to model the close environment of stars in 3D. Methods. The code solves for the statistical equilibrium and radiative transfer equations using the Multilevel Accelerated Lambda Iteration (MALI) method. We tested MCFOST-art on spherically symmetric models of stellar photospheres as well as on a standard model of the solar atmosphere. We computed atomic level populations and outgoing fluxes and compared these values with the results of the TURBOspectrum and RH codes. Calculations including expansion and rotation of the atmosphere were also performed. We tested both the pure local thermodynamic equilibrium and the out-of-equilibrium problems. Results. In all cases, the results from all codes agree within a few percent at all wavelengths and reach the sub-percent level between RH and MCFOST-art. We still note a few marginal discrepancies between MCFOST-art and TURBOspectrum as a result of different treatments of background opacities at some critical wavelength ranges.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
