Planar Three-Body Problem: theoretical predictions and simulation results

La gravità è la grande coreografa del cosmo, che detta come stelle, pianeti e buchi neri si muovano nel vasto vuoto dello spazio. Quando due oggetti massicci interagiscono, il loro percorso è prevedibile e stabile, come un pianeta che orbita attorno a una stella in un ciclo perfetto e ripetitivo. Ma se si introduce un terzo oggetto, le regole cambiano completamente. Il problema dei tre corpi è una sfida fondamentale della fisica in cui tre corpi di massa simile si attraggono a vicenda, creando una danza caotica in cui minuscole variazioni iniziali portano a esiti drasticamente diversi in seguito. Poiché questo caos rende impossibile prevedere il percorso futuro esatto di qualsiasi gruppo di tre oggetti, gli scienziati si sono rivolti alla statistica. Invece di chiedere dove finirà un particolare trio, chiedono cosa accada in media attraverso milioni di scenari. Questo approccio statistico aiuta a spiegare come i pianeti si formino in dischi rotanti di gas o come i buchi neri collidano nei centri delle galassie, eventi che spesso coinvolgono disposizioni quasi piatte e bidimensionali piuttosto che il caos tridimensionale completo dello spazio profondo.

In uno studio recente, i ricercatori si sono posti l'obiettivo di testare una specifica teoria statistica contro la realtà di queste interazioni piatte e bidimensionali. Si sono concentrati su uno scenario in cui tre corpi, partendo da fermi, cadono l'uno verso l'altro sotto l'effetto della gravità finché il sistema non si frammenta. In quasi tutti i casi, questa interazione caotica risulta nella formazione di una coppia stretta e legata, mentre il terzo oggetto viene scagliato via nello spazio. Il team voleva sapere se una teoria moderna, che prevede gli esiti basandosi sul flusso di stati possibili attraverso il tempo, potesse descrivere accuratamente ciò che accade quando il moto è confinato in un singolo piano, come avviene spesso nei reali dischi di materia. Per scoprirlo, hanno eseguito un milione di simulazioni al computer per otto diverse combinazioni di masse, tracciando ogni dettaglio di come i sistemi si siano evoluti e infine si siano frantumati.

I ricercatori hanno iniziato costruendo un laboratorio digitale dove potessero controllare le condizioni iniziali con precisionza. Hanno posizionato una coppia di corpi in un'orbita circolare e un terzo corpo, più leggero, a una certa distanza, il tutto disposto su una superficie piana. Hanno poi lasciato che la gravità prendesse il sopravvento, osservando mentre il corpo distante cadeva verso la coppia, innescando uno scompiglio caotico. Il computer ha tracciato il sistema fino a quando non si è stabilizzato in uno stato finale: una coppia binaria e un viaggiatore in fuga. Per garantire di studiare la parte veramente caotica del processo, il team ha filtrato i rari espulsioni rapide che avvenivano quasi immediatamente, concentrandosi solo sui sistemi che trascorrevano abbastanza tempo in uno stato caotico e mescolato prima di separarsi. Ciò ha permesso loro di confrontare direttamente i risultati simulati con le previsioni teoriche, che assumono che il sistema abbia dimenticato il proprio punto di partenza e abbia raggiunto uno stato di equilibrio statistico.

Uno dei risultati più sorprendenti riguardava la probabilità di quale corpo sarebbe stato quello in fuga. La teoria prevedeva che, in un mondo piatto e bidimensionale, le probabilità di fuga seguirebbero un modello specifico basato sulle masse coinvolte. Tuttavia, le simulazioni hanno rivelato qualcosa di inaspettato. I risultati delle simulazioni piatte non corrispondevano alle nuove previsioni teoriche per le due dimensioni. Al contrario, corrispondevano alle previsioni per un mondo tridimensionale con una precisione quasi perfetta, differendo di circa l'uno percento. Ciò suggerisce che, anche quando il moto è fisicamente vincolato a un piano piatto, il comportamento statistico del sistema imita il caso più complesso e tridimensionale. I ricercatori hanno notato che questa era una deviazione significativa rispetto a quanto suggerito dalla teoria attuale per i sistemi planari, suggerendo che la teoria potrebbe aver bisogno di aggiustamenti per tenere conto di come l'energia e la quantità di moto vengono condivise in questi ambienti piatti.

Lo studio ha anche esaminato i dettagli dell'espulsione stessa, in particolare i momenti in cui il terzo corpo riesce a malapena a liberarsi. In queste fughe "marginali", la teoria prevedeva una distribuzione specifica per il momento angolare del corpo in fuga. Le simulazioni hanno confermato questa previsione con alta precisione, mostrando che i dati si allineavano strettamente con l'aspettativa matematica. Allo stesso modo, i ricercatori hanno esaminato la forma delle orbite lasciate dalla coppia rimanente. Hanno scoperto che la distribuzione di quanto queste nuove orbite fossero allungate o circolari corrispondeva al modello teorico in modo sorprendente, anche se il modello si basava su un'ipotesi che i ricercatori non potevano dimostrare pienamente per questo caso specifico. Questo accordo suggerisce che le idee centrali della teoria siano robuste, anche se alcuni dei dettagli più fini riguardanti l'"emissività" del sistema rimangono un mistero.

Un altro ambito chiave di investigazione è stato quanto durassero queste interazioni caotiche. I ricercatori hanno scoperto che la durata del sistema seguiva un modello prevedibile: per le interazioni molto lunghe, il numero di sistemi sopravvissuti diminuiva secondo una specifica legge di potenza, un risultato che corrispondeva alle aspettative teoriche derivate dalle leggi del moto orbitale. Tuttavia, le fasi iniziali dell'interazione mostravano un mix di comportamenti, con alcuni sistemi che si separavano rapidamente e altri che indugiavano in uno stato di escursioni di sub-fuga, in cui i corpi si allontanano ma vengono poi richiamati insieme prima di fuggire definitivamente. Questo mix di comportamenti ha creato una distribuzione complessa di durate che le simulazioni hanno catturato in dettaglio, fornendo un quadro chiaro di come il caos e la stabilità temporanea coesistano in questi sistemi a tre corpi.

Infine, il team ha indagato se la direzione della fuga importasse. In un sistema piatto, la coppia rimanente può orbitare nella stessa direzione del spin totale del sistema (prograda) o nella direzione opposta (retrograda). La teoria prevedeva che questi due esiti dovessero essere ugualmente probabili, indipendentemente dalle masse coinvolte. Le simulazioni, tuttavia, hanno raccontato una storia diversa. I ricercatori hanno osservato che, all'aumentare della differenza di massa tra i corpi, la probabilità di un'espulsione prograda cresceva significativamente rispetto a una retrograda. Questo bias sistematico contraddiceva la previsione teorica di probabilità uguale, suggerendo che la distribuzione di massa gioca un ruolo più critico nel determinare la direzione dell'orbita finale di quanto il modello attuale tenga conto.

Gli autori concludono che, sebbene la loro teoria statistica abbia previsto con successo molti aspetti del problema dei tre corpi, come la distribuzione delle forme orbitali e la tempistica delle interazioni a lunga durata, essa non è riuscita a prevedere accuratamente le probabilità di fuga e la direzione dell'orbita finale nel caso planare. Il fatto che le simulazioni piatte abbiano corrisposto così da vicino alle previsioni tridimensionali è stata una grande sorpresa, indicando che i vincoli di un piano piatto non semplificano il risultato statistico tanto quanto precedentemente ipotizzato. I ricercatori suggeriscono che il fallimento dell'approssimazione "emissività-cieca", che assume che il sistema dimentichi rapidamente il proprio passato, possa essere la causa di queste discrepanze. Pianificano di perfezionare la loro teoria per tenere meglio conto di questi effetti, con l'obiettivo di creare una descrizione statistica più completa che possa essere applicata agli ambienti astrofisici reali, dalla formazione dei pianeti nel passato del nostro sistema solare alle violente collisioni dei buchi neri nei centri di galassie lontane.