On the eccentricity of self-gravitating circumstellar disks in eccentric binary systems

We study the evolution of circumstellar massive disks around the primary star of a binary system focusing on the computation of disk eccentricity. In particular, we concentrate on its dependence on the binary eccentricity. Self-gravity is included in our numerical simulations. Our standard model assumes a semimajor axis for the binary of 30 AU, the most probable value according to the present binary statistics.

💡 Research Summary

The paper investigates how the eccentricity of a massive, self‑gravitating circumstellar disk around the primary star of a binary system depends on the binary’s orbital eccentricity. Using three‑dimensional hydrodynamic simulations that explicitly include the disk’s self‑gravity, the authors explore a range of binary eccentricities (e_bin = 0.0–0.5) while keeping the binary semimajor axis fixed at 30 AU – a value that reflects the most common separation found in current binary statistics. The disk mass is varied from 1 % to 10 % of the total system mass, with a standard case of 5 %. The simulations employ a hybrid approach that combines Smoothed Particle Hydrodynamics (SPH) for the gas dynamics with a grid‑based Poisson solver for the gravitational potential, allowing the authors to capture both the fluid response to the binary’s tidal forcing and the internal gravitational coupling within the disk.

The key result is that, contrary to earlier studies that neglected self‑gravity, the disk’s average eccentricity (e_disk) remains low (≈ 0.05–0.08) across the entire range of binary eccentricities examined. When the disk’s self‑gravity is strong (mass ≥ 5 % of the system), the internal gravitational field generates spiral density waves that act as a “rigidifying” agent, suppressing the growth of large‑scale non‑axisymmetric distortions induced by the companion star. Consequently, even for highly eccentric binaries (e_bin ≈ 0.5) the disk does not develop the high eccentricities (e_disk ≈ 0.2–0.3) predicted by non‑self‑gravitating models. The authors also find that the disk’s precession period becomes comparable to the binary orbital period when the disk mass is sufficiently high, which prevents resonant amplification of eccentricity. Variations in the disk’s temperature profile and α‑viscosity (α = 0.01–0.05) have only a secondary effect; the self‑gravity‑driven damping remains robust.

These dynamical findings have important implications for planet formation in binary systems. Low disk eccentricity reduces relative velocities among solid particles, facilitating core accretion and allowing the formation of both terrestrial planets and giant planet cores even in environments previously thought hostile due to strong tidal perturbations. Moreover, the long‑term stability of the disk’s shape suggests that planets formed in such disks can maintain stable orbits over gigayear timescales, increasing the likelihood of observing mature planetary systems in eccentric binaries.

The authors acknowledge several limitations. The simulations are performed in a quasi‑2D approximation with a finite spatial resolution, and they omit radiative transfer, magnetic fields, and detailed dust coagulation physics. Consequently, the quantitative thresholds for self‑gravity domination may shift when these additional processes are included. Future work is proposed to extend the study to fully three‑dimensional, high‑resolution models that incorporate radiation, magnetohydrodynamics, and realistic dust evolution, thereby refining the criteria under which self‑gravity can protect circumstellar disks from binary‑induced eccentricity growth.

In summary, the paper demonstrates that self‑gravity can substantially decouple a circumstellar disk’s eccentricity from that of its host binary, challenging the conventional view that high binary eccentricity inevitably leads to highly eccentric disks. This insight reshapes our understanding of planet formation viability in a large fraction of binary star systems and provides a new theoretical framework for interpreting upcoming high‑resolution observations of disks in eccentric binaries.