We report on time-resolved optical imaging of the X-ray binary SAX J1808.4-3658 during its quiescent state and 2008 outburst. The binary, containing an accretion-powered millisecond pulsar, has a large sinusoidal-like modulation in its quiescent optical emission. We employ a Markov chain Monte Carlo technique to fit our multi-band light curve data in quiescence with an irradiated star model, and derive a tight constraint of $50^{+6}_{-5}$ deg on the inclination angle $i$ of the binary system. The pulsar and its companion are constrained to have masses of $0.97^{+0.31}_{-0.22} M_{\sun}$ and $0.04^{+0.02}_{-0.01} M_{\sun}$ (both 1$\sigma$ ranges), respectively. The dependence of these results on the measurements of the companion's projected radial velocity is discussed. We also find that the accretion disk had nearly constant optical fluxes over a $\sim$500 day period in the quiescent state our data covered, but started brightening 1.5 months before the 2008 outburst. Variations in modulation during the outburst were detected in our four observations made 7-12 days after the start of the outburst, and a sinusoidal-like modulation with 0.2 mag amplitude changed to have a smaller amplitude of 0.1 mag. The modulation variations are discussed. We estimate the albedo of the companion during its quiescence and the outburst, which was approximately 0 and 0.8 (for isotropic emission), respectively. This large difference probably provides additional evidence that the neutron star in the binary turns on as a radio pulsar in quiescence.

Galactic interacting double white dwarfs (DWD) are guaranteed gravitational wave (GW) sources for the GW detector LISA, with more than 10^4 binaries expected to be detected over the mission's lifetime. Part of this population is expected to be eccentric, and here we investigate the potential for constraining the white dwarf (WD) properties through apsidal precession in these binaries. We analyze the tidal, rotational, and general relativistic contributions to apsidal precession by using detailed He WD models, where the evolution of the star's interior is followed throughout the cooling phase. In agreement with previous studies of zero-temperature WDs, we find that apsidal precession in eccentric DWDs can lead to a detectable shift in the emitted GW signal when binaries with cool (old) components are considered. This shift increases significantly for hot (young) WDs. We find that apsidal motion in hot (cool) DWDs is dominated by tides at orbital frequencies above ~10^{-4}Hz (10^{- 3}$Hz). The analysis of apsidal precession in these sources while ignoring the tidal component would lead to an extreme bias in the mass determination, and could lead us to misidentify WDs as neutron stars or black holes. We use the detailed WD models to show that for older, cold WDs, there is a unique relationship that ties the radius and apsidal precession constant to the WD masses, therefore allowing tides to be used as a tool to constrain the source masses.

We have developed a new method for calculating common envelope (CE) events based on explicit consideration of the donor star's structural response to adiabatic mass loss. In contrast to existing CE prescriptions, which specify a priori the donor's remnant mass, we determine this quantity self-consistently and find it depends on binary and CE parameters. This aspect of our model is particularly important to realistic modeling for upper main sequence star donors without strongly degenerate cores (and hence without a clear core/envelope boundary). We illustrate the central features of our method by considering CE events involving 10 solar mass donors on or before their red giant branch. For such donors, the remnant core mass can be as much as 30% larger than the star's He-core mass. Applied across a population of such binaries, our methodology results in a significantly broader remnant mass and final orbital separation distribution and a 20% increase in CE survival rates as compared to previous prescriptions for the CE phase.

We study tidal interactions in white dwarf binaries in the limiting case of quasi-static tides. The formalism is valid for arbitrary orbital eccentricities and therefore applicable to white dwarf binaries in the Galactic disk as well as globular clusters. In the quasi-static limit, the total perturbation of the gravitational potential shows a phase shift with respect to the position of the companion, the magnitude of which is determined primarily by the efficiency of energy dissipation through convective damping. We determine rates of secular evolution of the orbital elements and white dwarf rotational angular velocity for a 0.3 solar mass helium white dwarf in binaries with orbital frequencies in the LISA gravitational wave frequency band and companion masses ranging from 0.3 to 10^5 solar masses. The resulting tidal evolution time scales for the orbital semi-major axis are longer than a Hubble time, so that convective damping of quasi-static tides need not be considered in the construction of gravitational wave templates of white dwarf binaries in the LISA band. Spin-up of the white dwarf, on the other hand, can occur on time scales of less than 10Myr, provided that the white dwarf is initially rotating with a frequency much smaller than the orbital frequency. For semi-detached white dwarf binaries spin-up can occur on time scales of less than 1Myr. Nevertheless, the time scales remain longer than the orbital inspiral time scales due to gravitational radiation, so that the degree of asynchronism in these binaries increases. As a consequence, tidal forcing eventually occurs at forcing frequencies beyond the quasi-static tide approximation. For the shortest period binaries, energy dissipation is therefore expected to take place through dynamic tides and resonantly excited g-modes.

This white paper briefly describes the astrophysics of ultra-compact binaries, with emphasis of the challenges and opportunities in the next decade.

We observed the accreting millisecond pulsar SAX J1808.4-3658 with Gemini-South in g' and i' bands, nearly simultaneous with XMM-Newton observations. A clear periodic flux modulation on the system's orbital period is present, consistent with the varying aspect of the donor star's heated face. We model the contributions of a disk and donor star to these optical bands. To produce the observed modulation amplitudes, we conclude that the donor must be irradiated by an external flux 2 orders of magnitude greater than provided by the measured X-ray luminosity. A possible explanation for this irradiation is that the radio pulsar mechanism becomes active during the quiescent state as suggested by Burderi et al., with relativistic particles heating the donor's day-side face. Our modelling constrains the binary inclination to be between 36 and 67 degrees. We obtain estimates for the pulsar mass of >2.2 solar masses (although this limit is sensitive to the source's distance), consistent with the accelerated NS cooling in this system indicated by X-ray observations. We also estimate the donor mass to be in the range of 0.07-0.11 solar masses, providing further indications that the system underwent non-standard binary evolution to reach its current state.

We observed SAX J1808.4-3658 (1808), the first accreting millisecond pulsar, in deep quiescence with XMM-Newton and (near-simultaneously) Gemini-South. The X-ray spectrum of 1808 is similar to that observed in quiescence in 2001 and 2006, describable by an absorbed power-law with photon index 1.74+-0.11 and unabsorbed X-ray luminosity L_X=7.9+-0.7*10^{31} ergs/s, for N_H=1.3*10^{21} cm^{-2}. Fitting all the quiescent XMM-Newton X-ray spectra with a power-law, we constrain any thermally emitting neutron star with a hydrogen atmosphere to have a temperature less than 30 eV and L_{NS}(0.01-10 keV)<6.2*10^{30} ergs/s. A thermal plasma model also gives an acceptable fit to the continuum. Adding a neutron star component to the plasma model produces less stringent constraints on the neutron star; a temperature of 36^{+4}_{-8} eV and L_{NS}(0.01-10 keV)=1.3^{+0.6}_{-0.8}*10^{31} ergs/s. In the framework of the current theory of neutron star heating and cooling, the constraints on the thermal luminosity of 1808 and 1H 1905+000 require strongly enhanced cooling in the cores of these neutron stars. We compile data from the literature on the mass transfer rates and quiescent thermal flux of the largest possible sample of transient neutron star LMXBs. We identify a thermal component in the quiescent spectrum of the accreting millisecond pulsar IGR J00291+5934, which is consistent with the standard cooling model. The contrast between the cooling rates of IGR J00291+5934 and 1808 suggests that 1808 may have a significantly larger mass. This can be interpreted as arising from differences in the binary evolution history or initial neutron star mass in these otherwise similar systems.

The accreting millisecond pulsar SAX J1808.4-3658 may be a transition object between accreting X-ray binaries and millisecond radio pulsars. We have constrained the thermal radiation from its surface through XMM-Newton X-ray observations, providing strong evidence for neutrino cooling processes from the neutron star core. We have also undertaken simultaneous X-ray and optical (Gemini) observations, shedding light on whether the strong heating of the companion star in quiescence may be due to X-ray irradiation, or to a radio pulsar turning on when accretion stops.

I review the evolutionary connection between low-mass X-ray binaries (LMXBs) and pulsars with binary companions (bPSRs) from a stellar binary evolution perspective. I focus on the evolution of stellar binaries with end-states consisting of a pulsar with a low-mass (<1.0 solar mass) companion, starting at the point the companion's progenitor first initiates mass transfer onto the neutron star. Whether this mass transfer is stable and the physics driving ongoing mass transfer partitions the phase space of the companions's initial mass and initial orbital period into five regions. The qualitative nature of the mass-transfer process and the binary's final end-state differ between systems in each region; four of these regions each produce a particular class of LMXBs. I compare the theoretical expectations to the populations of galactic field LMXBs with companion-mass constraints and field bPSRs. I show that the population of accreting millisecond pulsars are all identified with only two of the four LMXB classes and that these systems do not have readily identifiable progeny in the bPSR population. I discuss which sub-populations of bPSRs can be explained by binary evolution theory and those that currently are not. Finally I discuss some outstanding questions in this field.

We report on the discovery by the Swift Gamma-Ray Burst Explorer of the eighth known transient accretion-powered millisecond pulsar, SWIFT J1756.9-2508, as part of routine observations with the Swift Burst Alert Telescope hard X-ray transient monitor. The pulsar was subsequently observed by both the X-Ray Telescope on Swift and the Rossi X-Ray Timing Explorer Proportional Counter Array. It has a spin frequency of 182 Hz (5.5 ms) and an orbital period of 54.7 minutes. The minimum companion mass is between 0.0067 and 0.0086 solar masses, depending on the mass of the neutron star, and the upper limit on the mass is 0.030 solar masses (95% confidence level). Such a low mass is inconsistent with brown dwarf models, and comparison with white dwarf models suggests that the companion is a He-dominated donor whose thermal cooling has been at least modestly slowed by irradiation from the accretion flux. No X-ray bursts, dips, eclipses or quasi-periodic oscillations were detected. The current outburst lasted approximately 13 days and no earlier outbursts were found in archival data.

(Abridged) We calculate the full stellar-structural evolution of donors in AM CVn systems formed through the WD channel coupled to the binary's evolution. Contrary to assumptions made in prior modelling, these donors are not fully convective over much of the AM CVn phase and do not evolve adiabatically under mass loss indefinitely. Instead, we identify three distinct phases of evolution: a mass transfer turn-on phase (during which the orbital period continues to decrease even after contact, the donor contracts, and the mass transfer rate accelerates to its maximum), a phase in which the donor expands adiabatically in response to mass loss, and a cooling phase beginning at orbital periods of approximately 45--55 minutes during which the donor contracts. The physics that determines the behaviour in the first and third phases, both of which are new outcomes of this study, are discussed in some detail. We find the overall duration of the turn-on phase to be between $\sim 10^4$-$10^6$ yrs, significantly longer than prior estimates. We predict the donor's luminosity and effective temperature. During the adiabatic expansion phase (ignoring irradiation effects), the luminosity is approximately $10^{-6}$--$10^{-4} L_\odot$ and the effective temperature is approximately 1000--1800 K. However, the flux generated in the accretion flow dominates the donor's intrinsic light at all times. The impact of irradiation on the donor extends the phase of adiabatic expansion to longer orbital periods and alters the donor's observational characteristics. Irradiated donors during the adiabatic phase can attain a surface luminosity of up to $\approx10^{-2} L_\odot$. We argue that the turn-on and cooling phases both will leave significant imprints on the AM CVn population's orbital period distribution.

We discuss results of our study on AM CVn binaries formed with donors that never ignited He before contact. For the first time, we treat the donor's in these systems in the context of a full stellar structure evolution theory and find that the binary's evolution can described in terms of 3 phases: contact, adiabatic donor expansion, and late-time donor cooling. Details of the first and third phase are new results from this study and we focus on generally characterizing these two phases. Finally, we present our predictions for the donor's light in these systems.

We report on evolutionary calculations of the onset of mass transfer in AM CVn binaries, treating the donor's evolution in detail. We show that during the early contact phase, while the mass transfer rate, $\Mdot$, is increasing, gravity wave (GW) emission continues to drive the binary to shorter orbital period, $\Porb$. We argue that the phase where $\Mdot > 0$ and $\nudot > 0$ ($\nu = 1/\Porb$) can last between $10^3$ and $10^6$ yrs, significantly longer than previously estimated. These results are applied to RX J0806+1527 ($\Porb = 321$ s) and RX J914+2456 ($\Porb=569$ s), both of which have measured $\nudot > 0$. \emph{Thus, a $\nudot > 0$ does not select between the unipolar inductor and accretion driven models proposed as the source of X-rays in these systems}. For the accretion model, we predict for RX J0806 that $\ddot{\nu} \approx \ee{1.0-1.5}{-28}$ Hz s$^{-2}$ and argue that timing observations can probe $\ddot{\nu}$ at this level with a total $\approx 20$ yr baseline. We also place constraints on each system's initial parameters given current observational data.

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

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.

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.

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