Planar Three-Body Problem: theoretical predictions and simulation results

중력은 우주의 위대한 안무가로서, 별, 행성, 그리고 블랙홀이 광활한 우주의 빈 공간을 통해 어떻게 움직이는지를 결정한다. 두 개의 거대한 천체가 상호작용할 때, 그 경로는 마치 행성이 별 주위를 완벽하고 반복적인 궤도로 공전하는 것처럼 예측 가능하고 안정적이다. 하지만 세 번째 천체가 도입되면 규칙은 완전히 바뀐다. 삼체 문제(three-body problem)는 비슷한 질량을 가진 세 물체가 서로를 끌어당기며, 초기 단계의 미세한 변화가 나중에 완전히 다른 결과를 초래하는 혼돈의 춤을 만들어내는 물리학의 근본적인 과제이다. 이러한 혼돈으로 인해 특정 세 물체 그룹의 정확한 미래 경로를 예측하는 것이 불가능하기 때문에, 과학자들은 통계학에 의존하게 되었다. 특정 세 쌍의 조합이 정확히 어디로 갈 것인지 묻는 대신, 그들은 수백만 개의 시나리오 전반에 걸쳐 평균적으로 어떤 일이 일별되는지를 묻는다. 이러한 통계적 접근 방식은 가스 디스크 내에서 행성이 어떻게 형성되는지, 혹은 은하 중심에서 블랙홀이 어떻게 충돌하는지와 같은 현상을 설명하는 데 도움을 주는데, 이러한 사건들은 대개 심우주의 완전한 3차원적 혼돈보다는 거의 평평한 2차원적 배열을 포함하는 경우가 많다.

최근 한 연구에서 연구진은 이러한 평평한 2차원적 상호작용의 현실에 대해 특정 통계 이론을 검증하고자 했다. 그들은 세 물체가 정지 상태에서 출발하여 중력에 의해 서로를 향해 낙하하다가 시스템이 해체될 때까지의 시나리오에 집중했다. 거의 모든 경우에서, 이러한 혼돈스러운 상호작용은 두 물체가 긴밀하게 결합된 쌍을 형성하고, 세 번째 물체는 우주 멀리 튕겨 나가는 결과로 이어진다. 연구팀은 현대 이론—시간에 따른 가능한 상태들의 흐름을 바탕으로 결과를 예측하는 이론—이 실제 물질의 디스크에서 흔히 발생하는 것처럼 운동이 단일 평면에 국한될 때도 이를 정확하게 설명할 수 있는지 알고 싶었다. 이를 확인하기 위해 그들은 여덟 가지 질량 조합에 대해 100만 번의 컴퓨터 시뮬레이션을 실행하며, 시스템이 어떻게 진화하고 결국 산산조각 나는지에 대한 모든 세부 사항을 추적했다.

연구진은 정밀하게 초기 조건을 제어할 수 있는 디지털 실험실을 구축하며 시작했다. 그들은 두 물체를 원형 궤도에 배치하고, 세 번째로 더 가벼운 물체를 일정 거리 떨어진 곳에 배치한 뒤 모두 동일한 평면 위에 놓았다. 그런 다음 중력이 작용하도록 하여, 멀리 있던 물체가 쌍을 향해 떨어지면서 촉발되는 혼돈스러운 소동을 관찰했다. 컴퓨터는 시스템이 최종 상태, 즉 하나의 이진 쌍(binary pair)과 탈출하는 여행자로 정리될 때까지 추적했다. 진정으로 혼돈스러운 과정만을 연구하기 위해, 연구팀은 거의 즉각적으로 발생하는 드문 빠른 탈출 사례들을 걸러내고, 해체되기 전까지 충분한 시간 동안 혼돈스럽고 뒤섞인 상태를 유지한 시스템들에만 집중했다. 이를 통해 그들은 시스템이 시작점을 잊고 통계적 평형 상태에 도달했다고 가정하는 이론적 예측과 시뮬레이션 결과를 직접 비교할 수 있었다.

가장 놀라운 발견 중 하나는 어떤 물체가 탈출하게 될 것인가에 대한 가능성에 관한 것이었다. 이론은 평평한 2차원 세계에서 탈출 확률이 관련된 질량에 따라 특정한 패턴을 따를 것이라고 예측했다. 그러나 시뮬레이션은 예상치 못한 결과를 보여주었다. 평면 시뮬레이션의 결과는 2차원을 위한 새로운 이론적 예측과 일치하지 않았다. 대신, 그것은 3차원 세계에 대한 예측과 거의 완벽하게 일치했으며, 차이는 약 1%에 불과했다. 이는 운동이 물리적으로 평면에 국한되어 있음에도 불구하고, 시스템의 통계적 행동이 더 복잡한 3차원 사례를 모방한다는 것을 시사한다. 연구진은 이것이 평면 시스템에 대한 현재 이론의 제안으로부터 크게 벗어난 점이라고 언급하며, 이 평면 환경에서 에너지와 운동량이 어떻게 공유되는지를 설명하기 위해 이론이 조정될 필요가 있음을 시사했다.

연구는 또한 탈출의 세부 사항, 특히 세 번째 물체가 간신히 자유를 얻는 순간에 대해서도 조사했다. 이러한 "한계적(marginal)" 탈출에서, 이론은 탈출하는 물체의 각운동량에 대한 특정 분포를 예측했다. 시뮬레이션은 데이터가 수학적 기대치와 밀접하게 일치함을 보여줌으로써 이 예측을 높은 정밀도로 확인했다. 마찬가지로, 연구진은 남겨진 쌍이 남긴 궤도의 형태를 살펴보았다. 그들은 이 새로운 궤도가 얼마나 길게 늘어지거나 원형인지에 대한 분포가 이론적 모델과 놀라울 정도로 잘 일치한다는 것을 발견했는데, 이는 연구진이 이 특정 사례에 대해 완전히 증명할 수는 없었던 가설에 기반한 모델이었음에도 그러했다. 이러한 일치는 비록 '방사성(emissivity)'에 관한 세부적인 부분들이 여전히 미스터리로 남아 있을지라도, 이론의 핵심 아이디어가 견고하다는 것을 시사한다.

또 다른 주요 조사 영역은 이러한 혼돈스러운 상호작용이 얼마나 지속되는가였다. 연구진은 시스템의 수명이 예측 가능한 패턴을 따른다는 것을 발견했다. 매우 긴 상호작용의 경우, 생존하는 시스템의 수는 궤도 운동 법칙에서 유도된 특정 멱법칙(power law)에 따라 감소했는데, 이는 이론적 기대치와 일치하는 결과였다. 그러나 상호작용의 초기 단계는 다양한 행동을 보였는데, 어떤 시스템은 빠르게 분해되었고 어떤 시스템은 몸체들이 서로 떨어졌다가 마침내 탈출하기 전에 다시 끌려 들어오는 '서브 탈출 순회(sub-escape excursions)' 상태에서 머물렀다. 이러한 다양한 행동의 혼합은 복잡한 수명 분포를 만들어냈으며, 시뮬레이션은 이를 상세하게 포착하여 혼돈과 일시적 안정성이 이 삼체 시스템에서 어떻게 공존하는지에 대한 명확한 그림을 제공했다.

마지막으로, 연구팀은 탈출 방향이 중요한지를 조사했다. 평면 시스템에서 남겨진 쌍은 시스템의 전체 스핀과 같은 방향(순행, prograde)으로 궤도를 돌거나 반대 방향(역행, retrograde)으로 돌 수 있다. 이론은 질량에 관계없이 이 두 결과가 동일하게 발생할 것이라고 예측했다. 그러나 시뮬레이션은 다른 이야기를 들려주었다. 연구진은 물체 간의 질량 차이가 커질수록 순행 탈출의 가능성이 역행보다 현저히 높아지는 것을 관찰했다. 이러한 체계적인 편향은 두 결과의 확률이 동일할 것이라는 이론적 예측에 반하는 것으로, 질량 분포가 현재 모델이 고려하는 것보다 최종 궤도의 방향을 결정하는 데 더 결정적인 역할을 한다는 것을 시사한다.

저자들은 통계적 이론이 궤도 형태의 분포나 장기적 상호작용의 타이밍과 같은 삼체 문제의 많은 측면을 성공적으로 예측했지만, 평면형의 경우 탈출 확률과 최종 궤도의 방향을 정확하게 예측하는 데는 실패했다고 결론지었다. 평면 시뮬레이션이 3차원 예측과 매우 밀접하게 일치했다는 사실은 큰 놀라움이었으며, 이는 평면이라는 제약이 통계적 결과를 예상만큼 단순화하지 않는다는 것을 나타낸다. 연구진은 시스템이 과거를 빠르게 잊는다고 가정하는 '방사성 무관(emissivity-blind)' 근사의 붕괴가 이러한 불일치의 원인일 수 있다고 제안한다. 그들은 우리 태양계의 과거 속 행성 형성부터 먼 은하 중심의 격렬한 블랙홀 충돌에 이르기까지, 실제 천문학적 환경에 적용할 수 있는 더 완전한 통계적 설명을 만들기 위해 이론을 개선할 계획이다.