Enhancement of optically thick to thin line intensities in solar and stellar coronal plasmas through radiative transfer effects: an angularly resolved study.

An analysis of radiative transfer effects present in the Fe XV ion stage of solar and stellar coronal plasmas provides a general explanation of line radiation intensity enhancement above the optically thin limit. Full linearization radiation transfer is compared with the escape factor method and found to be in good agreement at the lower column densities. An angular study of the enhancement shows that symmetry factors are of great importance. This gives a possible reason for the indeterminate status of opacity in relation to coronal lines of distant stellar sources, where only emission integrated across the whole surface is detected.

Multidimensional simulations of radiative transfer in Type Ia supernovae

A three-dimensional Monte Carlo code for modelling radiation transport in Type Ia supernovae is described. In addition to tracking Monte Carlo quanta to follow the emission, scattering and deposition of radiative energy, a scheme involving volume-based Monte Carlo estimators is used to allow properties of the emergent radiation field to be extracted for specific viewing angles in a multidimensional structure. This eliminates the need to compute spectra or light curves by angular binning of emergent quanta. The code is applied to two test problems to illustrate consequences of multidimensional structure on the modelling of light curves. First, elliptical models are used to quantify how large-scale asphericity can introduce angular dependence to light curves. Secondly, a model which incorporates complex structural inhomogeneity, as predicted by modern explosion models, is used to investigate how such structure may affect light-curve properties. © 2006 RAS.

Line-driven Disk Winds in Activ Galactic Nuclei: The Critical Importance of Ionization and Radiative Transfer

Accretion disk winds are thought to produce many of the characteristic features seen in the spectra of active galactic nuclei (AGNs) and quasi-stellar objects (QSOs). These outflows also represent a natural form of feedback between the central supermassive black hole and its host galaxy. The mechanism for driving this mass loss remains unknown, although radiation pressure mediated by spectral lines is a leading candidate. Here, we calculate the ionization state of, and emergent spectra for, the hydrodynamic simulation of a line-driven disk wind previously presented by Proga & Kallman. To achieve this, we carry out a comprehensive Monte Carlo simulation of the radiative transfer through, and energy exchange within, the predicted outflow. We find that the wind is much more ionized than originally estimated. This is in part because it is much more difficult to shield any wind regions effectively when the outflow itself is allowed to reprocess and redirect ionizing photons. As a result, the calculated spectrum that would be observed from this particular outflow solution would not contain the ultraviolet spectral lines that are observed in many AGN/QSOs. Furthermore, the wind is so highly ionized that line driving would not actually be efficient. This does not necessarily mean that line-driven winds are not viable. However, our work does illustrate that in order to arrive at a self-consistent model of line-driven disk winds in AGN/QSO, it will be critical to include a more detailed treatment of radiative transfer and ionization in the next generation of hydrodynamic simulations.

Testing quasar unification: radiative transfer in clumpy winds

Various unification schemes interpret the complex phenomenology of quasars and luminous active galactic nuclei (AGN) in terms of a simple picture involving a central black hole, an accretion disc and an associated outflow. Here, we continue our tests of this paradigm by comparing quasar spectra to synthetic spectra of biconical disc wind models, produced with our state-of-the-art Monte Carlo radiative transfer code. Previously, we have shown that we could produce synthetic spectra resembling those of observed broad absorption line (BAL) quasars, but only if the X-ray luminosity was limited to 1043 erg s-1. Here, we introduce a simple treatment of clumping, and find that a filling factor of ˜0.01 moderates the ionization state sufficiently for BAL features to form in the rest-frame UV at more realistic X-ray luminosities. Our fiducial model shows good agreement with AGN X-ray properties and the wind produces strong line emission in, e.g., Lyα and C IV 1550 Å at low inclinations. At high inclinations, the spectra possess prominent LoBAL features. Despite these successes, we cannot reproduce all emission lines seen in quasar spectra with the correct equivalent-width ratios, and we find an angular dependence of emission line equivalent width despite the similarities in the observed emission line properties of BAL and non-BAL quasars. Overall, our work suggests that biconical winds can reproduce much of the qualitative behaviour expected from a unified model, but we cannot yet provide quantitative matches with quasar properties at all viewing angles. Whether disc winds can successfully unify quasars is therefore still an open question.

Scaling of magneto-quantum-radiative hydrodynamic equations: From laser-produced plasmas to astrophysics

We introduce the equations of magneto-quantum-radiative hydrodynamics. By rewriting them in a dimensionless form, we obtain a set of parameters that describe scale-dependent ratios of characteristic hydrodynamic quantities. We discuss how these dimensionless parameters relate to the scaling between astrophysical observations and laboratory experiments.

Radiative hydrodynamic modelling and observations of the X-class solar flare on 2011 March 9

Aims. We investigated the response of the solar atmosphere to non-thermal electron beam heating using the radiative transfer and hydrodynamics modelling code RADYN. The temporal evolution of the parameters that describe the non-thermal electron energy distribution were derived from hard X-ray observations of a particular flare, and we compared the modelled and observed parameters.

Methods. The evolution of the non-thermal electron beam parameters during the X1.5 solar flare on 2011 March 9 were obtained from analysis of RHESSI X-ray spectra. The RADYN flare model was allowed to evolve for 110 s, after which the electron beam heating was ended, and was then allowed to continue evolving for a further 300 s. The modelled flare parameters were compared to the observed parameters determined from extreme-ultraviolet spectroscopy.

Results. The model produced a hotter and denser flare loop than that observed and also cooled more rapidly, suggesting that additional energy input in the decay phase of the flare is required. In the explosive evaporation phase a region of high-density cool material propagated upward through the corona. This material underwent a rapid increase in temperature as it was unable to radiate away all of the energy deposited across it by the non-thermal electron beam and via thermal conduction. A narrow and high-density (ne ≤ 1015 cm-3) region at the base of the flare transition region was the source of optical line emission in the model atmosphere. The collision-stopping depth of electrons was calculated throughout the evolution of the flare, and it was found that the compression of the lower atmosphere may permit electrons to penetrate farther into a flaring atmosphere compared to a quiet Sun atmosphere.

The Physical and Chemical Structure of Hot Molecular Cores

We have made self-consistent models of the density and temperature profiles of the gas and dust surrounding embedded luminous objects using a detailed radiative transfer model together with observations of the spectral energy distribution of hot molecular cores. Using these profiles we have investigated the hot core chemistry which results when grain mantles are evaporated, taking into account the different binding energies of the mantle molecules, as well a model in which we assume that all molecules are embedded in water ice and have a common binding energy. We find that most of the resulting column densities are consistent with those observed toward the hot core G34.3+0.15 at a time around 10^4 years after central luminous star formation. We have also investigated the dependence of the chemical structure on the density profile which suggests an observational possibility of constraining density profiles from determination of the source sizes of line emission from desorbed molecules.

Chemical processes in protoplanetary disks. II. On the importance of photochemistry and X-ray ionization

We investigate the impact of photochemistry and X-ray ionization on the molecular composition of, and ionization fraction in, a protoplanetary disk surrounding a typical T Tauri star. We use a sophisticated physical model, which includes a robust treatment of the radiative transfer of UV and X-ray radiation, and calculate the time-dependent chemical structure using a comprehensive chemical network. In previous work, we approximated the photochemistry and X-ray ionization; here, we recalculate the photoreaction rates using the explicit UV wavelength spectrum and wavelength-dependent reaction cross sections. We recalculate the X-ray ionization rate using our explicit elemental composition and X-ray energy spectrum. We find that photochemistry has a larger influence on the molecular composition than X-ray ionization. Observable molecules sensitive to the photorates include OH, HCO+, N2H+, H2O, CO2, and CH3OH. The only molecule significantly affected by the X-ray ionization is N2H+, indicating that it is safe to adopt existing approximations of the X-ray ionization rate in typical T Tauri star-disk systems. The recalculation of the photorates increases the abundances of neutral molecules in the outer disk, highlighting the importance of taking into account the shape of the UV spectrum in protoplanetary disks. A recalculation of the photoreaction rates also affects the gas-phase chemistry due to the adjustment of the H/H2 and C+/C ratios. The disk ionization fraction is not significantly affected by the methods adopted to calculate the photochemistry and X-ray ionization. We determine that there is a probable "dead zone" where accretion is suppressed, present in a layer, Z/R lsim 0.1-0.2, in the disk midplane, within R ˜ 200 AU.

A search for line intensity enhancements in the far-UV spectra of active late-type stars arising from opacity

Context. Radiative transfer calculations have predicted intensity enhancements for optically thick emission lines, as opposed to the normal intensity reductions, for astrophysical plasmas under certain conditions. In particular, the results are predicted to be dependent both on the geometry of the emitting plasma and the orientation of the observer. Hence in principle the detection of intensity enhancement may provide a way of determining the geometry of an unresolved astronomical source.
Aims. To investigate such enhancements we have analysed a sample of active late-type stars observed in the far ultraviolet spectral region.
Methods. Emission lines of O vi in the FUSE satellite spectra of ϵ Eri, II Peg and Prox Cen were searched for intensity enhancements due to opacity.
Results. We have found strong evidence for line intensity enhancements due to opacity during active or flare-like activity for all three stars. The O vi 1032/1038 line intensity ratios, predicted to have a value of 2.0 in the optically thin case, are found to be up to ~30% larger during several orbital phases.
Conclusions. Our measurements, combined with radiative transfer models, allow us to constrain both the geometry of the O vi emitting regions in our stellar sources and the orientation of the observer. A spherical emitting plasma can be ruled out, as this would lead to no intensity enhancement. In addition, the theory tells us that the line-of-sight to the plasma must be close to perpendicular to its surface, as observations at small angles to the surface lead to either no intensity enhancement or the usual line intensity decrease over the optically thin value. For the future, we outline a laboratory experiment, that could be undertaken with current facilities, which would provide an unequivocal test of predictions of line intensity enhancement due to opacity, in particular the dependence on plasma geometry.

Double-detonation sub-Chandrasekhar supernovae: Synthetic observables for minimum helium shell mass models

In the double-detonation scenario for Type Ia supernovae, it is suggested that a detonation initiates in a shell of helium-rich material accreted from a companion star by a sub-Chandrasekhar-mass white dwarf. This shell detonation drives a shock front into the carbon-oxygen white dwarf that triggers a secondary detonation in the core. The core detonation results in a complete disruption of the white dwarf. Earlier studies concluded that this scenario has difficulties in accounting for the observed properties of Type Ia supernovae since the explosion ejecta are surrounded by the products of explosive helium burning in the shell. Recently, however, it was proposed that detonations might be possible for much less massive helium shells than previously assumed (Bildsten et al.). Moreover, it was shown that even detonations of these minimum helium shell masses robustly trigger detonations of the carbon-oxygen core (Fink et al.). Therefore, it is possible that the impact of the helium layer on observables is less than previously thought. Here, we present time-dependent multi-wavelength radiative transfer calculations for models with minimum helium shell mass and derive synthetic observables for both the optical and ? -ray spectral regions. These differ strongly from those found in earlier simulations of sub-Chandrasekhar-mass explosions in which more massive helium shells were considered. Our models predict light curves that cover both the range of brightnesses and the rise and decline times of observed Type Ia supernovae. However, their colors and spectra do not match the observations. In particular, their B - V colors are generally too red. We show that this discrepancy is mainly due to the composition of the burning products of the helium shell of the Fink et al. models which contain significant amounts of titanium and chromium. Using a toy model, we also show that the burning products of the helium shell depend crucially on its initial composition. This leads us to conclude that good agreement between sub-Chandrasekhar-mass explosions and observed Type Ia supernovae may still be feasible but further study of the shell properties is required.