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.

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.

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.

(Abridged) Using the Zwicky Transient Facility alert stream, we are conducting a large campaign to spectroscopically classify all transients occurring in galaxies in the Census of the Local Universe (CLU) catalog. The aim of the experiment is to construct a spectroscopically complete, volume-limited sample of transients coincident within 100" of CLU galaxies out to 200 Mpc, and to a depth of 20 mag. We describe the survey design and spectroscopic completeness from the first 16 months of operations. We present results from a systematic search for Calcium rich gap transients in the sample of 22 low luminosity (peak absolute magnitude $M > -17$), hydrogen poor events found in the experiment (out of 754 spectroscopically classified SNe). We report the detection of eight Calcium rich gap transients, and constrain their volumetric rate to be at least $\approx 15\pm5$% of the SN Ia rate. Combining this sample with ten events from the literature, we find a likely continuum of spectroscopic properties ranging from events with SN Ia-like features (Ca-Ia objects) to SN Ib/c-like features (Ca-Ib/c objects) at peak light. Within the Ca-Ib/c events, we find two populations of events distinguished by their red ($g - r \approx 1.5$ mag) or green ($g - r \approx 0.5$ mag) spectral colors at $r$-band peak, wherein redder events show strong line blanketing signatures, slower light curves, weaker He lines and lower [Ca II]/[O I] in the nebular phase. Together, we find that the spectroscopic continuum, volumetric rates and striking old environments are consistent with the explosive burning of He shells on low mass white dwarfs. We posit that Ca-Ia and red Ca-Ib/c objects are consistent with the double detonation of He shells with high He burning efficiency, while green Ca-Ib/c objects could arise from less efficient He burning scenarios such as detonations in low density He shells or He shell deflagrations.

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.

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.

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.

Astronomy, as a field, has long encouraged the development of free, open digital infrastructure (e.g., National Research Council 2010, 2011). Examples range from simple scripts that enable individual scientific research, through software instruments for entire communities, to data reduction pipelines for telescope operations at national facilities. As with the digital infrastructure of our larger society today (e.g., Eghbal 2016), nearly all astronomical research relies on free, open source software (FOSS) written and maintained by a small number of developers. And like the physical infrastructure of roads or bridges, digital infrastructure needs regular upkeep and maintenance (e.g., Eghbal 2016). In astronomy, financial support for maintaining existing digital infrastructure is generally much harder to secure than funding for developing new digital infrastructures that promise new science. Sustaining astronomy's digital infrastructure is a new topic for many, the sustainability challenges are not always widely known...

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.

General relativity predicts that short orbital period binaries emit significant gravitational radiation, and the upcoming Laser Interferometer Space Antenna (LISA) is expected to detect tens of thousands of such systems; however, few have been identified, and only one is eclipsing--the double white dwarf binary SDSS J065133.338+284423.37, which has an orbital period of 12.75 minutes. Here, we report the discovery of an eclipsing double white dwarf binary system with an orbital period of only 6.91 minutes, ZTF J153932.16+502738.8. This system has an orbital period close to half that of SDSS J065133.338+284423.37 and an orbit so compact that the entire binary could fit within the diameter of the planet Saturn. The system exhibits a deep eclipse, and a double-lined spectroscopic nature. We observe rapid orbital decay, consistent with that expected from general relativity. ZTF J153932.16+502738.8 is a significant source of gravitational radiation close to the peak of LISA's sensitivity, and should be detected within the first week of LISA observations.

We present advances in modeling Type IIP supernovae using MESA for evolution to shock breakout coupled with STELLA for generating light and radial velocity curves. Explosion models and synthetic light curves can be used to translate observable properties of supernovae (such as the luminosity at day 50 and the duration of the plateau, as well as the observable quantity $ET$, defined as the time-weighted integrated luminosity that would have been generated if there was no ${\rm ^{56}Ni}$ in the ejecta) into families of explosions which produce the same light curve and velocities on the plateau. These predicted families of explosions provide a useful guide towards modeling observed SNe, and can constrain explosion properties when coupled with other observational or theoretical constraints. For an observed supernova with a measured ${\rm ^{56}Ni}$ mass, breaking the degeneracies within these families of explosions (ejecta mass, explosion energy, and progenitor radius) requires independent knowledge of one parameter. We expect the most common case to be a progenitor radius measurement for a nearby supernova. We show that ejecta velocities inferred from the Fe II$\lambda$ 5169 line measured during the majority of the plateau phase provide little additional information about explosion characteristics. Only during the initial shock cooling phase can photospheric velocity measurements potentially aid in unraveling light curve degeneracies.

Using high-cadence observations from the Zwicky Transient Facility at low Galactic latitudes, we have discovered a new class of pulsating, hot, compact stars. We have found four candidates, exhibiting blue colors ($g-r\leq-0.1$ mag), pulsation amplitudes of $>5\%$, and pulsation periods of $200 - 475$ sec. Fourier transforms of the lightcurves show only one dominant frequency. Phase-resolved spectroscopy for three objects reveals significant radial velocity, T$_{\rm eff}$ and log(g) variations over the pulsation cycle, consistent with large amplitude radial oscillations. The mean T$_{\rm eff}$ and log(g) for these stars are consistent with hot subdwarf B (sdB) effective temperatures and surface gravities. We calculate evolutionary tracks using MESA and adiabatic pulsations using GYRE for low-mass helium-core pre-white dwarfs and low mass helium-burning stars. Comparison of low-order radial oscillation mode periods with the observed pulsation periods show better agreement with the pre-white dwarf models. Therefore, we suggest that these new pulsators and Blue Large-Amplitude Pulsators (BLAPs) could be members of the same class of pulsators, composed of young $\approx0.25-0.35$ M$_\odot$ helium-core pre-white dwarfs.

We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). RSP is a new functionality in MESAstar that models the non-linear radial stellar pulsations that characterize RR Lyrae, Cepheids, and other classes of variable stars. We significantly enhance numerical energy conservation capabilities, including during mass changes. For example, this enables calculations through the He flash that conserve energy to better than 0.001 %. To improve the modeling of rotating stars in MESA, we introduce a new approach to modifying the pressure and temperature equations of stellar structure, and a formulation of the projection effects of gravity darkening. A new scheme for tracking convective boundaries yields reliable values of the convective-core mass, and allows the natural emergence of adiabatic semiconvection regions during both core hydrogen- and helium-burning phases. We quantify the parallel performance of MESA on current generation multicore architectures and demonstrate improvements in the computational efficiency of radiative levitation. We report updates to the equation of state and nuclear reaction physics modules. We briefly discuss the current treatment of fallback in core-collapse supernova models and the thermodynamic evolution of supernova explosions. We close by discussing the new MESA Testhub software infrastructure to enhance source-code development.

Many isolated white dwarfs (WDs) show spectral evidence of atmospheric metal pollution. Since heavy element sedimentation timescales are short, this most likely indicates ongoing accretion. Accreted metals encounter a variety of mixing processes at the WD surface: convection, gravitational sedimentation, overshoot, and thermohaline instability. We present MESA WD models that explore each of these processes and their implications for inferred accretion rates. We provide diffusion timescales for many individual metals, and we quantify the regimes in which thermohaline mixing dominates over gravitational sedimentation in setting the effective settling rate of the heavy elements. We build upon and confirm earlier work finding that accretion rates as high as $10^{13} \, \rm g \, s^{-1}$ are needed to explain observed pollution in DA WDs for $T_{\rm eff} > 15,000 \, \rm K$, and we provide tabulated results from our models that enable accretion rate inferences from observations of polluted DA WDs. If these rates are representative of young WDs, we estimate that the total mass of planetesimal material accreted over a WD lifetime may be as high as $10^{28} \, \rm g$, though this estimate is susceptible to potential selection biases and uncertainties about the nature of disk processes that supply accretion to the WD surface. We also find that polluted DB WDs experience much less thermohaline mixing than DA WDs, and we do not expect thermohaline instability to be active for polluted DB WDs with $T_{\rm eff} < 18,000 \, \rm K$.

The detonation of a helium shell on a white dwarf has been proposed as a possible explosion triggering mechanism for Type Ia supernovae. Here, we report ZTF 18aaqeasu (SN 2018byg/ATLAS 18pqq), a peculiar Type I supernova, consistent with being a helium-shell double-detonation. With a rise time of $\approx 18$ days from explosion, the transient reached a peak absolute magnitude of $M_R \approx -18.2$ mag, exhibiting a light curve akin to sub-luminous SN 1991bg-like Type Ia supernovae, albeit with an unusually steep increase in brightness within a week from explosion. Spectra taken near peak light exhibit prominent Si absorption features together with an unusually red color ($g-r \approx 2$ mag) arising from nearly complete line blanketing of flux blue-wards of 5000 \AA. This behavior is unlike any previously observed thermonuclear transient. Nebular phase spectra taken at and after $\approx 30$ days from peak light reveal evidence of a thermonuclear detonation event dominated by Fe-group nucleosynthesis. We show that the peculiar properties of ZTF 18aaqeasu are consistent with the detonation of a massive ($\approx 0.15$ M$_\odot$) helium shell on a sub-Chandrasekhar mass ($\approx 0.75$ M$_\odot$) white dwarf after including mixing of $\approx 0.2$ M$_\odot$ of material in the outer ejecta. These observations provide evidence of a likely rare class of thermonuclear supernovae arising from detonations of massive helium shells.

(Abridged) Stars more massive than $20-30M_{\odot}$ are so luminous that the radiation force on the cooler, more opaque outer layers can balance or exceed the force of gravity. These near or super-Eddington outer envelopes represent a long standing challenge for calculating the evolution of massive stars in one dimension, a situation that limits our understanding of the stellar progenitors of some of the most exciting and energetic explosions in the universe. In particular, the proximity to the Eddington limit has been the suspected cause for the variability, large mass loss rate and giant eruptions of an enigmatic class of massive stars: the luminous blue variables (LBVs). When in quiescence, LBVs are usually found on the hot ($T_{eff} \approx 2 - 4 \times 10^4$ K) S Dor instability strip. While in outburst, most LBVs stay on the cold S Dor instability strip with a $T_{eff} \approx 9000$ K. Here we show that physically realistic three dimensional global radiation hydrodynamic simulations of radiation dominated massive stars with the largest supercomputers in the world naturally reproduce many observed properties of LBVs, specifically their location in the Hertzsprung-Russell (HR) diagram and their episodic mass loss with rates of $10^{-7}-10^{-5} M_{\odot}/yr$. The helium opacity peak is found to play an important role to determine these properties, which is not realized in the traditional one dimensional models of massive stars. The simulations also predict that convection causes irregular envelope oscillations and 10-30% brightness variations on a typical timescale of a few days. The variability is more prominent in our models that are on the cool part of the S Dor instability. These calculations pave the way to a quantitative understanding of the structure, stability and the dominant mode of mass loss of massive stars.

Many isolated, old white dwarfs (WDs) show surprising evidence of metals in their photospheres. Given that the timescale for gravitational sedimentation is astronomically short, this is taken as evidence for ongoing accretion, likely of tidally disrupted planetesimals. The rate of such accretion, $\dot M_{\rm acc}$, is important to constrain, and most modeling of this process relies on assuming an equilibrium between diffusive sedimentation and metal accretion supplied to the WD's surface convective envelope. Building on earlier work of Deal and collaborators, we show that high $\dot M_{\rm acc}$ models with only diffusive sedimentation are unstable to thermohaline mixing and that models which account for the enhanced mixing from the active thermohaline instability require larger accretion rates, sometimes reaching $\dot M_{\rm acc} \approx 10^{13} \, \rm g \, s^{-1}$ to explain observed Calcium abundances. We present results from a grid of MESA models that include both diffusion and thermohaline mixing. These results demonstrate that both mechanisms are essential for understanding metal pollution across the range of polluted WDs with hydrogen atmospheres. Another consequence of active thermohaline mixing is that the observed metal abundance ratios are identical to accreted material.

Most massive stars end their lives as Red Supergiants (RSGs), a short-lived evolution phase when they are known to pulsate with varying amplitudes. The RSG period-luminosity (PL) relation has been measured in the Milky Way, the Magellanic Clouds and M33 for about 120 stars in total. Using over 1500 epochs of R-band monitoring from the Palomar Transient Factory (PTF) survey over a five-year period, we study the variability of 255 spectroscopically cataloged RSGs in M31. We find that all RGSs brighter than M_K~ -10 mag (log(L/L_sun)>4.8) are variable at dm_R>0.05 mag. Our period analysis finds 63 with significant pulsation periods. Using the periods found and the known values of M_K for these stars, we derive the RSG PL relation in M31 and show that it is consistent with those derived earlier in other galaxies of different metallicities. We also detect, for the first time, a sequence of likely first-overtone pulsations. Comparison to stellar evolution models from MESA confirms the first overtone hypothesis and indicates that the variable stars in this sample have 12 M_sun<M<24 M_sun. As these RSGs are the immediate progenitors to Type II-P core-collapse supernovae (SNe), we also explore the implication of their variability in the initial-mass estimates for SN progenitors based on archival images of the progenitors. We find that this effect is small compared to the present measurement errors.

After an optical peak, a classical or recurrent nova settles into a brief (days to years) period of quasi-stable thermonuclear burning in a compact configuration nearly at the white dwarf (WD) radius. During this time, the underlying WD becomes visible as a strong emitter of supersoft X-rays. Observations during this phase have revealed oscillations in the X-ray emission with periods on the order of tens of seconds. A proposed explanation for the source of these oscillations are internal gravity waves excited by nuclear reactions at the base of the hydrogen-burning layer. In this work, we present the first models exhibiting unstable surface $g$-modes with periods similar to oscillation periods found in galactic novae. However, when comparing mode periods of our models to the observed oscillations of several novae, we find that the modes which are excited have periods shorter than that observed.

We update the capabilities of the software instrument Modules for Experiments in Stellar Astrophysics (MESA) and enhance its ease of use and availability. Our new approach to locating convective boundaries is consistent with the physics of convection, and yields reliable values of the convective core mass during both hydrogen and helium burning phases. Stars with $M<8\,{\rm M_\odot}$ become white dwarfs and cool to the point where the electrons are degenerate and the ions are strongly coupled, a realm now available to study with MESA due to improved treatments of element diffusion, latent heat release, and blending of equations of state. Studies of the final fates of massive stars are extended in MESA by our addition of an approximate Riemann solver that captures shocks and conserves energy to high accuracy during dynamic epochs. We also introduce a 1D capability for modeling the effects of Rayleigh-Taylor instabilities that, in combination with the coupling to a public version of the STELLA radiation transfer instrument, creates new avenues for exploring Type II supernovae properties. These capabilities are exhibited with exploratory models of pair-instability supernova, pulsational pair-instability supernova, and the formation of stellar mass black holes. The applicability of MESA is now widened by the capability of importing multi-dimensional hydrodynamic models into MESA. We close by introducing software modules for handling floating point exceptions and stellar model optimization, and four new software tools -- MESAWeb, MESA-Docker, pyMESA, and mesastar.org -- to enhance MESA's education and research impact.

We study the evolution and final outcome of long-lived (${\approx}10^5$ years) remnants from the merger of a He white dwarf (WD) with a more massive C/O or O/Ne WD. Using Modules for Experiments in Stellar Astrophysics ($\texttt{MESA}$), we show that these remnants have a red giant configuration supported by steady helium burning, adding mass to the WD core until it reaches $M_{\rm core}\approx 1.12-1.20 M_\odot$. At that point, the base of the surface convection zone extends into the burning layer, mixing the helium burning products (primarily carbon and magnesium) throughout the convective envelope. Further evolution depletes the convective envelope of helium, and dramatically slows the mass increase of the underlying WD core. The WD core mass growth re-initiates after helium depletion, as then an uncoupled carbon burning shell is ignited and proceeds to burn the fuel from the remaining metal-rich extended envelope. For large enough initial total merger masses, O/Ne WD cores would experience electron-capture triggered collapse to neutron stars (NSs) after growing to near Chandrasekhar mass ($M_{\rm Ch}$). Massive C/O WD cores could suffer the same fate after a carbon-burning flame converts them to O/Ne. The NS formation would release ${\approx}10^{50}$ ergs into the remaining extended low mass envelope. Using the STELLA radiative transfer code, we predict the resulting optical light curves from these exploded envelopes. Reaching absolute magnitudes of $M_V\approx -17$, these transients are bright for about one week, and have many features of the class of luminous, rapidly evolving transients studied by Drout and collaborators.

Every supernova hitherto observed has been considered to be the terminal explosion of a star. Moreover, all supernovae with absorption lines in their spectra show those lines decreasing in velocity over time, as the ejecta expand and thin, revealing slower moving material that was previously hidden. In addition, every supernova that exhibits the absorption lines of hydrogen has one main light-curve peak, or a plateau in luminosity, lasting approximately 100 days before declining. Here we report observations of iPTF14hls, an event that has spectra identical to a hydrogen-rich core-collapse supernova, but characteristics that differ extensively from those of known supernovae. The light curve has at least five peaks and remains bright for more than 600 days; the absorption lines show little to no decrease in velocity; and the radius of the line-forming region is more than an order of magnitude bigger than the radius of the photosphere derived from the continuum emission. These characteristics are consistent with a shell of several tens of solar masses ejected by the star at supernova-level energies a few hundred days before a terminal explosion. Another possible eruption was recorded at the same position in 1954. Multiple energetic pre-supernova eruptions are expected to occur in stars of 95-130 solar masses, which experience the pulsational pair instability. That model, however, does not account for the continued presence of hydrogen, or the energetics observed here. Another mechanism for the violent ejection of mass in massive stars may be required.

We report the discovery of the ultracompact hot subdwarf (sdOB) binary OW J074106.0-294811.0 with an orbital period of P$_{\rm orb}=44.66279\pm1.16\times10^{-4}$ min, making it the most compact hot subdwarf binary known. Spectroscopic observations using the VLT, Gemini and Keck telescopes revealed a He-sdOB primary with an intermediate helium abundance, T$_{\rm eff}=39 400\pm500$ K and log(g)=$5.74\pm0.09$. High signal-to-noise ratio lightcurves show strong ellipsoidal modulation resulting in a derived sdOB mass $M_{\rm sdOB}=0.23\pm0.12$ M$_\odot$ with a WD companion ($M_{\rm WD}=0.72\pm0.17$ M$_\odot$). The mass ratio was found to be $q = M_{\rm sdOB}/M_{\rm WD}=0.32\pm0.10$. The derived mass for the He-sdOB is inconsistent with the canonical mass for hot sbudwarfs of $\approx0.47$ M$_\odot$. To put constraints on the structure and evolutionary history of the sdOB star we compared the derived T$_{\rm eff}$, log(g) and sdOB mass to evolutionary tracks of helium stars and helium white dwarfs calculated with Modules for Experiments in Stellar Astrophysics (MESA). We find that the best fitting model is a helium white dwarf with a mass of $0.320$ M$_\odot$, which left the common envelope ${\approx}1.1$ Myr ago, is consistent with the observations. As a helium white dwarf with a massive white dwarf companion the object will reach contact in 17.6 Myr at an orbital period of 5 min. Depending on the spin-orbit synchronization timescale the object will either merge to form an R CrB star or end up as a stably accreting AM CVn-type system with a helium white dwarf donor.

The bright subdwarf-O star (sdO), HD 49798, is in a 1.55 day orbit with a compact companion that is spinning at 13.2 seconds. Using the measurements of the effective temperature ($T_{\rm eff}$), surface gravity ($\log g$), and surface abundances of the sdO, we construct models to study the evolution of this binary system using Modules for Experiments in Stellar Astrophysics ($\texttt{MESA}$). Previous studies of the compact companion have disagreed on whether it is a white dwarf (WD) or a neutron star (NS). From the published measurements of the companion's spin and spin-up rate, we agree with Mereghetti and collaborators that a NS companion is more likely. However, since there remains the possibility of a WD companion, we use our constructed $\texttt{MESA}$ models to run simulations with both WD and NS companions that help us constrain the past and future evolution of this system. If it presently contains a NS, the immediate mass transfer evolution upon Roche lobe (RL) filling will lead to mass transfer rates comparable to that implied in ultraluminous X-ray sources (ULXs). Depending on the rate of angular momentum extraction via a wind, the fate of this system is either a wide ($P_{\rm orb}{\approx} 3$ day) intermediate mass binary pulsar (IMPB) with a relatively rapidly spinning NS (${\approx} 0.3$ s) and a high mass WD (${\approx} 0.9 M_\odot$), or a solitary millisecond pulsar (MSP).

We study the evolution of accreting oxygen-neon (ONe) white dwarfs (WDs), with a particular emphasis on the effects of the presence of the carbon-burning products $\mathrm{^{23}Na}$ and $\mathrm{^{25}Mg}$. These isotopes lead to substantial cooling of the WD via the $\mathrm{^{25}Mg}$-$\mathrm{^{25}Na}$, $\mathrm{^{23}Na}$-$\mathrm{^{23}Ne}$, and $\mathrm{^{25}Na}$-$\mathrm{^{25}Ne}$ Urca pairs. We derive an analytic formula for the peak Urca-process cooling rate and use it to obtain a simple expression for the temperature to which the Urca process cools the WD. Our estimates are equally applicable to accreting carbon-oxygen WDs. We use the Modules for Experiments in Stellar Astrophysics (MESA) stellar evolution code to evolve a suite of models that confirm these analytic results and demonstrate that Urca-process cooling substantially modifies the thermal evolution of accreting ONe WDs. Most importantly, we show that MESA models with lower temperatures at the onset of the $\mathrm{^{24}Mg}$ and $\mathrm{^{24}Na}$ electron captures develop convectively unstable regions, even when using the Ledoux criterion. We discuss the difficulties that we encounter in modeling these convective regions and outline the potential effects of this convection on the subsequent WD evolution. For models in which we do not allow convection to operate, we find that oxygen ignites around a density of $\log(\rho_{\rm c}/\rm g\,cm^{-3}) \approx 9.95$, very similar to the value without Urca cooling. Nonetheless, the inclusion of the effects of Urca-process cooling is an important step in producing progenitor models with more realistic temperature and composition profiles which are needed for the evolution of the subsequent oxygen deflagration and hence for studies of the signature of accretion-induced collapse.

White dwarfs (WDs) that accrete helium at rates $\sim 10^{-8} \, M_\odot \, \rm yr^{-1}$, such as those in close binaries with sdB stars, can accumulate large ($\gtrsim 0.1 \, M_\odot$) helium envelopes which are likely to detonate. We perform binary stellar evolution calculations of sdB+WD binary systems with MESA, incorporating the important reaction chain $^{14}{\rm N}(e^-, \nu) {^{14}{\rm C}} ( \alpha, \gamma) {^{18}{\rm O}}$ (NCO), including a recent measurement for the ${^{14}{\rm C}} ( \alpha, \gamma) {^{18}{\rm O}}$ rate. In large accreted helium shells, the NCO reaction chain leads to ignitions at the dense base of the freshly accreted envelope, in contrast to $3\alpha$ ignitions which occur away from the base of the shell. In addition, at these accretion rates, the shells accumulate on a timescale comparable to their thermal time, leading to an enhanced sensitivity of the outcome on the accretion rate history. Hence, time dependent accretion rates from binary stellar evolution are necessary to determine the helium layer mass at ignition. We model the observed sdB+WD system ${\rm CD}\,-30^\circ 11223$ and find that the inclusion of these effects predicts ignition of a $0.153 \, M_\odot$ helium shell, nearly a factor of two larger than previous predictions. A shell with this mass will ignite dynamically, a necessary condition for a helium shell detonation.

Accretion-induced collapse (AIC) occurs when an O/Ne white dwarf (WD) grows to nearly the Chandrasekhar mass ($M_{\rm Ch}$), reaching central densities that trigger electron captures in the core. Using Modules for Experiments in Stellar Astrophysics ($\texttt{MESA}$), we present the first true binary simulations of He star + O/Ne WD binaries, focusing on a $1.5 M_\odot$ He star in a 3 hour orbital period with $1.1-1.3 M_\odot$ O/Ne WDs. The helium star fills its Roche lobe after core helium burning is completed and donates helium on its thermal timescale to the WD, $\dot{M}\approx3\times10^{-6} M_\odot$/yr, a rate high enough that the accreting helium burns stably on the WD. The accumulated carbon/oxygen ashes from the helium burning undergo an unstable shell flash that initiates an inwardly moving carbon burning flame. This flame is only quenched when it runs out of carbon at the surface of the original O/Ne core. Subsequent accumulation of fresh carbon/oxygen layers also undergo thermal instabilities, but no mass loss is triggered, allowing $M_{\rm WD}\rightarrow M_{\rm Ch}$, triggering the onset of AIC. We also discuss the scenario of accreting C/O WDs that experience shell carbon ignitions to become O/Ne WDs, and then, under continuing mass transfer, lead to AIC. Studies of the AIC event rate using binary population synthesis should include all of these channels, especially this latter channel, which has been previously neglected but might dominate the rate.

We substantially update the capabilities of the open-source software instrument Modules for Experiments in Stellar Astrophysics (MESA). MESA can now simultaneously evolve an interacting pair of differentially rotating stars undergoing transfer and loss of mass and angular momentum, greatly enhancing the prior ability to model binary evolution. New MESA capabilities in fully coupled calculation of nuclear networks with hundreds of isotopes now allow MESA to accurately simulate advanced burning stages needed to construct supernova progenitor models. Implicit hydrodynamics with shocks can now be treated with MESA, enabling modeling of the entire massive star lifecycle, from pre-main sequence evolution to the onset of core collapse and nucleosynthesis from the resulting explosion. Coupling of the GYRE non-adiabatic pulsation instrument with MESA allows for new explorations of the instability strips for massive stars while also accelerating the astrophysical use of asteroseismology data. We improve treatment of mass accretion, giving more accurate and robust near-surface profiles. A new MESA capability to calculate weak reaction rates "on-the-fly" from input nuclear data allows better simulation of accretion induced collapse of massive white dwarfs and the fate of some massive stars. We discuss the ongoing challenge of chemical diffusion in the strongly coupled plasma regime, and exhibit improvements in MESA that now allow for the simulation of radiative levitation of heavy elements in hot stars. We close by noting that the MESA software infrastructure provides bit-for-bit consistency for all results across all the supported platforms, a profound enabling capability for accelerating MESA's development.

We use three dimensional radiation magneto-hydrodynamic simulations to study the effects of magnetic fields on the energy transport and structure of radiation pressure dominated main sequence massive star envelopes at the region of the iron opacity peak. We focus on the regime where the local thermal timescale is shorter than the dynamical timescale, corresponding to inefficient convective energy transport. We begin with initially weak magnetic fields relative to the thermal pressure, from 100-1000G in differing geometries. The unstable density inversion amplifies the magnetic field, increasing the magnetic energy density to values close to equipartition with the turbulent kinetic energy density. By providing pressure support, the magnetic field's presence significantly increases the density fluctuations in the turbulent envelope, thereby enhancing the radiative energy transport by allowing photons to diffuse out through low density regions. Magnetic buoyancy brings small scale magnetic fields to the photosphere and increases the vertical energy transport with the energy advection velocity proportional to the Alfv\'en velocity, although in all cases we study photon diffusion still dominates the energy transport. The increased radiative and advective energy transport causes the stellar envelope to shrink by several scale heights. We also find larger turbulent velocity fluctuations compared to the purely hydrodynamic case, reaching $\approx$ 100 km/s at the stellar photosphere. The photosphere also shows vertical oscillations with similar averaged velocities and periods of a few hours. The increased turbulent velocity and oscillations will have strong impacts on the line broadening and periodic signals in massive stars.

We present the discovery of the hot subdwarf B star (sdB) binary PTF1 J082340.04+081936.5. The system has an orbital period P$_{\rm orb}=87.49668(1)$ min (0.060761584(10) days), making it the second-most compact sdB binary known. The lightcurve shows ellipsoidal variations. Under the assumption that the sdB primary is synchronized with the orbit, we find a mass $M_{\rm sdB}=0.45^{+0.09}_{-0.07}$ M$_\odot$, a companion white dwarf mass $M_{\rm WD}=0.46^{+0.12}_{-0.09}$ M$_\odot$ and a mass ratio $q = \frac{M_{\rm WD}}{M_{\rm sdB}}=1.03^{+0.10}_{-0.08}$. The future evolution was calculated using the MESA stellar evolution code. Adopting a canonical sdB mass of $M_{\rm sdB}=0.47$ M$_\odot$, we find that the sdB still burns helium at the time it will fill its Roche lobe if the orbital period was less than 106 min at the exit from the last common envelope phase. For longer common envelope exit periods the sdB will have stopped burning helium and turned into a C/O white dwarf at the time of contact. Comparing the spectroscopically derived log(g) and $T_{\rm eff}$ with our MESA models, we find that an sdB model with a hydrogen envelope mass of $5\times10^{-4} M_\odot$ matches the measurements at a post-common envelope age of 94 Myr, corresponding to a post-common envelope orbital period of 109 min which is close to the limit to start accretion while the sdB is still burning helium.

A hybrid C/O/Ne white dwarf (WD) -- an unburned C/O core surrounded by an O/Ne/Na mantle -- can be formed if the carbon flame is quenched in a super-AGB (SAGB) star or white dwarf merger remnant. We show that this segregated hybrid structure becomes unstable to rapid mixing within 2,000 years of the onset of WD cooling. Carbon burning includes a weak reaction that removes electrons, resulting in a lower electron-to-baryon ratio ($Y_{\rm e}$) in the regions processed by carbon burning compared to the unburned C/O core, making the O/Ne mantle denser than the C/O core as the WD cools. This is unstable to efficient mixing. We use the results of $\texttt{MESA}$ models with different size C/O cores to quantify the rate at which the cores mix with the mantle as they cool. In all cases, we find that the WDs undergo significant core/mantle mixing on timescales shorter than the time available to grow the WD to the Chandrasekhar mass ($M_{\rm Ch}$) by accretion. As a result, hybrid WDs that reach $M_{\rm Ch}$ due to later accretion will have lower central carbon fractions than assumed thus far. We briefly discuss the implications of these results for the possibility of Type Ia supernovae from hybrid WDs.

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, $\kappa$, i.e., when $D_{\rm t}>\kappa$. However, we find $D_{\rm t}<\kappa$ 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.

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 $\rho_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.

High-quality collections of Type II supernova (SN) light curves are scarce because they evolve for hundreds of days, making follow-up observations time consuming and often extending over multiple observing seasons. In light of these difficulties, the diversity of SNe II is not fully understood. Here we present ultraviolet and optical photometry of 12 SNe II monitored by the Las Cumbres Observatory Global Telescope Network (LCOGT) during 2013-2014, and compare them with previously studied SNe having well-sampled light curves. We explore SN II diversity by searching for correlations between the slope of the linear light-curve decay after maximum light (historically used to divide SNe II into IIL and IIP) and other measured physical properties. While SNe IIL are found to be on average more luminous than SNe IIP, SNe IIL do not appear to synthesize more 56Ni than SNe IIP. Finally, optical nebular spectra obtained for several SNe in our sample are found to be consistent with models of red supergiant progenitors in the 12-16 Msun range. Consequently, SNe IIL appear not to account for the deficit of massive red supergiants as SN II progenitors.

Non-radial pulsations have been identified in a number of accreting white dwarfs in cataclysmic variables. These stars offer insight into the excitation of pulsation modes in atmospheres with mixed compositions of hydrogen, helium, and metals, and the response of these modes to changes in the white dwarf temperature. Among all pulsating cataclysmic variable white dwarfs, GW Librae stands out by having a well-established observational record of three independent pulsation modes that disappeared when the white dwarf temperature rose dramatically following its 2007 accretion outburst. Our analysis of HST ultraviolet spectroscopy taken in 2002, 2010 and 2011, showed that pulsations produce variations in the white dwarf effective temperature as predicted by theory. Additionally in May~2013, we obtained new HST/COS ultraviolet observations that displayed unexpected behaviour: besides showing variability at ~275s, which is close to the post-outburst pulsations detected with HST in 2010 and 2011, the white dwarf exhibits high-amplitude variability on a ~4.4h time-scale. We demonstrate that this variability is produced by an increase of the temperature of a region on white dwarf covering up to ~30 per cent of the visible white dwarf surface. We argue against a short-lived accretion episode as the explanation of such heating, and discuss this event in the context of non-radial pulsations on a rapidly rotating star

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.

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.

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)

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.

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.

Magnetic fields play a role in almost all stages of stellar evolution. Most low-mass stars, including the Sun, show surface fields that are generated by dynamo processes in their convective envelopes. Intermediate-mass stars do not have deep convective envelopes, although 10% exhibit strong surface fields that are presumed to be residuals from the stellar formation process. These stars do have convective cores that might produce internal magnetic fields, and these might even survive into later stages of stellar evolution, but information has been limited by our inability to measure the fields below the stellar surface. Here we use asteroseismology to study the occurrence of strong magnetic fields in the cores of low- and intermediate-mass stars. We have measured the strength of dipolar oscillation modes, which can be suppressed by a strong magnetic field in the core, in over 3600 red giant stars observed by Kepler. About 20% of our sample show mode suppression but this fraction is a strong function of mass. Strong core fields only occur in red giants above 1.1 solar masses (1.1Msun), and the occurrence rate is at least 60% for intermediate-mass stars (1.6--2.0Msun), indicating that powerful dynamos were very common in the convective cores of these stars.

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.

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}$.

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.

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, $\tau_0$, to $\tau_c=c/c_g$, where $c_g \approx$ 50 km/s is the isothermal sound speed in the gas alone. When $\tau_0 \gg \tau_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 $\tau_0 \ll \tau_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.

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.

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.

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.

Type-Iax supernovae (SN Iax) are stellar explosions that are spectroscopically similar to some type-Ia supernovae (SN Ia) at maximum light, except with lower ejecta velocities. They are also distinguished by lower luminosities. At late times, their spectroscopic properties diverge from other SN, but their composition (dominated by iron-group and intermediate-mass elements) suggests a physical connection to normal SN Ia. These are not rare; SN Iax occur at a rate between 5 and 30% of the normal SN Ia rate. The leading models for SN Iax are thermonuclear explosions of accreting carbon-oxygen white dwarfs that do not completely unbind the star, implying they are "less successful" cousins of normal SN Ia, where complete disruption is observed. Here we report the detection of the luminous, blue progenitor system of the type-Iax SN 2012Z in deep pre-explosion imaging. Its luminosity, colors, environment, and similarity to the progenitor of the Galactic helium nova V445 Puppis, suggest that SN 2012Z was the explosion of a white dwarf accreting from a helium-star companion. Observations in the next few years, after SN 2012Z has faded, could test this hypothesis, or alternatively show that this supernova was actually the explosive death of a massive star.

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.

We search for the Fe K${\alpha}$ 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${\alpha}$ 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${\alpha}$ 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${\alpha}$ 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}$.

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).

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.

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.

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).

We present optical and infrared monitoring data of SN 2012hn collected by the Public ESO Spectroscopic Survey for Transient Objects (PESSTO). We show that SN 2012hn has a faint peak magnitude (MR ~ -15.7) and shows no hydrogen and no clear evidence for helium in its spectral evolution. Instead, we detect prominent Ca II lines at all epochs, which relates this transient to previously described 'Ca-rich' or 'gap' transients. However, the photospheric spectra (from -3 to +32 d with respect to peak) of SN 2012hn show a series of absorption lines which are unique, and a red continuum that is likely intrinsic rather than due to extinction. Lines of Ti II and Cr II are visible. This may be a temperature effect, which could also explain the red photospheric colour. A nebular spectrum at +150d shows prominent CaII, OI, CI and possibly MgI lines which appear similar in strength to those displayed by core-collapse SNe. To add to the puzzle, SN 2012hn is located at a projected distance of 6 kpc from an E/S0 host and is not close to any obvious starforming region. Overall SN 2012hn resembles a group of faint H-poor SNe that have been discovered recently and for which a convincing and consistent physical explanation is still missing. They all appear to explode preferentially in remote locations offset from a massive host galaxy with deep limits on any dwarf host galaxies, favouring old progenitor systems. SN 2012hn adds heterogeneity to this sample of objects. We discuss potential explosion channels including He-shell detonations and double detonations of white dwarfs as well as peculiar core-collapse SNe.

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.

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.

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.

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]

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.

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.

Various lines of evidence suggest that very massive stars experience extreme mass-loss episodes shortly before they explode as a supernova. Interestingly, several models predict such pre-explosion outbursts. Establishing a causal connection between these mass-loss episodes and the final supernova explosion will provide a novel way to study pre-supernova massive-star evolution. Here we report on observations of a remarkable mass-loss event detected 40 days prior to the explosion of the Type IIn supernova SN 2010mc (PTF 10tel). Our photometric and spectroscopic data suggest that this event is a result of an energetic outburst, radiating at least 6x10^47 erg of energy, and releasing about 0.01 Solar mass at typical velocities of 2000 km/s. We show that the temporal proximity of the mass-loss outburst and the supernova explosion implies a causal connection between them. Moreover, we find that the outburst luminosity and velocity are consistent with the predictions of the wave-driven pulsation model and disfavor alternative suggestions.

Of the more than 150000 targets followed by the Kepler Mission, about 10% were selected as red giants. Due to their high scientific value, in particular for Galaxy population studies and stellar structure and evolution, their Kepler light curves were made public in late 2011. More than 13000 (over 85%) of these stars show intrinsic flux variability caused by solar-like oscillations making them ideal for large scale asteroseismic investigations. We automatically extracted individual frequencies and measured the period spacings of the dipole modes in nearly every red giant. These measurements naturally classify the stars into various populations, such as the red giant branch, the low-mass (M/Msol < 1.8) helium-core-burning red clump, and the higher-mass (M/Msol > 1.8) secondary clump. The period spacings also reveal that a large fraction of the stars show rotationally induced frequency splittings. This sample of stars will undoubtedly provide an extremely valuable source for studying the stellar population in the direction of the Kepler field, in particular when combined with complementary spectroscopic surveys.

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.

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.

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.

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.

There is a consensus that Type-Ia supernovae (SNe Ia) arise from the thermonuclear explosion of white dwarf stars that accrete matter from a binary companion. However, direct observation of SN Ia progenitors is lacking, and the precise nature of the binary companion remains uncertain. A temporal series of high-resolution optical spectra of the SN Ia PTF 11kx reveals a complex circumstellar environment that provides an unprecedentedly detailed view of the progenitor system. Multiple shells of circumsteller are detected and the SN ejecta are seen to interact with circumstellar material (CSM) starting 59 days after the explosion. These features are best described by a symbiotic nova progenitor, similar to RS Ophiuchi.

We construct a photometrically calibrated catalog of non-variable sources from the Palomar Transient Factory (PTF) observations. The first version of this catalog presented here, the PTF photometric catalog 1.0, contains calibrated R_PTF-filter magnitudes for about 21 million sources brighter than magnitude 19, over an area of about 11233 deg^2. The magnitudes are provided in the PTF photometric system, and the color of a source is required in order to convert these magnitudes into other magnitude systems. We estimate that the magnitudes in this catalog have typical accuracy of about 0.02 mag with respect to magnitudes from the Sloan Digital Sky Survey. The median repeatability of our catalog's magnitudes for stars between 15 and 16 mag, is about 0.01 mag, and it is better than 0.03 mag for 95% of the sources in this magnitude range. The main goal of this catalog is to provide reference magnitudes for photometric calibration of visible light observations. Subsequent versions of this catalog, which will be published incrementally online, will be extended to a larger sky area and will also include g_PTF-filter magnitudes, as well as variability and proper motion information.

(Abridged). The optical light curve of some SNe may be powered by the outward diffusion of the energy deposited by the explosion shock in optically thick circumstellar matter (CSM). Recently, it was shown that the radiation-mediated and -dominated shock in an optically thick wind must transform into a collisionless shock and can produce hard X-rays. The X-rays are expected to peak at late times, relative to maximum visible light. Here we report on a search, using Swift and Chandra, for X-ray emission from 28 SNe that belong to classes whose progenitors are suspected to be embedded in dense CSM (IIn/Ibn/SLSN-I). Two SNe in our sample have X-ray properties that are roughly consistent with the expectation for X-rays from a collisionless shock in optically thick CSM. Therefore, we suggest that their optical light curves are powered by shock breakout in CSM. We show that two other events were too X-ray bright during the SN maximum optical light to be explained by the shock breakout model. We conclude that the light curves of some, but not all, type-IIn/Ibn SNe are powered by shock breakout in CSM. For the rest of the SNe in our sample, including all the SLSN-I events, our X-ray limits are not deep enough and were typically obtained at too early times to conclude about their nature. We argue that the optical light curves of SNe, for which the X-ray emission peaks at late times, are likely powered by the diffusion of shock energy from a dense CSM. We comment about the possibility to detect some of these events in radio.

The nearby Type Ia supernova SN 2011fe in M101 (cz=241 km s^-1) provides a unique opportunity to study the early evolution of a "normal" Type Ia supernova, its compositional structure, and its elusive progenitor system. We present 18 high signal-to-noise spectra of SN 2011fe during its first month beginning 1.2 days post-explosion and with an average cadence of 1.8 days. This gives a clear picture of how various line-forming species are distributed within the outer layers of the ejecta, including that of unburned material (C+O). We follow the evolution of C II absorption features until they diminish near maximum light, showing overlapping regions of burned and unburned material between ejection velocities of 10,000 and 16,000 km s^-1. This supports the notion that incomplete burning, in addition to progenitor scenarios, is a relevant source of spectroscopic diversity among SNe Ia. The observed evolution of the highly Doppler-shifted O I 7774 absorption features detected within five days post-explosion indicate the presence of O I with expansion velocities from 11,500 to 21,000 km s^-1. The fact that some O I is present above C II suggests that SN 2011fe may have had an appreciable amount of unburned oxygen within the outer layers of the ejecta.

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)

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.

The Palomar Transient Factory (PTF) provides multiple epoch imaging for a large fraction of the sky. Here we describe the photometric calibration of the PTF data products that relates the PTF magnitudes to other mag systems. The calibration process utilizes SDSS r~16 mag point source objects as photometric standards. During photometric conditions, this allows us to solve for the extinction coefficients and color terms, and to estimate the camera illumination correction. This also enables the calibration of fields that are outside the SDSS footprint. We test the precision and repeatability of the PTF photometric calibration. Given that PTF is observing in a single filter each night, we define a PTF calibrated magnitude system for R-band and g-band. We show that, in this system, 59% (47%) of the photometrically calibrated PTF R-band (g-band) data achieve a photometric precision of 0.02-0.04 mag. Given the objects' color, the PTF magnitude system can be converted to other systems. A night-by-night comparison of the calibrated magnitudes of individual stars observed on multiple nights shows that they are consistent to a level of ~0.02 mag. Most of the data that were taken under non-photometric conditions can be calibrated relative to other epochs of the same sky footprint obtained during photometric conditions. We provide a guide describing the use of the PTF photometric calibration data products, as well as the transformations between the PTF magnitude system and the SDSS and Johnson-Cousins systems. (abridged).

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.

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.

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.

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.

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.

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.

Supernovae (SNe) are stellar explosions driven by gravitational or thermonuclear energy, observed as electromagnetic radiation emitted over weeks or more. In all known SNe, this radiation comes from internal energy deposited in the outflowing ejecta by either radioactive decay of freshly-synthesized elements (typically 56Ni), stored heat deposited by the explosion shock in the envelope of a supergiant star, or interaction between the SN debris and slowly-moving, hydrogen-rich circumstellar material. Here we report on a new class of luminous SNe whose observed properties cannot be explained by any of these known processes. These include four new SNe we have discovered, and two previously unexplained events (SN 2005ap; SCP 06F6) that we can now identify as members. These SNe are all ~10 times brighter than SNe Ia, do not show any trace of hydrogen, emit significant ultra-violet (UV) flux for extended periods of time, and have late-time decay rates which are inconsistent with radioactivity. Our data require that the observed radiation is emitted by hydrogen-free material distributed over a large radius (~10^15 cm) and expanding at high velocities (>10^4 km s^-1). These long-lived, UV-luminous events can be observed out to redshifts z>4 and offer an excellent opportunity to study star formation in, and the interstellar medium of, primitive distant galaxies.

PTF09dav is a peculiar subluminous type Ia supernova (SN) discovered by the Palomar Transient Factory (PTF). Spectroscopically, it appears superficially similar to the class of subluminous SN1991bg-like SNe, but it has several unusual features which make it stand out from this population. Its peak luminosity is fainter than any previously discovered SN1991bg-like SN Ia (M_B -15.5), but without the unusually red optical colors expected if the faint luminosity were due to extinction. The photospheric optical spectra have very unusual strong lines of Sc II and Mg I, with possible Sr II, together with stronger than average Ti II and low velocities of ~6000 km/s. The host galaxy of PTF09dav is ambiguous. The SN lies either on the extreme outskirts (~41kpc) of a spiral galaxy, or in an very faint (M_R>-12.8) dwarf galaxy, unlike other 1991bg-like SNe which are invariably associated with massive, old stellar populations. PTF09dav is also an outlier on the light-curve-width--luminosity and color--luminosity relations derived for other sub-luminous SNe Ia. The inferred 56Ni mass is small (0.019+/-0.003Msun), as is the estimated ejecta mass of 0.36Msun. Taken together, these properties make PTF09dav a remarkable event. We discuss various physical models that could explain PTF09dav. Helium shell detonation or deflagration on the surface of a CO white-dwarf can explain some of the features of PTF09dav, including the presence of Sc and the low photospheric velocities, but the observed Si and Mg are not predicted to be very abundant in these models. We conclude that no single model is currently capable of explaining all of the observed signatures of PTF09dav.

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.

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.

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.

Type-IIn supernovae (SNe), which are characterized by strong interaction of their ejecta with the surrounding circumstellar matter (CSM), provide a unique opportunity to study the mass-loss history of massive stars shortly before their explosive death. We present the discovery and follow-up observations of a Type IIn SN, PTF 09uj, detected by the Palomar Transient Factory (PTF). Serendipitous observations by GALEX at ultraviolet (UV) wavelengths detected the rise of the SN light curve prior to the PTF discovery. The UV light curve of the SN rose fast, with a time scale of a few days, to a UV absolute AB magnitude of about -19.5. Modeling our observations, we suggest that the fast rise of the UV light curve is due to the breakout of the SN shock through the dense CSM (n~10^10 cm^-3). Furthermore, we find that prior to the explosion the progenitor went through a phase of high mass-loss rate (~0.1 solar mass per year) that lasted for a few years. The decay rate of this SN was fast relative to that of other SNe IIn.

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.

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/).

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 $\alpha$ capturing material, accounting for thermal, mass defect, and Coulomb corrections on the stratification of material with the same zeroth order $\mu_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.

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.

Supernovae (SNe) are thought to arise from two different physical processes. The cores of massive, short-lived stars undergo gravitational core collapse and typically eject a few solar masses during their explosion. These are thought to appear as as type Ib/c and II SNe, and are associated with young stellar populations. A type Ia SN is thought to arise from the thermonuclear detonation of a white dwarf star composed mainly of carbon and oxygen, whose mass approaches the Chandrasekhar limit. Such SNe are observed in both young and old stellar environments. Here we report our discovery of the faint type Ib SN 2005E in the halo of the nearby isolated galaxy, NGC 1032. The lack of any trace of recent star formation near the SN location (Fig. 1), and the very low derived ejected mass (~0.3 M_sun), argue strongly against a core-collapse origin for this event. Spectroscopic observations and the derived nucleosynthetic output show that the SN ejecta have high velocities and are dominated by helium-burning products, indicating that SN 2005E was neither a subluminous nor a regular SN Ia (Fig. 2). We have therefore found a new type of stellar explosion, arising from a low-mass, old stellar system, likely involving a binary with a primary white dwarf and a helium-rich secondary. The SN ejecta contain more calcium than observed in any known type of SN and likely additional large amounts of radioactive 44Ti. Such SNe may thus help resolve fundamental physical puzzles, extending from the composition of the primitive solar system and that of the oldest stars, to the Galactic production of positrons.

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.

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.

Nonradial pulsations in the primary white dwarfs of cataclysmic variables can now potentially allow us to explore the stellar interior of these accretors using stellar seismology. In this context, we conducted a multi-site campaign on the accreting pulsator SDSS J161033.64-010223.3 (V386 Ser) using seven observatories located around the world in May 2007 over a duration of 11 days. We report the best fit periodicities here, which were also previously observed in 2004, suggesting their underlying stability. Although we did not uncover a sufficient number of independent pulsation modes for a unique seismological fit, our campaign revealed that the dominant pulsation mode at 609s is an evenly spaced triplet. The even nature of the triplet is suggestive of rotational splitting, implying an enigmatic rotation period of about 4.8 days. There are two viable alternatives assuming the triplet is real: either the period of 4.8 days is representative of the rotation period of the entire star with implications for the angular momentum evolution of these systems, or it is perhaps an indication of differential rotation with a fast rotating exterior and slow rotation deeper in the star. Investigating the possibility that a changing period could mimic a triplet suggests that this scenario is improbable, but not impossible. Using time-series spectra acquired in May 2009, we determine the orbital period of SDSS J161033.64-010223.3 to be 83.8 +/- 2.9 min. Three of the observed photometric frequencies from our May 2007 campaign appear to be linear combinations of the 609s pulsation mode with the first harmonic of the orbital period at 41.5min. This is the first discovery of a linear combination between nonradial pulsation and orbital motion for a variable white dwarf.

We have measured solar-like oscillations in red giants using time-series photometry from the first 34 days of science operations of the Kepler Mission. The light curves, obtained with 30-minute sampling, reveal clear oscillations in a large sample of G and K giants, extending in luminosity from the red clump down to the bottom of the giant branch. We confirm a strong correlation between the large separation of the oscillations (Delta nu) and the frequency of maximum power (nu_max). We focus on a sample of 50 low-luminosity stars (nu_max > 100 muHz, L <~ 30 L_sun) having high signal-to-noise ratios and showing the unambiguous signature of solar-like oscillations. These are H-shell-burning stars, whose oscillations should be valuable for testing models of stellar evolution and for constraining the star-formation rate in the local disk. We use a new technique to compare stars on a single echelle diagram by scaling their frequencies and find well-defined ridges corresponding to radial and non-radial oscillations, including clear evidence for modes with angular degree l=3. Measuring the small separation between l=0 and l=2 allows us to plot the so-called C-D diagram of delta nu_02 versus Delta nu. The small separation delta nu_01 of l=1 from the midpoint of adjacent l=0 modes is negative, contrary to the Sun and solar-type stars. The ridge for l=1 is notably broadened, which we attribute to mixed modes, confirming theoretical predictions for low-luminosity giants. Overall, the results demonstrate the tremendous potential of Kepler data for asteroseismology of red giants.

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.

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).

The Palomar Transient Factory (PTF) is a fully-automated, wide-field survey aimed at a systematic exploration of the optical transient sky. The transient survey is performed using a new 8.1 square degree camera installed on the 48-inch Samuel Oschin telescope at Palomar Observatory; colors and light curves for detected transients are obtained with the automated Palomar 60-inch telescope. PTF uses eighty percent of the 1.2-m and fifty percent of the 1.5-m telescope time. With an exposure of 60-s the survey reaches a depth of approximately 21.3 in g' and 20.6 in R (5 sigma, median seeing). Four major experiments are planned for the five-year project: 1) a 5-day cadence supernova search; 2) a rapid transient search with cadences between 90 seconds and 1 day; 3) a search for eclipsing binaries and transiting planets in Orion; and 4) a 3-pi sr deep H-alpha survey. PTF provides automatic, realtime transient classification and follow up, as well as a database including every source detected in each frame. This paper summarizes the PTF project, including several months of on-sky performance tests of the new survey camera, the observing plans and the data reduction strategy. We conclude by detailing the first 51 PTF optical transient detections, found in commissioning data.

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.

AM Canum Venaticorum (AM CVn) binaries consist of a degenerate helium donor and a helium, C/O, or O/Ne WD accretor, with accretion rates of Mdot = 1e-13 - 1e-5 Msol/yr. For accretion rates < 1e-6 Msol/yr, the accreted helium ignites unstably, resulting in a helium flash. As the donor mass and Mdot decrease, the ignition mass increases and eventually becomes larger than the donor mass, yielding a "last-flash" ignition mass of < 0.1 Msol. Bildsten et al. (2007) predicted that the largest outbursts of these systems will lead to dynamical burning and thermonuclear supernovae. In this paper, we study the evolution of the He-burning shells in more detail. We calculate maximum achievable temperatures as well as the minimum envelope masses that achieve dynamical burning conditions, finding that AM CVn systems with accretors > 0.8 Msol will undergo dynamical burning. Triple-alpha reactions during the hydrostatic evolution set a lower limit to the 12C mass fraction of 0.001 - 0.05 when dynamical burning occurs, but core dredge-up may yield 12C, 16O, and/or 20Ne mass fractions of ~ 0.1. Accreted 14N will likely remain 14N during the accretion and convective phases, but regardless of 14N's fate, the neutron-to-proton ratio at the beginning of convection is fixed until the onset of dynamical burning. During explosive burning, the 14N will undergo 14N(a,g)18F(a,p)21Ne, liberating a proton for the subsequent 12C(p,g)13N(a,p)16O reaction, which bypasses the relatively slow alpha-capture onto 12C. Future hydrodynamic simulations must include these isotopes, as the additional reactions will reduce the Zel'dovich-von Neumann-Doring (ZND) length, making the propagation of the detonation wave more likely.

Fundamental mysteries remain regarding the physics of Type Ia supernovae (SNIa) and their stellar progenitors. We argue here that important clues to these questions may emerge by the identification of those SNIa that occur in extragalactic globular clusters--stellar systems with well defined ages and metallicities. We estimate an all-sky rate of approximately 0.1 eta (D/100 Mpc)^3 per year for SNIa in globular clusters within a distance D, where eta is the rate enhancement per unit mass as a result of dynamical production channels that are inaccessible in the galactic field. If eta is approximately 2-10, as suggested by observations and theory, the combined efforts of accurate supernova astrometry and deep follow-up imaging should identify the > 1% of nearby (D < 100 Mpc) SNIa that occur in globular clusters.

We report the search for low-mass X-ray binaries in quiescence (qLMXBs) in the globular cluster NGC 6304 using XMM observations. We present the spectral analysis leading to the identification of three candidate qLMXBs within the field of this globular cluster (GC), each consistent with the X-ray spectral properties of previously identified qLMXBs in the field and in other globular clusters -- specifically, with a hydrogen atmosphere neutron star with radius between 5--20\km. One (source 4, with R=11.7^{+8.3}_{-0.4} (D/5.97 kpc) km and kT_eff=117^{+59}_{-44} eV) is located within one core radius (r_c) of the centre of NGC 6304. This candidate also presents a spectral power-law component contributing 49 per cent of the 0.5-10 keV flux. A second one (source 9 with R=15.3^{+11.2}_{-6.5} (D/5.97 kpc) km and kT_eff=100^{+24}_{-19} eV) is found well outside the optical core (at 32 r_c) but still within the tidal radius. From spatial coincidence, we identify a bright 2MASS infrared counterpart which, at the distance of NGC 6304, seems to be a post-asymptotic giant branch star. The third qLMXB (source 5 with R=23^{+38}_{-14} (D/5.97 kpc) km and kT_eff=70^{+28}_{-20} eV) is a low signal-to-noise candidate for which we also identify from spatial coincidence a bright 2MASS infrared counterpart, with 99.916 per cent confidence. Three qLMXBs from this GC is marginally consistent with that expected from the encounter rate of NGC 6304. We also report a low signal-to-noise source with an unusually hard photon index (\alpha=-2.0^{+1.2}_{-2.2}). Finally, we present an updated catalogue of the X-ray sources lying in the field of NGC 6304, and compare this with the previous catalogue compiled from ROSAT observations.

Measuring g-mode pulsations of isolated white dwarfs can reveal their interior properties to high precision. With a spectroscopic mass of ~0.51 M_{\odot} (log g = 7.82), the DAV white dwarf HS 1824+6000 is near the transition between carbon/oxygen core and helium core white dwarfs, motivating our photometric search for additional pulsations from the Palomar 60-inch telescope. We confirmed (with much greater precision) the three frequencies: 2.751190 +/- 0.000010 mHz (363.479 sec), 3.116709 +/- 0.000006 mHz (320.851 sec), 3.495113 +/- 0.000009 mHz (286.114 sec), previously found by B. Voss and collaborators, and found an additional pulsation at 4.443120 +/- 0.000012 mHz (225.067 sec). These observed frequencies are similar to those found in other ZZ Ceti white dwarfs of comparable mass (e.g. log g < 8). We hope that future observations of much lower mass ZZ Ceti stars (< 0.4 M_{\odot}) will reveal pulsational differences attributable to a hydrogen covered helium core.

We present multi-band photometric and optical spectroscopic observations of SN2007ax, the faintest and reddest Type Ia supernova (SNIa) yet observed. With M_B = -15.9 and (B-V)max = 1.2, this SN is over half a magnitude fainter at maximum light than any other SNIa. Similar to subluminous SN2005ke, SN2007ax also appears to show excess in UV emission at late time. Traditionally, Delta-m_15(B) has been used to parameterize the decline rate for SNeIa. However, the B-band transition from fast to slow decline occurs sooner than 15 days for faint SNeIa. Therefore we suggest that a more physically motivated parameter, the time of intersection of the two slopes, be used instead. Only by explaining the faintest (and the brightest) supernovae, we can thoroughly understand the physics of thermonuclear explosions. We suggest that future surveys should carefully design their cadence, depth, pointings and follow-up to find an unbiased sample of extremely faint members of this subclass of faint SNeIa.

Detection of gravitational waves from accreting neutron stars (NSs) in our galaxy, due to ellipticity or internal oscillation, would be a breakthrough in our understanding of compact objects and explain the absence of NSs rotating near the break-up limit. Direct detection, however, poses a formidable challenge. Using the current data available on the properties of the accreting NSs in Low Mass X-Ray Binaries (LMXBs), we quantify the detectability for the known accreting NSs, considering various emission scenarios and taking into account the negative impact of parameter uncertainty on the data analysis process. Only a few of the persistently bright NSs accreting at rates near the Eddington limit are detectable by Advanced LIGO if they are emitting gravitational waves at a rate matching the torque from accretion. A larger fraction of the known population is detectable if the spin and orbital parameters are known in advance, especially with the narrow-band Advanced LIGO. We identify the most promising targets, and list specific actions that would lead to significant improvements in detection probability. These include astronomical observations (especially for unknown orbital periods), improvements in data analysis algorithms and capabilities, and further detector development.

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.

We investigate the crystallization rate of a one-component plasma (OCP) in the context of classical nucleation theory. From our derivation of the free energy of an arbitrary distribution of solid clusters embedded in a liquid phase, we derive the steady-state nucleation rate of an OCP as a function of the Coulomb coupling parameter. Our result for the rate is in accord with recent molecular dynamics simulations, but it is greater than that of previous analytical estimates by many orders of magnitude. Further molecular dynamics simulations of the nucleation rate of a supercooled liquid OCP for several values of the coupling parameter would clarify the physics of this process.

We have discovered a partially eclipsing white dwarf, low-mass M dwarf binary (3 hour 54.41 second orbital period), SDSS J143547.87+373338.5, from May 2007 observations at the WIYN telescope. Here we present blue band photometry of three eclipses. Eclipse fitting gives main sequence solutions to the M dwarf companion of $M_S=0.15-0.35 M_{\odot}$ and $R_S=0.17-0.32 R_{\odot}$. Analysis of the SDSS spectrum constrains the M dwarf further to be of type M4-M6 with $M_S=0.11-0.20 M_{\odot}$. Once full radial velocity curves are measured, high precision determinations of the masses and radii of both components will be easily obtained without any knowledge of stellar structure or evolution. ZZ Ceti pulsations from the white dwarf were not found at our 4 mmag detection limit.

During accretion, a neutron star (NS) is spun up as angular momentum is transported through its liquid surface layers. We study the resulting differentially rotating profile, focusing on the impact this has for type I X-ray bursts. The viscous heating is found to be negligible, but turbulent mixing can be activated. Mixing has the greatest impact when the buoyancy at the compositional discontinuity between accreted matter and ashes is overcome. This occurs preferentially at high accretion rates or low spin frequencies and may depend on the ash composition from the previous burst. We then find two new regimes of burning. The first is ignition in a layer containing a mixture of heavier elements with recurrence times as short as ~5-30 minutes, similar to short recurrence time bursts. When mixing is sufficiently strong, a second regime is found where accreted helium mixes deep enough to burn stably, quenching X-ray bursts altogether. The carbon-rich material produced by stable helium burning would be important for triggering and fueling superbursts.

Prior to the incineration of a white dwarf (WD) that makes a Type Ia supernova (SN Ia), the star "simmers" for ~1000 years in a convecting, carbon burning region. We have found that weak interactions during this time increase the neutron excess by an amount that depends on the total quantity of carbon burned prior to the explosion. This contribution is in addition to the metallicity (Z) dependent neutronization through the 22Ne abundance (as studied by Timmes, Brown, & Truran). The main consequence is that we expect a floor to the level of neutronization that dominates over the metallicity contribution when Z/Z_\odot<2/3, and it can be important for even larger metallicities if substantial energy is lost to neutrinos via the convective Urca process. This would mask any correlations between SN Ia properties and galactic environments at low metallicities. In addition, we show that recent observations of the dependences of SNe Ia on galactic environments make it clear that metallicity alone cannot provide for the full observed diversity of events.

The strong degeneracy of the 12C ignition layer on an accreting neutron star results in a hydrodynamic thermonuclear runaway, in which the nuclear heating time becomes shorter than the local dynamical time. We model the resulting combustion wave during these superbursts as an upward propagating detonation. We solve the reactive fluid flow and show that the detonation propagates through the deepest layers of fuel and drives a shock wave that steepens as it travels upward into lower density material. The shock is sufficiently strong upon reaching the freshly accreted H/He layer that it triggers unstable 4He burning if the superburst occurs during the latter half of the regular Type I bursting cycle; this is likely the origin of the bright Type I precursor bursts observed at the onset of superbursts. The cooling of the outermost shock-heated layers produces a bright, ~0.1s, flash that precedes the Type I burst by a few seconds; this may be the origin of the spike seen at the burst onset in 4U 1820-30 and 4U 1636-54, the only two bursts observed with RXTE at high time resolution. The dominant products of the 12C detonation are 28Si, 32S, and 36Ar. Gupta et al. showed that a crust composed of such intermediate mass elements has a larger heat flux than one composed of iron-peak elements and helps bring the superburst ignition depth into better agreement with values inferred from observations.

The study of the observational properties of uncommonly long bursts from low luminosity sources with extended decay times up to several tens of minutes is important when investigating the transition from a hydrogen-rich bursting regime to a pure helium regime and from helium burning to carbon burning as predicted by current burst theories. IGR J17254-3257 is a recently discovered X-ray burster of which only two bursts have been recorded: an ordinary short type I X-ray burst, and a 15 min long burst. An upper limit to its distance is estimated to about 14.5 kpc. The broad-band spectrum of the persistent emission in the 0.3-100 keV energy band obtained using contemporaneous INTEGRAL and XMM-Newton data indicates a bolometric flux of 1.1x10^-10 erg/cm2/s corresponding, at the canonical distance of 8 kpc, to a luminosity about 8.4x10^35 erg/s between 0.1-100 keV, which translates to a mean accretion rate of about 7x10^-11 solar masses per year. The low X-ray persistent luminosity of IGR J17254-3257 seems to indicate the source may be in a state of low accretion rate usually associated with a hard spectrum in the X-ray range. The nuclear burning regime may be intermediate between pure He and mixed H/He burning. The long burst is the result of the accumulation of a thick He layer, while the short one is a prematurate H-triggered He burning burst at a slightly lower accretion rate.

Helium that accretes onto a Carbon/Oxygen white dwarf in the double white dwarf AM Canum Venaticorum (AM CVn) binaries undergoes unstable thermonuclear flashes when the orbital period is in the 3.5-25 minute range. At the shortest orbital periods (and highest accretion rates, Mdot > 10^-7 Msol/yr), the flashes are weak and likely lead to the Helium equivalent of classical nova outbursts. However, as the orbit widens and Mdot drops, the mass required for the unstable ignition increases, leading to progressively more violent flashes up to a final flash with Helium shell mass ~ 0.02-0.1 Msol. The high pressures of these last flashes allow the burning to produce the radioactive elements 48Cr, 52Fe, and 56Ni that power a faint (M_V in the range of -15 to -18) and rapidly rising (few days) thermonuclear supernova. Current galactic AM CVn space densities imply one such explosion every 5,000-15,000 years in 10^11 Msol of old stars (~ 2-6% of the Type Ia rate in E/SO galaxies). These ".Ia" supernovae (one-tenth as bright for one-tenth the time as a Type Ia supernovae) are excellent targets for deep (e.g. V=24) searches with nightly cadences, potentially yielding an all-sky rate of 1,000 per year.

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.

During accretion a neutron star (NS) is spun up as angular momentum is transported through its surface layers. We study the resulting differentially rotating profile, focusing on the impact this has for type I X-ray bursts. The predominant viscosity is likely provided by the Tayler-Spruit dynamo. The radial and azimuthal magnetic field components have strengths of ~10^5 G and ~10^10 G, respectively. This leads to nearly uniform rotation at the depths of interest for X-ray bursts. A remaining small shear transmits the accreted angular momentum inward to the NS interior. Though this shear gives little viscous heating, it can trigger turbulent mixing. Detailed simulations will be required to fully understand the consequences of mixing, but our models illustrate some general features. Mixing has the greatest impact when the buoyancy at the compositional discontinuity between accreted matter and ashes is overcome. This occurs at high accretion rates, at low spin frequencies, or may depend on the ashes from the previous burst. We then find two new regimes of burning. The first is ignition in a layer containing a mixture of heavier elements from the ashes. If ignition occurs at the base of the mixed layer, recurrence times as short as ~5-30 minutes are possible. This may explain the short recurrence time of some bursts, but incomplete burning is still needed to explain these bursts' energetics. When mixing is sufficiently strong, a second regime is found where accreted helium mixes deep enough to burn stably, quenching X-ray bursts. We speculate that the observed change in X-ray burst properties near one-tenth the Eddington accretion rate is from this mechanism. The carbon-rich material produced by stable helium burning would be important for triggering and fueling superbursts. (abridged)

We present results from our 830 ksec observation of the globular cluster 47 Tucanae with the Chandra X-ray Observatory's High Resolution Camera-S. We limit our analysis here to the 19 previously known, localized millisecond pulsars (MSPs) in the cluster. This work more than doubles the sample of X-ray-detected MSPs observed with sensitivity to rotational variability; it is also the first survey of a large group of radio-discovered MSPs for which no previous X-ray pulsations have been detected and is therefore an unbiased survey of the X-ray properties of radio-discovered MSPs. We find that only 47 Tuc D, O and R show significant pulsations at the >~ 4-sigma level, but there is statistical evidence for rotational variability in five additional MSPs. Furthermore, we constrain the pulsed magnetospheric emission of 7 more MSPs using Monte Carlo simulations. The result is that the majority of the 47 Tuc MSPs are characterized by low pulsed fractions, <~ 50%. In cases where larger pulsed fractions are measured the folded pulse profiles show relatively large duty cycles. When considered with previous spectroscopic studies, this suggests that the X-ray emission arises from the neutron star's heated polar caps, and in some cases, from intra-binary shocks, but generally not directly from the star's magnetosphere. We discuss the impact of these results on our understanding of high energy emission from MSPs.

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.

The Advanced Compton Telescope (ACT), the next major step in gamma-ray astronomy, will probe the fires where chemical elements are formed by enabling high-resolution spectroscopy of nuclear emission from supernova explosions. During the past two years, our collaboration has been undertaking a NASA mission concept study for ACT. This study was designed to (1) transform the key scientific objectives into specific instrument requirements, (2) to identify the most promising technologies to meet those requirements, and (3) to design a viable mission concept for this instrument. We present the results of this study, including scientific goals and expected performance, mission design, and technology recommendations.

(Abridged) The Cataclysmic Variable (CV) population harbors a diverse range of donor stars and accreting white dwarfs (WDs). A range of WD masses is expected, from low mass Helium core WDs, to massive WDs which have previously accreted at rates high enough for Hydrogen to burn steadily. Furthermore, a wide range of Helium enrichment is expected in the accreted material depending on the degree to which the donor star is evolved. We investigate the impact of this diversity on the range of effective temperatures ($T_{\rm eff}$) for which g-modes are unstable. The critical $T_{\rm eff}$ below which modes are unstable ("blue edge") depends on both surface gravity, $g$, and He abundance, $Y$. The Hydrogen/first Helium ionization instability strip is more sensitive to $g$ than $Y$. We find that (for solar composition envelopes), relative to a fiducial WD mass $0.6 M_\odot$, the blue edge for a $0.4 M_\odot$ He core WD shifts downward by $\approx 1000 {\rm K}$, while that for a massive $\approx 1.2 M_\odot$ WD shifts upward by $\approx 2000 {\rm K}$. The second Helium ionization instability strip exhibits strong dependences on both $g$ and $Y$. Surprisingly, increasing $Y$ by only 10% relative to solar creates an instability strip near $15,000 {\rm K}$. Hence CV's below the period gap with evolved donor stars of $Y\ga 0.4$ may have an "intermediate" instability strip well outside of the isolated DA and DB variables. This "intermediate" instability strip also occurs for low mass He WD with solar composition envelopes. The lack of pulsations in CV's with $T_{\rm eff}$ in the pure Hydrogen ZZ Ceti instability strip is also easily explained.

We consider the thermal structure and radii of strongly irradiated gas giant planets over a range in mass and irradiating flux. The cooling rate of the planet is sensitive to the surface boundary condition, which depends on the detailed manner in which starlight is absorbed and energy redistributed by fluid motion. We parametrize these effects by imposing an isothermal boundary condition $T \equiv T_{\rm deep}$ below the photosphere, and then constrain $T_{\rm deep}$ from the observed masses and radii. We compute the dependence of luminosity and core temperature on mass, $T_{\rm deep}$ and core entropy, finding that simple scalings apply over most of the relevant parameter space. These scalings yield analytic cooling models which exhibit power-law behavior in the observable age range $0.1-10 {\rm Gyr}$, and are confirmed by time-dependent cooling calculations. We compare our model to the radii of observed transiting planets, and derive constraints on $T_{\rm deep}$. Only HD 209458 has a sufficiently accurate radius measurement that $T_{\rm deep}$ is tightly constrained; the lower error bar on the radii for other planets is consistent with no irradiation. More accurate radius and age measurements will allow for a determination of the correlation of $T_{\rm deep}$ with the equilibrium temperature, informing us about both the greenhouse effect and day-night asymmetries.

We calculate the heating and cooling of the accreting white dwarf (WD) in the ultracompact AM Canum Venaticorum (AM CVn) binaries and show that the WD can contribute significantly to their optical and ultraviolet emission. We estimate the WD's effective temperature, Teff, using the optical continuum for a number of observed binaries, and show that it agrees well with our theoretical calculations. Driven by gravitational radiation losses, the time averaged accretion rate, <Mdot>, decreases monotonically with increasing Porb, covering six orders of magnitude. If the short period (Porb<10 min) systems accrete at a rate consistent with gravitational radiation via direct impact, we predict their unpulsed optical/UV light to be that of the Teff>50,000$ K accreting WD. At longer Porb we calculate the Teff and absolute visual magnitude, M_V, that the accreting WD will have during low accretion states, and find that the WD naturally crosses the pulsational instability strip. Discovery and study of pulsations could allow for the measurement of the accumulated helium mass on the accreting WD, as well as its rotation rate. Accretion heats the WD core, but for Porb>40 minutes, the WD's Teff is set by its cooling as <Mdot> plummets. For the two long period AM CVn binaries with measured parallaxes, GP Com and CE 315, we show that the optical broadband colors and intensity are that expected from a pure helium atmosphere WD. This confirms that the WD brightness sets the minimum light in wide AM CVn binaries, allowing for meaningful constraints on their population density from deep optical searches, both in the field and in Globular Clusters.

We solve for the evolution of the vertical extent of the convective region of a neutron star atmosphere during a Type I X-ray burst. The convective region is well-mixed with ashes of nuclear burning and its extent determines the rise time of the burst light curve. Using a full nuclear reaction network, we show that the maximum vertical extent of the convective region during photospheric radius expansion (RE) bursts can be sufficiently great that: (1) some ashes of burning are ejected by the radiation driven wind during the RE phase and, (2) some ashes of burning are exposed at the neutron star surface following the RE phase. We find that ashes with mass number A ~ 30 - 60 are mixed in with the ejected material. We calculate the expected column density of ejected and surface ashes in hydrogen-like states and determine the equivalent widths of the resulting photoionization edges from both the wind and neutron star surface. We find that these can exceed 100 eV and are potentially detectable. A detection would probe the nuclear burning processes and might enable a measurement of the neutron star gravitational redshift. In addition, we find that in bursts with pure helium burning layers, protons from (alpha, p) reactions cause a rapid onset of the 12C(p, gamma)13N(alpha, p)16O reaction sequence. The sequence bypasses the relatively slow 12C(alpha, gamma)16O reaction and leads to a sudden surge in energy production that is directly observable as a rapid (~ ms) increase in flux during burst rise.

The discovery of the first gravitationally redshifted spectral line from a neutron star (NS) by Cottam, Paerels and Mendez has triggered theoretical studies of the physics of atomic line formation in NS atmospheres. Chang, Bildsten and Wasserman showed that the hydrogenic Fe H$\alpha$ line formed above the photosphere of a bursting NS is intrinsically broad. We now include rotational broadening within general relativity and compare the resulting profile to that observed during Type I bursts from EXO 0748-676. We show that the fine structure splitting of the line precludes a meaningful constraint on the radius. Our fitting of the data show that the line forming Fe column is ${\rm log}_{10} (N_{\rm Fe, n=2}/{\rm cm^{-2}})=17.9_{-0.42}^{+0.27}$ and gravitational redshift $1+z =1.345_{-0.008}^{+0.005}$ with 95% confidence. We calculate the detectability of this spectral feature for a large range of spins and inclinations assuming that the emission comes from the entire surface. We find that at 300 (600) Hz only 10-20% (5-10%) of NSs would have spectral features as deep as that seen in EXO 0748-676.

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*10^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. Large area surveys should also turn up large numbers of quiescent X-ray binaries. (Modified)

We calculate the photon energy dependence of the pulsed amplitude of neutron star (NS) surface modes. Simple approximations demonstrate that it depends most strongly on the bursting NS surface temperature. This result compares well with full integrations that include Doppler shifts from rotation and general relativistic corrections to photon propagation. We show that the energy dependence of type I X-ray burst oscillations agrees with that of a surface mode, lending further support to the hypothesis that they originate from surface waves. The energy dependence of the pulsed emission is rather insensitive to the NS inclination, mass and radius, or type of mode, thus hindering constraints on these parameters. We also show that, for this energy-amplitude relation, the majority of the signal (relative to the noise) comes in the 2-25 keV band, so that the current burst oscillation searches with the Rossi X-Ray Timing Explorer are close to optimal. The critical test of the mode hypothesis for X-ray burst oscillations would be a measurement of the energy dependence of burst oscillations from an accreting millisecond pulsar.

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.

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

Accreting neutron stars (NSs) often show coherent modulations during type I X-ray bursts, called burst oscillations. We consider whether a nonradial mode can serve as an explanation for burst oscillations from those NSs which are not magnetic. We find that a surface wave in the shallow burning layer transitions into a crustal interface wave as the envelope cools, a new and previously uninvestigated phenomenon. The surface modulations decrease dramatically as the mode switches, explaining why burst oscillations often disappear before burst cooling ceases. When we include rotational modifications, we find mode frequencies and drifts consistent with those observed. The large NS spin ($\approx270-620 {\rm Hz}$) needed to make this match implies that accreting NSs are spinning at frequencies $\approx4 {\rm Hz}$ above the burst oscillation. Since the asymptotic frequency is set by the crustal interface wave, the observed late time frequency drifts are a probe of the composition and temperature of NS crusts. We compare our model with the observed drifts and persistent luminosities of X-ray burst sources, and find that NSs with a higher average accretion rate show smaller drifts, as we predict. Furthermore, the drift sizes are consistent with crusts composed of iron-like nuclei, as expected for the ashes of the He-rich bursts that are exhibited by these objects.

Motivated by the measurement of redshifted Fe H$\alpha$ lines during type I X-ray bursts on EXO 0748-676 (Cottam, Paerels & Mendez), we study the formation of atomic Fe lines above the photosphere of a bursting neutron star ($k_BT_{\rm eff} \approx 1-2 {\rm keV}$). We discuss the effects of Stark broadening, resonant scattering and NLTE (level population) on the formation of hydrogenic Fe H$\alpha$, Ly$\alpha$ and P$\alpha$ lines. From the observed equivalent width of the Fe H$\alpha$ line, we find an implied Fe column of $1-3 \times 10^{20} {\rm cm}^{-2}$, which is 3-10 times larger than the Fe column calculated from the accretion/spallation model of Bildsten, Chang & Paerels. We also estimate that the implied Fe column is about a factor of 2-3 larger than a uniform solar metallicity atmosphere. We discuss the effects of rotational broadening and find that the rotation rate of \EXO must be slow, as confirmed by the recent measurement of a 45 Hz burst oscillation by Villarreal & Strohmayer. We also show that the Fe Ly$\alpha$ EW $\approx$ 15-20 eV (redshifted 11-15 eV) and the P$\alpha$ EW $\approx$ 4-7 eV (redshifted 3-5 eV) when the H$\alpha$ EW is 10 eV (redshifted 8 eV). The Ly$\alpha$ line is rotationally broadened to a depth of $\approx 10%$, making it difficult to observe with {\it Chandra}. We also show that radiative levitation can likely support the Fe column needed to explain the line.

We follow the time dependent thermal evolution of a white dwarf (WD) undergoing sudden accretion in a dwarf nova outburst, using both simulations and analytic estimates. The post-outburst lightcurve clearly separates into early times when the WD flux is high, and late times when the flux is near the quiescent level. The break between these two regimes, occurring at a time of order the outburst duration, corresponds to a thermal diffusion wave reaching the base of the freshly accreted layer. Our principal result is that long after the outburst, the fractional flux perturbation about the quiescent flux decays as a power law with time (and {\it not} as an exponential). We use this result to construct a simple fitting formula that yields estimates for both the quiescent flux and the accreted column, i.e. the total accreted mass divided by WD surface area. The WD mass is not well constrained by the late time lightcurve alone, but it can be inferred if the accreted mass is known from observations. We compare our work with the well-studied outburst of WZ Sge, finding that the cooling is well described by our model, giving an effective temperature $T_{\rm eff}=14,500 {\rm K}$ and accreted column $\Delta y\approx10^6 {\rm g cm^{-2}}$, in agreement with the modeling of Godon et al. To reconcile this accreted column with the accreted mass inferred from the bolometric accretion luminosity, a large WD mass $\gtrsim1.1M_\odot$ is needed. Our power law result is a valuable tool for making quick estimates of the outburst properties. We show that fitting the late time lightcurve with this formula yields a predicted column within 20% of that estimated from our full numerical calculations.

Classical Novae (CNe) are the brightest manifestation of mass transfer onto a white dwarf in a cataclysmic variable (CV). As such, they are probes of the mass transfer rate, Mdot, and WD mass, Mwd, in these interacting binaries. Our calculations of the dependence of the CN ignition mass, Mign, on Mdot and Mwd yields the recurrence times of these explosions. We show that the observed CNe orbital period distribution is consistent with the interrupted magnetic braking evolutionary scenario, where at orbital periods Porb > 3 hr mass transfer is driven by angular momentum loss via a wind from the companion star and at Porb < 3 hr by gravitational radiation. About 50% of CNe occur in binaries accreting at Mdot ~= 10^{-9} Msun/yr with Porb = 3-4 hr, with the remaining 50% split evenly between Porb longer (higher Mdot) and shorter (lower Mdot) than this. This resolution of the relative contribution to the CN rate from different CVs tells us that 3(9)x10^5 CVs with WD mass 1.0(0.6)Msun are needed to produce one CN per year. Using the K-band specific CN rate measured in external galaxies, we find a CV birthrate of 2(4)x10^{-4}/yr per 10^{10}Lsun,K, very similar to the luminosity specific Type Ia supernova rate in elliptical galaxies. Likewise, we predict that there should be 60-180 CVs for every 10^6Lsun,K in an old stellar population, similar to the number of X-ray identified CVs in the globular cluster 47 Tuc, showing no overabundance relative to the field. Using a two-component steady state model of CV evolution we show that the fraction of CVs which are magnetic (22%) implies a birthrate of 8% relative to non-magnetic CVs, similar to the fraction of strongly magnetic field WDs. (abridged)

We present the first simultaneous, multi-wavelength observations of an L dwarf, the L3.5 candidate brown dwarf 2MASS J00361617+1821104, conducted with the Very Large Array, the Chandra X-ray Observatory, and the Kitt Peak 4-m telescope. We detect strongly variable and periodic radio emission (P=3 hr) with a fraction of about 60% circular polarization. No X-ray emission is detected to a limit of L_X/L_{bol}<2e-5, several hundred times below the saturation level observed in early M dwarfs. Similarly, we do not detect H-alpha emission to a limit of L_{H-alpha}/L_{bol}<2e-7, the deepest for any L dwarf observed to date. The ratio of radio to X-ray luminosity is at least four orders of magnitude in excess of that observed in a wide range of active stars (including M dwarfs) providing the first direct confirmation that late-M and L dwarfs violate the radio/X-ray correlation. The radio emission is due to gyrosynchrotron radiation in a large-scale magnetic field of about 175 G, which is maintained on timescales longer than three years. The detected 3-hour period may be due to (i) the orbital motion of a companion at a separation of about five stellar radii, similar to the configuration of RS CVn systems, (ii) an equatorial rotation velocity of about 37 km/s and an anchored, long-lived magnetic field, or (iii) periodic release of magnetic stresses in the form of weak flares. In the case of orbital motion, the magnetic activity may be induced by the companion, possibly explaining the unusual pattern of activity and the long-lived signal. We conclude that fully convective stars can maintain a large-scale and stable magnetic field, but the lack of X-ray and H-alpha emission indicates that the atmospheric conditions are markedly different than in early-type stars and even M dwarfs. [abridged]

We apply the Deloye & Bildsten (2003) isentropic models for donors in ultracompact low-mass X-ray binaries to the AM CVn population of ultracompact, interacting binaries. The mass-radius relations of these systems' donors in the mass range of interest ($M_2<0.1 \msun$) are not single-valued, but parameterized by the donor's specific entropy. This produces a range in the relationships between system observables, such as orbital period, $\Porb$, and mass transfer rate, $\Mdot$. For a reasonable range in donor specific entropy, $\Mdot$ can range over several orders of magnitude at fixed $\Porb$. We determine the unique relation between $\Mdot$ and $M_2$ in the AM CVn systems with known donor to accretor mass ratios, $q=M_2/M_1$. We use structural arguments, as well as each system's photometric behavior, to place limits on $\Mdot$ and $M_2$ in each. Most systems allow a factor of about 3 variation in $\Mdot$, although V803 Cen, if the current estimates of its $q$ are accurate, is an exception and must have $M_2 \approx 0.02 \msun$ and $\Mdot \approx 10^{-10} \msun$ yr$^{-1}$. Our donor models also constrain each donor's core temperature, $T_c$, range and correlate $T_c$ with $M_2$. We examine how variations in donor specific entropy across the white dwarf family \citep{nele01a} of AM CVn systems affects this population's current galactic distribution. Allowing for donors that are not fully degenerate produces a shift in systems towards longer $\Porb$ and higher $\Mdot$ increasing the parameter space in which these systems can be found. This shift increases the fraction of systems whose $\Porb$ is long enough that their gravity wave (GW) signal is obscured by the background of detached double white dwarf binaries that dominate the GW spectrum below a frequency $\approx 2$ mHz.

We describe recent theoretical work on the final stage of accretion when material passes from an accretion disk onto a white dwarf surface. Our calculations focus on understanding the latitudinal spreading and differentially rotating profile of this material, so we call it the ``spreading layer'' (SL) model. The SL typically extends to an angle of 0.01-0.1 radians, with respect to the equator. At low accretion rates ($\dot{M}\la10^{18} {\rm g s}^{-1}$) the amount of spreading is small, so that the dissipated energy is radiated back into the accretion disk. When the accretion rate is high, such as in dwarf novae, symbiotic binaries, and supersoft sources, the material spreads to higher latitudes to be more easily observed. The SL may contain shallow surface modes, and we propose that such waves could produce dwarf nova oscillations (DNOs). This hypothesis naturally explains many key properties of DNOs, including their frequency range, sinusoidal nature, sensitivity to accretion rate, and why they are only seen during outburst.

Dwarf nova outbursts often show coherent ($Q\sim10^4-10^6$) sinusoidal oscillations with the largest pulsed fraction in the extreme ultraviolet. Called dwarf nova oscillations (DNOs), they have periods of $P\approx3-40 {\rm s}$ and scale with luminosity as $P\propto L^{-\beta}$ with $\beta\approx0.1-0.2$. We propose that DNOs may be produced by nonradial oscillations in a thin hydrostatic layer of freshly accreted material, the ``spreading layer'' (SL), at the white dwarf (WD) equator. This would naturally explain a number of key properties of DNOs, including their frequency range, sinusoidal nature, sensitivity to accretion rate, and why they are only seen during outburst. In support of this hypothesis we construct a simple model that treats the SL as a cavity containing shallow surface waves, each with the same radial structure, but split into three different modes denoted by their azimuthal wavenumber, $m$. The $m=0$ latitudinally propagating mode best matches the periods and scalings associated with most DNOs, and DNOs with periods shorter than the WD Keplerian period are explained by the $m=-1$ prograde mode. We also predict a third set of oscillations, produced by the $m=1$ retrograde mode, and show its expected dependence on accretion rate.

We compute the rate of diffusive nuclear burning for hydrogen on the surface of a "magnetar" (Soft Gamma-Ray Repeater or Anomalous X-Ray Pulsar). We find that hydrogen at the photosphere will be burned on an extremely rapid timescale of hours to years, depending on composition of the underlying material. Improving on our previous studies, we explore the effect of a maximally thick "inert" helium layer, previously thought to slow down the burning rate. Since hydrogen diffuses faster in helium than through heavier elements, we find this helium buffer actually increases the burning rate for magnetars. We compute simple analytic scalings of the burning rate with temperature and magnetic field for a range of core temperature. We conclude that magnetar photospheres are very unlikely to contain hydrogen. This motivates theoretical work on heavy element atmospheres that are needed to measure effective temperature from the observed thermal emission and constrains models of AXPs that rely on magnetar cooling through thick light element envelopes.

Pulsation modes have recently been observed in a handful of white dwarf (WD) primaries of cataclysmic variables, allowing an interesting new probe into the structure of accreting WD's. We briefly discuss the seismology of these objects, how stellar properties may be inferred from the observed mode frequencies, and mode driving mechanisms. For one pulsator, GW Lib, we have shown that a WD mass $M=1.05 M_\odot$ and accreted envelope mass $M_{\rm env}=0.4\times 10^{-4} M_\odot$ give the best match to the observed pulsation periods. A first exploration of mode driving favors $T_{\rm eff} = 14 000 K$ and a massive WD, but more work is necessary.

The eigenfrequencies of nonradial oscillations are a powerful probe of a star's interior structure. This is especially true when there exist discontinuities such as at the neutron star (NS) ocean/crust boundary, as first noted by McDermott, Van Horn & Hansen. The interface mode associated with this boundary has subsequently been neglected in studies of stellar nonradial oscillations. We revisit this mode, investigating its properties both analytically and numerically for a simple NS envelope model. We find that it acts like a shallow surface ocean wave, but with a large radial displacement at the ocean/crust boundary due to flexing of the crust with shear modulus $\mu\ll P$, the pressure. This displacement lowers the mode's frequency by a factor of $\sim(\mu/P)^{1/2}\sim0.1$ in comparison to a shallow surface wave frequency on a hard surface. The interface mode may be excited on accreting or bursting NSs and future work on nonradial oscillations should consider this mode. Our work also implies an additional mode on massive and/or cold white dwarfs with crystalline cores, which may have a frequency between the f-mode and g-modes, an otherwise empty part of the frequency domain.

We present a first analysis of the g-mode oscillation spectrum for the white dwarf (WD) primary of GW Lib, a faint cataclysmic variable (CV). Stable periodicities have been observed from this WD for a number of years, but their interpretation as stellar pulsations has been hampered by a lack of theoretical models appropriate to an accreting WD. Using the results of Townsley and Bildsten, we construct accreting models for the observed effective temperature and approximate mass of the WD in GW Lib. We compute g-mode frequencies for a range of accreted layer masses, Macc, and long term accretion rates, <Mdot>. If we assume that the observed oscillations are from l=1 g-modes, then the observed periods are matched when M ~= 1.02 Msun, Macc ~= 0.31 x 10^-4 Msun and <Mdot> ~= 7.3 x 10^-11 Msun/yr. Much more sensitive observations are needed to discover more modes, after which we will be able to more accurately measure these parameters and constrain or measure the WD's rotation rate.

Recently we have made measurements of thermonuclear burst energetics and recurrence times which are unprecedented in their precision, largely thanks to the sensitivity of the Rossi X-ray Timing Explorer. In the "Clocked Burster", GS 1826-24, hydrogen burns during the burst via the rapid-proton (rp) process, which has received particular attention in recent years through theoretical and modelling studies. The burst energies and the measured variation of alpha (the ratio of persistent to burst flux) with accretion rate strongly suggests solar metallicity in the neutron star atmosphere, although this is not consistent with the corresponding variation of the recurrence time. Possible explanations include extra heating between the bursts, or a change in the fraction of the neutron star over which accretion takes place. I also present results from 4U 1746-37, which exhibits regular burst trains which are interrupted by "out of phase" bursts.

Chandra observations of distant elliptical galaxies have revealed large numbers of Low Mass X-ray Binaries (LMXBs) accreting at rates in excess of 10^{-8} solar masses per year. The majority of these LMXBs reside in globular clusters (GCs) and it has been suggested that many of the field LMXBs also originated in GCs. We show here that ultracompact binaries with orbital periods of 8-10 minutes and He or C/O donors of 0.06-0.08 solar masses naturally provide the observed accretion rates from gravitational radiation losses alone. Such systems are predicted to be formed in the dense GC environment, a hypothesis supported by the 11.4 minute binary 4U 1820-30, the brightest persistent LMXB in a Galactic GC. These binaries have short enough lifetimes (less then 3 Myr) while bright that we calculate their luminosity function under a steady-state approximation. This yields a luminosity function slope in agreement with that observed for luminosities in the range of 6E37 ergs/sec to 5E38 ergs/sec, encouraging us to use the observed numbers of LMXBs per GC mass to calculate the accumulated number of ultracompact binaries. For a constant birthrate over 8 Gyrs, the number of ultracompact binaries which have evolved through this bright phase is nearly 4000 in a 10 million solar mass GC, consistent with dynamical interaction calculations. Perhaps most importantly, if all ultracompacts become millisecond radio pulsars, then the observed normalization agrees with the inferred number of millisecond radio pulsars in 47 Tuc and Galactic GCs in general.

When a white dwarf (WD) is weakly magnetized and its accretion disk is thin, accreted material first reaches the WD's surface at its equator. This matter slows its orbit as it comes into co-rotation with the WD, dissipating kinetic energy into thermal energy and creating a hot band of freshly accreted material around the equator. Radiating in the extreme ultraviolet and soft X-rays, this material moves toward the pole as new material piles behind it, eventually becoming part of the WD once it has a temperature and rotational velocity comparable with the surface. We present a set of solutions which describe the properties of this ``spreading layer'' in the steady state limit based on the conservation equations derived by Inogamov & Sunyaev (1999) for accreting neutron stars. Our analysis and subsequent solutions show that the case of WDs is qualitatively different. We investigate example solutions of the spreading layer for a WD of mass $M=0.6M_\odot$ and radius $R=9\times10^{8}{\rm cm}$. These solutions show that the spreading layer typically extends to an angle of $\theta_{\rm SL}\approx0.01-0.1$ (with respect to the equator), depending on accretion rate and the magnitude of the viscosity. At low accretion rates, $\dot{M}\lesssim10^{18}{\rm g s}^{-1}$, the amount of spreading is negligible and most of the dissipated energy is radiated back into the accretion disk. When the accretion rate is high, such as in dwarf novae, symbiotic binaries, and supersoft sources, the material may spread to latitudes high enough to be directly visible above the accretion disk.

We extend our initial study of diffusive nuclear burning (DNB) for neutron stars (NSs) with Hydrogen atmospheres and an underlying layer of proton capturing nuclei. Our initial study showed that DNB can alter the photospheric abundance of Hydrogen on surprisingly short timescales ($10^{2-4}\yrs$). Significant composition evolution impacts the radiated thermal spectrum from the NS as well as its overall cooling rate. In this paper, we consider the case when the rate limiting step for the H consumption is diffusion to the burning layer, rather than the local nuclear timescale. This is relevant for NSs with surface temperatures in excess of $10^6 {\rm K}$, such as young ($<10^5$ yr) radio pulsars and accreting NSs in quiescence. When downward diffusion is the limiting rate in DNB, the rate of H consumption is suppressed by 1-2 orders of magnitude compared to a DNB estimate that assumes diffusive equilibrium. In order to apply our ongoing study to young neutron stars, we also include the important effects of strong magnetic fields ($B \sim 10^{12} {\rm G}$). In this initial study of magnetic modifications to DNB, we find that the H burning time is lengthened by 2-3 orders of magnitude for a $10^{12} {\rm G}$ field. However, even for NSs with dipole field strengths of $10^{12}$ G, we find that all of the H can be burned before the pulsar reaches an age of $\sim 10^5 \ {\rm yr}$, thus potentially revealing the underlying proton-capturing elements. Finally, we conclude by providing an overview of what can be learned about fallback and pulsar winds from measuring the surface composition of a young NS.

Material accreted onto a neutron star can stably burn in steady state only when the accretion rate is high (typically super-Eddington) or if a large flux from the neutron star crust permeates the outer atmosphere. For such situations we have analyzed the stability of nonradial oscillations, finding one unstable mode for pure helium accretion. This is a shallow surface wave which resides in the helium atmosphere above the heavier ashes of the ocean. It is excited by the increase in the nuclear reaction rate during the oscillations, and it grows on the timescale of a second. For a slowly rotating star, this mode has a frequency of approximately 20-30 Hz (for l=1), and we calculate the full spectrum that a rapidly rotating (>>30 Hz) neutron star would support. The short period X-ray binary 4U 1820--30 is accreting helium rich material and is the system most likely to show this unstable mode,especially when it is not exhibiting X-ray bursts. Our discovery of an unstable mode in a thermally stable atmosphere shows that nonradial perturbations have a different stability criterion than the spherically symmetric thermal perturbations that generate type I X-ray bursts.

We analyze 24 type I X-ray bursts from GS 1826-24 observed by the Rossi X-ray Timing Explorer between 1997 November and 2002 July. The bursts observed between 1997-98 were consistent with a stable recurrence time of 5.74 +/- 0.13 hr. The persistent intensity of GS 1826-24 increased by 36% between 1997-2000, by which time the burst interval had decreased to 4.10 +/- 0.08 hr. In 2002 July the recurrence time was shorter again, at 3.56 +/- 0.03 hr. The bursts within each epoch had remarkably identical lightcurves over the full approx. 150 s burst duration; both the initial decay timescale from the peak, and the burst fluence, increased slightly with the rise in persistent flux. The decrease in the burst recurrence time was proportional to Mdot^(-1.05+/-0.02) (where Mdot is assumed to be linearly proportional to the X-ray flux), so that the ratio alpha between the integrated persistent and burst fluxes was inversely correlated with Mdot. The average value of alpha was 41.7 +/- 1.6. Both the alpha value, and the long burst durations indicate that the hydrogen is burning during the burst via the rapid-proton (rp) process. The variation in alpha with Mdot implies that hydrogen is burning stably between bursts, requiring solar metallicity (Z ~ 0.02) in the accreted layer. We show that solar metallicity ignition models naturally reproduce the observed burst energies, but do not match the observed variations in recurrence time and burst fluence. Low metallicity models (Z ~ 0.001) reproduce the observed trends in recurrence time and fluence, but are ruled out by the variation in alpha. We discuss possible explanations, including extra heating between bursts, or that the fraction of the neutron star covered by the accreted fuel increases with Mdot.

White dwarfs experience a thermal renaissance when they receive mass from a stellar companion in a binary. For accretion rates < 10^-8 Msun/yr, the freshly accumulated hydrogen/helium envelope ignites in a thermally unstable manner that results in a classical novae (CN) outburst and ejection of material. We have undertaken a theoretical study of the impact of the accumulating envelope on the thermal state of the underlying white dwarf (WD). This has allowed us to find the equilibrium WD core temperatures (T_c), the classical nova ignition masses (M_ign) and the thermal luminosities for WDs accreting at rates of 10^-11 - 10^-8 Msun/yr. These accretion rates are most appropriate to WDs in cataclysmic variables (CVs) of P_orb <~ 7 hr, many of which accrete sporadically as dwarf novae. We have included ^3He in the accreted material at levels appropriate for CVs and find that it significantly modifies the CN ignition mass. We compare our results with several others from the CN literature and find that the inclusion of ^3He leads to lower M_ign for <Mdot> >~ 10^-10 Msun/yr, and that for <Mdot> below this the particular author's assumption concerning T_c, which we calculate consistently, is a determining factor. Initial comparisons of our CN ignition masses with measured ejected masses find reasonable agreement and point to ejection of material comparable to that accreted.

Our previous theoretical study of the impact of an accreting envelope on the thermal state of an underlying white dwarf (WD) has yielded equilibrium core temperatures, classical nova ignition masses and thermal luminosities for WDs accreting at time averaged rates of <Mdot> = 10^-11 - 10^-8 Msun/yr. These <Mdot>'s are appropriate to WDs in cataclysmic variables (CVs) of P_orb <~ 7 hr, many of which accrete sporadically as Dwarf Novae. Approximately thirty nonmagnetic Dwarf Novae have been observed in quiescence, when the accretion rate is low enough for spectral detection of the WD photosphere, and a measurement of T_eff. We use our theoretical work to translate the measured T_eff's into local time-averaged accretion rates, confirming the factor of ten drop in <Mdot> predicted for CV's as they transit the period gap. For DN below the period gap, we show that if <Mdot> is that given by gravitational radiation losses alone, then the WD masses are > 0.8 Msun. An alternative conclusion is that the masses are closer to 0.6 Msun and <Mdot> is 3-4 times larger than that expected from gravitational radiation losses. In either case, it is very plausible that a subset of CVs with P_orb < 2 hours will have T_eff's low enough for them to become non-radial pulsators, as discovered by van Zyl and collaborators in GW Lib.

We discuss the mass-radius (M-R) relations for low-mass (M<0.1 Msun) white dwarfs (WDs) of arbitrary degeneracy and evolved (He, C, O) composition. We do so with both a simple analytical model and models calculated by integration of hydrostatic balance using a modern equation of state valid for fully ionized plasmas. The M-R plane is divided into three regions where either Coulomb physics, degenerate electrons or a classical gas dominate the WD structure. For a given M and central temperature, T_c, the M-R relation has two branches differentiated by the model's entropy content. We present the M-R relations for a sequence of constant entropy WDs of arbitrary degeneracy parameterized by M and T_c for pure He, C, and O. We discuss the applications of these models to the recently discovered accreting millisecond pulsars. We show the relationship between the orbital inclination for these binaries and the donor's composition and T_c. In particular we find from orbital inclination constraints that the probability XTE J1807-294 can accommodate a He donor is approximately 15% while for XTE J0929-304, it is approximately 35%. We argue that if the donors in ultracompact systems evolve adiabatically, there should be 60-160 more systems at orbital periods of 40 min than at orbital periods of 10 min, depending on the donor's composition.

The gravitational redshift measured by Cottam, Paerels and Mendez for the neutron star (NS) in the low-mass X-ray binary EXO 0748-676 depends on the identification of an absorption line during a type I burst as the H$\alpha$ line from hydrogenic Fe. We show that Fe is present above the photosphere as long as $\dot M>4\times 10^{-13}M_\odot {\rm yr^{-1}}$ during the burst. In this limit, the total Fe column is $N_{\rm Fe}\approx 3\times 10^{19}{\rm cm^{-2}}$ for incident material of solar abundances and only depends on the nuclear physics of the proton spallation. The Fe destruction creates many heavy elements with $Z<26$ which may imprint photo-ionization edges on the NS spectra during a radius expansion event or in a burst cooling tail. Detecting these features in concert with those from Fe would confirm a redshift measurement. We also begin to address the radiative transfer problem, and find that a concentrated Fe layer with $kT=1.2-1.4 {\rm keV}$ and column $N_{\rm Fe}= 7-20 \times 10^{20} {\rm cm}^{-2}$ (depending on the line depth) above the hotter continuum photosphere is required to create the H$\alpha$ line of the observed strength. This estimate must be refined by considerations of non-LTE effects as well as resonant line transport. Until these are carried out, we cannot say whether the Fe column from accretion and spallation is in conflict with the observations. We also show that hydrogenic Fe might remain in the photosphere due to radiative levitation from the high burst flux.

Since the advent of powerful new X-ray observatories, NASA's Rossi X-ray Timing Explorer (RXTE), the Italian - Dutch BeppoSAX mission, XMM-Newton and Chandra, a number of entirely new phenomena associated with thermonuclear burning on neutron stars have been discovered. These include: (i) the discovery of millisecond (300 - 600 Hz) oscillations during bursts, so called ``burst oscillations'', (ii) a new regime of nuclear burning on neutron stars which manifests itself through the generation of hours long flares about once a decade, now referred to as ``superbursts'',(iii) discoveries of bursts from low accretion rate neutron stars, and (iv) new evidence for discrete spectral features from bursting neutron stars. In this article we review our current understanding of thermonuclear bursts on neutron stars, with a focus on these new phenomena.

Recent theoretical work has made it plausible for neutron stars (NSs) to lose angular momentum via gravitational radiation on long timescales (around Myrs) while actively accreting. The gravitational waves (GWs) can either be emitted via the excitation of r-modes or from a deformed crust. GW emission can thus intervene to slow-down or halt the otherwise relentless spin-up from accretion. Prior to this theoretical work (and the measurements of NS rotation rates in LMXBs) the community was rather confident that an accreting NS would be spun-up to rotation rates near breakup, motivating searches for sub-millisecond objects. After only briefly describing the physics of the GW processes, I argue that the limiting spin frequency might be appreciably lower than the breakup frequency. Millisecond radio pulsar observers would likely discover the impact of GW emission as a dropoff in the number of pulsars beyond 600 Hz, and I show here that the millisecond pulsar inventory in 47 Tuc might already exhibit such a cutoff. These theoretical ideas will be tested by GW searches with ground-based interferometers, such as the advanced LIGO instrument proposed for operation by 2008.

We present a new mode of hydrogen burning on neutron stars (NSs) called diffusive nuclear burning (DNB). In DNB, the burning occurs in the exponentially suppressed tail of hydrogen that extends to the hotter regions of the envelope where protons are readily captured. Diffusive nuclear burning changes the compositional structure of the envelope on timescales $\sim 10^{2-4} {\rm yrs}$, much shorter than otherwise expected. This mechanism is applicable to the physics of young pulsars, millisecond radio pulsars (MSPs) and quiescent low mass X-ray binaries (LMXBs).

We calculate the rate of hydrogen burning for neutron stars (NSs) with hydrogen atmospheres and an underlying reservoir of nuclei capable of proton capture. This burning occurs in the exponentially suppressed diffusive tail of H that extends to the hotter depths of the envelope where protons are rapidly captured. This process, which we call diffusive nuclear burning (DNB), can change the H abundance at the NS photosphere on timescales as short as $10^{2-4}$ years. In the absence of diffusion, the hydrogen at the photosphere (where $T\approx 10^6 {\rm K}$ and $\rho\sim 0.1 {\rm g cm^{-2}}$) would last for far longer than a Hubble time. Our work impacts the understanding of the evolution of surface abundances of isolated NSs, which is important to their thermal spectrum and their effective temperature-core temperature relation. In this paper, we calculate the rate of H burning when the overall consumption rate is controlled by the nuclear timescales, rather than diffusion timescales. The immediate application is for H burning on millisecond radio pulsars and in quiescence for the accreting NS Cen X-4. We will apply this work to young radio pulsars and magnetars once we have incorporated the effects of strong $B>10^{12} {\rm G}$ magnetic fields.

We assess the impact of the trace element Ne-22 on the cooling and seismology of a liquid C/O white dwarf (WD). Due to this element's neutron excess, it sinks towards the interior as the liquid WD cools. The gravitational energy released slows the WD's cooling by 0.5-1.6 Gyr. In addition the Ne-22 abundance gradient changes the periods of the high radial order g-modes at the 1% level.

We explain the physics of compressional heating of the deep interior of an accreting white dwarf (WD) at accretion rates low enough so that the accumulated hydrogen burns unstably and initiates a classical nova (CN). In this limit, the WD core temperature (T_c) reaches an equilibrium value (T_c,eq) after accreting an amount of mass much less than the WD's mass. Once this equilibrium is reached, the compressional heating from within the envelope exits the surface. This equilibrium yields useful relations between the WD surface temperature, accretion rate and mass that can be employed to measure accretion rates from observed WD effective temperatures, thus testing binary evolution models for cataclysmic variables.

Our recent theoretical work (Townsley and Bildsten 2002) on the thermal state of white dwarfs (WDs) in low mass transfer rate binaries allows us to predict the broadband colors of the binary from those of the WD and companion when the disk is dim. The results based on standard CV evolution are presented h ere. These will aid the discovery of such objects in field surveys and proper-motion selected globular cluster surveys with HST; especially for the largely unexplored post period minim um Cataclysmic Variables (CVs) with the lowest accretion rates and degenerate companions. We have also calculated the fraction of time that the WD resides in the ZZ Ceti instability strip thus clarifying that we expect many accreting WDs to exhibit non-radial oscillations. The study of these will provide new insights into the rotational and thermal structure of an actively accreting WD.

The discovery of two accreting millisecond X-ray pulsars in binaries with 43 minute orbital periods allows for a new probe of the donor's structure. For XTE J1751-305, only a hot white dwarf (WD) can fill the Roche Lobe. A cold He WD is a possible solution for XTE J0929-314, though I will show that evolutionary arguments make a hot WD more likely. In addition to being larger than the T=0 models, these finite entropy, low-mass (<0.03 solar masses) WDs have a minimum mass for a fixed core temperature. If they remain hot as they lose mass and expand, they can ``evaporate'' to leave an isolated millisecond radio pulsar. They also adiabatically expand upon mass loss at a rate faster than the growth of the Roche radius if the angular momentum deposited in the disk is not returned to the donor. If the timescale of the resulting runaway mass transfer is shorter than the viscous timescale in the outer disk, then the mass transfer instability of Ruderman and Shaham for He WDs would be realized. However, my estimates of these timescales still makes the instability unlikely for adiabatic responses. I close by noting the possible impact of finite T WDs on our understanding of AM CVn binaries.

We assess the impact of the trace element ^{22}Ne on the cooling and seismology of a liquid C/O white dwarf (WD). Due to this elements' neutron excess, it sinks towards the interior as the liquid WD cools. The subsequent gravitational energy released slows the cooling of the WD by 0.25--1.6 Gyrs by the time it has completely crystallized, depending on the WD mass and the adopted sedimentation rate. The effects will make massive WDs or those in metal rich clusters (such as NGC 6791) appear younger than their true age. Our diffusion calculations show that the ^{22}Ne mass fraction in the crystallized core actually increases outwards. The stability of this configuration has not yet been determined. In the liquid state, the settled ^{22}Ne enhances the internal buoyancy of the interior and changes the periods of the high radial order g-modes by approximately 1%. Though a small adjustment, this level of change far exceeds the accuracy of the period measurements. A full assessment and comparison of mode frequencies for specific WDs should help constrain the still uncertain ^{22}Ne diffusion coefficient for the liquid interior.

Chandra observations of elliptical galaxies have revealed large numbers of Low Mass X-ray Binaries (LMXBs) accreting at rates above 10^-9 solar masses per year. One scenario which generates this transfer rate from an old stellar population is nuclear driven mass transfer onto a neutron star or black hole from a Roche lobe filling red giant. However, in our Galaxy, most of these systems accrete sporadically as transients due to a thermal instability in the accretion disk. Using the common criterion for disk instability (including irradiation), we find that this mode of mass transfer leads to transient accretion for at least 75% of the binary's life. Repeated Chandra observations of elliptical galaxies should reveal this population. The recurrence times might exceed 1000 years (depending on the orbital period at the onset of mass transfer), and outbursts might last for 1-100 years. Mass transfering binaries can also be formed in old populations via interactions in dense stellar environments, such as globular clusters. These tend to have shorter orbital periods and are more likely stable accretors, making them apparently a large fraction of the elliptical's LMXB population.

We obtained four Chandra/ACIS-S observations beginning two weeks after the end of the November 2000 outburst of the neutron star (NS) transient Aql X-1. Over the five month span in quiescence, the X-ray spectra are consistent with thermal emission from a NS with a pure hydrogen photosphere and R_{\infty}=15.9+{0.8}-{2.9} (d/5 kpc) km at the optically implied X-ray column density. We also detect a hard power-law tail during two of the four observations. The intensity of Aql X-1 first decreased by 50+/-4% over three months, then increased by 35+/-5% in one month, and then remained constant (<6% change) over the last month. These variations in the first two observations cannot be explained by a change in the power-law spectral component, nor in the X-ray column density. Presuming that R_{\infty} is not variable and a pure hydrogen atmosphere, the long-term changes can only be explained by variations in the NS effective temperature, from kT_{eff, \infty}=130+3-5 eV, down to 113+3-4 eV, finally increasing to 118+9-4 eV for the final two observations. During one of these observations, we observe two phenomena which were previously suggested as indicators of quiescent accretion onto the NS: short-timescale (<1e4 sec) variability (at 32+8-6% rms), and a possible absorption feature near 0.5 keV. The possible absorption feature can potentially be explained as due to a time-variable response in the ACIS detector. Even so, such a feature has not been detected previously from a NS, and if confirmed and identified, can be exploited for simultaneous measurements of the photospheric redshift and NS radius.

The composition of the outer 100 m of a neutron star sets the heat flux that flows outwards from the core. For an accreting neutron star in an X-ray transient, the thermal quiescent flux depends sensitively on the amount of hydrogen and helium remaining on the surface after an accretion outburst and on the composition of the underlying ashes of previous H/He burning. Because H/He has a higher thermal conductivity, a larger mass of H/He implies a shallower thermal gradient through the low density envelope and hence a higher effective temperature for a given core temperature. The mass of residual H and He varies from outburst to outburst, so the thermal quiescent flux is variable even though the core temperature is constant for timescales < 10 000 yr. Heavy elements settle from a H/He envelope in a few hours; we therefore model the quiescent envelope as two distinct layers, H/He over heavier elements, and treat the mass of H/He as a free parameter. We find that the emergent thermal quiescent flux can vary by a factor of 2 to 3 between different quiescent epochs. The variation is more pronounced at lower interior temperatures, making systems with low quiescent luminosities and frequent outbursts, such as SAX J1808.4-3658, ideal candidates from which to observe this effect. We compute, for different ash compositions, the interior temperatures of Cen X-4, Aql X-1, and SAX J1808.4-3658. In the case of Aql X-1, the inferred high interior temperature suggests that neutrino cooling contributes to the neutron star's thermal balance.

In a photometric variability survey of 131 stars with B < 14 mag, we have found four new Delta Scuti stars. We were sensitive to oscillation amplitudes as low as a few mmag. The detection rate of short-period (P < 0.1 d) pulsating variable stars, which may be relevant for planned large variability surveys such as GAIA, POI, and even the LSST, was therefore 3%. All four new variable stars have low pulsation amplitude (tens of mmag), and one has a very short period (0.0198 d). This star is one of the fastest known Delta Scuti pulsators. The short period of this variable star makes it observationally tractable, and it may therefore be a particularly good candidate for asteroseismological studies. All four new variable stars will add to the cadre of low-amplitude and relatively short-period Delta Scuti stars that are potentially useful for learning about the structure of stars on or near the main sequence, slightly more massive than the Sun.

Current evolutionary models imply that most Cataclysmic Variables (CVs) have P_orb<2 hours and are Dwarf Nova (DN) systems that are quiescent most of the time. Observations of nearby quiescent DN find that the UV spectrum is dominated by the hot white dwarf (WD), indicating that it provides a significant fraction of the optical light in addition to the quiescent disk and main sequence companion. Hence, identifying a faint, quiescent CV in either the field or a globular cluster (GC) from broadband colors depends on our ability to predict the WD contribution in quiescence. We are undertaking a theoretical study of the compressional heating of WDs, extending down to very low time averaged accretion rates, <Mdot>~10^{-11}M_sun/yr, which allows us to self-consistently find the T_eff} of the WD. We demonstrate here that most of the compressional heating occurs in the freshly accreted envelope and that the WD core temperature reaches a fixed value on a timescale less than typical evolutionary times. Since nuclear burning is unstable at these <Mdot>'s, we have incorporated the recurrent heating and cooling of the WD core throughout the classical novae limit cycle in order to find the T_eff-<Mdot> relations. Comparing to observations of field DN confirms the <Mdot>-P_orb relation of disrupted magnetic braking. We also predict broad-band colors of a quiescent CV as a function of <Mdot> and companion mass and show that this leads to the identification of what may be many CVs in deep HST images of GCs.

In recent years several Dwarf Novae (DN) systems have been observed in quiescence, when the accretion rate is low and the WD photosphere can be directly detected. The WDs are observed to cool after the DN outburst from high effective temperatures to lower effective temperatures (T_eff) thought to be indicative of the thermal state of the deep interior of the WD. Sion has argued that the most likely energy source for this quiescent luminosity is the gravitational compression of the WD interior, which rejuvenates an otherwise cold WD into a much hotter state. We are undertaking a theoretical study of the compressional heating of WD's, extending down to the very low time averaged accretion rates, <Mdot>~10^{-11}M_sun/yr, applicable to the post-turnaround CV's (the ``TOADS''). Nuclear burning is unstable at these <Mdot>'s, so we have incorporated the recurrent heating and cooling of the WD throughout the classical novae limit cycle. In addition to self-consistently finding the range of T_eff as a function of <Mdot> during the cycle, we also self-consistently find the ignition masses. Comparing these theoretical masses to the observed ejected masses will tell us whether the WD mass in CV's is secularly increasing or decreasing. We close by comparing our results to the accumulated observations of quiescent DN and making predictions for the colors of low <Mdot> CV's in quiescence that are applicable to searches for faint CVs in the field and galactic globular clusters.

(Abridged). The type-I X-ray bursting low mass X-ray binary KS 1731-260 was recently detected for the first time in quiescence by Wijnands et al., following an approximately 13 yr outburst which ended in Feb 2001. Unlike all other known transient neutron stars, the duration of this recent outburst is as long as the thermal diffusion time of the crust. The large amount of heat deposited by reactions in the crust will have heated the crust to temperatures much higher than the equilibrium core temperature. As a result, the thermal luminosity currently observed from the neutron star is dominated not by the core, but by the crust. Moreover, the level and the time evolution of quiescent luminosity is determined mostly by the amount of heat deposited in the crust during the most recent outburst. Using estimates of the outburst mass accretion rate, our calculations of the quiescent flux immediately following the end of the outburst agree with the observed quiescent flux to within a factor of a few. We present simulations of the evolution of the quiescent lightcurve for different scenarios of the crust microphysics, and demonstrate that monitoring observations (with currently flying instruments) spanning from 1--30 yr can measure the crust cooling timescale and the total amount of heat stored in the crust. These quantities have not been directly measured for any neutron star.

We present calculations of the spin-down of a neutron star atmosphere due to hydrostatic expansion during a Type I X-ray burst. We show that (i) Cumming and Bildsten overestimated the spin-down of rigidly-rotating atmospheres by a factor of two, and (ii) general relativity has a small (5-10%) effect on the angular momentum conservation law. We rescale our results to different neutron star masses, rotation rates and equations of state, and present some detailed rotational profiles. Comparing with recent observations of large frequency shifts in MXB 1658-298 and 4U 1916-053, we find that the spin-down expected if the atmosphere rotates rigidly is a factor of two to three less than the observed values. If differential rotation is allowed to persist, we find that the upper layers of the atmosphere spin down by an amount comparable to the observed values; however, there is no compelling reason to expect the observed spin frequency to be that of only the outermost layers. We conclude that hydrostatic expansion and angular momentum conservation alone cannot account for the largest frequency shifts observed during Type I bursts.

We show that burning of a small mass fraction of carbon in a neutron star ocean is thermally unstable at low accumulated masses when the ocean contains heavy ashes from the hydrogen burning rapid proton (rp) process. The key to early unstable ignition is the low thermal conductivity of a heavy element ocean. The instability requires accretion rates in excess of one-tenth the Eddington limit when the carbon mass fraction is 0.1 or less. The unstable flashes release 10^{42} to 10^{43} ergs over hours to days, and are likely the cause of the recently discovered large Type I X-ray bursts (so-called ``superbursts'') from six Galactic low mass X-ray binaries. In addition to explaining the energetics, recurrence times, and durations of the superbursts, these mixed carbon/heavy element flashes have an accretion rate dependence of unstable burning similar to that observed. Though the instability is present at accretion rates near Eddington, there is less contrast with the accretion luminosity there, explaining why most detections are made at accretion rates between 0.1 and 0.3 Eddington. Future comparisons of time dependent calculations with observations will provide new insights into the rp process.

Using the Chandra/ACIS-I detector, we have identified an X-ray source (CXOU 132619.7-472910.8) in the globular cluster NGC 5139 with a thermal spectrum identical to that observed from transiently accreting neutron stars in quiescence. The absence of intensity variability on timescales as short as 4 seconds (< 25% rms variability) and as long as 5 years (<50% variability) supports the identification of this source as a neutron star, most likely maintained at a high effective temperature (approximately 1e6 K) by transient accretion from a binary companion. The ability to spectrally identify quiescent neutron stars in globular clusters (where the distance and interstellar column densities are known) opens up new opportunities for precision neutron star radius measurements.

We report on the quiescent spectrum measured with Chandra/ACIS-S of the transient, type-I X-ray bursting neutron star Aql X-1, immediately following an accretion outburst. The neutron star radius, assuming a pure hydrogen atmosphere and hard power-law spectrum, is $R_\infty$=13.4{+5}{-4} (d/5 \kpc) km. Based on the historical outburst record of RXTE/ASM, the quiescent luminosity is consistent with that predicted by Brown, Bildsten and Rutledge from deep crustal heating, lending support to this theory for providing a minimum quiescent luminosity of transient neutron stars. While not required by the data, the hard power-law component can account for 18+/-8% of the 0.5-10 keV thermal flux. Short-timescale intensity variability during this observation is less than 15% rms (3 sigma; 0.0001-1 Hz, 0.2-8 keV). Comparison between the Chandra spectrum and three X-ray spectral observations made between Oct 1992 and Oct 1996 find all spectra consistent with a pure H atmosphere, but with temperatures ranging from 145--168 eV, spanning a factor of 1.87+/-0.21 in observed flux. The source of variability in the quiescent luminosity on long timescales (greater than years) remains a puzzle. If from accretion, then it remains to be explained why the quiescent accretion rate provides a luminosity so nearly equal to that from deep crustal heating.

We present observations of the Rapid Burster (RB, also known as MXB 1730-335) using the Chandra High Energy Transmission Grating Spectrometer. The average interval between type II (accretion) bursts was about 40 s. There was one type I (thermonuclear flash) burst and about 20 "mini-bursts" which are probably type II bursts whose peak flux is 10-40% of the average peak flux of the other type II bursts. The time averaged spectra of the type II bursts are well fit by a blackbody with a temperature of kT = 1.6 keV, a radius of 8.9 km for a distance of 8.6 kpc, and an interstellar column density of 1.7e22 per sq. cm. No narrow emission or absorption lines were clearly detected. The 3 sigma upper limits to the equivalent widths of any features are < 10 eV in the 1.1-7.0 keV band and as small as 1.5 eV near 1.7 keV. We suggest that Comptonization destroys absorption features such as the resonance line of Fe XXVI.

BeppoSAX has revealed a population of faint neutron star X-ray transients in the Galactic bulge. King conjectured that these neutron stars are accreting from brown dwarfs (BD) with a time-averaged mass transfer rate <\dot M>\approx 10^{-11} M_sun/year. We show that the orbital parameters of the 401 Hz accreting millisecond pulsar SAX J1808.4-3658 support this hypothesis. A main-sequence mass donor requires a nearly face-on inclination and a higher <\dot M> than observed. However, the range of allowed inclinations is substantially relaxed and the predicted mass transfer rate is consistent with that observed if a hot 0.05 M_sun BD is the donor. We then explain the BD radius required (0.13 R_\sun) to fill the Roche lobe. All transiently accreting neutron stars have a minimum luminosity in quiescence set by the time-averaged mass transfer rate onto the neutron star. We show that the constant heating of the BD by this emission can maintain the higher entropy implied by a 0.13 R_\sun radius. Our work bolsters the case that SAX J1808.4-3658 is a progenitor to compact ms radio pulsar binaries (e.g. like those found in 47 Tuc). The very low <\dot M>'s implies that the progenitors to these radio pulsars are long-lived (Gyrs) transient systems, rather than short-lived (Myrs) Eddington-limited accretors. Hence, the accreting progenitor population to millisecond radio pulsars in 47 Tuc could still be present and found in quiescence with Chandra.(abridged)

We report on our survey for rapid (time scale of minutes) photometric variability in symbiotic binaries. These binaries are becoming an increasingly important place to study accretion onto white dwarfs since they are candidate Type Ia supernovae progenitors. Unlike in most cataclysmic variables, the white dwarfs in symbiotics typically accrete from a wind, at rates greater than or equal to 10^{-9} solar masses per year. In order to elucidate the differences between symbiotics and other white dwarf accretors, as well as search for magnetism in symbiotic white dwarfs, we have studied 35 primarily northern symbiotic binaries via differential optical photometry. Our study is the most comprehensive to date of rapid variability in symbiotic binaries. We have found one magnetic accretor, Z And, previously reported by Sokoloski & Bildsten (1999). In four systems (EG And, BX Mon, CM Aql, and BF Cyg), some evidence for flickering at a low level (roughly 10 mmag) is seen for the first time. These detections are, however, marginal. For 25 systems, we place tight upper limits (order of mmag) on both aperiodic and periodic variability, highlighting a major difference between symbiotics and cataclysmic variables. The remaining five of the objects included in our sample (the 2 recurrent novae RS Oph and T CrB, plus CH Cyg, o Ceti, and MWC 560) had previous detections of large-amplitude optical flickering, and we present our extensive observations of these systems in a separate paper. We discuss the impact of our results on the ``standard'' picture of wind-fed accretion, and speculate on the possibility that in most symbiotics, light from quasi-steady nuclear burning on the surface of the white dwarf hides the fluctuating emission from accretion.

We calculate the rapid proton (rp) capture process of hydrogen burning on the surface of an accreting neutron star with an updated reaction network that extends up to Xe, far beyond previous work. In both steady-state nuclear burning appropriate for rapidly accreting neutron stars (such as the magnetic polar caps of accreting X-ray pulsars) and unstable burning of Type I X-ray bursts, we find that the rp process ends in a closed SnSbTe cycle. This prevents the synthesis of elements heavier than Te and has important consequences for X-ray burst profiles, the composition of accreting neutron stars, and potentially galactic nucleosynthesis of light p nuclei.

I show that the recent calculation of Moffatt's regarding the viscous dissipation of a spinning coin overlooked the importance of the finite width of the viscous boundary layer. My new estimates are more in accord with that observed. I also point out that the frequency ``chirp'' of the specially designed toy ``Euler's Disk'' is similar to that expected during the last few minutes of the life of a coalescing binary of two neutron stars. As such, this toy is an excellent desktop demonstration for the expected phenomena.

We investigate whether the magnetic field of an accreting neutron star may be diamagnetically screened by the accreted matter. We assume the freshly accumulated material is unmagnetized, and calculate the rate at which the intrinsic stellar magnetic flux is transported into it by Ohmic diffusion. We calculate the one-dimensional steady-state magnetic field profiles, and show that the magnetic field strength decreases as one moves up through the outer crust and ocean by roughly (Mdot/0.02 Mdot_Edd) orders of magnitude, where Mdot is the accretion rate and Mdot_Edd the Eddington accretion rate. We show that buoyancy instabilities set a limit to the strength of any buried field of roughly 10^10-10^11 G. Our results show that magnetic screening is ineffective for Mdot<0.01 Mdot_Edd, so that, no matter how the accreted material joins onto the star, the underlying stellar field should always be evident. In this respect, we point out the only known persistently-pulsing accreting X-ray millisecond pulsar, SAX J1808.4-3658, has an accretion rate of 10^-3 Mdot_Edd, far below the regime where magnetic screening can play a role. Most steadily accreting neutron stars in low-mass X-ray binaries in our Galaxy accrete at rates where screening would be effective if the simplified magnetic and accretion geometry we adopt were correct. If screened, then the underlying field will emerge after accretion halts, on a timescale of only 100--1000 years, set by the Ohmic diffusion time across the outer crust. It thus seems unlikely that screening alone can explain the low magnetic fields of the millisecond radio pulsars.

The nuclear reactions that occur in the stellar progenitors of white dwarfs (WDs) lead to an internal composition of 12C, 16O, and a ``contaminant'' nucleus, 22Ne. The 22Ne is produced by helium captures on 14N left from hydrogen burning via the CNO cycle. By virtue of its two excess neutrons (relative to the predominant A=2Z nuclei), a downward force of 2m_pg is exerted on 22Ne in the WD interior. This biases its diffusive equilibrium, forcing 22Ne to settle towards the center of the WD. We discuss the physics of the gravitational settling when the WD is in the liquid state and the luminosity generated by it. This modifies the cooling of WD's with masses in excess of a solar mass. The current uncertainties in the microphysics even allow for solutions where a 1.2M_\odot WD remains mostly liquid for a few Gyrs due to the internal heating from 22Ne sedimentation. This highlights the need for an accurate calculation of the inter-diffusion coefficient, especially in the quantum liquid regime relevant for high mass WD's. There is also time in old, liquid WD's (such as those found in cataclysmic variables and possibly in accreting Type Ia progenitors) for partial settling.

We report on spectral and intensity variability analysis from a Chandra/ACIS-S observation of the transient, type-I X-ray bursting low-mass X-ray binary Cen X-4. The quiescent X-ray spectrum during this observation is statistically identical to one observed previously with Beppo/SAX, and close, but not identical, to one observed previously with ASCA. The X-ray spectrum is best described as a pure Hydrogen atmosphere thermal spectrum plus a power-law component that dominates the spectrum above 2 keV. The best-fit radius of the neutron star is r=12.9+/-2.6 (d/1.2 kpc) km if the interstellar absorption is fixed at the value implied by the optical reddening. Allowing the interstellar absorption to be a free parameter yields r=19+45-10 (d/1.2 kpc) km (90% confidence). The thermal spectrum from the neutron star surface is inconsistent with a solar metallicity. We find a 3sigma upper-limit of root-mean-square variability <18% (0.2-2.0 keV; 0.0001-1 Hz) during the observation. On the other hand, the 0.5-10.0 keV luminosity decreased by 40+/-8% in the 4.9 years between the Asca and Chandra observations. This variability can be attributed to the power-law component. Moreover, we limit the variation in thermal temperature to <10% over these 4.9 years. The stability of the thermal temperature and emission area radius supports the interpretation that the quiescent thermal emission is due to the hot neutron star core.

We present an analysis of the nine type I X-ray bursts that were observed from KS 1731-260 with RXTE. We find that the bursts divide naturally into two populations: ``fast bursts'' occur on the Banana Branch when the accretion rate is high and exhibit short decay times, high peak fluxes, and radius expansion episodes. ``Slow bursts'' occur in the Island State at lower accretion rates, have lower peak fluxes, higher fluences, longer decay times, and show no evidence of radius expansion. All five of the fast bursts, and none of the four slow bursts, show coherent oscillations near 524 Hz. We perform in-burst phase connection of the burst pulsations, which allows us to unambiguously characterize their frequency evolution. That evolution exhibits a variety of behaviors, including a sharp spin-down during one burst. Applying our phase models, we find that the pulsations are spectrally harder than the burst emission, with the strength of the pulsations increasing monotonically with photon energy. Coherently summing all of the burst signals, we find upper limits of ~5% relative to the strength of the main pulse on any modulation at 1/2 or 3/2 the main pulse frequency. We discuss our results in the context of models in which the burst pulsations are produced by a hot spot on the surface of a spinning neutron star.

We have searched the rising portion of type I X-ray bursts observed from the Rapid Burster with the Rossi X-ray Timing Explorer for the presence of periodicities. The 95 per cent confidence upper limit on the average root-mean-square variation of near coherent pulsations with a width of <1 Hz (in 60--2048 Hz) during the first second of the bursts is <8.8 per cent. We find a possible detection (>98 per cent significance) at 306.5 Hz.

Motivated by the narrow range of spin frequencies of nearly 20 accreting neutron stars, Bildsten (1998) conjectured that their spin-up had been halted by the emission of gravitational waves. He also pointed out that small nonaxisymmetric temperature variations in the accreted crust will lead to "wavy" electron capture layers, whose horizontal density variations naturally create a mass quadrupole moment. We present a full calculation of the crust's elastic adjustment to these density perturbations and find that the elastic response of the crust reduces Bildsten's original estimate of the quadrupole moment in the thin outer crust by a factor of 20-50. However, this basic picture, when applied to capture layers in the deep inner crust, can generate quadrupoles in the necessary range as long as there are ~5% lateral temperature variations in the inner crust. By calculating the thermal flow throughout the core and the crust, we find that temperature gradients this large are easily maintained by asymmetric heat sources or lateral composition gradients in the crust. We also derive a general relation between the stresses and strains in the crust and the maximum quadrupole moment they can generate. We show under quite general conditions that maintaining a quadrupole of the magnitude necessary to balance the accretion torque requires dimensionless strains close to 0.01 at near-Eddington accretion rates, of order the breaking strain of conventional materials.

The rotation rates of six weakly-magnetic neutron stars accreting in low-mass X-ray binaries have most likely been measured by Type I X-ray burst observations with RXTE. The nearly coherent oscillations detected during the few seconds of thermonuclear burning are most simply understood as rotational modulation of brightness asymmetries on the neutron star surface. We show that, as suggested by Strohmayer and colleagues, the frequency changes of 1-2 Hz observed during bursts are consistent with angular momentum conservation as the burning shell hydrostatically expands and contracts. We calculate how vertical heat propagation through the radiative outer layers of the atmosphere and convection affect the coherence of the oscillation. We show that the evolution of the rotational profile depends strongly on whether the burning layers are composed of pure helium or mixed hydrogen/helium. Our results help explain the absence (presence) of oscillations from hydrogen-burning (helium-rich) bursts that was found by Muno and collaborators. We investigate angular momentum transport within the burning layers and the recoupling of the burning layers with the star. We show that the Kelvin-Helmholtz instability is quenched by the strong stratification, and that mixing between the burning fuel and underlying ashes by the baroclinic instability does not occur. However, the baroclinic instability may have time to operate within the differentially rotating burning layer, potentially bringing it into rigid rotation.

We have detected a bright X-ray flare from the nearby brown dwarf LP 944-20 with the Chandra/ACIS-S. This is an old (500 Myr), rapidly rotating, lithium-bearing M9 object, with a bolometric luminosity of 6e29 ergs cm-2 s-1. It was only detected by Chandra during an X-ray flare of duration 1-2 hours near the end of a 12.1 hour observation. The peak X-ray luminosity was 1.2(+0.5/-0.3)e26 ergs cm-2 s-1 in the brightest ~550 seconds, corresponding to Lx/Lbol ~ 2e-4. A total of 2e29 ergs was released during the 43,773 sec observation, giving a time-averaged Lx/Lbol ~ 7e-6. LP 944-20 was not detected before the flare, with a 3 sigma upper limit on the emission at Lx/Lbol < 2e-6 (Lx<1e24 ergs cm-2 s-1). This is faint for a rapidly rotating late-type star, and establishes a record lower limit to the quiescent flux about an order of magnitude below the flux limit (and a factor of 5 below the Lx/Lbol limit) placed on quiescent X-ray emission from the M8 dwarf VB 10. The inferred flaring duty cycle is comparable to that measured via variable H-alpha emission for other late M-type, fully convective stars.

We summarize the quiescent X-ray observations of transient low-mass X-ray binaries. These observations show that, in quiescence, binaries containing black holes are fainter than those containing neutron stars. This has triggered a number of theoretical ideas about what causes the quiescent X-ray emission. For black hole binaries, the options are accretion onto the black hole or coronal emission from the rapidly rotating stellar companion. There are more possibilities for the neutron stars: accretion, thermal emission from the surface or non-thermal emission from a ``turned-on'' radio pulsar. We review recent theoretical work on these mechanisms and note where current observations can distinguish between them. We highlight the re-analysis of the quiescent neutron star emission by Rutledge and collaborators that showed thermal emission to be a predominant contributor in many of these systems. Our knowledge of these binaries is bound to dramatically improve now that the Chandra and XMM satellites are operating successfully,

Observations of soft X-ray transients (SXTs) in quiescence have found that the binaries harboring black holes are fainter than those that contain a neutron star. Narayan and collaborators postulated that the faint X-ray emission from black hole binaries was powered by an advection dominated accretion flow (ADAF). We explore an alternative explanation for the quiescent X-ray emission from the black hole systems: coronal emission from the rapidly rotating optical companion. This is commonly observed and well studied in other tidally locked binaries, such as the RS CVns. We show that two of the three X-ray detected black hole binaries (A0620-00 and GRO J1655-40) exhibit X-ray fluxes entirely consistent with coronal emission. The X-ray spectra of these objects should be best fit with thermal Raymond-Smith models rich in lines when coronal emission predominates. One black hole system (V404 Cyg) is too X-ray bright to be explained as coronal emission. The quiescent X-ray emission from the neutron star binaries is far too bright for coronal emission. It might be that all SXT's have variable accretion rates in quiescence and that the basal quiescent X-ray flux is set by either coronal emission from the companion or -- when present -- by thermal emission from the neutron star. We also show that the lithium abundances in the black hole systems are comparable to those in the RS CVns, reducing the need for production mechanisms that involve the compact object.

I review our understanding of the thermonuclear instabilities on accreting neutron stars that produce Type I X-Ray bursts. I emphasize those observational and theoretical aspects that should interest the broad audience of this meeting. The easily accessible timescales of the bursts (durations of tens of seconds and recurrence times of hours to days) allow for a very stringent comparison to theory. The largest discrepancy (which was found with EXOSAT observations) is the accretion rate dependence of the Type I burst properties. Bursts become less frequent and energetic as the global accretion rate increases, just the opposite of what the spherical theory predicts. I present a resolution of this issue by taking seriously the observed dependence of the burning area on the global accretion rate, which implies that as the accretion rate increases, the accretion rate per unit area decreases. This resurrects the unsolved problem of knowing where the freshly accreted material accumulates on the star, equally relevant to the likely signs of rotation during the bursts summarized by Swank at this meeting. I close by highlighting the Type I bursts from GS 1826-238 that were found with BeppoSAX and RXTE. Their energetics, recurrence times and temporal profiles clearly indicate that hydrogen is being burned during these bursts, most likely by the rapid-proton (rp) process.

We summarize the observations of the spin periods of rapidly accreting neutron stars. If gravitational radiation is responsible for balancing the accretion torque at the observed spin frequencies of ~300 Hz, then the brightest of these systems make excellent gravitational wave sources for LIGO-II and beyond. We review the recent theoretical progress on two mechanisms for gravitational wave emission: mass quadrupole radiation from deformed neutron star crusts and current quadrupole radiation from r-mode pulsations in neutron star cores.

A deep all-sky imaging hard x-ray survey and wide-field monitor is needed to extend soft (ROSAT) and medium (ABRIXAS2) x-ray surveys into the 10-100 keV band (and beyond) at comparable sensitivity (~0.05 mCrab). This would enable discovery and study of >3000 obscured AGN, which probably dominate the hard x-ray background; detailed study of spectra and variability of accreting black holes and a census of BHs in the Galaxy; Gamma-ray bursts and associated massive star formation (PopIII) at very high redshift and Soft Gamma-ray Repeaters throughout the Local Group; and a full galactic survey for obscured supernova remnants. The Energetic X-ray Imaging Survey Telescope (EXIST) is a proposed array of 8 x 1m^2 coded aperture telescopes fixed on the International Space Station (ISS) with 160deg x 40deg field of view which images the full sky each 90 min orbit. EXIST has been included in the most recent NASA Strategic Plan as a candidate mission for the next decade. An overview of the science goals and mission concept is presented.

Recent work has raised the exciting possibility that r-modes (Rossby waves) in rotating neutron star cores might be strong gravitational wave sources. We estimate the effect of a solid crust on their viscous damping rate and show that the dissipation rate in the viscous boundary layer between the oscillating fluid and the nearly static crust is >10^5 times higher than that from the shear throughout the interior. This increases the minimum frequency for the onset of the gravitational r-mode instability to at least 500 Hz when the core temperature is less than 10^10 K. It eliminates the conflict of the r-mode instability with the accretion-driven spin-up scenario for millisecond radio pulsars and makes it unlikely that the r-mode instability is active in accreting neutron stars. For newborn neutron stars, the formation of a solid crust shortly after birth affects their gravitational wave spin-down and hence detectability by ground-based interferometric gravitational wave detectors.

We fit hydrogen atmosphere models to the X-ray data for four neutron stars (three from a previous paper, plus 4U 2129+47) and six black hole candidates (A0620-00, GS 2000+25, GS 1124-68, GS 2023+33, GRO J1655-40, and GRO J0422+32). While the neutron stars are similar in their intrinsic X-ray spectra (similar effective temperatures and emission area radii ~10 km), the spectra of two black hole candidates are significantly different, and the spectra of the remaining four are consistent with a very large parameter space that includes the neutron stars. The spectral differences between the neutron stars and black hole candidates favors the interpretation that the quiescent neutron star emission is predominantly thermal emission from the neutron star surface. Our work suggests that an X-ray spectral comparison in quiescence provides an additional means for distinguishing between neutron stars and black holes. The faint X-ray sources in globular clusters are also a class of objects which can be investigated in this manner.

We calculate the nucleosynthesis during stable nuclear burning on an accreting neutron star. This is appropriate for weakly magnetic neutron stars accreting at near-Eddington rates in low mass X-ray binaries, and for most accreting X-ray pulsars. We show that the nuclear burning proceeds via the rapid proton capture process (rp process), and makes nuclei far beyond the iron group. The final mixture of nuclei consists of elements with a range of masses between approximately A=60 and A=100. The average nuclear mass of the ashes is set by the extent of helium burning via (alpha,p) reactions, and depends on the local accretion rate. Our results imply that the crust of these accreting neutron stars is made from a complex mixture of heavy nuclei, with important implications for its thermal, electrical and structural properties. A crustal lattice as impure as our results suggest will have a conductivity set mostly by impurity scattering, allowing more rapid Ohmic diffusion of magnetic fields than previously estimated.

We describe the evolutionary progression of an outburst of the Rapid Burster. Four outbursts have been observed with the Rossi X-Ray Timing Explorer between February 1996 and May 1998, and our observations are consistent with a standard evolution over the course of each. An outburst can be divided into two distinct phases: Phase I is dominated by type I bursts, with a strong persistent emission component; it lasts for 15-20 days. Phase II is characterized by type II bursts, which occur in a variety of patterns. The light curves of time-averaged luminosity for the outbursts show some evidence for reflares, similar to those seen in soft X-ray transients. The average recurrence time for Rapid Burster outbursts during this period has been 218 days, in contrast with an average ~180 day recurrence period observed during 1976-1983.

We present observations of a series of 10 outbursts of pulsed hard X-ray flux from the transient 10.6 mHz accreting pulsar GS 1843-02, using the Burst and Transient Source Experiment on the Compton Gamma Ray Observatory. These outbursts occurred regularly every 242 days, coincident with the ephemeris of the periodic transient GRO J1849-03 (Zhang et al. 1996), which has recently been identified with the SAS 3 source 2S 1845-024 (Soffitta et al. 1998). Our pulsed detection provides the first clear identification of GS 1843-02 with 2S 1845-024. We present a pulse timing analysis which shows that the 2S 1845-024 outbursts occur near the periastron passage of the neutron star's highly eccentric (e = 0.88+-0.01) 242.18+-0.01 day period binary orbit about a high mass (M > 7 solar masses) companion. The orbit and transient outburst pattern strongly suggest the pulsar is in a binary system with a Be star. Our observations show a long-term spin-up trend, with most of the spin-up occurring during the outbursts. From the measured spin-up rates and inferred luminosities we conclude that an accretion disk is present during the outbursts.

We report the first result from our survey of rapid variability in symbiotic binaries: the discovery of a persistent oscillation at P=1682.6 +- 0.6 s in the optical emission from the prototype symbiotic, Z Andromedae. The oscillation was detected on all 8 occasions on which the source was observed over a timespan of nearly a year, making it the first such persistent periodic pulse found in a symbiotic binary. The amplitude was typically 2 - 5 mmag, and it was correlated with the optical brightness during a relatively small outburst of the system. The most natural explanation is that the oscillation arises from the rotation of an accreting, magnetic (B_S > 10^5 G) white dwarf. This discovery constrains the outburst mechanisms, since the oscillation emission region near the surface of the white dwarf was visible during the outburst.

We re-analyze the available X-ray spectral data of the type I bursting neutron star transients Aql X-1, Cen X-4, and 4U 1608-522 using realistic hydrogen atmosphere models. Previous spectral fits assumed a blackbody spectrum; because the free-free dominated photospheric opacity decreases with increasing frequency, blackbody spectral fits overestimate the effective temperature and underestimate, by as much as two orders of magnitude, the emitting area. Hydrogen atmosphere spectral models, when fit to the available observational data, imply systematically larger emission area radii, consistent with the canonical 10 km radius of a neutron star. This suggests that a substantial fraction of the quiescent luminosity is thermal emission from the surface of the neutron star. The magnitude of the equivalent hydrogen column density toward these systems, however, presents a considerable systematic uncertainty, which can only be eliminated by high signal-to-noise X-ray spectral measurements (e.g., with AXAF or XMM) which would permit simultaneous determination of the equivalent hydrogen column density, emission area, and thermal temperature.

Nuclear reactions occurring deep in the crust of a transiently accreting neutron star efficiently maintain the core at a temperature >5e7 K. When accretion halts, the envelope relaxes to a thermal equilibrium set by the flux from the hot core, as if the neutron star were newly born. For the time-averaged accretion rates typical of low-mass X-ray transients, standard neutrino cooling is unimportant and the core thermally re-radiates the deposited heat. The resulting luminosity has the same magnitude as that observed from several transient neutron stars in quiescence. Confirmation of this mechanism would strongly constrain rapid neutrino cooling mechanisms for neutron stars. Thermal emission had previously been dismissed as a predominant source of quiescent emission since blackbody spectral fits implied an emitting area much smaller than a neutron star's surface. However, as with thermal emission from radio pulsars, fits with realistic emergent spectra will imply a substantially larger emitting area. Other emission mechanisms, such as accretion or a pulsar shock, can also operate in quiescence and generate intensity and spectral variations over short timescales. Indeed, quiescent accretion may produce gravitationally redshifted metal photoionization edges in the quiescent spectra (detectable with AXAF and XMM). We discuss past observations of Aql~X-1 and note that the low luminosity X-ray sources in globular clusters and the Be star/X-ray transients are excellent candidates for future study.

We investigate hydrogen electron capture in the oceans of neutron stars accreting at rates appropriate for most Low-Mass X-Ray Binaries. These stars burn the accreted hydrogen and helium unstably in the upper atmosphere and accumulate material which usually contains some small amount of hydrogen (mass fractions are typically 10%) mixed in with the heavier iron group ashes. The subsequent evolution of this matter is determined by compression towards higher densities until electron capture on the hydrogen occurs. We construct steady-state models of the electron captures and the subsequent neutron recombinations onto the heavy nuclei. The density discontinuity from these captures gives rise to a new g-mode (much like a surface wave), which has a lowest order (l=1) frequency of 35 Hz when X=0.1 on a slowly rotating star. We also discuss, for the first time, a new set of non-radial g-modes unique to these neutron stars. These modes have most of their nodes in the finite thickness layer where the electron captures are occurring. The lowest order mode frequencies are in the 1-10 Hz range for a few radial nodes on a slowly rotating star. We conclude by discussing how the dispersion relations for these modes are modified for a rapidly rotating neutron star. Whether any of these modes are observable depends on their excitation and damping mechanisms and the ability to excite them, issues we will address in a future paper.

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.

Several classes of stars (most notably O and B main-sequence stars, as well as accreting white dwarfs and neutron stars) rotate quite rapidly, at spin frequencies greater than the typical g-mode frequencies. We discuss how rapid rotation modifies the $\kappa$-mechanism excitation and observability of g-mode oscillations. We find that, by affecting the timescale match between the mode period and the thermal time at the driving zone, rapid rotation stabilizes some of the g-modes that are excited in a non-rotating star, and, conversely, excites g-modes that are damped in absence of rotation. The fluid velocities and temperature perturbations are strongly concentrated near the equator for most g-modes in rapidly rotating stars, which means that a favorable viewing angle may be required to observe the pulsations. Moreover, the stability of modes of the same $l$ but different $m$ is affected differently by rotation. We illustrate this by considering g-modes in Slowly Pulsating B-type stars as a function of the rotation rate.

We investigate the atmosphere, ocean, and crust of neutron stars accreting at rates sufficiently high (typically in excess of the local Eddington limit) to stabilize the burning of accreted hydrogen and helium. For hydrogen-rich accretion at global rates in excess of 10^-8 solar masses per year (typical of a few neutron stars), we discuss the thermal state of the deep ocean and crust and their coupling to the neutron star core, which is heated by conduction (from the crust) and cooled by neutrino emission. We estimate the Ohmic diffusion time in the hot, deep crust and find that it is noticeably shortened (to less than 10^8 yr) from the values characteristic of the colder crusts in slowly accreting neutron stars. We speculate on the implications of these calculations for magnetic field evolution in the bright accreting X-ray sources. We also explore the consequences of rapid compression at local accretion rates exceeding ten times the Eddington rate. This rapid accretion heats the atmosphere/ocean to temperatures of order 10^9 K at relatively low densities; for stars accreting pure helium, this causes unstable ignition of the ashes (mostly carbon) resulting from stable helium burning. This unstable burning can re-occur on timescales as short as hours to days, and might be the cause of some flares on helium accreting pulsars, in particular 4U~1626--67. Such rapid local accretion rates are common on accreting X-ray pulsars, where the magnetic field focuses the accretion flow onto a small fraction of the stellar area.

I review the rotation measurements of accreting neutron stars. Many of the highly magnetic accreting X-ray pulsars have been continuously observed with the Burst and Transient Source Experiment (BATSE) aboard the Compton Gamma-Ray Observatory (CGRO) since April 1991. These observations show that the accretion torque exerted on many disk-fed accreting X-ray pulsars changes sign on a monthly to yearly timescale. This results in alternating periods of spin-up and spin-down with nearly the same torques, leading to little net angular momentum gained by accretion. I also summarize recent discoveries with the Rossi X-Ray Timing Explorer (RXTE) of periodicities during Type I X-ray bursts. These seem to indicate that many of the rapidly accreting and weakly magnetic neutron stars in our galaxy are rotating at frequencies greater than 250 Hertz. Most remarkable is that they all rotate within a rather narrow range of frequencies.

Observations of quasi-periodic oscillations (QPOs) in the luminosity from many accreting neutron stars (NS) have led us to investigate a source of periodicity prevalent in other stars: non-radial oscillations. After summarizing the structure of the atmosphere and ocean of an accreting NS, we discuss the various low l g-modes with frequencies in the 1-100 Hz range. Successful identification of a non-radial mode with an observed frequency would yield new information about the thermal and compositional makeup of the NS, as well as its radius. We close by discussing how rapid rotation changes the g-mode frequencies.

We present an analytic calculation of the thermonuclear depletion of the light elements lithium, beryllium, and boron in fully convective, low-mass stars. Under the presumption that the pre--main-sequence star is always fully mixed during contraction, we find that the burning of these rare light elements can be computed analytically, even when the star is degenerate. Using the effective temperature as a free parameter, we constrain the properties of low-mass stars from observational data, independently of the uncertainties associated with modeling their atmospheres and convection. Our analytic solution explains the dependence of the age at a given level of elemental depletion on the stellar effective temperature, nuclear cross sections, and chemical composition. Most importantly, our results allow observers to translate lithium non-detections in young cluster members into a model-independent minimum age for that cluster. Using this procedure, we have found lower limits to the ages of the Pleiades (100 Myr) and Alpha Persei (60 Myr) clusters. Recent experimental work on the low energy resonance in the ^10B(p,\alpha)^7Be reaction has greatly enhanced estimates of the destruction rate of ^10B, making it possible for stars with M>0.1 M_sun to deplete both ^10B and ^11B before reaching the main sequence. Moreover, there is an interesting range of masses, 0.085 M_sun < M < 0.13 M_sun, where boron depletion occurs on the main sequence in less than a Hubble time, providing a potential ``clock'' for dating low-mass stars.

Neutron stars in mass-transferring binaries are accreting the hydrogen and helium rich matter from the surfaces of their companions. This article simply explains the physics associated with how that material eventually fuses to form heavier nuclei and the observations of the time dependent phenomena (such as Type I X-ray bursts) associated with the thermally unstable thermonuclear reactions. We explain how the outcome depends on the composition of the accreting matter, the accretion rate and the mass, radius and thermal state of the neutron star. We also introduce many new analytic relations that are convenient for comparisons to both observations and computational results. After explaining nuclear burning for spherically symmetric accretion onto neutron stars, we discuss the possibility of asymmetric burning. In particular, we discuss some of the mysteries from EXOSAT observations of Type I X-Ray bursts and how the solution to these puzzles may lie in considering the lateral propagation of nuclear burning fronts around the star. Fully understanding this problem requires knowledge of parameters previously neglected such as the distribution of fresh fuel on the star, the magnetic field strength, and the stellar rotation. Recent RXTE observations of bursters may finally tell us some of these parameters.

Dramatic torque reversals between spin up and spin down have been observed in half of the persistent X-ray pulsars monitored by the BATSE all-sky monitor on CGRO. Theoretical models developed to explain early pulsar timing data can explain spin down torques via a disk-magnetosphere interaction if the star nearly corotates with the inner accretion disk. To produce the observed BATSE torque reversals, however, these equilibrium models require the disk to alternate between two mass accretion rates, with $\dot M_{\pm}$ producing accretion torques of similar magnitude, but always of opposite sign. Moreover, in at least one pulsar (GX 1+4) undergoing secular spin down the neutron star spins down faster during brief ($\sim 20$ day) hard X-ray flares -- this is opposite the correlation expected from standard theory, assuming BATSE pulsed flux increases with mass accretion rate. The $10$ day to 10 yr intervals between torque reversals in these systems are much longer than any characteristic magnetic or viscous time scale near the inner disk boundary and are more suggestive of a global disk phenomenon. We discuss possible explanations of the observed torque behavior. Despite the preferred sense of rotation defined by the binary orbit, the BATSE observations are surprisingly consistent with an earlier suggestion by Makishima \etal (1988) for GX~1+4: the disks in these systems somehow alternate between episodes of prograde and retrograde rotation. We are unaware of any mechanism that could produce a stable retrograde disk in a binary undergoing Roche-lobe overflow, but such flip-flop behavior does occur in numerical simulations of wind-fed systems. One possibility is that the disks in some of these binaries are fed by an X-ray excited wind.

We summarize five years of continuous monitoring of accretion-powered pulsars with the Burst and Transient Source Experiment (BATSE) on the Compton Gamma Ray Observatory. Our 20-70 keV observations have determined or refined the orbital parameters of 13 binaries, discovered 5 new transient accreting pulsars, measured the pulsed flux history during outbursts of 12 transients (GRO J1744-28, 4U 0115+634, GRO J1750-27, GS 0834-430, 2S 1417-624, GRO J1948+32, EXO 2030+375, GRO J1008-57, A 0535+26, GRO J2058+42, 4U 1145-619 and A 1118-616), and also measured the accretion torque history of during outbursts of 6 of those transients whose orbital parameters were also known. We have also continuously measured the pulsed flux and spin frequency for eight persistently accreting pulsars (Her X-1, Cen X-3, Vela X-1, OAO 1657-415, GX 301-2, 4U 1626-67, 4U 1538-52, and GX 1+4). Because of their continuity and uniformity over a long baseline, BATSE observations have provided new insights into the long-term behavior of accreting magnetic stars. We have found that all accreting pulsars show stochastic variations in their spin frequencies and luminosities, including those displaying secular spin-up or spin-down on long time scales, blurring the conventional distinction between disk-fed and wind-fed binaries. Pulsed flux and accretion torque are strongly correlated in outbursts of transient accreting pulsars, but uncorrelated, or even anticorrelated, in persistent sources.

Over five years of daily hard X-ray (>20 keV) monitoring of the 2-min accretion-powered pulsar GX 1+4 with the Compton Gamma Ray Observatory/BATSE large-area detectors has found nearly continuous rapid spin-down, interrupted by a bright 200-d spin-up episode. During spin-down, the torque becomes more negative as the luminosity increases (assuming that the 20-60 keV pulsed flux traces bolometric luminosity), the opposite of what is predicted by standard accretion torque theory. No changes in the shape of the 20-100 keV pulsed energy spectrum were detected, so that a very drastic change in the spectrum below 20 keV or the pulsed fraction would be required to make the 20-60 keV pulsed flux a poor luminosity tracer. These are the first observations which flatly contradict standard magnetic disk accretion theory, and they may have important implications for understanding the spin evolution of X-ray binaries, cataclysmic variables, and protostars. We briefly discuss the possibility that GX 1+4 may be accreting from a retrograde disk during spin-down, as previously suggested.

We present an analytic calculation of the thermonuclear depletion of lithium in contracting, fully convective, pre-main sequence stars of mass M < 0.5 M_sun. Previous numerical work relies on still-uncertain physics (atmospheric opacities and convection, in particular) to calculate the effective temperature as a unique function of stellar mass. We assume that the star's effective temperature, T_eff, is fixed during Hayashi contraction and allow its actual value to be a free parameter constrained by observation. Using this approximation, we compute lithium burning analytically and explore the dependence of lithium depletion on T_eff, M, and composition. Our calculations yield the radius, age, and luminosity of a pre-main sequence star as a function of lithium depletion. This allows for more direct comparisons to observations of lithium depleted stars. Our results agree with those numerical calculations that explicitly determine stellar structure during Hayashi contraction. In agreement with Basri, Marcy, and Graham (1996), we show that the absence of lithium in the Pleiades star HHJ 3 implies that it is older than 100 Myr. We also suggest a generalized method for dating galactic clusters younger than 100 Myr (i.e., those with contracting stars of M > 0.08 M_sun) and for constraining the masses of lithium depleted stars.

We investigate the thermal stability of nuclear burning on the accreting X-ray pulsar GRO J1744-28. The neutron star's dipolar magnetic field is <3\times 10^{11} G if persistent spin-up implies that the magnetospheric radius is less than the co-rotation radius. After inferring the properties of the neutron star, we study the thermal stability of hydrogen/helium burning and show that thermonuclear instabilities are unlikely causes of the hourly bursts seen at very high accretion rates. We then discuss how the stability of the thermonuclear burning depends on both the global accretion rate and the neutron star's magnetic field strength. We emphasize that the appearance of the instability (i.e., whether it looks like a Type I X-ray burst or a flare lasting a few minutes) will yield crucial information on the neutron star's surface magnetic field and the role of magnetic fields in convection. We suggest that a thermal instability in the accretion disk is the origin of the long (~300 days) outburst and that the recurrence time of these outbursts is >50 years. We also discuss the nature of the binary and point out that a velocity measurement of the stellar companion (most likely a Roche-lobe filling giant with m_K>17) will constrain the neutron star mass.

The spin-induced quadrupole moment of a rapidly rotating star changes the orbital dynamics in a binary system, giving rise to advance (or regression) of periastron and precession of the orbital plane. We show that these effects are important in the recently discovered radio pulsar/main sequence star binary system PSR J0045$-$7319, and can reliably account for the observed peculiar timing residuals. Precise measurements of the apsidal motion and orbital plane precession can yield valuable information on the internal structure and rotation of the star. The detection of orbital precession implies that the spin of the companion star is not aligned with the orbital angular momentum, and suggests that the supernova gave the pulsar a kick out of the original orbital plane. Tidal excitation of g-mode oscillations in the PSR J0045$-$7319 system induces an orbital period change of order $|\Delta P_{\rm orb}/P_{\rm orb}|\sim 10^{-6}$ at each periastron passage, but the secular trend depends on the radiative damping time of the g-modes. We also discuss the spin-orbit coupling effects for the accreting X-ray pulsars and the other known radio pulsar/main sequence binary, PSR B1259$-$63.

Gravitational-wave interferometers are expected to monitor the last three minutes of inspiral and final coalescence of neutron star and black hole binaries at distances approaching cosmological, where the event rate may be many per year. Because the binary's accumulated orbital phase can be measured to a fractional accuracy $\ll 10^{-3}$ and relativistic effects are large, the waveforms will be far more complex, carry more information, and be far harder to model theoretically than has been expected. Theorists must begin now to lay a foundation for extracting the waves' information.