[论文解读] sCOLA: The N-body COLA Method Extended to the Spatial Domain

Svetlin Tassev, Daniel J. Eisenstein+2 more|arXiv (Cornell University)|Feb 26, 2015

Galaxies: Formation, Evolution, Phenomena参考文献 27被引用 31

一句话总结

sCOLA 通过将近场 N 体动力学与通过拉格朗日微扰理论（LPT）建模的远场引力影响解耦，将 N 体 COLA 方法扩展到空间域，实现了快速、精确的 zoom-in 模拟以及极易并行化的宇宙学模拟。它使小体积 N 体模拟能够重现大规模模拟的结果，通过微扰方法捕捉长程模式，显著降低计算成本，同时保持大尺度结构研究的高精度。

ABSTRACT

We present sCOLA -- an extension of the N-body COmoving Lagrangian Acceleration (COLA) method to the spatial domain. Similar to the original temporal-domain COLA, sCOLA is an N-body method for solving for large-scale structure in a frame that is comoving with observers following trajectories calculated in Lagrangian Perturbation Theory. Incorporating the sCOLA method in an N-body code allows one to gain computational speed by capturing the gravitational potential from the far field using perturbative techniques, while letting the N-body code solve only for the near field. The far and near fields are completely decoupled, effectively localizing gravity for the N-body side of the code. Thus, running an N-body code for a small simulation volume using sCOLA can reproduce the results of a standard N-body run for the same small volume embedded inside a much larger simulation. We demonstrate that sCOLA can be safely combined with the original temporal-domain COLA. sCOLA can be used as a method for performing zoom-in simulations. It also allows N-body codes to be made embarrassingly parallel, thus allowing for efficiently tiling a volume of interest using grid computing. Moreover, sCOLA can be useful for cheaply generating large ensembles of accurate mock halo catalogs required to study galaxy clustering. Surveys that will benefit the most are ones with large aspect ratios, such as pencil-beam surveys, where sCOLA can easily capture the effects of large-scale transverse modes without the need to substantially increase the simulated volume. As an illustration of the method, we present proof-of-concept zoom-in simulations using a freely available sCOLA-based N-body code.

研究动机与目标

• 开发一种在宇宙学模拟中将近场 N 体动力学与远场引力势分离的方法。
• 通过微扰理论建模长程引力效应，实现在小体积区域内的精确、高速模拟。
• 通过将大体积划分为独立的小模拟，实现高效、极易并行的模拟。
• 支持生成大量精确的模拟晕系度量表，用于星系聚类与协方差矩阵估计。
• 通过在不增加模拟体积的情况下捕捉横向大尺度模式，提高对铅笔束与高长宽比调查的模拟效率。

提出的方法

• sCOLA 通过将拉格朗日微扰理论（LPT）应用于空间模式，将时间域的 COLA 方法扩展，将大尺度引力势视为微扰修正。
• 该方法将模拟分为局部 N 体盒子（近场）与全体积 LPT 计算（远场），远场通过解析势能修正影响近场。
• 使用两点或四点有限差分规则计算一阶和二阶 LPT 位移场，以实现引力势的高精度近似。
• N 体代码中的力计算使用 PM 网格与云-在-细胞赋值方法，在由 LPT 场偏移的欧拉位置上评估势能。
• 该方法采用位移场方程的修改形式，引入参数 α 和 β，以确保与 ΛCDM 和物质主导极限的一致性。
• sCOLA 框架已集成到公开的 pyCOLA 代码中，允许用户在标准 COLA 与 sCOLA 模式之间切换，实现混合模拟。

实验结果

研究问题

• RQ1 是否可以通过微扰理论在不运行完整大体积模拟的情况下，精确建模远处物质对局部模拟体积的引力影响？
• RQ2 sCOLA 在小体积区域内通过解析方式捕捉远场效应，能在多大程度上重现完整 N 体模拟的结果？
• RQ3 sCOLA 与原始时间域 COLA 的结合如何提升宇宙学模拟中的计算效率与精度？
• RQ4 sCOLA 是否能够实现高效、可扩展且极易并行的模拟，适用于高长宽比的大尺度结构调查？
• RQ5 sCOLA 在生成用于星系聚类研究的模拟晕系度量表方面，其性能与精度如何？

主要发现

• sCOLA 通过 LPT 建模远场引力效应，在小体积 N 体模拟中成功重现了完整模拟的大尺度结构。
• 该方法可在不扩展模拟体积的情况下实现精确的 zoom-in 模拟，显著降低计算成本。
• sCOLA 通过将感兴趣区域划分为独立的局部 N 体模拟并由远场 LPT 修正，实现了极易并行的模拟。
• sCOLA 与原始时间域 COLA 的结合，为宇宙学模拟提供了一个灵活、高精度且计算高效的框架。
• sCOLA 在铅笔束与高长宽比调查中特别有效，可在不增加模拟体积的情况下捕捉横向大尺度模式。
• 基于公开的 sCOLA 版 pyCOLA 代码的原理验证模拟表明，该方法在保持计算速度提升的同时实现了高精度。