[논문 리뷰] Fast ray-tracing algorithm for circumstellar structures (FRACS) I. Algorithm description and parameter-space study for mid-IR interferometry of B[e] stars

G. Niccolini, Philippe Bendjoya외 1명|arXiv (Cornell University)|2010. 09. 28.

Astrophysics and Star Formation Studies참고 문헌 27인용 수 4

한 줄 요약

이 논문은 축대칭 원환성 성간 원반을 매개변수화된 발광도와 대칭성을 사용하여 빠르게 강도 맵과 기울기(visibility)를 계산할 수 있도록 하는 빠른 레이 트레이싱 알고리즘인 FRACS를 소개한다. 중간 적외선 데이터가 기하학적 매개변수(내부 반경, 위치각, 기울기)와 온도 구조(내부 온도 및 경사도)를 잘 제약할 수 있음을 보여주며, 밀도 매개변수는 덜 제약받지만 상한선과 상관관계를 통해 범위를 좁힐 수 있다.

ABSTRACT

The physical interpretation of spectro-interferometric data is strongly model-dependent. On one hand, models involving elaborate radiative transfer solvers are too time consuming in general to perform an automatic fitting procedure and derive astrophysical quantities and their related errors. On the other hand, using simple geometrical models does not give sufficient insights into the physics of the object. We propose to stand in between these two extreme approaches by using a physical but still simple parameterised model for the object under consideration. Based on this philosophy, we developed a numerical tool optimised for mid-infrared (mid-IR) interferometry, the fast ray-tracing algorithm for circumstellar structures (FRACS) which can be used as a stand-alone model, or as an aid for a more advanced physical description or even for elaborating observation strategies. FRACS is based on the ray-tracing technique without scattering, but supplemented with the use of quadtree meshes and the full symmetries of the axisymmetrical problem to significantly decrease the necessary computing time to obtain e.g. monochromatic images and visibilities. We applied FRACS in a theoretical study of the dusty circumstellar environments (CSEs) of B[e] supergiants (sgB[e]) in order to determine which information (physical parameters) can be retrieved from present mid-IR interferometry (flux and visibility). From a set of selected dusty CSE models typical of sgB[e] stars we show that together with the geometrical parameters (position angle, inclination, inner radius), the temperature structure (inner dust temperature and gradient) can be well constrained by the mid-IR data alone. Our results also indicate that the determination of the parameters characterising the CSE density structure is more challenging but, in some cases, upper limits as well as correlations on the parameters characterising the mass loss can be obtained. Good constraints for the sgB[e] central continuum emission (central star and inner gas emissions) can be obtained whenever its contribution to the total mid-IR flux is only as high as a few percents. Ray-tracing parameterised models such as FRACS are thus well adapted to prepare and/or interpret long wavelengths (from mid-IR to radio) observations at present (e.g. VLTI/MIDI) and near-future (e.g. VLTI/MATISSE, ALMA) interferometers.

연구 동기 및 목표

• 복잡한 방사열전달 코드가 너무 느려서 사용이 불가능한 상황에서 중간 적외선 간섭계 데이터를 해석하기 위한 계산적으로 효율적인 방법을 개발하기 위해.
• 간단한 기하학적 모델과 복잡한 자기 일관성 있는 방사열전달 시뮬레이션 사이의 격차를 메우기 위해 매개변수화된 물리적 모델을 사용하기 위해.
• 현재 중간 적외선 간섭계 데이터(예: VLTI/MIDI)로 B[e] 항성 원환성 원반의 어떤 물리적 매개변수를 제약할 수 있는지 규명하기 위해.
• 더 복잡한 몬테카를로 코드를 사용한 세부 모델링을 안내하기 위해 초기 모델 피팅과 매개변수 공간 탐색을 위한 빠르고 영리한 도구를 제공하기 위해.
• MATISSE와 ALMA와 같은 차세대 간섭계를 사용한 관측의 계획 및 해석을 지원하기 위해.

제안 방법

• FRACS는 산산이 흩트리는 것을 고려하지 않고 축대칭성과 완전한 회전 대칭성을 활용하여 계산 부담을 줄인다.
• 원천의 구조에 따라 공간 영역을 적응적으로 해상도를 높이는 4분면 트리 메쉬를 사용하여 효율성을 향상시킨다.
• 원환성 매질의 매개변수화된 발광도 함수로부터 단파장 강도 맵과 기울기 곡선을 계산한다.
• 내부 반경, 온도 경사도, 밀도 구조 등의 물리적 매개변수는 발광도 함수에 포함되어 체계적으로 변화시킨다.
• 중복 계산을 피하기 위해 대칭성을 활용하여 기울기 및 복사율 프로파일 계산을 크게 가속화한다.
• 결과는 몬테카를로 방사열전달 시뮬레이션과 인공 간섭계 데이터와의 비교를 통해 검증되었다.

실험 결과

연구 질문

• RQ1 중간 적외선 간섭계 데이터만으로 B[e] 항성 원환성 원반의 어떤 물리적 매개변수를 신뢰성 있게 복원할 수 있는가?
• RQ2 중간 적외선 기울기 및 복사율 측정을 통해 기하학적 및 열적 매개변수(예: 내부 반경, 온도 경사도)는 얼마나 정확하게 제약될 수 있는가?
• RQ3 밀도 구조와 질량 손실 매개변수는 어느 정도까지 제약될 수 있으며, 상한선이나 상관관계를 통해 어떤 정보를 도출할 수 있는가?
• RQ4 기울기 반응은 특히 내부 원반 영역에서 온도 구조에 따라 어떻게 달라지는가?
• RQ5 FRACS는 더 복잡한 몬테카를로 방사열전달 피팅을 위한 효과적인 초기 모델이 될 수 있는가?

주요 결과

• 중간 적외선 간섭계 데이터는 B[e] 항성 원반의 내부 반경, 위치각, 기울기를 15% 이내의 정확도로 제약할 수 있다.
• 내부 먼지 온도와 온도 경사도(γ)는 각각 20% 이하의 상대 오차로 결정될 수 있다.
• 중앙 소스가 총 중간 적외선 복사율에 몇 퍼센트 정도 기여할 경우, 총 복사율에 대한 기여도는 30% 이내로 제약될 수 있다.
• 밀도 구조 매개변수(n_in, m, A2)는 덜 잘 제약되며, 일부 경우에서 n_in은 최대 약 50%의 정확도로만 도출될 수 있다.
• 외부 온도 경사도(γ')는 신뢰성 있게 결정되지 않으며, 테스트된 구성에서 진짜 값과의 상대적 차이가 28%를 초과한다.
• 대칭성과 적응형 4분면 트리 메쉬를 활용하여 알고리즘은 빠른 계산(매핑당 몇 초)을 달성하여 매개변수 공간 탐색을 효율적으로 가능하게 한다.