A Potential Signature of HD 7977's Passage Among Observed Long-Period Comet Orbits

Technical Summary: A Potential Signature of HD 7977's Passage Among Observed Long-Period Comet Orbits

Problem Statement
The prevailing model for the production of long-period comets (LPCs) posits that the Galactic tide is the dominant perturbation driving bodies from the Oort cloud into the inner solar system. This mechanism is expected to produce a distinct anisotropy in the arguments of perihelion ($\omega$) of dynamically new LPCs (semimajor axes $a > 10^4$ au), specifically a bias toward positive values of $\sin(2\omega_G)$ (denoted as $s2\omega$). However, recent observations of dynamically new LPCs reveal a distribution that is significantly more isotropic than predicted by tide-dominated simulations. Conversely, observed "returning" LPCs ($10^3 < a < 10^4$ au) exhibit anisotropy consistent with tidal models. This discrepancy challenges the assumption that the Galactic tide is the sole driver of recent LPC production. Furthermore, Gaia data suggests a close encounter with the Sun-like star HD 7977 approximately 2.5 Myrs ago, with a median impact parameter of 13,000 au but a non-negligible probability of a much closer approach (4,000 au). This paper investigates whether such a stellar passage could explain the observed isotropy of new LPCs while preserving the anisotropy of returning LPCs.

Methodology
The authors utilized the SCATR N-body code to simulate LPC production under two distinct scenarios:

1. TIDE Simulation: A 2 Gyr integration of $10^6$ test particles representing the Oort cloud, perturbed by the Sun, giant planets, the Galactic disk tide, and "weak" field star passages. This simulation models the epoch where the Galactic tide dominates Oort cloud dynamics.
2. Passage Simulations: A set of 13 shorter (3.5 Myr) simulations modeling the Oort cloud response to a hypothetical passage of HD 7977. These simulations utilized a co-added sample of $\sim 10^7$ particles derived from the TIDE simulation to ensure sufficient statistics for the short post-encounter timescale. The simulations varied the impact parameter of HD 7977 logarithmically from 3,162 au to 50,119 au.

The study compared the simulated distributions of $s2\omega$ and original semimajor axes ($a_{orig}$) against observed LPCs. The observed sample was restricted to comets discovered after 1989 to ensure high-quality orbital data and uniform sky coverage. Dynamically new LPCs were defined as $a > 10^4$ au, and young returning LPCs as $10^3 < a < 10^4$ au. Statistical comparisons employed Kolmogorov-Smirnov (K-S) and Anderson-Darling tests, with particle weights adjusted using Whipple's fading law ($\Phi = n^{-0.6}$) to account for comet survival probabilities.

Key Results

• TIDE Simulation Discrepancy: The tide-dominated simulation successfully reproduced the anisotropic $s2\omega$ distribution of observed returning LPCs but failed to match the observed dynamically new LPCs. The simulation predicted a strong bias toward $s2\omega > 0$ (only 33% of orbits with $s2\omega < 0$), whereas observed new LPCs are nearly isotropic (48.6% with $s2\omega < 0$). The K-S test rejected the null hypothesis that observed new LPCs and the TIDE simulation share the same distribution ($p = 0.0039$).
• HD 7977 Passage Match: Simulations of HD 7977 passages with impact parameters between 6,000 and 10,000 au provided the best fit to observations. In these scenarios:
  • New LPCs ($a > 10^4$ au): The stellar perturbation is strong enough to inject a large population of distant comets ($a \gtrsim 2 \times 10^4$ au) into the loss cone. These bodies have long orbital periods and undergo sufficient tidal evolution to randomize their $s2\omega$ values before reaching the inner solar system 2–3 Myrs post-encounter, resulting in a nearly isotropic distribution.
  • Returning LPCs ($10^3 < a < 10^4$ au): These bodies are primarily derived from the inner Oort cloud ($a \sim 1.7 \times 10^4$ au) injected earlier (within the first 1–2 Myrs). Due to their shorter orbital periods and origin in a region less affected by the full isotropizing effect of the flyby, they retain the anisotropic signature of the Galactic tide, matching observed returning LPC distributions.
• Flux Enhancement: For impact parameters in the 6,000–10,000 au range, the simulations indicate that the current LPC flux is approximately twice the long-term rate expected from Galactic tide dominance alone. Over the past 2.5 Myrs, the flux was elevated by a factor of 7–12.
• Semimajor Axis and Fading: The best-fitting HD 7977 simulations initially showed a discrepancy in the distribution of original semimajor axes compared to observations. However, the authors note that accounting for non-gravitational forces (which typically shift inferred $a_{orig}$ by $\sim 10^{-5}$ au$^{-1}$) or reweighting the Oort cloud density profile to be less centrally concentrated can reconcile these distributions, though a single model simultaneously matching both $a$ and $s2\omega$ remains elusive.

Significance and Claims
The paper argues that the observed isotropy of dynamically new LPCs, combined with the anisotropy of returning LPCs, provides a potential signature of a recent, powerful stellar encounter. The authors conclude that if HD 7977 passed within 6,000–10,000 au of the Sun ~2.5 Myrs ago, the solar system is currently in the latter stages of a comet shower. This scenario implies:

1. Revised Population Estimates: Because the current LPC flux is temporarily elevated by a factor of ~2, previous estimates of the total Oort cloud population (based on current flux rates) should be revised downward by a similar factor.
2. Future Predictions: The analysis predicts that upcoming Gaia data releases will favor an HD 7977 impact parameter within the 6,000–10,000 au range.
3. Historical Context: If confirmed, this encounter would represent the most powerful stellar passage the solar system has experienced in the past 100–300 Myrs.

The authors remain cautious, noting that alternative explanations involving a flatter Oort cloud density profile or specific comet fading laws could also contribute to the observed distributions, but the HD 7977 scenario offers a coherent explanation for the simultaneous presence of isotropic new LPCs and anisotropic returning LPCs.