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

R. Mor, A. C. Robin외 2명|arXiv (Cornell University)|2018. 09. 10.

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

한 줄 요약

이 논문은 약간의 근사 베이지안 계산(ABC)을 기반으로 한 BGM FASt를 소개한다. BGM FASt는 사전에 샘플링된 시뮬레이션과 재가중치화 기법을 활용하여 태양계 부근에서 초기 질량 함수(IMF), 은하계 형성 역사(SFH), 별 밀도를 동시에 다중 매개변수 추론할 수 있는 빠르고 확장 가능한 프레임워크이다. 표준 BGM에 비해 약 10⁴배 빠른 성능을 달성하면서도, 색-등급도와 같은 주요 관측 결과를 정확히 재현한다. 이로써 얇은 디스크의 SFH는 시간이 지남에 따라 감소하는 경향을 보이며, 현재의 별 형성률은 1.2±0.2 M⊙/yr, 총 별 질량 밀도는 0.051+0.002−0.005 M⊙/pc³이며, 0.5에서 1.53 M⊙ 사이의 복합 IMF 기울기는 α₂=2.1+0.1−0.3이다.

ABSTRACT

We develop a new theoretical framework to generate Besan\c{c}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 $\alpha_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 $\alpha_3={2.9}_{-0.2}^{+0.2}$ and $\alpha_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.

연구 동기 및 목표

• 현대 베이지안 추론을 사용할 때 대규모 매개변수 공간을 탐색할 때 전통적인 Milky Way 모델링의 계산적 한계를 해결하기 위해.
• Gaia와 같은 대규모 설문 조사 데이터를 처리할 수 있도록 Apache Spark 및 Hadoop 같은 대용량 데이터 플랫폼과 호환되는 빠르고 확장 가능한 시뮬레이션 프레임워크를 개발하기 위해.
• Tycho-2의 관측 데이터를 활용해 초기 질량 함수(IMF), SFH, 별 질량 밀도와 같은 여러 핵심 은하 매개변수를 동시에 추론할 수 있도록 하기 위해.
• 표준 BGM 시뮬레이션과의 성능를 검증하면서도 역학적 일관성을 유지하기 위해.
• 실제 광학 및 천체측량 데이터를 활용해 국지적 IMF, SFH, 밀도를 제약하는 데 프레임워크의 유용성을 입증하기 위해.

제안 방법

• BGM FASt 프레임워크는 표준 Besançon Galaxy Model(BGM Std)를 기반으로 한 사전에 샘플링된 시뮬레이션(모더 시뮬레이션)을 기초로 한다.
• 전체 시뮬레이션을 다시 실행하지 않고도 새로운 매개변수 값(예: IMF, SFH, 밀도 법칙)에 따라 시뮬레이션된 별 집단을 재가중치화하는 기법을 적용한다.
• 은하계 내 별의 분포 함수가 분석적 표현으로 기술될 수 있다고 가정함으로써, 효율적인 매개변수 공간 탐색이 가능해진다.
• 합성 색-등급도(CMD)와 관측된 Tycho-2 데이터를 비교하기 위해 약간의 근사 베이지안 계산(ABC) 알고리즘과 통합하여 자동으로 매개변수 추론을 수행한다.
• Apache Spark 및 Hadoop에 구현되어 Gaia 데이터 분석 프레임워크(GDAF)와 같은 대용량 데이터 인프라에서 분산 고성능 계산을 가능하게 한다.
• 역학적 일관성을 유지하기 위해 분석적 제약 조건을 통해 국소 역학 통계 평형(LDSE)이 유지되도록 하며, 이는 전체 BGM Std 시뮬레이션과의 비교를 통해 검증된다.

실험 결과

연구 질문

• RQ1 대규모 설문 조사 데이터를 활용해 은하 매개변수의 다중 매개변수 추론을 효율적으로 수행할 수 있는 빠르고 확장 가능한 시뮬레이션 프레임워크를 개발할 수 있는가?
• RQ2 BGM FASt는 표준 BGM에 비해 주요 관측 결과인 CMD와 별 분포를 얼마나 정확하게 재현할 수 있는가?
• RQ3 Tycho-2 데이터를 활용해 태양계 부근에서의 초기 질량 함수(IMF), 은하계 형성 역사(SFH), 별 질량 밀도의 추정치는 무엇인가?
• RQ4 특히 흡수 지도의 시스템적 불확실성은 고질량에서의 IMF 기울기에 어떤 영향을 미치는가?
• RQ5 비모수적 SFH, 금속성에 의존하는 IMF, 또는 시간에 따라 변화하는 디스크 구조를 포함하는 데에 이 프레임워크를 얼마나 넓게 확장할 수 있는가?

주요 결과

• BGM FASt는 표준 BGM 코드 대비 약 10⁴배 빠른 성능을 달성하면서도 시뮬레이션된 색-등급도와 별 분포에서 뛰어난 일치를 유지한다.
• 태양계 부근에서 추정된 별 형성 역사는 시간이 지남에 따라 감소하는 경향을 보이며, 현재의 별 형성률은 1.2±0.2 M⊙/yr이다.
• 태양계 부근의 총 별 질량 밀도는 0.051+0.002−0.005 M⊙/pc³이며, 국소 암흑 물질 밀도는 0.012±0.001 M⊙/pc³이다.
• 복합 IMF의 경우 0.5–1.53 M⊙ 범위에서의 기울기는 α₂=2.1+0.1−0.3이며, Salpeter 기울기와 일치한다.
• 고질량 영역(4 M⊙)에서의 IMF 기울기는 흡수 지도에 매우 민감하다: Drimmel 지도를 사용할 경우 α₃=3.7+0.2−0.2, Stilism 지도를 사용할 경우 α₃=2.9+0.2−0.2이며, 둘 다 ~3σ에서 Salpeter 기울기를 기각한다.
• 얇은 디스크의 반경 스케일 길이는 잘 제약되지 않았다(예: Stilism 지도를 사용할 경우 2151+421−247 pc), SFH의 형태는 흡수 지도 선택과 상관관계가 있으며, 현재의 별 형성률은 여전히 강건하다.