[논문 리뷰] $\mathit{Chempy}$: A flexible chemical evolution model for abundance fitting - Do the Sun's abundances alone constrain chemical evolution models?

J. Rybizki, A. Just외 1명|arXiv (Cornell University)|2017. 02. 28.

Stellar, planetary, and galactic studies참고 문헌 98인용 수 16

한 줄 요약

Chempy는 유연한 베이지안 화학진화 모델로, 매개변수화된 일점 모델 프레임워크를 사용해 별의 원소 농도에 맞추며, 모델 파라미터와 핵합성 수확 집합에 대한 전체 사후 분포 샘플링을 가능하게 한다. 이 모델은 IMF의 고질량 기울기를 α_IMF = −2.42⁺⁰.⁰⁶₋₀.⁰⁵로 제약하며(11.6⁺².¹₋₁.₆%의 IMF가 중심붕괴 초신성로 폭발함), 초신성 Ia의 발생률을 1000 M⊙당 0.5–1.4로 제약하면서도, 여전히 존재하는 농도 불일치를 드러내어 누락된 핵합성 채널을 시사한다.

ABSTRACT

Elemental abundances of stars are the result of the complex enrichment history of their galaxy. Interpretation of observed abundances requires flexible modeling tools to explore and quantify the information about Galactic chemical evolution (GCE) stored in such data. Here we present Chempy, a newly developed code for GCE modeling, representing a parametrized open one-zone model within a Bayesian framework. A Chempy model is specified by a set of 5-10 parameters that describe the effective galaxy evolution along with the stellar and star-formation physics: e.g. the star-formation history, the feedback efficiency, the stellar initial mass function (IMF) and the incidence of supernova type Ia (SN Ia). Unlike established approaches, Chempy can sample the posterior probability distribution in the full model parameter space and test data-model matches for different nucleosynthetic yield sets. We extend Chempy to a multi-zone scheme. As an illustrative application, we show that interesting parameter constraints result from only the ages and elemental abundances of Sun, Arcturus and the present-day interstellar medium (ISM). For the first time, we use such information to infer IMF parameter via GCE modeling, where we properly marginalize over nuisance parameters and account for different yield sets. We find that of the IMF $11.6_{-1.6}^{+2.1}$ % explodes as core-collapse SN, compatible with Salpeter 1955. We also constrain the incidence of SN Ia per 10^3 Msun to 0.5-1.4. At the same time, this Chempy application shows persistent discrepancies between predicted and observed abundances for some elements, irrespective of the chosen yield set. These cannot be remedied by any variations of Chempy's parameters and could be an indication for missing nucleosynthetic channels. Chempy should be a powerful tool to confront predictions from stellar nucleosynthesis with far more complex abundance data sets.

연구 동기 및 목표

• 복잡한 별 농도 데이터를 베이지안 프레임워크 내에서 맞출 수 있는 유연하고 오픈소스의 화학진화 모델을 개발하는 것.
• 은하 형성력, IMF, 피드백 효율, 초신성 Ia 발생 빈도와 같은 모델 파라미터에 대한 전체 사후 분포 샘플링을 가능하게 하는 것.
• 다양한 핵합성 수확 집합을 동시에 테스트하고 비교하여 불확실성과 체계적 오차를 고려하는 것.
• 태양, 아크투루스, 국지 B형 항성의 관측 농도로부터 IMF의 고질량 기울기를 제약하며, 부수적 파라미터를 통합하여 분석하는 것.
• 예측된 농도와 관측된 농도 사이의 지속적인 불일치를 식별하여, 누락된 핵합성 과정이 존재할 가능성을 밝혀내는 것.

제안 방법

• Chempy는 은하 형성력, 피드백, IMF, 초신성 Ia 발생 빈도를 묘사하는 5–10개 자유도 파라미터를 가진 매개변수화된 일점 화학진화 모델을 사용한다.
• 모델 파라미터에 대한 전체 사후 분포 샘플링을 가능하게 하여 불확실성 정량화와 통합을 가능하게 하는 베이지안 추론 프레임워크를 사용한다.
• 다양한 핵합성 수확 집합을 지원하여, 관측된 농도를 재현하는 데 있어 각 집합의 능력을 직접 비교할 수 있다.
• 다중영역 확장은 서로 다른 항성에 대해 독립적인 은하간 매질(ISM) 진화 역사 모델링을 가능하게 하여 공통 물리적 파라미터에 대한 제약 강화를 이룬다.
• 항성 연령 정보를 통합하고 농도 측정값과 수확 불확실성 간의 공분산을 고려한다.
• 모델은 공개되어 있으며, SSP 피드백 출력을 통해 n-체 시뮬레이션과 통합될 수 있다.

실험 결과

연구 질문

• RQ1 적절히 부수적 파라미터와 수확 불확실성을 통합할 때, 태양의 원소 농도만으로 초기 질량 함수(IMF)의 고질량 기울기를 제약할 수 있는가?
• RQ2 1000 M⊙의 별 형성당 초신성 Ia의 발생률은 얼마이며, 수확 집합에 따라 어떻게 달라지는가?
• RQ3 태양, 아크투루스, 국지 B형 항성의 관측 농도는 은하수의 은하 형성력과 피드백 효율을 어느 정도 제약하는가?
• RQ4 모델 파라미터를 최적화한 후에도 나트륨과 같은 특정 원소에 대해 여전히 예측 농도와 관측 농도 사이의 불일치가 지속되는 이유는 무엇인가?
• RQ5 Chempy의 다중영역 확장은 두꺼운 디스크와 얇은 디스크와 같은 서로 다른 항성 집단에 대해 별개의 화학적 풍화 역사를 드러낼 수 있는가?

주요 결과

• IMF의 고질량 기울기는 α_IMF = −2.42⁺⁰.⁰⁶₋₀.⁰⁵로 제약되며, 이는 IMF의 11.6⁺².¹₋₁.₆%가 중심붕괴 초신성로 폭발함을 의미한다.
• 초신성 Ia의 발생률은 1000 M⊙의 별 형성당 0.5–1.4건으로 제약되며, 이는 이전 메타분석과 일치하지만 더 정밀한 결과를 제공한다.
• 초신성 Ia의 지연 시간은 3–6 Gyr로 추정되어, 아크투루스의 α-과잉을 설명하기 위해 장수명의 원천 채널이 필요함을 시사한다.
• 모델 파라미터와 수확 집합의 철저한 설정에도 불구하고 나트륨은 일관되게 과잉 예측되며, 이는 누락된 핵합성 채널이나 비균일한 혼합이 존재할 가능성을 시사한다.
• 다중영역 확장은 태양, 아크투루스, B형 항성에 대해 별개의 ISM 풍화 역사를 성공적으로 회복하여, 이들이 서로 다른 은하 성분과 연관되어 있음을 지지한다.
• Chempy는 최소한의 데이터로부터 IMF와 초신성 Ia 파라미터에 대한 강력하고 불확실성 인식 기반 추론을 가능하게 하여, 농도 조사 해석에 있어 그 힘을 입증한다.