[论文解读] 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)|Oct 2, 2014

Astrophysics and Star Formation Studies参考文献 33被引用 6

一句话总结

本文提出一种快速、基于树结构的算法，用于在自适应矩量网格（AMR）代码中估算柱密度，实现分子云模拟过程中高效的实时辐射转移计算。该方法将计算成本降低至O(N log N)，且处理器间通信极少，使消光误差低于10%，柱密度估算误差低于50%。研究发现，紫外辐射屏蔽会优先冷却致密气体，降低 Jeans 质量并增强碎片化，导致包含屏蔽效应的模拟中形成更紧凑、高密度的团块。

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代码中实现基于树结构的算法，利用其分层树结构高效计算柱密度。
• 在树遍历过程中采用汇聚方法，累积来自距离逐渐增加的单元的贡献，近似完整4π立体角。
• 引入预计算模块，存储几何项以加快评估速度，从而在不损失精度的前提下减少计算时间。
• 通过加权求和树中单元的贡献来估算柱密度，权重包括角大小和距离，且贡献随分辨率变化。
• 利用估算的柱密度计算消光，并建模尘埃对远紫外辐射的屏蔽效应。
• 在模拟中应用该方法，对比包含与不包含紫外屏蔽的模拟结果，分析气体温度、密度PDF及团块质量函数的差异。

实验结果

研究问题

• RQ1 如何在高分辨率AMR模拟中加速柱密度估算，同时保持较高精度？
• RQ2 紫外辐射屏蔽对分子云热结构与动力学结构有何影响？
• RQ3 包含紫外屏蔽后，对致密团块的形成及结构质量谱有何影响？
• RQ4 紫外屏蔽在湍流自引力气体中在多大程度上改变Jeans质量与碎片化效率？

主要发现

• 基于树的算法在柱密度估算中相对误差最高达50%，但消光图的误差低于10%。
• 该方法在并行计算中极为高效，当树结构完全线程化时，无需处理器间通信。
• 紫外屏蔽显著降低了冷中性介质（CNM）的温度，降幅达50%，而暖气体（WNM）温度基本不受影响。
• 由于紫外屏蔽导致的冷却作用，致密区域的Jeans质量降低，更有利于小尺度结构的引力坍缩。
• 包含紫外屏蔽的模拟产生更多紧凑、高密度的团块，尤其在密度阈值≥2500 cm⁻³时，质量谱中表现明显。
• 当包含紫外屏蔽时，气体及每密度区间温度的概率密度函数（PDF）在致密区域表现出更强的偏离。