Common Envelope Evolution Research Articles

High-velocity string of knots in the outburst of the planetary nebula Hb4

The bipolar collimated outflows of the Hb4 Planetary Nebula (PN) exhibit an evident decrease in their expansion velocity with respect to the distance from the central star. So far, similar velocity law has also been found in Herbig-Haro objects. The interpretation of this peculiar velocity law and the classification of the outflows is the main focal point of this paper. High dispersion long-slit echelle spectra along with high resolution images from Hubble Space Telescope (HST) are applied in the astronomical code SHAPE in order to reproduce a three-dimensional morpho-kinematical model for the core and the bipolar outflows. Its central part shows a number of low-ionization filamentary structures (knots and jets) indicative of common-envelope PNe evolution and it is reconstructed assuming a toroidal structure. The high-resolution HST [N II] image of Hb4 unveils the fragmented structure of outflows. The northern and southern outflows are composed by four and three knots, respectively, and each knot moves outwards with its own expansion velocity. They are reconstructed as string of knots rather than jets.This string of knots is formed by ejection events repeated every 200- 250 years. Hb4 displays several indirect evidence for a binary central system with a [WR] companion evolved through the common envelopes channel.The observed deceleration of the knots is likely associated with shock collisions between the knots and the interstellar medium or nebular material.

Constraints on the common-envelope evolution process from wide triple systems

Common envelope (CE) is an important phase in the evolution of interacting evolved binary systems. The interaction of the binary components during the CE evolution (CEE) stage gives rise to orbital inspiral and the formation of a short-period binary or a merger, on the expense of extending and/or ejecting the envelope. CEE is not well understood, as hydrodynamical simulations show that only a fraction of the CE-mass is ejected during the dynamical inspiral, in contrast with observations of post-CE binaries. Different CE models suggest different timescales are involved in the CE-ejection, and hence a measurement of the CE-ejection timescale could provide direct constraints on the CEE-process. Here we propose a novel method for constraining the mass-loss timescale of the CE, using post-CE binaries which are part of wide-orbit triple systems. The orbit/existence of a third companion constrains the CE mass-loss timescale, since rapid CE mass-loss may disrupt the triple system, while slower CE mass-loss may change the orbit of the third companion without disrupting it. As first test-cases we examine two observed post-CE binaries in wide triples, Wolf-1130 and GD-319. We follow their evolution due to mass-loss using analytic and numerical tools, and consider different mass-loss functions. We calculate a wide grid of binary parameters and mass-loss timescales in order to determine the most probable mass-loss timescale leading to the observed properties of the systems. We find that mass-loss timescales of the order of $10^{3}-10^{5}{\rm yr}$ are the most likely to explain these systems. Such long timescales are in tension with most of the CE mass-loss models, which predict shorter, dynamical timescales, but are potentially consistent with the longer timescales expected from the dust-driven winds model for CE ejection.

Evidence for bimodal orbital separations of white dwarf–red dwarf binary stars

We present the results of a radial velocity survey of 20 white dwarf plus M dwarf binaries selected as a follow up to a \textit{Hubble Space Telescope} study that aimed to spatially resolve suspected binaries. Our candidates are taken from the list of targets that were spatially unresolved with \textit{Hubble}. We have determined the orbital periods for 16 of these compact binary candidates. The period distribution ranges from 0.14 to 9.16\,d and peaks near 0.6\,d. The original sample therefore contains two sets of binaries, wide orbits ($\approx100-1000$\,au) and close orbits ($\lesssim1-10$\,au), with no systems found in the $\approx10-100$\,au range. This observational evidence confirms the bimodal distribution predicted by population models and is also similar to results obtained in previous studies. We find no binary periods in the months to years range, supporting the post common envelope evolution scenario. One of our targets, WD\,1504+546, was discovered to be an eclipsing binary with a period of 0.93\,d.

Signs of accretion in the white dwarf + brown dwarf binary NLTT5306

We present new XSHOOTER spectra of NLTT5306, a 0.44 $\pm$ 0.04\msun white dwarf in a short period (101\,min) binary system with a brown dwarf companion that is likely to have previously undergone common envelope evolution. We have confirmed the presence of H$\alpha$ emission and discovered Na I absorption associated with the white dwarf. These observations are indicative of accretion. Accretion is typically evidenced by high energy emission in the UV and X-ray regime. However our \textit{Swift} observations covering the full orbital period in three wavebands (uvw1, uvm2, uvw2) revealed no UV excess or modulation. We used the X-ray non-detection to put an upper limit on the accretion rate of 2$\times$10$^{-15}$\msun yr$^{-1}$. We compare NLTT5306 to similar accreting binaries with brown dwarf donors and suggest the inferred accretion rate could be from wind accretion or accretion from a debris/dust disk. The lack of evidence for a disk implies NLTT5306 is magnetically funnelling a weak wind from a potentially low gravity brown dwarf. The upper limit on the accretion rate suggests a magnetic field as low as 0.45\,kG would be sufficient to achieve this. If confirmed this would constitute the first detection of a brown dwarf wind and could provide useful constraints on mass loss rates.

On the formation of neutron stars via accretion-induced collapse in binaries

We investigate evolutionary pathways leading to neutron star formation through the collapse of oxygen-neon white dwarf (ONe WD) stars in interacting binaries. We consider (1) non-dynamical mass transfer where an ONe WD approaches the Chandrasekhar mass leading to accretion-induced collapse (AIC) and (2) dynamical timescale merger-induced collapse (MIC) between an ONe WD and another WD. We present rates, delay times, and progenitor properties for two different treatments of common envelope evolution. We show that AIC neutron stars are formed via many different channels and the most dominant channel depends on the adopted common envelope physics. Most AIC and MIC neutron stars are born shortly after star formation, though some have delay times >10 Gyr. The shortest delay time (25-50 Myr) AIC neutron stars have stripped-envelope, compact, helium-burning star donors, though many prompt AIC neutron stars form via wind-accretion from an asymptotic giant branch star. The longest delay time AIC neutron stars, which may be observed as young milli-second pulsars among globular clusters, have a red giant or main sequence donor at the time of NS formation and will eventually evolve into NS + helium WD binaries. We discuss AIC & MIC binaries as potential gravitational wave sources for LISA. Neutron stars created via AIC undergo a LMXB phase, offering an electromagnetic counterpart for those shortest orbital period sources that LISA could identify. The formation of neutron stars from interacting WDs in binaries is likely to be a key mechanism for the production of LIGO/Virgo gravitational wave sources (NS-NS and BH-NS mergers) in globular clusters.

On the requirement for a high resolution spectrograph to investigate the multiplicity of massive stars with the 3.6-m Devasthal Optical Telescope

The investigation of the multiplicity of massive stars constitutes a key aspect of their understanding. On the one hand, the evolution of short period systems can be tremendously influenced by the companion, for instance through mass exchange or common envelope evolution. On the other hand, the wind-wind interaction region in massive multiple systems is the seat of a high level physics, including strong thermal X-ray emission, particle acceleration and non-thermal emission processes. The description and understanding of these processes require the accurate determination of orbital parameters, notably through spectroscopic studies. In this context, it is timely to anticipate some relevant applications of a high resolution spectrograph (HRS) to be mounted on the Devasthal Optical Telescope, emphasizing its potential complementarity with other Indian facilities such as the GMRT and Astrosat.

Variable jets at the termination of the common envelope evolution

I propose that at the termination of the common envelope evolution (CEE) the companion to the giant star might launch jets that have variable directions and intensities, hence the jets shape the inner zones of the descendant nebula causing it to lose any type of symmetry. This might account for some of the chaotic departures from axisymmetric structures in some planetary nebulae. I base my study on the assumption that at the termination of the CEE the binary interaction forms a circumbinary thick disk that feeds the companion via an accretion disk that launches opposite jets. I consider two processes that might lead to the development of Rayleigh-Taylor instabilities that lead to chaotic accretion flow on to the companion. The chaotic accretion flow leads to a variable accretion disk around the companion, hence to variable jets. In the first process the jets inflate low-density hot bubbles in the inner regions of the envelope. In the second process the inner boundary of the thick circumbinary disk, that faces the center, cools by radiation thus leading to a temperature, and hence also a pressure, increase with radius. The opposite pressure and density gradients lead to Rayleigh-Taylor instabilities. This study further underlines the importance of exploring the termination phase of the CEE to the understanding of the entire CEE and its outcomes.

On the Role of Supernova Kicks in the Formation of Galactic Double Neutron Star Systems

Abstract In this work we focus on a group of Galactic double neutron star (DNS) systems with long orbital periods of ≳1 day and low eccentricities of ≲0.4. The feature of these orbital parameters is used to constrain the evolutionary processes of progenitor binaries and the supernova (SN) kicks of the second born NSs. Adopting that the mass transfer during primordial binary evolution is highly nonconservative (rotation-dependent), the formation of DNS systems involves a double helium star binary phase, the common envelope (CE) evolution initiates before the first NS formation. During the CE evolution the binary orbital energy is obviously larger when using a helium star rather than an NS to expel the donor envelope, this can help explain the formation of DNS systems with long periods. SN kicks at NS birth can lead to eccentric orbits and even the disruption of binary systems, and the low eccentricities require that the DNSs receive a small natal kick at the second collapse. Compared with the overall distribution of orbital parameters for observed DNS binaries, we propose that the second born NSs in most DNS systems are subject to small natal kicks with the Maxwellian dispersion velocity of less than 80 km s−1, which can provide some constraints on the SN explosion processes. The mass distribution of DNS binaries is also briefly discussed. We suggest that the rotation-dependent mass transfer mode and our results about SN kicks should be applied to massive binary evolution and population synthesis studies.

The demographics of neutron star – white dwarf mergers

Context. The mergers of neutron stars (NSs) and white dwarfs (WDs) could give rise to explosive transients, potentially observable with current and future transient surveys. However, the expected properties and distribution of such events is not well understood. Aims. Here we characterise the rates of such events, their delay-time distributions, their progenitors, and the distribution of their properties. Methods. We use binary population synthesis models and consider a wide range of initial conditions and physical processes. In particular we consider different common-envelope evolution models and different NS natal kick distributions. We provide detailed predictions arising from each of the models considered. Results. We find that the majority of NS–WD mergers are born in systems in which mass-transfer played an important role, and the WD formed before the NS. For the majority of the mergers the WDs have a carbon-oxygen composition (60−80%) and most of the rest are with oxygen-neon WDs. The time-integrated rates of NS–WD mergers are in the range of 3−15% of the type Ia supernovae (SNe) rate. Their delay-time distribution is very similar to that of type Ia SNe, but is slightly biased towards earlier times. They typically explode in young 100 Myr < τ < 1 Gyr environments, but have a tail distribution extending to long, gigayear-timescales. Models including significant kicks give rise to relatively wide offset distribution extending to hundreds of kiloparsecs. Conclusions. The demographic and physical properties of NS–WD mergers suggest they are likely to be peculiar type Ic-like SNe, mostly exploding in late-type galaxies. Their overall properties could be related to a class of recently observed rapidly evolving SNe, while they are less likely to be related to the class of Ca-rich SNe.

Planetary Nebulae Shaped by Common Envelope Evolution

The morphologies of planetary nebula have long been believed to be due to wind shaping processes in which a “fast wind” from the central star impacts a previously ejected envelope. It is assumed that asymmetries existing in the “slow wind” envelope would lead to inertial confinement, shaping the resulting interacting wind flow. We present new results demonstrating the effectiveness of Common Envelope Evolution (CEE) at producing aspherical envelopes which, when impinged upon by a spherical fast stellar wind, produce highly bipolar, jet-like outflows. We have run two simple cases using the output of a single PHANTOM SPH CEE simulation. Our work uses the Adaptive Mesh Refinement code AstroBEAR to track the interaction of the fast wind and CEE ejecta allows us to follow the morphological evolution of the outflow lobes at high resolution in 3-D. Our two models bracket low and high momentum output fast winds. We find the interaction leads to highly collimated bipolar outflows. In addition, the bipolar morphology depends on the fast wind momentum injection rate. With this dependence comes the initiation of significant symmetry breaking between the top and bottom bipolar lobes. Our simulations, though simplified, confirm the long-standing belief that CEE can plan a major role in PPN and PN shaping. These simulations are intended as an initial exploration of the post-CE/PPN flow patterns that can be expected from central source outflows and CE ejecta.

Formation of hot subdwarf B stars with neutron star components

Context. Binary population synthesis predicts the existence of subdwarf B stars (sdBs) with neutron star (NS) or black hole (BH) companions. Several works have been dedicated to finding such systems, but none has been confirmed yet. Theoretically, the formation of sdBs with white dwarf (WD) and main sequence (MS) companions has been well investigated, while those with NS or BH companions remain to be explored further. Aims. We systematically investigate the formation of sdB+NS binaries from binary evolution and aim to obtain some clues for a search for such systems. Methods. We started from a series of MS+NS systems and determined the parameter spaces for producing sdB+NS binaries from the stable Roche-lobe overflow (RLOF) channel and from the common envelope (CE) ejection channel. The parameters for sdB+NS binaries were obtained from detailed binary evolution calculation with the code called modules for experiments in stellar astrophysics (MESA), and the CE parameters were given by the standard energy budget for CE evolution. The MS star had an initial mass ranging from 0.8 to 5 M⊙. Various NS accretion efficiencies and NS masses were examined to investigate the effects they have. We show the characteristics of the produced sdB+NS systems, such as the mass of components, orbital period, the semi-amplitude of the radial velocity (K), and the spin of the NS component. Results. sdB+NS binaries can be produced either from stable RLOF or from CE ejection. In the stable RLOF channel, sdBs can be formed when the donor starts mass transfer close to the tip of the giant branch if the donor has an initial mass ≤2.0 M⊙. For more massive donors, sdBs can be formed when the donor starts mass transfer during the Hertzsprung gap or near the end of the MS. The orbital period of sdB+NS binaries produced in this way ranges from several days to more than 1000 days and moves toward the short-period (∼hr) side with increasing initial MS mass. The highest K is about 150 km s−1 for an MS star of initially 5 M⊙. However, the sdB+NS systems that result from CE ejection have very short orbital periods and then high values of K (up to 800 km s−1). Such systems are born in very young populations (younger than 0.3 Gyr) and are potential gravitational wave sources that might be resolved by the Laser Interferometer Space Antenna (LISA) in the future. Gravitational wave radiation may again bring them into contact on a timescale of only ∼Myr. As a consequence, they are rare and hard to discover. The pulsar signal is likely a feature of sdB+NS systems caused by stable RLOF, and some NS components in sdB binaries may be millisecond pulsars. Various NS accretion efficiencies and NS masses change some properties of sdB+NS binaries, but not our general results.

Origin of magnetic fields in cataclysmic variables

In a series of recent papers, it has been proposed that high field magnetic white dwarfs are the result of close binary interaction and merging. Population synthesis calculations have shown that the origin of isolated highly magnetic white dwarfs is consistent with the stellar merging hypothesis. In this picture, the observed fields are caused by an alpha-Omega dynamo driven by differential rotation. The strongest fields arise when the differential rotation equals the critical break-up velocity and result from the merging of two stars (one of which has a degenerate core) during common envelope evolution or from the merging of two white dwarfs. We now synthesise a population of binary systems to investigate the hypothesis that the magnetic fields in the magnetic cataclysmic variables also originate during stellar interaction in the common envelope phase. Those systems that emerge from common envelope more tightly bound form the cataclysmic variables with the strongest magnetic fields. We vary the common envelope efficiency parameter and compare the results of our population syntheses with observations of magnetic cataclysmic variables. We find that common envelope interaction can explain the observed characteristics of these magnetic systems if the envelope ejection efficiency is low.

The scatter of the M dwarf mass–radius relationship

M dwarfs are prime targets in the hunt for habitable worlds around other stars. This is due to their abundance as well as their small radii and low masses and temperatures, which facilitate the detection of temperate, rocky planets in orbit around them. However, the fundamental properties of M dwarfs are difficult to constrain, often limiting our ability to characterise the planets they host. Here we test several theoretical relationships for M dwarfs by measuring 23 high precision, model-independent masses and radii for M dwarfs in binaries with white dwarfs. We find a large scatter in the radii of these low-mass stars, with 25 per cent having radii consistent with theoretical models while the rest are up to 12 per cent over-inflated. This scatter is seen in both partially- and fully-convective M dwarfs. No clear trend is seen between the over-inflation and age or metallicity, but there are indications that the radii of slowly rotating M dwarfs are more consistent with predictions, albeit with a similar amount of scatter in the measurements compared to more rapidly rotating M dwarfs. The sample of M dwarfs in close binaries with white dwarfs appears indistinguishable from other M dwarf samples, implying that common envelope evolution has a negligible impact on their structure. We conclude that theoretical and empirical mass-radius relationships lack the precision and accuracy required to measure the fundamental parameters of M dwarfs well enough to determine the internal structure and bulk composition of the planets they host.

Radiating the Hydrogen Recombination Energy during Common Envelope Evolution

Abstract Using the stellar evolution code Modules for Experiments in Stellar Astrophysics (MESA), we show that most of the hydrogen recombination energy that is released as the envelope expands during a regular common envelope evolution—namely, the initial dynamical phase or plunge-in phase—is radiated, and hence substantially increases the stellar luminosity. Only about 10% of the hydrogen recombination energy might be used to remove the envelope. We show that the key property of energy transport is that when convection becomes inefficient in the outer parts of the envelope, where the ionization degree of hydrogen falls below about 30%, photon diffusion becomes very efficient and removes the recombination energy. The expanding envelope absorbs most of the gravitational energy that is released by the spiraling-in process of the secondary star inside the common envelope, and so it is the hydrogen recombination energy that is responsible for most of the luminosity increase of the system. The recombination energy of hydrogen adds only a small fraction of the energy required to remove the common envelope, and hence does not play a significant role in the ejection of the envelope.

Are there pre-main-sequence/black hole X-ray binaries?

A large fraction of black-hole low-mass X-ray binaries (BHLMXBs) are in short period orbits, which require strong orbital angular momentum loss during the previous evolutionary stage. Ivanova put forward the possibility that some of the BHLMXBs may harbour pre-main-sequence (pre-MS) donor stars, in order to explain the Li-overabundances in three BHLMXBs. In this work we investigate the evolution of low-mass pre-MS stars in binaries with a BH companion. We calculate the evolution of the spin and orbital periods, the stellar radius and the Roche lobe radius during the pre-MS stage with the stellar evolution code mesa. We find that, because of the relatively slow rotation of the pre-MS star after the common envelope evolution and the long turnover time in the pre-MS stars, tidal torque is not always able to synchronize the pre-MS star, so the spin periods are generally longer than the orbital periods. Mass transfer can occur only for stars with mass larger than ∼1.2 M , which experience expansion due to nuclear reaction after the Hayashi contraction phase. The effective temperatures and orbital periods of these systems do not match the observations of BHLMXBs. Our results show that the observed BHLMXBs with lithium overabundances are unlikely to host pre-MS donor stars.

Oxygen-neon-rich merger during common envelope evolution

We conduct a population synthesis study of the common envelope evolution (CEE) of a white dwarf (WD) and an asymptotic giant branch (AGB) star and find that the potential number of type Ia supernovae (SNe Ia) from the core degenerate (CD) and from the double degenerate (DD) scenarios are of the same order of magnitude. For the CD scenario we consider cases where a carbon oxygen rich (CO) AGB core and a CO WD merge during the CEE and leave a WD remnant with a mass of >1.35Mo, and for the DD scenario we count surviving CO WD binary systems that merge within a time of 10 Gyr. When either the AGB core or the WD are oxygen neon rich (ONe) we assume that the outcome might be a peculiar SN Ia. We find that the number of potential peculiar SNe Ia in the channels we study, that do not include peculiar SNe Ia that involve helium WDs or helium donors, is non-negligible, but less than the number of normal SNe Ia. If a SN Ia or a peculiar SN Ia explosion takes place within about million year after CEE, whether in the CD scenario or in the DD scenario, a massive circumstellar matter is present at explosion time. Our results are compatible with the suggestion that Chandrasekhar-mass SNe Ia mostly come from the CD scenario, and sub-Chandrasekhar SNe Ia mostly come from the DD scenario.

Surveying Planetary Nebulae Central Stars for Close Binaries: Constraining Evolution of Central Stars Based on Binary Parameters

The increase in discovered close binary central stars of planetary nebulae is leading to a sufficiently large sample to begin to make broader conclusions about the effect of close binary stars on common envelope evolution and planetary nebula formation. Herein I review some of the recent results and conclusions specifically relating close binary central stars to nebular shaping, common envelope evolution off the red giant branch, and the total binary fraction and double degenerate fraction of central stars. Finally, I use parameters of known binary central stars to explore the relationship between the proto-planetary nebula and planetary nebula stages, demonstrating that the known proto-planetary nebulae are not the precursors of planetary nebulae with close binary central stars.

Accretion in common envelope evolution

AbstractCommon envelope evolution (CEE) occurs in some binary systems involving asymptotic giant branch (AGB) or red giant branch (RGB) stars, and understanding this process is crucial for understanding the origins of various transient phenomena. CEE has been shown to be highly asymmetrical and global 3D simulations are needed to help understand the dynamics. We perform and analyze hydrodynamic CEE simulations with the adaptive mesh refinement (AMR) code AstroBEAR, and focus on the role of accretion onto the companion star. We bracket the range of accretion rates by comparing a model that removes mass and pressure using a subgrid accretion prescription with one that does not. Provided a pressure-release valve, such as a bipolar jet, is available, super-Eddington accretion could be common. Finally, we summarize new results pertaining to the energy budget, and discuss the overall implications relating to the feasibility of unbinding the envelope in CEE simulations.

Constraining the progenitor evolution of GW 150914

AbstractOne of the largest surprises from the LIGO results regarding the first gravitational wave detection (GW 150914) was the fact the black holes (BHs) were “heavy”, of order 30 - 40 Mȯ. The most promising explanation for this obesity is that the BH-BH merger occurred at low metallicity (Z): when the iron (Fe) contents is lower this is expected to result in weaker mass loss during the Wolf-Rayet (WR) phase. We therefore critically evaluate the claims for the reasons of heavy BHs as a function of Z in the literature. Furthermore, weaker stellar winds might lead to more rapid stellar rotation, allowing WR and BH progenitor evolution in a chemically homogeneous manner. However, there is as yet no empirical evidence for more rapid rotation amongst WR stars in the low Z environment of the Magellanic Clouds. Due to the intrinsic challenge of determining WR rotation rates from emission lines, the most promising avenue to constrain rotation-rate distributions amongst various WR subgroups is through the utilisation of their emission lines in polarised light. We thus provide an overview of linear spectro-polarimetry observations of both single and binary WRs in the Galaxy, as well as the Large and Small Magellanic Clouds, at 50% and 20% of solar Z, respectively. Initial results suggest that the route of chemically homogeneous evolution (CHE) through stellar rotation is challenging, whilst the alternative of a post-LBV or common envelope evolution is more likely.

High-Mass X-ray Binaries: progenitors of double compact objects

AbstractA summary is given of the present state of our knowledge of High-Mass X-ray Binaries (HMXBs), their formation and expected future evolution. Among the HMXB-systems that contain neutron stars, only those that have orbital periods upwards of one year will survive the Common-Envelope (CE) evolution that follows the HMXB phase. These systems may produce close double neutron stars with eccentric orbits. The HMXBs that contain black holes do not necessarily evolve into a CE phase. Systems with relatively short orbital periods will evolve by stable Roche-lobe overflow to short-period Wolf-Rayet (WR) X-ray binaries containing a black hole. Two other ways for the formation of WR X-ray binaries with black holes are identified: CE-evolution of wide HMXBs and homogeneous evolution of very close systems. In all three cases, the final product of the WR X-ray binary will be a double black hole or a black hole neutron star binary.