Common-envelope (CE) evolution in massive binary systems is thought to be one of the most promising channels for the formation of compact binary mergers. In the case of merging binary black holes (BBHs), the essential CE phase takes place at a stage when the first BH is already formed and the companion star expands as a supergiant. We study which BH binaries with supergiant companions will evolve through and potentially survive a CE phase. To this end, we compute envelope binding energies from detailed massive stellar models at different evolutionary stages and metallicities. We make multiple physically extreme choices of assumptions that favor easier CE ejection as well as account for recent advancements in mass transfer stability criteria. We find that even with the most optimistic assumptions, a successful CE ejection in BH (and also NS) binaries is only possible if the donor is a massive convective-envelope giant, a red supergiant (RSG). In other words, pre-CE progenitors of BBH mergers are BH binaries with RSG companions. We find that due to its influence on the radial expansion of massive giants, metallicity has an indirect but a very strong effect on the envelope structure and binding energies of RSGs. Our results suggest that merger rates from population synthesis models could be severely overestimated, especially at low metallicity. Additionally, the lack of observed RSGs with luminosities above log($L/L_{\odot}$) = 5.6-5.8, corresponding to stars with $M > 40 M_{\odot}$, puts into question the viability of the CE channel for the formation of the most massive BBH mergers. Either such RSGs elude detection due to very short lifetimes, or they do not exist and the CE channel can only produce BBH systems with masses $< 50 M_{\odot}$. Finally, we discuss an alternative CE scenario, in which a partial envelope ejection is followed by a phase of possibly long and stable mass transfer.

Metallicity is known to significantly affect the radial expansion of a massive star: the lower the metallicity, the more compact the star, especially during its post-MS evolution. We study this effect in the context of binary evolution. Using the stellar-evolution code MESA, we computed evolutionary tracks of stars at different metallicities, exploring variations of factors known to affect the radial expansion (e.g. semiconvection, overshooting, rotation). We find observational support for an evolution in which already at metallicity $0.2Z_{\odot}$ massive stars remain relatively compact during the Hertzprung-Gap (HG) phase and most of their expansion occurs during core-helium burning (CHeB). Consequently, we show that metallicity has a strong influence on the type of mass transfer evolution in binary systems. At solar metallicity, a case-B mass transfer is initiated shortly after the end of MS, and a giant donor is almost always a rapidly expanding HG star. At lower metallicity, the parameter space for mass transfer from a more evolved CHeB star increases dramatically. This means that envelope stripping and formation of helium stars in low-metallicity environments occurs later in the evolution of the donor, implying a much shorter duration of the Wolf-Rayet phase (even by an order of magnitude) and higher final core masses. This metallicity effect is independent of the impact of metallicity-dependent stellar winds. At very low metallicities, a significant fraction of massive stars in binaries engages in the first episode of mass transfer very late into their evolution, when they already have a well-developed CO core. The remaining lifetime ($< 10^4$ yr) is unlikely to be enough to strip the entire H-rich envelope. We also briefly discuss the extremely small parameter space for mass transfer from massive convective-envelope donors in the context of binary black hole merger formation.

Binaries harbouring millisecond pulsars enable a unique path to determine neutron star masses: radio pulsations reveal the motion of the neutron star, while that of the companion can be characterised through studies in the optical range. PSR J1012+5307 is a millisecond pulsar in a 14.5-h orbit with a helium-core white dwarf companion. In this work we present the analysis of an optical spectroscopic campaign, where the companion star absorption features reveal one of the lightest known white dwarfs. We determine a white dwarf radial velocity semi-amplitude of K_2 = 218.9 +- 2.2 km/s, which combined with that of the pulsar derived from the precise radio timing, yields a mass ratio of q=10.44+- 0.11. We also attempt to infer the white dwarf mass from observational constraints using new binary evolution models for extremely low-mass white dwarfs, but find that they cannot reproduce all observed parameters simultaneously. In particular, we cannot reconcile the radius predicted from binary evolution with the measurement from the photometric analysis (R_WD=0.047+-0.003 Rsun). Our limited understanding of extremely low-mass white dwarf evolution, which results from binary interaction, therefore comes as the main factor limiting the precision with which we can measure the mass of the white dwarf in this system. Our conservative white dwarf mass estimate of M_WD = 0.165 +- 0.015 Msun, along with the mass ratio enables us to infer a pulsar mass of M_NS = 1.72 +- 0.16 Msun. This value is clearly above the canonical 1.4 Msun, therefore adding PSR J1012+5307 to the growing list of massive millisecond pulsars.

White dwarfs are the burnt out cores of Sun-like stars and are the final fate of 97% of all stars in our Galaxy. The internal structure and composition of white dwarfs are hidden by their high gravities, which causes all elements, apart from the lightest ones, to settle out of their atmospheres. The most direct method to probe the inner structure of stars and white dwarfs in detail is via asteroseismology. Here we present the first known pulsating white dwarf in an eclipsing binary system, enabling us to place extremely precise constraints on the mass and radius of the white dwarf from the light curve, independent of the pulsations. This 0.325M$_\odot$ white dwarf --- one member of SDSS J115219.99+024814.4 --- will serve as a powerful benchmark to constrain empirically the core composition of low-mass stellar remnants and investigate the effects of close binary evolution on the internal structure of white dwarfs.

We report on the discovery of the companion star to the millisecond pulsar J1342+2822B in the globular cluster M3. We exploited a combination of near-ultraviolet and optical observations acquired with the Hubble Space Telescope in order to search for the optical counterparts to the known millisecond pulsars in this cluster. At a position in excellent agreement with that of the radio pulsar J1342+2822B (M3B), we have identified a blue and faint object (mF275W \approx 22.45) that, in the color-magnitude diagram of the cluster, is located in the region of He core white dwarfs. From the comparison of the observed magnitudes with theoretical cooling tracks we have estimated the physical properties of the companion star: it has a mass of only 0.19 \pm 0.02 Msun, a surface temperature of 12 \pm 1 x 10^3 K and a cooling age of 1.0\pm0.2 Gyr. Its progenitor was likely a ~ 0.84 M star and the bulk of the mass-transfer activity occurred during the sub-giant branch phase. The companion mass, combined with the pulsar mass function, implies that this system is observed almost edge-on and that the neutron star has a mass of 1.1 \pm 0.3 Msun, in agreement with the typical values measured for recycled neutron stars in these compact binary systems. We have also identified a candidate counterpart to the wide and eccentric binary millisecond pulsar J1342+2822D. It is another white dwarf with a He core and a mass of 0.22 \pm 0.2 Msun, implying that the system is observed at a high inclination angle and hosts a typical NS with a mass of 1.3 \pm 0.3 Msun. At the moment, the large uncertainty on the radio position of this millisecond pulsar prevents us from robustly concluding that the detected star is its optical counterpart.

We use Hubble Space Telescope multicolour photometry of the globular cluster 47 Tucanae to uncover a population of 24 objects with no previous classification that are outliers from the single-star model tracks in the colour-magnitude diagram and yet are likely cluster members. By comparing those sources with evolutionary models and X-ray source catalogues, we were able to show that the majority of those sources are likely binary systems that do not have any X-ray source detected nearby, most possibly formed by a white dwarf and a main-sequence star and a small number of possible double-degenerate systems.

PSR J0636+5128 is a millisecond pulsar in one of the most compact pulsar binaries known, with a 96\,min orbital period. The pulsar mass function suggests a very low-mass companion, similar to that seen in so-called "black widow" binaries. Unlike in most of those, however, no radio eclipses by material driven off from the companion were seen leading to the possibility that the companion was a degenerate remnant of a carbon-oxygen white dwarf. We report the discovery of the optical counterpart of its companion in images taken with the Gemini North and Keck~I telescopes. The companion varies between $r=25$ and $r=23$ on the 96\,min orbital period of the binary, caused by irradiation from the pulsar's energetic wind. We modeled the multi-color lightcurve using parallax constraints from pulsar timing and determine a companion mass of $(1.71\pm0.23)\times 10^{-2}\,M_\odot$, a radius of $(7.6\pm1.4)\times 10^{-2}\,R_\odot$, and a mean density of $54\pm26\,{\rm g\,cm}^{-3}$, all for an assumed neutron star mass of $1.4\,M_\odot$. This makes the companion to PSR J0636+5128 one of the densest of the "black widow" systems. Modeling suggests that the composition is not predominantly hydrogen, perhaps due to an origin in an ultra-compact X-ray binary.

Previous optical and radio observations of the binary millisecond pulsar PSR J1640+2224 have come to inconsistent conclusions about the identity of its companion, with some observations suggesting the companion is a low-mass helium-core (He-core) white dwarf (WD), while others indicate it is most likely a high-mass carbon-oxygen (CO) WD. Binary evolution models predict PSR J1640+2224 most likely formed in a low-mass X-ray binary (LMXB) based on the pulsar's short spin period and long-period, low-eccentricity orbit, in which case its companion should be a He-core WD with mass about $0.35 - 0.39 \, M_\odot$, depending on metallicity. If it is instead a CO WD, that would suggest the system has an unusual formation history. In this paper we present the first astrometric parallax measurement for this system from observations made with the Very Long Baseline Array (VLBA), from which we determine the distance to be $1520^{+170}_{-150}\,\mathrm{pc}$. We use this distance and a reanalysis of archival optical observations originally taken in 1995 with the Wide Field Planetary Camera 2 (WFPC2) on the Hubble Space Telescope (HST) in order to measure the WD's mass. We also incorporate improvements in calibration, extinction model, and WD cooling models. We find that the existing observations are not sufficient to tightly constrain the companion mass, but we conclude the WD mass is $>0.4\,M_\odot$ with $>90\%$ confidence. The limiting factor in our analysis is the low signal-to-noise ratio of the original HST observations.

We present an analysis of the evolutionary and pulsation properties of the extremely low-mass white dwarf precursor (B) component of the double-lined eclipsing system WASP 0247-25. Given that the fundamental parameters of that star have been obtained previously at a unique level of precision, WASP 0247-25B represents the ideal case for testing evolutionary models of this newly-found category of pulsators. Taking into account the known constraints on the mass, orbital period, effective temperature, surface gravity, and atmospheric composition, we present a model that is compatible with these constraints and show pulsation modes that have periods very close to the observed values. Importantly, these modes are predicted excited. Although the overall consistency remains perfectible, the observable properties of WASP 0247-25B are closely reproduced. A key ingredient of our binary evolutionary models is represented by rotational mixing as the main competitor against gravitational settling. Depending on assumptions made about the values of the degree index l for the observed pulsation modes, we found three possible seismic solutions. We discuss two tests, rotational splitting and multicolor photometry, that should readily identify the modes and discriminate between these solutions. However, this will require improved temporal resolution and higher S/N observations than currently available.

In 2012, five high school students involved in the Pulsar Search Collaboratory discovered the millisecond pulsar PSR J1400$-$1431 and initial timing parameters were published in Rosen et al. (2013) a year later. Since then, we have obtained a phase-connected timing solution spanning five years, resolving a significant position discrepancy and measuring $\dot{P}$, proper motion, parallax, and a monotonic slope in dispersion measure over time. Due to PSR J1400$-$1431's proximity and significant proper motion, we use the Shklovskii effect and other priors to determine a 95% confidence interval for PSR J1400$-$1431's distance, $d=270^{+130}_{-80}$ pc. With an improved timing position, we present the first detection of the pulsar's low-mass white dwarf (WD) companion using the Goodman Spectrograph on the 4.1-m SOAR telescope. Deeper imaging suggests that it is a cool DA-type WD with $T_{\rm eff}=3000\pm100$ K and $R/R_\odot=(2.19\pm0.03)\times10^{-2}\,(d/270\,{\rm pc})$. We show a convincing association between PSR J1400$-$1431 and a $\gamma$-ray point source, 3FGL J1400.5$-$1437, but only weak (3.3-$\sigma$) evidence of pulsations after folding $\gamma$-ray photons using our radio timing model. We detect an X-ray counterpart with XMM-Newton but the measured X-ray luminosity ($1\times10^{29}$ ergs s$^{-1}$) makes PSR J1400$-$1431 the least X-ray luminous rotation-powered millisecond pulsar (MSP) detected to date. Together, our findings present a consistent picture of a nearby ($d\approx230$ pc) MSP in a 9.5-day orbit around a cool, $\sim$0.3 M$_\odot$ WD companion, with orbital inclination, $i\gtrsim60^\circ$.

Tight binaries of helium white dwarfs (He WDs) orbiting millisecond pulsars (MSPs) will eventually "merge" due to gravitational damping of the orbit. The outcome has been predicted to be the production of long-lived ultra-compact X-ray binaries (UCXBs), in which the WD transfers material to the accreting neutron star (NS). Here we present complete numerical computations, for the first time, of such stable mass transfer from a He WD to a NS. We have calculated a number of complete binary stellar evolution tracks, starting from pre-LMXB systems, and evolved these to detached MSP+WD systems and further on to UCXBs. The minimum orbital period is found to be as short as 5.6 minutes. We followed the subsequent widening of the systems until the donor stars become planets with a mass of ~0.005 Msun after roughly a Hubble time. Our models are able to explain the properties of observed UCXBs with high helium abundances and we can identify these sources on the ascending or descending branch in a diagram displaying mass-transfer rate vs. orbital period.

We exploit the recent discovery of pulsations in mixed-atmosphere (He/H), extremely low-mass white dwarf precursors (ELM proto-WDs) to test the proposition that rotational mixing is a fundamental process in the formation and evolution of low-mass helium core white dwarfs. Rotational mixing has been shown to be a mechanism able to compete efficiently against gravitational settling, thus accounting naturally for the presence of He, as well as traces of metals such as Mg and Ca, typically found in the atmospheres of ELM proto-WDs. Here we investigate whether rotational mixing can maintain a sufficient amount of He in the deeper driving region of the star, such that it can fuel, through HeII-HeIII ionization, the observed pulsations in this type of stars. Using state-of-the-art evolutionary models computed with MESA, we show that rotational mixing can indeed explain qualitatively the very existence and general properties of the known pulsating, mixed-atmosphere ELM proto-WDs. Moreover, such objects are very likely to pulsate again during their final WD cooling phase.

A large number of extremely low-mass helium white dwarfs (ELM WDs) have been discovered in recent years. The majority of them are found in close binary systems suggesting they are formed either through a common-envelope phase or via stable mass transfer in a low-mass X-ray binary (LMXB) or a cataclysmic variable (CV) system. Here, we investigate the formation of these objects through the LMXB channel with emphasis on the proto-WD evolution in environments with different metallicities. We study, for the first time, the combined effects of rotational mixing and element diffusion (e.g. gravitational settling, thermal and chemical diffusion) on the evolution of proto-WDs and on the cooling properties of the resulting WDs. We present state-of-the-art binary stellar evolution models computed with MESA for metallicities between Z=0.0002 and Z=0.02, producing WDs with masses between 0.16-0.45 M$_{\odot}$. Our results confirm that element diffusion plays a significant role in the evolution of proto-WDs that experience hydrogen shell flashes. The occurrence of these flashes produces a clear dichotomy in the cooling timescales of ELM WDs, which has important consequences e.g. for the age determination of binary millisecond pulsars. Rotational mixing is found to counteract the effect of gravitational settling in the surface layers of young, bloated ELM proto-WDs and therefore plays a key role in determining their surface chemical abundances. We predict that these proto-WDs have helium-rich envelopes through a significant part of their lifetime, a crucial ingredient for understanding the newly observed ELM proto-WD pulsators. The hydrogen envelope at detachment, although small compared to the total mass of the WD, contains enough angular momentum such that the spin frequency of the resulting WD on the cooling track is well above the orbital frequency.

A significant fraction of the envelope of low- and intermediate-mass stars is unstable to convection, leading to sub-surface turbulent motion. Here, we consider and include the effects of turbulence pressure in our stellar evolution calculations. In search of an observational signature, we compare the fractional contribution of turbulent pressure to the observed macroturbulent velocities in stars at different evolutionary stages. We find a strong correlation between the two quantities, similar to what was previously found for massive OB stars. We therefore argue that turbulent pressure fluctuations of finite amplitude may excite high-order, high-angular degree stellar oscillations, which manifest themselves at the surface an additional broadening of the spectral lines, i.e., macroturbulence, across most of the HR diagram. When considering the locations in the HR diagram where we expect high-order oscillations to be excited by stochastic turbulent pressure fluctuations, we find a close match with the observational $\gamma$ Doradus instability strip, which indeed contains high-order, non-radial pulsators. We suggest that turbulent pressure fluctuations on a percentual level may contribute to the $\gamma$ Dor phenomenon, calling for more detailed theoretical modelling in this direction.

Gaia will identify several 1e5 white dwarfs, most of which will be in the solar neighborhood at distances of a few hundred parsecs. Ground-based optical follow-up spectroscopy of this sample of stellar remnants is essential to unlock the enormous scientific potential it holds for our understanding of stellar evolution, and the Galactic formation history of both stars and planets.

We present Hubble Space Telescope UV spectra of the 4.6 h period double white dwarf SDSS J125733.63+542850.5. Combined with Sloan Digital Sky Survey optical data, these reveal that the massive white dwarf (secondary) has an effective temperature T2 = 13030 +/- 70 +/- 150 K and a surface gravity log g2 = 8.73 +/- 0.05 +/- 0.05 (statistical and systematic uncertainties respectively), leading to a mass of M2 = 1.06 Msun. The temperature of the extremely low-mass white dwarf (primary) is substantially lower at T1 = 6400 +/- 37 +/- 50 K, while its surface gravity is poorly constrained by the data. The relative flux contribution of the two white dwarfs across the spectrum provides a radius ratio of R1/R2 = 4.2, which, together with evolutionary models, allows us to calculate the cooling ages. The secondary massive white dwarf has a cooling age of about 1 Gyr, while that of the primary low-mass white dwarf is likely to be much longer, possibly larger than 5 Gyrs, depending on its mass and the strength of chemical diffusion. These results unexpectedly suggest that the low-mass white dwarf formed long before the massive white dwarf, a puzzling discovery which poses a paradox for binary evolution.

Recently, a large number of low-mass (<0.30 M_sun) helium white dwarfs (He WDs) have been discovered as a result of several surveys campaigns such as WASP, ELM, Kepler or SDSS. The far majority of them have as companion another compact object. There appears to be discrepancies between current theoretical modelling of low-mass He WDs and a number of key observational cases indicating that some details of their formation scenario yet remain to be understood.

Millisecond pulsars (MSPs) are generally believed to be old neutron stars (NSs) which have been spun up to high rotation rates via accretion of matter from a companion star in a low-mass X-ray binary (LMXB). However, many details of this recycling scenario remain to be understood. Here we investigate binary evolution in close LMXBs to study the formation of radio MSPs with low-mass helium white dwarf companions (He WDs) in tight binaries with orbital periods P_orb = 2-9 hr. In particular, we examine: i) if such observed systems can be reproduced from theoretical modelling using standard prescriptions of orbital angular momentum losses (i.e. with respect to the nature and the strength of magnetic braking), ii) if our computations of the Roche-lobe detachments can match the observed orbital periods, and iii) if the correlation between WD mass and orbital period (M_WD, P_orb) is valid for systems with P_orb < 2 days. Numerical calculations with a detailed stellar evolution code were used to trace the mass-transfer phase in ~ 400 close LMXB systems with different initial values of donor star mass, NS mass, orbital period and the so-called gamma-index of magnetic braking. Subsequently, we followed the orbital and the interior evolution of the detached low-mass (proto) He WDs, including stages with residual shell hydrogen burning. We find that a severe fine-tuning is necessary to reproduce the observed MSPs in tight binaries with He WD companions of mass < 0.20 M_sun, which suggests that something needs to be modified or is missing in the standard input physics of LMXB modelling. We demonstrate that the theoretically calculated (M_WD, P_orb)-relation is in general also valid for systems with P_orb < 2 days, although with a large scatter in He WD masses between 0.15-0.20 M_sun. The results of the thermal evolution of the (proto) He WDs are reported in a follow-up paper (Paper II).

A large number of low-mass (< 0.20 M_sun) helium white dwarfs (He WDs) have recently been discovered. The majority of these are orbiting another WD or a millisecond pulsar (MSP) in a close binary system; a few examples are found to show pulsations or to have a main-sequence star companion. There appear to be discrepancies between the current theoretical modelling of such low-mass He WDs and a number of key observed cases, indicating that their formation scenario remains to be fully understood. Here we investigate the formation of detached proto-He WDs in close-orbit low-mass X-ray binaries (LMXBs). Our prime focus is to examine the thermal evolution and the contraction phase towards the WD cooling track and investigate how this evolution depends on the WD mass. Our calculations are then compared to the most recent observational data. Numerical calculations with a detailed stellar evolution code were used to trace the mass-transfer phase in a large number of close-orbit LMXBs. Subsequently, we followed the evolution of the detached low-mass proto-He WDs, including stages with residual shell hydrogen burning and vigorous flashes caused by unstable CNO burning. We find that the time between Roche-lobe detachment until the low-mass proto-He WD reaches the WD cooling track is typically Delta_t_proto = 0.5 - 2 Gyr, depending systematically on the WD mass and therefore on its luminosity. The lowest WD mass for developing shell flashes is ~0.21 M_sun for progenitor stars of mass M2 <= 1.5 M_sun (and ~0.18 M_sun for M2 = 1.6 M_sun). The long timescale of low-mass proto-He WD evolution can explain a number of recent observations, including some MSP systems hosting He WD companions with very low surface gravities and high effective temperatures. We find no evidence for Delta_t_proto to depend on the occurrence of flashes and thus question the suggested dichotomy in thermal evolution of proto-WDs.