ExoplaNeT accRetion mOnitoring sPectroscopic surveY (ENTROPY) - II. Time series of Balmer line profiles of Delorme 1(AB)b

기술 요약: ENTROPY II. Delorme 1(AB)b의 발머 선 프로파일 시계열 분석

문제 및 배경
행성 질량 영역에서의 강착 과정은 행성 및 주변 행성 원반(CPD)의 형성 및 진화에 결정적인 역할을 함에도 불구하고 여전히 제약이 많은 상태이다. 활발한 강착을 나타내는 방출선들이 수십 개의 행성 질량 천체(PMO)에서 검출되었으나, 이러한 강착 흐름을 구동하는 메커즘—구체적으로 이것이 자기권 강착 퍼널(magnetospheric accretion funnels)에서 기인하는지 아니면 국소적인 강착 충격(localized accretion shocks)에서 기인하는지—은 아직 완전히 이해되지 않았다. 또한, 프로토플래닛 탐색 전략에 정보를 제공할 수 있는 이들 선의 변동성 시간 척도 역시 행성 질량 동반체에 대해서는 제대로 규명되지 않았다. PDS 70 b, c 및 GQ Lup b와 같은 천체에 대한 기존 연구들은 변동성을 보여주었으나, 이 질량 영역에서 고해외 분해된 선 프로파일에 대한 고해상도 다중 에포크 모니터링은 부족한 실정이다.

방법론
본 연구는 VLT/UVES 기기를 사용하여 2022년 10월부터 2023년 1월 사이 31개 에포크 동안 수행된, 약 13 $M_{Jup}$의 동반성 Delorme 1 (AB)b에 대한 고해상도(R ≈ 50,000) 광학 분광 모니터링 캠페인을 제시한다. 관측은 330–680 nm 파장 범위를 커버하며, 시간 단위에서 년 단위에 이르는 변동성을 조사한다.

데이터 감소 및 분석에는 다음과 같은 몇 가지 핵심적인 기술적 단계가 포함되었다:

1. 추출 및 오염 제거: 호스트 쌍성의 근접성(분리각 ~1.77″)으로 인해, 호스트 별의 헤일로를 제거하기 위한 2단계 프로세스가 개발되었다. 여기에는 2D 스펙트럼 트레이스 내의 주 신호에 대한 다중 가우시안 피팅이 포함되었으며, 이후 마스터 스펙트럼과 파장 의존적 오염 계수를 사용하여 잔류 항성 오염을 모델링하고 제거하는 후처리 단계가 이어졌다.
2. 선 분해(Line Decomposition): 가우시안 프로파일을 가정하는 대신, L-BFGS-B 최적화 방법을 기반으로 한 커스텀 알고리즘을 사용하여 발머 방출선(H$\alpha$부터 H10까지)을 가변 진폭을 가진 두 개의 정적 성분(고정된 형태)으로 분해하였다. 이러한 경험적 접근 방식은 특정 물리적 프로파일을 사전에 가정하지 않고 복잡한 선 형태를 처리하기 위해 선택되었다.
3. 모델링: 분해된 선 성분들을 두 가지 이론적 프레임워크와 비교하였다:
   • 자기권 강착 모델: Thanathibodee 등(2019)의 퍼널 흐름 모델 사용.
   • 충격 모델: Aoyama 등(2018)의 강착 충격 모델 사용.
4. UV 과잉 분석: 자외선 연속체 과잉(약 3700 Å 미만)은 국소 열평형(LTE) 상태의 수소 방출에 대한 슬랩(slab) 모델을 사용하여 모델링되었으며, 이를 통해 강착 광도와 강착률을 도출하여 선 플럭스로부터 도출된 값들과 비교하였다.

주요 결과

• 검출 및 변동성: 연구는 H$\alpha$부터 H10까지의 발머 선과 가변적인 UV 연속체 과잉을 검출하여 지속적인 강착을 확인하였다. 모든 선과 UV 과잉은 변동성을 보였다.
• 두 성분 구조: 방출선은 성공적으로 두 개의 뚜렷한 성분으로 분해되었다:
  • 깊은 흡수 특징과 확장된 속도 날개(최대 ±500 km/s)를 특징으로 하는 넓은 "날개(wings)" 성분.
  • 더 대칭적이고 고원(plateau) 형태의 프로파일을 가진 더 좁은 "핵(core)" 성분.
• 물리적 기원:
  • 날개 성분은 UV 과잉과 강한 상관관계를 보이며, 자기권 강착 퍼널 모델에 의해 질적으로 재현된다. 이 프로파일의 특징(비대칭성, 넓은 날개)은 단순한 충격 모델과는 일치하지 않는다.
  • 핵 성분은 UV 과잉과 상관관계가 상대적으로 낮다. 이는 1차원 강착 충격 모델로 재현될 수 있으나, 프로파일의 유사성은 활성 M형 왜성에서도 나타나는 특징으로, 이 또한 흑점 활동(chromospheric activity)에서 기인했을 가능성을 시사한다.
• 변동성 시간 척도:
  • 시간 단위에서는 선들이 낮은 상대적 플럭스 변동성(∼8–10%)을 보인다.
  • 주 단위에서 년 단위의 시간 척도에서는 변동성이 크게 증가하여 최대 ∼100%에 달한다.
  • 넓은 날개 성분은 가장 높은 상대적 변동성을 보이며, 특히 2022년 10월 관측된 강착 폭발(outburst) 동안 선 프로파일 형태의 변화를 주도한다. 좁은 핵 성분은 절대적 플럭스 변동의 대부분을 주도한다.
• 강착률: UV 과잉 모델링 결과, 대부분의 에포크에서 강착률은 $\dot{M} \approx (1–3) \times 10^{-12} M_{\odot}/yr$로 나타났으며, 2022년 10월 폭발 중에는 $\sim 6.6 \times 10^{-12} M_{\odot}/yr$의 정점에 도달했다. 이 비율들은 날개 성분 모델링과는 일치하지만, 핵 성분의 충격 모델링으로부터 도출된 비율보다는 현저히 낮다.

의의 및 결론
본 논문은 넓은 날개 성분의 특성이 Delorme 1 (AB)b에 대한 자기권 강착 시나리오를 강력하게 뒷받침하며, 이는 행성 질량 동반체에서 이 메커니즘에 대한 최초의 분해된 선 프로파일 증거를 제공한다고 결론짓는다. 좁은 핵 성분의 기원은 모호하게 남아 있으며, 강착 충격과 흑점 활동이 혼합되어 있을 가능성이 있다.

이 연구는 해당 천체가 오래되었음에도 불구하고(30–45 Myr), 질량 저장고(아마도 "피터 팬" 원반)를 유지하며 계속해서 강착하고 있음을 강조한다. 시간 단위의 낮은 변동성과 장기적인 큰 변화를 특징으로 하는 변동성 거동은 T Tauri 항성과 유사하다. 저자들은 시간 단위의 중간 정도의 변동성이 프로토플래닛 H$\alpha$ 이미징 탐사의 낮은 검출률을 설명하기에는 불충분할 수 있으며, 이는 장기적인 변동성이나 다른 요인들이 더 관련이 있을 수 있음을 시사한다고 언급했다. 마지막으로, 본 연구는 행성 질량 천체에서 작동하는 복잡한 강착 기하학 및 메커니즘을 구분하기 위해 다중 에포크, 다파장 분석이 필수적임을 강조한다.