Mass and Spin Growth of Very Massive Stars in Star Clusters Potentially Associated with Little Red Dots

기술 요약: 리틀 레드 닷(Little Red Dots)과 연관될 가능성이 있는 성단 내 초거대 질량 별의 질량 및 스핀 성장

문제 제기
본 연구는 밀집 성단 내에서 초기 질량이 $\sim 10^2\text{--}10^4 \, M_\odot$에 달하는 초거대 질량 별(Very Massive Stars, VMSs)의 형성과 진화를 조사한다. 이 성단들은 최근 제임스 웹 우주 망원경(JWST)에 의해 식별된 붉고 거대한 압축 천체 집단인 "리틀 레드 닷(LRDs)"의 전구체로 제안된다. LRD는 넓은 방출선을 특징으로 하지만, 종종 검출 가능한 X선이나 라디오 방출이 결여되어 있다. LRD에 대해서는 두 가지 주요 모델이 존재한다: "BH 엔벨로프 모델"(밀집 가스에 둘러싸인 거대 블랙홀)과 "별-전용 모델"(극도로 밀집된 항성계). 두 시나리오 모두에서 충돌에 의한 질량 폭주(runaway stellar collisions)는 VMS를 형성하는 유력한 메커니즘이며, 이는 이후 중간 질량 블랙홀(IMBH)로 붕괴할 수 있다.

이전 수치 시뮬레이션의 결정적인 공백은 VMS의 충돌 후 상태에 대한 처리이다. 별의 충돌은 팽창된 반경과 낮은 유효 온도($\sim 5000$ K)를 가진 하야시 궤도(Hayashi track) 상에 머무는 부풀어 오른(bloated) 천체를 생성할 것으로 예상된다. 이러한 부풀어 오른 상태는 충돌 단면적을 크게 증가시켜, 추가적인 폭주 성장을 용이하게 할 가능성이 있다. 그러나 이전 연구들은 종종 더 높은 온도를 가정하거나 이러한 지속적인 팽창을 충분히 고려하지 않았다. 또한, 이러한 특정 성단 조건 내에서 VMS의 스핀 각운동량—그 붕괴의 성격과 잠재적인 중력파(GW) 신호를 결정하는 데 중요한 요소—은 포괄적으로 추적되지 않았다.

방법론
저자들은 네 가지 서로 다른 밀집 성단 모델의 진화를 모델링하기 위해 PETER 코드를 사용하여 중력 $N$-체 시뮬레이션을 수행하였다. 시뮬레이션에는 다음과 같은 물리적 및 수치적 가정이 포함되었다:

• 초기 조건: 성단은 질량 분리나 프랙탈 구성이 없는 King 모델($W_0=12$)을 따랐다. 이들은 금속 함량이 $Z = 2 \times 10^{-3}$ ($0.1 \, Z_\odot$)인 고립계로 취급되었다.
• 항성 집단: 두 가지 초기 질량 함수(IMF)가 테스트되었다: 표준 Kroupa IMF와 "상부 편중형(top-heavy)" IMF ($10 \le m/M_\odot \le 150$ 범위에서 $dN/dm \propto m^{-1}$).
• 항성 진화: 시뮬레이션은 질량이 최대 $1280 \, M_\odot$에 이르는 별들을 위한 커스텀 단일 항성 진화 모델("L 모델")을 활용하였다. 항성풍은 뜨거운 거대 질량 별 및 로 luminous blue variables(LBVs)를 위한 BSE 패키지 처방을 사용하여 모델링되었다.
• 충돌 물리: 충도는 질량 손실이 없는 완벽하게 끈적한(perfectly sticky) 충돌로 처리되었다. 병합 시, 충돌 생성물은 전구 성의 질량비에 따라 의존적인 "재활성화(rejuvenation)" 과정을 거쳤다.
• 부풀어 오른 상태 구현: 본 연구의 핵심 방법론적 혁신은 충돌 생성물($M > 150 \, M_\odot$)에 대한 부풀어 오른 상태의 명시적 모델링이다. 별의 반경은 유효 온도 $T_{\text{eff}} = 5000$ K를 유지하도록 고정되었으며, 광도는 부풀지 않은 상태의 광도를 보존하였다. 이는 각 성단 모델에 대해 부풀어 오른 상태를 포함하는 시뮬레이션("Y" 접미사)과 포함하지 않는 시뮬레이션("N" 접미사) 두 세트로 적용되었다.
• 사후 처리 검증: 저자들은 VMS의 강착률이 임계값($\dot{M}_{\text{crit}} = 3 \times 10^{-2} \, M_\odot \, \text{yr}^{-1}$)을 초과하는지 확인하고, 충돌 사이의 시간이 켈빈-헬름홀츠(KH) 시간 timescale(부풀어 오른 반경과 부풀지 않은 반경 모두에 대해 계산됨)보다 짧은지 체크함으로써 부풀어 오른 상태 가정의 자기 일관성을 평가하였다.

주요 결과
시뮬레이션은 모든 모델에서 질량 폭주가 발생하여 $1000 \, M_\odot$을 초과하는 질량을 가진 VMS 형성을 입증하였다.

• 질량 성장: 부풀어 오른 상태의 포함은 VMS 성장을 현저히 강화했다. 가장 극단적인 경우(모델 M6D7H-Y), VMS 질량은 1 Myr 이내에 $\sim 9.5 \times 10^4 \, M_\odot$(전체 성단 질량의 약 $10\%$)에 도달했다. 이는 부풀지 않은 시나리오(모델 M6D7H-N, $\sim 2.4 \times 10^4 \, M_\odot$)에서 달성된 질량보다 약 3배 높다. 강착률이 0.5 Myr 이후에 임계값 아래로 떨어지는 모델에서도, 초기 부풀어 오른 단계 덕분에 VMS는 비부풀림(non-bloated) 사례보다 훨씬 높은 질량에 도달할 수 있었다.
• 스핀 진화: 생성된 VMS들의 무차원 스핀 파라미터($J c / G M^2$)는 모든 경우에서 10을 초과하며 매우 높게 나타났다. 이러한 빠른 회전은 충돌하는 별들의 궤도 각운동량이 VMS의 스핀 각운동량으로 전환되는 것에 기인한다. 저자들은 시뮬레이션 결과와 일치하는 분석적 추정식(Eq. 15)을 도출하였으며, 스핀이 충격 매개변수(impact parameter)의 제곱근과 충돌 횟수에 비례함을 보여주었다.
• 성단 의존성: VMS 질량비는 초기 반질량 이완 시간(half-mass relaxation time)과 역상관 관계를 보였다. 더 높은 초기 질량 밀도, 상부 편중형 IMF, 그리고 (특정 구성에서의) 더 적은 총 별의 개수가 더 효율적인 질량 폭주를 촉진했다.

의의 및 주장
본 논문은 LRD 전구체와 유사한 밀집 성단에서 형성되는 VMS가 자연스럽게 높은 스핀을 가지며, $10^3\text{--}10^4 \, M_\odot$(또는 특정 부풀림 시나리오에서는 그 이상)의 질량 범위에 도달할 수 있다고 주장한다. 저자들은 이러한 빠르게 회전하는 VMS가 붕괴할 때, 거대한 강착 원반에 둘러싸인 블랙홀을 형성할 것으로 예상한다고 상정한다. 이 BH-원반 시스템은 다음과 같은 현상의 동력이 될 수 있는 plausable 엔진이다:

1. 전자기적 과도 현상: 점성 각운동량 수송을 통한 폭발을 일으키고 잠재적으로 상대론적 제트(초장기 감마선 폭발)를 발사함.
2. 중력파: 거대 강착 원반의 비축퇴적 불안정성(non-axisymmetric instabilities)으로 인해 버스트 중력파를 방출하며, 이는 쌍 불안정성 질량 간극(pair-instability mass gap)에 해당하는 잔류 질량을 포함하는 GW190521 또는 GW231123과 같은 사건을 설명할 수 있음.

저자들은 자신들의 결과가 상한치(upper limits)를 나타낸다는 점을 명시적으로 언급하였다. 시뮬레이션은 충돌 중의 질량 및 각운동량 손실, 주계열 이후의 진화, 그리고 특정 VMS 항성풍의 질량 손실 메커니즘을 생략하였다. 또한, 각운동량 수송(예: 자기장을 통한)이 비효식적일 경우 강체 회전(rigid rotation) 가정은 과대평가일 수 있다. 이러한 주의사항에도 불구하고, 본 연구는 부풀어 오른 상태를 고려하는 것이 LRD와 유사한 환경에서 VMS 성장을 정확하게 모델링하는 데 필수적이며, 생성된 천체들이 다중 신호 천문학적 과도 현상(multi-messenger transients)의 강력한 후보라는 결론을 내린다.