skip to main content

Title: Jets from main sequence and white dwarf companions during common envelope evolution

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 more » 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.

« less
; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 3041-3057
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this

    We present small-scale 3D hydrodynamical simulations of the evolution of a 0.3 M⊙ main-sequence (MS) star that launches two perpendicular jets within the envelope of a 0.88 M⊙ red giant (RG). Based on previous large-scale simulations, we study the dynamics of the jets either when the secondary star is grazing, when it has plunged-in, or when it is well within the envelope of the RG (in each stage for ∼11 d). The dynamics of the jets through the common envelope (CE) depend on the conditions of the environment as well as on their powering. In the grazing stage and the commencement of the plunge self-regulated jets need higher efficiencies to break out of the envelope of the RG. Deep inside the CE, on the time-scales simulated, jets are choked independently of whether they are self-regulated or constantly powered. Jets able to break out of the envelope of the RG in large-scale simulations, are choked in our small-scale simulations. The accreted angular momentum on to the secondary star is not large enough to form a disc. The mass accretion on to the MS star is 1–10 per cent of the Bondi–Hoyle–Littleton rate (∼10−3–10−1 M⊙ yr−1). High-luminosity emission, from X-rays to ultraviolet and optical, is expected ifmore »the jets break out of the CE. Our simulations illustrate the need for inclusion of more realistic accretion and jet models in the dynamical evolution of the CEs.

    « less
  2. Abstract Long-duration γ -ray bursts (GRBs) accompany the collapse of massive stars and carry information about the central engine. However, no 3D models have been able to follow these jets from their birth via black hole (BH) to the photosphere. We present the first such 3D general-relativity magnetohydrodynamic simulations, which span over six orders of magnitude in space and time. The collapsing stellar envelope forms an accretion disk, which drags inwardly the magnetic flux that accumulates around the BH, becomes dynamically important, and launches bipolar jets. The jets reach the photosphere at ∼10 12 cm with an opening angle θ j ∼ 6° and a Lorentz factor Γ j ≲ 30, unbinding ≳90% of the star. We find that (i) the disk–jet system spontaneously develops misalignment relative to the BH rotational axis. As a result, the jet wobbles with an angle θ t ∼ 12°, which can naturally explain quiescent times in GRB lightcurves. The effective opening angle for detection θ j + θ t suggests that the intrinsic GRB rate is lower by an order of magnitude than standard estimates. This suggests that successful GRBs are rarer than currently thought and emerge in only ∼0.1% of supernovae Ib/c, implyingmore »that jets are either not launched or choked inside most supernova Ib/c progenitors. (ii) The magnetic energy in the jet decreases due to mixing with the star, resulting in jets with a hybrid composition of magnetic and thermal components at the photosphere, where ∼10% of the gas maintains magnetization σ ≳ 0.1. This indicates that both a photospheric component and reconnection may play a role in the prompt emission.« less
  3. Abstract We compute the forces, torque and rate of work on the companion-core binary due to drag in global simulations of common envelope (CE) evolution for three different companion masses. Our simulations help to delineate regimes when conventional analytic drag force approximations are applicable. During and just prior to the first periastron passage of the in-spiral phase, the drag force is reasonably approximated by conventional analytic theory and peaks at values proportional to the companion mass. Good agreement between global and local 3D “wind tunnel” simulations, including similar net drag force and flow pattern, is obtained for comparable regions of parameter space. However, subsequent to the first periastron passage, the drag force is up to an order of magnitude smaller than theoretical predictions, quasi-steady, and depends only weakly on companion mass. The discrepancy is exacerbated for larger companion mass and when the inter-particle separation reduces to the Bondi-Hoyle-Lyttleton accretion radius, creating a turbulent thermalized region. Greater flow symmetry during this phase leads to near balance of opposing gravitational forces in front of and behind the companion, hence a small net drag. The reduced drag force at late times helps explain why companion-core separations necessary for envelope ejection are not reachedmore »by the end of limited duration CE simulations.« less

    Tidal evolution of eccentric binary systems containing at least one massive main-sequence (MS) star plays an important role in the formation scenarios of merging compact-object binaries. The dominant dissipation mechanism in such systems involves tidal excitation of outgoing internal gravity waves at the convective-radiative boundary and dissipation of the waves at the stellar envelope/surface. We have derived analytical expressions for the tidal torque and tidal energy transfer rate in such binaries for arbitrary orbital eccentricities and stellar rotation rates. These expressions can be used to study the spin and orbital evolution of eccentric binaries containing massive MS stars, such as the progenitors of merging neutron star binaries. Applying our results to the PSR J0045-7319 system, which has a massive B-star companion and an observed, rapidly decaying orbit, we find that for the standard radius of convective core based on non-rotating stellar models, the B-star must have a significant retrograde and differential rotation in order to explain the observed orbital decay rate. Alternatively, we suggest that the convective core may be larger as a result of rapid stellar rotation and/or mass transfer to the B-star in the recent past during the post-MS evolution of the pulsar progenitor.

  5. ABSTRACT Common envelope (CE) evolution, which is crucial in creating short-period binaries and associated astrophysical events, can be constrained by reverse modelling of such binaries’ formation histories. Through analysis of a sample of well-constrained white dwarf (WD) binaries with low-mass primaries (seven eclipsing double WDs, two non-eclipsing double WDs, one WD-brown dwarf), we estimate the CE energy efficiency αCE needed to unbind the hydrogen envelope. We use grids of He- and CO-core WD models to determine the masses and cooling ages that match each primary WD’s radius and temperature. Assuming gravitational wave-driven orbital decay, we then calculate the associated ranges in post-CE orbital period. By mapping WD models to a grid of red giant progenitor stars, we determine the total envelope binding energies and possible orbital periods at the point CE evolution is initiated, thereby constraining αCE. Assuming He-core WDs with progenitors of 0.9–2.0 M⊙, we find αCE ∼ 0.2–0.4 is consistent with each system we model. Significantly higher values of αCE are required for higher mass progenitors and for CO-core WDs, so these scenarios are deemed unlikely. Our values are mostly consistent with previous studies of post-CE WD binaries, and they suggest a nearly constant and low envelope ejection efficiencymore »for CE events that produce He-core WDs.« less