General Relativistic Radiative Transfer

We present a general method to calculate radiative transfer including scattering in the continuum as well as in lines in spherically symmetric systems that are influenced by the effects of general relativity (GR). We utilize a comoving wavelength ansatz that allows to resolve spectral lines throughout the atmosphere. The used numerical solution is an operator splitting (OS) technique that uses a characteristic formal solution. The bending of photon paths and the wavelength shifts due to the effects of GR are fully taken into account, as is the treatment of image generation in a curved spacetime. We describe the algorithm we use and demonstrate the effects of GR on the radiative transport of a two level atom line in a neutron star like atmosphere for various combinations of continuous and line scattering coefficients. In addition, we present grey continuum models and discuss the effects of different scattering albedos on the emergent spectra and the determination of effective temperatures and radii of neutron star atmospheres

Conservative Formulations of General Relativistic Radiative Transfer

Accurate accounting of particle number and 4-momentum in radiative transfer may be facilitated by the use of transport equations that allow transparent conversion between volume and surface integrals in both spacetime and momentum space. Such conservative formulations of general relativistic radiative transfer in multiple spatial dimensions are presented, and their relevance to core-collapse supernova simulations described.Comment: 4 pages. Talk presented at MG10, Rio de Janeiro, Brazil, 20-26 July 2003. To be published in Proceedings of the Tenth Marcel Grossmann Meeting on General Relativity, edited by M. Novello, S. Perez-Bergliaffa, and R. Ruffini, World Scientific, Singapore, 200

Radiative Transfer Along Rays in Curved Spacetimes

Radiative transfer in curved spacetimes has become increasingly important to understanding high-energy astrophysical phenomena and testing general relativity in the strong field limit. The equations of radiative transfer are physically equivalent to the Boltzmann equation, where the latter has the virtue of being covariant. We show that by a judicious choice of the basis of the phase space, it is generally possible to make the momentum derivatives in the Boltzmann equation vanish along an arbitrary (including nongeodesic) path, thus reducing the problem of radiative transfer along a ray to a path integral in coordinate space.Comment: To be published in MNRAS Letter

Radiative Heat Transfer in Anisotropic Many-Body Systems: Tuning and Enhancement

A general formalism for calculating the Radiative Heat Transfer in many body systems with anisotropic component is presented. Our scheme extends the theory of radiative heat transfer in isotropic many body systems to anisotropic cases. In addition, the radiative heating of the particles by the thermal bath is taken into account in our formula. It is shown that the radiative heat exchange (HE) between anisotropic particles and their radiative cooling/heating (RCH) could be enhanced several order of magnitude than that of isotropic particles. Furthermore, we demonstrate that both the HE and RCH can be tuned dramatically by particles relative orientation in many body systems

Optical depth in polarised Monte Carlo radiative transfer

Context: The Monte Carlo method is the most widely used method to solve radiative transfer problems in astronomy, especially in a fully general 3D geometry. A crucial concept in any Monte Carlo radiative transfer code is the random generation of the next interaction location. In polarised Monte Carlo radiative transfer with aligned non-spherical grains, the nature of dichroism complicates the concept of optical depth. Aims: We investigate, in detail, the relation between optical depth and the optical properties and density of the attenuating medium in polarised Monte Carlo radiative transfer codes that take dichroic extinction into account. Methods: Based on solutions for the radiative transfer equation, we discuss the optical depth scale in polarised radiative transfer with spheroidal grains. We compare the dichroic optical depth to the extinction and total optical depth scale. Results: In a dichroic medium, the optical depth is not equal to the usual extinction optical depth, nor to the total optical depth. For representative values of the optical properties of dust grains, the dichroic optical depth can differ from the extinction or total optical depth by several tens of percent. A closed expression for the dichroic optical depth cannot be given, but it can be derived efficiently through an algorithm that is based on the analytical result corresponding to elongated grains with a uniform grain alignment. Conclusions: Optical depth is more complex in dichroic media than in systems without dichroic attenuation, and this complexity needs to be considered when generating random free path lengths in Monte Carlo radiative transfer simulations. There is no benefit in using approximations instead of the dichroic optical depth

Covariant Compton Scattering Kernel in General Relativistic Radiative Transfer

A covariant scattering kernel is a core component in any self-consistent general relativistic radiative transfer formulation in scattering media. An explicit closed-form expression for a covariant Compton scattering kernel with a good dynamical energy range has unfortunately not been available thus far. Such an expression is essential to obtain numerical solutions to the general relativistic radiative transfer equations in complicated astrophysical settings where strong scattering effects are coupled with highly relativistic flows and steep gravitational gradients. Moreover, this must be performed in an efficient manner. With a self-consistent covariant approach, we have derived a closed-form expression for the Compton scattering kernel for arbitrary energy range. The scattering kernel and its angular moments are expressed in terms of hypergeometric functions, and their derivations are shown explicitly in this paper. We also evaluate the kernel and its moments numerically, assessing various techniques for their calculation. Finally, we demonstrate that our closed-form expression produces the same results as previous calculations, which employ fully numerical computation methods and are applicable only in more restrictive settings.Comment: 29 pages, 10 figures, 2 tables; Accepted for publication in MNRA

Mechanical relations between conductive and radiative heat transfer

We present a general nonequilibrium Green's function formalism for modeling heat transfer in systems characterized by linear response that establishes the formal algebraic relationships between phonon and radiative conduction, and reveals how upper bounds for the former can also be applied to the latter. We also propose an extension of this formalism to treat systems susceptible to the interplay of conductive and radiative heat transfer, which becomes relevant in atomic systems and at nanometric and smaller separations where theoretical descriptions which treat each phenomenon separately may be insufficient. We illustrate the need for such coupled descriptions by providing predictions for a low-dimensional system of carbyne wires in which the total heat transfer can differ from the sum of its radiative and conductive contributions. Our framework has ramifications for understanding heat transfer between large bodies that may approach direct contact with each other or that may be coupled by atomic, molecular, or interfacial film junctions.Comment: 16 pages, 2 figures, 1 table, 2 appendice

A class of Galerkin schemes for time-dependent radiative transfer

The numerical solution of time-dependent radiative transfer problems is challenging, both, due to the high dimension as well as the anisotropic structure of the underlying integro-partial differential equation. In this paper we propose a general framework for designing numerical methods for time-dependent radiative transfer based on a Galerkin discretization in space and angle combined with appropriate time stepping schemes. This allows us to systematically incorporate boundary conditions and to preserve basic properties like exponential stability and decay to equilibrium also on the discrete level. We present the basic a-priori error analysis and provide abstract error estimates that cover a wide class of methods. The starting point for our considerations is to rewrite the radiative transfer problem as a system of evolution equations which has a similar structure like first order hyperbolic systems in acoustics or electrodynamics. This analogy allows us to generalize the main arguments of the numerical analysis for such applications to the radiative transfer problem under investigation. We also discuss a particular discretization scheme based on a truncated spherical harmonic expansion in angle, a finite element discretization in space, and the implicit Euler method in time. The performance of the resulting mixed PN-finite element time stepping scheme is demonstrated by computational results