[논문 리뷰] Time stepping N-body simulations

Thomas Quinn, Neal Katz외 2명|arXiv (Cornell University)|1997. 10. 03.

Numerical methods for differential equations참고 문헌 2인용 수 14

한 줄 요약

이 논문은 개별 입자마다 별도의 시간단계를 허용하면서도 단일 입자 시험에서는 심플렉틱성을 유지하는 N-체 시뮬레이션을 위한 시간에 따라 적응하는 루프프로 정적분기법을 소개한다. 전체 시뮬레이션에서는 시간단계 선택으로 인해 심플렉틱성이 손실되지만, 저자들은 심플렉틱성을 기각하지만 높은 정확도를 유지하는 수정된 정적분기법을 개발하여 대규모 천체역학 시뮬레이션에서 표준 방법을 능가한다. 또한, 천체역학적 N-체 코드에서 오랫동안 존재해 온 결함을 보완하기 위해 공변 좌표계에서 심플렉틱성을 갖는 루프프로 형식을 유도한다.

ABSTRACT

Leapfrog integration has been the method of choice in N-body simulations owing to its low computational cost for a symplectic integrator with second order accuracy. We introduce a new leapfrog integrator that allows for variable timesteps for each particle in large N-body simulations. Tests with single particles in fixed potentials show that it behaves as a symplectic integrator. We then examine the results of both standard leapfrog and our temporally adaptive leapfrog on full N-body integrations of clusters and large scale structure establishing accuracy criteria for both methods. The adaptive method shows significant speed-ups over single step integrations---but the integrator no longer appears to be symplectic or, in the case of large scale structure simulations, accurate. This loss of accuracy appears to be caused by the way that the timestep is chosen, not by the integrator itself. We present a related integration technique that does retain sufficient accuracy. Although it is not symplectic, it is apparently better than previous implementations and is our current integrator of choice for large astrophysical simulations. We also note that the standard leapfrog difference equations used in cosmological N-body integrations in comoving coordinates are not symplectic. We derive an implementation of leapfrog that is in comoving canonical coordinates to correct for this deficiency.

연구 동기 및 목표

• 대규모 N-체 시뮬레이션에서 개별 입자 시간단계를 허용하는 시간에 따라 적응하는 심플렉틱 정적분기법을 개발하는 것.
• 적응형 루프프로 방법을 전체 N-체 시스템에 적용했을 때 심플렉틱성과 정확도가 손실되는 이유를 다루는 것.
• 천체역학적 시뮬레이션에서 널리 사용되는 공변 좌표계에서의 표준 루프프로 방법이 비심플렉틱한 점을 수정하는 것.
• 대규모 천체물리학적 시뮬레이션에 실용적이고 정확하며 효율적인 정적분기법을 제공하는 것.

제안 방법

• 단일 입자 시험에서 심플렉틱성을 유지하면서도 각 입자마다 가변 시간단계를 허용하는 새로운 루프프로 정적분기법을 제안한다.
• 입자가 지역적인 천체역학적 시간스케일에 따라 개별 시간단계로 진행되는 계층적 시간단계 기법을 사용한다.
• 라그랑지안에 캐논성 변환를 적용하여 공변 캐논ical 좌표계에서 심플렉틱 정적분기법을 유도함으로써 에너지 및 위상공간 보존을 보장한다.
• 카우프만-암스트라트(KAM) 정리를 활용하여 해밀토니안 시스템에서 심플렉틱 정적분기법의 안정성을 정당화한다.
• 전체 N-체 시뮬레이션에서 심플렉틱성을 기각하지만 높은 정확도와 성능를 유지하는 수정된 정적분기법을 구현한다.
• 비평형 및 임계 상태의 물질 우주에서 시간 적분의 해석적 표현을 유도하여 공변 좌표계에서 정확한 통합을 가능하게 한다.

실험 결과

연구 질문

• RQ1 적응형 루프프로 정적분기법은 대규모 N-체 시뮬레이션에서 심플렉틱성을 유지하면서도 시간에 따라 적응할 수 있는가?
• RQ2 단일 입자에서는 잘 작동하지만 전체 시스템에 적용했을 때 적응형 루프프로 방법이 정확도와 심플렉틱성을 잃는 이유는 무엇인가?
• RQ3 적응형 시간단계 사용으로 인해 정확도가 떨어지는 원인와 성능을 희생시키지 않고 이를 보완할 수 있는가?
• RQ4 공변 좌표계에서 사용되는 표준 루프프로 방법은 진정으로 심플렉틱한가? 만약 그렇지 않다면 어떻게 수정할 수 있는가?
• RQ5 비심플렉틱이지만 매우 정확한 정적분기법을 개발하여 대규모 천체역학 시뮬레이션에서 표준 방법을 능가할 수 있는가?

주요 결과

• 제안된 적응형 루프프로 정적분기법은 고정된 잠재력 내에서 단일 입자 시험에서 심플렉틱성과 2차 정확도를 유지한다.
• 전체 N-체 시뮬레이션에서 적응형 루프프로 방법은 정적분기법의 구조가 아니라 시간단계 선택 알고리즘으로 인해 심플렉틱성과 정확도를 상실한다.
• 심플렉틱성을 기각하지만 개선된 시간단계 기준(η = 0.3)을 사용하는 수정된 정적분기법은 높은 정확도를 달성하며 이제 저자들이 대규모 천체물리학 시뮬레이션에 가장 선호하는 정적분기법이 되었다.
• 공변 좌표계에서 사용되는 표준 루프프로 방법은 비심플렉틱하지만, 캐논성 변환와 적절한 공변 캐논ical 좌표계에서의 시간 적분을 적용한 수정된 버전은 심플렉틱함을 입증하였다.
• 공변 좌표계에서의 새로운 심플렉틱 정적분기법은 임계 및 비평형 우주에서 시간 적분의 해석적 해를 사용하여 정확하고 효율적인 장기 통합을 가능하게 한다.
• 저자들은 시간에 따라 적응하는 기법이 자연 시간단계의 역수에 비례하는 작업을 기반으로 하면 대규모 N-체 시뮬레이션에서 최대 약 50배의 성능 향상을 이룰 수 있음을 입증하였다.