skip to main content


Title: Common envelope evolution on the asymptotic giant branch: unbinding within a decade?
ABSTRACT Common envelope (CE) evolution is a critical but still poorly understood progenitor phase of many high-energy astrophysical phenomena. Although 3D global hydrodynamic CE simulations have become more common in recent years, those involving an asymptotic giant branch (AGB) primary are scarce, due to the high computational cost from the larger dynamical range compared to red giant branch (RGB) primaries. But CE evolution with AGB progenitors is desirable to simulate because such events are the likely progenitors of most bi-polar planetary nebulae (PNe), and prominent observational testing grounds for CE physics. Here we present a high-resolution global simulation of CE evolution involving an AGB primary and 1-$\mathrm{M_\odot }$ secondary, evolved for 20 orbital revolutions. During the last 16 of these orbits, the envelope unbinds at an almost constant rate of about 0.1–$0.2\, \mathrm{M_\odot \, yr^{-1}}$. If this rate were maintained, the envelope would be unbound in less than $10\, {\rm yr}$. The dominant source of this unbinding is consistent with inspiral; we assess the influence of the ambient medium to be subdominant. We compare this run with a previous run that used an RGB phase primary evolved from the same 2-$\mathrm{M_\odot }$ main-sequence star to assess the influence of the evolutionary state of the primary. When scaled appropriately, the two runs are quite similar, but with some important differences.  more » « less
Award ID(s):
1813298
PAR ID:
10163568
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
495
Issue:
4
ISSN:
0035-8711
Page Range / eLocation ID:
4028 to 4039
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    It has long been speculated that jet feedback from accretion on to the companion during a common envelope (CE) event could affect the orbital evolution and envelope unbinding process. We present global 3D hydrodynamical simulations of CE evolution (CEE) that include a jet subgrid model and compare them with an otherwise identical model without a jet. Our binary consists of a 2-M⊙ red giant branch primary and a 1- or 0.5-M⊙ main sequence (MS) or white dwarf (WD) secondary companion modelled as a point particle. We run the simulations for 10 orbits (40 d). Our jet model adds mass at a constant rate $\dot{M}_\mathrm{j}$ of the order of the Eddington rate, with maximum velocity vj of the order of the escape speed, to two spherical sectors with the jet axis perpendicular to the orbital plane. We explore the influence of the jet on orbital evolution, envelope morphology and envelope unbinding, and assess the dependence of the results on the jet mass-loss rate, launch speed, companion mass, opening angle, and accretion rate. In line with our theoretical estimates, jets are choked around the time of first periastron passage and remain choked thereafter. Subsequent to choking, but not before, jets efficiently transfer energy to bound envelope material. This leads to increases in unbound mass of up to $\sim 10{{\ \rm per\ cent}}$, as compared to the simulations without jets. We also estimate the cumulative effects of jets over a full CE phase, finding that jets launched by MS and WD companions are unlikely to dominate envelope unbinding.

     
    more » « less
  2. Abstract

    About ten percent of Sun-like (1–2M) stars will engulf a 1–10MJplanet as they expand during the red giant branch (RGB) or asymptotic giant branch (AGB) phase of their evolution. Once engulfed, these planets experience a strong drag force in the star’s convective envelope and spiral inward, depositing energy and angular momentum. For these mass ratios, the inspiral takes ∼10–102yr (∼102–103orbits); the planet undergoes tidal disruption at a radius of ∼1R. We use the Modules for Experiments in Stellar Astrophysics (MESA) software instrument to track the stellar response to the energy deposition while simultaneously evolving the planetary orbit. For RGB stars, as well as AGB stars withMp≲ 5MJplanets, the star responds quasi-statically but still brightens measurably on a timescale of years. In addition, asteroseismic indicators, such as the frequency spacing or rotational splitting, differ before and after engulfment. For AGB stars, engulfment of anMp≳ 5MJplanet drives supersonic expansion of the envelope, causing a bright, red, dusty eruption similar to a “luminous red nova.” Based on the peak luminosity, color, duration, and expected rate of these events, we suggest that engulfment events on the AGB could be a significant fraction of low-luminosity red novae in the Galaxy. We do not find conditions where the envelope is ejected prior to the planet’s tidal disruption, complicating the interpretation of short-period giant planets orbiting white dwarfs as survivors of common envelope evolution.

     
    more » « less
  3. ABSTRACT

    Binaries consisting of a hot subdwarf star and an accreting white dwarf (WD) are sources of gravitational wave radiation at low frequencies and possible progenitors of Type Ia supernovae if the WD mass is large enough. Here, we report the discovery of the third binary known of this kind: It consists of a hot subdwarf O (sdO) star and a WD with an orbital period of 3.495 h and an orbital shrinkage of 0.1 s in 6 yr. The sdO star overfills its Roche lobe and likely transfers mass to the WD via an accretion disc. From spectroscopy, we obtain an effective temperature of $T_{\mathrm{eff}}=54\, 240\pm 1840$ K and a surface gravity of log g = 4.841 ± 0.108 for the sdO star. From the light curve analysis, we obtain an sdO mass of MsdO = 0.55 M⊙ and a mass ratio of q = MWD/MsdO = 0.738 ± 0.001. Also, we estimate that the disc has a radius of $\sim\!0.41\ \mathrm{R}_\odot$ and a thickness of $\sim\!0.18\ \mathrm{R}_\odot$. The origin of this binary is probably a common envelope ejection channel, where the progenitor of the sdO star is either a red giant branch star or, more likely, an early asymptotic giant branch star; the sdO star will subsequently evolve into a WD and merge with its WD companion, likely resulting in an R Coronae Borealis (R CrB) star. The outstanding feature in the spectrum of this object is strong Ca H&K lines, which are blueshifted by ∼200 km s−1 and likely originate from the recently ejected common envelope, and we estimated that the remnant common envelope (CE) material in the binary system has a density $\sim\!6\times 10^{-10}\ {\rm g\, cm}^{-3}$.

     
    more » « less
  4. ABSTRACT

    We develop a rapid algorithm for the evolution of stable, circular, circumbinary discs suitable for parameter estimation and population synthesis modelling. Our model includes disc mass and angular momentum changes, accretion on to the binary stars, and binary orbital eccentricity pumping. We fit our model to the post-asymptotic giant branch (post-AGB) circumbinary disc around IRAS 08544−4431, finding reasonable agreement despite the simplicity of our model. Our best-fitting disc has a mass of about $0.01\, \mathrm{M}_{\odot }$ and angular momentum $2.7\times 10^{52}\, \mathrm{g}\, \mathrm{cm}^{2}\, \mathrm{s}^{-1}\simeq 9 \,\mathrm{M}_{\odot }\, \mathrm{km}\, \mathrm{s}^{-1}\, \mathrm{au}$, corresponding to 0.0079 and 0.16 of the common-envelope mass and angular momentum, respectively. The best-fitting disc viscosity is αdisc = 5 × 10−3 and our tidal torque algorithm can be constrained such that the inner edge of the disc Rin ∼ 2a. The inner binary eccentricity reaches about 0.13 in our best-fitting model of IRAS 08544−4431, short of the observed 0.22. The circumbinary disc evaporates quickly when the post-AGB star reaches a temperature of $\sim \! 6\times 10^4\, \mathrm{K}$, suggesting that planetismals must form in the disc in about $10^{4}\, \mathrm{yr}$ if secondary planet formation is to occur, while accretion from the disc on to the stars at ∼10 times the inner-edge viscous rate can double the disc lifetime.

     
    more » « less
  5. ABSTRACT

    We conduct a long-time-scale ($5000\,$ d) 3D simulation of a common-envelope event with a $2\, {\rm M}_{\odot }$ red giant and a $1\, {\rm M}_{\odot }$ main-sequence companion, using the moving-mesh hydrodynamic solver manga. Starting with an orbital radius of $52\, \mathrm{ R}_{\odot }$, our binary shrinks to an orbital radius of $5\, \mathrm{ R}_{\odot }$ in $200\,$ d. We show that over a time-scale of about $1500\,$ d, the envelope is completely ejected, while 80 per cent is ejected in about $400\,$ d. The complete ejection of the envelope is solely powered by the orbital energy of the binary, without the need for late-time reheating from recombination or jets. Motivated by recent theoretical and observational results, we also find that the envelope enters a phase of homologous expansion about $550\, \rm d$ after the start of our simulation. We also run a simplified 1D model to show that heating from the central binary in the envelope at late times does not influence the ejection. This homologous expansion of the envelope would likely simplify calculations of the observational implications such as light curves.

     
    more » « less