HPF로 측정한 저질량 M왜성 이중성의 궤도와 물리량 정밀 분석

초록

Habitable‑Zone Planet Finder(HPF) 고해상도 근적외선 분광기를 이용해 M dwarf 이중성 3개(LSPM J0515+5911, NLTT 43564, NLTT 45468)의 방사속도와 광도 데이터를 수년간 축적하였다. LSPM J0515+5911은 이중선 스펙트럼(SB2)으로 126.948 d 주기의 궤도를, 나머지 두 시스템은 단일선(SB1)으로 각각 1877 d와 9.686 d 주기를 구했다. 최소 질량과 1차 질량을 SED 피팅과 모델링을 결합해 추정했으며, 관측된 질량‑반지름 관계가 기존 저질량 별 모델과 차이를 보임을 확인하였다.

상세 분석

본 연구는 HPF의 R≈53 000, 0.8–1.3 µm 파장대 고해상도 분광 데이터를 이용해 M dwarf 이중성의 동역학적 매개변수를 정밀하게 측정한 점이 가장 큰 강점이다. 데이터는 2019년부터 2025년까지 수십 회에 걸쳐 수집되었으며, 각 관측은 배경 별 오염을 최소화하기 위해 전략적으로 궤도 위상에 맞춰 진행되었다.

RV 추출은 두 단계로 이루어졌다. 먼저 LFC(Laser Frequency Comb) 기반의 파장 보정과 고정밀 LSF 모델링을 적용해 절대 파장을 재구성했으며, 이후 telluric absorption을 제거하기 위해 물, O₂, CO₂, CH₄의 라인‑바이‑라인 전이 모델을 HPF LSF와 컨볼루션한 복합 텔루어 모델을 사용하였다. 이 과정은 특히 NIR에서 강하게 나타나는 수증기 라인에 대한 보정 정확도를 크게 향상시켰다.

SB2인 LSPM J0515+5911에서는 두 성분의 스펙트럼을 교차상관법(CCF)과 다중 가우시안 피팅으로 분리해 각각의 RV 시계열을 얻었다. 최소 질량 M sin³i는 0.058 M⊙(1차)와 0.046 M⊙(2차)로, i가 알려지지 않아 실제 질량은 이보다 크다. 반면 NLTT 43564와 NLTT 45468은 SB1이므로 1차 질량을 사전 광도·색 지표와 Gaia DR3 파라미터를 이용해 추정하고, 질량 함수와 관측된 RV 반주기를 결합해 최소 2차 질량을 제한하였다.

광도와 색 정보를 활용한 SED 피팅은 BT‑Settl 및 PARSEC 모델을 기반으로 수행했으며, 각 별의 유효 온도, 반지름, 금속성을 동시에 추정했다. 특히 NLTT 45468의 경우 짧은 9.686 d 주기에 비해 1차 질량 0.35 M⊙(±0.02)와 반지름 0.36 R⊙(±0.03)으로, 기존 저질량 별 진화 모델이 예측하는 값보다 약 5–10 % 큰 편차를 보였다.

이 연구는 관측된 질량‑반지름 관계가 현재 모델(예: Baraffe 2015, Dartmouth)과 일관되지 않음을 재확인한다. 특히 0.05–0.07 M⊙ 범위의 초저질량 별에서 모델이 과소평가하는 경향이 뚜렷하며, 이는 대류 효율, 자기장 억제, 혹은 금속성 차이와 연관될 수 있다.

한계점으로는 세 시스템 모두 식이(전이) 현상이 없으며, i를 직접 측정할 수 없다는 점이다. 따라서 실제 질량은 최소 질량보다 크게 될 가능성이 있다. 또한, SB1 시스템에서 2차 성분의 스펙트럼이 검출되지 않아 온도·반지름을 직접 구할 수 없으며, 이는 모델 의존도가 높아지는 원인이다. 향후 고해상도 광학/근적외선 인터페라미터(예: CHARA, VLTI)와 장기 광도 변이 관측을 결합하면 궤도 기울기와 실제 질량을 직접 구할 수 있을 것이다.