Source author record

Lars Bildsten

Lars Bildsten appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.SR astro-ph.HE astro-ph.CO astro-ph astro-ph.GA astro-ph.IM gr-qc physics.flu-dyn

Catalog footprint

What is connected

69works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

3D Hydrodynamics of Pre-supernova Outbursts in Convective Red Supergiant Envelopes

Eruptive mass loss likely produces the energetic outbursts observed from some massive stars before they undergo core-collapse supernovae (CCSNe). The resulting dense circumstellar medium (CSM) may also cause the subsequent SNe to be observed as Type IIn events. The leading hypothesis of the cause of these outbursts is the response of the envelope of the red supergiant (RSG) progenitor to energy deposition in the months to years prior to collapse. Early theoretical studies of this phenomena were limited to 1D, leaving the 3D convective RSG structure unaddressed. Using FLASH's hydrodynamic capabilities, we explore the 3D outcomes by constructing convective RSG envelope models and depositing energies less than the envelope binding energies on timescales shorter than the envelope dynamical time deep within them. We confirm the 1D prediction of an outward moving acoustic pulse steepening into a shock, unbinding the outermost parts of the envelope. However, we find that the initial 2-4 km/s convective motions seed the intrinsic convective instability associated with the high entropy material deep in the envelope, enabling gas from deep within the envelope to escape, increasing the amount of ejected mass compared to an initially "quiescent" envelope. The 3D models reveal a rich density structure, with column densities varying by 10x along different lines of sight. Our work highlights that the 3D convective nature of RSG envelopes impacts our ability to reliably predict the outburst dynamics, the amount, and the spatial distribution of the ejected mass associated with deep energy deposition.

preprint2022arXiv

Discovery of a double detonation thermonuclear supernova progenitor

We present the discovery of a new double detonation progenitor system consisting of a hot subdwarf B (sdB) binary with a white dwarf companion with an P=76.34179(2) min orbital period. Spectroscopic observations are consistent with an sdB star during helium core burning residing on the extreme horizontal branch. Chimera light curves are dominated by ellipsoidal deformation of the sdB star and a weak eclipse of the companion white dwarf. Combining spectroscopic and light curve fits we find a low mass sdB star, $M_{\rm sdB}=0.383\pm0.028$ M$_\odot$ with a massive white dwarf companion, $M_{\rm WD}=0.725\pm0.026$ M$_\odot$. From the eclipses we find a blackbody temperature for the white dwarf of 26,800 K resulting in a cooling age of $\approx$25 Myrs whereas our MESA model predicts an sdB age of $\approx$170 Myrs. We conclude that the sdB formed first through stable mass transfer followed by a common envelope which led to the formation of the white dwarf companion $\approx$25 Myrs ago. Using the MESA stellar evolutionary code we find that the sdB star will start mass transfer in $\approx$6 Myrs and in $\approx$60 Myrs the white dwarf will reach a total mass of $0.92$ M$_\odot$ with a thick helium layer of $0.17$ M$_\odot$. This will lead to a detonation that will likely destroy the white dwarf in a peculiar thermonuclear supernova. PTF1 2238+7430 is only the second confirmed candidate for a double detonation thermonuclear supernova. Using both systems we estimate that at least $\approx$1% of white dwarf thermonuclear supernovae originate from sdB+WD binaries with thick helium layers, consistent with the small number of observed peculiar thermonuclear explosions.

preprint2022arXiv

Numerical Simulations of Convective 3-Dimensional Red Supergiant Envelopes

We explore the three-dimensional properties of convective, luminous ($L\approx10^{4.5}-10^{5}L_\odot$), Hydrogen-rich envelopes of Red Supergiants (RSGs) based on radiation hydrodynamic simulations in spherical geometry using $\texttt{Athena++}$. These computations comprise $\approx30\%$ of the stellar volume, include gas and radiation pressure, and self-consistently track the gravitational potential for the outer $\approx 3M_\odot$ of the simulated $M\approx15M_\odot$ stars. This work reveals a radius, $R_\mathrm{corr}$, around which the nature of the convection changes. For $r>R_\mathrm{corr}$, though still optically thick, diffusion of photons dominates the energy transport. Such a regime is well-studied in less luminous stars, but in RSGs, the near- (or above-) Eddington luminosity (due to opacity enhancements at ionization transitions) leads to the unusual outcome of denser regions moving outwards rather than inward. This region of the star also has a large amount of turbulent pressure, yielding a density structure much more extended than 1D stellar evolution predicts. This "halo" of material will impact predictions for both shock breakout and early lightcurves of Type II-P supernovae. Inside of $R_\mathrm{corr}$, we find a nearly flat entropy profile as expected in the efficient regime of mixing-length-theory (MLT). Radiation pressure provides $\approx1/3$ of the support against gravity in this region. Our comparisons to MLT suggest a mixing length of $α=3-4$, consistent with the sizes of convective plumes seen in the simulations. The temporal variability of these 3D models is mostly on the timescale of the convective plume lifetimes ($\approx300$ days), with amplitudes consistent with those observed photometrically.

preprint2022arXiv

Shock Breakout in 3-Dimensional Red Supergiant Envelopes

Using Athena++, we perform 3D Radiation-Hydrodynamic calculations of the radiative breakout of the shock wave in the outer envelope of a red supergiant (RSG) which has suffered core collapse and will become a Type IIP supernova. The intrinsically 3D structure of the fully convective RSG envelope yields key differences in the brightness and duration of the shock breakout (SBO) from that predicted in a 1D stellar model. First, the lower-density `halo' of material outside of the traditional photosphere in 3D models leads to a shock breakout at lower densities than 1D models. This would prolong the duration of the shock breakout flash at any given location on the surface to $\approx$1-2 hours. However, we find that the even larger impact is the intrinsically 3D effect associated with large-scale fluctuations in density that cause the shock to break out at different radii at different times. This substantially prolongs the SBO duration to $\approx$3-6 hours and implies a diversity of radiative temperatures, as different patches across the stellar surface are at different stages of their radiative breakout and cooling at any given time. These predicted durations are in better agreement with existing observations of SBO. The longer durations lower the predicted luminosities by a factor of 3-10 ($L_\mathrm{bol}\sim10^{44}\mathrm{erg\ s^{-1}}$), and we derive the new scalings of brightness and duration with explosion energies and stellar properties. These intrinsically 3D properties eliminate the possibility of using observed rise times to measure the stellar radius via light-travel time effects.

preprint2021arXiv

Still Brighter than Pre-Explosion, SN 2012Z Did Not Disappear: Comparing Hubble Space Telescope Observations a Decade Apart

Type Iax supernovae represent the largest class of peculiar white-dwarf supernovae. The type Iax SN~2012Z in NGC 1309 is the only white dwarf supernova with a detected progenitor system in pre-explosion observations. Deep \textit{Hubble Space Telescope} images taken before SN~2012Z show a luminous, blue source that we have interpreted as a helium-star companion (donor) to the exploding white dwarf. We present here late-time \textit{HST} observations taken $\sim$1400 days after the explosion to test this model. We find the SN light curve can empirically be fit by an exponential decay model in magnitude units. The fitted asymptotic brightness is within $10\%$ of our latest measurements and approximately twice the brightness of the pre-explosion source. The decline of the light curve is too slow to be powered by $^{56}$Co or $^{57}$Co decay: if radioactive decay is the dominate power source, it must be from longer half-life species like $^{55}$Fe. Interaction with circumstellar material may contribute to the light curve, as may shock heating of the companion star. Companion-star models underpredict the observed flux in the optical, producing most of their flux in the UV at these epochs. A radioactively-heated bound remnant, left after only a partial disruption of the white dwarf, is also capable of producing the observed excess late-time flux. Our analysis suggests that the total ejecta + remnant mass is consistent with the Chandrasekhar mass for a range of type Iax supernovae.

preprint2021arXiv

Stochastic Low Frequency Variability in 3-Dimensional Radiation Hydrodynamical Models of Massive Star Envelopes

Increasing main sequence stellar luminosity with stellar mass leads to the eventual dominance of radiation pressure in stellar envelope hydrostatic balance. As the luminosity approaches the Eddington limit, additional instabilities (beyond conventional convection) can occur. These instabilities readily manifest in the outer envelopes of OB stars, where the opacity increase associated with iron yields density and gas pressure inversions in 1D models. Additionally, recent photometric surveys (e.g. TESS) have detected excess broadband low frequency variability in power spectra of OB star lightcurves, called stochastic low frequency variability (SLFV). This motivates our novel 3D Athena++ radiation hydrodynamical (RHD) simulations of two 35$\,$M$_\odot$ star envelopes (the outer $\approx$15$\%$ of the stellar radial extent), one on the zero-age main sequence and the other in the middle of the main sequence. Both models exhibit turbulent motion far above and below the conventional iron opacity peak convection zone (FeCZ), obliterating any ``quiet" part of the near-surface region and leading to velocities at the photosphere of 10-100$\,$km$\,$s$^{-1}$, directly agreeing with spectroscopic data. Surface turbulence also produces SLFV in model lightcurves with amplitudes and power-law slopes that are strikingly similar to those of observed stars. The characteristic frequencies associated with SLFV in our models are comparable to the thermal time in the FeCZ ($\approx$3-7$\,$days$^{-1}$). These simulations, which have no free parameters, are directly validated by observations and, though more models are needed, we remain optimistic that 3D RHD models of main sequence O star envelopes exhibit SLFV originating from the FeCZ.

preprint2020arXiv

A Massive Star's Dying Breaths: Pulsating Red Supergiants and Their Resulting Type IIP Supernovae

Massive stars undergo fundamental-mode and first-overtone radial pulsations with periods of 100-1000 days as Red Supergiants (RSGs). At large amplitudes, these pulsations substantially modify the outer envelope's density structure encountered by the outgoing shock wave from the eventual core collapse of these $M>9M_\odot$ stars. Using Modules for Experiments in Stellar Astrophysics (MESA), we model the effects of fundamental-mode and first-overtone pulsations in the RSG envelopes, and the resulting Type IIP supernovae (SNe) using MESA+STELLA. We find that, in the case of fundamental mode pulsations, SN plateau observables such as the luminosity at day 50, $L_{50}$, time-integrated shock energy $ET$, and plateau duration $t_{\rm p}$ are consistent with radial scalings derived considering explosions of non-pulsating stars. Namely, most of the effect of the pulsation is consistent with the behavior expected for a star of a different size at the time of explosion. However, in the case of overtone pulsations, the Lagrangian displacement is not monotonic. Therefore, in such cases, excessively bright or faint SN emission at different times reflects the underdense or overdense structure of the emitting region near the SN photosphere.

preprint2020arXiv

A new class of Roche lobe-filling hot subdwarf binaries

We present the discovery of the second binary with a Roche lobe-filling hot subdwarf transferring mass to a white dwarf (WD) companion. This 56 minute binary was discovered using data from the Zwicky Transient Facility. Spectroscopic observations reveal an He-sdOB star with an effective temperature of $T_{\rm eff}=33,700\pm1000$ K and a surface gravity of $log(g)=5.54\pm0.11$. The GTC+HiPERCAM light curve is dominated by the ellipsoidal deformation of the He-sdOB star and shows an eclipse of the He-sdOB by an accretion disk as well as a weak eclipse of the WD. We infer a He-sdOB mass of $M_{\rm sdOB}=0.41\pm0.04$ M$_\odot$ and a WD mass of $M_{\rm WD}=0.68\pm0.05$ M$_\odot$. The weak eclipses imply a WD black-body temperature of $63,000\pm10,000$ K and a radius $R_{\rm WD}=0.0148\pm0.0020$ M$_\odot$ as expected for a WD of such high temperature. The He-sdOB star is likely undergoing hydrogen shell burning and will continue transferring mass for $\approx1$ Myrs at a rate of $10^{-9} M_\odot {\rm yr}^{-1}$ which is consistent with the high WD temperature. The hot subdwarf will then turn into a WD and the system will merge in $\approx30$ Myrs. We suggest that Galactic reddening could bias discoveries towards preferentially finding Roche lobe-filling systems during the short-lived shell burning phase. Studies using reddening corrected samples should reveal a large population of helium core-burning hot subdwarfs with $T_{\rm eff}\approx25,000$ K in binaries of 60-90 minutes with WDs. Though not yet in contact, these binaries would eventually come into contact through gravitational wave emission and explode as a sub-luminous thermonuclear supernova or evolve into a massive single WD.

preprint2020arXiv

Comparing Moment-Based and Monte Carlo Methods of Radiation Transport Modeling for Type II-Plateau Supernova Light Curves

Time-dependent electromagnetic signatures from core-collapse supernovae are the result of detailed transport of the shock-deposited and radioactively-powered radiation through the stellar ejecta. Due to the complexity of the underlying radiative processes, considerable approximations are made to simplify key aspects of the radiation transport problem. We present a systematic comparison of the moment-based radiation hydrodynamical code STELLA and the Monte Carlo radiation transport code Sedona in the 1D modeling of Type II-Plateau supernovae. Based on explosion models generated from the Modules for Experiments in Stellar Astrophysics (MESA) instrument, we find remarkable agreements in the modeled light curves and the ejecta structure thermal evolution, affirming the fidelity of both radiation transport modeling approaches. The radiative moments computed directly by the Monte Carlo scheme in Sedona also verify the accuracy of the moment-based scheme. We find that the coarse resolutions of the opacity tables and the numerical approximations in STELLA have insignificant impact on the resulting bolometric light curves, making it an efficient tool for the specific task of optical light curve modeling.

preprint2020arXiv

The first ultracompact Roche lobe-filling hot subdwarf binary

We report the discovery of the first short period binary in which a hot subdwarf star (sdOB) fills its Roche lobe and started mass transfer to its companion. The object was discovered as part of a dedicated high-cadence survey of the Galactic Plane named the Zwicky Transient Facility and exhibits a period of $P_{\rm orb}=39.3401(1)$ min, making it the most compact hot subdwarf binary currently known. Spectroscopic observations are consistent with an intermediate He-sdOB star with an effective temperature of $T_{\rm eff}=42,400\pm300$ K and a surface gravity of $\log(g)=5.77\pm0.05$. A high-signal-to noise GTC+HiPERCAM light curve is dominated by the ellipsoidal deformation of the sdOB star and an eclipse of the sdOB by an accretion disk. We infer a low-mass hot subdwarf donor with a mass $M_{\rm sdOB}=0.337\pm0.015$ M$_\odot$ and a white dwarf accretor with a mass $M_{\rm WD}=0.545\pm0.020$ M$_\odot$. Theoretical binary modeling indicates the hot subdwarf formed during a common envelope phase when a $2.5-2.8$ M$_\odot$ star lost its envelope when crossing the Hertzsprung Gap. To match its current $P_{\rm orb}$, $T_{\rm eff}$, $\log(g)$, and masses, we estimate a post-common envelope period of $P_{\rm orb}\approx150$ min, and find the sdOB star is currently undergoing hydrogen shell burning. We estimate that the hot subdwarf will become a white dwarf with a thick helium layer of $\approx0.1$ M$_\odot$ and will merge with its carbon/oxygen white dwarf companion after $\approx17$ Myr and presumably explode as a thermonuclear supernova or form an R CrB star.

preprint2020arXiv

The Value of Progenitor Radius Measurements for Explosion Modeling of Type II-Plateau Supernovae

Using Modules for Experiments in Stellar Astrophysics (MESA)+STELLA, we show that very different physical models can adequately reproduce a specific observed Type II-Plateau Supernova (SN). We consider SN2004A, SN2004et, SN2009ib, SN2017eaw, and SN2017gmr, Nickel-rich ($M_\mathrm{Ni}>0.03M_\odot$) events with bolometric lightcurves and a well-sampled decline from the plateau. These events also have constraints on the progenitor radius, via a progenitor image, or, in the case of SN2017gmr, a radius from fitting shock-cooling models. In general, many explosions spanning the parameter space of progenitors can yield excellent lightcurve and Fe line velocity agreement, demonstrating the success of scaling laws in motivating models which match plateau properties for a given radius and highlighting the degeneracy between plateau luminosity and velocity in models and observed events, which can span over 50% in ejecta mass, radius, and explosion energy. This can help explain disagreements in explosion properties reported for the same event using different model calculations. Our calculations yield explosion properties when combined with pre-explosion progenitor radius measurements or a robust understanding of the outermost $<0.1\,M_\odot$ of material that quantifies the progenitor radius from SN observations a few days after explosion.

preprint2020arXiv

Variability of massive stars in M31 from the Palomar Transient Factory

Using data from the (intermediate) Palomar Transient Factory (iPTF), we characterize the time variability of ~500 massive stars in M31. Our sample is those stars which are spectrally typed by Massey and collaborators, including Luminous Blue Variables, Wolf-Rayets, and warm and cool supergiants. We use the high-cadence, long-baseline (~5 years) data from the iPTF survey, coupled with data-processing tools that model complex features in the light curves. We find widespread photometric (R-band) variability in the upper Hertzsprung Russell diagram (or CMD) with an increasing prevalence of variability with later spectral types. Red stars (V-I>1.5) exhibit larger amplitude fluctuations than their bluer counterparts. We extract a characteristic variability timescale, tch, via wavelet transformations that are sensitive to both continuous and localized fluctuations. Cool supergiants are characterized by longer timescales (>100 days) than the hotter stars. The latter have typical timescales of tens of days but cover a wider range, from our resolution limit of a few days to longer than 100 days timescales. Using a 60-night block of data straddling two nights with a cadence of around 2 minutes, we extracted tch in the range 0.1--10 days with amplitudes of a few percent for 13 stars. Though there is broad agreement between the observed variability characteristics in the different parts of the upper CMD with theoretical predictions, detailed comparison requires models with a more comprehensive treatment of the various physical processes operating in these stars such as pulsation, subsurface convection, and the effect of binary companions.

preprint2019arXiv

Remnants of Subdwarf Helium Donor Stars Ejected from Close Binaries with Thermonuclear Supernovae

Some binary systems composed of a white dwarf (WD) and a hot subdwarf (sdB) helium star will make contact within the helium burning lifetime of the sdB star. The accreted helium on the WD inevitably undergoes a thermonuclear instability, causing a detonation that is expected to transition into the WD core and lead to a thermonuclear supernova while the donor orbits nearby with high velocity. Motivated by the recent discovery of fast-moving objects that occupy unusual locations on the HR diagram, we explore the impact of the thermonuclear supernovae on the donors in this specific double detonation scenario. We use MESA to model the binary up to the moment of detonation, then 3D Athena++ to model the hydrodynamic interaction of the supernova ejecta with the donor star, calculating the amount of mass that is stripped and the entropy deposited in the deep stellar interior by the strong shock that traverses it. We show that these donor remnants are ejected with velocities primarily set by their orbital speeds: $700-900\ \rm km\ s^{-1}$. We model the long-term thermal evolution of remnants by introducing the shock entropy into MESA models. In response to this entropy change, donor remnants expand and brighten for timescales ranging from $10^6-10^8$ years, giving ample time for these runaway stars to be observed in their inflated state before they leave the galaxy. Even after surface layers are stripped, some donors retain enough mass to resume core helium burning and further delay fading for more than $10^8$ years.

preprint2016arXiv

Asteroseismic signatures of evolving internal stellar magnetic fields

Recent asteroseismic analyses have revealed the presence of strong (B $\gtrsim 10^5$ G) magnetic fields in the cores of many red giant stars. Here, we examine the implications of these results for the evolution of stellar magnetic fields, and we make predictions for future observations. Those stars with suppressed dipole modes indicative of strong core fields should exhibit moderate but detectable quadrupole mode suppression. The long magnetic diffusion times within stellar cores ensure that dynamo-generated fields are confined to mass coordinates within the main sequence convective core, and the observed sharp increase in dipole mode suppression rates above $1.5 \, M_\odot$ may be explained by the larger convective core masses and faster rotation of these more massive stars. In clump stars, core fields of $\sim10^5 \, {\rm G}$ can suppress dipole modes, whose visibility should be equal to or less than the visibility of suppressed modes in ascending red giants. High dipole mode suppression rates in low-mass ($M \lesssim 2 \, M_\odot$) clump stars would indicate that magnetic fields generated during the main sequence can withstand subsequent convective phases and survive into the compact remnant phase. Finally, we discuss implications for observed magnetic fields in white dwarfs and neutron stars, as well as the effects of magnetic fields in various types of pulsating stars.

preprint2016arXiv

Carbon Shell or Core Ignitions in White Dwarfs Accreting from Helium Stars

White dwarfs accreting from helium stars can stably burn at the accreted rate and avoid the challenge of mass loss associated with unstable Helium burning that is a concern for many Type Ia supernovae scenarios. We study binaries with helium stars of mass $1.25 M_\odot\le M_{\rm{He}} \le 1.8 M_\odot$, which have lost their hydrogen rich envelopes in an earlier common envelope event and now orbit with periods ($P_{\rm orb}$) of several hours with non-rotating $0.84$ and $1.0 M_\odot$ C/O WDs. The helium stars fill their Roche lobes (RLs) after exhaustion of central helium and donate helium on their thermal timescales (${\sim}10^5$yr). As shown by others, these mass transfer rates coincide with the steady helium burning range for WDs, and grow the WD core up to near the Chandrasekhar mass ($M_{\rm Ch}$) and a core carbon ignition. We show here, however, that many of these scenarios lead to an ignition of hot carbon ashes near the outer edge of the WD and an inward going carbon flame that does not cause an explosive outcome. For $P_{\rm orb} = 3$ hours, $1.0 M_\odot$ C/O WDs with donor masses $M_{\rm He}\gtrsim1.8 M_\odot$ experience a shell carbon ignition, while $M_{\rm He}\lesssim1.3 M_\odot$ will fall below the steady helium burning range and undergo helium flashes before reaching core C ignition. Those with $1.3 M_\odot \lesssim M_{\rm He} \lesssim 1.7 M_\odot$ will experience a core C ignition. We also calculate the retention fraction of accreted helium when the accretion rate leads to recurrent weak helium flashes.

preprint2016arXiv

Late-time Spectroscopy of Type Iax Supernovae

We examine the late-time (t > 200 days after peak brightness) spectra of Type Iax supernovae (SNe Iax), a low-luminosity, low-energy class of thermonuclear stellar explosions observationally similar to, but distinct from, Type Ia supernovae. We present new spectra of SN 2014dt, resulting in the most complete published late-time spectral sequence of a SN Iax. At late times, SNe Iax have generally similar spectra, all with a similar continuum shape and strong forbidden-line emission. However, there is also significant diversity where some late-time SN Iax spectra display narrow P-Cygni features and a continuum indicative of a photosphere in addition to strong narrow forbidden lines, while others have no obvious P-Cygni features, strong broad forbidden lines, and weak narrow forbidden lines. Finally, some SNe Iax have spectra intermediate to these two varieties with weak P-Cygni features and broad/narrow forbidden lines of similar strength. We find that SNe Iax with strong broad forbidden lines also tend to be more luminous and have higher-velocity ejecta at peak brightness. We estimate blackbody and kinematic radii of the late-time photosphere, finding the latter an order of magnitude larger than the former. We propose a two-component model that solves this discrepancy and explains the diversity of the late-time spectra of SNe Iax. In this model, the broad forbidden lines originate from the SN ejecta, while the photosphere, P-Cygni lines, and narrow forbidden lines originate from a wind launched from the remnant of the progenitor white dwarf and is driven by the radioactive decay of newly synthesized material left in the remnant. The relative strength of the two components accounts for the diversity of late-time SN Iax spectra. This model also solves the puzzle of a long-lived photosphere and slow late-time decline of SNe Iax. (Abridged)

preprint2016arXiv

Rapidly Rising Transients in the Supernova - Superluminous Supernova Gap

We present observations of four rapidly rising (t_{rise}~10d) transients with peak luminosities between those of supernovae (SNe) and superluminous SNe (M_{peak}~-20) - one discovered and followed by the Palomar Transient Factory (PTF) and three by the Supernova Legacy Survey (SNLS). The light curves resemble those of SN 2011kl, recently shown to be associated with an ultra-long-duration gamma ray burst (GRB), though no GRB was seen to accompany our SNe. The rapid rise to a luminous peak places these events in a unique part of SN phase space, challenging standard SN emission mechanisms. Spectra of the PTF event formally classify it as a Type II SN due to broad Halpha emission, but an unusual absorption feature, which can be interpreted as either high velocity Halpha (though deeper than in previously known cases) or Si II (as seen in Type Ia SNe), is also observed. We find that existing models of white dwarf detonations, CSM interaction, shock breakout in a wind (or steeper CSM) and magnetar spindown can not readily explain the observations. We consider the possibility that a "Type 1.5 SN" scenario could be the origin of our events. More detailed models for these kinds of transients and more constraining observations of future such events should help better determine their nature.

preprint2016arXiv

Seismology of Rapidly Rotating Accreting White Dwarfs

A number of White Dwarfs (WDs) in cataclysmic binaries have shown brightness variations consistent with non-radial oscillations as observed in isolated WDs. A few objects have been well-characterized with photometric campaigns in the hopes of gleaning information about the mass, spin, and possibly internal structural characteristics. The novel aspect of this work is the possiblity to measure or constrain the interior structure and spin rate of WDs which have spent gigayears accreting material from their companion, undergoing thousands of nova outbursts in the process. In addition, variations in the surface temperature affect the site of mode driving, and provide unique and challenging tests for mode driving theories previously applied to isolated WD's. Having undergone long-term accretion, these WDs are expected to have been spun up. Spin periods in the range 60-100 seconds have been measured by other means for two objects, GW Lib and V455 And. Compared to typical mode frequencies, the spin frequency may be similar or higher, and the Coriolis force can no longer be treated as a small perturbation on the fluid motions. We present the results of a non-perturbative calculation of the normal modes of these WDs, using interior thermal structures appropriate to accreting systems. This includes a discussion of the surface brightness distributions, which are strongly modified from the non-rotating case. Using the measured spin period of approximately 100 seconds, we show that the observed pulsations from GW Lib are consistent with the three lowest azimuthal order rotationally modified modes that have the highest frequency in the stellar frame. The high frequencies are needed for the convective driving, but are then apparently shifted to lower frequencies by a combination of their pattern motion and the WD rotation.

preprint2016arXiv

Thermal Runaway During the Evolution of ONeMg Cores towards Accretion-Induced Collapse

We study the evolution of degenerate electron cores primarily composed of the carbon burning products oxygen, neon, and magnesium (hereafter ONeMg cores) that are undergoing compression. Electron capture reactions on A=20 and A=24 isotopes reduce the electron fraction and heat the core. We develop and use a new capability of the Modules for Experiments in Stellar Astrophysics (MESA) stellar evolution code that provides a highly accurate implementation of these key reactions. These new accurate rates and the ability of MESA to perform extremely small spatial zoning demonstrates a thermal runaway in the core triggered by the temperature and density sensitivity of the Ne-20 electron capture reactions. Both analytics and numerics show that this thermal runaway does not trigger core convection, but rather leads to a centrally concentrated (r < km) thermal runaway that will subsequently launch an oxygen deflagration wave from the center of the star. We use MESA to perform a parameter study that quantifies the influence of the magnesium mass fraction, the central temperature, the compression rate, and uncertainties in the electron capture reaction rates on the ONeMg core evolution. This allows us to establish a lower limit on the central density at which the oxygen deflagration wave initiates of $ρ_c > 8.5 \times 10^9\, \textrm{g cm}^{-3}$. Based on previous work and order-of-magnitude calculations, we expect objects which ignite oxygen at or above these densities to collapse and form a neutron star. Calculations such as these are an important step in producing more realistic progenitor models for studies of the signature of accretion-induced collapse.

preprint2016arXiv

Turbulent Chemical Diffusion in Convectively Bounded Carbon Flames

It has been proposed that mixing induced by convective overshoot can disrupt the inward propagation of carbon deflagrations in super-asymptotic giant branch stars. To test this theory, we study an idealized model of convectively bounded carbon flames with 3D hydrodynamic simulations of the Boussinesq equations using the pseudospectral code Dedalus. Because the flame propagation timescale is much longer than the convection timescale, we approximate the flame as fixed in space, and only consider its effects on the buoyancy of the fluid. By evolving a passive scalar field, we derive a {\it turbulent} chemical diffusivity produced by the convection as a function of height, $D_{\rm t}(z)$. Convection can stall a flame if the chemical mixing timescale, set by the turbulent chemical diffusivity, $D_{\rm t}$, is shorter than the flame propagation timescale, set by the thermal diffusivity, $κ$, i.e., when $D_{\rm t}>κ$. However, we find $D_{\rm t}<κ$ for most of the flame because convective plumes are not dense enough to penetrate into the flame. Extrapolating to realistic stellar conditions, this implies that convective mixing cannot stall a carbon flame and that "hybrid carbon-oxygen-neon" white dwarfs are not a typical product of stellar evolution.

preprint2015arXiv

AM Canum Venaticorum Progenitors with Helium Star Donors and the resultant Explosions

We explore the outcome of mass transfer via Roche lobe overflow (RLOF) of $M_{\rm He}\lesssim0.51 M_\odot$ pure helium burning stars in close binaries with white dwarfs (WDs). The evolution is driven by the loss of angular momentum through gravitational wave radiation (GWR), and both stars are modeled using Modules for Experiments in Stellar Astrophysics (MESA). The donors have masses of $M_{\rm He}=0.35, 0.4, \&\ 0.51M_\odot$ and accrete onto WDs of mass $M_{\rm WD}$ from $0.6M_\odot$ to $1.26M_\odot$. The initial orbital periods ($P_{\rm{orb}}$) span 20 to 80 minutes. For all cases, the accretion rate onto the WD is below the stable helium burning range, leading to accumulation of helium followed by unstable ignition. The mass of the convective core in the donors is small enough so that the WD accretes enough helium-rich matter to undergo a thermonuclear runaway in the helium shell before any carbon-oxygen enriched matter is transferred. The mass of the accumulated helium shell depends on $M_{\rm WD}$ and the accretion rate. We show that for $M_{\rm He}\gtrsim0.4 M_\odot$ and $M_{\rm WD}\gtrsim0.8 M_\odot$, the first flash is likely vigorous enough to trigger a detonation in the helium layer. These thermonuclear runaways may be observed as either faint and fast .Ia SNe, or, if the carbon in the core is also detonated, Type Ia SNe. Those that survive the first flash and eject mass will have a temporary increase in orbital separation, but GWR drives the donor back into contact, resuming mass transfer and triggering several subsequent weaker flashes.

preprint2015arXiv

Asteroseismology Can Reveal Strong Internal Magnetic Fields in Red Giant Stars

Internal stellar magnetic fields are inaccessible to direct observations and little is known about their amplitude, geometry and evolution. We demonstrate that strong magnetic fields in the cores of red giant stars can be identified with asteroseismology. The fields can manifest themselves via depressed dipole stellar oscillation modes, which arises from a magnetic greenhouse effect that scatters and traps oscillation mode energy within the core of the star. The Kepler satellite has observed a few dozen red giants with depressed dipole modes which we interpret as stars with strongly magnetized cores. We find field strengths larger than $\sim\! 10^5 \,{\rm G}$ may produce the observed depression, and in one case we infer a minimum core field strength of $\approx \! \! 10^7 \,{\rm G}$.

preprint2015arXiv

Local Radiation Hydrodynamic Simulations of Massive Star Envelopes at the Iron Opacity Peak

We perform three dimensional radiation hydrodynamic simulations of the structure and dynamics of radiation dominated envelopes of massive stars at the location of the iron opacity peak. One dimensional hydrostatic calculations predict an unstable density inversion at this location, whereas our simulations reveal a complex interplay of convective and radiative transport whose behavior depends on the ratio of the photon diffusion time to the dynamical time. The latter is set by the ratio of the optical depth per pressure scale height, $τ_0$, to $τ_c=c/c_g$, where $c_g \approx$ 50 km/s is the isothermal sound speed in the gas alone. When $τ_0 \gg τ_c$, convection reduces the radiation acceleration and removes the density inversion. The turbulent energy transport in the simulations agrees with mixing length theory and provides its first numerical calibration in the radiation dominated regime. When $τ_0 \ll τ_c$, convection becomes inefficient and the turbulent energy transport is negligible. The turbulent velocities exceed $c_g$, driving shocks and large density fluctuations that allow photons to preferentially diffuse out through low-density regions. However, the effective radiation acceleration is still larger than the gravitational acceleration so that the time average density profile contains a modest density inversion. In addition, the simulated envelope undergoes large-scale oscillations with periods of a few hours. The turbulent velocity field may affect the broadening of spectral lines and therefore stellar rotation measurements in massive stars, while the time variable outer atmosphere could lead to variations in their mass loss and stellar radius.

preprint2015arXiv

Magnetar Driven Shock Breakout and Double Peaked Supernova Light Curves

The light curves of some luminous supernovae are suspected to be powered by the spindown energy of a rapidly rotating magnetar. Here we describe a possible signature of the central engine: a burst of shock breakout emission occurring several days after the supernova explosion. The energy input from the magnetar inflates a high-pressure bubble that drives a shock through the pre-exploded supernova ejecta. If the magnetar is powerful enough, that shock will near the ejecta surface and become radiative. At the time of shock breakout, the ejecta will have expanded to a large radius (~10^{14} cm) so that the radiation released is at optical/ultraviolet wavelengths (T ~ 20,000 K) and lasts for several days. The luminosity and timescale of this magnetar driven shock breakout are similar to the first peak observed recently in the double-peaked light curve of SNLSQ14BDQ. However, for a large region of model parameter space, the breakout emission is predicted to be dimmer than the diffusive luminosity from direct magnetar heating. A distinct double peaked light curve may therefore only be conspicuous if thermal heating from the magnetar is suppressed at early times. We describe how such a delay in heating may naturally result from inefficient dissipation and thermalization of the pulsar wind magnetic energy. Without such suppression, the breakout may only be noticeable as a small bump or kink in the early luminosity or color evolution, or as a small but abrupt rise in the photospheric velocity. A similar breakout signature may accompany other central engines in supernovae, such as a black hole accreting fallback material.

preprint2015arXiv

Photoionization Heating of Nova Ejecta by the Post-Outburst Supersoft Source

The expanding ejecta from a classical nova remains hot enough ($\sim10^{4}\, {\rm K}$) to be detected in thermal radio emission for up to years after the cessation of mass loss triggered by a thermonuclear instability on the underlying white dwarf (WD). Nebular spectroscopy of nova remnants confirms the hot temperatures observed in radio observations. During this same period, the unstable thermonuclear burning transitions to a prolonged period of stable burning of the remnant hydrogen-rich envelope, causing the WD to become, temporarily, a super-soft X-ray source. We show that photoionization heating of the expanding ejecta by the hot WD maintains the observed nearly constant temperature of $(1-4)\times10^4\mathrm{~K}$ for up to a year before an eventual decline in temperature due to either the cessation of the supersoft phase or the onset of a predominantly adiabatic expansion. We simulate the expanding ejecta using a one-zone model as well as the Cloudy spectral synthesis code, both incorporating the time-dependent WD effective temperatures for a range of masses from $0.60\ M_{\odot}$ to $1.10\ M_{\odot}$. We show that the duration of the nearly isothermal phase depends most strongly on the velocity and mass of the ejecta and that the ejecta temperature depends on the WD's effective temperature, and hence its mass.

preprint2015arXiv

Population of post-nova supersoft X-ray sources

Novae undergo a supersoft X-ray phase of varying duration after the optical outburst. Such transient post-nova supersoft X-ray sources (SSSs) are the majority of the observed SSSs in M31. In this paper, we use the post-nova evolutionary models of Wolf et al. to compute the expected population of post-nova SSSs in M31. We predict that depending on the assumptions about the WD mass distribution in novae, at any instant there are about 250-600 post-nova SSSs in M31 with (unabsorbed) 0.2-1.0 keV luminosity L_x>10^36 erg/s. Their combined unabsorbed luminosity is of the order of ~10^39 erg/s. Their luminosity distribution shows significant steepening around log(L_x)~37.7-38 and becomes zero at L_x~2x10^38 erg/s, the maximum L_x achieved in the post-nova evolutionary tracks. Their effective temperature distribution has a roughly power law shape with differential slope of ~4-6 up to the maximum temperature of T_eff~1.5x10^6 K. We compare our predictions with the results of the XMM-Newton monitoring of the central field of M31 between 2006 and 2009. The predicted number of post-nova SSSs exceed the observed number by a factor of ~2-5, depending on the assumed WD mass distribution in novae. This is good agreement, considering the number and magnitude of uncertainties involved in calculations of the post-nova evolutionary models and their X-ray output. Furthermore, only a moderate circumstellar absorption, with hydrogen column density of the order of ~10^21 cm^-2, will remove the discrepancy.

preprint2014arXiv

An Accreting White Dwarf near the Chandrasekhar Limit in the Andromeda Galaxy

The iPTF detection of the most recent outburst of the recurrent nova system RX J0045.4+4154 in the Andromeda Galaxy has enabled the unprecedented study of a massive ($M>1.3\ M_\odot$) accreting white dwarf (WD). We detected this nova as part of the near daily iPTF monitoring of M31 to a depth of $R\approx 21$\,mag and triggered optical photometry, spectroscopy and soft X-ray monitoring of the outburst. Peaking at an absolute magnitude of $M_R=-6.6$ mag, and with a decay time of 1 mag per day, it is a faint and very fast nova. It shows optical emission lines of He/N and expansion velocities of 1900 to 2600 km s$^{-1}$ 1--4 days after the optical peak. The {\it Swift} monitoring of the X-ray evolution revealed a supersoft source (SSS) with $kT_{\rm eff}\approx 90-110\ {\rm eV}$ that appeared within 5 days after the optical peak, and lasted only 12 days. Most remarkably, this is not the first event from this system, rather it is a recurrent nova with a time between outbursts of approximately 1 year, the shortest known. Recurrent X-ray emission from this binary was detected by ROSAT in 1992 and 1993, and the source was well characterized as a $M>1.3\ M_\odot$ WD SSS. Based on the observed recurrence time between different outbursts, the duration and effective temperature of the SS phase, MESA models of accreting WDs allow us to constrain the accretion rate to $\dot M>1.7\times10^{-7}\ {M_{\odot}\ {\rm yr}}^{-1}$ and WD mass $>1.30\ M_{\odot}$. If the WD keeps $30\%$ of the accreted material, it will take less than a Myr to reach core densities high enough for carbon ignition (if made of C/O) or electron capture (if made of O/Ne) to end the binary evolution.

preprint2014arXiv

Angular momentum transport within evolved low-mass stars

Asteroseismology of 1.0-2.0 Msun red giants by the Kepler satellite has enabled the first definitive measurements of interior rotation in both first ascent red giant branch (RGB) stars and those on the Helium burning clump. The inferred rotation rates are 10-30 days for the ~0.2Msun He degenerate cores on the RGB and 30-100 days for the He burning core in a clump star. Using the MESA code we calculate state-of-the-art stellar evolution models of low mass rotating stars from the zero-age main sequence to the cooling white dwarf (WD) stage. We include transport of angular momentum due to rotationally induced instabilities and circulations, as well as magnetic fields in radiative zones (generated by the Tayler-Spruit dynamo). We find that all models fail to predict core rotation as slow as observed on the RGB and during core He burning, implying that an unmodeled angular momentum transport process must be operating on the early RGB of low mass stars. Later evolution of the star from the He burning clump to the cooling WD phase appears to be at nearly constant core angular momentum. We also incorporate the adiabatic pulsation code, ADIPLS, to explicitly highlight this shortfall when applied to a specific Kepler asteroseismic target, KIC8366239. The MESA inlist adopted to calculate the models in this paper can be found at \url{https://authorea.com/1608/} (bottom of the document).

preprint2014arXiv

Identification of the Optical Counterpart of Fermi Black Widow Millisecond Pulsar PSR J1544+4937

We report the optical identification of the companion to the {\it Fermi} black widow millisecond pulsar PSR J1544+4937. We find a highly variable source on Keck LRIS images at the nominal pulsar position, with 2 magnitude variations over orbital period in the B, g, R, and I bands. The nearly achromatic light curves are difficult to explain with a simply irradiated hemisphere model, and suggest that the optical emission is dominated by a nearly isothermal hot patch on the surface of the companion facing the pulsar. We roughly constrain the distance to PSR J1544+4937 to be between 2 and 5 kpc. A more reliable distance measurement is needed in order to constrain the composition of the companion.

preprint2014arXiv

Possible Detection of the Stellar Donor or Remnant for the Type Iax Supernova 2008ha

Type Iax supernovae (SNe Iax) are thermonuclear explosions that are related to SNe Ia, but are physically distinct. The most important differences are that SNe Iax have significantly lower luminosity (1% - 50% that of typical SNe Ia), lower ejecta mass (~0.1 - 0.5 M_sun), and may leave a bound remnant. The most extreme SN Iax is SN 2008ha, which peaked at M_V = -14.2 mag, about 5 mag below that of typical SNe Ia. Here, we present Hubble Space Telescope (HST) images of UGC 12682, the host galaxy of SN 2008ha, taken 4.1 years after the peak brightness of SN 2008ha. In these deep, high-resolution images, we detect a source coincident (0.86 HST pixels; 0.043"; 1.1 sigma) with the position of SN 2008ha with M_F814W = -5.4 mag. We determine that this source is unlikely to be a chance coincidence, but that scenario cannot be completely ruled out. If this source is directly related to SN 2008ha, it is either the luminous bound remnant of the progenitor white dwarf or its companion star. The source is consistent with being an evolved >3 M_sun initial mass star, and is significantly redder than the SN Iax 2012Z progenitor system, the first detected progenitor system for a thermonuclear SN. If this source is the companion star for SN 2008ha, there is a diversity in SN Iax progenitor systems, perhaps related to the diversity in SN Iax explosions. If the source is the bound remnant of the white dwarf, it must have expanded significantly. Regardless of the nature of this source, we constrain the progenitor system of SN 2008ha to have an age of <80 Myr.

preprint2014arXiv

The Ignition of Carbon Detonations via Converging Shock Waves in White Dwarfs

The progenitor channel responsible for the majority of Type Ia supernovae is still uncertain. One emergent scenario involves the detonation of a He-rich layer surrounding a C/O white dwarf, which sends a shock wave into the core. The quasi-spherical shock wave converges and strengthens at an off-center location, forming a second, C-burning, detonation that disrupts the whole star. In this paper, we examine this second detonation of the double detonation scenario using a combination of analytic and numeric techniques. We perform a spatially resolved study of the imploding shock wave and outgoing detonation and calculate the critical imploding shock strengths needed to achieve a core C detonation. We find that He detonations in recent two-dimensional simulations yield converging shock waves that are strong enough to ignite C detonations in high-mass C/O cores, with the caveat that a truly robust answer requires multi-dimensional detonation initiation calculations. We also find that convergence-driven detonations in low-mass C/O cores and in O/Ne cores are harder to achieve and are perhaps unrealized in standard binary evolution.

preprint2014arXiv

X-ray diagnostics of chemical composition of the accretion disk and donor star in UCXBs II: XMM-Newton observations

We search for the Fe K$α$ line in spectra of Ultra Compact X-ray Binaries (UCXBs). For this purpose we have analyzed XMM-Newton observations of five confirmed UCXBs. We find that the object 2S 0918-549 - whose optical spectrum bears tentative signatures of a C/O accretion disk - is devoid of any emission features in the 6-7 keV range, with an upper limit of less than 10 eV for the equivalent width (EW) of the iron line. 4U 1916-05 - whose optical spectrum is consistent with reflection from a He-rich accretion disk - exhibits a bright broad iron emission line. This behavior is in agreement with the theoretical predictions presented in Koliopanos, Gilfanov and Bildsten (2013). Namely, we expect strong suppression of the Fe K$α$ emission line in spectra originating in moderately bright (LogLx less than $\approx$ 37.5) UCXBs with C/O or O/Ne/Mg-rich donors. On the other hand the EW of the iron line in spectra from UCXBs with He-rich donors is expected to retain its nominal value of $\approx$ 100 eV. Our analysis also reveals a strong Fe K$α$ line in the spectrum of 4U 0614+091. This detection points towards a He-rich donor and seems to be at odds with the source's classification as C/O-rich. Nevertheless, a He-rich donor would explain the bursting activity reported for this system. Lastly, based on our theoretical predictions, we attribute the lack of a strong iron emission line - in the two remaining UCXB sources in our sample (XTE J1807-294, 4U 0513-40) - as an indication of a C/O or O/Ne/Mg white dwarf donor. From the upper limits of the Fe K$α$ line EW in 4U 0513-40, 2S 0918-549 and XTE J1807-294 we obtain a lower limit on the oxygen-to-iron ratio,O/Fe$\ge$ 10$\times$[O/Fe]$_{\odot}$.

preprint2013arXiv

Direct evidence for a supernova interacting with a large amount of hydrogen-free circumstellar material

We present our observations of SN 2010mb, a Type Ic SN lacking spectroscopic signatures of H and He. SN 2010mb has a slowly-declining light curve ($\sim600\,$days) that cannot be powered by $^{56}$Ni/$^{56}$Co radioactivity, the common energy source for Type Ic SNe. We detect signatures of interaction with hydrogen-free CSM including a blue quasi-continuum and, uniquely, narrow oxygen emission lines that require high densities ($\sim10^9$cm$^{-3}$). From the observed spectra and light curve we estimate that the amount of material involved in the interaction was $\sim3$M$_{\odot}$. Our observations are in agreement with models of pulsational pair-instability SNe described in the literature.

preprint2013arXiv

Hydrogen Burning on Accreting White Dwarfs: Stability, Recurrent Novae, and the Post-Novae Supersoft Source

We examine the properties of white dwarfs (WDs) accreting hydrogen-rich matter in and near the stable burning regime of accretion rates as modeled by time-dependent calculations done with Modules for Experiments in Stellar Astrophysics (MESA). We report the stability boundary for WDs of masses between 0.51 solar masses and 1.34 solar masses as found via time-dependent calculations. We also examine recurrent novae that are accreting at rates close to, but below, the stable burning limit and report their recurrence times and ignition masses. Our dense grid in accretion rates finds the expected minimum possible recurrence times as a function of the WD mass. This enables inferences to be made about the minimum WD mass possible to reach a specific recurrence time. We compare our computational models of post-outburst novae to the stably burning WDs and explicitly calculate the duration and effective temperature (Teff) of the post-novae WD in the supersoft phase. We agree with the measured turnoff time - Teff relation in M31 by Henze and collaborators, infer WD masses in the 1.0-1.3 solar masses range, and predict ejection masses consistent with those observed. We close by commenting on the importance of the hot helium layer generated by stable or unstable hydrogen burning for the short- and long-term evolution of accreting white dwarfs.

preprint2013arXiv

Modules for Experiments in Stellar Astrophysics (MESA): Giant Planets, Oscillations, Rotation, and Massive Stars

We substantially update the capabilities of the open source software package Modules for Experiments in Stellar Astrophysics (MESA), and its one-dimensional stellar evolution module, MESA Star. Improvements in MESA Star's ability to model the evolution of giant planets now extends its applicability down to masses as low as one-tenth that of Jupiter. The dramatic improvement in asteroseismology enabled by the space-based Kepler and CoRoT missions motivates our full coupling of the ADIPLS adiabatic pulsation code with MESA Star. This also motivates a numerical recasting of the Ledoux criterion that is more easily implemented when many nuclei are present at non-negligible abundances. This impacts the way in which MESA Star calculates semi-convective and thermohaline mixing. We exhibit the evolution of 3-8 Msun stars through the end of core He burning, the onset of He thermal pulses, and arrival on the white dwarf cooling sequence. We implement diffusion of angular momentum and chemical abundances that enable calculations of rotating-star models, which we compare thoroughly with earlier work. We introduce a new treatment of radiation-dominated envelopes that allows the uninterrupted evolution of massive stars to core collapse. This enables the generation of new sets of supernovae, long gamma-ray burst, and pair-instability progenitor models. We substantially modify the way in which MESA Star solves the fully coupled stellar structure and composition equations, and we show how this has improved MESA's performance scaling on multi-core processors. Updates to the modules for equation of state, opacity, nuclear reaction rates, and atmospheric boundary conditions are also provided. We describe the MESA Software Development Kit (SDK) that packages all the required components needed to form a unified and maintained build environment for MESA. [Abridged]

preprint2013arXiv

Properties of the Eclipsing Double-White Dwarf Binary NLTT 11748

We present high-quality ULTRACAM photometry of the eclipsing detached double-white dwarf binary NLTT 11748. This system consists of a carbon/oxygen white dwarf and an extremely-low mass (< 0.2 Msun) helium-core white dwarf in a 5.6 hr orbit. To date such extremely-low mass WDs, which can have thin, stably-burning outer layers, have been modeled via poorly-constrained atmosphere and cooling calculations where uncertainties in the detailed structure can strongly influence the eventual fates of these systems when mass-transfer begins. With precise (individual precision ~1%) high-cadence (~2 s) multi-color photometry of multiple primary and secondary eclipses spanning >1.5 yr, we constrain the masses and radii of both objects in the NLTT 11748 system to a statistical uncertainty of a few percent. However, we find that overall uncertainty in the thickness of the envelope of the secondary carbon/oxygen white dwarf leads to a larger (~13%) systematic uncertainty in the primary He WD's mass. Over the full range of possible envelope thicknesses we find that our primary mass (0.136-0.162 Msun) and surface gravity (log(g)=6.32-6.38; radii are 0.0423-0.0433 Rsun) constraints do not agree with previous spectroscopic determinations. We use precise eclipse timing to detect the Romer delay at 7 sigma significance, providing an additional weak constraint on the masses and limiting the eccentricity to e*cos(omega)= -4e-5 +/- 5e-5. Finally, we use multi-color data to constrain the secondary's effective temperature (7600+/-120 K) and cooling age (1.6-1.7 Gyr).

preprint2013arXiv

R Coronae Borealis Stars in M31 from the Palomar Transient Factory

We report the discovery of R Coronae Borealis (RCB) stars in the Andromeda galaxy (M31) using the Palomar Transient Factory (PTF). RCB stars are rare hydrogen-deficient, carbon-rich supergiant variables, most likely the merger products of two white dwarfs. These new RCBs, including two confirmed ones and two candidates, are the first to be found beyond the Milky Way and the Magellanic Clouds. All of M31 RCBs showed >1.5 mag irregular declines over timescales of weeks to months. Due to the limiting magnitude of our data (R~21-22 mag), these RCB stars have R~19.5 to 20.5 mag at maximum light, corresponding to M_R= -4 to -5, making them some of the most luminous RCBs known. Spectra of two objects show that they are warm RCBs, similar to the Milky Way RCBs RY Sgr and V854 Cen. We consider these results, derived from a pilot study of M31 variables, as an important proof-of-concept for the study of rare bright variables in nearby galaxies with the PTF or other synoptic surveys.

preprint2013arXiv

The Effects of Curvature and Expansion on Helium Detonations on White Dwarf Surfaces

Accreted helium layers on white dwarfs have been highlighted for many decades as a possible site for a detonation triggered by a thermonuclear runaway. In this paper, we find the minimum helium layer thickness that will sustain a steady laterally propagating detonation and show that it depends on the density and composition of the helium layer, specifically C12 and O16. Detonations in these thin helium layers have speeds slower than the Chapman-Jouget (CJ) speed from complete helium burning, v_cj=1.5e9 cm/s. Though gravitationally unbound, the ashes still have unburned helium (~ 80% in the thinnest cases) and only reach up to heavy elements such as Ca40, Ti44, Cr48, and Fe52. It is rare for these thin shells to generate large amounts of Ni56. We also find a new set of solutions that can propagate in even thinner helium layers when O16 is present at a minimum mass fraction of ~0.07. Driven by energy release from alpha captures on O16 and subsequent elements, these slow detonations only create ashes up to Si28 in the outer detonated He shell. We close by discussing how the unbound helium burning ashes may create faint and fast ".Ia" supernovae as well as events with virtually no radioactivity, and speculate on how the slower helium detonation velocities impact the off-center ignition of a carbon detonation that could cause a Type Ia supernova in the double detonation scenario.

preprint2013arXiv

X-ray diagnostics of chemical composition of the accretion disk and donor star in ultra-compact X-ray binaries

Non-solar composition of the donor star in ultra-compact X-ray binaries may have a pronounced effect on the fluorescent lines appearing in their spectra due to reprocessing of primary radiation by the accretion disk and the white dwarf surface. We show that the most dramatic and easily observable consequence of the anomalous C/O abundance, is the significant, by more than an order of magnitude, attenuation of the Ka line of iron. It is caused by screening of the presence of iron by oxygen - in the C/O dominated material the main interaction process for a E ~ 7keV photon is absorption by oxygen rather than by iron, contrary to the solar composition case. Ionization of oxygen at high mass accretion rates adds a luminosity dependence to this behavior - the iron line is significantly suppressed only at low luminosity, log(LX) less than 37-37.5, and should recover its nominal strength at higher luminosity. The increase of the EW of the Ka lines of carbon and oxygen, on the other hand, saturates at rather moderate values. Screening by He is less important, due to its low ionization threshold and because in the accretion disk it is mostly ionized. Consequently, in the case of the He-rich donor, the iron line strength remains close to its nominal value, determined by the iron abundance in the accretion disk. This opens the possibility of constraining the nature of donor stars in UCXBs by means of X-ray spectroscopy with moderate energy resolution.

preprint2012arXiv

Circumstellar Shell Formation in Symbiotic Recurrent Novae

We present models of spherically symmetric recurrent nova shells interacting with circumstellar material in a symbiotic system composed of a red giant expelling a wind, and a white dwarf accreting from this material. Recurrent nova eruptions periodically eject material at high velocities ($\gtrsim 10^3$ km/s) into the red giant wind profile, creating a decelerating shock wave as circumstellar material is swept up. High circumstellar material densities cause the shocked wind and ejecta to have very short cooling times of days to weeks. Thus, the late time evolution of the shell is determined by momentum conservation instead of energy conservation. We compute and show evolutionary tracks of shell deceleration, as well as post-shock structure. After sweeping up all the red giant wind, the shell coasts at a velocity $\sim 100$ km/s, depending on system parameters. These velocities are similar to those measured in blue-shifted circumstellar material from the symbiotic nova RS Oph, as well as a few Type Ia supernovae that show evidence of circumstellar material, such as 2006X, 2007le, and PTF 11kx. Supernovae occurring in such systems may not show circumstellar material interaction until the inner nova shell gets hit by the supernova ejecta, days to months after the explosion.

preprint2012arXiv

Laterally Propagating Detonations in Thin Helium Layers on Accreting White Dwarfs

Theoretical work has shown that intermediate mass (0.01Msun<M_He<0.1Msun) Helium shells will unstably ignite on the accreting white dwarf (WD) in an AM CVn binary. For more massive (M>0.8Msun) WDs, these helium shells can be dense enough (5x10^5 g/cc) that the convectively burning region runs away on a timescale comparable to the sound travel time across the shell; raising the possibility for an explosive outcome. The nature of the explosion (i.e. deflagration or detonation) remains ambiguous. In the case of detonation, this causes a laterally propagating front whose properties in these geometrically thin and low density shells we begin to study here. Our calculations show that the radial expansion time of <0.1 s leads to incomplete helium burning, in agreement with recent work by Sim and collaborators, but that the nuclear energy released is still adequate to realize a self-sustaining detonation propagating laterally at slower than the Chapman-Jouguet speed. Our simulations resolve the subsonic region behind the front and are consistent with a direct computation of the reaction structure from the shock strength. The ashes are typically He rich, and consist of predominantly Ti-44, Cr-48, along with a small amount of Fe-52, with very little Ni-56 and with significant Ca-40 in carbon-enriched layers. If this helium detonation results in a Type Ia Supernova, its spectral signatures would appear for the first few days after explosion. (abridged)

preprint2012arXiv

MESA Models of Classical Nova Outbursts: The Multicycle Evolution and Effects of Convective Boundary Mixing

Novae are cataclysmic variables driven by accretion of H-rich material onto a white-dwarf (WD) star from its low-mass main-sequence binary companion. New time-domain observational capabilities, such as the Palomar Transient Factory and Pan-STARRS, have revealed a diversity of their behaviour that should be theoretically addressed. Nova outbursts depend sensitively on nuclear physics data, and more readily available nova simulations are needed in order to effectively prioritize experimental effort in nuclear astrophysics. In this paper we use the MESA stellar evolution code to construct multicycle nova evolution sequences with CO WD cores. We explore a range of WD masses and accretion rates as well as the effect of different cooling times before the onset of accretion. In addition, we study the dependence on the elemental abundance distribution of accreted material and convective boundary mixing at the core-envelope interface. Models with such convective boundary mixing display an enrichment of the accreted envelope with C and O from the underlying white dwarf that is commensurate with observations. We compare our results with the previous work and investigate a new scenario for novae with the 3He-triggered convection.

preprint2012arXiv

Orbital Evolution of Compact White Dwarf Binaries

The new-found prevalence of extremely low mass (ELM, Mhe<0.2 Msun) helium white dwarfs (WDs) in tight binaries with more massive WDs has raised our interest in understanding the nature of their mass transfer. Possessing small (Menv~1e-3 Msun) but thick hydrogen envelopes, these objects have larger radii than cold WDs and so initiate mass transfer of H-rich material at orbital periods of 6-10 minutes. Building on the original work of D'Antona et al., we confirm the 1e6 yr period of continued inspiral with mass transfer of H-rich matter and highlight that the inspiraling direct-impact double WD binary HM Cancri likely has an ELM WD donor. The ELM WDs have less of a radius expansion under mass loss, thus enabling a larger range of donor masses that can stably transfer matter and become a He mass transferring AM CVn binary. Even once in the long-lived AM CVn mass transferring stage, these He WDs have larger radii due to their higher entropy from the prolonged H burning stage.

preprint2012arXiv

The Long-Term Evolution of Double White Dwarf Mergers

In this paper, we present a model for the long-term evolution of the merger of two unequal mass C/O white dwarfs (WDs). After the dynamical phase of the merger, magnetic stresses rapidly redistribute angular momentum, leading to nearly solid-body rotation on a viscous timescale of 1e-4 to 1 yr, long before significant cooling can occur. Due to heating during the dynamical and viscous phases, the less massive WD is transformed into a hot, slowly rotating, and radially extended envelope supported by thermal pressure. Following the viscous phase of evolution, the maximum temperature near the envelope base may already be high enough to begin off-center convective carbon-burning. If not, Kelvin-Helmholtz contraction of the inner region of the envelope on a thermal timescale of 1e3-1e4 yr compresses the base of the envelope, again yielding off-center burning. As a result, the long-term evolution of the merger remnant is similar to that seen in previous calculations: the burning shell diffuses inwards over ~1e4 yr, eventually yielding a high-mass O/Ne WD or a collapse to a neutron star. During the cooling and shell-burning phases, the merger remnant radiates near the Eddington limit. Given the double WD merger rate of a few per 1000 yr, tens of these ~1e38 erg/s sources should exist in a Milky Way-type galaxy. While the end result is similar to that of previous studies, the physical picture and the dynamical state of the matter in our model differ from previous work. Furthermore, remaining uncertainties related to the convective structure near the photosphere and mass loss during the thermal evolution may significantly affect our conclusions. Thus, future work within the context of the physical model presented here is required to better address the eventual fate of double WD mergers, including those for which one or both of the components is a He WD.

preprint2012arXiv

The Observational Signatures of Convectively Excited Internal Gravity Modes in Main Sequence Stars

We predict the flux and surface velocity perturbations produced by convectively excited gravity modes (g-modes) in main sequence stars. Core convection in massive stars can excite g-modes to sufficient amplitudes to be detectable with high precision photometry by Kepler and CoRoT, if the thickness of the convective overshoot region is < 30 per cent of a pressure scale height. The g-modes manifest as excess photometric variability, with amplitudes of ~ 10 micromagnitudes at frequencies ~ 10 microHz (0.8 1/d) near the solar metallicity zero-age main sequence. The flux variations are largest for stars with M > 5 solar masses, but are potentially detectable down to M ~ 2 - 3 solar masses. During the main sequence evolution, radiative damping decreases such that ever lower frequency modes reach the stellar surface and flux perturbations reach up to ~ 100 micromagnitudes at the terminal-age main sequence. Using the same convective excitation model, we confirm previous predictions that solar g-modes produce surface velocity perturbations of < 0.3 mm/s. This implies that stochastically excited g-modes are more easily detectable in the photometry of massive main sequence stars than in the Sun.

preprint2011arXiv

A Search for Pulsations in Helium White Dwarfs

The recent plethora of sky surveys, especially the Sloan Digital Sky Survey, have discovered many low-mass (M < 0.45 Msun) white dwarfs that should have cores made of nearly pure helium. These WDs come in two varieties; those with masses 0.2 < M < 0.45 Msun and H envelopes so thin that they rapidly cool, and those with M < 0.2 Msun (often called extremely low mass, ELM, WDs) that have thick enough H envelopes to sustain 10^9 years of H burning. In both cases, these WDs evolve through the ZZ Ceti instability strip, Teff = 9,000-12,000 K, where g-mode pulsations always occur in Carbon/Oxygen WDs. This expectation, plus theoretical work on the contrasts between C/O and He core WDs, motivated our search for pulsations in 12 well characterized helium WDs. We report here on our failure to find any pulsators amongst our sample. Though we have varying amplitude limits, it appears likely that the theoretical expectations regarding the onset of pulsations in these objects requires closer consideration. We close by encouraging additional observations as new He WD samples become available, and speculate on where theoretical work may be needed.

preprint2011arXiv

Acoustic Signatures of the Helium Core Flash

All evolved stars with masses M <2 solar masses undergo an initiating off-center helium core flash in their 0.48 solar mass He core as they ascend the red giant branch (RGB). This off-center flash is the first of a few successive helium shell subflashes that remove the core electron degeneracy over 2 Myrs, converting the object into a He burning star. Though characterized by Thomas over 40 years ago, this core flash phase has yet to be observationally probed. Using the Modules for Experiments in Stellar Astrophysics (MESA) code, we show that red giant asteroseismology enabled by space-based photometry (i.e. Kepler and CoRoT) can probe these stars during the flash. The rapid (< 100,000 years) contraction of the red giant envelope after the initiating flash dramatically improves the coupling of the p-modes to the core g-modes, making the detection of l=1 mixed modes possible for these 2 Myrs. This duration implies that 1 in 35 stars near the red clump in the HR diagram will be in their core flash phase. During this time, the star has a g-mode period spacing of Delta P_g = 70-100 seconds, lower than the Delta P_g=250 seconds of He burning stars in the red clump, but higher than the RGB stars at the same luminosity. This places them in an underpopulated part of the large frequency spacing (Delta nu) vs. Delta P_g diagram that should ease their identification amongst the thousands of observed red giants.

preprint2011arXiv

Calcium-rich gap transients in the remote outskirts of galaxies

From the first two seasons of the Palomar Transient Factory, we identify three peculiar transients (PTF09dav, PTF10iuv, PTF11bij) with five distinguishing characteristics: peak luminosity in the gap between novae and supernovae (M_R = 15.5 to -16.5), rapid photometric evolution (rise-time ~12--15 days), large photospheric velocities (~6000 to 11000 km/s), early spectroscopic evolution into nebular phase (~1 to 3 months) and peculiar nebular spectra dominated by Calcium. We also culled the extensive decade-long Lick Observatory Supernova Search database and identified an additional member of this group, SN 2007ke. Our choice of photometric and spectroscopic properties was motivated by SN 2005E (Perets et al. 2010). To our surprise, as in the case of SN 2005E, all four members of this group are also clearly offset from the bulk of their host galaxy. Given the well-sampled early and late-time light curves, we derive ejecta masses in the range of 0.4--0.7 Msun. Spectroscopically, we find that there may be a diversity in the photospheric phase, but the commonality is in the unusual nebular spectra. Our extensive follow-up observations rule out standard thermonuclear and standard core-collapse explosions for this class of "Calcium-rich gap" transients. If the progenitor is a white dwarf, we are likely seeing a detonation of the white dwarf core and perhaps, even shockfront interaction with a previously ejected nova shell. In the less likely scenario of a massive star progenitor, a very non-standard channel specific to a low-metallicity environment needs to be invoked (e.g., ejecta fallback leading to black hole formation). Detection (or lack thereof) of a faint underlying host (dwarf galaxy, cluster) will provide a crucial and decisive diagnostic to choose between these alternatives.

preprint2011arXiv

Constraints on the Progenitor System of the Type Ia Supernova SN 2011fe/PTF11kly

Type Ia supernovae (SNe) serve as a fundamental pillar of modern cosmology, owing to their large luminosity and a well-defined relationship between light-curve shape and peak brightness. The precision distance measurements enabled by SNe Ia first revealed the accelerating expansion of the universe, now widely believed (though hardly understood) to require the presence of a mysterious "dark" energy. General consensus holds that Type Ia SNe result from thermonuclear explosions of a white dwarf (WD) in a binary system; however, little is known of the precise nature of the companion star and the physical properties of the progenitor system. Here we make use of extensive historical imaging obtained at the location of SN 2011fe/PTF11kly, the closest SN Ia discovered in the digital imaging era, to constrain the visible-light luminosity of the progenitor to be 10-100 times fainter than previous limits on other SN Ia progenitors. This directly rules out luminous red giants and the vast majority of helium stars as the mass-donating companion to the exploding white dwarf. Any evolved red companion must have been born with mass less than 3.5 times the mass of the Sun. These observations favour a scenario where the exploding WD of SN 2011fe/PTF11kly, accreted matter either from another WD, or by Roche-lobe overflow from a subgiant or main-sequence companion star.

preprint2011arXiv

PTF10iya: A short-lived, luminous flare from the nuclear region of a star-forming galaxy

We present the discovery and characterisation of PTF10iya, a short-lived (dt ~ 10 d, with an optical decay rate of ~ 0.3 mag per d), luminous (M_g ~ -21 mag) transient source found by the Palomar Transient Factory. The ultraviolet/optical spectral energy distribution is reasonably well fit by a blackbody with T ~ 1-2 x 10^4 K and peak bolometric luminosity L_BB ~ 1-5 x 10^44 erg per s (depending on the details of the extinction correction). A comparable amount of energy is radiated in the X-ray band that appears to result from a distinct physical process. The location of PTF10iya is consistent with the nucleus of a star-forming galaxy (z = 0.22405 +/- 0.00006) to within 350 mas (99.7 per cent confidence radius), or a projected distance of less than 1.2 kpc. At first glance, these properties appear reminiscent of the characteristic "big blue bump" seen in the near-ultraviolet spectra of many active galactic nuclei (AGNs). However, emission-line diagnostics of the host galaxy, along with a historical light curve extending back to 2007, show no evidence for AGN-like activity. We therefore consider whether the tidal disruption of a star by an otherwise quiescent supermassive black hole may account for our observations. Though with limited temporal information, PTF10iya appears broadly consistent with the predictions for the early "super-Eddington" phase of a solar-type star disrupted by a ~ 10^7 M_sun black hole. Regardless of the precise physical origin of the accreting material, the large luminosity and short duration suggest that otherwise quiescent galaxies can transition extremely rapidly to radiate near the Eddington limit; many such outbursts may have been missed by previous surveys lacking sufficient cadence.

preprint2011arXiv

Supernova 2011fe from an Exploding Carbon-Oxygen White Dwarf Star

Type Ia supernovae (SNe Ia) have been used empirically as standardized candles to reveal the accelerating universe even though fundamental details, such as the nature of the progenitor system and how the star explodes, remained a mystery. There is consensus that a white dwarf star explodes after accreting matter in a binary system, but the secondary could be anything from a main sequence star to a red giant, or even another white dwarf. The uncertainty stems from the fact that no recent SN Ia has been discovered close enough to detect the stars before explosion. Here we report early observations of SN 2011fe (PTF11kly) in M101 at a distance of 6.4 Mpc, the closest SN Ia in the past 25 years. We find that the exploding star was likely a carbon-oxygen white dwarf, and from the lack of an early shock we conclude that the companion was most likely a main sequence star. Early spectroscopy shows high-velocity oxygen that varies on a time scale of hours and extensive mixing of newly synthesized intermediate mass elements in the outermost layers of the supernova. A companion paper uses pre-explosion images to rule out luminous red giants and most helium stars as companions.

preprint2010arXiv

A ground-based measurement of the relativistic beaming effect in a detached double WD binary

We report on the first ground-based measurement of the relativistic beaming effect (aka Doppler boosting). We observed the beaming effect in the detached, non-interacting eclipsing double white dwarf (WD) binary NLTT 11748. Our observations were motivated by the system's high mass ratio and low luminosity ratio, leading to a large beaming-induced variability amplitude at the orbital period of 5.6 hr. We observed the system during 3 nights at the 2.0m Faulkes Telescope North with the SDSS-g' filter, and fitted the data simultaneously for the beaming, ellipsoidal and reflection effects. Our fitted relative beaming amplitude is (3.0 +/- 0.4) x 10^(-3), consistent with the expected amplitude from a blackbody spectrum given the photometric primary radial velocity amplitude and effective temperature. This result is a first step in testing the relation between the photometric beaming amplitude and the spectroscopic radial velocity amplitude in NLTT 11748 and similar systems. We did not identify any variability due to the ellipsoidal or reflection effects, consistent with their expected undetectable amplitude for this system. Low-mass, helium-core WDs are expected to reside in binary systems where in some of those systems the binary companion is a faint C/O WD and the two stars are detached and non-interacting, as in the case of NLTT 11748. The beaming effect can be used to search for the faint binary companion in those systems using wide-band photometry.

preprint2010arXiv

Clearing the Gas from Globular Clusters & Dwarf Spheroidals with Classical Novae

Observations of the intra-cluster medium (ICM) in galactic globular clusters (GCs) show a systematic deficiency in ICM mass as compared to that expected from accumulation of stellar winds in the time available between galactic plane crossings. In this paper, we reexamine the original hypothesis of Scott and Durisen that hydrogen-rich explosions on accreting white dwarfs, classical novae (CNe), will sweep out the ICM from the cluster more frequently than galactic plane crossings. From the CNe rate and stellar mass-loss rate, this clearing mechanism predicts that ~ 0.03 M_sun should be present in <= 10^5 M_sun GCs. We model the expanding remnant made from the 10^-4 M_sun nova ejecta and show that it escapes long before it has cooled. We discuss the few positive ICM measurements and use a Monte-Carlo simulation of the accumulation and CNe recurrence times to reveal the possible variance in the ICM masses for the higher mass (> 5x10^5 M_sun) GCs. We find that nova shells are effective at clearing the ICM in low-mass GCs (<= 10^5 M_sun), whereas higher-mass clusters may experience a quiescent time between novae long enough to prevent the next nova shell from escaping. The nova clearing mechanism will also operate in ultra-faint Milky Way satellites, where many upper limits on gas masses are available.

preprint2010arXiv

Diffusive Nuclear Burning of Helium on Neutron Stars

Diffusive nuclear burning of H by an underlying material capable of capturing protons can readily consume H from the surface of neutron stars (NSs) during their early cooling history. In the absence of subsequent accretion, it will be depleted from the photosphere. We now extend diffusive nuclear burning to He, motivated by the recent observation by Ho \& Heinke of a carbon atmosphere on the NS in the Cassiopeia A supernova remnant. We calculate the equilibrium structure of He on an underlying $α$ capturing material, accounting for thermal, mass defect, and Coulomb corrections on the stratification of material with the same zeroth order $μ_e = A/Z$. We show that Coulomb corrections dominate over thermal and mass defect corrections in the highly degenerate part of the envelope. We also show that the bulk of the He sits deep in the envelope rather than near the surface. Thus, even if the photospheric He abundance is low, the total He column could be substantially larger than the photospheric column, which may have implications for rapid surface evolution ($\approx 1$ yr timescales) of neutron stars. When nuclear reactions are taken into account, we find that for base temperatures $\gtrsim 1.6 \times 10^8$ K, He is readily captured onto C. As these high temperatures are present during the early stages of NS evolution, we expect that the primordial He is completely depleted from the NS surface like the case for primordial H. We also find that magnetic fields $\lesssim 10^{12}$ G do not affect our conclusions. Armed with the results of this work and our prior efforts, we expect that primordial H and He are depleted, and so any observed H or He on the surfaces of these NS must be due to subsequent accretion (with or without spallation). If this subsequent accretion can be prevented, the underlying mid-Z material would be exposed.

preprint2010arXiv

Discovery of the Eclipsing Detached Double White Dwarf Binary NLTT 11748

We report the discovery of the first eclipsing detached double white dwarf (WD) binary. In a pulsation search, the low-mass helium core WD NLTT 11748 was targeted for fast (approx 1 minute) differential photometry with the Las Cumbres Observatory's Faulkes Telescope North. Rather than pulsations, we discovered approx 180 s 3%-6% dips in the photometry. Subsequent radial velocity measurements of the primary white dwarf from the Keck telescope found variations with a semi-amplitude K_1 = 271 +/- 3 km/s, and confirmed the dips as eclipses caused by an orbiting WD with a mass M_2 = 0.648-0.771 M_sun for M_1 = 0.1-0.2 M_sun. We detect both the primary and secondary eclipses during the P_orb = 5.64 hr orbit and measure the secondary's brightness to be 3.5% +/- 0.3% of the primary at SDSS-g'. Assuming that the secondary follows the mass-radius relation of a cold C/O WD and including the effects of microlensing in the binary, the primary eclipse yields a primary radius of R_1 = 0.043-0.039 R_sun for M_1 = 0.1-0.2 M_sun, consistent with the theoretically expected values for a helium core WD with a thick, stably burning hydrogen envelope. Though nearby (at approx 150 pc), the gravitational wave strain from NLTT 11748 is likely not adequate for direct detection by the Laser Interferometer Space Antenna. Future observational efforts will determine M_1, yielding accurate WD mass-radius measurement of both components, as well as a clearer indication of the binary's fate once contact is reached.

preprint2010arXiv

Evidence for the White Dwarf Nature of Mira B

The nature of the accreting companion to Mira --- the prototypical pulsating asymptotic giant branch star --- has been a matter of debate for more than 25 years. Here we use a quantitative analysis of the rapid optical brightness variations from this companion, Mira B, which we observed with the Nickel telescope at Lick Observatory, to show that it is a white dwarf (WD). The amplitude of aperiodic optical variations on time scales of minutes to tens of minutes (approximately 0.2 mag) is consistent with that of accreting WDs in cataclysmic variables on these same time scales. It is significantly greater than that expected from an accreting main-sequence star. With Mira B identified as a WD, its ultraviolet (UV) and optical luminosities, along with constraints on the WD effective temperature from the UV, indicate that it accretes at ~1e-10 solar masses per year. We do not find any evidence that the accretion rate is higher than predicted by Bondi-Hoyle theory. The accretion rate is high enough, however, to explain the weak X-ray emission, since the accretion-disk boundary layer around a low-mass WD accreting at this rate is likely to be optically thick and therefore to emit primarily in the far or extreme UV. Furthermore, the finding that Mira B is a WD means that it has experienced, and will continue to experience nova explosions, roughly every million years. It also highlights the similarity between Mira AB and other jet-producing symbiotic binaries such as R Aquarii, CH Cygni, and MWC 560, and therefore raises the possibility that Mira B launched the recently discovered bipolar streams from this system.

preprint2010arXiv

Modules for Experiments in Stellar Astrophysics (MESA)

Stellar physics and evolution calculations enable a broad range of research in astrophysics. Modules for Experiments in Stellar Astrophysics (MESA) is a suite of open source libraries for a wide range of applications in computational stellar astrophysics. A newly designed 1-D stellar evolution module, MESA star, combines many of the numerical and physics modules for simulations of a wide range of stellar evolution scenarios ranging from very-low mass to massive stars, including advanced evolutionary phases. MESA star solves the fully coupled structure and composition equations simultaneously. It uses adaptive mesh refinement and sophisticated timestep controls, and supports shared memory parallelism based on OpenMP. Independently usable modules provide equation of state, opacity, nuclear reaction rates, and atmosphere boundary conditions. Each module is constructed as a separate Fortran 95 library with its own public interface. Examples include comparisons to other codes and show evolutionary tracks of very low mass stars, brown dwarfs, and gas giant planets; the complete evolution of a 1 Msun star from the pre-main sequence to a cooling white dwarf; the Solar sound speed profile; the evolution of intermediate mass stars through the thermal pulses on the He-shell burning AGB phase; the interior structure of slowly pulsating B Stars and Beta Cepheids; evolutionary tracks of massive stars from the pre-main sequence to the onset of core collapse; stars undergoing Roche lobe overflow; and accretion onto a neutron star. Instructions for downloading and installing MESA can be found on the project web site (http://mesa.sourceforge.net/).

preprint2010arXiv

Pulsations in Hydrogen Burning Low Mass Helium White Dwarfs

Helium core white dwarfs (WDs) with mass M <~ 0.20 Msun undergo several Gyrs of stable hydrogen burning as they evolve. We show that in a certain range of WD and hydrogen envelope masses, these WDs may exhibit g-mode pulsations similar to their passively cooling, more massive carbon/oxygen core counterparts, the ZZ Cetis. Our models with stably burning hydrogen envelopes on helium cores yield g-mode periods and period spacings longer than the canonical ZZ Cetis by nearly a factor of 2. We show that core composition and structure can be probed using seismology since the g-mode eigenfunctions predominantly reside in the helium core. Though we have not carried out a fully nonadiabatic stability analysis, the scaling of the thermal time in the convective zone with surface gravity highlights several low-mass helium WDs that should be observed in search of pulsations: NLTT 11748, SDSS J0822+2753, and the companion to PSR J1012+5307. Seismological studies of these He core WDs may prove especially fruitful, as their luminosity is related (via stable hydrogen burning) to the hydrogen envelope mass, which eliminates one model parameter.

preprint2010arXiv

Rapidly Decaying Supernova 2010X: A Candidate ".Ia" Explosion

We present the discovery, photometric and spectroscopic follow-up observations of SN 2010X (PTF 10bhp). This supernova decays exponentially with tau_d=5 days, and rivals the current recordholder in speed, SN 2002bj. SN 2010X peaks at M_r=-17mag and has mean velocities of 10,000 km/s. Our light curve modeling suggests a radioactivity powered event and an ejecta mass of 0.16 Msun. If powered by Nickel, we show that the Nickel mass must be very small (0.02 Msun) and that the supernova quickly becomes optically thin to gamma-rays. Our spectral modeling suggests that SN 2010X and SN 2002bj have similar chemical compositions and that one of Aluminum or Helium is present. If Aluminum is present, we speculate that this may be an accretion induced collapse of an O-Ne-Mg white dwarf. If Helium is present, all observables of SN 2010X are consistent with being a thermonuclear Helium shell detonation on a white dwarf, a ".Ia" explosion. With the 1-day dynamic-cadence experiment on the Palomar Transient Factory, we expect to annually discover a few such events.

preprint2010arXiv

Thermonuclear .Ia Supernovae from Helium Shell Detonations: Explosion Models and Observables

During the early evolution of an AM CVn system, helium is accreted onto the surface of a white dwarf under conditions suitable for unstable thermonuclear ignition. The turbulent motions induced by the convective burning phase in the He envelope become strong enough to influence the propagation of burning fronts and may result in the onset of a detonation. Such an outcome would yield radioactive isotopes and a faint rapidly rising thermonuclear ".Ia" supernova. In this paper, we present hydrodynamic explosion models and observable outcomes of these He shell detonations for a range of initial core and envelope masses. The peak UVOIR bolometric luminosities range by a factor of 10 (from 5e41 - 5e42 erg/s), and the R-band peak varies from M_R,peak = -15 to -18. The rise times in all bands are very rapid (<10 d), but the decline rate is slower in the red than the blue due to a secondary near-IR brightening. The nucleosynthesis primarily yields heavy alpha-chain elements (40Ca through 56Ni) and unburnt He. Thus, the spectra around peak light lack signs of intermediate mass elements and are dominated by CaII and TiII features, with the caveat that our radiative transfer code does not include the non-thermal effects necessary to produce He features.

preprint2009arXiv

Exploring the Optical Transient Sky with the Palomar Transient Factory

The Palomar Transient Factory (PTF) is a wide-field experiment designed to investigate the optical transient and variable sky on time scales from minutes to years. PTF uses the CFH12k mosaic camera, with a field of view of 7.9 deg^2 and a plate scale of 1 asec/pixel, mounted on the the Palomar Observatory 48-inch Samuel Oschin Telescope. The PTF operation strategy is devised to probe the existing gaps in the transient phase space and to search for theoretically predicted, but not yet detected, phenomena, such as fallback supernovae, macronovae, .Ia supernovae and the orphan afterglows of gamma-ray bursts. PTF will also discover many new members of known source classes, from cataclysmic variables in their various avatars to supernovae and active galactic nuclei, and will provide important insights into understanding galactic dynamics (through RR Lyrae stars) and the Solar system (asteroids and near-Earth objects). The lessons that can be learned from PTF will be essential for the preparation of future large synoptic sky surveys like the Large Synoptic Survey Telescope. In this paper we present the scientific motivation for PTF and describe in detail the goals and expectations for this experiment.

preprint2009arXiv

Helium White Dwarfs in Cataclysmic Variables

Binary evolution predicts a population of helium core (M < 0.5 Msol) white dwarfs (WDs) that are slowly accreting hydrogen-rich material from low mass main sequence or brown dwarf donors with orbital periods less than four hours. Four binaries are presently known in the Milky Way that will reach such a mass-transferring state in a few Gyr. Despite these predictions and observations of progenitor binaries, there are still no secure cases of helium core WDs among the mass-transferring cataclysmic variables (CVs). This led us to calculate the fate of He WDs once accretion begins at a rate Mdot < 1e-10 Msol/yr set by angular momentum losses. We show here that the cold He core temperatures (T_c < 1e7 K) and low Mdot result in ~ 1e-3 Msol of accumulated H-rich material at the onset of the thermonuclear runaway. Shara and collaborators noted that these large accumulated masses may lead to exceptionally long classical nova (CN) events. For a typical donor star of 0.2 Msol, such binaries will only yield a few hundred CNe, making these events rare amongst all CNe. We calculate the reheating of the accreting WD, allowing a comparison to the measured WD effective temperatures in quiescent dwarf novae and raising the possibility that WD seismology may be the best way to confirm the presence of a He WD. We also find that a very long (> 1000 yr) stable burning phase occurs after the CN outburst, potentially explaining enigmatic short orbital period supersoft sources like RX J0537-7034 (P_orb = 3.5 hr) and 1E 0035.4-7230 (P_orb = 4.1 hr).

preprint2009arXiv

Supernova Light Curves Powered by Young Magnetars

We show that energy deposited into an expanding supernova remnant by a highly magnetic (B ~ 5 x 10^14 G) neutron star spinning at an initial period of P ~ 2-20 ms can substantially brighten the light curve. For magnetars with parameters in this range, the rotational energy is released on a timescale of days to weeks, which is comparable to the effective diffusion time through the supernova remnant. The late time energy injection can then be radiated without suffering overwhelming adiabatic expansion losses. The magnetar input also produces a central bubble which sweeps ejecta into an internal dense shell, resulting in a prolonged period of nearly constant photospheric velocity in the observed spectra. We derive analytic expressions for the light curve rise time and peak luminosity as a function of B, P and the properties of the supernova ejecta that allow for direct inferences about the underlying magnetar in bright supernovae. We perform numerical radiation hydrodynamical calculations of a few specific instances and compare the resulting light curves to observed events. Magnetar activity is likely to impact more than a few percent of all core collapse supernovae, and may naturally explain some of the brightest events ever seen (e.g., SN 2005ap and SN 2008es) at L > 10^44 ergs/s.

preprint2008arXiv

The Effect of Composition on Nova Ignitions

The accretion of hydrogen-rich matter onto C/O and O/Ne white dwarfs in binary systems leads to unstable thermonuclear ignition of the accreted envelope, triggering a convective thermonuclear runaway and a subsequent classical, recurrent, or symbiotic nova. Prompted by uncertainties in the composition at the base of the accreted envelope at the onset of convection, as well as the range of abundances detected in nova ejecta, we examine the effects of varying the composition of the accreted material. For high accretion rates and carbon mass fractions < 0.002, we find that carbon, which is usually assumed to trigger the runaway via proton captures, is instead depleted and converted to 14N. Additionally, we quantify the importance of 3He, finding that convection is triggered by 3He+3He reactions for 3He mass fractions > 0.002. These different triggering mechanisms, which occur for critical abundances relevant to many nova systems, alter the amount of mass that is accreted prior to a nova, causing the nova rate to depend on accreted composition. Upcoming deep optical surveys such as Pan-STARRS-1, Pan-STARRS-4, and the Large Synoptic Survey Telescope may allow us to detect the dependence of nova rates on accreted composition. Furthermore, the burning and depletion of 3He with a mass fraction of 0.001, which is lower than necessary for triggering convection, still has an observable effect, resulting in a pre-outburst brightening in disk quiescence to > Lsun and an increase in effective temperature to 6.5e4 K for a 1.0 Msun white dwarf accreting at 1e-8 Msun/yr.

preprint2007arXiv

Thermally Stable Nuclear Burning on Accreting White Dwarfs

One of the challenges to increasing the mass of a white dwarf through accretion is the tendency for the accumulating hydrogen to ignite unstably and potentially trigger mass loss. It has been known for many years that there is a narrow range of accretion rates for which the hydrogen can burn stably, allowing for the white dwarf mass to increase as a pure helium layer accumulates. We first review the physics of stable burning, providing a clear explanation for why radiation pressure stabilization leads to a narrow range of accretion rates for stable burning near the Eddington limit, confirming the recent work of Nomoto and collaborators. We also explore the possibility of stabilization due to a high luminosity from beneath the burning layer. We then examine the impact of the beta-decay-limited ''hot'' CNO cycle on the stability of burning. Though this plays a significant role for accreting neutron stars, we find that for accreting white dwarfs, it can only increase the range of stably-burning accretion rates for metallicities below 0.01 solar metallicity.

preprint2006arXiv

Hydrodynamic Thermonuclear Runaways in Superbursts

We calculate the thermal and dynamical evolution of the surface layers of an accreting neutron star during the rise of a superburst. For the first few hours following unstable 12C ignition, the nuclear energy release is transported by convection. However, as the base temperature rises, the heating time becomes shorter than the eddy turnover time and convection becomes inefficient. This results in a hydrodynamic nuclear runaway, in which the heating time becomes shorter than the local dynamical time. Such hydrodynamic burning can drive shock waves into the surrounding layers and may be the trigger for the normal X-ray burst found to immediately precede the onset of the superburst in both cases where the Rossi X-Ray Timing Explorer was observing.

preprint2005arXiv

The discovery of X-ray binaries in the Sculptor Dwarf Spheroidal Galaxy

We report the results of a deep Chandra survey of the Sculptor dwarf spheroidal galaxy. We find five X-ray sources with $L_X$ of at least $6\times10^{33}$ ergs/sec with optical counterparts establishing them as members of Sculptor. These X-ray luminosities indicate that these sources are X-ray binaries, as no other known class of Galactic point sources can reach 0.5-8 keV luminosities this high. Finding these systems proves definitively that such objects can exist in an old stellar population without stellar collisions. Three of these objects have highly evolved optical counterparts (giants or horizontal branch stars), as do three other sources whose X-ray luminosities are in the range which includes both quiescent low mass X-ray binaries and the brightest magnetic cataclysmic variables. We predict that large area surveys of the Milky Way should also turn up large numbers of quiescent X-ray binaries.

preprint2005arXiv

The Type Ia Supernova Rate

We explore the idea that the Type Ia supernovae (SNe Ia) rate consists of two components: a prompt piece that is proportional to the star formation rate (SFR) and an extended piece that is proportional to the total stellar mass. We fit the parameters of this model to the local observations of Mannucci and collaborators and then study its impact on three important problems. On cosmic scales, the model reproduces the observed SNe Ia rate density below z=1, and predicts that it will track the measured SFR density at higher redshift, reaching a value of 1-3.5 X 10^-4 per yr per Mpc^3 at z=2. In galaxy clusters, a large prompt contribution helps explain the iron content of the intracluster medium. Within the Galaxy, the model reproduces the observed stellar [O/Fe] abundance ratios if we allow a short (approximately 0.7 Gyr) delay in the prompt component. Ongoing medium-redshift SN surveys will yield more accurate parameters for our model

preprint1998arXiv

Gravitational Radiation and Rotation of Accreting Neutron Stars

Recent discoveries by the Rossi X-Ray Timing Explorer indicate that most of the rapidly accreting and weakly magnetic neutron stars in the Galaxy are rotating at spin frequencies greater than 250 Hz. Remarkably, they all rotate in a narrow range of frequencies. I suggest that these stars rotate fast enough so that, on average, the angular momentum added by accretion is lost to gravitational radiation. The strong spin frequency dependence of the angular momentum loss rate from gravitational radiation then provides a natural reason for similar spin frequencies. Provided that the interior temperature has a large scale asymmetry misaligned from the spin axis, then the temperature sensitive electron captures in the deep crust can provide the quadrupole needed to reach this limiting situation at 300 Hz. This quadrupole is only present during accretion and makes it difficult to form radio pulsars spinning more rapidly than 600-800 Hz via rapid accretion. The gravity wave strength is <10^{-26} from most neutron stars and >2 X 10^{-26} for Sco X-1. Prior knowledge of the position, spin frequency and orbital periods will allow for deep searches for these periodic signals with gravitational wave interferometers (LIGO, VIRGO and the dual-recycled GEO 600 detector) and experimenters need to take such sources into account. Sco X-1 will most likely be detected first.

Lars Bildsten

What is connected

Connect this record

See the researcher in context

Building this map preview

69 published item(s)

3D Hydrodynamics of Pre-supernova Outbursts in Convective Red Supergiant Envelopes

Discovery of a double detonation thermonuclear supernova progenitor

Numerical Simulations of Convective 3-Dimensional Red Supergiant Envelopes

Shock Breakout in 3-Dimensional Red Supergiant Envelopes

Still Brighter than Pre-Explosion, SN 2012Z Did Not Disappear: Comparing Hubble Space Telescope Observations a Decade Apart

Stochastic Low Frequency Variability in 3-Dimensional Radiation Hydrodynamical Models of Massive Star Envelopes

A Massive Star's Dying Breaths: Pulsating Red Supergiants and Their Resulting Type IIP Supernovae

A new class of Roche lobe-filling hot subdwarf binaries

Comparing Moment-Based and Monte Carlo Methods of Radiation Transport Modeling for Type II-Plateau Supernova Light Curves

The first ultracompact Roche lobe-filling hot subdwarf binary

The Value of Progenitor Radius Measurements for Explosion Modeling of Type II-Plateau Supernovae

Variability of massive stars in M31 from the Palomar Transient Factory

Remnants of Subdwarf Helium Donor Stars Ejected from Close Binaries with Thermonuclear Supernovae

Asteroseismic signatures of evolving internal stellar magnetic fields

Carbon Shell or Core Ignitions in White Dwarfs Accreting from Helium Stars

Late-time Spectroscopy of Type Iax Supernovae

Rapidly Rising Transients in the Supernova - Superluminous Supernova Gap

Seismology of Rapidly Rotating Accreting White Dwarfs

Thermal Runaway During the Evolution of ONeMg Cores towards Accretion-Induced Collapse

Turbulent Chemical Diffusion in Convectively Bounded Carbon Flames

AM Canum Venaticorum Progenitors with Helium Star Donors and the resultant Explosions

Asteroseismology Can Reveal Strong Internal Magnetic Fields in Red Giant Stars

Local Radiation Hydrodynamic Simulations of Massive Star Envelopes at the Iron Opacity Peak

Magnetar Driven Shock Breakout and Double Peaked Supernova Light Curves

Photoionization Heating of Nova Ejecta by the Post-Outburst Supersoft Source

Population of post-nova supersoft X-ray sources

An Accreting White Dwarf near the Chandrasekhar Limit in the Andromeda Galaxy

Angular momentum transport within evolved low-mass stars

Identification of the Optical Counterpart of Fermi Black Widow Millisecond Pulsar PSR J1544+4937

Possible Detection of the Stellar Donor or Remnant for the Type Iax Supernova 2008ha

The Ignition of Carbon Detonations via Converging Shock Waves in White Dwarfs

X-ray diagnostics of chemical composition of the accretion disk and donor star in UCXBs II: XMM-Newton observations

Direct evidence for a supernova interacting with a large amount of hydrogen-free circumstellar material

Hydrogen Burning on Accreting White Dwarfs: Stability, Recurrent Novae, and the Post-Novae Supersoft Source

Modules for Experiments in Stellar Astrophysics (MESA): Giant Planets, Oscillations, Rotation, and Massive Stars

Properties of the Eclipsing Double-White Dwarf Binary NLTT 11748

R Coronae Borealis Stars in M31 from the Palomar Transient Factory

The Effects of Curvature and Expansion on Helium Detonations on White Dwarf Surfaces

X-ray diagnostics of chemical composition of the accretion disk and donor star in ultra-compact X-ray binaries

Circumstellar Shell Formation in Symbiotic Recurrent Novae

Laterally Propagating Detonations in Thin Helium Layers on Accreting White Dwarfs

MESA Models of Classical Nova Outbursts: The Multicycle Evolution and Effects of Convective Boundary Mixing

Orbital Evolution of Compact White Dwarf Binaries

The Long-Term Evolution of Double White Dwarf Mergers

The Observational Signatures of Convectively Excited Internal Gravity Modes in Main Sequence Stars

A Search for Pulsations in Helium White Dwarfs

Acoustic Signatures of the Helium Core Flash

Calcium-rich gap transients in the remote outskirts of galaxies

Constraints on the Progenitor System of the Type Ia Supernova SN 2011fe/PTF11kly

PTF10iya: A short-lived, luminous flare from the nuclear region of a star-forming galaxy

Supernova 2011fe from an Exploding Carbon-Oxygen White Dwarf Star

A ground-based measurement of the relativistic beaming effect in a detached double WD binary

Clearing the Gas from Globular Clusters & Dwarf Spheroidals with Classical Novae

Diffusive Nuclear Burning of Helium on Neutron Stars

Discovery of the Eclipsing Detached Double White Dwarf Binary NLTT 11748

Evidence for the White Dwarf Nature of Mira B

Modules for Experiments in Stellar Astrophysics (MESA)

Pulsations in Hydrogen Burning Low Mass Helium White Dwarfs

Rapidly Decaying Supernova 2010X: A Candidate ".Ia" Explosion

Thermonuclear .Ia Supernovae from Helium Shell Detonations: Explosion Models and Observables

Exploring the Optical Transient Sky with the Palomar Transient Factory

Helium White Dwarfs in Cataclysmic Variables

Supernova Light Curves Powered by Young Magnetars

The Effect of Composition on Nova Ignitions

Thermally Stable Nuclear Burning on Accreting White Dwarfs

Hydrodynamic Thermonuclear Runaways in Superbursts

The discovery of X-ray binaries in the Sculptor Dwarf Spheroidal Galaxy

The Type Ia Supernova Rate

Gravitational Radiation and Rotation of Accreting Neutron Stars