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 of Z = 0.02, 0.01, 0.001 and 0.0002, producing WDs with masses between ∼0.16−0.45 M . 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. In addition, we confirm that the threshold mass at which this dichotomy occurs depends on metallicity. 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, i.e. the observed presence of metals in their atmospheres. We predict that these proto-WDs have helium-rich envelopes through a significant part of their lifetime. This is of great importance as helium is a crucial ingredient in the driving of the κ-mechanism suggested for the newly observed ELM proto-WD pulsators. However, we find that the number of hydrogen shell flashes and, as a result, the hydrogen envelope mass at the beginning of the cooling track, are not influenced significantly by rotational mixing. In addition to being dependent on proto-WD mass and metallicity, the hydrogen envelope mass of the newly formed proto-WDs depends on whether or not the donor star experiences a temporary contraction when the H-burning shell crosses the hydrogen discontinuity left behind by the convective envelope. 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.

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 aim to decipher the kinds of BH binaries with supergiant companions that could potentially evolve through and 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 binaries is only possible if the donor is a massive convective-envelope giant, namely a red supergiant (RSG). The same is true for neutron-star binaries with massive companions. In other words, pre-CE progenitors of BBH mergers are BH binaries with RSG companions. We find that because of its influence on the radial expansion of massive giants, metallicity has an indirect but a very strong effect on the chemical profile, density structure, and the binding energies of RSG envelopes. 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⊙) ≈ 5.6 − 5.8, corresponding to stars with M ≳ 40 M⊙, 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 total mass ≲50 M⊙. 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.

Context. A large number of low-mass ( < ∼ 0.20 M ) 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 lowmass He WDs and a number of key observed cases, indicating that their formation scenario yet remains to be fully understood. Aims. 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. Methods. Numerical calculations with a detailed stellar evolution code were used to trace the mass-transfer phase in a large number of close-orbit LMXBs with different initial values of donor star mass, neutron star mass, orbital period, and strength of magnetic braking. 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. Results. We find that the time between Roche-lobe detachment until the low-mass proto-He WD reaches the WD cooling track is typically Δ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 for progenitor stars of mass M 2 ≤ 1.5 M (and ∼0.18 M for M 2 = 1.6 M ). Conclusions. 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 Δt proto to depend on the occurrence of flashes and thus question the suggested dichotomy in thermal evolution of proto-WDs.

Context. Millisecond pulsars (MSPs) are generally believed to be old neutron stars (NSs) that have been spun up to high rotation rates via accretion of matter from a companion star in a low-mass X-ray binary (LMXB). This scenario has been strongly supported by various pieces of observational evidence. However, many details of this recycling scenario remain to be understood. Aims. 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 h. In particular, we examine i) if the observed systems can be reproduced by 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. Methods. 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 γ-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. Results. We find that severe fine-tuning is necessary to reproduce the observed MSPs in tight binaries with He WD companions of mass <0.20 M , which suggests that something needs to be modified or is missing in the standard input physics of LMXB modelling. Results from previous independent studies support this conclusion. 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 . The results of the thermal evolution of the (proto) He WDs are reported in a follow-up paper (Paper II).

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. Our goal is to study the impact of this effect in the context of binary evolution. Using the stellar-evolution code MESA, we compute evolutionary tracks of massive stars at six different metallicities between 1.0 Z and 0.01 Z . We explore variations of factors known to affect the radial expansion of massive stars (eg. semiconvection, overshooting, rotation). Using observational constraints, we find support for evolution in which already at metallicity Z ≈ 0.2 Z massive stars stay relatively compact (∼ 100 R ) during the Hertzprung-Gap phase (HG) and most of their expansion happens 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. However, at lower metallicity the parameter space for mass transfer from a more evolved, slowly-expanding CHeB star increases dramatically. This means that envelope stripping and formation of helium stars in low metallicity environments happens later during the evolution of the donor, implying a 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 metallicities Z ≤ 0.04 Z a significant fraction of massive stars in binaries with periods above 100 days engage in their 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. Cases of unstable mass transfer leading to a merger would produce CO cores that are fast spinning at the moment of collapse. We find that the parameter space for mass transfer from massive donors (> 40 M ) with outer convective envelopes is extremely small or even non-existent. We briefly discuss this finding in the context of formation of binary black hole mergers.

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

AM CVn binaries consist of a WD accreting from a hydrogen-deficient star (or WD) companion (Warner, 1995;Solheim, 2010). In their formation history (Fig. 1.6 and Section 1.3.1.1), AM CVns form after at least one CE phase of their progenitor system. The current RLO is initiated, due to orbital damping caused by GW radiation, at orbital periods of typically 5−20 min (depending on the nature and the temperature of the companion star), and the mass-transfer rate is determined

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