[논문 리뷰] Chemical abundances in metal-poor giants: limitations imposed by the use of classical 1D stellar atmosphere models

V. Dobrovolskas, A. Kučinskas외 4명|arXiv (Cornell University)|2010. 10. 12.

Stellar, planetary, and galactic studies참고 문헌 7인용 수 4

한 줄 요약

이 연구는 금속 빈도가 낮은 거성에서 대류 및 온도 비균일성의 영향을 평가하기 위해 3D 수열역학적(CO 5 BOLD) 및 1D 수직평형(LHD) 별대기모델을 비교한다. [M/H] = –3.0에서 중성 원자에 대해 3D–1D 원소 농도 보정치는 –0.8 dex에 이를 수 있으며, 이는 3D 모델에서 수평적 온도 변동이 더 강하고 평균 온도가 낮기 때문에 발생한다. 이는 1D 모델이 저금속성 별에 적용될 때 심각한 한계를 가짐을 시사한다.

ABSTRACT

In this work we have used 3D hydrodynamical (CO5BOLD) and 1D hydrostatic (LHD) stellar atmosphere models to study the importance of convection and horizontal temperature inhomogeneities in stellar abundance work related to late-type giants. We have found that for a number of key elements, such as Na, Mg, Si, Ca, Ti, Fe, Ni, Zn, Ba, Eu, differences in abundances predicted by 3D and 1D models are typically minor (< 0.1 dex) at solar metallicity. However, at [M/H] = -3 they become larger and reach to -0.5...-0.8 dex. In case of neutral atoms and fixed metallicity, the largest abundance differences were obtained for the spectral lines with lowest excitation potential, while for ionized species the largest 3D-1D abundance differences were found for lines of highest excitation potential. The large abundance differences at low metallicity are caused by large horizontal temperature fluctuations and lower mean temperature in the outer layers of the 3D hydrodynamical model compared with its 1D counterpart.

연구 동기 및 목표

• 금속 빈도가 낮은 붉은 거성에서 3D 대류 운동 및 수평적 온도 비균일성이 스펙트럼 선 형성에 미치는 영향을 평가하기 위해.
• 3D 수열역학적(CO 5 BOLD) 및 1D 수직평형(LHD) 별대기 모델에서 유도된 원소 농도의 차이를 정량화하기 위해.
• 특히 중성 및 이온화된 이종에 대해 고금속성 별 농도 분석에서 전통적인 1D 모델의 신뢰성 평가하기 위해.
• Fe, Mg, Ca, Eu와 같은 주요 원소에 대해 3D–1D 농도 보정치가 금속 농도, 이온화 잠재력, 파장에 따라 어떻게 변화하는지 조사하기 위해.
• 3D 및 1D 모델 간 일관된 미세물리학(흡수율, 상태방정식)을 사용함으로써 이전 연구와의 모순을 해결하기 위해.

제안 방법

• 동일한 대기 조건, 화학 조성, 흡수율, 상태방정식을 가진 3D 수열역학적 CO 5 BOLD 모델과 1D 수직평형 LHD 모델을 사용하였다.
• 태양 및 [M/H] = –3.0 금속 농도에서 3D 및 1D 모델 간의 온도 구배 및 수평적 온도 변동을 비교하였다.
• LINFOR3D를 사용하여 LTE 조건 하에서 선 강도를 계산하였으며, 미세난류 효과를 최소화하기 위해 약한 선(W < 5 mÅ)에 집중하였다.
• 고정된 등가폭에서 3D 및 1D 모델 간의 등가폭 유도 농도 차이를 통해 3D–1D 농도 보정치를 계산하였다.
• 이온화 잠재력 0에서 6 eV까지 2 eV 간격으로 λ = 4000 Å 및 8500 Å에서 허구의 선을 합성하여 이온화 잠재력의 영향을 분리하였다.
• Na I, Mg I, Mg II, Si I, Si II, Ca I, Ca II, Ti I, Ti II, Fe I, Fe II, Ni I, Ni II, Zn I, Zn II, Ba II, Eu II 등 다양한 원소에 대해 보정치를 분석하였다.

실험 결과

연구 질문

• RQ1 금속 빈도가 낮은 거성에서 3D 수열역학적 모델은 1D 모델에 비해 예측하는 원소 농도에 어떤 영향을 미치는가?
• RQ2 중성 및 이온화된 이종에 대해 3D–1D 농도 보정치가 금속 농도, 이온화 잠재력, 파장에 따라 어떻게 달라지는가?
• RQ3 왜 3D와 1D 모델 간의 농도 차이가 저금속성에서 더 크며, 이러한 차이를 일으키는 물리적 메커니즘은 무엇인가?
• RQ4 외부 대기층에서의 수평적 온도 변동 및 평균 온도 차이가 선 강도 및 농도 결정에 어떤 영향을 미치는가?
• RQ5 연구 간의 3D–1D 보정치 차이를 설명하는 데 있어 모델의 미세물리학(예: 흡수율, 상태방정식) 차이가 어느 정도 기여하는가?

주요 결과

• [M/H] = –3.0에서 중성 원자에 대한 3D–1D 농도 보정치는 최대 –0.80 dex에 이르며, 주로 λ = 8500 Å에서 저이온화 잠재력의 선에서 발생한다.
• 일차 이온화 이종의 경우 3D–1D 농도 보정치는 더 작은 편이며, [M/H] = –3.0에서도 일반적으로 ±0.1 dex를 초과하지 않는다.
• 중성 이종에 대해 가장 큰 3D–1D 농도 보정치는 저이온화 잠재력(예: 0 eV)에서 발생하며, 이는 대기 외부의 더 차가우며 난류가 심한 층에서 선이 형성되기 때문이다.
• 이온화 이종의 경우 가장 음의 보정치는 고이온화 잠재력에서 나타나며, 이는 깊은 층(로그 τ_Ross ≈ 0)에서 온도 비균일성이 더 강하게 작용함을 시사한다.
• 이 연구와 이전 연구([4]) 간의 모순은 [4]에서 사용된 3D–1D 모델 간 외부 층에서 더 큰 온도 차이가 있었기 때문이며, 이는 미세물리학의 비일관성에서 기인한다.
• 이 연구는 3D 모델에서 수평적 온도 변동과 낮은 평균 온도가 저금속성에서 원소 농도의 큰 차이를 유발하는 주요 원인임을 확인한다. 특히 중성 원자에 대해 그렇다.