[논문 리뷰] Redback: A Bayesian inference software package for electromagnetic transients

Nikhil Sarin, M. T. Hübner외 9명|arXiv (Cornell University)|2023. 08. 24.

Gamma-ray bursts and supernovaePhysics and Astronomy참고 문헌 130인용 수 3

한 줄 요약

Redback는 전자기 투과성 현상에서 베이지안 추론을 위한 오픈소스, 객체지향적 파이썬 패키지로, 100개 이상의 모델과 12개 이상의 샘플러를 사용하여 카일로노비아, 초신성, 감마선파르트 후광, 기타 투과성 현상에 대한 종단 간 모델링, 시뮬레이션 및 피팅을 가능하게 한다. 주요 카탈로그에서 데이터 수신, 실제 설계된 설문 조사 시뮬레이션(예: LSST, ZTF), 그리고 다중 메신저 추론을 지원하며, Bilby 및 SNCosmo와의 통합을 위한 모듈식 설계를 통해 확장성과 유연성을 제공한다.

ABSTRACT

Fulfilling the rich promise of rapid advances in time-domain astronomy is only possible through confronting our observations with physical models and extracting the parameters that best describe what we see. Here, we introduce {\sc Redback}; a Bayesian inference software package for electromagnetic transients. {\sc Redback} provides an object-orientated {\sc python} interface to over 12 different samplers and over 100 different models for kilonovae, supernovae, gamma-ray burst afterglows, tidal disruption events, engine-driven transients among other explosive transients. The models range in complexity from simple analytical and semi-analytical models to surrogates built upon numerical simulations accelerated via machine learning. {\sc Redback} also provides a simple interface for downloading and processing data from various catalogs such as extit{Swift} and Fink. The software can also serve as an engine to simulate transients for telescopes such as the Zwicky Transient Facility and Vera Rubin with realistic cadences, limiting magnitudes, and sky-coverage or a hypothetical user-constructed survey or a generic transient for target-of-opportunity observations with different telescopes. As a demonstration of its capabilities, we show how {\sc Redback} can be used to jointly fit the spectrum and photometry of a kilonova, enabling a more powerful, holistic probe into the properties of a transient. We also showcase general examples of how {\sc Redback} can be used as a tool to simulate transients for realistic surveys, fit models to real, simulated, or private data, multi-messenger inference with gravitational waves, and serve as an end-to-end software toolkit for parameter estimation and interpreting the nature of electromagnetic transients.

연구 동기 및 목표

• 시간 영역 광도 측정 데이터를 해석하기 위한 유연하고 확장 가능하며 견고한 소프트웨어 툴킷에 대한 증가하는 수요를 해결하기 위해.
• 기존 소프트웨어 패키지의 한계를 극복하기 위해, 예를 들어 제한된 모델 라이브러리, 유연하지 못한 인터페이스, 또는 시뮬레이션과 추론 간의 통합 부족을 해결하기 위해.
• 다양한 투과성 현상에 대해 종단 간 베이지안 추론 워크플로우(모델 피팅, 매개변수 추정, 설문 조사 최적화 포함)를 가능하게 하기 위해.
• GW170817와 같은 중력파 사건과의 공동 추론을 가능하게 하여 다중 메신저 천문학을 지원하기 위해.
• 최신 소프트웨어 개발 관행을 기반으로 한 현대적이고 지속적으로 유지 관리되며 커뮤니티 기반의 오픈소스 플랫폼을 제공하여 천체물리학 애플리케이션에 적합한 솔루션을 제공하기 위해.

제안 방법

• Redback는 모델, 가능도, 사전분포 및 샘플러를 분리하는 객체지향 파이썬 인터페이스를 구현하여, 소스 코드를 수정하지 않고도 사용자가 커스텀 구성 요소를 통합할 수 있도록 한다.
• Bilby와의 통합을 통해 12개 이상의 오픈소스 마르코프 체인 몬테카를로 및 네스티드 샘플링 샘플러를 노출하여 상호 검증과 고성능 추론을 가능하게 한다.
• LSST, ZTF 및 사용자 정의 설문 조사에 대해 실제 관측 조건(주기, 한계 등급, 하늘 덮개)을 재현하는 설문 조사 시뮬레이션 엔진을 포함한다.
• Swift, Fink, Lasair, BATSE 및 오픈 액세스 카탈로그를 포함한 주요 공개 카탈로그에서 데이터를 다운로드하고 처리할 수 있는 내장 인터페이스를 제공한다.
• 모델은 분석적 및 반분석적 표현에서부터 수치 시뮬레이션에 기반한 기계학습 가속 서rogate에 이르기까지 다양하다.
• 표준 과학용 파이썬 라이브러리(astropy, SNCosmo, matplotlib, corner)를 활용하며, 베이지안 계층적 모델링, 모델 선택 및 중요도 샘플링을 지원한다.

Figure 1 : Flowchart showcasing the different modules of Redback and how they interact for different workflows.

실험 결과

연구 질문

• RQ1 모듈식 오픈소스 소프트웨어 패키지가 전자기 투과성 현상 천문학에서의 베이지안 추론의 신뢰성과 재현 가능성에 어떻게 기여할 수 있는가?
• RQ2 Redback는 실제 설문 조사 관측 조건 하에서 얼마나 정확하게 투과성 광도 곡선을 시뮬레이션할 수 있으며, 이를 통해 정확한 설문 조사 설계 및 선택 함수 모델링가 가능할 수 있는가?
• RQ3 Redback의 민첩한 인터페이스는 카일로노비아 및 중력파와 같은 다중 투과성 유형과의 공동 추론을 얼마나 잘 지원할 수 있는가?
• RQ4 CI/CD, 단위 테스트, 활발한 개발과 같은 현대 소프트웨어 관행과의 통합이 기존 패키지에 비해 장기적 유지보수성을 어떻게 향상시키는가?
• RQ5 Redback는 다양한 투과성 유형(예: tidal disruption events, 자석성 초신성 등)에 대해 실제 및 시뮬레이션 데이터를 피팅하는 데 있어 성능과 정확도가 어느 정도인가?

주요 결과

• Redback는 카일로노비아, 초신성, 감마선파르트 후광, 조리개 파괴 사건 등을 포함한 100개 이상의 전자기 투과성 현상 모델에 대해 종단 간 베이지안 추론을 가능하게 하며, 사용자 정의 모델 및 사전분포 정의를 완전히 지원한다.
• LSST 및 ZTF에 대해 실제 설문 조사 조건(주기, 깊이, 하늘 덮개)을 반영한 투과성 시뮬레이션을 성공적으로 수행하여 정확한 설문 조사 설계 및 선택 함수 모델링이 가능하다.
• Redback는 실질적인 연구에 활용되어 Sarin & Metzger(2023)의 조리개 파괴 사건, Schulze 등(2023)의 SN 2018ibb, Levan 등(2023)의 공동 GRB-카일로노비아 사건에 대한 매개변수 추정에 사용되었으며, 실제 응용 가능성에 대한 증거를 제공한다.
• Bilby와의 통합을 통해 다양한 샘플러의 원활한 사용이 가능해져 상호 검증과 효율적인 추론이 가능하며, 병렬 처리 및 고급 샘플링 기법 지원이 가능하다.
• 소프트웨어의 모듈식 아키텍처와 활발한 개발 팀 덕분에 장기적인 유지보수성과 확장성이 보장되며, 커뮤니티 기여를 통해 새로운 모델과 데이터 유형에 대한 지원이 가능하다.
• 모든 결과, 예제 스크립트, 플로팅 노트북은 GitHub 및 Zenodo(DOI: 10.5281/zenodo.8273145)를 통해 공개되어 있으며, 패키지는 PyPI에 배포되고, 자세한 문서는 https://redback.readthedocs.io/en/latest/ 에서 제공된다.

Figure 2 : (Left) Simulated kilonova (one component kilonova model) in ZTF and in LSST (Right). We emphasize that aesthetic features such as the colours of the data points, axes limits etc can all be modified by passing in relevant keyword arguments to the plotting methods.