[논문 리뷰] BGM FASt: Besançon Galaxy Model for Big Data. Simultaneous inference of the IMF, SFH and density in the Solar Neighbourhood

R. Mor, A. C. Robin외 2명|HAL (Le Centre pour la Communication Scientifique Directe)|2018. 09. 10.

Stellar, planetary, and galactic studies참고 문헌 58인용 수 11

한 줄 요약

이 논문은 베사종 은하모델 기반으로, 사전에 샘플링된 시뮬레이션의 재가중치화와 근사 베이지안 계산을 활용하여 태양계 부근의 초기 질량 함수(IMF), 은하형성 역사(SFH), 항성 밀도를 동시에 추론하는 빠르고 대용량 데이터 기반의 프레임워크인 BGM FASt를 소개한다. 표준 BGM 대비 약 10⁴배 빠른 성능을 달성하며, 현재의 은하형성률은 1.2±0.2 M⊙/yr, 총 항성 밀도는 0.051⁻⁰.⁰⁰⁵⁺⁰.⁰⁰² M⊙/pc³, 0.5–1.53 M⊙ 범위 내 IMF 기울기는 α₂=2.1⁻⁰.³⁺⁰.¹로 추정된다.

ABSTRACT

We develop a new theoretical framework to generate Besançon Galaxy Model fast approximate simulations (BGM FASt) to address fundamental questions of the Galactic structure and evolution performing multi-parameter inference. As a first application of our strategy we simultaneously infer the IMF, the star formation history and the stellar mass density in the Solar Neighbourhood. The BGM FASt strategy is based on a reweighing scheme, that uses a specific pre-sampled simulation, and on the assumption that the distribution function of the generated stars in the Galaxy can be described by an analytical expression. To validate BGM FASt we execute a set of tests. Finally, we use BGM FASt with an approximate Bayesian computation algorithm to obtain the posterior PDF of the inferred parameters, by comparing synthetic versus Tycho-2 colour-magnitude diagrams. Results: The validation shows a very good agreement between BGM FASt and the standard BGM, with BGM FASt being $\approx 10^4$ times faster. By analysing Tycho-2 data we obtain a thin disc star formation history decreasing in time and a present rate of $1.2 \pm 0.2 M_\odot/yr$. The resulting total stellar mass density in the Solar Neighbourhood is $0.051_{-0.005}^{+0.002} M_\odot/pc^3$ and the local dark matter density is $0.012 \pm 0.001 M_\odot/pc^3$. For the composite IMF we obtain a slope of $α_2={2.1}_{-0.3}^{+0.1}$ in the mass range between $0.5 M_\odot$ and $1.53M_\odot$. The results of the slope at the high mass range are trustable up to $4M_\odot$ and highly depend on the choice of the extinction map (obtaining $α_3={2.9}_{-0.2}^{+0.2}$ and $α_3={3.7}_{-0.2}^{+0.2}$ respectively, for two different extinction maps). Systematic uncertainties are not included. Conclusions: The good performance of BGM FASt demonstrates that it is a very valuable tool to perform multi-parameter inference using Gaia data releases.

연구 동기 및 목표

• 기존 은하 시뮬레이션 코드로 타원은하 모델링에서 큰 매개변수 공간을 탐색할 때 발생하는 계산적 병목 현상을 해결하기 위해.
• 대용량 데이터 플랫폼에 적합한 확장성 있고 고성능인 빠른 근사 시뮬레이션을 위한 프레임워크를 개발하기 위해.
• 관측 데이터를 활용해 태양계 부근의 초기 질량 함수(IMF), 은하형성 역사(SFH), 항성 질량 밀도를 동시에 추론할 수 있도록 하기 위해.
• 표준 BGM 시뮬레이션과 실제 Tycho-2 데이터를 기준으로 프레임워크의 정확성과 효율성을 검증하기 위해.

제안 방법

• BGM FASt 프레임워크는 사전에 샘플링된 표준 BGM 시뮬레이션에 재가중치화 기법을 적용하여 빠른 근사 시뮬레이션을 생성한다.
• 은하 내 항성 분포 함수가 해석적 표현으로 기술될 수 있다고 가정하여, 효율적인 매개변수 공간 탐색이 가능하다.
• 분산 컴퓨팅을 위해 Apache Spark와 Hadoop를 활용하여 대용량 데이터 인프라에서의 확장성을 확보한다.
• 합성 색-등급도(CMD)를 관측된 Tycho-2 데이터와 비교하기 위해 근사 베이지안 계산(ABC)을 사용하여 사후 분포를 추론한다.
• IMF, SFH, 밀도 법칙에 대한 민감한 매개변수화를 가능하게 하며, 비모수적 형태도 포함한다.
• 식 (35)를 사용해 새로운 모델 매개변수에 기반해 사전에 샘플링된 시뮬레이션을 재가중시키며, 전체 재시뮬레이션을 피한다.

실험 결과

연구 질문

• RQ1 대용량 데이터 플랫폼을 활용해 은하 모델의 다중 매개변수 추론을 수행할 수 있는 빠르고 확장 가능한 프레임워크를 개발할 수 있는가?
• RQ2 Tycho-2 데이터를 기반으로 태양계 부근의 초기 질량 함수(IMF), 은하형성 역사(SFH), 항성 질량 밀도에 대한 사후 제약 조건은 무엇인가?
• RQ3 3차원 먼지 투과도 지도의 선택이 고질량 영역에서 추론된 IMF 기울기에 어떤 영향을 미치는가?
• RQ4 표준 BGM 대비 10⁴배 빠른 성능을 달성하면서도 정확도를 얼마나 유지할 수 있는가?
• RQ5 프레임워크를 비모수적 SFH 및 더 복잡한 IMF 매개변수화로 확장할 수 있는가?

주요 결과

• BGM FASt 프레임워크는 표준 BGM 코드 대비 약 10⁴배 빠른 성능을 달성하면서도 검증 테스트에서 뛰어난 일치를 보였다.
• 태양계 부근의 현재 은하형성률은 1.2±0.2 M⊙/yr로 추정되며, 사용된 먼지 투과도 지도에 관계없이 동일하다.
• 태양계 부근의 총 항성 체적 질량 밀도는 0.051⁻⁰.⁰⁰⁵⁺⁰.⁰⁰² M⊙/pc³이며, 지역 암흑물질 밀도는 0.012±0.001 M⊙/pc³이다.
• 0.5–1.53 M⊙ 범위 내 IMF 기울기는 α₂=2.1⁻⁰.³⁺⁰.¹로 추정되며, 이는 살페터 IMF와 3σ 수준에서 불일치한다.
• 고질량 영역(4 M⊙)에서의 IMF 기울기는 먼지 투과도 지도에 매우 민감하여, Stilism 지도에서는 α₃=2.9⁻⁰.²⁺⁰.², Drimmel 지도에서는 α₃=3.7⁻⁰.²⁺⁰.²로 도출되었으며, 둘 다 살페터 기울기보다 더 기울어져 있다.
• 박층의 반경 스케일 길이는 잘 제약되지 않았다(2151⁻²⁴⁷⁺⁴²¹ pc), 반면 SFH 감쇠 시간 상수 γ는 Drimmel 기준 γ=0.09⁻⁰.²⁺⁰.³, Stilism 기준 γ=0.14⁻⁰.²⁺⁰.³로 도출되었으며, 형태가 다름에도 불구하고 현재의 SFR은 일관되게 유사하게 추정된다.