Self-Consistent Thermal Accretion Disk Corona Models for Compact Objects: I. Properties of the Corona and the Spectrum of Escaping Radiation

We present the properties of accretion disk corona (ADC) models, where the radiation field, the temperature, and the total opacity of the corona are determined self-consistently. We use a non-linear Monte Carlo code to perform the calculations. As an example, we discuss models where the corona is situated above and below a cold accretion disk with a plane-parallel (slab) geometry, similar to the model of Haardt and Maraschi. By Comptonizing the soft radiation emitted by the accretion disk, the corona is responsible for producing the high-energy component of the escaping radiation. Our models include the reprocessing of radiation in the accretion disk. Here, the photons either are Compton reflected or photo-absorbed, giving rise to fluorescent line emission and thermal emission. The self-consistent coronal temperature is determined by balancing heating (due to viscous energy dissipation) with Compton cooling, determined using the fully relativistic, angle-dependent cross-sections. The total opacity is found by balancing pair productions with annihilations. We find that, for a disk temperature kT_bb \lta 200 eV, these coronae are unable to have a self-consistent temperature higher than \sim 120 keV if the total optical depth is \gta 0.2, regardless of the compactness parameter of the corona and the seed opacity. This limitation corresponds to the angle-averaged spectrum of escaping radiation having a photon index \gta 1.8 within the 5 keV - 30 keV band. Finally, all models that have reprocessing features also predict a large thermal excess at lower energies. These constraints make explaining the X-ray spectra of persistent black hole candidates with ADC models very problematic.Comment: 15 pages, Latex, 9 .eps figures, uses emulateapj.sty (included). To be published in ApJ, October 1, 1997, Vol. 48

Radiative transfer in very optically thick circumstellar disks

In this paper we present two efficient implementations of the diffusion approximation to be employed in Monte Carlo computations of radiative transfer in dusty media of massive circumstellar disks. The aim is to improve the accuracy of the computed temperature structure and to decrease the computation time. The accuracy, efficiency and applicability of the methods in various corners of parameter space are investigated. The effects of using these methods on the vertical structure of the circumstellar disk as obtained from hydrostatic equilibrium computations are also addressed. Two methods are presented. First, an energy diffusion approximation is used to improve the accuracy of the temperature structure in highly obscured regions of the disk, where photon counts are low. Second, a modified random walk approximation is employed to decrease the computation time. This modified random walk ensures that the photons that end up in the high-density regions can quickly escape to the lower density regions, while the energy deposited by these photons in the disk is still computed accurately. A new radiative transfer code, MCMax, is presented in which both these diffusion approximations are implemented. These can be used simultaneously to increase both computational speed and decrease statistical noise. We conclude that the diffusion approximations allow for fast and accurate computations of the temperature structure, vertical disk structure and observables of very optically thick circumstellar disks.Comment: Accepted for publication in A&

An algorithm for Monte-Carlo time-dependent radiation transfer

A new Monte-Carlo algorithm for calculating time-dependent radiative-transfer under the assumption of LTE is presented. Unlike flux-limited diffusion the method is polychromatic, includes scattering, and is able to treat the optically thick and free-streaming regimes simultaneously. The algorithm is tested on a variety of 1-d and 2-d problems, and good agreement with benchmark solutions is found. The method is used to calculate the time-varying spectral energy distribution from a circumstellar disc illuminated by a protostar whose accretion luminosity is varying. It is shown that the time lag between the optical variability and the infrared variability results from a combination of the photon travel time and the thermal response in the disc, and that the lag is an approximately linear function of wavelength.Comment: 10 pages, 10 figures, accepted for publication by MNRA

RADAMESH: Cosmological Radiative Transfer for Adaptive Mesh Refinement Simulations

We present a new three-dimensional radiative transfer (RT) code, RADAMESH, based on a ray-tracing, photon-conserving and adaptive (in space and time) scheme. RADAMESH uses a novel Monte Carlo approach to sample the radiation field within the computational domain on a "cell-by-cell" basis. Thanks to this algorithm, the computational efforts are now focused where actually needed, i.e. within the Ionization-fronts (I-fronts). This results in an increased accuracy level and, at the same time, a huge gain in computational speed with respect to a "classical" Monte Carlo RT, especially when combined with an Adaptive Mesh Refinement (AMR) scheme. Among several new features, RADAMESH is able to adaptively refine the computational mesh in correspondence of the I-fronts, allowing to fully resolve them within large, cosmological boxes. We follow the propagation of ionizing radiation from an arbitrary number of sources and from the recombination radiation produced by H and He. The chemical state of six species (HI, HII, HeI, HeII, HeIII, e) and gas temperatures are computed with a time-dependent, non-equilibrium chemistry solver. We present several validating tests of the code, including the standard tests from the RT Code Comparison Project and a new set of tests aimed at substantiating the new characteristics of RADAMESH. Using our AMR scheme, we show that properly resolving the I-front of a bright quasar during Reionization produces a large increase of the predicted gas temperature within the whole HII region. Also, we discuss how H and He recombination radiation is able to substantially change the ionization state of both species (for the classical Stroemgren sphere test) with respect to the widely used "on-the-spot" approximation.Comment: 19 pages, 24 figures; accepted for publication in MNRAS, version with high-resolution figures is avalaible at http://www.ast.cam.ac.uk/~cantal/Papers/CP10.pd

Fast and accurate frequency-dependent radiation transport for hydrodynamics simulations in massive star formation

Context: Radiative feedback plays a crucial role in the formation of massive stars. The implementation of a fast and accurate description of the proceeding thermodynamics in pre-stellar cores and evolving accretion disks is therefore a main effort in current hydrodynamics simulations. Aims: We introduce our newly implemented three-dimensional frequency dependent radiation transport algorithm for hydrodynamics simulations of spatial configurations with a dominant central source. Methods: The module combines the advantage of the speed of an approximate Flux Limited Diffusion (FLD) solver with the high accuracy of a frequency dependent first order ray-tracing routine. Results: We prove the viability of the scheme in a standard radiation benchmark test compared to a full frequency dependent Monte-Carlo based radiative transfer code. The setup includes a central star, a circumstellar flared disk, as well as an envelope. The test is performed for different optical depths. Considering the frequency dependence of the stellar irradiation, the temperature distributions can be described precisely in the optically thin, thick, and irradiated transition regions. Resulting radiative forces onto dust grains are reproduced with high accuracy. The achievable parallel speedup of the method imposes no restriction on further radiative (magneto-) hydrodynamics simulations. Conclusions: The proposed approximate radiation transport method enables frequency dependent radiation hydrodynamics studies of the evolution of pre-stellar cores and circumstellar accretion disks around an evolving massive star in a highly efficient and accurate manner.Comment: 16 pages, 11 figure

Monte Carlo Radiative Transfer

I outline methods for calculating the solution of Monte Carlo Radiative Transfer (MCRT) in scattering, absorption and emission processes of dust and gas, including polarization. I provide a bibliography of relevant papers on methods with astrophysical applications.Comment: To appear in the Chandra Centennial issue of the Bulletin of the Astronomical Society of India, volume 39 (2011), eds D.J. Saikia and Virginia Trimble; 27 pages, 1 figur

HYPERION: An open-source parallelized three-dimensional dust continuum radiative transfer code

HYPERION is a new three-dimensional dust continuum Monte-Carlo radiative transfer code that is designed to be as generic as possible, allowing radiative transfer to be computed through a variety of three-dimensional grids. The main part of the code is problem-independent, and only requires an arbitrary three-dimensional density structure, dust properties, the position and properties of the illuminating sources, and parameters controlling the running and output of the code. HYPERION is parallelized, and is shown to scale well to thousands of processes. Two common benchmark models for protoplanetary disks were computed, and the results are found to be in excellent agreement with those from other codes. Finally, to demonstrate the capabilities of the code, dust temperatures, SEDs, and synthetic multi-wavelength images were computed for a dynamical simulation of a low-mass star formation region. HYPERION is being actively developed to include new features, and is publicly available (http://www.hyperion-rt.org).Comment: Accepted for publication in Astronomy & Astrophysics. HYPERION is being prepared for release at the start of 2012, but you can already sign up to the mailing list at http://www.hyperion-rt.org to be informed once it is available for downloa

Monte-Carlo methods for NLTE spectral synthesis of supernovae

We present JEKYLL, a new code for modelling of supernova (SN) spectra and lightcurves based on Monte-Carlo (MC) techniques for the radiative transfer. The code assumes spherical symmetry, homologous expansion and steady state for the matter, but is otherwise capable of solving the time-dependent radiative transfer problem in non-local-thermodynamic-equilibrium (NLTE). The method used was introduced in a series of papers by Lucy, but the full time-dependent NLTE capabilities of it have never been tested. Here, we have extended the method to include non-thermal excitation and ionization as well as charge-transfer and two-photon processes. Based on earlier work, the non-thermal rates are calculated by solving the Spencer-Fano equation. Using a method previously developed for the SUMO code, macroscopic mixing of the material is taken into account in a statistical sense. In addition, a statistical Markov-chain model is used to sample the emission frequency, and we introduce a method to control the sampling of the radiation field. Except for a description of JEKYLL, we provide comparisons with the ARTIS, SUMO and CMFGEN codes, which show good agreement in the calculated spectra as well as the state of the gas. In particular, the comparison with CMFGEN, which is similar in terms of physics but uses a different technique, shows that the Lucy method does indeed converge in the time-dependent NLTE case. Finally, as an example of the time-dependent NLTE capabilities of JEKYLL, we present a model of a Type IIb SN, taken from a set of models presented and discussed in detail in an accompanying paper. Based on this model we investigate the effects of NLTE, in particular those arising from non-thermal excitation and ionization, and find strong effects even on the bolometric lightcurve. This highlights the need for full NLTE calculations when simulating the spectra and lightcurves of SNe.Comment: Accepted for publication by Astronomy & Astrophysic