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.
- Award ID(s):
- 1813298
- PAR ID:
- 10163568
- 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
-
ABSTRACT -
Abstract About ten percent of Sun-like (1–2
M ⊙) stars will engulf a 1–10M Jplanet 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 withM p≲ 5M Jplanets, 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 anM p≳ 5M Jplanet 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. -
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}$.
-
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.
-
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.