베타 피키스 b의 위성 및 이중 행성 탐색을 위한 최초 천체측정 제한

초록

β 피키스 시스템의 직접 영상 행성 β 피키스 b와 c를 대상으로 GRAVITY와 SPHERE의 고정밀 천체측정 데이터를 분석하였다. 위성(또는 이중 행성)의 존재가 행성의 궤도에 미치는 미세한 흔들림을 모델링하고, 3σ 신뢰수준에서 다양한 궤도 주기에 대한 질량 상한을 제시한다. 50 일 주기에서는 180 M⊕(0.6 MJup) 이하, 300 일에서는 65 M⊕(0.2 MJup) 이하, 700 일에서는 50 M⊕(0.15 MJup) 이하, 1 100 일에서는 30 M⊕(0.1 MJup) 이하의 위성만 허용된다. 이는 β 피키스 시스템에서 위성 존재 가능성을 최초로 천체측정적으로 제한한 결과이다.

상세 분석

본 연구는 β 피키스 b와 c의 장기 천체측정 기록을 활용해 위성 존재 가능성을 정량적으로 평가하였다. 먼저 SPHERE와 GRAVITY에서 얻은 34개의 SPHERE 측정치와 12개의 GRAVITY 측정치를 정밀 보정하여 RA·DEC 좌표와 오차를 확보하였다. 데이터 전처리 단계에서 NACO 등 낮은 정밀도 데이터를 제외하고, 각 관측기의 시스템 오차를 재조정(j‑factor)함으로써 실제 측정 불확실성을 반영하였다.

위성 검출을 위한 핵심은 행성 중심(photocenter)과 질량 중심(center of mass) 사이의 미세한 이동을 모델링하는 것이었다. 저자는 β 피키스 b의 광도에 비해 위성이 차지하는 K‑밴드 밝기 비율을 ATMO 진화 모델을 이용해 추정했으며, 2 MJup 위성조차 전체 밝기의 5% 미만에 불과하므로 광도 효과는 무시할 수 있다고 결론지었다. 따라서 위성의 존재는 순수히 중력에 의한 궤도 흔들림으로만 나타난다.

두 행성(β 피키스 b와 c)의 궤도는 Orbitize! 패키지를 기반으로 한 7차원(반장축, 이심률, 경사, 근점인자, 평균이동, 질량 등) 파라미터 공간을 MCMC 방식으로 탐색하였다. 로그우도 함수는 관측값과 모델값의 차이를 가우시안 오차로 가정해 구성했으며, 각 파라미터에 대한 사전분포는 기존 문헌(La​grange et al. 2020)에서 제시된 값들을 사용하였다.

위성 효과는 β 피키스 b의 궤도에 추가적인 3차원 진동을 도입하는 형태로 구현하였다. 위성 질량(Mmoon)과 궤도 주기(Pmoon)를 자유 파라미터로 두고, 각 조합에 대해 잔차(Residual) 제곱합을 계산해 3σ 신뢰구간을 도출했다. 결과는 주기별 질량 상한이 비선형적으로 감소함을 보여준다. 짧은 주기(≈50 일)에서는 위성의 궤도 반지름이 작아 행성의 중심 이동이 크게 나타나므로 질량 제한이 상대적으로 높으며, 장기 주기(≈700–1 100 일)에서는 위성의 궤도 반지름이 커져 행성 중심에 미치는 가속도가 감소해 더 낮은 질량까지 탐지 가능해진다.

또한, 위성 질량‑주기 관계를 기존 이론(달 형성 모델, 포획 시나리오)과 비교했을 때, 현재 제한은 목성형 행성 주변에 존재할 수 있는 대형 위성(수십 M⊕ 수준)까지 배제한다는 점에서 의미가 크다. 특히 β 피키스 b가 약 10 MJup의 질량을 가지고 있음에도 불구하고, 0.1–0.6 MJup 수준의 위성은 관측적으로 허용되지 않는다. 이는 행성‑위성 시스템 형성 이론에 중요한 제약을 제공한다.

마지막으로 저자는 향후 GRAVITY+와 ELT‑MICADO 같은 차세대 인터페라메트리 장비가 10 μas 이하의 정밀도를 달성하면, 지구질량 위성까지도 탐지 가능할 것으로 전망한다. 현재 연구는 이러한 미래 관측 전략의 기반을 마련하고, β 피키스와 같은 젊은 별 주변의 위성 탐색에 천체측정 방법이 유효함을 입증한다.