[논문 리뷰] 3-D radiative transfer in clumped hot star winds I. Influence of clumping on the resonance line formation

B. Šurlan, W.‐R. Hamann외 3명|arXiv (Cornell University)|2012. 02. 21.

Astrophysics and Star Formation Studies참고 문헌 42인용 수 47

한 줄 요약

이 연구는 구름으로 이루어진 뜨거운 항성의 흡수선 형성에 대해 처음으로 전체 3차원 몽테카를로 복사전달 모델을 제시한다. 이 모델은 현실적인 3차원 밀도 구조, 비단조선형 속도장, 그리고 구름 간 매질(ICM)을 고려한다. 연구 결과, 매크로클러스터링은 특히 구름 간 거리가 클수록 선 투과도를 감소시키며, 이는 구름 구조를 忽시할 경우 질량 손실률이 체계적으로 과소평가됨을 의미한다. 따라서 이는 ρ 기반(예: Hα) 및 ρ² 기반(예: 자외선 흡수선) 진단법 간의 모순을 해결한다.

ABSTRACT

The true mass-loss rates from massive stars are important for many branches of astrophysics. For the correct modeling of the resonance lines, which are among the key diagnostics of stellar mass-loss, the stellar wind clumping turned out to be very important. In order to incorporate clumping into radiative transfer calculation, 3-D models are required. Various properties of the clumps may have strong impact on the resonance line formation and, therefore, on the determination of empirical mass-loss rates. We incorporate the 3-D nature of the stellar wind clumping into radiative transfer calculations and investigate how different model parameters influence the resonance line formation. We develop a full 3-D Monte Carlo radiative transfer code for inhomogeneous expanding stellar winds. The number density of clumps follows the mass conservation. For the first time, realistic 3-D models that describe the dense as well as the tenuous wind components are used to model the formation of resonance lines in a clumped stellar wind. At the same time, non-monotonic velocity fields are accounted for. The 3-D density and velocity wind inhomogeneities show very strong impact on the resonance line formation. The different parameters describing the clumping and the velocity field results in different line strengths and profiles. We present a set of representative models for various sets of model parameters and investigate how the resonance lines are affected. Our 3-D models show that the line opacity is reduced for larger clump separation and for more shallow velocity gradients within the clumps. Our new model demonstrates that to obtain empirically correct mass-loss rates from the UV resonance lines, the wind clumping and its 3-D nature must be taken into account.

연구 동기 및 목표

• 1D 또는 2D 근사의 한계를 넘어 현실적인 항성 풍속 클러스터링을 모델링할 수 있는 전체 3차원 복사전달 프레임워크를 개발하는 것.
• 3차원 밀도 비균일성, 비단조선형 속도장, 그리고 구름 간 매질(ICM)이 뜨거운 항성 풍속에서 흡수선 형성에 미치는 영향을 조사하는 것.
• ρ 기반(예: Hα) 및 ρ² 기반(예: 자외선 흡수선) 진단법으로 유도된 질량 손실률 간 오랜 기간 지속된 모순을 현실적인 클러스터링 효과를 통합함으로써 해결하는 것.
• 클러스터링 매개변수(예: 구름 간 거리, 시작 반경, 속도 분산)가 P-Cygni 프로파일 형상과 선 강도에 미치는 영향을 정량화하는 것.
• 향후 관측 스펙트럼 분석을 위한 융통성 있고 물리적으로 일관된 모델링 도구를 제공하여 실증적 질량 손실률 추정을 향상시키는 것.

제안 방법

• 임의의 3차원 클러스터 기하구조를 가진 비균일하고 팽창하는 항성 풍속을 통과하는 광자의 전파를 시뮬레이션하는 전체 3차원 몽테카를로 복사전달 코드를 개발하였다.
• 질량 보존 원리를 적용하여 클러스터의 수밀도가 매끄러운 풍속과 동일한 질량 손실률을 따르도록 하였으며, 클러스터 간 거리는 매개변수 $ L_0 $ 로 제어하였다.
• 속도 법에 매개변수 $ m $ 를 도입하여 클러스터 내 비단조선형 속도장과 속도 분산을 구현하였다.
• 구름 간 매질(ICM)을 고려하기 위해 공백이 아닌 배경을 도입하고, 투과도 $ d $ 를 설정함으로써 현실적인 선 포화 행동을 가능하게 하였다.
• 초기형 항성의 자외선 스펙트럼에서 핵심 진단 도구인 双선선 프로파일(예: $ ext{C}^{4+} $, $ ext{Si}^{4+} $)을 사용하여 모델을 시험하였다.
• 모델 간 비교를 위해 $ D $, $ L_0 $, $ r_{\text{cl}} $, $ m $, $ d $ 를 다양하게 설정하였으며, 일관성을 확보하기 위해 기본 매개변수 세트(표 1)를 사용하였다.

실험 결과

연구 질문

• RQ1 3차원 풍속 클러스터링은 뜨거운 항성에서 자외선 흡수선 형성에 어떤 영향을 미치는가? 특히 선 투과도와 프로파일 형상 측면에서.
• RQ2 구름 간 거리($ L_0 $)와 클러스터링 시작 반경($ r_{\text{cl}} $)은 P-Cygni 프로파일의 강도와 구조에 어떤 영향을 미치는가?
• RQ3 클러스터 내 비단조선형 속도장과 속도 분산은 흡수선의 파랑 끝부분과 유도된 최종 속도($ v_{\infty} $)에 어떤 영향을 미치는가?
• RQ4 공백이 아닌 구름 간 매질(ICM)이 존재할 경우, 선 포화 및 포화/비포화 선 간 일관성에 어떤 영향을 미치는가?
• RQ5 매크로클러스터링 효과는 ρ 기반 및 ρ² 기반 진단법으로 유도된 질량 손실률 간 관측된 모순을 설명할 수 있는가?

주요 결과

• 매크로클러스터링은 특히 큰 구름 간 거리($ L_0 $)일수록 효과적인 선 투과도를 감소시키며, 이는 $ v_{\infty} $ 근처에서 특징적인 흡수 균열이 나타나는 약한 흡수를 초래한다.
• 고정된 클러스터링 인자 $ D $ 를 가정할 경우, 투과도 감소 효과는 $ L_0 $ 가 클수록 가장 뚜렷하며, 이는 구름을 忽시할 경우 자외선 흡수선에서 유도된 실증적 질량 손실률이 체계적으로 과소평가됨을 의미한다.
• 클러스터링의 시작($ r_{\text{cl}} $)은 선 프로파일에 강한 영향을 미친다: 조기에 시작할수록 선 중심부에 더 두드러진 흡수 균열이 나타나며, 이는 클러스터링 시작을 진단하는 데 유용할 수 있다.
• 투과도 $ d = 0.05 $ 인 공백이 아닌 구름 간 매질(ICM)이 반드시 필요하며, 관측 결과에서 요구하는 포화 및 비포화 선을 동시에 재현할 수 있다.
• 클러스터 내 속도 분산은 $ m = 0.2 $ 로 모델링되었으며, 이는 흡수의 파랑 끝부분을 $ v_{\infty} $ 초과 영역으로 연장시킨다. 이는 비단조선형 유동을 고려하지 않을 경우 최종 속도가 과대평가될 수 있음을 시사한다.
• 밀도 비균일성과 비단조선형 속도장의 상호작용은 실제 클러스터링 풍속에서 본질적으로 연관되어 있으므로, 이 둘을 동시에 모델링해야 한다.