The quest for Magrathea planets. II. Orbital stability of exoplanets formed around double white dwarfs

기술 요약: 이중 백색왜성 주위에서 형성된 외계 행성의 궤도 안정성

문제 제기
이론적 모델은 공통 외피(common-envelope) 단계에서 방출된 물질로 생성된 공통 회전 원반(circumbinary discs)으로부터 "2세대" 행성이 형성될 수 있음을 시사하지만, 이러한 시스템의 장기적인 역학적 안정성은 아직 검증되지 않았다. 기존 연구들은 이러한 환경에서 거대 행성 형성의 가능성을 입증했으나(Ledda et al. 2023), 다행성 공통 회전 구조에서 발생하는 혼돈스러운 중력 상호작용 속에서 이러한 시스템이 수백만 년의 시간 척도 동안 생존할 수 있는지는 불분명하다. 또한, 이러한 시스템의 안정성은 중력파(GW) 주파수 변조를 통해 이중 백색왜성(DWD) 주위의 외계 행성을 탐지하고자 하는 차세대 레이저 간섭계 우주 안테나(LISA) 임무의 탐지 가능성을 평가하는 데 매우 중요하다.

연구 방법론
저자들은 LISA의 DWD 모집단에서 추출된 세 가지 특정 DWD 구성(DWD2, DWD3, DWD4)을 중심으로 하는 공통 회전 행성계에 대해 $N$-체 시뮬레이션을 수행하였다. 본 연구는 장기적인 케플러 진화와 근접 조우/충돌을 모두 처리하기 위해 이중 성계에 특화되어 조정된 하이브리드 심플렉틱 적분기(Chambers et al. 2002)를 사용하였다.

시뮬레이션은 초기 조건에 따라 네 가지 별개의 모집단을 다루었다:

• 모집단 A: 저질량 쌍성계($M_{bin} < 0.6 M_\odot$) 주변의 저질량 행성($M_p \leq 1.2 M_J$)을 가진 비공명 시스템.
• 모집단 B: 고질량 쌍성계($M_{bin} \sim 1 M_\odot$) 주변의 고질량 행성($M_p \geq 2.5 M_J$)을 가진 비공명 시스템.
• 모집단 C: 모집단 B의 공명 대응체로, 행성들이 초기에 2:1 평균 운동 공명(mean-motion resonances)에 묶여 있는 시스템.
• 모집단 D: DWD*4 주변의 특정 "NMS" 공명(4:3 및 5:3 공명의 조합)을 가진 5개 행성 시스템.

초기 조건에는 Ledda et al. (2023)의 분포를 기반으로 한 행성 질량, 임계 안정 한계($a_{crit}$)와 1 au 사이의 장반경(또는 공명 사례를 위한 스케일링된 값), 그리고 낮은 초기 이심률($e < 0.01$) 및 궤도 경사각($1^\circ-3^\circ$)이 포함되었다. 시뮬레이션은 $3 \times 10^6$에서 $10^7$년 동안 실행되었다.

결과 분석을 위해 저자들은 다음과 같은 지표를 활용하였다:

1. 정규화된 각운동량 결핍(Normalized Angular Momentum Deficit, NAMD): 시스템의 역학적 흥분도(온도)를 정량화하기 위함.
2. 재앙적 사건 추적(Catastrophic Event Tracking): 행성 방출(100 au 너머로 산란), 충돌, 쌍성계와의 근접 접근을 기록.
3. 궤도 간격 및 질량 집중도: 궤도 간격($S_s$)과 가장 무거운 행성의 질량 분율($S_m$)의 변화를 분석.
4. 질량 중심(CoM) 이동: 어떤 행성이 소실되었는지 추론하기 위해 행성계의 질량 중심 변위를 계산.

주요 결과

• 안정성 및 다중성: DWD 주위의 행성계는 연구된 시간 척도 동안 안정적으로 유지될 수 있으나, 안정성은 초기 다중성에 강하게 의존한다. 두 개의 행성으로 시작하는 시스템이 변형 없이 생존할 가능성이 가장 높다. 반면, 세 개, 네 개 또는 다섯 개의 행성을 가진 시스템은 빈번하게 재앙적 불안정성을 경험한다.
• 모집단별 결과:
  • 모집단 A (저질량, 비공명): 82.5%가 초기 행성 수를 유지하였고, 17.5%가 재앙적 사건을 경험하였다.
  • 모집단 B (고질량, 비공명): 49%만이 초기 다중성을 유지하였으며, 51%는 행성 손실을 겪었다.
  • 모집단 C (고질량, 공명): 65.6%가 초기 다중성을 유지하였다. 공명은 저다중성 시스템의 안정성을 향상시켰으나, 네 개의 행성을 가진 공명 시스템은 비공명 대응체보다 파괴되기 더 쉬웠다.
  • 모집단 D (다섯 행성 공명): 100%의 시스템이 재앙적 사건을 경험하였다. 대다수(58%)는 정확히 두 개의 행성을 잃었으며, 17.2%는 완전히 파괴되었다(모든 행성 소실).
• 역학적 진화: 불안정성을 겪은 시스템은 일반적으로 "역학적으로 뜨거운" 상태(NAMD > $1.3 \times 10^{-3}$)로 진화하였으며, 이는 증가된 궤도 간격과 생존한 가장 무거운 행성의 높은 질량 집중도를 특징으로 한다. 최종 다중성과 NAMD 사이에는 강한 역상관 관계가 관찰되었다: 더 적은 수의 행성으로 끝난 시스템은 일반적으로 더 높은 역학적 흥분을 보였다.
• 공명 효과: 초기 2:1 공명이 저다중성(2-3개 행성) 시스템을 보존하는 데 도움이 되었지만, 고다중성 체인(4-5개 행성)의 붕괴를 막지는 못했으며, 이들은 거의 예외 없이 공명을 깨뜨리고 행성을 잃었다.
• LISA 탐지 가능성:
  • 단일 행성 시스템: 전체 시뮬레이션 시스템의 약 7%가 단일 행성 시스템으로 종료되었다. 이 생존자 대부분은 유리한 하늘 위치와 거리가 확보된다면 LISA의 탐지 범위 내에 해당하는 질량과 주기를 가진다.
    간의 전역 중력파 주파수 변동 진폭($\Delta f$)을 추정하기 위한 공식(식 15)을 도출하였다. 대부분의 시뮬레이션 시스템에서 이 변동은 LISA의 주파수 해상도($\sim 7.92 \times 10^{-9}$ Hz, 4년 임무 기준) 미만이다. 오직 일부 공명 시스템(모집단 C)만이 8년 임무 내에 탐지 가능한 잠재적 변동을 보였다.
• 질량 중심(CoM) 가속도: 행성 방출이 쌍성의 질량 중심에 측정 가능한 가속도(특정 사례에서 최대 $\sim 0.9$ au/yr$^2$)를 유발할 수 있음을 주목하였으나, LISA를 위한 구체적인 탐지 임계값은 아직 정의되지 않았다.

의의 및 주장
본 논문은 DWD 주변의 2세대 행성계가 역학적으로 관측 가능한 유망한 후보이지만, 그 최종 구조는 초기 역학적 불안정성에 의해 크게 형성된다고 주장한다. 저자들은 다음을 강조한다:

1. 생존은 흔하지만 수정된다: 완전한 파괴는 드물지만(전체 시스템의 7.7%), 행성 손실은 흔하며, 이는 고다중성 조상으로부터 유래한 2개 행성 시스템의 인구를 122% 증가시킨다.
2. 진단 도구로서의 NAMD: NAMD 지표는 저다중성으로 형성된 시스템과 역학적으로 흥분되어 행성을 잃은 시스템을 구분하는 중요한 도구이다. 관측된 2개 행성 시스템에서 높은 NAMD 값은 불안정성의 역사를 나타낼 수 있다.
3. LISA 잠재력: 본 연구는 "마그라테아(Magrathea)" 행성(DWD 주변의 2세대 행성)을 LISA로 탐지하는 것이 가능하다는 점을 뒷받침하며, 특히 단일 행성 생존자의 경우 그러하다. 다만, 다중 행성 시스템을 주파수 변조를 통해 탐지하는 것은 작은 주파수 변동으로 인해 어려울 수 있으므로, 특정 파형 연구와 질량 중심 가속도 신호에 대한 연구가 필요하다.

저자들은 본 시뮬레이션이 초기 역학적 진화($\sim 10$ Myr까지)를 다루고 있지만, 그 결과가 미래의 LISA 데이터를 해석하기 위한 기초적인 틀을 제공한다고 결론짓는다. 그들은 자신들의 탐지 가능성 주장이 예시적인 것이며, 관측 전망을 확정하기 위해서는 전용 파형 분석이 필요함을 명시적으로 밝히고 있다.