3D deflagration simulations leaving bound remnants:A model for 2002cx-like Type Ia supernovae

2002cx-like supernovae are a sub-class of sub-luminous Type Ia supernovae (SNe). Their light curves and spectra are characterized by distinct features that indicate strong mixing of the explosion ejecta. Pure turbulent deflagrations have been shown to produce such mixed ejecta. Here, we present hydrodynamics, nucleosynthesis and radiative-transfer calculations for a 3D full-star deflagration of a Chandrasekhar-mass white dwarf. Our model is able to reproduce the characteristic observational features of SN 2005hk (a prototypical 2002cx-like supernova), not only in the optical, but also in the near-infrared. For that purpose we present, for the first time, five near-infrared spectra of SN 2005hk from -0.2 to 26.6 d with respect to B-band maximum. Since our model burns only small parts of the initial white dwarf, it fails to completely unbind the white dwarf and leaves behind a bound remnant of ~1.03Mconsisting mainly of unburned carbon and oxygen, but also enriched by some amount of intermediate-mass and iron-group elements from the explosion products that fall back on the remnant.We discuss possibilities for detecting this bound remnant and how it might influence the late-time observables of 2002cx-like SNe. © 2013 The Authors Published by Oxford University Press on behalf of the Royal 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.

Constraining type Ia supernova models:SN 2011fe as A test case

The nearby supernova SN 2011fe can be observed in unprecedented detail. Therefore, it is an important test case for Type Ia supernova (SN Ia) models, which may bring us closer to understanding the physical nature of these objects. Here, we explore how available and expected future observations of SN 2011fe can be used to constrain SN Ia explosion scenarios. We base our discussion on three-dimensional simulations of a delayed detonation in a Chandrasekhar-mass white dwarf and of a violent merger of two white dwarfs (WDs) - realizations of explosion models appropriate for two of the most widely discussed progenitor channels that may give rise to SNe Ia. Although both models have their shortcomings in reproducing details of the early and near-maximum spectra of SN 2011fe obtained by the Nearby Supernova Factory (SNfactory), the overall match with the observations is reasonable. The level of agreement is slightly better for the merger, in particular around maximum, but a clear preference for one model over the other is still not justified. Observations at late epochs, however, hold promise for discriminating the explosion scenarios in a straightforward way, as a nucleosynthesis effect leads to differences in the Co production. SN 2011fe is close enough to be followed sufficiently long to study this effect. © © 2012 The American Astronomical Society. All rights reserved.

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.

nero- a post-maximum supernova radiation transport code

The interpretation of supernova (SN) spectra is essential for deriving SN ejecta properties such as density and composition, which in turn can tell us about their progenitors and the explosion mechanism. A very large number of atomic processes are important for spectrum formation. Several tools for calculating SN spectra exist, but they mainly focus on the very early or late epochs. The intermediate phase, which requires a non-local thermodynamic equilibrium (NLTE) treatment of radiation transport has rarely been studied. In this paper, we present a new SN radiation transport code, nero, which can look at those epochs. All the atomic processes are treated in full NLTE, under a steady-state assumption. This is a valid approach between roughly 50 and 500days after the explosion depending on SN type. This covers the post-maximum photospheric and the early and the intermediate nebular phase. As a test, we compare nero to the radiation transport code of Jerkstrand, Fransson & Kozma and to the nebular code of Mazzali et al. All three codes have been developed independently and a comparison provides a valuable opportunity to investigate their reliability. Currently, nero is one-dimensional and can be used for predicting spectra of synthetic explosion models or for deriving SN properties by spectral modelling. To demonstrate this, we study the spectra of the 'normal' Type Ia supernova (SN Ia) 2005cf between 50 and 350 days after the explosion and identify most of the common SN Ia line features at post-maximum epochs. © 2011 The Authors Monthly Notices of the Royal Astronomical Society © 2011 RAS.

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

Highly ionized fast accretion disc winds have been suggested as an explanation for a variety of observed absorption and emission features in the X-ray spectra of active galactic nuclei. Simple estimates have suggested that these flows may be massive enough to carry away a significant fraction of the accretion energy and could be involved in creating the link between supermassive black holes and their host galaxies. However, testing these hypotheses, and quantifying the outflow signatures, requires high-quality theoretical spectra for comparison with observations. Here, we describe extensions of our Monte Carlo radiative transfer code that allow us to generate realistic theoretical spectra for a much wider variety of disc wind models than that was possible in our previous work. In particular, we have expanded the range of atomic physics simulated by the code so that L- and M-shell ions can now be included. We have also substantially improved our treatment of both ionization and radiative heating such that we are now able to compute spectra for outflows containing far more diverse plasma conditions. We present example calculations that illustrate the variety of spectral features predicted by parametrized outflow models and demonstrate their applicability to the interpretation of data by comparison with observations of the bright quasar PG1211+143. We find that the major features in the observed 2-10 keV spectrum of this object can be well reproduced by our spectra, confirming that it likely hosts a massive outflow. © 2010 The Authors. Journal compilation © 2010 RAS.

Late-time supernova light curves:The effect of internal conversion and auger electrons

Energy release from radioactive decays contributes significantly to supernova light curves. Previous works, which considered the energy deposited by ?-rays and positrons produced by Ni, Co, Ni, Co, Ti and Sc, have been quite successful in explaining the light curves of both core collapse and thermonuclear supernovae. We point out that Auger and internal conversion electrons, together with the associated X-ray cascade, constitute an additional heat source. When a supernova is transparent to ?-rays, these electrons can contribute significantly to light curves for reasonable nucleosynthetic yields. In particular, the electrons emitted in the decay of Co, which are largely due to internal conversion from a fortuitously low-lying 3/2 state in the daughter Fe, constitute an additional significant energy-deposition channel. We show that when the heating by these electrons is accounted for, a slow-down in the light curve of SN 1998bw is naturally obtained for typical hypernova nucleosynthetic yields. Additionally, we show that for generic Type Ia supernova yields, the Auger electrons emitted in the ground-state to ground-state electron capture decay of Fe exceed the energy released by the Ti decay chain for many years after the explosion. © 2009 RAS.

Nebular emission-line profiles of Type Ib/c supernovae - Probing the ejecta asphericity

In order to assess qualitatively the ejecta geometry of stripped-envelope core-collapse supernovae (SNe), we investigate 98 late-time spectra of 39 objects, many of them previously unpublished. We perform a Gauss-fitting of the [O ] ??6300, 6364 feature in all spectra, with the position, full width at half maximum and intensity of the ?6300 Gaussian as free parameters, and the ?6364 Gaussian added appropriately to account for the doublet nature of the [O ] feature. On the basis of the best-fitting parameters, the objects are organized into morphological classes, and we conclude that at least half of all Type Ib/c SNe must be aspherical. Bipolar jet models do not seem to be universally applicable, as we find too few symmetric double-peaked [O ] profiles. In some objects, the [O ] line exhibits a variety of shifted secondary peaks or shoulders, interpreted as blobs of matter ejected at high velocity and possibly accompanied by neutron-star kicks to assure momentum conservation. At phases earlier than ~200 d, a systematic blueshift of the [O ] ??6300, 6364 line centroids can be discerned. Residual opacity provides the most convincing explanation of this phenomenon, photons emitted on the rear side of the SN being scattered or absorbed on their way through the ejecta. Once modified to account for the doublet nature of the oxygen feature, the profile of Mg i] ?4571 at sufficiently late phases generally resembles that of [O ] ??6300, 6364, suggesting negligible contamination from other lines and confirming that O and Mg are similarly distributed within the ejecta. © 2009 RAS.

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.

Corneal biomechanical properties in primary open angle glaucoma and normal tension glaucoma

PURPOSE: To determine and compare the corneal biomechanical properties between eyes with primary open angle glaucoma (POAG) and eyes with normal tension glaucoma (NTG). PATIENTS AND METHODS: Prospective cross-sectional study. Consecutive eligible POAG and NTG patients attending the Glaucoma Clinic had assessment of their corneal biomechanical properties-corneal hysteresis (CH) and corneal resistance factor (CRF)-using the Ocular Response Analyzer by an observer masked to the diagnosis. Exclusion criteria included previous intraocular surgery, corneal pathology, inflammatory connective tissue disease, and refraction of 5-dimensional or over. If both eyes were eligible, then the right eye was used for analysis. The main outcome measures were corneal hysteresis and CRF measurements. Data analysis was performed using the t test and general linear model. RESULTS: Eighty-one patients (80 whites) were analyzed. Forty had NTG, whereas 41 had POAG. Thirty-five were females. There was a statistically significant difference in mean CH (NTG 9.6±1.3 mm Hg; POAG 9.0±1.4 mm Hg; P=0.01), but not in mean CRF (NTG 9.9±1.4; POAG 10.8±1.7; P=0.06). The highest recorded Goldmann applanation intraocular pressure (IOP) was statistically significantly associated with lower CH (P=0.01) and higher CRF (P=0.02). CONCLUSIONS: There was a small but statistically significant difference in the mean CH between POAG and NTG (CH was higher in NTG). The highest recorded Goldmann applanation IOP was also statistically significantly correlated with lower CH and higher CRF, suggesting that alterations to the corneal biomechanical properties may occur as a result of chronic raised IOP in POAG. © 2008 by Lippincott Williams & Wilkins.

Reversal of optic disc cupping after glaucoma surgery analyzed with a scanning laser tomograph

Objective: To detect and quantitate changes in optic nerve morphology after glaucoma surgery using the Heidelberg Retina Tomograph (HRT, Heidelberg Instruments, Heidelberg, Germany). Design: Nonconsecutive observational case series. Participants and Intervention: The authors prospectively enrolled 21 adult patients undergoing incisional glaucoma surgery for progressive glaucoma damage. Quantitative analysis of the optic nerve head by scanning laser tomography and automated perimetry were performed before and after glaucoma surgery. Main Outcome Measures: Changes in optic nerve parameters were subjected to linear regression analysis with respect to percent of postoperative reduction of intraocular pressure (IOP), as well as with respect to age, refraction, preoperative cup:disc ratio, and change in visual field parameters. Results: Seventeen patients had pre- and postoperative images suitable for analysis. Mean IOP at the time of image acquisition before surgery was 30.5 ± 12 mmHg, and after surgery 11.8 ± 5.2 mmHg (mean follow-up, 26 ± 7 weeks). Eleven of 13 (85%) patients having IOP reduction of greater than 40% showed improvement in optic disc parameters. All four patients with less than 25% reduction in IOP showed worsening of most parameters. Changes in optic disc parameters were highly correlated with percent IOP reduction and with age. The parameters in which change most strongly correlated with percent change of IOP were cup area, rim area, cup:disc ratio, and mean cup depth (each, P <0.005). The age of the patient correlated highly with change in maximum cup depth (P <0.005). Refraction and clinically determined cup:disc ratio correlated poorly with changes in measured optic disc parameters. Clinical improvement in visual fields was correlated with the degree of improvement of cup:disc ratio (P = 0.025). Conclusion: Most patients showing a 40% lowering of IOP after glaucoma surgery show improved optic nerve morphology as measured by the HRT. The amount of improvement correlated highly with the percent reduction of IOP.

Galactic Bulge PNe:Carbon molecules in oxygen-rich environments

Galactic bulge planetary nebulae show evidence of mixed chemistry with emission from both silicate dust and PAHs. This mixed chemistry is unlikely to be related to carbon dredge up, as third dredge-up is not expected to occur in the low mass bulge stars. We show that the phenomenon is widespread, and is seen in 30 nebulae out of our sample of 40. A strong correlation is found between strength of the PAH bands and morphology, in particular, the presence of a dense torus. A chemical model is presented which shows that hydrocarbon chains can form within oxygen-rich gas through gas-phase chemical reactions. We conclude that the mixed chemistry phenomenon occurring in the galactic bulge planetary nebulae is best explained through hydrocarbon chemistry in an UV-irradiated, dense torus. © 2012 International Astronomical Union.

A filamentation instability for streaming cosmic rays

We demonstrate that cosmic rays form filamentary structures in the precursors of supernova remnant shocks due to their self-generated magnetic fields. The cosmic ray filamentation results in the growth of a long-wavelength instability, and naturally couples the rapid non-linear amplification on small scales to larger length-scales. Hybrid magnetohydrodynamics-particle simulations are performed to confirm the effect. The resulting large-scale magnetic field may facilitate the scattering of high-energy cosmic rays as required to accelerate protons beyond the knee in the cosmic ray spectrum at supernova remnant shocks. Filamentation far upstream of the shock may also assist in the escape of cosmic rays from the accelerator.

Cosmic-ray acceleration and escape from supernova remnants

Galactic cosmic-ray (CR) acceleration to the knee in the spectrum at a few PeV is only possible if the magnetic field ahead of a supernova remnant (SNR) shock is strongly amplified by CRs escaping the SNR. A model formulated in terms of the electric charge carried by escaping CRs predicts the maximum CR energy and the energy spectrum of CRs released into the surrounding medium. We find that historical SNRs such as Cas A, Tycho and Kepler may be expanding too slowly to accelerate CRs to the knee at the present time.

Universal behaviour of shock precursors in the presence of efficient cosmic ray acceleration

The self-consistent interaction between energetic particles and self-generated hydromagnetic waves in a cosmic ray pressure dominated plasma is considered. Using a three-dimensional hybrid magnetohydrodynamics (MHD)-kinetic code, which utilizes a spherical harmonic expansion of the Vlasov-Fokker-Planck equation, high-resolution simulations of the magnetic field growth including feedback on the cosmic rays are carried out. It is found that for shocks with high cosmic ray acceleration efficiency, the magnetic fields become highly disorganized, resulting in near isotropic diffusion, independent of the initial orientation of the ambient magnetic field. The possibility of sub-Bohm diffusion is demonstrated for parallel shocks, while the diffusion coefficient approaches the Bohm limit from below for oblique shocks. This universal behaviour suggests that Bohm diffusion in the root-mean-squared field inferred from observation may provide a realistic estimate for the maximum energy acceleration time-scale in young supernova remnants. Although disordered, the magnetic field is not self-similar suggesting a non-uniform energy-dependent behaviour of the energetic particle transport in the precursor. Possible indirect radiative signatures of cosmic ray driven magnetic field amplification are discussed.