[논문 리뷰] An accelerated Monte Carlo method to solve two-dimensional radiative transfer and molecular excitation

M. R. Hogerheijde, F. van der Tak|arXiv (Cornell University)|2000. 08. 10.

Atmospheric Ozone and Climate참고 문헌 1인용 수 12

한 줄 요약

이 논문은 축대칭 천체물리 모델에서 2차원 복사전달과 분자 진동을 해결하기 위한 가속화된 몬테카를로 방법을 제시한다. 셀 기반의 수식을 사용하여 고광학 깊이(τ ~ 10^4)에서 국소 및 외부 복사장 간의 분리를 실현함으로써 정확도와 속도를 향상시킨다. 이 방법은 문제에 특화된 캘리브레이션 없이 별 형성 환경에서 분자 선 방출과-dust 연속 방출을 효율적이고도 영리하게 모델링할 수 있게 한다.

ABSTRACT

We present a numerical method and computer code to calculate the radiative transfer and excitation of molecular lines. Formulating the Monte Carlo method from the viewpoint of cells rather than photons allows us to separate local and external contributions to the radiation field. This separation is critical to accurate and fast performance at high optical depths (tau>100). The random nature of the Monte Carlo method serves to verify the independence of the solution to the angular, spatial, and frequency sampling of the radiation field. These features allow use of our method in a wide variety of astrophysical problems without specific adaptations: in any axially symmetric source model and for all atoms or molecules for which collisional rate coefficients are available. Continuum emission and absorption by dust is explicitly taken into account but scattering is neglected. We illustrate these features in calculations of (i) the HCO+ J=1-0 and 3-2 emission from a flattened protostellar envelope with infall and rotation, (ii) the CO, HCO+, CN and HCN emission from a protoplanetary disk and (iii) HCN emission from a high-mass young stellar object, where infrared pumping is important. The program can be used for optical depths up to 1000-10,000, depending on source model. We expect this program to be an important tool in analysing data from present and future infrared and (sub) millimetre telescopes.

연구 동기 및 목표

• 고광학 깊이(τ ≥ 100)에서 전통적인 몬테카를로 방법의 계산 비효율성을 해결하여 실제 천체물리 시뮬레이션에의 활용을 제한하는 문제를 해결한다.
• 축대칭 기하학에서 분자 선 방출과 먼지 연속 방출을 포함한 다양한 천체물리 문제를 다룰 수 있는 영리하고 일반적인 수치 도구를 개발한다.
• 형성 중인 별주위 환경, 행성계 형성 원반, 고질량 젊은 별주위 환경에서 복사전달과 비-LTE 분자 진동을 정확하게 모델링할 수 있도록 한다.
• 기존 몬테카를로 코드의 한계를 극복하기 위해 국소 복사장 해법과 전반적 전달을 분리함으로써 수렴성과 안정성을 향상시킨다.
• 향후 인프라적 및 밀리미터파 망원경(예: ALMA, SOFIA)의 관측 데이터 해석을 지원하기 위해 견고하고 재사용 가능한 시뮬레이션 프레임워크를 제공한다.

제안 방법

• 개별 광자를 중심으로 하는 몬테카를로 방법이 아닌 공간 셀의 관점에서 수식을 재구성함으로써 국소 및 외부 복사장 기여를 분리할 수 있도록 한다.
• 이중 단계 반복 알고리즘을 구현한다: 첫 번째 단계는 고정된 각도 및 주파수 샘플링을 사용하여 노이즈를 최소화하면서 인구 분포에 수렴한다. 두 번째 단계는 오차 기준을 초과하는 셀에서 레이 수를 두 배로 늘여 정확도를 향상시킨다.
• 가속화된 라미다 반복(ACI) 기반 기법을 사용하여 국소 진동 및 복사장 해를 자성적으로 해결함으로써 고광학 깊이에서 수렴성을 향상시킨다.
• 먼지 연속 방출 및 흡수를 명시적으로 포함하지만 산란은 무시하여 계산 복잡도를 감소시킨다.
• 최종 해법에 광선 추적을 적용하여 주어진 소스 기울기 및 거리에 따른 천구 평균 강도 분포를 계산한다.
• 동일한 충돌률 계수와 기하학을 사용하여 베른스의 CO 구름 및 최치 등(1995)의 B335 모델을 기준으로 방법을 검증한다.

실험 결과

연구 질문

• RQ1 셀 기반 몬테카를로 수식은 고광학 깊이(τ ≥ 100)에서 복사전달 및 분자 진동 계산의 정확도와 효율성을 향상시키는가?
• RQ2 구형 대칭, 고광학 깊이의 구름에서 CO 및 CS 선 진동의 기존 기준 결과를 얼마나 잘 재현하는가?
• RQ3 유동, 도는 형성 중인 별주위 환경 및 행성계 형성 원반과 같은 복잡한 축대칭 소스 구조를 얼마나 잘 처리할 수 있는가?
• RQ4 국소 및 외부 복사장 기여를 분리함으로써 고다이내믹 레인지 시뮬레이션에서 더 빠른 수렴성과 낮은 노이즈를 달성하는가?
• RQ5 고질량 젊은 별주위 환경에서 복사 포획이 강하게 영향을 미치는 조건에서, 이 방법은 적외선 펌프 효과를 얼마나 정확하게 모델링할 수 있는가?

주요 결과

• 코드는 소스 구조에 따라 광학 깊이 10^3~10^4까지 복사전달 및 분자 진동을 문제에 특화된 캘리브레이션 없이 성공적으로 모델링한다.
• 원래 계산과 현재 계산의 정확도 범위 내에서 베른스(1979)의 CO J=1–0 및 J=2–1 진동 온도를 재현하여 복사 포획 처리의 정확성을 검증한다.
• B335 모델에 대해 코드는 최치 등(1995)의 결과와 선형형태에서 10% 이내로 일치시키며, 속도장 및 광학 깊이 효과 처리의 정확성을 확인한다.
• H2CO 선 진동의 불일치는 코드 오류가 아니라 충돌률 계수의 차이로 기인하며, 탈출 선형형태에서 약 10% 수준의 편차가 일관되게 나타난다.
• 실린더 대칭 구현은 구형 대칭 버전과 일관된 결과를 도출하며, 미세한 수치적 진동은 실린더 격자에서 기하학적 이산화로 기인한 것으로 기인한다.
• 기준 CO 구름에 대해 표준 워크스테이션에서 약 5분 내로 코드가 실행되어 복잡한 천체물리 모델에 대한 계산 가능성을 입증한다.