Steps for solving the radiative transfer equation for arbitrary flows in stationary space–times

We derive the radiative transfer equation for arbitrary stationary relativistic flows in stationary space–times, i.e. for steady-state transfer problems. We show how the standard characteristics method of solution developed by Mihalas and used throughout the radiative transfer community can be adapted to multidimensional applications with isotropic sources. Because the characteristics always coincide with geodesics and can always be specified by constants, direct integration of the characteristics derived from the transfer equation as commonly done in 1D applications is not required. The characteristics are known for a specified metric from the geodesics. We give details in both flat and static spherically symmetric space–times. This work has direct application in 3D simulations of supernovae, gamma-ray bursts, and active galactic nuclei, as well as in modelling neutron star atmospheres.     

Albedo-Shifting Method: Scattering in a Line with Complete Frequency Redistribution

The albedo-shifting method is used to solve the problem of radiative transfer at line frequencies. Radiative transfer with complete frequency redistribution in a plane-parallel, semi-infinite atmosphere is considered. It is shown that the method works well in this case and enables one to considerably improve the convergence in an iterative solution of the equation for the source function.     

An exact linearization and decoupling of the integral equations satisfied by time-dependent X- and Y-functions

We discuss a simple method of linearization and decoupling of the integral equations satisfied by time-dependentX - andY -functions which play an important rôle in the study of non-stationary radiative transfer problems.     

A relaxation method for the computation of model stellar atmospheres

A general Monte Carlo relaxation method has been formulated for the computation of physically self-consistent model stellar atmospheres. The local physical state is obtained by solving simultaneously the equations of statistical equilibrium for the atomic and ionic level populations, the kinetic energy balance equation for the electron gas to obtain the electron temperature, and the equation of radiative transfer. Anisotropic Thomson scattering is included in the equation of transfer and radiation pressure effects are included in the hydrostatic equation. The constraints of hydrostatic and radiative equilibrium are enforced. Local thermodynamic equilibrium (L.T.E.) is assumed as a boundary condition deep in the atmosphere. Elsewhere in the atmosphere L.T.E. is not assumed.The statistical equilibrium equations are solved with no assumptions made concerning detailed balance for the bound-bound radiative processes. The source function is formulated in microscopic detail. All atomic processes contributing to the absorption and emission of radiation are included. The kinetic energy balance equation for the electron gas is formulated in detail. All atomic processes by which kinetic energy is gained and lost by the electron gas are included.The method has been applied to the computation of a model atmosphere for a pure hydrogen early-type star. An idealized model of the hydrogen atom with five bound levels and the continuum was adopted. The results of the trial calculation are discussed with reference to stability, accuracy, and convergence of the solution.Contribution No. 385 from the Kitt Peak National Observatory.Operated by the Association of Universities for Research in Astronomy, Inc., under contract with the National Science Foundation.     

The polarization-free approximation applied to multi-level non-LTE radiative transfer

The polarization-free (POF) approximation (Trujillo Bueno and Landi Degl'Innocenti, 1996) is capable of accounting for the approximate influence of the magnetic field on the statistical equilibrium, without actually solving the full Stokes vector radiative transfer equation. The method introduces the Zeeman splitting or broadening of the line absorption profile I in the scalar radiative transfer equation, but the coupling between Stokes I and the other Stokes parameters is neglected. The expected influence of the magnetic field is largest for strongly-split strong lines and the effect is greatly enhanced by gradients in the magnetic field strength. Formally the interaction with the other Stokes parameters may not be neglected for strongly-split strong lines, but it turns out that the error in Stokes I obtained through the POF approximation to a large extent cancels the neglect of interaction with the other Stokes parameters, so that the resulting line source functions and line opacities are more accurate than those obtained with the field-free approach. Although its merits have so far only been tested for a two-level atom, we apply the POF approximation to multi-level non-LTE radiative transfer problems on the premise that there is no essential difference between these two cases. Final verification of its validity in multi-level cases still awaits the completion of a non-LTE Stokes vector transfer code.For two realistic multi-level cases (CaII and MgI in the solar atmosphere) it is demonstrated that the POF method leads to small changes, with respect to the field-free method, in the line source functions and emergent Stokes vector profiles (much smaller than for a two-level atom). Real atoms are dominated by strong ultraviolet lines (only weakly split) and continua, and most lines with large magnetic splitting (in the red and the infrared) are at higher excitation energies, i.e. they are relatively weak and unable to produce significant changes in the statistical equilibrium. We find that it is generally unpredictable by how much the POF results will differ from the field-free results, so that it is nearly always necessary to confirm predictions by actual computations.The POF approximation provides more reliable results than the field-free approximation without significantly complicating the radiative transfer problem, i.e. without solving any extra equations and without excessive computational resource requirements, so that it is to be preferred over the field-free approximation.     

交互投影迭代算法的原理与实践：用线性化分离方法求解Non—LTE恒星大气模型

讨论交互投影的迭代算法在线性化分离求解Non-LTE恒星大气模型中的运用和实践，交互投影的迭代算法可以实现降维和病态分离，它的解可以收敛至原来解的最小二乘解，实践证实了这一点，因此线性化分离求解大气模型的方法可以代替完全线性化方法。     

Radiative heat transfer to hydromagnetic flow of a slightly rarefied binary gas in a vertical channel

The paper considers the fully-developed slip flow in a vertical channel with radiative heat transfer and mass transfer in the presence of an externally applied magnetic field. The problem is modelled by the compressible Navier-Stokes equations, so that the gas is only slightly rarefied. Invoking the exact integral equation for radiation, the problem is reduced to a set of ordinary integro-differential equations. By realistic assumptions, the set is linearized and the temperature is reduced to a mixed Fredholm-Volterra integral equation which is solved by standard iterative procedure. Thereafter the concentration equation is solved by the WKB approximation while the velocity is obtained by the finite difference scheme. These solutions are discussed qualitatively.     

A discrete spherical harmonics method for radiative transfer analysis in inhomogeneous polarized planar atmosphere

A discrete spherical harmonics method is developed for the radiative transfer problem in inhomogeneous polarized planar atmosphere illuminated at the top by a collimated sunlight while the bottom reflects the radiation. The method expands both the Stokes vector and the phase matrix in a finite series of generalized spherical functions and the resulting vector radiative transfer equation is expressed in a set of polar directions. Hence, the polarized characteristics of the radiance within the atmosphere at any polar direction and azimuthal angle can be determined without linearization and/or interpolations. The spatial dependent of the problem is solved using the spectral Chebyshev method. The emergent and transmitted radiative intensity and the degree of polarization are predicted for both Rayleigh and Mie scattering. The discrete spherical harmonics method predictions for optical thin atmosphere using 36 streams are found in good agreement with benchmark literature results. The maximum deviation between the proposed method and literature results and for polar directions \(\vert \mu \vert \geq0.1 \) is less than 0.5% and 0.9% for the Rayleigh and Mie scattering, respectively. These deviations for directions close to zero are about 3% and 10% for Rayleigh and Mie scattering, respectively.     

Angular distribution of radiation in an isotropically scattering slab

The equation of radiative transfer in an isotropically scattering slab subject to general boundary conditions is solved. The Padé approximation technique is used to calculate the reflected and transmitted angular distributions. Numerical results for angular distributions through and at the boundaries of a finite slab are calculated by the Padé approximation technique. The results for a Padé approximation of [0/1] are compared with results obtained by the Galerkin method.     

Solution of a Riemann–Hilbert problem for a new expression of Chandrasekhar’s H-function in radiative transfer

The linear singular integral equation derived from the nonlinear integral equation of Chandrasekhar’s H-function in radiative transfer is considered here to develop a new form of H-function as a solution of a Riemann–Hilbert problem using Plemelj and Cauchy integral formulae for complex domain. This new form of H-function is a simple integral of known functions. Forms of H-function both for conservative and nonconservative cases are obtained. Their numerical evaluations are made by Simpson’s one-third rule to arrive at an accuracy to ninth places of decimals.     

The solution of a time-dependent equation of transfer in nonconservative finite media

A time-dependent nonconservative radiative transfer equation in a finite media has been solved. Firstly, the time-dependent function is converted into a time-independent function by the use of a Laplace transform and the transformed equation is then solved by a method using linear operator theory.     

Radiation hydrodynamics using characteristics on adaptive decomposed domains for massively parallel star formation simulations

We present an algorithm for solving the radiative transfer problem on massively parallel computers using adaptive mesh refinement and domain decomposition. The solver is based on the method of characteristics which requires an adaptive raytracer that integrates the equation of radiative transfer. The radiation field is split into local and global components which are handled separately to overcome the non-locality problem. The solver is implemented in the framework of the magneto-hydrodynamics code FLASH and is coupled by an operator splitting step. The goal is the study of radiation in the context of star formation simulations with a focus on early disc formation and evolution. This requires a proper treatment of radiation physics that covers both the optically thin as well as the optically thick regimes and the transition region in particular. We successfully show the accuracy and feasibility of our method in a series of standard radiative transfer problems and two 3D collapse simulations resembling the early stages of protostar and disc formation.     

Exact solution of a key equation in a finite stellar atmosphere by the method of Laplace transform and linear singular operators

The key equation which commonly appears for radiative transfer in a finite stellar atmosphere having ground reflection according to Lambert's law is considered in this paper. The exact solution of this equation is obtained for surface quantities in terms of theX-Y equations of Chandrasekhar by the method of Laplace transform and linear singular operators. This exact method is widely applicable for obtaining the solution for surface quantities in a finite atmosphere.     

Solution of the equation of transfer with Rayleigh's phase function in a thin atmosphere

The equation of radiative transfer with scattering according to Rayleigh's phase function has been solved in a thin atmosphere by use of a modification of the spherical-harmonic method suggested by Wanet al. (1986).     

Spectra and light curves of GRB afterglows

We performed accurate numerical calculations of angle-, time-, and frequency-dependent radiative transfer for the relativistic motion of matter in gamma-ray burst (GRB) models. Our technique for solving the transfer equation, which is based on the method of characteristics, can be applied to the motion of matter with a Lorentz factor up to 1000. The effect of synchrotron self-absorption is taken into account. We computed the spectra and light curves from electrons with a power-law energy distribution in an expanding relativistic shock and compare them with available analytic estimates. The behavior of the optical afterglows from GRB 990510 and GRB 000301c is discussed qualitatively.     

On nonstationary radiative transfer in a spectral line in stellar atmospheres

We consider nonstationary radiative transfer in a line in stellar atmospheres simulated as a stationary semi-infinite plane-parallel medium. We assume complete frequency redistribution in the elementary act of scattering. We assume that the time a photon spends in the medium is determined only by the mean time spent in the absorbed state. We obtain an explicit expression for the resolvent of the nonstationary integral equation of transfer, which is a bilinear expansion with respect to the eigenfunctions found in [12] for the corresponding stationary transfer equation.Translated fromAstrofizika, Vol. 37, No. 3, 1994.I am grateful to the American Astronomical Society for financial support of this work.     

Polarization in spectral lines

The general formalism, presented in a previous paper of this series (Landi Degl'Innocenti, 1983a), is particularized to deduce the radiative transfer equations for polarized radiation and the statistical equilibrium equations for a multi-level atom in the zero-magnetic field, collisionless regime. The formulae are developed both in the standard representation and in the representation of the statistical tensors. For resonance scattering in a two-level atom, in the limiting case of complete depolarization of the ground level, we recover the classical results for Rayleigh scattering and we derive the expression of the phase matrix in terms of ordinary rotation matrices. The law of scattering is then generalized to take properly into account the influence of the anisotropy of the radiation field on the atomic polarization of the ground level (depopulation pumping).     

Radiative transfer

Chandrasekhar’s work in radiative transfer theory began in 1944 and culminated with the publication of his influential treatiseRadiative Transfer in 1950. In this review his major contributions to radiative transfer will be recounted and evaluated. These include his development of the discrete ordinates method, the invariance principles, and his formulation and solution of the transfer equation for polarized light.     

A new moment method for continuum radiative transfer in cosmological re-ionization

We introduce a new code for computing time-dependent continuum radiative transfer and non-equilibrium ionization states in static density fields with periodic boundaries. Our code solves the moments of the radiative transfer equation, closed by an Eddington tensor computed using a long characteristics (LC) method. We show that traditional short characteristics and the optically thin approximation are inappropriate for computing Eddington factors for the problem of cosmological re-ionization. We evolve the non-equilibrium ionization field via an efficient and accurate (errors <1 per cent) technique that switches between fully implicit or explicit finite differencing depending on whether the local time-scales are long or short compared to the time-step. We tailor our code for the problem of cosmological re-ionization. In tests, the code conserves photons, accurately treats cosmological effects and reproduces analytic Strömgren sphere solutions. Its chief weakness is that the computation time for the LC calculation scales relatively poorly compared to other techniques  ( t _LC∝ N ^∼1.5_cells)  ; however, we mitigate this by only recomputing the Eddington tensor when the radiation field changes substantially. Our technique makes almost no physical approximations, so it provides a way to benchmark faster but more approximate techniques. It can readily be extended to evolve multiple frequencies, though we do not do so here. Finally, we note that our method is generally applicable to any problem involving the transfer of continuum radiation through a periodic volume.