Gas and grain chemistry in a protoplanetary disk

The chemistry in a protoplanetary accretion disk is modelled between a radius of 100 and 0.1 AU of the central object. We find that interaction of the gas with the dust grains is very important, both by removing a large fraction of the material from the gas in the outer regions and through the chemical reactions which can occur on the dust grain surfaces. In addition, collision with grains neutralises gaseous ions effectively and keeps the ionization fraction low. This results in a chemistry which is dominated by neutral-neutral reactions, even if ionization is provided by cosmic rays or by the decay of radioactive isotopes. We model the effects of two desorption processes with very different efficiencies and find that while these produce similar results over much of the disk for many species, some molecules are extremely sensitive to the nature of the desorption and may one day be used as an observational test for the desorption process.

Roche tomography of cataclysmic variables - V. A high-latitude star-spot on RU Pegasi

We present Roche tomograms of the secondary star in the dwarf nova system RU Pegasi derived from blue and red arm ISIS data taken on the 4.2-m William Herschel Telescope. We have applied the entropy landscape technique to determine the system parameters and obtained component masses of M1 = 1.06 Msun, M2 = 0.96 Msun, an orbital inclination angle of i = 43 degrees, and an optimal systemic velocity of gamma = 7 km/s. These are in good agreement with previously published values. Our Roche tomograms of the secondary star show prominent irradiation of the inner Lagrangian point due to illumination by the disc and/or bright spot, which may have been enhanced as RU Peg was in outburst at the time of our observations.We find that this irradiation pattern is axi-symmetric and confined to regions of the star which have a direct view of the accretion regions. This is in contrast to previous attempts to map RU Peg which suggested that the irradiation pattern was non-symmetric and extended beyond the terminator. We also detect additional inhomogeneities in the surface distribution of stellar atomic absorption that we ascribe to the presence of a large star-spot. This spot is centred at a latitude of about 82 degrees and covers approximately 4 per cent of the total surface area of the secondary. In keeping with the high latitude spots mapped on the cataclysmic variables AE Aqr and BV Cen, the spot on RU Peg also appears slightly shifted towards the trailing hemisphere of the star. Finally, we speculate that early mapping attempts which indicated non-symmetric irradiation patterns which extended beyond the terminator of CV donors could possibly be explained by a superposition of symmetric heating and a large spot.

The geometry and ionization structure of the wind in the eclipsing nova-like variables RW Tri and UX UMa

The UV spectra of nova-like variables are dominated by emission from the accretion disk, modified by scattering in a wind emanating from the disk. Here, we model the spectra of RW Tri and UX UMa, the only two eclipsing nova-like variables which have been observed with the Hubble Space Telescope in the far-ultraviolet, in an attempt to constrain the geometry and the ionization structure of their winds. Using our Monte Carlo radiative transfer code, we computed spectra for simply parameterized axisymmetric biconical outflow models and were able to find plausible models for both systems. These reproduce the primary UV resonance lines-N v, Si iv, and C iv-in the observed spectra in and out of eclipse. The distribution of these ions in the wind models is similar in both cases as is the extent of the primary scattering regions in which these lines are formed. The inferred mass-loss rates are 6%-8% of the mass accretion rates for the systems. We discuss the implication of our point models for our understanding of accretion disk winds in cataclysmic variables. © 2010. The American Astronomical Society. All rights reserved.

Modeling the FeK line profiles in type i active galactic nuclei with a compton-thick disk wind

We have modeled a small sample of Seyfert galaxies that were previously identified as having simple X-ray spectra with little intrinsic absorption. The sources in this sample all contain moderately broad components of FeK-shell emission and are ideal candidates for testing the applicability of a Compton-thick accretion disk wind model to active galactic nucleus (AGN) emission components. Viewing angles through the wind allow the observer to see the absorption signature of the gas, whereas face-on viewing angles allow the observer to see the scattered light from the wind. We find that the FeK emission line profiles are well described with a model of a Compton-thick accretion disk wind of solar abundances, arising tens to hundreds of gravitational radii from the central black hole. Further, the fits require a neutral component of FeKa emission that is too narrow to arise from the inner part of the wind, and likely comes from a more distant reprocessing region. Our study demonstrates that a Compton-thick wind can have a profound effect on the observed X-ray spectrum of an AGN, even when the system is not viewed through the flow. © 2012. The American Astronomical Society. All rights reserved..

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.

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.

On the brightness distribution of type ia supernovae from violent white dwarf mergers

We investigate the brightness distribution expected for thermonuclear explosions that might result from the ignition of a detonation during the violent merger of white dwarf (WD) binaries. Violent WD mergers are a subclass of the canonical double degenerate scenario where two carbon-oxygen (CO) WDs merge when the larger WD fills its Roche lobe. Determining their brightness distribution is critical for evaluating whether such an explosion model could be responsible for a significant fraction of the observed population of Type Ia supernovae (SNe Ia). We argue that the brightness of an explosion realized via the violent merger model is mainly determined by the mass of Ni produced in the detonation of the primary COWD. To quantify this link, we use a set of sub-Chandrasekhar mass WD detonation models to derive a relationship between primary WD mass (m) and expected peak bolometric brightness (M). We use this m-M relationship to convert the masses of merging primary WDs from binary population models to a predicted distribution of explosion brightness. We also investigate the sensitivity of our results to assumptions about the conditions required to realize a detonation during violent mergers ofWDs. We find a striking similarity between the shape of our theoretical peak-magnitude distribution and that observed for SNe Ia: our model produces a M distribution that roughly covers the range and matches the shape of the one observed for SNe Ia. However, this agreement hinges on a particular phase of mass accretion during binary evolution: the primary WD gains ~0.15-0.35M? from a slightly evolved helium star companion. In our standard binary evolution model, such an accretion phase is predicted to occur for about 43 per cent of all binary systems that ultimately give rise to binary CO WD mergers. We also find that with high probability, violent WD mergers involving the most massive primaries (?1.3M?, which should produce bright SNe) have delay times ?500 Myr. © 2012 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society.

On the relativistic iron line and soft excess in the Seyfert 1 galaxy Markarian 335

We report a 133-ks XMM-Newton observation of the Seyfert 1 galaxy Markarian 335. The 0.4-12 keV spectrum contains an underlying power-law continuum, a soft excess below 2 keV, and a double-peaked iron emission feature in the 6-7 keV range. We investigate the possibility that the double-peaked emission might represent the characteristic signature of the accretion disc. Detailed investigations show that a moderately broad accretion disc line is most likely present, but that the peaks may be due to narrower components from more distant material. The peaks at 6.4 and 7 keV can be identified, respectively, with the molecular torus in active galactic nucleus unification schemes, and very highly ionized, optically thin gas filling the torus. The X-ray variability spectra on both long (~100 ks) and short (~1 ks) time-scales do not support the recent suggestion that the soft excess is an artefact of variable, moderately ionized absorption. © 2007 The Authors. Journal compilation © 2007 RAS.

A simple disc wind model for broad absorption line quasars

Approximately 20 per cent of quasi-stellar objects (QSOs) exhibit broad, blue-shifted absorption lines in their ultraviolet spectra. Such features provide clear evidence for significant outflows from these systems, most likely in the form of accretion disc winds. These winds may represent the ‘quasar’ mode of feedback that is often invoked in galaxy formation/evolution models, and they are also key to unification scenarios for active galactic nuclei (AGN) and QSOs. To test these ideas, we construct a simple benchmark model of an equatorial, biconical accretion disc wind in a QSO and use a Monte Carlo ionization/radiative transfer code to calculate the ultraviolet spectra as a function of viewing angle. We find that for plausible outflow parameters, sightlines looking directly into the wind cone do produce broad, blue-shifted absorption features in the transitions typically seen in broad absorption line (BAL) QSOs. However, our benchmark model is intrinsically X-ray weak in order to prevent overionization of the outflow, and the wind does not yet produce collisionally excited line emission at the level observed in non-BAL QSOs. As a first step towards addressing these shortcomings, we discuss the sensitivity of our results to changes in the assumed X-ray luminosity and mass-loss rate, Ṁwind. In the context of our adopted geometry, Ṁwind ∼ Ṁacc is required in order to produce significant BAL features. The kinetic luminosity and momentum carried by such outflows would be sufficient to provide significant feedback.

Variability of the high-velocity outflow in the quasar PDS 456

We present a comparison of two Suzaku X-ray observations of the nearby (z = 0.184), luminous (L ∼ 10 erg s) type I quasar, PDS 456. A new 125 ks Suzaku observation in 2011 caught the quasar during a period of low X-ray flux and with a hard X-ray spectrum, in contrast with a previous 190 ks Suzaku observation in 2007 when the quasar appeared brighter and had a steep (Γ > 2) X-ray spectrum. The 2011 X-ray spectrum contains a pronounced trough near 9 keV in the quasar rest frame, which can be modeled with blueshifted iron K-shell absorption, most likely from the He- and H-like transitions of iron. The absorption trough is observed at a similar rest-frame energy as in the earlier 2007 observation, which appears to confirm the existence of a persistent high-velocity wind in PDS 456, at an outflow velocity of 0.25-0.30c. The spectral variability between 2007 and 2011 can be accounted for by variations in a partial covering absorber, increasing in covering fraction from the brighter 2007 observation to the hard and faint 2011 observation. Overall, the low-flux 2011 observation can be explained if PDS 456 is observed at relatively low inclination angles through a Compton-thick wind, originating from the accretion disk, which significantly attenuates the X-ray flux from the quasar. © 2014. The American Astronomical Society. All rights reserved.

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.

The spotty donor star in the X-ray transient Cen X-4

We accurately determine the fundamental system parameters of the neutron star X-ray transient Cen X-4 solely using phase-resolved high-resolution UV-Visual Echelle Spectrograph spectroscopy. We first determine the radial-velocity curve of the secondary star and then model the shape of the phase-resolved absorption line profiles using an X-ray binary model. The model computes the exact rotationally broadened, phase-resolved spectrum and does not depend on assumptions about the rotation profile, limb-darkening coefficients and the effects of contamination from an accretion disc. We determine the secondary star-to-neutron star binary mass ratio to be 0.1755 ± 0.0025, which is an order of magnitude more accurate than previous estimates. We also constrain the inclination angle to be 32^{+8}_{-2} degrees. Combining these values with the results of the radial-velocity study gives a neutron star mass of 1.94^{+0.37}_{-0.85}M⊙ consistent with previous estimates. Finally, we perform the first Roche tomography reconstruction of the secondary star in an X-ray binary. The tomogram reveals surface inhomogeneities that are due to the presence of cool starspots. A large cool polar spot, similar to that seen in Doppler images of rapidly rotating isolated stars, is present on the Northern hemisphere of the K7 secondary star and we estimate that ~4 percent of the total surface area of the donor star is covered with spots.This evidence for starspots supports the idea that magnetic braking plays an important role in the evolution of low-mass X-ray binaries.

Water in protoplanetary disks

Infrared water line emission from protoplanetary disks, recently observed by the Spitzer and Herschel space telescopes, is thought to trace the surface layer of the inner to outer regions of the disks. We have modelled the water abundance profile and line emission, especially focusing on the effects of dust size growth and turbulent mixing. Comparison between model calculations and observations suggests a small grain model with turbulent mixing is preferred. Copyright © International Astronomical Union 2014.

Roche tomography of cataclysmic variables - VI. Differential rotation of AE Aqr - not tidally locked!

We present Roche tomograms of the K4V secondary star in the cataclysmic variable AE Aqr, reconstructed from two data sets taken 9 d apart, and measure the differential rotation of the stellar surface. The tomograms show many large, cool starspots, including a large high-latitude spot and a prominent appendage down the trailing hemisphere. We find two distinct bands of spots around 22° and 43° latitude, and estimate a spot coverage of 15.4-17 per cent on the Northern hemisphere. Assuming a solar-like differential rotation law, the differential rotation of AE Aqr was measured using two different techniques. The first method yields an equator-pole lap time of 269 d and the second yields a lap time of 262 d. This shows that the star is not fully tidally locked, as was previously assumed for CVs, but has a co-rotation latitude of ˜40°. We discuss the implications that these observations have on stellar dynamo theory, as well as the impact that spot traversal across the L1 point may have on accretion rates in CVs as well as some of their other observed properties. The entropy landscape technique was applied to determine the system parameters of AE Aqr. For the two independent data sets, we find M1 = 1.20 and 1.17 M⊙, M2 = 0.81 and 0.78 M⊙, and orbital inclinations of 50° to 51° at optimal systemic velocities of γ = -64.7 and -62.9 km s-1.