Kick Velocities and Mass Function of Free-Floating Planets from Dynamical Ejection in Hierarchical Three-Body Systems

Technische Samenvatting: Kick-snelheden en de Massafunctie van Vrijzwevende Planeten door Dynamische Ejectie in Hiërarchische Drie-Lichamenstelsels

Probleemstelling
Vrijzwevende planeten (FFP's), of rogue planeten, zijn substellaire objecten die niet gebonden zijn aan een gastster. Hun snelheid en massaverdeling dienen als kinematische fossielen van hun vormings- en ejectiegeschiedenis. Hoewel gravitationele microlensing een populatie van FFP's heeft gedetecteerd die varieert van aardmassa tot Jupitermassa, zijn de dynamische mechanismen die de ejectiesnelheden en de resulterende massafunctie beheersen, nog onvolledig gekarakteriseerd. Specifiek is er een behoefte om de uitkomsten van de ejectie van hiërarchische drie-lichamenstelsels (een gastster, een massieve reusachtige verstoring en een lichtere binnenste planeet) in kaart te brengen naar de observeerbare eigenschappen van de Galactische FFP-populatie, inclusief hun asymptotische snelheden ($v_\infty$) en de tijdschalen van ejectie ($t_{ejec}$).

Methodologie
De auteurs presenteren een reeks directe $N$-body simulaties met behulp van de rebound integratiepackage met de IAS15 integrator (15e-orde adaptieve tijdstap) om hoge nauwkeurigheid tijdens nauwe ontmoetingen te garanderen. De studie richt zich op het beperkte drie-lichamenprobleem, waarbij de binnenste planeet wordt behandeld als een testdeeltje ten opzichte van een zonachtige gastster ($M_\star = 1 M_\odot$) en een massieve reusachtige verstoring.

Het simulatiegrid is gecontroleerd en systematisch:

• Vaste Parameters: Gastmassa ($1 M_\odot$), semi-grote halve as van de reusachtige planeet ($a_J = 5,2$ AU, een Jupiter-analoog) en initiële inclinaties ($0^\circ$–$3^\circ$).
• Variabele Parameters:
  • Geëjecteerde planetenmassa ($m_p$): $10^{-5}$ tot $10^{-1} M_{Jupiter}$.
  • Reusachtige planetenmassa ($M_J$): $1$ tot $75 M_{Jupiter}$.
  • Geëjecteerde planeet semi-grote halve as ($a_p$): $0,7$ tot $1,5 a_J$.
  • Excentriciteiten van beide planeten ($e_p$ en $e_J$): $0,0$ tot $0,9$.
• Realisaties: Elke configuratie wordt gesimuleerd met $N=10.000$ onafhankelijke realisaties met gedeelde random seeds voor de gemiddelde anomalieën om de effecten van parameters te isoleren.
• Metrieken: De studie meet de asymptotische ejectiesnelheid ($v_\infty$), de eindige-straal kick-snelheid ($v_{kick}$) op het moment van bevestiging van de ejectie, en de ejectietijdsschaal ($t_{ejec}$). Ejectie wordt gedefinieerd door positieve specifieke energie, een afstand van de gast groter dan $10 \times$ de maximale semi-grote halve as van het systeem, en een naar buiten gerichte radiale snelheid die tien orbitale perioden aanhoudt.

Analytisch Kader
De numerieke resultaten worden geïnterpreteerd via een semi-analytisch kader dat de volgende elementen combineert:

1. Energiebudget: Het energie-increment dat nodig is voor ejectie wordt afgeleid van het verschil tussen de gebonden orbitale energie en de asymptotische kinetische energie.
2. Gravitationele Slingshot: De maximale kick-snelheid wordt begrensd door het slingshot-mechanisme in het ruststelsel van de reus, afhankelijk van de ontmoetingssnelheid $u$ en de orbitale snelheid van de reus.
3. Hill-schaal Verstrooiing: De mediane ejectiesnelheid wordt beheerst door typische ontmoetingen binnen de Hill-sfeer van de reus, schalend als $v_H \propto M_J^{1/3}$.
4. Tisserand-parameter: Gebruikt als een geconserveerde diagnostische methode om orbitale elementen te relateren aan ontmoetingssnelheden.
5. Seculaire Evolutie: Tijdschalen voor instabiliteit zijn gekoppeld aan seculaire pompwerking en de overlap van gemiddelde bewegingsresonantie (MMR).

Belangrijkste Resultaten

1. Testdeeltje Limiet: De massa van de geëjecteerde planeet ($m_p$) heeft een verwaarloosbare impact op de ejectiestatistieken over vier ordes van grootte ($10^{-5}$ tot $10^{-1} M_{Jupiter}$). De ejectiefractie varieert met minder dan 2,3%, en de gemiddelde snelheden blijven stabiel. Dit bevestigt dat voor $m_p/M_J \lesssim 10^{-2}$, de dynamica wordt gedomineerd door de reusachtige verstoring, en het geëjecteerde lichaam zich gedraagt als een testdeeltje.
2. Schaal van de Reusachtige Planetenmassa: De massa van de reusachtige planeet ($M_J$) bepaalt de schaal en tijd van de ejectie.
   • De gemiddelde ejectiesnelheid schaalt als $\langle v_\infty \rangle \propto M_J^{0,42}$.
   • De ejectietijdsschaal schaalt als $\langle t_{ejec} \rangle \propto M_J^{-1,09}$, consistent met de seculaire voorspelling $t_{sec} \propto M_J^{-1}$.
   • De ejectiefractie is niet-monotoon, met een piek nabij $10 M_{Jupiter}$ en een daling bij zeer lage en zeer hoge massa's vanwege respectievelijk seculaire tijdsschaalbeperkingen en Hill-sfeer vangst.
3. Afhankelijkheid van de Semi-grote Halve As: Voor configuraties binnen de resonantie-overlapzone ($|a_p - a_J|/a_J \lesssim 0,40$ voor een $10 M_{Jupiter}$ verstoring) zijn de ejectietijden vlak ($\sim 1,9 \times 10^4$ jaar), beheerst door de snelle chaotische overlap-tijdsschaal in plaats van seculaire evolutie. Gemiddelde snelheden pieken nabij $a_p = 0,9 a_J$.
4. Excentriciteitseffecten (Twee-snelheidsrespons):
   • Binnenste Planeet Excentriciteit ($e_p$): Beïnvloedt primair de hoog-snelheidstaart. Het verhogen van $e_p$ van 0 naar 0,9 verhoogt de maximale kick-snelheid met een factor $\sim 3$ (van $\sim 7,8$ naar $\sim 23$ km/s), maar verhoogt de gemiddelde snelheid slechts met een factor $\sim 1,7$.
   • Reusachtige Planeet Excentriciteit ($e_J$): Heeft een aanzienlijk sterker effect op de staart. Het verhogen van $e_J$ van 0 naar 0,9 verhoogt de maximale kick-snelheid met een factor $\sim 10$ (van $\sim 7,8$ naar $\sim 80$ km/s). Dit wordt gedreven door de niet-lineaire groei van de periastraal-snelheid van de reus ($v_{peri} \propto \sqrt{(1+e_J)/(1-e_J)}$), wat zeldzame, hoog-energetische slingshot-ontmoetingen mogelijk maakt.
5. Analytische Bovenste Enveloppe: De auteurs leiden een analytische formule af voor de maximale ejectiesnelheid op basis van het slingshot-mechanisme en de Tisserand-parameter. Deze enveloppe voorspelt accuraat de bovengrens van de gesimuleerde snelheidverdelingen, met name voor $e_p \gtrsim 0,5$.

Betekenis en Implicaties
Het artikel stelt dat deze resultaten een directe link bieden tussen de architectuur van planetaire stelsels en de observeerbare eigenschappen van de Galactische FFP-populatie:

• Massafunctie: Omdat de ejectiefractie onafhankelijk is van de planeetmassa, volgt de huidige FFP-massafunctie direct de onderliggende planeet-opkomstfrequentie ($dN_{FFP}/dm_p \approx \bar{f}_{eject} \cdot dN_{planet}/dm_p$). Dit impliceert dat als de opkomstfrequentie stijgt naar lagere massa's (zoals gesuggereerd door Kepler-data), de sub-Aarde FFP-populatie numeriek dominant zal zijn.
• Microlensing Signaturen: De getransformeerde microlensing-tijdsschaalverdeling ($dN/d\log t_E$) wordt voorspeld als vlak in $\log t_E$ als de planeet-opkomstfrequentie vlak is, of versterkt bij korte tijdschalen als de opkomstfrequentie stijgt naar lagere massa's. Dit biedt een testbare voorspelling voor de Roman Wide Field Survey en Euclid.
• Kinematica: De bulk van de geëjecteerde planeten heeft asymptotische snelheden van slechts enkele km/s, waardoor ze kinematisch ononderscheidbaar zijn van hun ouderlijke sterpopulatie (het verandert de 3D-snelheidsdispersie met $\sim 1\%$). Echter, systemen met zeer excentrische reusachtige verstoringen ($e_J \gtrsim 0,5$) kunnen een hoog-snelheidstaart produceren ($v_\infty \sim 80$ km/s), wat potentieel hoog-snelheidsuitbijters creëert die detecteerbaar zijn door Gaia of Roman.
• Ejectietijdsschalen: De meeste ejecties uit massieve of excentrische systemen vinden snel plaats ($10^3$–$10^4$ jaar), vaak tijdens de pre-hoofdreeks fase. Bijgevolg is de ruimtelijke verdeling van de algemene FFP-populatie effectief ononderscheidbaar van het veld, waarbij alleen jonge, recent geëjecteerde objecten een geheugen behouden van hun geboortemilieu.

Beperkingen
De auteurs merken op dat de beperking tot exact drie lichamen de primaire beperking is, aangezien echte systemen meerdere planeten kunnen bevatten die leiden tot cascade-instabiliteiten. Daarnaast bevatten de simulaties geen protoplanetaire schijven, sterrenmassaverlies of een volledig scala aan inclinaties (specifiek het Kozai-Lidov mechanisme), wat de ejectiestatistieken zou kunnen veranderen. Toekomstig werk wordt geïdentificeerd als noodzakelijk om deze factoren te integreren in een volledig populatiesynthesemodel.