Towards an agnostic algorithm for sampling empirical structure models: The case of Uranus and Neptune

Probleemstelling

Het bepalen van de interne structuur van planeten, zoals Uranus en Neptunus, is een fundamenteel maar sterk gedegenereerd probleem. De beschikbare waarnemingsdata (totale massa, straal en gravitationele momenten $J_{2n}$) zijn beperkt (minder dan tien parameters), terwijl een dichtheidsprofiel uit duizenden datapunten bestaat. Dit betekent dat er oneindig veel dichtheidsprofielen bestaan die allemaal consistent zijn met de waarnemingen.

Traditionele methoden, zoals Markov Chain Monte Carlo (MCMC), worden vaak gebruikt om deze ruimte te verkennen. Echter, om de rekentijd haalbaar te maken, vereenvoudigen eerdere studies het probleem vaak door dichtheidsprofielen te parametriseren met polynomen, sigmoid-functies of polytropen. Deze aannames introduceren onbedoelde vooroordelen (bias) en sluiten mogelijk geldige oplossingen uit, vooral in de diepe binnenkant van de planeet. Er is behoefte aan een "agnostische" aanpak die de volledige ruimte van mogelijke oplossingen verkent zonder onnodige aannames over de vorm van het profiel.

Methodologie

De auteurs presenteren een nieuw, op optimalisatie gebaseerd algoritme dat de MCMC-methode vervangt door een gradiëntafdalingsbenadering (gradient descent) voor computerefficiëntie.

1. Agnostische Parametrisatie:

   • In plaats van een vooraf gedefinieerde functie te gebruiken, wordt het dichtheidsprofiel beschreven door $N=1024$ discrete dichtheidswaarden op gelijke afstanden van niveau-oppervlakken.
   • De enige aannames zijn: geen dichtheidsplateaus of -inversies (dichtheid moet monotoon toenemen naar het centrum) en de maximale dichtheid mag de dichtheid van een fysisch haalbaar materiaal niet overschrijden ($\rho_{max} = 20.000$ kg/m³).
   • De dichtheidswaarden worden afgeleid van $N-1$ vrije parameters $p_i$ via een exponentiële som, wat zorgt voor een monotoon stijgend profiel.

2. Optimalisatie-algoritme:

   • Het algoritme gebruikt de Theory of Figures (ToF) om de gravitationele momenten ($J_2, J_4$) te berekenen op basis van een gegeven dichtheidsprofiel.
   • Een kostenfunctie (cost function) wordt geminimaliseerd die de afwijking van de waargenomen gravitationele momenten en de totale massa meet.
   • De optimalisatie verloopt via de Adam-gradiëntafdalingsmethode (Kingma & Ba, 2017).
   • Na elke stap wordt een projectie toegepast om de minimale toename van de dichtheid per stap ($\Delta_{min}$) af te dwingen.

3. Onbevooroordeelde Startgeneratie:

   • Om te voorkomen dat het algoritme vastloopt in lokale minima of bepaalde profielen prefereert, wordt een speciale binair-subdivisie-algoritme gebruikt om startprofielen te genereren.
   • Dit algoritme zorgt ervoor dat de kans om een profiel te genereren evenredig is met het volume van de oplossingruimte, waardoor de startverdeling onbevooroordeeld is.

4. Simulatie-omgeving:

   • De simulaties zijn uitgevoerd voor Uranus en Neptunus met gebruikmaking van de meest recente waarnemingsdata (o.a. French et al. 2024 voor Uranus).
   • Er zijn honderdduizenden runs uitgevoerd om de statistische verdeling van de oplossingen te bepalen.

Belangrijkste Bijdragen

• Nieuwe Agnostische Algorithm: Een efficiënt algoritme dat de volledige ruimte van empirische dichtheidsprofielen verkent zonder de beperkingen van parametrische modellen (zoals polynomen of polytropen).
• Statistische Relatie voor Discontinuïteiten: De auteurs bieden een statistische leidraad voor modellers die werken met een vast aantal lagen. Ze kwantificeren het verband tussen de scherpte van een discontinuïteit (dichtheidsgradiënt) en het verwachte aantal discontinuïteiten in een oplossing.
• Locatie van Discontinuïteiten: Identificatie van waar in de planeet discontinuïteiten het meest waarschijnlijk voorkomen.
• Vergelijking met Bestaande Modellen: Een uitgebreide vergelijking met eerdere studies (o.a. Nettelmann et al., Helled & Fortney, Movshovitz & Fortney) om te tonen hoe hun aannames de oplossingruimte hebben beperkt.

Resultaten

1. Dichtheidsprofielen:

   • De oplossingsruimte is het minst beperkt in het uiterste binnenste van de planeten.
   • De mediaan van de dichtheidsprofielen toont een geleidelijke toename, consistent met compositiesgradiënten, zonder dat het profiel vlak wordt in het centrum (in tegenstelling tot sommige eerdere modellen).
   • De centrale dichtheid voor Uranus ligt rond de $\sim 6$ g/cm³ en voor Neptunus rond de $\sim 7.5$ g/cm³.
   • De 95% betrouwbaarheidsintervallen voor de centrale dichtheid zijn groot ($\Delta\rho \approx 11$ g/cm³ voor Uranus en $13$ g/cm³ voor Neptunus), wat aantoont dat de data de interne structuur niet uniek bepaalt.

2. Vergelijking met Literatuur:

   • Modellen met scherpe discontinuïteiten uit de literatuur neigen naar een hoge centrale dichtheid (boven de 84e percentiel).
   • Modellen zonder discontinuïteiten neigen naar het piek van de verdeling, maar missen de variatie in de diepe binnenkant.
   • De nieuwe methode toont een bredere spreiding, vooral voor hogere centrale dichtheden, die door eerdere methoden vaak werd gemist.

3. Discontinuïteiten:

   • Als discontinuïteiten worden gedefinieerd als een dichtheidsgradiënt $> 0.02$ kg m⁻⁴, dan hebben de meeste oplossingen maximaal één dergelijke discontinuïteit.
   • Minder scherpe discontinuïteiten (gradiënt 0.002 - 0.02 kg m⁻⁴) komen veel vaker voor.
   • Locatie: Als er discontinuïteiten zijn, zijn deze geconcentreerd rond een genormaliseerde straal van $r/R \approx 0.65$ voor Uranus en $r/R \approx 0.7$ voor Neptunus. Dit zou kunnen corresponderen met de overgang tussen de waterstof/helium-rijke buitenlaag en de diepere lagen van water/rots.

4. Afgeleide Grootheden:

   • De berekende traagheidsmomenten en afplattingen komen overeen met eerdere schattingen, maar tonen de onzekerheid in deze waarden aan.

Betekenis en Conclusie

De studie benadrukt dat het interpreteren van planetaire interne structuren niet alleen afhangt van precisie-data, maar ook van de methodologie die wordt gebruikt om die data te interpreteren. Door de volledige oplossingsruimte te verkennen zonder vooroordelen, onthult dit onderzoek dat er veel meer variatie in mogelijke interne structuren is dan eerder werd aangenomen.

De conclusie is dat de interne structuur van Uranus en Neptunus niet uniek bepaald kan worden met de huidige data; er bestaat een breed spectrum van fysisch mogelijke modellen. Het algoritme biedt een robuust kader voor toekomstig onderzoek, inclusief de analyse van chemische samenstelling en thermische structuren, en kan worden toegepast op andere planeten met gravitationele veldgegevens. De auteurs waarschuwen dat het aannemen van een vast aantal lagen of specifieke functievormen in modellen een significant deel van de mogelijke realiteit kan uitsluiten.