Three-dimensional pure deflagration models with nucleosynthesis and synthetic observables for type ia supernovae

We investigate whether pure deflagration models ofChandrasekhar-mass carbon-oxygen white dwarf stars can account for one or more subclass of the observed population of Type Ia supernova (SN Ia) explosions. We compute a set of 3D full-star hydrodynamic explosion models, in which the deflagration strength is parametrized using the multispot ignition approach. For each model, we calculate detailed nucleosynthesis yields in a post-processing step with a 384 nuclide nuclear network. We also compute synthetic observables with our 3D Monte Carlo radiative transfer code for comparison with observations. For weak and intermediate deflagration strengths (energy release E {less-than or approximate} 1.1 × 10 erg), we find that the explosion leaves behind a bound remnant enriched with 3 to 10 per cent (by mass) of deflagration ashes. However, we do not obtain the large kick velocities recently reported in the literature. We find that weak deflagrations with E ~ 0.5 × 10 erg fit well both the light curves and spectra of 2002cx-like SNe Ia, and models with even lower explosion energies could explain some of the fainter members of this subclass. By comparing our synthetic observables with the properties of SNe Ia, we can exclude the brightest, most vigorously ignited models as candidates for any observed class of SN Ia: their B-V colours deviate significantly from both normal and 2002cx-like SNe Ia and they are too bright to be candidates for other subclasses.

A spectral synthesis code for rapid modelling of supernovae

We present TARDIS-an open-source code for rapid spectral modelling of supernovae (SNe). Our goal is to develop a tool that is sufficiently fast to allow exploration of the complex parameter spaces of models for SN ejecta. This can be used to analyse the growing number of highquality SN spectra being obtained by transient surveys. The code uses Monte Carlo methods to obtain a self-consistent description of the plasma state and to compute a synthetic spectrum. It has a modular design to facilitate the implementation of a range of physical approximations that can be compared to assess both accuracy and computational expediency. This will allow users to choose a level of sophistication appropriate for their application. Here, we describe the operation of the code and make comparisons with alternative radiative transfer codes of differing levels of complexity (SYN++, PYTHON and ARTIS). We then explore the consequence of adopting simple prescriptions for the calculation of atomic excitation, focusing on four species of relevance to Type Ia SN spectra-Si II, SII, MgII and Ca II. We also investigate the influence of three methods for treating line interactions on our synthetic spectra and the need for accurate radiative rate estimates in our scheme.

Type Ia supernova bolometric light curves and ejected mass estimates from the nearby supernova factory

We present a sample of normal Type Ia supernovae (SNe Ia) from the Nearby Supernova Factory data set with spectrophotometry at sufficiently late phases to estimate the ejected mass using the bolometric light curve.Wemeasure Ni masses from the peak bolometric luminosity, then compare the luminosity in the Co-decay tail to the expected rate of radioactive energy release from ejecta of a given mass. We infer the ejected mass in a Bayesian context using a semi-analytic model of the ejecta, incorporating constraints from contemporary numerical models as priors on the density structure and distribution of Ni throughout the ejecta. We find a strong correlation between ejected mass and light-curve decline rate, and consequently Ni mass, with ejected masses in our data ranging from 0.9 to 1.4 M. Most fast-declining (SALT2 x <-1) normal SNe Ia have significantly sub-Chandrasekhar ejected masses in our fiducial analysis.

The effect of helium accretion efficiency on rates of Type Ia supernovae: Double detonations in accreting binaries

he double-detonation explosion scenario of Type Ia supernovae (SNe Ia) has gained increased support from the SN Ia community as a viable progenitor model, making it a promising candidate alongside the well-known single degenerate and double degenerate scenarios. We present delay times of double-detonation SNe, in which a sub-Chandrasekhar mass carbon–oxygen white dwarf (WD) accretes non-dynamically from a helium-rich companion. One of the main uncertainties in quantifying SN rates from double detonations is the (assumed) retention efficiency of He-rich matter. Therefore, we implement a new prescription for the treatment of accretion/accumulation of He-rich matter on WDs. In addition, we test how the results change depending on which criteria are assumed to lead to a detonation in the helium shell. In comparing the results to our standard case (Ruiter et al.), we find that regardless of the adopted He accretion prescription, the SN rates are reduced by only ∼25 per cent if low-mass He shells (≲0.05 M⊙) are sufficient to trigger the detonations. If more massive (0.1 M⊙) shells are needed, the rates decrease by 85 per cent and the delay time distribution is significantly changed in the new accretion model – only SNe with prompt (<500 Myr) delay times are produced. Since theoretical arguments favour low-mass He shells for normal double-detonation SNe, we conclude that the rates from double detonations are likely to be high, and should not critically depend on the adopted prescription for accretion of He.

Synthetic light curves and spectra for three-dimensional delayed-detonation models of type ia supernovae

In a companion paper, Seitenzahl et al. have presented a set of three-dimensional delayed detonation models for thermonuclear explosions of near-Chandrasekhar-mass white dwarfs (WDs). Here,we present multidimensional radiative transfer simulations that provide synthetic light curves and spectra for those models. The model sequence explores both changes in the strength of the deflagration phase (which is controlled by the ignition configuration in our models) and the WD central density. In agreement with previous studies, we find that the strength of the deflagration significantly affects the explosion and the observables. Variations in the central density also have an influence on both brightness and colour, but overall it is a secondary parameter in our set of models. In many respects, the models yield a good match to the observed properties of normal Type Ia supernovae (SNe Ia): peak brightness, rise/decline time-scales and synthetic spectra are all in reasonable agreement. There are, however, several differences. In particular, the models are systematically too red around maximum light, manifest spectral line velocities that are a little too high and yield I-band light curves that do not match observations. Although some of these discrepancies may simply relate to approximations made in the modelling, some pose real challenges to the models. If viewed as a complete sequence, our models do not reproduce the observed light-curve width- luminosity relation (WLR) of SNe Ia: all our models show rather similar B-band decline rates, irrespective of peak brightness. This suggests that simple variations in the strength of the deflagration phase in Chandrasekhar-mass deflagration-to-detonation models do not readily explain the observed diversity of normal SNe Ia. This may imply that some other parameter within the Chandrasekhar-mass paradigm is key to the WLR, or that a substantial fraction of normal SNe Ia arise from an alternative explosion scenario.

'Super-Chandrasekhar' Type Ia Supernovae at nebular epochs

We present a first systematic comparison of superluminous Type Ia supernovae (SNe Ia) at late epochs, including previously unpublished photometric and spectroscopic observations of SN 2007if, SN 2009dc and SNF20080723-012. Photometrically, the objects of our sample show a diverse late-time behaviour, some of them fading quite rapidly after a light-curve break at ∼ 150-200 d. The latter is likely the result of flux redistribution into the infrared, possibly caused by dust formation, rather than a true bolometric effect. Nebular spectra of superluminous SNe Ia are characterized by weak or absent [Fe III] emission, pointing at a low ejecta ionization state as a result of high densities. To constrain the ejecta and Ni masses of superluminous SNe Ia, we compare the observed bolometric light curve of SN 2009dc with synthetic model light curves, focusing on the radioactive tail after ∼60 d. Models with enough Ni to explain the light-curve peak by radioactive decay, and at the same time sufficient mass to keep the ejecta velocities low, fail to reproduce the observed light-curve tail of SN 2009dc because of too much γ -ray trapping.We instead propose a model with ∼1M of Ni and ∼2 M of ejecta, which may be interpreted as the explosion of a Chandrasekhar-mass white dwarf (WD) enshrouded by 0.6-0.7 M of C/O-rich material, as it could result from a merger of two massive C/O WDs. This model reproduces the late light curve of SN 2009dc well. A flux deficit at peak may be compensated by light from the interaction of the ejecta with the surrounding material.

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.

A window on exoplanet dynamical histories: Rossiter-McLaughlin observations of WASP-13b and WASP-32b

We present Rossiter-McLaughlin observations of WASP-13b and WASP-32b and determine the sky-projected angle between the normal of the planetary orbit and the stellar rotation axis (λ). WASP-13b and WASP-32b both have prograde orbits and are consistent with alignment with measured sky-projected angles of λ =8°^{+13}_{-12} and λ =-2°^{+17}_{-19}, respectively. Both WASP-13 and WASP-32 have Teff < 6250 K, and therefore, these systems support the general trend that aligned planetary systems are preferentially found orbiting cool host stars. A Lomb-Scargle periodogram analysis was carried out on archival SuperWASP data for both systems. A statistically significant stellar rotation period detection (above 99.9 per cent confidence) was identified for the WASP-32 system with Prot =11.6 ± 1.0 days. This rotation period is in agreement with the predicted stellar rotation period calculated from the stellar radius, R*, and vsin i if a stellar inclination of i* =90° is assumed. With the determined rotation period, the true 3D angle between the stellar rotation axis and the planetary orbit, ψ, was found to be ψ = 11° ± 14°. We conclude with a discussion on the alignment of systems around cool host stars with Teff < 6150 K by calculating the tidal dissipation time-scale. We find that systems with short tidal dissipation time-scales are preferentially aligned and systems with long tidal dissipation time-scales have a broad range of obliquities.

A high resolution study of complex organic molecules in hot cores

We present the results of a line identification analysis using data from the IRAM Plateau de Bure Inferferometer, focusing on six massive star-forming hot cores: G31.41+0.31, G29.96-0.02, G19.61-0.23, G10.62-0.38, G24.78+0.08A1 and G24.78+0.08A2. We identify several transitions of vibrationally excited methyl formate (HCOOCH$_3$) for the first time in these objects as well as transitions of other complex molecules, including ethyl cyanide (C$_2$H$_5$CN), and isocyanic acid (HNCO). We also postulate a detection of one transition of glycolaldehyde (CH$_2$(OH)CHO) in two new hot cores. We find G29.96-0.02, G19.61-0.23, G24.78+0.08A1 and 24.78+0.08A2 to be chemically very similar. G31.41+0.31, however, is chemically different: it manifests a larger chemical inventory and has significantly larger column densities. We suggest that it may represent a different evolutionary stage to the other hot cores in the sample, or it may surround a star with a higher mass. We derive column densities for methyl formate in G31.41+0.31, using the rotation diagram method, of $\times$10$^{17}$ cm$^{-2}$ and a T$_{rot}$ of $\sim$170 K. For G29.96-0.02, G24.78+0.08A1 and G24.78+0.08A2, glycolaldehyde, methyl formate and methyl cyanide all seem to trace the same material and peak at roughly the same position towards the dust emission peak. For G31.41+0.31, however, glycolaldehyde shows a different distribution to methyl formate and methyl cyanide and seems to trace the densest, most compact inner part of hot cores.

New detections of HC5N towards hot cores associated with 6.7 GHz methanol masers

We present new detections of cyanodiacetylene (HC5N) toward hot molecular cores, observed with the Tidbinbilla 34 m radio telescope (DSS–34). In a sample of 79 hot molecular cores, HC5N was detected towards 35. These results are counter to the expectation that long chain cyanopolyynes, such as HC5N, are not typically found in hot molecular cores, unlike their shorter chain counterpart HC3N. However it is consistent with recent models which suggest HC5N may exist for a limited period during the evolution of hot molecular cores.

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.

Defocused transmission spectroscopy: a potential detection of sodium in the atmosphere of WASP-12b

We report on a pilot study of a novel observing technique, defocussed transmission spectroscopy, and its application to the study of exoplanet atmospheres using ground-based platforms. Similar to defocussed photometry, defocussed transmission spectroscopy has an added advantage over normal spectroscopy in that it reduces systematic errors due to flat-fielding, PSF variations, slit-jaw imperfections and other effects associated with ground-based observations. For one of the planetary systems studied, WASP-12b, we report a tentative detection of additional Na absorption of 0.12+/-0.03[+0.03]% during transit using a 2A wavelength mask. After consideration of a systematic that occurs mid-transit, it is likely that the true depth is actually closer to 0.15%. This is a similar level of absorption reported in the atmosphere of HD209458b (0.135+/-0.017%, Snellen et al. 2008). Finally, we outline methods that will improve the technique during future observations, based on our findings from this pilot study.

Spitzer Space Telescope spectra of post-AGB stars in the large magellanic Cloud-polycyclic aromatic hydrocarbons at low metallicities

This paper reports variations of polycyclic aromatic hydrocarbons (PAHs) features that were found in Spitzer Space Telescope spectra of carbon-rich post-asymptotic giant branch (post-AGB) stars in the Large Magellanic Cloud (LMC). The paper consists of two parts. The first part describes our Spitzer spectral observing programme of 24 stars including post-AGB candidates. The latter half of this paper presents the analysis of PAH features in 20 carbon-rich post-AGB stars in the LMC, assembled from the Spitzer archive as well as from our own programme.We found that five post-AGB stars showed a broad feature with a peak at 7.7 μm, that had not been classified before. Further, the 10-13 μm PAH spectra were classified into four classes, one of which has three broad peaks at 11.3, 12.3 and 13.3 μm rather than two distinct sharp peaks at 11.3 and 12.7 μm, as commonly found in HII regions. Our studies suggest that PAHs are gradually processed while the central stars evolve from post-AGB phase to planetary nebulae, changing their composition before PAHs are incorporated into the interstellar medium. Although some metallicity dependence of PAH spectra exists, the evolutionary state of an object is more significant than its metallicity in determining the spectral characteristics of PAHs for LMC and Galactic post-AGB stars. © 2014 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society.

Superluminous supernovae from PESSTO

We present optical spectra and light curves for three hydrogen-poor superluminous supernovae followed by the Public ESO Spectroscopic Survey of Transient Objects (PESSTO). Time series spectroscopy from a fewdays aftermaximum light to 100 d later shows them to be fairly typical of this class, with spectra dominated by Ca II, MgII, FeII, and Si II, which evolve slowly over most of the post-peak photospheric phase. We determine bolometric light curves and apply simple fitting tools, based on the diffusion of energy input by magnetar spin-down, Ni-56 decay, and collision of the ejecta with an opaque circumstellar shell. We investigate how the heterogeneous light curves of our sample (combined with others from the literature) can help to constrain the possible mechanisms behind these events. We have followed these events to beyond 100-200 d after peak, to disentangle host galaxy light from fading supernova flux and to differentiate between the models, which predict diverse behaviour at this phase. Models powered by radioactivity require unrealistic parameters to reproduce the observed light curves, as found by previous studies. Both magnetar heating and circumstellar interaction still appear to be viable candidates. A large diversity is emerging in observed tail-phase luminosities, with magnetar models failing in some cases to predict the rapid drop in flux. This would suggest either that magnetars are not responsible, or that the X-ray flux from the magnetar wind is not fully trapped. The light curve of one object shows a distinct rebrightening at around 100 d after maximum light. We argue that this could result either from multiple shells of circumstellar material, or from a magnetar ionization front breaking out of the ejecta.

Line-driven radiative outflows in luminous quasars

An analysis of ≃19 500 narrow (≲200 km s^-1) CIV λλ1548.2,1550.8 absorbers in ≃34 000 Sloan Digital Sky Survey quasar spectra is presented. The statistics of the number of absorbers as a function of outflow velocity shows that in approximately two-thirds of outflows, with multiple C IV absorbers present, absorbers are line-locked at the 500 km s^-1 velocity separation of the C IV absorber doublet; appearing as 'triplets' in the quasar spectra. Line-locking is an observational signature of radiative line-driving in outflowing material, where the successive shielding of 'clouds' of material in the outflow locks the clouds together in outflow velocity. Line-locked absorbers are seen in both broad absorption line (BAL) quasars and non-BAL quasars with comparable frequencies and with velocities out to at least 20 000 km s^-1. There are no detectable differences in the absorber properties and the dust content of single C IV doublets and line-locked C IV doublets. The gas associated with both single and line-locked CIV absorption systems includes material with a wide range of ionization potential (14-138 eV). Both single and line-locked CIV absorber systems show strong systematic trends in their ionization as a function of outflow velocity, with ionization decreasing rapidly with increasing outflow velocity. Initial simulations, employing CLOUDY, demonstrate that a rich spectrum of line-locked signals at various velocities may be expected due to significant opacities from resonance lines of Li-, He- and H-like ions of O, C and N, along with contributions from He II and HI resonance lines. The simulations confirm that line-driving can be the dominant acceleration mechanism for clouds with N(H I) ≃ 10¹⁹ cm^-2.