The stakes could hardly be higher. Gravitational wave detections of merging black holes and neutron stars have made it abundantly clear that some mechanism must shrink stellar orbits by orders of magnitude before the final, feeble gravitational-wave-driven inspiral can begin. For most systems, that mechanism is the common envelope. The concept dates back to 1976, when Bohdan Paczyński proposed it to explain V 471 Tauri, a compact white dwarf and K dwarf binary whose tight orbit could only be understood if the two stars had once shared an envelope. Since then, the common envelope has become the linchpin of compact binary formation theory, invoked to explain cataclysmic variables, binary pulsars, short gamma-ray burst progenitors, and the close white dwarf pairs that will one day merge as Type Ia supernovae. Practically all massive stars are born in multiple systems, and up to seventy percent of them experience binary interaction during their lives, making this phase not an exotic curiosity but a routine chapter in stellar history.

The challenge for theorists is that the common envelope phase is fast, asymmetric, and wildly multi-scale. The review organizes the interaction into three stages: a pre-envelope phase in which unstable mass transfer begins, a dynamical inspiral in which the companion plunges into the giant’s envelope, and a post-inspiral phase in which the system either ejects the envelope or settles toward merger. The dynamical inspiral, the heart of the problem, unfolds on timescales of roughly fifty days for the envelopes of red giant and asymptotic giant branch stars. Yet the cores of those same stars, with sound speeds exceeding a thousand kilometers per second, evolve on dynamical timescales of just twenty seconds. This five-order-of-magnitude temporal gap, combined with spatial scales spanning up to eight orders of magnitude when companions range from planets to black holes, means that fully resolved three-dimensional simulations remain out of reach for the foreseeable future.

For decades, the field relied on parametric shortcuts. The energy formalism, introduced by van den Heuvel and Webbink, balances the binding energy of the envelope against the orbital energy lost during inspiral, mediated by an efficiency parameter called alpha. Population synthesis models that predict gravitational wave event rates and Type Ia supernova progenitor channels depend critically on this parameter, yet its value remains deeply uncertain. Observational calibrations based on post-common-envelope white dwarf binaries have yielded conflicting results, with some studies favoring a constant alpha and others a value that varies with system parameters. A parallel angular momentum formalism, the so-called gamma prescription, suffers from its own lack of predictive power. Compounding the problem, the very definition of envelope binding energy is ambiguous: whether recombination energy released as ionized hydrogen and helium recombine can be tapped to eject the envelope is one of the most contested questions in the field, and the answer may push the effective efficiency above unity.

One-dimensional mechanical models attempted something more physical, integrating an equation of motion for the inspiralling companion under gravity and an assumed drag force while feeding the released energy into a stellar evolution code. The drag itself traces back to classical work on gravitational focusing: Hoyle and Lyttleton’s treatment of accretion onto a moving point mass, refined by Bondi and Hoyle into the famous accretion column picture, and complemented by Bondi’s spherical limit. Dynamical friction, first analyzed by Chandrasekhar for collisionless systems and extended to gaseous media by Dokuchaev and later Ostriker, produces the backward pull that shrinks the orbit. These analytic expressions, involving the Mach number of the companion’s motion and a Coulomb logarithm whose integration limits remain arbitrary, capture the functional dependencies but not the full nonlinear reality of a turbulent, shock-laden envelope.

Three-dimensional hydrodynamic simulations, begun tentatively in the late 1980s and maturing through waves of increasingly capable codes, have transformed the picture. They revealed just how violently non-spherical the interaction is: the envelope deforms into a toroidal shape, spiral shocks issue from the core binary, and shear instabilities shred the outflow. Crucially, they showed that ejecting the envelope using orbital energy alone is extraordinarily difficult. The inspiral tends to stall while a significant fraction of the envelope remains bound, a phenomenon attributed to several conspiring effects: the expansion and dilution of gas around the cores weakens the drag, the gas can be dragged into co-rotation with the binary until the velocity contrast vanishes, and the growing gravitational attraction between the approaching cores demands ever larger forces to sustain the inspiral. Simulations that include recombination energy through realistic equations of state eject substantially more material, though whether that energy thermalizes or radiates away remains contested, making radiation transport an urgent addition to the modeling toolkit.

The numerical techniques themselves form a fascinating landscape. Smoothed particle hydrodynamics, a Lagrangian particle method, excels at conserving angular momentum and avoiding advection errors, making it naturally suited to two orbiting stars embedded in vacuum, but it struggles to resolve shocks and low-density flows without enormous particle counts. Eulerian grid-based finite-volume schemes capture shocks and instabilities with superior accuracy through Riemann solvers, but suffer advection errors on fixed grids and require adaptive mesh refinement to concentrate resolution where it matters. Moving-mesh codes such as arepo, which evolved from cosmological simulation technology, thread the needle by advecting an unstructured Voronoi mesh with the flow, combining Lagrangian flexibility with Godunov-type accuracy. The first three-dimensional magnetohydrodynamic simulations of common envelope evolution, performed with arepo, showed that magnetic fields are amplified enormously by the magnetorotational instability in the accretion flow around the companion, though they remain dynamically subdominant during the inspiral itself.

Equally delicate is the art of setting up these simulations. The core of a giant star cannot be resolved without triggering fatal timestep restrictions, so modelers replace it with a gravitating point particle and reconstruct the envelope in hydrostatic equilibrium using a modified Lane-Emden equation. Gravitational softening must be applied to prevent the point masses from generating unphysical singular forces, yet the choice of softening length measurably alters the inspiral rate and envelope unbinding. Mapping one-dimensional stellar evolution models onto three-dimensional grids introduces noise and spurious velocities that must be damped away over several dynamical timescales before the companion is released. Even the pseudo-vacuum used to fill empty grid cells can contaminate predictions of observable light curves and colors. Energy and angular momentum conservation, monitored to the ten-percent level in early work, remains a critical diagnostic, because the envelope is so loosely bound that numerical errors can masquerade as physical ejection.

What emerges from the review is a field in rapid ascent. The past decade has seen a proliferation of global simulations, wind-tunnel experiments probing drag forces at high resolution, and growing recognition that the pre-envelope and post-inspiral phases, which unfold on thermal timescales of thousands of years, must eventually be coupled to the dynamical calculations. Observational anchors abound: post-common-envelope binaries provide statistical constraints, planetary nebulae bear the morphological fingerprints of ejected envelopes, and luminous red novae such as V1309 Sco, the merger that was caught in pre-outburst survey data, offer direct glimpses of the interaction in action. The ultimate prize is a predictive connection between the parameters of a binary entering the common envelope phase and the properties of the remnant that emerges, whether a tight double white dwarf, an X-ray binary, or a future gravitational wave source. Röpke and De Marco conclude that with rapidly advancing computational power, refined numerical techniques, and a clearer grasp of the relevant physics, one of the last fundamental unsolved problems of stellar astrophysics may finally yield.

Subject of Research: Three-dimensional hydrodynamic simulations of common-envelope evolution in binary stellar systems

Article Title: Simulations of common-envelope evolution in binary stellar systems: physical models and numerical techniques

Article References: Simulations of common-envelope evolution in binary stellar systems: physical models and numerical techniques. (n.d.). https://doi.org/10.1007/s41115-023-00017-x

Image Credits: AI Generated

DOI: 10.1007/s41115-023-00017-x

Keywords: common envelope evolution, binary stars, stellar astrophysics, hydrodynamic simulations, gravitational waves, compact objects, smoothed particle hydrodynamics, adaptive mesh refinement, recombination energy, stellar mergers, computational astrophysics, binary evolution