Synthetic X-ray spectra for simulations of the dynamics of an accretion flow irradiated by a quasar

Ultraviolet and X-ray observations show evidence of outflowing gas around many active galactic nuclei. It has been proposed that some of these outflows are driven off gas infalling towards the central supermassive black hole. We perform radiative transfer calculations to compute the gas ionization state and the emergent X-ray spectra for both two- and three-dimensional (3D) hydrodynamical simulations of this outflow-from-inflow scenario. By comparison with observations, our results can be used to test the theoretical models and guide future numerical simulations. We predict both absorption and emission features, most of which are formed in a polar funnel of relatively dense (10 -10 g cm ) outflowing gas. This outflow causes strong absorption for observer orientation angles of ?35°. Particularly in 3D, the strength of this absorption varies significantly for different lines of sight owing to the fragmentary structure of the gas flow. Although infalling material occupies a large fraction of the simulation volume, we do not find that it imprints strong absorption features in the X-ray spectra since the ionization state is predicted to be very high. Thus, an absence of observed inflow absorption features does not exclude the models. The main spectroscopic consequence of the infalling gas is a Compton-scattered continuum component that partially re-fills the absorption features caused by the outflowing polar funnel. Fluorescence and scattering in the outflow are predicted to give rise to several emission features including a multicomponent Fe Ka emission complex for all observer orientations. For the hydrodynamical simulations considered, we predict both ionization states and column densities for the outflowing gas that are too high to be quantitatively consistent with well-observed X-ray absorption systems. Nevertheless, our results are qualitatively encouraging and further exploration of the model parameter space is warranted. Higher resolution hydrodynamic simulations are needed to determine whether the outflows fragment on scales unresolved in our current study, which may yield the denser lower ionization material that could reconcile the models and the observations. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 RAS.

Detonations in sub-chandrasekhar-mass C+O white dwarfs

Explosions of sub-Chandrasekhar-mass white dwarfs (WDs) are one alternative to the standard Chandrasekhar-mass model of Type Ia supernovae (SNe Ia). They are interesting since binary systems with sub-Chandrasekhar-mass primary WDs should be common and this scenario would suggest a simple physical parameter which determines the explosion brightness, namely the mass of the exploding WD. Here we perform one-dimensional hydrodynamical simulations, associated post-processing nucleosynthesis, and multi-wavelength radiation transport calculations for pure detonations of carbon-oxygen WDs. The light curves and spectra we obtain from these simulations are in good agreement with observed properties of SNe Ia. In particular, for WD masses from 0.97 to 1.15 Msun we obtain 56Ni masses between 0.3 and 0.8 Msun, sufficient to capture almost the complete range of SN Ia brightnesses. Our optical light curve rise times, peak colors, and decline timescales display trends which are generally consistent with observed characteristics although the range of B-band decline timescales displayed by our current set of models is somewhat too narrow. In agreement with observations, the maximum light spectra of the models show clear features associated with intermediate-mass elements and reproduce the sense of the observed correlation between explosion luminosity and the ratio of the Si II lines at ?6355 and ?5972. We therefore suggest that sub-Chandrasekhar-mass explosions are a viable model for SNe Ia for any binary evolution scenario leading to explosions in which the optical display is dominated by the material produced in a detonation of the primary WD. © 2010. The American Astronomical Society.

Three-dimensional delayed-detonation models with nucleosynthesis for type ia supernovae

We present results for a suite of 14 three-dimensional, high-resolution hydrodynamical simulations of delayed-detonation models of Type Ia supernova (SN Ia) explosions. This model suite comprises the first set of three-dimensional SN Ia simulations with detailed isotopic yield information. As such, it may serve as a data base for Chandrasekhar-mass delayed-detonation model nucleosynthetic yields and for deriving synthetic observables such as spectra and light curves. We employ aphysically motivated, stochastic model based on turbulent velocity fluctuations and fuel density to calculate in situ the deflagration-to-detonation transition probabilities. To obtain different strengths of the deflagration phase and thereby different degrees of pre-expansion, we have chosen a sequence of initial models with 1, 3, 5, 10, 20, 40, 100, 150, 200, 300 and 1600 (two different realizations) ignition kernels in a hydrostatic white dwarf with a central density of 2.9 × 10 g cm, as well as one high central density (5.5 × 10 g cm) and one low central density (1.0 × 10 g cm) rendition of the 100 ignition kernel configuration. For each simulation, we determined detailed nucleosynthetic yields by postprocessing10 tracer particles with a 384 nuclide reaction network. All delayed-detonation models result in explosions unbinding thewhite dwarf, producing a range of 56Ni masses from 0.32 to 1.11M. As a general trend, the models predict that the stableneutron-rich iron-group isotopes are not found at the lowest velocities, but rather at intermediate velocities (~3000×10 000 km s) in a shell surrounding a Ni-rich core. The models further predict relatively low-velocity oxygen and carbon, with typical minimum velocities around 4000 and 10 000 km s, respectively. © 2012 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society.

Plasma properties and Stokes profiles during the lifetime of a photospheric magnetic bright point

Aims: In this paper we aim to investigate the evolution of plasmaproperties and Stokes parameters in photospheric magnetic bright pointsusing 3D magneto-hydrodynamical simulations and radiative diagnostics ofsolar granulation.

Methods: Simulated time-dependent radiationparameters and plasma properties were investigated throughout theevolution of a bright point. Synthetic Stokes profiles for the FeI630.25 nm line were calculated, which also allowed the evolution of theStokes-I line strength and Stokes-V area and amplitude asymmetries to beinvestigated.

Results: Our results are consistent withtheoretical predictions and published observations describing convectivecollapse, and confirm this as the bright point formation process.Through degradation of the simulated data to match the spatialresolution of SOT, we show that high spatial resolution is crucial forthe detection of changing spectro-polarimetric signatures throughout amagnetic bright point's lifetime. We also show that the signaturedownflow associated with the convective collapse process tends towardszero as the radiation intensity in the bright point peaks, because ofthe magnetic forces present restricting the flow of material in the fluxtube.

Multidimensional modelling of X-ray spectra for AGN accretion disc outflows - III. Application to a hydrodynamical simulation

We perform multidimensional radiative transfer simulations to compute spectra for a hydrodynamical simulation of a line-driven accretion disc wind from an active galactic nucleus. The synthetic spectra confirm expectations from parametrized models that a disc wind can imprint a wide variety of spectroscopic signatures including narrow absorption lines, broad emission lines and a Compton hump. The formation of these features is complex with contributions originating from many of the different structures present in the hydrodynamical simulation. In particular, spectral features are shaped both by gas in a successfully launched outflow and in complex flows where material is lifted out of the disc plane but ultimately falls back. We also confirm that the strong Fe Ka line can develop a weak, red-skewed line wing as a result of Compton scattering in the outflow. In addition, we demonstrate that X-ray radiation scattered and reprocessed in the flow has a pivotal part in both the spectrum formation and determining the ionization conditions in the wind. We find that scattered radiation is rather effective in ionizing gas which is shielded from direct irradiation from the central source. This effect likely makes the successful launching of a massive disc wind somewhat more challenging and should be considered in future wind simulations. © 2010 The Authors. Journal compilation © 2010 RAS.

Two-dimensional Monte Carlo simulations of H i line formation in massive young stellar object disc winds

Massive young stellar objects (YSOs) are powerful infrared Hi line emitters. It has been suggested that these lines form in an outflow from a disc surrounding the YSO. Here, new two-dimensional Monte Carlo radiative transfer calculations are described which test this hypothesis. Infrared spectra are synthesized for a YSO disc wind model based on earlier hydrodynamical calculations. The model spectra are in qualitative agreement with the observed spectra from massive YSOs, and therefore provide support for a disc wind explanation for the Hi lines. However, there are some significant differences: the models tend to overpredict the Bra/Br? ratio of equivalent widths and produce line profiles which are slightly too broad and, in contrast to typical observations, are double-peaked. The interpretation of these differences within the context of the disc wind picture and suggestions for their resolution via modifications to the assumed disc and outflow structure are discussed. © 2005 RAS.

Multi-Dimensional Simulations of the Expanding Supernova Remnant of SN 1987A

The expanding remnant from SN 1987A is an excellent laboratory for investigating the physics of supernovae explosions. There is still a large number of outstanding questions, such as the reason for the asymmetric radio morphology, the structure of the pre-supernova environment, and the efficiency of particle acceleration at the supernova shock. We explore these questions using three-dimensional simulations of the expanding remnant between days 820 and 10,000 after the supernova. We combine a hydrodynamical simulation with semi-analytic treatments of diffusive shock acceleration and magnetic field amplification to derive radio emission as part of an inverse problem. Simulations show that an asymmetric explosion, combined with magnetic field amplification at the expanding shock, is able to replicate the persistent one-sided radio morphology of the remnant. We use an asymmetric Truelove & McKee progenitor with an envelope mass of 10 M-circle dot and an energy of 1.5 x 10(44) J. A termination shock in the progenitor's stellar wind at a distance of 0 ''.43-0 ''.51 provides a good fit to the turn on of radio emission around day 1200. For the H II region, a minimum distance of 0 ''.63 +/- 0 ''.01 and maximum particle number density of (7.11 +/- 1.78) x 10(7) m(-3) produces a good fit to the evolving average radius and velocity of the expanding shocks from day 2000 to day 7000 after explosion. The model predicts a noticeable reduction, and possibly a temporary reversal, in the asymmetric radio morphology of the remnant after day 7000, when the forward shock left the eastern lobe of the equatorial ring.

An investigation into the use of atomic datasets in simulations of the Ni-like gadolinium x-ray laser

The effect of differing the datasets used in the modelling of the Ni-like Gd x-ray laser (XRL) is examined through the 1.50 hydro-atomic code, EHYBRID. Two atomic datasets, including energy levels and radiative and collisional excitation rates, are used as input data for the code. It is found that the behaviour of the XRL is somewhat different than might be expected from superficial examination of the atomic data. The similarities in the gain profiles at low densities are found to have encouraging implications. in our attempts to model XRLs.