ALMA로 본 대질량 원시군집의 원시별 변동성 탐색

초록

**
본 연구는 ALMA Band 6(1.3 mm)에서 22개의 대질량 원시군집을 다중 에폭으로 관측하여 383개의 응축체를 식별하고, 정밀한 상대 플럭스 보정과 표준편차·차이 지도 분석을 통해 5개의 변동원을 발견하였다. 특히 I13111‑6228의 초소형 H II 영역이 1년 동안 68 % 증가한 것은 고질량 별 형성에서 폭발적 급증(burst‑mode) accretion이 작동함을 뒷받침한다.

**

상세 분석

**
이 논문은 대질량 원시군집 내 원시별의 밀리미터 파장 변동성을 최초로 통계적으로 조사한 점에서 의미가 크다. 먼저 22개의 군집을 선정하고, QUARKS와 MaMMOtH 두 ALMA 대규모 프로그램에서 얻은 Band 6 데이터를 0.3″ 해상도로 일관되게 재구성하였다. 데이터 파이프라인은 (1) 표준 CASA 캘리브레이션, (2) 라인‑프리 채널 선택 후 컨티뉴엄 이미지 생성, (3) 자체 보정(self‑calibration)으로 신호‑대‑잡음비를 크게 향상시켰으며, (4) 모든 에폭을 동일한 합성빔으로 스무딩해 해상도 차이를 제거하였다.

응축체 추출에는 astrodendro 알고리즘을 적용해 3σ 이상 신호를 갖는 383개의 구조를 정의하고, 동일한 마스크를 사용해 각 에폭에서 피크 강도를 측정하였다. 상대 플럭스 보정은 동일한 시점에 관측된 다중 스케일의 비정상적인 변동을 최소화하기 위해, 비변동성 소스(전역 평균 0.5 % 이하 변동)들을 기준으로 스케일 팩터를 도출하였다.

변동성 검출은 두 가지 독립적인 방법으로 수행되었다. 첫 번째는 각 응축체의 피크 강도 시계열에 대한 표준편차를 구하고, 전체 샘플의 평균 표준편차(≈2 %)보다 현저히 큰 값을 보이는 경우를 후보로 선정하였다. 두 번째는 에폭 간 차이 지도를 만든 뒤, 5σ 이상(σ는 이미지 노이즈)인 잔여 신호를 가진 픽셀 군집을 찾는 방식이다. 두 방법이 모두 일치한 5개의 소스가 최종 변동원으로 확정되었다.

가장 눈에 띄는 변동원은 I13111‑6228으로, 초소형 H II 영역이 1년 사이에 피크 강도가 68 % 증가했으며, 통계적 불확실성은 2 % 수준이다. 이는 기존에 보고된 고질량 별 폭발 사례(S255IR‑NIRS3, NGC 6334I‑MM1 등)와 유사한 급격한 질량 공급을 시사한다. 변동원 비율은 383개 중 5개, 즉 최소 1.3 %이며, 이는 저질량 구역에서 보고된 변동 비율(≈5 %)보다 낮지만, 관측 거리와 복잡한 배경을 고려하면 의미 있는 상한선이다.

논문은 또한 현재 한계점을 명확히 제시한다. 첫째, 시간 간격이 불규칙하고 일부 에폭은 몇 시간 차이만 존재해 장기 변동 탐지에 제약이 있다. 둘째, 플럭스 보정에 사용된 기준 소스가 충분히 많은지에 대한 검증이 필요하며, 특히 강한 자유‑자유 방출을 보이는 H II 영역은 자체 변동성을 내포할 수 있다. 셋째, 0.3″ 해상도는 1000 AU 이하 구조를 완전히 분리하지 못해, 복수 원시별이 하나의 응축체로 합쳐질 위험이 있다. 향후 고해상도(≤0.05″)와 더 긴 모니터링 기간을 결합하면 변동성의 발생 빈도와 물리적 메커니즘을 보다 정밀히 규명할 수 있을 것으로 기대된다.

**