[논문 리뷰] N-MODY: a code for collisionless N-body simulations in modified Newtonian dynamics

P. Londrillo, Carlo Nipoti|ArXiv.org|2008. 03. 31.

Scientific Research and Discoveries인용 수 3

한 줄 요약

N-MODY는 수정 뉴턴 역학(MOND) 하에서 충돌 없는 시스템을 시뮬레이션하기 위한 구면좌표계 기반의 병렬 입자-메쉬 N체 코드로, 구면 조화함수를 사용한 반복 방법을 통해 비선형 MOND 필드 방정식을 해결한다. 이는 뉴턴 중력과 다크 물질이 있는 경우와 비교할 때 폭력적 완화, 융합, 천체역학적 마찰과 같은 천체역학적 과정이 상당히 다름을 드러내며, 완화와 융합의 시간스케일이 더 길고 마찰력이 더 강하다.

ABSTRACT

We describe the numerical code N-MODY, a parallel particle-mesh code for collisionless N-body simulations in modified Newtonian dynamics (MOND). N-MODY is based on a numerical potential solver in spherical coordinates that solves the non-linear MOND field equation, and is ideally suited to simulate isolated stellar systems. N-MODY can be used also to compute the MOND potential of arbitrary static density distributions. A few applications of N-MODY indicate that some astrophysically relevant dynamical processes are profoundly different in MOND and in Newtonian gravity with dark matter.

연구 동기 및 목표

• 수정 뉴턴 역학(MOND)에서 은하 운동을 설명하는 데 다크 물질 대신 사용할 수 있는 충돌 없는 N체 시뮬레이션을 위한 수치 코드를 개발하기 위해.
• 구면좌표계에서 구면 조화함수 전개를 사용한 반복 방법을 통해 비선형 MOND 필드 방정식을 효율적으로 해결하기 위해.
• 고립된 별계 및 임의의 정적 밀도 분포를 가진 시스템에 대해 MOND 하에서의 시뮬레이션을 가능하게 하여 MOND를 뉴턴 역학과 대비해 시험할 수 있는 도구를 제공하기 위해.
• 폭력적 완화, 위상 혼합, 융합, 천체역학적 마찰과 같은 기본 천체역학적 과정이 다크 물질이 있는 뉴턴 중력과 비교해 MOND에서 어떻게 다르게 나타나는지 조사하기 위해.
• 공개 가능한 병렬 FORTRAN 90 코드를 제공하여 고성능 컴퓨팅 환경에서 순차적 및 MPI 병렬 실행을 모두 지원하기 위해.

제안 방법

• N-MODY는 입자-메쉬 기반의 구면좌표계를 사용하며, 격자 위에서 선형 또는 2차형태 함수를 통해 밀도장을 재구성한다.
• MOND 포텐셜은 비선형 필드 방정식 ∇·[μ(|∇ϕ|/a₀)∇ϕ] = 4πGρ를 반복 방법과 구면 조화함수 전개를 사용해 격자에서 해결함으로써 계산된다.
• 수치적 해를 얻기 위해 반복 계수(ω⁻¹ ≈ 0.3–0.5)와 수렴 정밀도(ε = 10⁻³)를 사용하여 ∥M̂(ϕ)∥∞ ≈ 0.1 및 ∥M̂(ϕ)∥rms ≈ 0.01의 정확도를 달성한다.
• 입자 가속도는 동일한 형태 함수를 사용해 격자에서 보간한 후, 구면좌표에서 카르테시안 성분으로 변환하여 시간 적분을 수행한다.
• 시간 적분은 2차 또는 4차 수준의 리프레그 방법을 사용하여 시뮬레이션에서 에너지 및 운동량 보존을 보장한다.
• 코드는 MPI를 통해 병렬화되며, 입자들이 프로세서 간에 균일하게 분할되고, θ좌표에 대해 도메인 분할을 적용하여 포텐셜 해소기의 부분적 병렬화를 실현한다.

실험 결과

연구 질문

• RQ1 MOND에서 폭력적 완화와 위상 혼합의 시간스케일은 다크 물질이 있는 뉴턴 중력과 비교해 어떻게 다를까?
• RQ2 MOND 시스템에서 천체역학적 마찰의 효율성은 뉴턴 중력 시스템과 비교해 어떻게 다를까?
• RQ3 은하의 융합 시간스케일은 MOND와 뉴턴 중력 간에 어떻게 다를까?
• RQ4 비선형 MOND 필드 방정식은 뉴턴 포텐셜과 비교해 고립된 별계의 구조와 진화에 얼마나 큰 영향을 미치는가?
• RQ5 구면좌표 기반의 입자-메쉬 코드는 천체물리학적으로 관련 있는 시스템에 대해 MOND 필드 방정식을 효율적이고 정확하게 해결할 수 있는가?

주요 결과

• N-MODY의 MOND 포텐셜 해소기는 ε = 10⁻³ 정밀도로 5~10회 반복 내에 수렴하며, ∥M̂(ϕ)∥∞ ≈ 0.1 및 ∥M̂(ϕ)∥rms ≈ 0.01의 정확도를 달성한다.
• 시리얼 모드에서 N ≈ 10Ng일 때 MOND 포텐셜 해소기는 총 CPU 시간의 약 10%를 차지하지만, 병렬 실행에서는 병렬화가 제한되어 있어 주요 성능 저하 요소가 된다.
• MOND 해소기의 병렬화 효율은 약 70%이며, 계산 비용은 격자 포인트 수 Ng = (Nr+1)NθNφ에 비례하여 증가한다.
• 16개 프로세서를 사용할 때도 N = 10⁷개 입자 및 Ng = 128×64×128일 경우, N-MODY의 입자당 타임스텝당 총 CPU 시간은 10⁻⁵초 이하로 유지된다.
• MOND 시뮬레이션에서는 폭력적 완화, 위상 혼합, 은하 융합이 뉴턴 중력보다 상당히 느리게 일어나며, 이는 천체역학적 과정이 억제됨을 시사한다.
• MOND 시스템에서는 다크 물질이 있는 동일한 뉴턴 시스템보다 천체역학적 마찰력이 더 강해, 질량이 큰 천체의 궤도 감쇠 속도가 증가함을 시사한다.