[논문 리뷰] Particle-Mesh code for cosmological simulations

Anatoly Klypin, Jon A. Holtzman|arXiv (Cornell University)|1997. 12. 17.

Scientific Research and Discoveries인용 수 55

한 줄 요약

이 논문은 천체역학 시뮬레이션을 위한 공개 가능하고 빠르며 단순한 입자-메쉬(Particle-Mesh, PM) N-body 코드를 제시한다. 이는 평탄한, 열린, 닫힌 우주를 포함한 다양한 천체역학 모델의 효율적 탐색을 가능하게 하며, 우주상수의 유무, 변하는 허블 상수, 뜨거운 중성자, 기울어진 초기 스펙트럼 등을 수용한다. 이 코드는 표준 워크스테이션에서 최대 $800^3$ 메쉬와 $256^3$ 입자를 지원하며, 파wer 스펙트럼 분석, 유사결집 밀도 최대값 기반의 은하단 탐지, 포괄적인 데이터 분석 루틴을 내장하고 있다.

ABSTRACT

Particle-Mesh (PM) codes are still very useful tools for testing predictions of cosmological models in cases when extra high resolution is not very important. We release for public use a cosmological PM N-body code. We provide a complete package of routines needed to set initial conditions, to run the code, and to analyze the results. The package allows you to simulate models with numerous combinations of parameters: open/flat/closed background, with or without the cosmological constant, different values of the Hubble constant, with or without hot neutrinos, tilted or non-tilted initial spectra, different amount of baryons. Routines are included to measure the power spectrum and the density distribution function in your simulations, and a bound-density-maxima code for halo finding. We also provide results of test runs. A simulation with 256^3 mesh and 128^3 particles can be done in a couple of days on a typical workstation (70Mb of RAM are needed). To run simulations with 800^3 mesh and 256^3 particles one needs a computer with 1Gb memory and 1Gb disk space. The code has been successfully tested on an HP workstation and on a Sun workstation running Solaris. Most of the files (not tests) can be obtained from ftp://astro.nmsu.edu/pub/aklypin/PMCODE The package can be downloaded from http://astro.nmsu.edu/~aklypin/PM/pmcode/index.html We provide this tool as a service to the astronomical community, but we cannot guarantee results or publications.

연구 동기 및 목표

• 천체역학 모델 테스트를 지원하기 위해 빠르고 접근 가능하며 확장 가능한 입자-메쉬(Particle-Mesh, PM) N-body 코드를 제공하는 것.
• 다양한 천체역학적 파라미터, 즉 변하는 허블 상수, 우주상수, 곡률, 복사 성분, 질량이 있는 중성자를 포함한 시뮬레이션을 가능하게 하는 것.
• 초기 조건 설정, 시뮬레이션 실행, 결과 분석을 위한 도구를 포함한 완전한 소프트웨어 패키지를 제공하는 것—파워 스펙트럼 및 은하단 카탈로그 생성 포함.
• 표준 워크스테이션에서 높은 입자 수(최대 $256^3$)와 큰 메쉬($800^3$)를 지원하여 대규모 구조 연구에 적합하게 하는 것.
• 천문학 공동체를 위해 안정적이고 테스트된, 이식 가능한 코드베이스를 제공하며, 문서화된 예제와 광범위한 분석 도구를 제공하는 것.

제안 방법

• 운동 방정식을 단순화하고 수치적 안정성을 향상시키기 위해 공변 좌표계와 무차원 변수를 사용한다.
• 고정 메쉬 위에서 힘 계산을 위해 7점 이산 라플라스 연산자와 클라우드-인-셀(Cloud-in-Cell, CIC) 밀도 할당을 사용한다.
• 무차원 밀도 대비 및 천체역학적 파rameter를 사용하여 중력 포텐셜에 대한 포아송 방정식을 해결한다.
• 확장 인자 $a$ 기반의 시간 스텝 스킴을 사용하여 입자 궤적을 통합하며, 운동 방정식을 입자 운동량 $\mathbf{p}$ 기반으로 재구성한다.
• 임의의 초기 스펙트럼을 가진 초기 조건을 생성하기 위한 루틴을 제공하며, 피팅 함수(예: cdm.fit) 또는 직접 코드 수정을 통해 가능하다.
• 구성 가능한 파rameter(구의 수, 셸 해상도, 반경 빈팅)를 가진 유사결집 밀도 최대값 알고리즘을 통해 은하단 탐지를 지원하며, 탐지 최적화를 위해 설계되어 있다.

실험 결과

연구 질문

• RQ1 다양한 파라미터를 가진 다양한 천체역학 모델을 지원할 수 있도록 빠르고 단순한 PM N-body 코드를 설계하는 방법은 무엇인가?
• RQ2 유사결집 밀도 최대값 기반의 은하단 탐지 알고리즘에 최적의 설정 파라미터는 무엇이며, 이는 대량 및 저질량 은하단 모두를 효과적으로 탐지할 수 있도록 보장하는가?
• RQ3 PM 방법이 $256^3$ 입자와 $800^3$ 메쉬를 가진 천체역학 시뮬레이션에서 파워 스펙트럼과 밀도 분포 함수를 얼마나 정확하게 재현할 수 있는가?
• RQ4 표준 워크스테이션에서 $800^3$ 메쉬와 $256^3$ 입자를 가진 대규모 구조 시뮬레이션을 수행하는 데 필요한 계산 비용과 메모리 요구량은 얼마인가?
• RQ5 질량이 있는 중성자를 포함한 비표준 모델인 CHDM의 경우 초기 조건는 어떻게 생성할 수 있으며, 그에 적용되는 적합한 적색편이 및 스펙트럼 형태 제약 조건은 무엇인가?

주요 결과

• 표준 워크스테이션에서 $256^3$ 메쉬와 $128^3$ 입자를 가진 시뮬레이션은 약 70MB 메모리로 며칠 내로 완료될 수 있다.
• $800^3$ 메쉬와 $256^3$ 입자를 가진 시뮬레이션은 약 1GB 메모리와 1GB 디스크 공간이 필요하며, 현대 워크스테이션에서 실행 가능하다.
• 이 코드는 평탄한, 열린, 닫힌 우주를 포함한 다양한 천체역학 모델을 성공적으로 지원하며, 우주상수의 유무, 변하는 허블 상수, 질량이 있는 중성자를 포함한다.
• 적절한 반경 빈팅과 최소 반경 제약 조건을 적용할 경우, 100,000~150,000개의 시드 구와 1개의 반경 셸에 5~6개의 입자를 사용하여 유사결집 밀도 최대값 알고리즘이 은하단을 효과적으로 탐지한다.
• PMpower 루틴을 통해 파워 스펙트럼과 밀도 분포 함수를 신뢰성 있게 측정할 수 있으며, 출력은 Spectrum.DAT 파일에 저장된다.
• 검증을 위해 $32^3$ 입자와 $128^3$ 메쉬, $128^3$ 입자와 $256^3$ 메쉬의 테스트 시뮬레이션을 제공하며, 전체 데이터와 분석 출력을 포함한다.