[논문 리뷰] New mass estimates for massive binary systems: a probabilistic approach using polarimetric radiative transfer

Andrew G. Fullard, Jack O’Brien외 5명|arXiv (Cornell University)|2022. 02. 21.

Gamma-ray bursts and supernovae참고 문헌 67인용 수 5

한 줄 요약

이 논문은 광학 두께가 두꺼운 바람을 가진 질량이 큰 월프-레인어 + O 이중성에서 질량을 추정하기 위한 새로운 확률적 방법을 제시한다. 이 방법은 신경망에 의한 몽테카를로 복사전달 모델을 사용하여 분석 모델의 편향을 극복한다. 이 접근법은 더 정확한 경사각 추정을 가능하게 하여, 특히 WR 153의 경우 이전 추정치와 상당한 불일치를 보이며, 질량 추정치를 개선하기 위해 더 나은 극화 및 스펙트로스코픽 데이터가 필요함을 시사한다. 이는 질량이 큰 이중성의 진화 모델을 정교화하는 데 기여한다.

ABSTRACT

Understanding the evolution of massive binary stars requires accurate estimates of their masses. This understanding is critically important because massive star evolution can potentially lead to gravitational wave sources such as binary black holes or neutron stars. For Wolf-Rayet stars with optically thick stellar winds, their masses can only be determined with accurate inclination angle estimates from binary systems which have spectroscopic $M \sin i$ measurements. Orbitally-phased polarization signals can encode the inclination angle of binary systems, where the Wolf-Rayet winds act as scattering regions. We investigated four Wolf-Rayet + O star binary systems, WR 42, WR 79, WR 127, and WR 153, with publicly available phased polarization data to estimate their masses. To avoid the biases present in analytic models of polarization while retaining computational expediency, we used a Monte Carlo radiative transfer model accurately emulated by a neural network. We used the emulated model to investigate the posterior distribution of parameters of our four systems. Our mass estimates calculated from the estimated inclination angles put strong constraints on existing mass estimates for three of the systems, and disagrees with the existing mass estimates for WR 153. We recommend a concerted effort to obtain polarization observations that can be used to estimate the masses of Wolf-Rayet binary systems and increase our understanding of their evolutionary paths.

연구 동기 및 목표

• 질량이 큰 이중성계, 특히 월프-레인어 + O 별에서 전통적인 방법이 편향과 불확실성으로 인해 어려움을 겪는 점을 해결하기 위해 질량 추정치를 향상시키는 것.
• 분석적 극화 모델의 한계를 해결하기 위해, 이는 고경사각에 대해 편향되며 광학적으로 얇은 산란 영역에 국한된다는 점이다.
• 단계별 극화 및 복사전달을 이용한 궤도 기울기와 질량을 추정하는 계산 효율성과 정확성이 높은 방법을 개발하는 것.
• 기울기, 광학 두께, 복사율 비율 등의 모델 매개변수에 대한 강력한 사후 분포를 제공하고, 불확실성을 정량화하는 것.
• 공개된 극화 데이터를 사용해 네 개의 WR + O 시스템을 재분석함으로써 기존 분석 방법에 비해 본 방법의 우수성을 입증하는 것.

제안 방법

• 광학 두께가 두꺼운 바람을 가진 이중성계에서 편광된 빛을 시뮬레이션하기 위해 몽테카를로 복사전달(MCRT) 코드인 slip을 사용한다.
• 전체 복사전달 모델을 근사하기 위해 신경망 에뮬레이터를 훈련시켜, 빠른 샘플링이 가능하도록 MCRT 시뮬레이션을 가속화한다.
• 에뮬레이터를 사용하여 베이지안 추론를 적용하여 기울기 각도, 광학 두께, 복사율 비율 등의 모델 매개변수의 사후 분포를 탐색한다.
• 시스템 매개변수에 대한 사전 지식을 통합하고, 신뢰구간을 추론하기 위해 마르코프 체인 몽테카를로(MCMC) 샘플링을 사용한다.
• 훈련 범위 내에서 정확도를 확보하기 위해 에뮬레이터를 전체 MCRT 시뮬레이션과 검증한다.
• 광학 두께(τ), 바람 비대칭성(zell), 복사율 비율(Ifrac) 등의 모델 매개변수를 함께 제약 조건을 둔다. 이는 강한 디그레시를 드러낸다.

실험 결과

연구 질문

• RQ1 신경망에 의해 에뮬레이트된 복사전달 모델은 기존 분석 모델에 비해 질량이 큰 이중성계에서 더 정확하고 편향이 없는 기울기 각도 추정치를 제공할 수 있는가?
• RQ2 극화 데이터로부터 재조정된 기울기 각도 추정치는 월프-레인어 + O 이중성계의 유도된 질량에 어떤 영향을 미치는가?
• RQ3 광학 두께, 바람 비대칭성, 복사율 비율 등의 핵심 매개변수의 사후 분포는 어떻게 되며, 이들 간의 상관관계는 어떠한가?
• RQ4 새로운 질량 추정치는 기존 추정치와 얼마나 다를까? 특히 WR 153의 경우 어떻게 되며, 이는 진화 모델에 어떤 함의를 갖는가?
• RQ5 이 방법은 타원도가 있는 궤도 또는 B형 별 동반성을 가진 다른 질량이 큰 이중성계로 일반화될 수 있는가?

주요 결과

• 이 방법은 단계별 극화 및 에뮬레이트된 복사전달을 이용해 고정된 기울기 각도를 정밀하게 추정함으로써 M sin i의 불확실성 문제를 해결한다.
• WR 153의 경우 새로운 질량 추정치가 이전 극화 추정치와 상당한 불일치를 보이며, 이는 이전 모델이 편향되었음을 시사한다.
• 모든 네 개의 시스템에서 오르막 절점의 장위치(Ω)에 대한 사후 분포는 이전 추정치와 완전히 불일치하며, 이는 적절한 불확실성 정량화가 부족했음을 나타낸다.
• 광학 두께(τ)는 일반적으로 낮은 수준(약 0.2)이며, 이는 밴드폭 극화가 월프-레인어 별 바람의 외부, 덜 조밀한 영역에서 발생함을 의미한다.
• τ와 zell 간, Ifrac와 τ 간의 강한 상관관계가 발견되어, 해석 시 고려해야 할 디그레시가 있음을 드러낸다.
• 이 방법은 모든 시스템에서 질량에 대한 강력한 상한선을 제공하고 불확실성을 정량화하여, 분석 모델에 비해 통계적으로 타당한 대안을 제공한다.