[논문 리뷰] Three-dimensional surface convection simulations of metal-poor stars -- The effect of scattering on the photospheric temperature stratification

R. Collet, Wolfgang Hayek외 4명|arXiv (Cornell University)|2011. 01. 17.

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

한 줄 요약

이 연구는 Bifrost 코드를 사용한 극저금속 항성의 3차원 복사유체역학 시뮬레이션에서 산산각산 처리의 영향을 조사한다. 세 가지 방법을 비교한다: 완전한 비편향 산산각산, 투명도가 낮은 층에서 산산각산을 무시하는 것, 산산각산을 순수 흡수로 취급하는 것. 주요 발견은 투명도가 낮은 층에서 산산각산 기여를 무시할 경우, 전체 비편향 산산각산 해법과 거의 동일한 온도 구조를 얻는다는 것이다. 반면 산산각산을 흡수로 취급할 경우 더 높은 온도, 더 낮은 정확도의 결과를 얻으며, 이는 이전의 냉각된 온도 구조가 산산각산 근사의 산물이 아니라는 것을 확인한다.

ABSTRACT

Three-dimensional (3D) radiative hydrodynamic model atmospheres of metal-poor late-type stars are characterized by cooler upper photospheric layers than their 1D counterparts. This property of 3D models can dramatically affect elemental abundances derived from temperature-sensitive spectral lines. We investigate whether the cool surface temperatures predicted by metal-poor 3D models can be ascribed to the approximated treatment of scattering in the radiative transfer. We use the Bifrost code to test three different ways to handle scattering in 3D model atmospheres of metal-poor stars. First, we solve self-consistently the radiative transfer equation for a source function with a coherent scattering term. Second, we solve the radiative transfer equation for a Planckian source function, neglecting the contribution of continuum scattering to extinction in the optically thin layers; this has been the default mode in previous models of ours. Third, we treat scattering as pure absorption everywhere, which is the standard case in CO5BOLD models. We find that the second approach produces temperature structures with cool upper photospheric layers very similar to the correct coherent scattering solution. In contrast, treating scattering as pure absorption leads to significantly hotter and shallower temperature stratifications. The main differences in temperature structure between our published models and those generated with the CO5BOLD code can be traced to the different treatments of scattering. Neglecting the contribution of continuum scattering to extinction in optically thin layers provides a good approximation to the full radiative transfer solution for metal-poor stars. Our results demonstrate that the cool temperature stratifications predicted for metal-poor late-type stellar atmospheres by previous models of ours are not an artifact of the approximated treatment of scattering.

연구 동기 및 목표

• 극저금속 항성의 3차원 대기 모델에서 관측된 냉각된 상부 광구 온도 구조가 산산각산 근사의 산물인지 여부를 확인하기 위해
• 복사전달에서 연속 산산각산의 세 가지 다른 처리 방법을 비교하기 위해: 완전한 비편향 산산각산, 투명도가 낮은 층에서 산산각산를 흡수에서 제외하는 것, 산산각산를 순수 흡수로 취급하는 것
• 이러한 방법들이 극저금속 항성의 3차원 복사유체역학 시뮬레이션에서 물리적·계산적 트레이드오프를 어떻게 초래하는지 평가하기 위해
• 이전 3차원 모델 결과의 강건성을 검증하기 위해, 냉각된 온도 구조가 물리적 원인인지 산산각산 근사의 결과인지 테스트하기 위해

제안 방법

• Bifrost 코드를 사용하여 극저금속 항성(두 개의 거대항성과 한 개의 전환점 항성)의 3차원 복사유체역학 시뮬레이션을 수행하였다.
• 세 가지 복사전달 처리 방법을 구현하였다: (1) 반복적 해법을 사용한 비편향 등방성 산산각산, (2) 산산각산를 투명도가 낮은 층에서 흡수에서 제외한 플랑크 광원 함수, (3) 산산각산를 순수 흡수로 취급한 플랑크 광원 함수.
• 산산각산를 흡수로 취급하는 경우 실제 조건을 반영하기 위해 평균 구조를 반복 조정하는 투과도 분할 기법을 사용하였다.
• 세 방법 간의 온도 및 밀도 구조를 비교하여 정확도와 계산 비용을 평가하였다.
• 상부 광구 온도를 조절하는 복사전달의 역할을 등온 조건과 비교하여 평가하였다.
• 온도 구조의 정확성에 기반해 스펙트럼 선 형성 예측을 암시적으로 평가하였다.

실험 결과

연구 질문

• RQ1 투명도가 낮은 층에서 산산각산 기여를 무시하는 것이 극저금속 항성의 3차원 시뮬레이션에서 전체 비편향 산산각산 해법에 대한 타당한 근사치를 제공하는가?
• RQ2 산산각산를 순수 흡수로 취급할 경우, 극저금속 항성의 상부 광구에서 예측되는 온도 구조에 어떤 영향을 미치는가?
• RQ3 이전 3차원 모델에서 관측된 냉각된 온도 구조는 산산각산 근사의 산물인지, 아니면 대류의 물리적 결과인지?
• RQ4 완전한 비편향 산산각산와 간단한 근사치 간의 계산 비용 차이는 무엇인가?
• RQ5 상부 광구의 온도 구조를 조절하는 데 복사전달이 등온 과정보다 얼마나 중요한가?

주요 결과

• 투명도가 낮은 층에서 산산각산를 제외하는 근사치는 전체 반복적 해법을 사용한 비편향 산산각산와 거의 구분되지 않는 온도 및 밀도 구조를 제공한다.
• 산산각산를 순수 흡수로 취급할 경우, 상부 광구에서 상당히 더 높고 얕은 온도 구조를 유도하며, [Fe/H] = -3 전환점 항성의 경우 최대 약 200 K의 편차가 발생한다.
• 따라서 이전 3차원 모델에서 관측된 냉각된 상부 광구 온도는 산산각산 근사의 산물이 아니며, 3차원 대류의 강력한 물리적 결과임을 확인한다.
• 산산각산를 얇은 층에서 무시하는 근사치는 정확도와 계산 효율성 사이의 좋은 균형을 제공하여, 낮은 비용으로 고정밀 시뮬레이션을 가능하게 한다.
• 산산각산가 없더라도 복사전달은 여전히 상부 광구 온도를 조절하는 데 필수적이며, 이는 등온 값보다 상당히 높은 온도를 유지함을 의미한다.
• 산산각산를 흡수로 취급하는 방법은 물리적으로 비현실적이며, 특히 [Fe/H] = -3 및 [Fe/H] = -2 금속성에서 온도 구조에 체계적인 오류를 유발한다.