[논문 리뷰] A family of models of partially relaxed stellar systems II. Comparison with the products of collisionless collapse

Michele Trenti, G. Bertin외 1명|University of Groningen research database (University of Groningen / Centre for Information Technology)|2004. 11. 18.

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

한 줄 요약

이 논문은 비충격적 붕괴의 N-body 시뮬레이션과 부분적으로 완화된 별성계를 묘사하기 위한 목적의 분석 모델 가족(f^(ν)) 간의 상세한 비교를 제시한다. f^(ν) 모델이 밀도, 이방성, 위상공간 구조를 매우 정확하게 재현함을 보여주며, 특히 위상공간에 덩어리가 존재하는 초기 조건일 경우에 그 정확도가 높아진다. 이는 모델이 붕괴의 최종 산물과 높은 정밀도로 일치함을 확인하며, 특히 경계 근처의 축방향 궤도 불안정성 조건에서도 성립함을 시사한다.

ABSTRACT

[Abridged] N-body simulations of collisionless collapse have offered important clues to the construction of realistic stellar dynamical models of elliptical galaxies. Understanding this idealized and relatively simple process, by which stellar systems can reach partially relaxed equilibrium configurations (characterized by isotropic central regions and radially anisotropic envelopes), is a prerequisite to more ambitious attempts at constructing physically justified models of elliptical galaxies in which the problem of galaxy formation is set in the generally accepted cosmological context of hierarchical clustering. In a previous paper, we have discussed the dynamical properties of a family of models of partially relaxed stellar systems (the f^{(ν)} models), designed to incorporate the qualitative properties of the products of collisionless collapse at small and at large radii. Here we revisit the problem of incomplete violent relaxation, by making a direct comparison between the detailed properties of such family of models and those of the products of collisionless collapse found in N-body simulations that we have run for the purpose. Surprisingly, the models thus identified are able to match the simulated density distributions over nine orders of magnitude and also to provide an excellent fit to the anisotropy profiles and a good representation of the overall structure in phase space.

연구 동기 및 목표

• 비충격적 붕괴를 통해 형성된 시스템의 구조적 및 운동학적 성질을 f^(ν) 모델 가족이 정량적으로 재현할 수 있는지 테스트하기 위해.
• 비충격적 붕괴에서 단일 입자 각운동량 혼합을 효율적으로 가능하게 하는 초기 위상공간 덩어리의 역할을 조사하기 위해.
• 붕괴 과정 중 양자수 Q의 수치적 보존성과 그 모델의 예측 능력에 대한 관련성을 평가하기 위해.
• 붕괴의 최종 산물이 축방향 궤도 불안정성의 임계점 근처의 이방성 수준을 가지는지 확인하기 위해.
• 우주론적 맥락에서 타원은하를 모델링하는 데 사용할 수 있도록 이상화된 시뮬레이션과 f^(ν) 모델 간의 정량적 다리를 구축하기 위해.

제안 방법

• 균일한(U-type) 및 덩어리가 있는(C1–C4.5) 초기 조건를 포함한 다양한 초기 조건으로 비충격적 붕괴의 광범위한 N-body 시뮬레이션을 수행하기 위해.
• 6차원 위상공간 덩어리 정도를 측정하기 위해 cl = ⟨ρ_local^(6)⟩ / ⟨ρ^(6)⟩ 를 사용하며, 국소 구의 반지름 rs는 전체 입자 수의 약 1/250를 포함하도록 고정한다.
• 밀도 프로파일, 이방성 프로파일, 전체 위상공간 구조를 비교하여 최종 시뮬레이션 시스템과 f^(ν) 모델 간의 일치를 평가하기 위해.
• 특히 다양한 초기 비르갈 비율 u에 대해 양자수 Q의 시간적 변화를 추적하여 붕괴 과정에서의 근사 보존성을 테스트하기 위해.
• 위치 및 속도 공간에서의 초기 덩어리 정도가 혼합 효율성과 최종 평형 구조에 미치는 영향을 분석하기 위해.
• 수치적 통합을 위해 트리 코드를 사용하고, 기존 기준(예: Dehnen 2000)과의 비교를 통해 결과를 검증하기 위해.

실험 결과

연구 질문

• RQ1 f^(ν) 모델은 N-body 시뮬레이션에서 비충격적 붕괴를 통해 형성된 시스템의 밀도 및 이방성 프로파일을 정량적으로 일치시킬 수 있는가?
• RQ2 초기 위상공간 덩어리가 비충격적 붕괴에서 효과적인 각운동량 혼합과 평형에 도달하는 데 어떤 역할을 하는가?
• RQ3 에너지와 질량 이외의 제3의 보존량으로서 도입된 Q가 붕괴 과정 중에 근사적으로 보존되는가?
• RQ4 붕괴의 최종 산물은 축방향 궤도 불안정성의 임계점 근처의 이방성 수준을 가지는가?
• RQ5 초기 조건—특히 덩어리의 수와 크기—는 시뮬레이션 시스템의 최종 구조적 및 운동학적 성질에 어떤 영향을 미치는가?

주요 결과

• f^(ν) 모델은 반지름이 9개 온스 오더에 걸쳐 시뮬레이션된 밀도 프로파일을 매우 정확하게 맞춘다.
• 모델은 또한 시뮬레이션된 이방성 프로파일과 높은 정확도로 일치하며, 이는 f^(ν) 가족이 외곽에서는 축방향 이방성, 핵심에서는 등방성을 잘 묘사하고 있음을 시사한다.
• 시뮬레이션의 최종 산물과 최적의 f^(ν) 모델은 모두 축방향 궤도 불안정성의 임계점 근처에 위치한 전반적인 이방성 파라미터를 공유한다.
• 특히 속도 공간에서의 위상공간 덩어리가 단일 입자 각운동량 혼합을 효율적으로 가능하게 하여 f^(ν) 모델과 일치하는 상태에 도달하는 데 핵심적인 역할을 한다.
• 덩어리 정도 측정치 cl 은 덩어리가 있는 초기 조건(예: C4.4는 10개의 작은 덩어리)에서 ξ = 1/250 에서 cl ≈ 300 을 보이고, 균일한 조건에서는 cl ≈ 0.72 를 보이며, 이는 비균일성의 중요성을 확인한다.
• 붕괴 정도가 높은 경우(u ≈ 0.1–0.3)에 Q 가 시간에 따라 근사적으로 보존되며, 이는 f^(ν) 모델 구축에서 Q 가 핵심 요소임을 뒷받침한다.