[논문 리뷰] A barred Milky Way surrogate from an N-body simulation

Thor Tepper-García, Joss Bland‐Hawthorn외 11명|arXiv (Cornell University)|2021. 11. 10.

Stellar, planetary, and galactic studies참고 문헌 8인용 수 9

한 줄 요약

이 논문은 약 20억 년 이내에 관측된 성질과 일치하는 안정적인 상자/아몬드 모양의 바를 자가일관적으로 형성하는 고해상도 N체 시뮬레이션을 제시한다. 이 성질들은 바의 질량(~10¹⁰ M⊙), 크기(5 kpc), 패턴 속도(40–45 km s⁻¹ kpc⁻¹)를 포함한다. 모델은 관측된 디스크 질량 분율(fₚ ≈ 0.3–0.6)을 반영한 자가일관된 初기 조건을 사용하며, 장기적으로 안정된 진동수 감쇠율(~−2 km s⁻¹ kpc⁻¹ Gyr⁻¹)을 보이며 최근 가이아 데이터와 일치한다.

ABSTRACT

We present an N-body model for the barred Milky Way (MW) galaxy that reproduces many of its properties, including the overall mass distribution, the disc kinematics, and the properties of the central bar. Our high-resolution (N ~ 10^8 particles) simulation, performed with the Ramses code, starts from an axisymmetric non-equilibrium configuration constructed within the AGAMA framework. This is a self-consistent dynamical model of the MW defined by the best available parameters for the dark matter halo, the stellar disc and the bulge. For the known (stellar and gas) disc mass (4.5 x 10^10 Msun) and disc mass fraction at R ~ 2.2 R_d (f_d ~ 0.3 - 0.6), the low mass limit does not yield a bar in a Hubble time. The high mass limit adopted here produces a box/peanut bar within about 2 Gyr with the correct mass (~10^10 Msun), size (~5 kpc) and peak pattern speed (~ 40-45 km/s/kpc). In agreement with earlier work, the bar formation timescale scales inversely with f_d (i.e. log [T/Gyr] ~ 0.60/f_d - 0.83 for 1 < f_d < 0.3). The disc radial heating is strong, but, in contrast to earlier claims, we find that disc vertical heating outside of the box/peanut bulge structure is negligible. The synthetic barred MW exhibits long-term stability, except for the slow decline (roughly -2 km/s/kpc/Gyr) of the bar pattern speed, consistent with recent estimates. If our model is indicative of the Milky Way, we estimate that the bar first emerged 3-4 Gyr ago.

연구 동기 및 목표

• 은하 중심 바와 디스크의 주요 관측 성질을 재현하는 자가일관적이고 고해상도 N체 모델을 구축하는 것.
• 관측된 디스크 질량 분율(fₚ ≈ 0.3–0.6)이 허블 시간 내에 현실적인 바 형성 시스템과 형태를 유도하는지 테스트하는 것.
• 특히 패턴 속도의 진화에 초점을 맞춰, 시뮬레이션된 바가 장기적으로 안정적인지 평가하는 것.
• 바의 중력 잠재력에 의해 유도되는 디스크 가열, 특히 수직 가열의 영향을 평가하는 것.
• 가이아 데이터 해석과 향후 시뮬레이션을 지원하기 위해 공개 가능한 물리적으로 일관된 N체 은하 대체 모델을 제공하는 것.

제안 방법

• 시뮬레이션은 약 10⁸ 개의 입자를 사용하는 ramses N체 코드를 활용하여 자가일관된 은하 유사 시스템을 모델링한다.
• 초기 조건은 관측된 암흑물질 허브, 항성 디스크, 부울지의 파arameters를 반영한 비평형 축대칭 구성으로 agama 프레임워크를 사용해 생성된다.
• 바 형성 시간 내에 성공적으로 형성되도록 하기 위해 fₚ를 관측 범위의 상한(0.3–0.6)으로 설정한다.
• 외부 편향 없이 중력만 작용하는 조건에서 시뮬레이션을 진행하여 내재적 안정성과 천천히 진화하는 성질을 평가한다.
• 표면 밀도, 속도장(Vₜ, V_z), 패턴 속도 진화 등의 합성 관측량을 계산하고 관측 결과와 비교한다.
• 바 성질(질량, 크기, 패턴 속도)과 디스크 운동학을 가이아 및 기타 관측 제약 조건과 비교하여 시뮬레이션을 검증한다.

실험 결과

연구 질문

• RQ1 관측된 디스크 질량 분율(fₚ ≈ 0.3–0.6)을 가진 자가일관된 N체 모델이 은하의 관측된 성질인 질량, 크기, 패턴 속도를 잘 재현하는가?
• RQ2 실제 질량 분율을 반영한 고해상도 자가일관 모델에서 바 형성의 시스템은 얼마나 걸리는가?
• RQ3 상자/아몬드 부울지 영역 외부의 디스크에서 바가 의미 있는 수직 가열을 유도하는가?
• RQ4 특히 패턴 속도의 진화에 초점을 맞춰, Gyr 단위의 시간 스케일에서 시뮬레이션된 바가 은하의 안정성은 어떠한가?
• RQ5 고해상도 N체 시뮬레이션은 은하 바와 내부 디스크 평면의 정렬성과 운동학적 일관성을 재현할 수 있는가?

주요 결과

• 시뮬레이션에서 바가 약 20억 년 이내에 형성되며, 질량은 약 ~10¹⁰ M⊙, 반장력은 약 ~5 kpc, 패턴 속도는 40–45 km s⁻¹ kpc⁻¹로 관측된 은하 바 성질과 일치한다.
• 바 형성 시간은 fₚ에 반비례하며, 0.3 ≤ fₚ ≤ 1 범위에서 log T_bar (Gyr) ≈ 0.60/fₚ − 0.83 의 관계를 따른다.
• 상자/아몬드 부울지 영역 외부의 디스크에서는 유의미한 수직 가열이 관측되지 않아, 이전에 바에 의한 강한 수직 가열을 주장한 주장과 정면으로 배치된다.
• 바 패턴 속도는 약 ~−2 km s⁻¹ kpc⁻¹ Gyr⁻¹ 의 느린 감쇠율을 보이며 최근 가이아 기반 추정치와 일치한다.
• 모델은 장기적으로 천체역학적 안정성을 유지하며, 가스와 디스크는 시간이 지남에 따라 초기 반경 분포를 유지한다.
• 시뮬레이션은 건조하고 자가일관된 은하 대체 모델을 생성하며, 바의 정렬성과 운동학적 일관성과 같은 관측 제약 조건과 잘 부합한다.