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.

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.

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.

High-contrast imaging search for planets and brown dwarfs around the most massive stars in the solar neighborhood

There has been a long-standing discussion in the literature as to whether core accretion or disk instability is the dominant mode of planet formation. Over the last decade, several lines of evidence have been presented showing that core accretion is most likely the dominant mechanism for the close-in population of planets probed by radial velocity and transits. However, this does not by itself prove that core accretion is the dominant mode for the total planet population, since disk instability might conceivably produce and retain large numbers of planets in the far-out regions of the disk. If this is a relevant scenario, then the outer massive disks of B-stars should be among the best places for massive planets and brown dwarfs to form and reside. In this study, we present high-contrast imaging of 18 nearby massive stars of which 15 are in the B2-A0 spectral-type range and provide excellent sensitivity to wide companions. By comparing our sensitivities to model predictions of disk instability based on physical criteria for fragmentation and cooling, and using Monte Carlo simulations for orbital distributions, we find that ~85% of such companions should have been detected in our images on average. Given this high degree of completeness, stringent statistical limits can be set from the null-detection result, even with the limited sample size. We find that

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.

Magma Evolution in the Primitive, Intra-oceanic Tonga Arc: Rapid Petrogenesis of Dacites at Fonualei Volcano

Fonualei is unusual amongst subaerial volcanoes in the Tonga arc because it has erupted dacitic vesicular lavas, tuffs and phreomagmatic deposits for the last 165 years. The total volume of dacite may approach 5 km(3) and overlies basal basaltic andesite and andesite lavas that are constrained to be less than a few millennia in age. All of the products are crystal-poor and formed from relatively low-viscosity magmas inferred to have had temperatures of 1100-1000 degrees C, 2-4 wt % H2O and oxygen fugacities 1-2 log units above the quartz-fayalite-magnetite buffer. Major and trace element data, along with Sr-Nd-Pb and U-Th-Ra isotope data, are used to assess competing models for the origin of the dacites. Positive correlations between Sc and Zr and Sr rule out evolution of the within-dacite compositional array by closed-system crystal fractionation of a single magma batch. An origin by partial melting of lower crustal amphibolites cannot reproduce these data trends or, arguably, any of the dacites either. Instead, we develop a model in which the dacites reflect mixing between two dacitic magmas, each the product of fractional crystallization of basaltic andesite magmas formed by different degrees of partial melting. Mixing was efficient because the two magmas had similar temperatures and viscosities. This is inferred to have occurred at shallow (2-6 km) depths beneath the volcano. U-Th-Ra disequilibria in the basaltic andesite and andesite indicate that the parental magmas had fluids added to their mantle source regions less than 8 kyr ago and that fractionation to the dacitic compositions took less than a few millennia. The 165 year eruption period for the dacites implies that mixing occurred on a similar timescale, possibly during ascent in conduits. The composition of the dacites renders them unsuitable candidates as contributors to average continental crust.

Rapid magmatic processes accompany arc-continent collision: the Western Bismarck arc, Papua New Guinea

New U-Th-Ra, major and trace element, and Sr-Nd-Pb isotope data are presented for young lavas from the New Britain and Western Bismarck arcs in Papua New Guinea. New Britain is an oceanic arc, whereas the latter is the site of an arc-continent collision. Building on a recent study of the Manus Basin, contrasts between the two arcs are used to evaluate the processes and timescales of magma generation accompanying arc-continent collision and possible slab detachment. All three suites share many attributes characteristic of arc lavas that can be ascribed to the addition of a regionally uniform subduction component derived from the subducting altered oceanic crust and sediment followed by dynamic melting of the modified mantle. However, the Western Bismarck arc lavas diverge from the Pb isotope mixing array formed by the New Britain and the Manus Basin lavas toward elevated Pb-208/Pb-204. We interpret this to reflect a second and subsequent addition of sediment melt at crustal depth during collision. U-238 and Ra-226 excesses are preserved in all of the lavas and are greatest in the Western Bismarck arc. High-Mg andesites with high Sr/Y ratios in the westernmost arc are attributed to recent shallow mantle flux melting at the slab edge. Data for two historical rhyolites are also presented. Although these rhyolites formed in quite different tectonic settings and display different geochemical and isotopic compositions, both formed from mafic parents within millennia.

Complex organic molecules in protoplanetary disks

Context. Protoplanetary disks are vital objects in star and planet formation, possessing all the material, gas and dust, which may form a planetary system orbiting the new star. Small, simple molecules have traditionally been detected in protoplanetary disks; however, in the ALMA era, we expect the molecular inventory of protoplanetary disks to significantly increase.

Aims. We investigate the synthesis of complex organic molecules (COMs) in protoplanetary disks to put constraints on the achievable chemical complexity and to predict species and transitions which may be observable with ALMA.

Methods. We have coupled a 2D steady-state physical model of a protoplanetary disk around a typical T Tauri star with a large gas-grain chemical network including COMs. We compare the resulting column densities with those derived from observations and perform ray-tracing calculations to predict line spectra. We compare the synthesised line intensities with current observations and determine those COMs which may be observable in nearby objects. We also compare the predicted grain-surface abundances with those derived from cometary comae observations.

Results. We find COMs are efficiently formed in the disk midplane via grain-surface chemical reactions, reaching peak grain-surface fractional abundances similar to 10(-6)-10(-4) that of the H nuclei number density. COMs formed on grain surfaces are returned to the gas phase via non-thermal desorption; however, gas-phase species reach lower fractional abundances than their grain-surface equivalents, similar to 10(-12)-10(-7). Including the irradiation of grain mantle material helps build further complexity in the ice through the replenishment of grain-surface radicals which take part in further grain-surface reactions. There is reasonable agreement with several line transitions of H2CO observed towards T Tauri star-disk systems. There is poor agreement with HC3(N) lines observed towards LkCa 15 and GO Tau and we discuss possible explanations for these discrepancies. The synthesised line intensities for CH3OH are consistent with upper limits determined towards all sources. Our models suggest CH3OH should be readily observable in nearby protoplanetary disks with ALMA; however, detection of more complex species may prove challenging, even with ALMA "Full Science" capabilities. Our grain-surface abundances are consistent with those derived from cometary comae observations providing additional evidence for the hypothesis that comets (and other planetesimals) formed via the coagulation of icy grains in the Sun's natal disk.