FORGE'd in the Early Universe: The Effect of Protostellar Outflows on Pop III Accretion

기술 요약: FORGE'd in FIRE 프레임워크를 이용한 Pop. III (원시)별의 시뮬레이션

문제 정의
첫 번째 세대의 별(Population III)의 형성과 진화는 직접적인 관측의 부재와 우주론적 구조 형성으로부터 상세한 원시별 물리학에 이르기까지 요구되는 거대한 동적 범위(dynamic range)로 인해 여전히 제약이 많다. 기존의 이론적 모델들은 종종 이상적인 초기 조건을 사용하거나, 자기장, 복사, 제트와 같은 핵심적인 물리학을 누락하거나, 금속이 풍부한 환경에서 관찰되는 실제적인 초기 질량 함수(IMF) 및 항성 다중성을 생성하는 데 실패했다. 특히, Pop. III 별이 원시별 제트를 분출할 수 있는지에 대한 구체적인 논쟁이 존재하는데, 이는 미량의 자기장 증폭과 회전하는 자화된 주변 성간 원반의 존재 여부에 달려 있다.

방법론
본 연구는 FORGE'd in FIRE 프레임워크를 사용하여 우주론적 "줌인(zoom-in)" 복사 자기 유체 역학(RMHD) 시뮬레이션을 제시한다. 이 시뮬레이션은 완전한 우주론적 맥락 내에서 적색편이 $z \sim 14$에서의 단일한 금속 결핍 Pop. III (원시)별의 형성과 진화를 추적한다.

• 시뮬레이션 프레임워크: 중력, 자기 유체 역학(MHD) 및 복사 전달 방정식을 결합하여 해결하기 위해 질량 보존 유한 체적(MFM) 방법을 사용하는 GIZMO 코드를 사용한다.
• 포함된 물리학: 시뮬레이션은 다음을 포함하는 FIRE-3 물리학 모델을 통합한다:
  • 복사: 가스와 먼지에 자가 일관적으로 결합된 다중 주파수 전달(광이온화, Lyman-Werner 등).
  • 열화학: $1\text{--}10^{10}$ K의 온도와 $10^{-8}\text{--}10^{24}$ cm$^{-3}$의 밀도 범위를 포함하며, H, He, HD 및 먼지 함량을 추적한다.
  • 자기장: 초기 우주 시드 자기장($\sim 10^{-15}$ G)으로 시작하여 자가 일관적으로 증폭된다. 비이데알(non-ideal) MHD 효과(Ohmic 저항, Hall 효과, ambipolar 확산)는 배제되었으며, 이는 포스트 프로세싱 시간 척도 비교를 통해 금속 결핍 가스의 역학적 시간 척도에 비해 무시할 수 있는 수준임을 확인했기 때문이다.
  • 원시별 진화: 싱크 입자(sink particles)는 (원시)별을 나타내며, 강착, 주계열 단계 및 피드백 메커니즘을 통해 진화한다.
• 제트 및 윈드 모델링: 원시별 제트는 Cunningham 등의 모델을 따라, 강착된 질량의 일부($f_W = 0.3$)를 케플러 속도의 일부($f_K = 0.3$)로 정렬된 양극 유출(bipolar outflows)로 전환하여 모델링한다. 금속 결핍 별의 경우 금속 함량에 따라 규모가 달라지므로 항성풍(stellar winds)은 고려하지 않는다.
• 해상도: 시뮬레이션은 "하이퍼-리파인먼트(hyper-refinement)"를 활용하여, (원시)별 근처의 공간 척도 $\sim 10^{-4}$ pc ($< 100$ au) 및 가스 셀 $\sim 10^{-3} M_\odot$까지 해상도를 높이는 동시에, $z \sim 100$부터의 우주론적 맥락을 유지한다.

주요 결과

1. 항성 성장 및 질량:

   • (원시)별은 초기 강착률 $\sim 5 \times 10^{-3} M_\odot \text{yr}^{-1}$에서 시작하여 31,000년 동안 최종 질량 $\sim 27 M_\odot$에 도달한다.
   • 성장률은 별이 "Shell D 연소" 단계에 진입함에 따라 안정화된다. 최종 질량은 주변 저장고의 고갈보다는 피드백(제트 및 복사)에 의해 조절된다.
   • 시뮬레이션 종료 시, 별의 반지름은 $\sim 10 R_\odot$, 유효 온도는 $\sim 3.5 \times 10^4$ K, 광도는 $\sim 1.3 \times 10^5 L_\odot$이다. 켈빈-헬름홀츠 시간 척도($t_{KH}$)가 강착 시간 척도($t_{acc}$)보다 짧아지는데, 이는 별이 추가적인 질량을 강착하기 전에 주계열에 도달할 것임을 나타낸다.

2. 주변 성간 원반:

   • 회전에 의해 지지되는 원반이 $\sim 0.01$ pc ($\sim 2000$ au)까지 형성된다.
   • 원반은 시뮬레이션 전반에 걸쳐 중력적으로 안정적이며, Toomre $Q$ 파라미터가 지속적으로 $>> 1$을 유지한다.
   • 안정성은 열적 지지(내부 영역, $T \sim 10^4$ K)와 난류 지지(외부 영역, $R > 1000$ au)의 결합에 의해 유지된다.
   • 표면 밀도 프로파일은 외부 원반의 $\Sigma \propto R^{-1}$에서 내부 영역의 $\alpha$-원반 프로파일로 전이되며, $\sim 20$ au에서 밀도가 급격히 상승한다.

3. 원시별 제트 및 유출:

   • 제트 분출: 제트는 Machida와 Doi가 정의한 대로, 제트 분출을 허용하는 자기 질량-자속비($\mu$) 영역($\mu \sim 10^{-1} \text{--} 10^1$) 내에 있기 때문에 형성된다. 이는 미량의 시드 자기장을 사용했음에도 불구하고 가능하다.
   • 역학: 제트는 $\sim 1$ pc까지 전파되며, 낙하하는 가스와 충돌하여 종단 충격파(termination shocks)를 형성한다.
   • 조절: 제트는 상당한 운동량을 전달하여 가스를 쓸어내며, 대규모 반경에서 순 강착률을 약 10배 정도 감소시킨다. 이러한 운동량 기반 피드백이 최종 항성 질량을 조절한다.
   • 유출률은 처음 1,000년 동안 $\sim 10^{-3} M_\odot \text{yr}^{-1}$로 정점을 찍은 후 감소하는 반면, 유입률은 약 $3\text{--}5 \times 10^{-4} M_\odot \text{yr}^{-1}$로 안정화된다.

4. 화학 및 열 구조:

   • 내부 원반에서는 높은 온도($T \sim 10^4$ K)로 인해 수소 분자($H_2$)가 억제되고 이온화된 수소가 선호되어, 냉각 기작이 $H_2$에서 Lyman 냉각으로 전환된다.
   • 원반은 주요 냉각 복사에 대해 광학적으로 얇으며, 이를 통해 IR 광자가 효율적으로 탈출할 수 있다.

의의 및 주장
본 논문은 완전한 우주론적 줌인 설정 내에서 개별적으로 해상된 Pop. III (원시)별에 FORGE'd in FIRE 프레임워크를 적용한 최초의 사례를 제시한다고 주장한다. 주요 의의는 다음과 같다:

• Pop. III 형성에 있어 제트의 생존 가능성: 미량의 우주론적 자기장이 있음에도 불구하고, 강력한 제트 분출을 위한 조건이 충족됨을 보여주며, 이는 제트가 금속 결핍 환경에서 존재하지 않을 수 있다는 기존의 가설에 도전한다.
• 피드백에 의한 질량 조절: 원시별 제트는 별이 단순히 주변 저장고를 모두 소비하는 것이 아니라, 운동량 기반 유출을 통해 강착을 조절함으로써 최종 질량을 결정하는 데 결정적인 역할을 한다.
• 원반 안정성: 열적 및 난류적 지지에 의해 유지되는, 파편화되지 않는 안정적인 주변 성간 원반의 형성이 가능하다.

이러한 결과는 초기 별 형성 영역에 대한 더 물리적으로 완전한 그림을 제공하며, Pop. III 초기 질량 함수(IMF)와 피드백 과정에 대한 통찰을 제공한다. 저자들은 향후 연구를 통해 이를 여러 개의 Pop. III 별 그룹으로 확장하여 예측된 IMF와 항성 다중성을 완전히 규명할 것이라고 언급했다.