Modified Einstein's gravity as a possible missing link between sub- and super-Chandrasekhar type Ia supernovae

We explore the effect of modification to Einstein's gravity in white dwarfs for the first time in the literature, to the best of our knowledge. This leads to significantly sub- and super-Chandrasekhar limiting masses of white dwarfs, determined by a single model parameter. On the other hand, type Ia supernovae (SNeIa), a key to unravel the evolutionary history of the universe, are believed to be triggered in white dwarfs having mass close to the Chandrasekhar limit. However, observations of several peculiar, under- and over-luminous SNeIa argue for exploding masses widely different from this limit. We argue that explosions of the modified gravity induced sub- and super-Chandrasekhar limiting mass white dwarfs result in under- and over-luminous SNeIa respectively, thus unifying these two apparently disjoint sub-classes and, hence, serving as a missing link. Our discovery raises two fundamental questions. Is the Chandrasekhar limit unique? Is Einstein's gravity the ultimate theory for understanding astronomical phenomena? Both the answers appear to be no!

Dynamical stability of nascent neutron stars

This thesis presents a study of the dynamical stability of nascent neutron stars resulting from the accretion induced collapse of rapidly rotating white dwarfs.

Chapter 2 and part of Chapter 3 study the equilibrium models for these neutron stars. They are constructed by assuming that the neutron stars have the same masses, angular momenta, and specific angular momentum distributions as the pre-collapse white dwarfs. If the pre-collapse white dwarf is rapidly rotating, the collapsed object will contain a high density central core of size about 20 km, surrounded by a massive accretion torus extending to hundreds of kilometers from the rotation axis. The ratio of the rotational kinetic energy to gravitational binding energy, β, of these neutron stars is all found to be less than 0.27.

Chapter 3 studies the dynamical stability of these neutron stars by numerically evolving the linearized hydrodynamical equations. A dynamical bar-mode instability is observed when the β of the star is greater than the critical value β_d ≈ 0.25. It is expected that the unstable mode will persist until a substantial amount of angular momentum is carried away by gravitational radiation. The detectability of these sources is studied and it is estimated that LIGO II is unlikely to detect them unless the event rate is greater than 10^-6/year/galaxy.

All the calculations on the structure and stability of the neutron stars in Chapters 2 and 3 are carried out using Newtonian hydrodynamics and gravity. Chapter 4 studies the relativistic effects on the structure of these neutron stars. New techniques are developed and used to construct neutron star models to the first post-Newtonian (1PN) order. The structures of the 1PN models are qualitatively similar to the corresponding Newtonian models, but the values of β are somewhat smaller. The maximum β for these 1PN neutron stars is found to be 0.24, which is 8% smaller than the Newtonian result (0.26). However, relativistic effects will also change the critical value β_d. A detailed post-Newtonian stability analysis has yet to be carried out to study the relativistic effects on the dynamical stability of these neutron stars.

Fully nonlinear ion-sound waves in a dense Fermi magnetoplasma

The nonlinear propagation of ion-sound waves in a collisionless dense electron-ion magnetoplasma is investigated. The inertialess electrons are assumed to follow a non-Boltzmann distribution due to the pressure for the Fermi plasma and the ions are described by the hydrodynamic (HD) equations. An energy balance-like equation involving a new Sagdeev-type pseudo-potential is derived in the presence of the quantum statistical effects. Numerical calculations reveal that the profiles of the Sagdeev-like potential and the ion-sound density excitations are significantly affected by the wave direction cosine and the Mach number. The present studies might be helpful to understand the excitation of nonlinear ion-sound waves in dense plasmas such as those in superdense white dwarfs and neutron stars as well as in intense laser-solid density plasma experiments.

A Brown Dwarf Mass Donor in an Accreting Binary

A long-standing and unverified prediction of binary star evolution theory is the existence of a population of white dwarfs accreting from substellar donor stars. Such systems ought to be common, but the difficulty of finding them, combined with the challenge of detecting the donor against the light from accretion, means that no donor star to date has a measured mass below the hydrogen burning limit. We applied a technique that allowed us to reliably measure the mass of the unseen donor star in eclipsing systems. We were able to identify a brown dwarf donor star, with a mass of 0.052 ± 0.002 solar mass. The relatively high mass of the donor star for its orbital period suggests that current evolutionary models may underestimate the radii of brown dwarfs.

ULTRACAM photometry of the eclipsing cataclysmic variable OU Vir

We present high-speed, three-colour photometry of the faint eclipsing cataclysmic variable OU Vir. For the first time in OU Vir, separate eclipses of the white dwarf and the bright spot have been observed. We use timings of these eclipses to derive a purely photometric model of the system, obtaining a mass ratio of q=0.175+/-0.025, an inclination of i=79.degrees2+/-0.degrees7 and a disc radius of R-d/a=0.2315+/-0.0150. We separate the white dwarf eclipse from the light curve and, by fitting a blackbody spectrum to its flux in each passband, obtain a white dwarf temperature of T=13900+/-600 K and a distance of D=51+/-17 pc. Assuming that the primary obeys the Nauenberg mass-radius relation for white dwarfs and allowing for temperature effects, we also find a primary mass M-w/M-circle dot=0.89+/-0.20, a primary radius R-w/R-circle dot=0.0097+/-0.0031 and an orbital separation a/R-circle dot=0.74+/-0.05.

Pulsating or not? A search for hidden pulsations below the red edge of the ZZ Ceti instability strip

The location of the red edge of the ZZ Ceti instability strip is defined observationally as being the lowest temperature for which a white dwarf with a H-rich atmosphere (DA) is known to exhibit periodic brightness variations. Whether this cut-off in flux variations is actually due to a cessation of pulsation or merely due to the attenuation of any variations by the convection zone, rendering them invisible, is not clear. The latter is a theoretical possibility because with decreasing effective temperature, the emergent flux variations become an ever smaller fraction of the amplitude of the flux variations in the interior. In contrast to the flux variations, the visibility of the velocity variations associated with the pulsations is not thought to be similarly affected. Thus, models imply that were it still pulsating, a white dwarf just below the observed red edge should show velocity variations. In order to test this possibility, we used time-resolved spectra of three DA white dwarfs that do not show photometric variability, but which have derived temperatures only slightly lower than the coolest ZZ Ceti variables. We find that none of our three targets show significant periodic velocity variations, and set 95% confidence limits on amplitudes of 3.0, 5.2, and 8.8 km s(-1). Thus, for two out of our three objects, we can rule out velocity variations as large as 5.4 km s(-1) observed for the strongest mode in the cool white dwarf pulsator ZZ Psc. In order to verify our procedures, we also examined similar data of a known ZZ Ceti, HL Tau 76. Applying external information from the light curve, we detect significant velocity variations for this object with amplitudes of up to 4 km s(-1). Our results suggest that substantial numbers of pulsators having large velocity amplitudes do not exist below the observed photometric red edge and that the latter probably reflects a real termination of pulsations.

WD 0346+246: A very low luminosity, cool degenerate in Taurus

We report the serendipitous discovery of a new, very low luminosity, cool degenerate in the region of Taurus. The object was found as a very high proper-motion star (mu = 1.'' 3 yr(-1)) on seven I-band UK Schmidt Telescope plates, dating from 1987 to 1994, via digitized scans from the new, fast, high-precision microdensitometer SuperCOSMOS. Photometry and spectrophotometry indicate that the object has a temperature comparable to those of the handful of coolest white dwarfs currently known (T similar to 3900 K). We discuss the relevance of this discovery to current research concerning Galactic structure and evolution.

Delay times and rates for Type Ia supernovae and thermonuclear explosions from double-detonation sub-Chandrasekhar mass models

We present theoretical delay times and rates of thermonuclear explosions that are thought to produce Type Ia supernovae (SNe Ia), including the double-detonation sub-Chandrasekhar mass model, using the population synthesis binary evolution code startrack. If detonations of sub-Chandrasekhar mass carbon-oxygen white dwarfs following a detonation in an accumulated layer of helium on the white dwarf's surface ('double-detonation' models) are able to produce thermonuclear explosions which are characteristically similar to those of SNe Ia, then these sub-Chandrasekhar mass explosions may account for at least some substantial fraction of the observed SN Ia rate. Regardless of whether all double-detonations look like 'normal' SNe Ia, in any case the explosions are expected to be bright and thus potentially detectable. Additionally, we find that the delay time distribution of double-detonation sub-Chandrasekhar mass SNe Ia can be divided into two distinct formation channels: the 'prompt' helium-star channel with delay times

Double-detonation sub-Chandrasekhar supernovae: Can minimum helium shell masses detonate the core?

The explosion of sub-Chandrasekhar mass white dwarfs via the double detonation scenario is a potential explanation for type Ia supernovae. In this scenario, a surface detonation in a helium layer initiates a detonation in the underlying carbon/oxygen core leading to an explosion. For a given core mass, a lower bound has been determined on the mass of the helium shell required for dynamical burning during a helium flash, which is a necessary prerequisite for detonation. For a range of core and corresponding minimum helium shell masses, we investigate whether an assumed surface helium detonation is capable of triggering a subsequent detonation in the core even for this limiting case. We carried out hydrodynamic simulations on a co-expanding Eulerian grid in two dimensions assuming rotational symmetry. The detonations are propagated using the level-set approach and a simplified scheme for nuclear reactions that has been calibrated with a large nuclear network. The same network is used to determine detailed nucleosynthetic abundances in a post-processing step. Based on approximate detonation initiation criteria in the literature, we find that secondary core detonations are triggered for all of the simulated models, ranging in core mass from 0.810 up to 1.385 M? with corresponding shell masses from 0.126 down to 0.0035 M?. This implies that, as soon as a detonation triggers in a helium shell covering a carbon/oxygen white dwarf, a subsequent core detonation is virtually inevitable.

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.

Modeling Type Ia supernova explosions

Despite their astrophysical significanceas a major contributor to cosmic nucleosynthesis and as distance indicators in observational cosmologyType Ia supernovae lack theoretical explanation. Not only is the explosion mechanism complex due to the interaction of (potentially turbulent) hydrodynamics and nuclear reactions, but even the initial conditions for the explosion are unknown. Various progenitor scenarios have been proposed. After summarizing some general aspects of Type Ia supernova modeling, recent simulations of our group are discussed. With a sequence of modeling starting (in some cases) from the progenitor evolution and following the explosion hydrodynamics and nucleosynthesis we connect to the formation of the observables through radiation transport in the ejecta cloud. This allows us to analyze several models and to compare their outcomes with observations. While pure deflagrations of Chandrasekhar-mass white dwarfs and violent mergers of two white dwarfs lead to peculiar events (that may, however, find their correspondence in the observed sample of SNe Ia), only delayed detonations in Chandrasekhar-mass white dwarfs or sub-Chandrasekhar-mass explosions remain promising candidates for explaining normal Type Ia supernovae. © 2011 Elsevier B.V. All rights reserved.

CI Aql: a Type Ia supernova progenitor?

If recurrent novae are progenitors of Type Ia supernovae, their white dwarfs must have masses close to the Chandrasekhar limit. The most reliable means of determining white dwarf masses in recurrent novae is dynamically, via radial-velocity and rotational-broadening measurements of the companion star. Such measurements require the system to be both eclipsing and to show absorption features from the secondary star. Prior to the work reported here, the only dynamical mass estimate of a recurrent nova was for U Sco, which has a white dwarf mass of 1.55 +/- 0.24 Msolar (Thoroughgood et al. 2001). We present new time-resolved, intermediate-resolution spectroscopy of the eclipsing recurrent nova CI Aquilae (CI Aql) during quiescence. We find the mass of the white dwarf to be 1.00 +/- 0.14 Msolar and the mass of the secondary star to be 2.32 +/- 0.19 Msolar. We estimate the radius of the secondary to be 2.07 +/- 0.06 Rsolar, implying that it is a slightly-evolved early A-type star. The high mass ratio of q = 2.35 +/- 0.24 and the high secondary-star mass implies that the mass transfer occurs on a thermal timescale. We suggest that CI Aql is rapidly evolving into a supersoft X-ray source, and ultimately may explode as a Type Ia supernova within 10 Myr.

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.

OGLE-2013-SN-079: A Lonely Supernova Consistent with a Helium Shell Detonation

We present observational data for a peculiar supernova discovered by the OGLE-IV survey and followed by the Public ESO Spectroscopic Survey for Transient Objects. The inferred redshift of z = 0.07 implies an absolute magnitude in the rest-frame I-band of M-1 similar to -17.6 mag. This places it in the luminosity range between normal Type Ia SNe and novae. Optical and near infrared spectroscopy reveal mostly Ti and Ca lines, and an unusually red color arising from strong depression of flux at rest wavelengths

Early ultraviolet emission in the Type Ia supernova LSQ12gdj: No evidence for ongoing shock interaction

We present photospheric-phase observations of LSQ12gdj, a slowly declining, UV-bright Type Ia supernova. Classified well before maximum light, LSQ12gdj has extinction-corrected absolute magnitude MB = -19.8, and pre-maximum spectroscopic evolution similar to SN 1991T and the super-Chandrasekhar-mass SN 2007if. We use ultraviolet photometry from Swift, ground-based optical photometry, and corrections from a near-infrared photometric template to construct the bolometric (1600-23 800 Å) light curve out to 45 d past B-band maximum light. We estimate that LSQ12gdj produced 0.96 ± 0.07 M· of 56Ni, with an ejected mass near or slightly above the Chandrasekhar mass. As much as 27 per cent of the flux at the earliest observed phases, and 17 per cent at maximum light, is emitted bluewards of 3300 Å. The absence of excess luminosity at late times, the cutoff of the spectral energy distribution bluewards of 3000 Å and the absence of narrow line emission and strong Na I D absorption all argue against a significant contribution from ongoing shock interaction. However, ~10 per cent of LSQ12gdj's luminosity near maximum light could be produced by the release of trapped radiation, including kinetic energy thermalized during a brief interaction with a compact, hydrogen-poor envelope (radius <1013 cm) shortly after explosion; such an envelope arises generically in double-degenerate merger scenarios.