Comparative Study of Two Luminous Red Novae I. Progenitor Modeling and Dust Formation

기술 요약: 두 개의 광도 높은 적색 신성(LRN)에 대한 비교 연구 (I. 전구체 모델링 및 먼지 형성)

문제 제기
광도 높은 적색 신성(Luminous Red Novae, LRN)은 쌍성계의 불안정한 질량 이동과 연관된 천체 현상으로, 흔히 공통 외층(common-envelope, CE) 진화 및 항성 병합으로 이어진다. 이러한 사건들은 컴팩트 쌍성계 및 중력파원 형성에도 연결되어 있으나, 전구체의 물리적 특성과 유출물의 질량 예산은 여전히 불분명한 상태로 남아 있다. 기존 연구들은 대개 전구체, 폭발, 그리고 잔해 단계를 개별적으로 다루거나 단일 항성 모델에 의존해 왔다. 전구체의 쌍성 특성과 관측된 과도 현상 및 잔해를 명시적인 질량 이동 불안정성 기준을 포함하여 연결하는 통합적인 분석은 부족한 실정이었다. 본 연구는 AT2021biy와 AT2021blu라는 구체적인 사례를 대상으로, 이들의 쌍성 전구체 질량과 궤도 매개변수를 제약하고, 방출된 가스 질량을 추정하며, 잔해 단계에서 형성된 먼지 질량을 추론하는 것을 목표로 한다.

방법론
저자들은 항성 진화 모델링, 에너지 예산 형식론, 광도 곡선 분석, 그리고 적외선 먼지 모델링을 결고합한 다각적인 접근 방식을 채택하였다:

1. 전구체 모델링 (MESA): MESA 코드(v11554)를 사용하여 두 개의 쌍성 진화 궤적 그리드를 생성하였다. 모델은 AT2021biy의 허블 우주 망원경 관측 데이터와 AT2021blu의 지상/ZTF/Pan-STARRS 데이터를 통해 제약되었다. 본 연구는 아카이브 데이터와 일치하는 시간 척도로 병합이 가능한 전구체를 선택하기 위해 질량 이동 불안정성 기준(준-정적, 외측 로브 넘침, 역학적)을 명시적으로 적용하였다. 모델은 비보존적 질량 이동을 가정하였으며, 질량의 절반은 등방성 재방출을 통해, 나머지 절반은 원반형 토로이드(circumbinary toroid)를 통해 손실되는 것으로 설정하였다.
2. 에너지 예산 추정: MESA 모델을 바탕으로, 저자들은 CE 단계 동안 방출된 외층 질량을 추정하였다. 이들은 궤도 에너지(내입 시작부터 동반성이 기증자의 헬륨 핵과 병합될 때까지의 과정)를 외층의 결합 에너지 및 유출물의 벌크 운동 에너지와 비교함으로써 하한값과 상한값을 모두 계산하였다.
3. 광도 곡선 비교: 에너지 예산 추정치를 다양한 해석적 및 반-해석적 광도 곡선 모델(예: Kasen & Woosley, Matsumoto & Metzger)과 비교하여, 다양한 물리적 가정(재결합, 충격파 상호작용)에 걸친 일관성을 평가하였다.
4. 먼지 질량 추론: NEOWISE 및 Spitzer 데이터를 복사 전달 코드인 dusty(via pydusty)를 사용하여 모델링함으로써, 폭발 후 약 3년까지의 먼지 질량과 온도를 도출하였다. 이는 유출된 총 가스 질량에 대한 독립적인 하한값을 제공하였다.

주요 결과

• 전구체 제약: 연구 결과, AT2021biy의 기증자 질량은 $M_d = 18 - 23 M_\odot$, AT2021blu의 경우 $M_d = 14 \pm 0.5 M_\odot$로 제한되었다. 생존 가능한 전구체들은 높은 질량비($q \gtrsim 3$ for AT2021biy; $q \gtrsim 10$ for AT2021blu)와 각각 $300-700 R_\odot$ 및 $200-300 R_\odot$의 초기 궤도 분리 거리를 갖는 것으로 나타났다.
• 유출물 질량 추정:
  • 에너지 예산: 방출된 외층 질량의 하한값 추정치는 AT2021biy의 경우 $0.03 - 2.98 M_\odot$, AT2021blu의 경우 $0.02 - 0.14 M_\odot$ 범위였다. 상한값은 각각 $1.72 - 7.82 M_\odot$ 및 $1.8 - 3.8 M_\odot$에 달했다.
  • 일관성: 광도 곡선 모델과의 비교 결과, 중간 질량비 영역($3 \lesssim q \lesssim 10$ for AT2021biy 및 $5 < q \lesssim 15$ for AT2021blu)에서 가장 높은 일관성을 보였다.
• 먼지 형성: 측정된 먼지 질량은 추정된 방출 외층 질량보다 1~5 자릿수 낮았으며, 이는 방출된 가스의 아주 작은 부분만이 먼지로 응축됨을 의미한다. 먼지 질량의 진화는 기존의 성간 물질(CSM)과의 충격파 상호작용과 일치하며, 이는 폭발 전 상당한 질량 손실이 있었다는 MESA의 예측을 뒷받침한다.
• 민감도: 추론된 유출물 질량은 가정된 유출 속도에 매우 민감하다. 저자들은 속도가 층별로 다른 프로파일(외층이 내층보다 더 빠르게 움직이는 구조)을 채택할 경우 벌크 운동 에너지가 감소하여, 단일 속도 가정에 비해 하한값 질량 추정치가 최대 2배까지 증가할 수 있음을 언급하였다.

의의 및 주장
본 논문은 질량 이동 불안정성 기준을 모델링에 명시적으로 포함함으로써, LRN의 쌍성 전구체, 폭발, 그리고 잔해를 연결하는 최초의 통합 연구를 제시한다고 주장한다. 불안정성 제약(특히 외측 로브 넘침 기준)을 통해 매개변수 공간을 좁힘으로써, 저자들은 기존의 단일 항성 또는 인구 통계학적 합성 연구보다 더 정밀한 전구체 질량 및 질량비 제약을 제공한다.

본 연구는 질량이 큰 기증자에서 관찰되는 "지연된 역학적 불안정성"(질량 이동이 일정 기간 준-안정적 넘침 이후 폭주하는 현상)이 이러한 사건의 주변에서 발견되는 느리게 움직이는 가스의 검출과 일치함을 강조한다. 또한, 거대한 방출 가스 질량과 상대적으로 작은 먼지 질량 사이의 불일치는 LRN이 효율적인 먼지 제조 공장이 아님을 강화하며, 관측된 먼지는 병합 유출물 자체보다는 기존의 CSM에서 형성되었을 가능성이 높다는 견해를 뒷받침한다.

저자들은 MESA의 구형 대칭 가정, 동반성에 대한 점질량 근사, 그리고 완벽한 에너지 전환 효율($\alpha_{CE}=1$) 가정 등의 한계를 인정한다. 이들은 자신들의 추정치가 이상적인 경우에 대한 지표 역할을 한다고 보며, 질량 추정치를 정교화하기 위해서는 3차원 유체역학 시뮬레이션과 자기 일관적인 동반성 진화를 포함하는 향후 연구가 필요할 것이라고 제언한다.