[논문 리뷰] The core helium flash revisited III. From Pop I to Pop III stars

Miroslav Mocák, S. W. Campbell외 2명|arXiv (Cornell University)|2010. 03. 18.

Stellar, planetary, and galactic studies참고 문헌 1인용 수 23

한 줄 요약

이 연구는 저질량 항성에서 핵 헬륨 폭발의 동역학을 다차원 시뮬레이션을 통해 분석하며, 기존 1D 별진화 모델의 한계를 해결하고자 한다. 금속 농도가 높은 Pop I 항성(1.25 M⊙, Z=0.02)에서 난류 유입이 단일 대류권의 급격한 성장을 유도하여 약 23일 후 수소가 헬륨 핵으로 주입되고, 최종적으로 이중 대류권이 형성됨을 밝혀냈다. 반면, 금속이 없는 Pop III 항성(0.85 M⊙, Z=0)에서는 이중 대류권이 내부 중력파에 의해 빠르게 붕괴되며, 3D 모델에서의 대류 속도는 1D 예측보다 2D보다 더 잘 일치한다. 이는 기존 별진화 모델의 가정에 도전하는 결과이다.

ABSTRACT

Degenerate ignition of helium in low-mass stars at the end of the red giant branch phase leads to dynamic convection in their helium cores. One-dimensional (1D) stellar modeling of this intrinsically multi-dimensional dynamic event is likely to be inadequate. Previous hydrodynamic simulations imply that the single convection zone in the helium core of metal-rich Pop I stars grows during the flash on a dynamic timescale. This may lead to hydrogen injection into the core, and a double convection zone structure as known from one-dimensional core helium flash simulations of low-mass Pop III stars. We perform hydrodynamic simulations of the core helium flash in two and three dimensions to better constrain the nature of these events. To this end we study the hydrodynamics of convection within the helium cores of a 1.25 \Msun metal-rich Pop I star (Z=0.02), and a 0.85 \Msun metal-free Pop III star (Z=0) near the peak of the flash. These models possess single and double convection zones, respectively. We use 1D stellar models of the core helium flash computed with state-of-the-art stellar evolution codes as initial models for our multidimensional hydrodynamic study, and simulate the evolution of these models with the Riemann solver based hydrodynamics code Herakles which integrates the Euler equations coupled with source terms corresponding to gravity and nuclear burning. The hydrodynamic simulation of the Pop I model involving a single convection zone covers 27 hours of stellar evolution, while the first hydrodynamic simulations of a double convection zone, in the Pop III model, span 1.8 hours of stellar life. We find differences between the predictions of mixing length theory and our hydrodynamic simulations. The simulation of the single convection zone in the Pop I model shows a strong growth of the size of the convection zone due to turbulent entrainment. Hence we predict that for the Pop I model a hydrogen injection phase (i.e. hydrogen injection into the helium core) will commence after about 23 days, which should eventually lead to a double convection zone structure known from 1D stellar modeling of low-mass Pop III stars. Our two and three-dimensional hydrodynamic simulations of the double (Pop III) convection zone model show that the velocity field in the convection zones is different from that predicted by stellar evolutionary calculations. The simulations suggest that the double convection zone decays quickly, the flow eventually being dominated by internal gravity waves.

연구 동기 및 목표

• 다차원 시뮬레이션을 활용하여 저질량 항성에서의 核 헬륨 폭발의 유체역학을 연구하고, 1D 별진화 모델의 한계를 해결하고자 한다.
• 1D 모델이 예측한 바와 같이 난류 유입이 금속 농도가 높은 Pop I 항성에서 수소가 헬륨 핵으로 주입되는지 여부를 규명하고자 한다.
• 특히 CNO 연소와 내부 중력파의 역할을 고려하여 금속이 없는 Pop III 항성에서의 이중 대류권의 안정성과 동역학을 조사하고자 한다.
• 혼합길이 이론과 1D 별진화 예측의 정확성을 고해상도 2D 및 3D 유체역학 시뮬레이션과 비교하여 평가하고자 한다.
• 수치 해상도와 핵반응의 영향을 고려하여 핵 헬륨 폭발에서 대류 흐름과 속도 분포의 발달 과정을 평가하고자 한다.

제안 방법

• 유체역학 시뮬레이션은 히라클레스 코드를 사용하여 중력과 핵반응 항목을 포함한 오일러 방정식을 해석하였다.
• 초기 모델은 1.25 M⊙ 금속 농도가 높은 Pop I 항성(Z=0.02)과 0.85 M⊙ 금속이 없는 Pop III 항성(Z=0)에 대해 최신 1D 별진화 코드에서 유도되었다.
• 난류 대류와 유입 과정을 해석하기 위해 리만 해법 기반의 2D 및 3D 시뮬레이션을 수행하였다.
• 시뮬레이션은 Pop I 항성의 경우 27시간, Pop III 항성의 경우 1.8시간 동안 대류권의 진화를 추적하였으며, 흐름 역학과 에너지 생성에 중점을 두었다.
• 1D 예측과의 비교를 위해 속도장과 운동에너지 분석을 수행하였으며, 해상도와 핵반응 활성화에 대한 민감도 테스트도 실시하였다.
• 고해상도 3D 시뮬레이션(heflpopIII.3d)을 저해상도 모델과 비교하여 수렴성과 수치적 오류를 평가하였다.

실험 결과

연구 질문

• RQ1 금속 농도가 높은 Pop I 항성에서 난류 유입이 대류권의 동적 성장을 유도하여 수소가 헬륨 핵으로 주입되는가?
• RQ2 핵 헬륨 폭발의 2D 및 3D 유체역학 시뮬레이션은 대류권 진화와 속도 구조 예측에서 어떻게 다를까?
• RQ3 왜 1D 별진화 모델은 Pop III 항성에서 이중 대류권을 예측하는가? 이는 다차원 유체역학 시뮬레이션에서도 유지되는가?
• RQ4 유체역학 시뮬레이션은 어느 정도 1D 별진화 모델이 예측한 대류 속도를 재현하는가?
• RQ5 핵반응은 대류를 유지하는 데 어떤 역할을 하는가? 만약 반응이 없이 온도 기울기만 존재한다면 대류는 어떻게 발생하는가?

주요 결과

• 1.25 M⊙ 금속 농도가 높은 Pop I 항성에서 난류 유입으로 인해 단일 대류권이 급격히 성장하여 약 23일 후에 수소가 헬륨 핵으로 주입된다.
• Pop I 항성의 3D 유체역학 시뮬레이션은 대류권에서 초기온도 기울기가 초기온도 기울기 이상으로 나타나며, 고해상도 역학과 이전 저해상도 모델보다 향상된 결과를 보인다.
• 0.85 M⊙ 금속이 없는 Pop III 항성에서는 이중 대류권이 빠르게 붕괴되며, 대류 발생 후 흐름이 내부 중력파에 의해 지배된다.
• Pop III 항성의 3D 시뮬레이션은 내부 헬륨 연소 영역에서 1D 예측보다 약 2배 높은 대류 속도를 보였지만, 이론적 기대와는 합리적인 일치를 보였다.
• 핵반응이 없더라도 시뮬레이션에서 대류가 자발적으로 발생함을 확인하여, 초기 온도 기울기만으로도 대류가 유도될 수 있으며, 핵반응은 주로 유지 역할을 한다는 점을 입증하였다.
• 이중 대류권의 2D 및 3D 시뮬레이션 결과, 대류 흐름이 급격히 감쇠됨을 확인하여, 1D 모델이 예측한 이중 대류권의 구조는 다차원 유체역학에서는 안정하지 않음을 시사한다.