거대한 충돌이 만든 거대 트랜스‑넵튠 천체 위성들의 초기 형성 메커니즘

초록

본 연구는 1000 km 이상 직경을 가진 대형 트랜스‑넵튠 천체(TNO)들의 위성들이 모두 거대한 충돌에 의해 형성되었음을, 수치 유체역학(SPH) 시뮬레이션과 장기 조석 진화를 결합해 증명한다. 충돌 파라미터에 따라 위성의 질량비(γsp)가 10⁻¹부터 10⁻⁴까지 다양하게 나타나며, 초기 유동체 상태에서 일정 기간(10⁴–10⁶ yr) 후 강체로 전이될 때 현재 관측되는 낮은 궤도 이심률과 동기 회전이 재현된다. 이는 위성 형성 시점에 원시 천체가 부분 혹은 전부가 용융된 상태였음을 시사한다.

상세 분석

이 논문은 대형 TNO(지름 > 1000 km)의 위성계가 모두 거대한 충돌에 의해 생성됐다는 가설을 검증하기 위해 434개의 SPH 시뮬레이션을 수행하였다. 충돌 대상과 충격체를 각각 차등(differentiated) 혹은 비차등(undifferentiated) 구조로 설정하고, 질량비(Mtar : Mimp = 2 : 1), 충돌 속도(vimp ≈ 1.0–1.4 vesc) 및 충돌각(θimp = 45°–75°)을 다양하게 조정했다. 결과는 크게 두 가지 위성 형성 메커니즘을 보여준다. 첫째, 파편 디스크에서 재축적되는 ‘디스크‑기원 위성’이며, 둘째, 충격체의 큰 파편이 그대로 살아남아 ‘완전한 위성(intact moon)’을 이루는 경우다. 특히 θimp ≥ 45°에서 충돌각이 크면 완전한 위성이 50 % 이상의 확률로 형성되며, γsp는 10⁻³–10⁻¹ 범위에 고르게 분포한다. 차등 구조와 비차등 구조 모두에서 이러한 경향이 나타났지만, 차등 구조에서는 γsp > 10⁻¹인 대형 위성이 더 자주 생성된다.

위성의 초기 궤도는 periapsis ≈ 3–4 Rp, 이심률 eini = 0–1을 보이며, 이는 관측된 저이심률 위성계와 차이가 있다. 이를 해소하기 위해 저자들은 장기 조석 진화 모델을 적용했으며, 물리적 강성(k)과 점성(Q) 값을 시간에 따라 변하게 하여 ‘유동‑상태(fluid‑like)’와 ‘강체‑상태(rigid)’ 전이를 모사했다. 시뮬레이션 결과, 위성이 즉시 강체 상태라면 대부분 높은 이심률(e > 0.1)로 남아 관측과 불일치한다. 반면, 충돌 직후 10⁴ ~ 10⁶ 년 동안 유동‑상태를 유지한 뒤 강체로 전이하면 거의 모든 위성이 원형 궤도로 감쇠하고, 특히 γsp ≤ 10⁻¹인 경우에는 거의 완전한 원형 궤도와 동기 회전을 달성한다. 이는 초기 충돌 후 원시 천체가 충분히 용융되어 낮은 전단 강성을 가졌다는 강력한 증거이다.

또한, 시뮬레이션에서 얻은 질량비와 얼음 함량(fice) 사이의 상관관계는 관측된 밀도 차와도 일치한다. γsp ≳ 10⁻²인 위성은 fice ≈ 0.5(얼음 + 암석)인 반면, γsp ≲ 10⁻²인 소형 위성은 거의 순수 얼음(fice ≈ 1)이다. 이는 플루토‑카론(γsp ≈ 0.12, 두 천체 모두 비슷한 밀도)과 하우메아‑히이아카(γsp ≈ 4.5×10⁻³, 주성은 암석‑얼음 혼합, 위성은 순수 얼음) 시스템을 자연스럽게 설명한다.

마지막으로, 충돌 빈도와 속도는 태양계 초기의 낮은 궤도 이심률(e, i < 10⁻³)과 연관된 ‘콜드 클래식’ TNO들의 형성 환경을 가정하면, 충돌 속도는 vimp ≈ vesc 정도가 되어 위에서 제시한 파라미터 범위와 일치한다. 따라서 거대한 충돌이 대형 TNO 위성의 기원이라는 결론은 동역학적, 열역학적, 관측적 증거가 모두 일관한다는 점에서 설득력을 가진다.