[논문 리뷰] A fast tree-based method for estimating column densities in Adaptive Mesh Refinement codes Influence of UV radiation field on the structure of molecular clouds

Valeska Valdivia, P. Hennebelle|arXiv (Cornell University)|2014. 10. 02.

Astrophysics and Star Formation Studies참고 문헌 33인용 수 6

한 줄 요약

이 논문은 적응 격자 분해(AMR) 코드에서 열량 밀도를 추정하기 위한 빠르고 트리 기반의 알고리즘을 제시하며, 분자의 구름 시뮬레이션 중 실시간 방사선 전달 계산을 가능하게 한다. 이 방법은 상호작용이 최소화된 병렬 처리를 통해 계산 비용을 O(N log N)으로 줄이며, 흡수율 오차는 10% 미만, 열량 밀도 추정 오차는 50% 미만을 달성한다. 이는 자외선 스크리닝이 밀도가 높은 가스를 선호적으로 냉각시켜, 잼 질량을 낮추고 분열을 증가시켜 시뮬레이션에서 더 밀도가 높고 컴팩트한 덩어리의 형성을 이끈다.

ABSTRACT

Context. Ultraviolet radiation plays a crucial role in molecular clouds. Radiation and matter are tightly coupled and their interplay influences the physical and chemical properties of gas. In particular, modeling the radiation propagation requires calculating column densities, which can be numerically expensive in high-resolution multidimensional simulations. Aims. Developing fast methods for estimating column densities is mandatory if we are interested in the dynamical influence of the radiative transfer. In particular, we focus on the effect of the UV screening on the dynamics and on the statistical properties of molecular clouds. Methods. We have developed a tree-based method for a fast estimate of column densities, implemented in the adaptive mesh refinement code RAMSES. We performed numerical simulations using this method in order to analyze the influence of the screening on the clump formation. Results. We find that the accuracy for the extinction of the tree-based method is better than 10%, while the relative error for the column density can be much more. We describe the implementation of a method based on precalculating the geometrical terms that noticeably reduces the calculation time. To study the influence of the screening on the statistical properties of molecular clouds we present the probability distribution function (PDF) of gas and the associated temperature per density bin and the mass spectra for different density thresholds. Conclusions. The tree-based method is fast and accurate enough to be used during numerical simulations since no communication is needed between CPUs when using a fully threaded tree. It is then suitable to parallel computing. We show that the screening for far UV radiation mainly affects the dense gas, thereby favoring low temperatures and affecting the fragmentation.

연구 동기 및 목표

• 분자의 구름에 대한 고해상도, 다차원 시뮬레이션에서 열량 밀도를 계산하는 데 발생하는 계산적 과제를 해결한다.
• 시뮬레이션 중 실시간으로 열량 밀도를 추정할 수 있는 빠르고 병렬 처리 가능한 방법을 개발한다.
• 자외선 복사 스크리닝이 분자의 구름의 역학적 및 열역학적 구조, 특히 밀도가 높은 영역에 미치는 영향을 조사한다.
• 자외선 스크리닝이 덩어리 형성과 기체의 통계적 특성(예: PDF 및 질량 스펙트럼)에 끼치는 영향을 평가한다.

제안 방법

• RAMSES AMR 코드에 계층적 트리 구조를 활용하여 열량 밀도를 효율적으로 계산하는 트리 기반 알고리즘을 구현한다.
• 트리 순회 중에 거리가 증가하는 셀의 기여를 누적하는 수거(gather) 방식을 사용하여 전체 4π 스테라디안의 입체각을 근사한다.
• 기하학적 항목을 사전 계산하여 저장하는 모듈을 도입함으로써 평가 시간을 단축시키며 정확도를 유지한다.
• 셀의 각도 크기와 거리에 가중치를 부여하여 트리 내 셀의 기여를 합산함으로써 열량 밀도를 추정한다. 이는 해상도에 따라 달라지는 기여를 포함한다.
• 추정된 열량 밀도를 사용하여 흡수율을 계산하고 먼지가 원거리 자외선 복사에 미치는 스크리닝 효과를 모델링한다.
• 이 방법을 시뮬레이션에 적용하여 자외선 스크리닝 유무에 따라 기체 온도, 밀도 PDF, 덩어리 질량 함수 등에서의 차이를 분석한다.

실험 결과

연구 질문

• RQ1 고해상도 AMR 시뮬레이션에서 정확도를 크게 떨어뜨리지 않고 열량 밀도 추정을 어떻게 가속화할 수 있는가?
• RQ2 자외선 복사 스크리닝은 분자의 구름의 열역학적 및 역학적 구조에 어떤 영향을 미치는가?
• RQ3 자외선 스크리닝을 포함할 경우 밀도가 높은 덩어리 형성과 그로 인한 구조의 질량 스펙트럼에 어떤 영향을 미치는가?
• RQ4 turbulent, 자기중력 작용이 있는 기체에서 자외선 스크리닝은 잼 질량과 분열 효율성에 얼마나 큰 영향을 미치는가?

주요 결과

• 트리 기반 방법은 상대 오차가 최대 50%까지 허용되지만, 흡수율 지도의 오차는 10% 이하로 유지된다.
• 이 방법은 병렬 처리에 매우 효율적이며, 트리가 완전히 스레드 처리된 경우 상호 프로세서 통신이 필요하지 않다.
• 자외선 스크리닝은 냉각 중성 매질(CNM)의 온도를 최대 50%까지 감소시키지만, 따뜻한 기체(WNM)에는 거의 영향을 주지 않는다.
• 자외선 스크리닝에 의한 냉각으로 인해 밀도가 높은 영역에서 잼 질량이 감소하여 더 작은 스케일의 구조가 중력 응축을 유도한다.
• 자외선 스크리닝을 포함한 시뮬레이션은 특히 밀도 임계치 ≥2500 cm⁻³에서 더 컴팩트하고 밀도가 높은 덩어리를 생성한다. 질량 스펙트럼에서 이를 확인할 수 있다.
• 자외선 스크리닝이 포함된 경우, 기체의 확률 밀도 함수(PDF)와 각 밀도 구간의 온도 분포에서 밀도가 높은 영역에서의 편차가 더 두드러진다.