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.

Multidimensional modelling of X-ray spectra for AGN accretion disc outflows

We use a multidimensional Monte Carlo code to compute X-ray spectra for a variety of active galactic nucleus (AGN) disc-wind outflow geometries. We focus on the formation of blueshifted absorption features in the Fe K band and show that line features similar to those which have been reported in observations are often produced for lines of sight through disc-wind geometries. We also discuss the formation of other spectral features in highly ionized outflows. In particular, we show that, for sufficiently high wind densities, moderately strong Fe K emission lines can form and that electron scattering in the flow may cause these lines to develop extended red wings. We illustrate the potential relevance of such models to the interpretation of real X-ray data by comparison with observations of a well-known AGN, Mrk 766. Journal compilation © 2008 RAS.

Detection of magnetic dipole lines of Fe XII in the ultraviolet spectrum of the dwarf star ε Eri

We report observations of the dwarf star e Eri (K2V) made with the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope. The high sensitivity of the STIS instrument has allowed us to detect the magnetic dipole transitions of Fe XII at 1242.00 and 1349 38 Å for the first time in a star other than the Sun. The width of the stronger line at 1242.00 Å has also been measured; such measurements are not possible for the permitted lines of Fe XII in the extreme-ultraviolet. To within the accuracy of the measurements the N v and the Fe XII lines occur at their rest wavelengths. Electron densities and linewidths have been measured from other transition region lines. Together, these can be used to investigate the non-thermal energy flux in the lower and upper transition regions, which is useful in constraining possible heating processes. The Fe XII lines are also present in archival STIS spectra of other G/K-type dwarfs.

Analytical approximations to numerical solutions of theoretical emission measure distributions

Emission line fluxes from cool stars are widely used to establish an apparent emission measure distribution, EmdApp(Te), between temperatures characteristic of the low transition region and the low corona. The true emission measure distribution, EmdTrue(Te), is determined by the energy balance and geometry adopted and, with a numerical model, can be used to predict EmdApp(Te), to guide further modelling. The scaling laws that exist between coronal parameters arise from the dimensions of the terms in the energy balance equation. Here, analytical approximations to numerical solutions for EmdTrue(Te) are presented, which show how the constants in the coronal scaling laws are determined. The apparent emission measure distributions show a minimum value at some T0 and a maximum at the mean coronal temperature Tc (although in some stars, emission from active regions can contribute). It is shown that, for the energy balance and geometry adopted, the analytical values of the emission measure and electron pressure at T0 and Tc depend on only three parameters: the stellar surface gravity and the values of T0 and Tc. The results are tested against full numerical solutions for e Eri (K2 V) and are applied to Procyon (a CMi, F5 IV/V). The analytical approximations can be used to restrict the required range of full numerical solutions, to check the assumed geometry and to show where the adopted energy balance may not be appropriate. © 2011 The Authors Monthly Notices of the Royal Astronomical Society © 2011 RAS.

The electron pressure in the outer atmosphere of ∈ Eri (K2 V)

Observations of ? Eri (K2 V) have been made with the Space Telescope Imaging Spectrograph on the Hubble Space Telescope. The spectra obtained show a number of emission lines which can be used to determine, or place limits on, the electron density and pressure. Values of the electron pressure are required in order to make quantitative models of the transition region and inner corona from absolute line fluxes, and to constrain semi-empirical models of the chromosphere. Using line flux ratios in Si II and O IV a mean electron pressure of P = NT = 4.8 × 10 cm K is derived. This value is compatible with the lower and upper limits to P found from flux ratios in C III, O V and Fe XII. Some inconsistencies which may be because of small uncertainties in the atomic data used are discussed.

Modelling the H Lyman lines in evolved late-type stars

The importance of partial redistribution (PRD) in the modelling of the Lyman a and Lyman ß emission lines of hydrogen in stellar atmospheres is examined using simple atmospheric models of a range of late-type stars. These models represent the subgiant Procyon (F5 IV-V), and the two giants ß Gem (K0 III) and a Tau (K5 III). These stars are selected to span a wide range of surface gravities: 1.25 <log g <4.00. The calculations are performed using the computer code MULTI with the modifications made by Hubeny & Lites. It is found that PRD effects are highly significant, both in the direct prediction of the Lyman line profiles and in the application of hydrostatic equilibrium to calculate the atmospheric electron density in static atmospheric models.

Modelling the chromosphere and transition region of ε Eri (K2 V)

Measurements of ultraviolet line fluxes from Space Telescope Imaging Spectrograph and Far-Ultraviolet Spectroscopic Explorer spectra of the K2-dwarf e Eri are reported. These are used to develop new emission measure distributions and semi-empirical atmospheric models for the chromosphere and lower transition region of the star. These models are the most detailed constructed to date for a main-sequence star other than the Sun. New ionization balance calculations, which account for the effect of finite density on dielectronic recombination rates, are presented for carbon, nitrogen, oxygen and silicon. The results of these calculations are significantly different from the standard Arnaud & Rothenflug ion balance, particularly for alkali-like ions. The new atmospheric models are used to place constraints on possible first ionization potential (FIP)-related abundance variations in the lower atmosphere and to discuss limitations of single-component models for the interpretation of certain optically thick line fluxes. © 2005 RAS.

On the filling factor of emitting material in the upper atmosphere of ε Eri (K2 V)

The emission measure distribution in the upper transition region and corona of e Eri is derived from observed emission-line fluxes. Theoretical emission measure distributions are calculated assuming that the radiation losses are balanced by the net conductive flux. We discuss how the area factor of the emitting regions as a function of temperature can be derived from a comparison between these emission measure distributions. It is found that the filling factor varies from ~0.2 in the mid-transition region to ~1.0 in the inner corona. The sensitivity of these results to the adopted ion fractions, the iron abundance and other parameters is discussed. The area factors found are qualitatively similar to the observed structure of the solar atmosphere, and can be used to constrain two-component models of the chromosphere. Given further observations, the method could be applied to investigate the trends in filling factors with indicators of stellar activity.

A study of velocity fields in the transition region of ε Eri (K2 V)

Analyses of the widths and shifts of optically thin emission lines in the ultraviolet spectrum of the active dwarf e Eri (K2 V) are presented. The spectra were obtained using the Space Telescope Imaging Spectrograph on the Hubble Space Telescope and the Far Ultraviolet Spectroscopic Explorer. The linewidths are used to find the non-thermal energy density and its variation with temperature from the chromosphere to the upper transition region. The energy fluxes that could be carried by Alfvén and acoustic waves are investigated, to test their possible roles in coronal heating. Acoustic waves do not appear to be a viable means of coronal heating. There is, in principle, ample flux in Alfvén waves, but detailed calculations of wave propagation are required before definite conclusions can be drawn concerning their viability. The high sensitivity and spectral resolution of the above instruments have allowed two-component Gaussian fits to be made to the profiles of the stronger transition region lines. The broad and narrow components that result share some similarities with those observed in the Sun, but in e Eri the broad component is redshifted relative to the narrow component and contributes more to the total line flux. The possible origins of the two components and the energy fluxes implied are discussed. On balance our results support the conclusion of Wood, Linsky & Ayres, that the narrow component is related to Alfvén waves reaching to the corona, but the origin of the broad component is not clear.

Mass-loss rates for hot luminous stars:The influence of line branching

The effect of photon frequency redistribution by line branching on mass-loss rates for hot luminous stars is investigated. Monte Carlo simulations are carried out for a range of OB star models which show that previous mass-loss calculations which neglect non-resonance line scattering overestimate mass-loss rates for luminous O stars by ~20 per cent. For luminous B stars the effect is somewhat larger, typically ~50 per cent. A Wolf-Rayet star model is used to investigate line branching in the strong wind limit. In this case the effect of line branching is much greater, giving mass-loss rates that are smaller by a factor ~3 from computations which neglect branching.

Modelling the X-ray spectra of high-velocity outflows from quasars

High-velocity outflows from supermassive black holes have been invoked to explain the recent identification of strong absorption features in the hard X-ray spectra of several quasars. Here, Monte Carlo radiative transfer calculations are performed to synthesize X-ray spectra from models of such flows. It is found that simple, parametric biconical outflow models with plausible choices for the wind parameters predict spectra that are in good qualitative agreement with observations in the 2-10 keV band. The influence on the spectrum of both the mass-loss rate and opening angle of the flow are considered: the latter is important since photon leakage plays a significant role in establishing an ionization gradient within the flow, a useful discriminant between spherical and conical outflow for this and other applications. Particular attention is given to the bright quasar PG 1211+143 for which constraints on the outflow geometry and mass-loss rate are discussed subject to the limitations of the currently available observational data.

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.

Evidence for relativistic features in the X-ray spectrum of Mrk 335

We present an analysis of hard X-ray features in the spectrum of the bright Sy 1 galaxy Mrk 335 observed by the XMM-Newton satellite. Our analysis confirms the presence of a broad, ionized Fe Ka emission line in the spectrum, first found by Gondoin et al. The broad line can be modelled successfully by relativistic accretion disc reflection models. This interpretation is unusually robust in the case of Mrk 335 because of the lack of any ionized ('warm') absorber and the absence a clear narrow core to the line. Partial covering by neutral gas cannot, however, be ruled out statistically as the origin of the broad residuals. Regardless of the underlying continuum we report, for the first time in this source, the detection of a narrow absorption feature at the rest frame energy of ~5.9 keV. If the feature is identified with a resonance absorption line of iron in a highly ionized medium, the redshift of the line corresponds to an inflow velocity of ~0.11-0.15c. We present a simple model for the inflow, accounting approximately for relativistic and radiation pressure effects, and use Monte Carlo methods to compute synthetic spectra for qualitative comparison with the data. This modelling shows that the absorption feature can plausibly be reproduced by infalling gas providing that the feature is identified with Fe xxvi. We require the inflowing gas to extend over a limited range of radii at a few tens of r to match the observed feature. The mass accretion rate in the flow corresponds to 60 per cent of the Eddington limit, in remarkable agreement with the observed rate. The narrowness of the absorption line tends to argue against a purely gravitational origin for the redshift of the line, but given the current data quality we stress that such an interpretation cannot be ruled out. © 2006 The Authors.

Light curves for off-centre ignition models of Type Ia supernovae

Motivated by recent models involving off-centre ignition of Type Ia supernova explosions, we undertake three-dimensional time-dependent radiation transport simulations to investigate the range of bolometric light-curve properties that could be observed from supernovae in which there is a lop-sided distribution of the products from nuclear burning. We consider both a grid of artificial toy models which illustrate the conceivable range of effects and a recent three-dimensional hydrodynamical explosion model. We find that observationally significant viewing angle effects are likely to arise in such supernovae and that these may have important ramifications for the interpretation of the observed diversity of Type Ia supernova and the systematic uncertainties which relate to their use as standard candles in contemporary cosmology. © 2007 RAS.

2D simulations of the double-detonation model for thermonuclear transients from low-mass carbon-oxygen white dwarfs

Thermonuclear explosions may arise in binary star systems in which a carbon-oxygen (CO) white dwarf (WD) accretes helium-rich material from a companion star. If the accretion rate allows a sufficiently large mass of helium to accumulate prior to ignition of nuclear burning, the helium surface layer may detonate, giving rise to an astrophysical transient. Detonation of the accreted helium layer generates shock waves that propagate into the underlying CO WD. This might directly ignite a detonation of the CO WD at its surface (an edge-lit secondary detonation) or compress the core of the WD sufficiently to trigger a CO detonation near the centre. If either of these ignition mechanisms works, the two detonations (helium and CO) can then release sufficient energy to completely unbind the WD. These 'double-detonation' scenarios for thermonuclear explosion of WDs have previously been investigated as a potential channel for the production of Type Ia supernovae from WDs of ~ 1 M . Here we extend our 2D studies of the double-detonation model to significantly less massive CO WDs, the explosion of which could produce fainter, more rapidly evolving transients. We investigate the feasibility of triggering a secondary core detonation by shock convergence in low-mass CO WDs and the observable consequences of such a detonation. Our results suggest that core detonation is probable, even for the lowest CO core masses that are likely to be realized in nature. To quantify the observable signatures of core detonation, we compute spectra and light curves for models in which either an edge-lit or compression-triggered CO detonation is assumed to occur. We compare these to synthetic observables for models in which no CO detonation was allowed to occur. If significant shock compression of the CO WD occurs prior to detonation, explosion of the CO WD can produce a sufficiently large mass of radioactive iron-group nuclei to significantly affect the light curves. In particular, this can lead to relatively slow post-maximum decline. If the secondary detonation is edge-lit, however, the CO WD explosion primarily yields intermediate-mass elements that affect the observables more subtly. In this case, near-infrared observations and detailed spectroscopic analysis would be needed to determine whether a core detonation occurred. We comment on the implications of our results for understanding peculiar astrophysical transients including SN 2002bj, SN 2010X and SN 2005E. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 RAS.