[논문 리뷰] On the evolution of shape in N-body simulations

Ch. Theis, Rainer Spurzem|arXiv (Cornell University)|1998. 10. 01.

Stellar, planetary, and galactic studies인용 수 4

한 줄 요약

이 논문은 특수한 GRAPE 및 HARP 컴퓨터를 사용하여 냉각되고 에너지 손실이 없는 N-body 시스템의 형상 진화를 조사하며, 세 가지 명확한 단계를 밝혀내었다: 빠른 붕괴로 인한 반경 방향 궤도 불안정성과 최대 삼축성, 그 후 폭발적이고 이중 입자 상호작용에 의한 완화로 축대칭으로 향하고 마지막으로 구형 대칭으로 수렴하는 과정이다. 주요 발견은 부드러움 조정을 정확히 반영할 경우, 천천히 형상 진화하는 시간 스케일이 이중 입자 상호작용 시간 스케일과 매우 밀접하게 일치하며, 핵 진동이 폭발적 진화의 원인임을 규명하였다.

ABSTRACT

A database on shape evolution of direct N-body models formed out of cold, dissipationless collapse is generated using GRAPE and HARP special purpose computers. Such models are important to understand the formation of elliptical galaxies. Three dynamically distinct phases of shape evolution were found, first a fast dynamical collapse which gives rise to the radial orbit instability (ROI) and generates at its end the maximal triaxiality of the system. Subsequently, two phases of violent and two-body shape relaxation occur, which drive the system first towards axisymmetry, finally to spherical symmetry (the final state, however, is still much more concentrated than the initial model). In a sequence of models the influence of numerical and physical parameters, like particle number, softening, initial virial ratio, timestep choice, different N-body codes, are examined. We find that an improper combination of softening and particle number can produce erroneous results. Selected models were evolved on the secular timescale until they became spherically symmetric again. The secular shape relaxation time scale is shown to agree very well with the two-body relaxation time, if softening is properly taken into account for the latter. Finally, we argue, that the intermediate phase of violent shape relaxation after collapse is induced by strong core oscillations in the centre, which cause potential fluctuations, dampening out the triaxiality.

연구 동기 및 목표

• 냉소성 은하 형성과 관련된 냉각되고 에너지 손실이 없는 N-body 시스템의 형상 진화 역학적 진화를 이해하기 위해.
• 삼축성에서 구형 대칭으로의 형상 진화를 이끄는 물리적 메커니즘을 규명하기 위해.
• 입자 수, 부드러움, 타임스텝과 같은 수치적 매개변수들이 형상 진화 시뮬레이션의 정확성과 신뢰성에 미치는 영향을 평가하기 위해.
• 적절한 부드러움 조건 하에서 천천히 형상 진화하는 시간 스케일과 이중 입자 상호작용 시간 스케일 간의 연결 고리를 검증하기 위해.

제안 방법

• 고정밀도로 냉각되고 에너지 손실이 없는 붕괴를 모델링하기 위해 GRAPE 및 HARP 전용 하드웨어를 사용하여 시뮬레이션을 수행하였다.
• 다양한 입자 수, 부드러움 길이, 초기 비르갈 비율, 타임스텝 선택지를 가진 N-body 모델 데이터베이스를 구축하였다.
• 다양한 구현 방식 간의 일관성과 수치적 강건성을 검증하기 위해 여러 N-body 코드를 사용하였다.
• 시간에 따른 주축 비율과 삼축성 척도를 통해 형상 진화를 추적하였다.
• 관측된 천천히 형상 진화하는 시간 스케일과 비교하기 위해, 적절한 부드러움 보정을 포함한 이중 입자 상호작용 시간을 계산하였다.
• 폭발적 형상 진화의 역할을 규명하기 위해 핵 진동과 포텐셜 변동을 분석하였다.

실험 결과

연구 질문

• RQ1 냉각되고 에너지 손실이 없는 N-body 시스템의 형상 진화에서 명백한 역학적 단계는 무엇인가?
• RQ2 붕괴 후 폭발적 형상 진화에 기여하는 핵 진동과 포텐셜 변동은 어떻게 작용하는가?
• RQ3 부드러움과 입자 수와 같은 수치적 매개변수가 형상 진화 시뮬레이션의 정확성에 얼마나 큰 영향을 미치는가?
• RQ4 부드러움을 정확히 반영할 경우, 천천히 형상 진화하는 시간 스케일이 이중 입자 상호작용 시간 스케일과 일치하는가?
• RQ5 이 시스템에서 삼축성에서 구형 대칭으로의 전이를 일으키는 원인는 무엇인가?

주요 결과

• 세 가지 명백한 형상 진화 단계가 확인되었다: 반경 방향 궤도 불안정성과 최대 삼축성을 유도하는 빠른 천체적 붕괴, 그 후 폭발적 및 이중 입자 상호작용에 의한 완화 단계.
• 시스템은 삼축성에서 축대칭으로 진화하여 궁극적으로는 구형 대칭으로 수렴하지만, 최종 상태는 초기 모델보다 훨씬 더 농축되어 있다.
• 부드러움과 입자 수의 부적절한 조합은 N-body 시뮬레이션에서 오류 결과를 초래할 수 있다.
• 부드러움을 정확히 계산에 포함할 경우, 천천히 형상 진화하는 시간 스케일이 이중 입자 상호작용 시간 스케일과 매우 잘 일치한다.
• 폭발적 형상 진화의 중간 단계는 삼축성을 감소시키는 잠재적 변동을 유도하는 강력한 핵 진동에 의해 이끌린다.
• 핵 진동은 붕괴 이후 폭발적 진화를 유도하는 주요 물리적 메커니즘으로 규명되었다.