[논문 리뷰] A code for optically thick and hot photoionized media

A. M. Dumont, A. Abrassart외 1명|arXiv (Cornell University)|2000. 03. 15.

Nuclear reactor physics and engineering인용 수 3

한 줄 요약

이 논문은 딱딱한 X선 연속 스펙트럼에 의해 조사되는 뜨겁고 광학적으로 두꺼운 매체를 위한 새로운 광이온화 코드인 titan을 제시한다. 이 코드는 연속과 선에 대한 전체 복사전달을 이중 흐름 근사법을 사용하여 해결하며 국소 탈출 확률 근사법을 피하고, Compton 산란을 정확히 모델링하기 위해 몬테카를로 코드(noar)와 결합한다. 이로 인해 기존의 cloudy와 같은 코드보다도 고 Thomson 광학 두께 조건에서 더 정확한 열, 이온화 및 스펙트럼 구조를 제공한다.

ABSTRACT

We describe a code designed for hot media {(T $\ge$} a few 10$^4$ K), optically thick to Compton scattering. It computes the structure of a plane-parallel slab of gas in thermal and ionization equilibrium, illuminated on one or on both sides by a given spectrum. Contrary to the other photoionization codes, it solves the transfer of the continuum and of the lines in a two stream approximation, without using the local escape probability formalism to approximate the line transfer. We stress the importance of taking into account the returning flux even for small column densities (10$^{22}$ cm$^{-2}$), and we show that the escape probability approximation can lead to strong errors in the thermal and ionization structure, as well as in the emitted spectrum, for a Thomson thickness larger than a few tenths. The transfer code is coupled with a Monte Carlo code which allows to take into account Compton and inverse Compton diffusions, and to compute the spectrum emitted up to MeV energies, in any geometry. Comparisons with {\sc{cloudy}} show that it gives similar results for small column densities. Several applications are mentioned.

연구 동기 및 목표

• 디지털 X선 연속 스펙트럼에 의해 조사되는 뜨겁고 광학적으로 두꺼운 매체에 대해 정확한 광이온화 코드의 부족을 해결한다.
• 두꺼운 매체에서 붕괴하는 국소 탈출 확률 근사법을 사용하는 기존 코드의 한계를 극복한다.
• 완전한 복사전달을 고려한 평면평행 판에서 에너지 균형, 이온화 및 통계적 평형의 자기 일관성 있는 처리를 개발한다.
• 몬테카를로 Compton 산란 코드와의 결합을 통해 MeV 에너지까지 반사 및 방출 스펙트럼을 정확하게 계산할 수 있도록 한다.
• AGN, 순환 디스크, X선 이진성과 같은 다양한 천체물리적 환경을 모델링하기 위한 유연하고 모듈식 프레임워크(titan, noar, pegas, iris)를 제공한다.

제안 방법

• 연속과 선에 대한 복사전달을 국소 탈출 확률 형식을 피하기 위해 이중 흐름 근사법을 사용하여 해결한다.
• 에너지, 이온화 및 통계적 평형 방정식의 결합 해법을 사용하여 평면평행 판 기하에서 열 및 이온화 구조를 계산한다.
• 직접 및 역방향 Compton 산란을 모델링하기 위해 몬테카를로 방법(noar)을 사용하여 MeV 에너지까지 정확한 스펙트럼 계산을 가능하게 한다.
• 모듈식 코드 구조를 활용: titan은 구조 및 전달를 담당하고, noar는 Compton 산란 및 스펙트럼 합성을 담당하며, pegas는 얇은 매체를 위해, iris는 고해상도 선 농도를 위해 사용된다.
• 최신 원자 데이터를 통합하고, L-껍질 철 선과 추가 원소를 포함한 다수준 원자 방정식을 해결한다.
• 표준 광이온화 코드에서 무시되는 반사되는 복사 플럭스를 일관되게 처리하며, 이는 정확한 결과에 매우 중요하다는 것이 입증되었다.

실험 결과

연구 질문

• RQ1 광학적으로 두꺼운 뜨거운 매체에서 광이온화 코드에서 반사되는 플럭스를 생략할 경우 열 및 이온화 구조에 어떤 영향을 미치는가?
• RQ2 Compton 두꺼운 매체에서 국소 탈출 확률 근사법이 얼마나 붕괴되는가? 이로 인해 발생하는 이온화 및 스펙트럼 구조 오차는 어느 정도인가?
• RQ3 이중 흐름 복사전달 접근법과 몬테카를로 Compton 산란을 조합하면, 기존의 탈출 확률 기반 코드인 cloudy와 비교해 두꺼운 뜨거운 조사된 매체에서 더 정확한 결과를 얻을 수 있는가?
• RQ4 동일한 물리적 설정에서 신규 코드(titan/noar)와 기존 모델(Zycki 등 1994) 간의 열 및 이온화 구조 및 반사 스펙트럼은 어떻게 비교되는가?
• RQ5 이러한 복잡한 광이온화 코드에서 수렴 속도 향상 및 원자 데이터 처리 향상을 위해 어떤 개선이 필요한가?

주요 결과

• 조직 두께가 10^22 cm^-2이상일지라도 반사 플럭스를 포함시켜야 하며, 이를 생략할 경우 열 및 이온화 구조에 심각한 오차가 발생한다.
• Thomson 광학 두께가 몇 분의 일 이상을 초과할 경우 탈출 확률 근사법이 열 및 이온화 구조와 방출 스펙트럼에 강력한 오차를 유발한다.
• ξ ≈ 200 erg cm s^-1인 판에서, 신규 코드의 스펙트럼에서는 OVIII 흡수 에지가 발산으로 나타나지만, Zycki 등(1994)의 모델에서는 흡수로 나타나며, 이는 표면 온도의 차이와 더 깊은 층의 냉각 차이 때문이므로, 열 구조의 정확성이 매우 중요하다.
• titan가 계산한 온도 프로파일은 Zycki 등(1994)의 모델보다 높은 표면 온도를 보이며, τ_T ≈ 0.7에서 블랙바디 온도에 도달하고, τ_T > 0.7에서는 그 이하로 떨어진다.
• titan와 Zycki 등(1994)이 계산한 반사 스펙트럼은 전반적인 형태는 유사하지만, OVIII 에지의 격차는 정확한 열 구조의 중요성을 부각시킨다.
• 전체 구조 및 스펙트럼에 대해 합리적인 시간 내에 수렴을 달성하였으며, 향후 주요 개선 과제로 Accelerated Lambda Iteration 방법을 통한 수렴 가속화가 필요하다고 확인되었다.