[논문 리뷰] sCOLA: The N-body COLA Method Extended to the Spatial Domain

Svetlin Tassev, Daniel J. Eisenstein외 2명|arXiv (Cornell University)|2015. 02. 26.

Galaxies: Formation, Evolution, Phenomena참고 문헌 27인용 수 31

한 줄 요약

sCOLA는 근접장 N-body 역학과 원거리 중력 영향을 라그랑주 편미분 이론(LPT)을 통해 모델링함으로써 시간 도메인의 N-body COLA 방법을 공간 도메인으로 확장한다. 이는 빠르고 정확한 줌인 시뮬레이션과 막연하게도 병렬화 가능한 천체물리학적 시뮬레이션을 가능하게 한다. 작은 부피의 N-body 시뮬레이션을 통해 매우 큰 시뮬레이션의 결과를 재현할 수 있도록 장거리 모드를 편미분 이론적으로 포괄함으로써 계산 비용을 크게 줄이고, 대규모 구조 연구에서 높은 정확도를 유지한다.

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-body 역학과 원거리 중력 잠재력을 분리하는 방법을 개발하는 것.
• 장거리 중력 효과를 편미분 이론을 사용해 모델링함으로써, 작은 부피 영역에 대해 정확하고 빠른 시뮬레이션을 가능하게 하는 것.
• 큰 부피를 독립적인 작은 시뮬레이션으로 타일링함으로써 효율적이고 막연하게도 병렬화 가능한 시뮬레이션을 지원하는 것.
• 은하 군집 및 공분산 행렬 추정을 위한 정확한 가짜 헬로 카탈로그 대용량 세트 생성을 촉진하는 것.
• 수평 비율이 높거나 펜슬빔 형태의 설문 조사에서 시뮬레이션의 효율성을 향상시키기 위해 횡방향 장거리 모드를 증가시키지 않고도 포착하는 것.

제안 방법

• sCOLA는 시간 도메인의 COLA 방법을 공간 모드에 적용하기 위해 라그랑주 편미분 이론(LPT)을 사용하여, 대규모 중력 잠재력을 편미분 보정으로 간주한다.
• 이 방법은 국소 N-body 상자(근접장)와 전체 부피의 LPT 계산(원거리장)으로 시뮬레이션을 분리하며, 원거리장이 분석적 잠재력 보정을 통해 근접장에 영향을 미친다.
• 1차 및 2차 LPT 이동장은 2점 또는 4점 유한 차분 규칙을 사용하여 중력 잠재력의 고정밀 근사치를 도출한다.
• N-body 코드의 힘 계산은 클러스터-인-셀 할당을 사용하는 PM 격자에 기반하며, LPT 장에 의해 이동된 유레지안 위치에서 잠재력을 평가한다.
• α 및 β 매개변수를 포함한 수정된 이동장 방정식 형태를 사용하여 ΛCDM 및 물질 우주론적 한계와의 일致성을 확보한다.
• sCOLA 프레임워크는 공개된 pyCOLA 코드에 통합되어 사용자가 하이브리드 시뮬레이션을 위해 표준 COLA와 sCOLA 모드를 전환할 수 있도록 한다.

실험 결과

연구 질문

• RQ1 큰 부피의 전체 시뮬레이션을 요구하지 않고도, 편미분 이론을 사용해 원거리 물질이 국소 시뮬레이션 부피에 미치는 중력 영향을 정확하게 모델링할 수 있는가?
• RQ2 원거리 효과를 분석적으로 포착함으로써 sCOLA가 작은 부피 영역에서 전체 N-body 시뮬레이션 결과를 어느 정도 재현할 수 있는가?
• RQ3 원래의 시간 도메인 COLA와 결합된 sCOLA가 천체물리학적 시뮬레이션의 계산 효율성과 정확도를 얼마나 향상시키는가?
• RQ4 sCOLA는 수평 비율이 높은 대규모 구조 설문 조사에 대해 효율적이고 확장 가능하며 막연하게도 병렬화 가능한 시뮬레이션을 가능하게 하는가?
• RQ5 sCOLA는 은하 군집 연구를 위한 가짜 헬로 카탈로그 생성에서 성능과 정확도가 어느 정도인가?

주요 결과

• sCOLA는 LPT를 통해 원거리 중력 영향을 모델링함으로써 작은 부피의 N-body 시뮬레이션 내에서 전체 시뮬레이션의 대규모 구조를 성공적으로 재현한다.
• 이 방법은 시뮬레이션 부피를 늘리지 않고도 정확한 줌인 시뮬레이션을 가능하게 하며, 계산 비용을 크게 줄인다.
• sCOLA는 관심 부피를 독립적인 국소 N-body 실행으로 타일링함으로써 막연하게도 병렬화 가능한 시뮬레이션을 가능하게 하며, 각 실행은 원거리 LPT에 의해 보정된다.
• sCOLA와 원래의 시간 도메인 COLA의 조합은 우주론적 시뮬레이션을 위한 민첩하고 정확하며 계산 효율성이 높은 프레임워크를 제공한다.
• sCOLA는 특히 수평 비율이 높거나 펜슬빔 형태의 설문 조사에서 횡방향 대규모 모드를 증가시키지 않고도 포착하는 데 매우 효과적이다.
• 공개된 sCOLA 기반 pyCOLA 코드를 사용한 개념 증명 시뮬레이션은 이 방법이 높은 정확도를 유지하면서도 계산 속도 향상을 달성함을 보여준다.